Integrative Analyses of Tumor and Peripheral Biomarkers in the Treatment of Advanced Renal Cell Carcinoma.

The phase III JAVELIN Renal 101 trial demonstrated prolonged progression-free survival (PFS) in patients (N = 886) with advanced renal cell carcinoma treated with first-line avelumab + axitinib (A+Ax) versus sunitinib. We report novel findings from integrated analyses of longitudinal blood samples and baseline tumor tissue. PFS was associated with elevated lymphocyte levels in the sunitinib arm and an abundance of innate immune subsets in the A+Ax arm. Treatment with A+Ax led to greater T-cell repertoire modulation and less change in T-cell numbers versus sunitinib. In the A+Ax arm, patients with tumors harboring mutations in ≥2 of 10 previously identified PFS-associated genes (double mutants) had distinct circulating and tumor-infiltrating immunologic profiles versus those with wild-type or single-mutant tumors, suggesting a role for non-T-cell-mediated and non-natural killer cell-mediated mechanisms in double-mutant tumors. We provide evidence for different immunomodulatory mechanisms based on treatment (A+Ax vs. sunitinib) and tumor molecular subtypes. SIGNIFICANCE: Our findings provide novel insights into the different immunomodulatory mechanisms governing responses in patients treated with avelumab (PD-L1 inhibitor) + axitinib or sunitinib (both VEGF inhibitors), highlighting the contribution of tumor biology to the complexity of the roles and interactions of infiltrating immune cells in response to these treatment regimens. This article is featured in Selected Articles from This Issue, p. 384.

HLA-A*03 and response to immune checkpoint blockade in cancer: an epidemiological biomarker study.

BACKGROUND: Predictive biomarkers could allow more precise use of immune checkpoint inhibitors (ICIs) in treating advanced cancers. Given the central role of HLA molecules in immunity, variation at the HLA loci could differentially affect the response to ICIs. The aim of this epidemiological study was to determine the effect of HLA-A*03 as a biomarker for predicting response to immunotherapy. METHODS: In this epidemiological study, we investigated the clinical outcomes (overall survival, progression free survival, and objective response rate) after treatment for advanced cancer in eight cohorts of patients: three observational cohorts of patients with various types of advanced tumours (the Memorial Sloan Kettering Integrated Mutation Profiling of Actionable Cancer Targets [MSK-IMPACT] cohort, the Dana-Farber Cancer Institute [DFCI] Profile cohort, and The Cancer Genome Atlas) and five clinical trials of patients with advanced bladder cancer (JAVELIN Solid Tumour) or renal cell carcinoma (CheckMate-009, CheckMate-010, CheckMate-025, and JAVELIN Renal 101). In total, these cohorts included 3335 patients treated with various ICI agents (anti-PD-1, anti-PD-L1, and anti-CTLA-4 inhibitors) and 10 917 patients treated with non-ICI cancer-directed therapeutic approaches. We initially modelled the association of HLA amino-acid variation with overall survival in the MSK-IMPACT discovery cohort, followed by a detailed analysis of the association between HLA-A*03 and clinical outcomes in MSK-IMPACT, with replication in the additional cohorts (two further observational cohorts and five clinical trials). FINDINGS: HLA-A*03 was associated in an additive manner with reduced overall survival after ICI treatment in the MSK-IMPACT cohort (HR 1·48 per HLA-A*03 allele [95% CI 1·20-1·82], p=0·00022), the validation DFCI Profile cohort (HR 1·22 per HLA-A*03 allele, 1·05-1·42; p=0·0097), and in the JAVELIN Solid Tumour clinical trial for bladder cancer (HR 1·36 per HLA-A*03 allele, 1·01-1·85; p=0·047). The HLA-A*03 effect was observed across ICI agents and tumour types, but not in patients treated with alternative therapies. Patients with HLA-A*03 had shorter progression-free survival in the pooled patient population from the three CheckMate clinical trials of nivolumab for renal cell carcinoma (HR 1·31, 1·01-1·71; p=0·044), but not in those receiving control (everolimus) therapies. Objective responses were observed in none of eight HLA-A*03 homozygotes in the ICI group (compared with 59 [26·6%] of 222 HLA-A*03 non-carriers and 13 (17·1%) of 76 HLA-A*03 heterozygotes). HLA-A*03 was associated with shorter progression-free survival in patients receiving ICI in the JAVELIN Renal 101 randomised clinical trial for renal cell carcinoma (avelumab plus axitinib; HR 1·59 per HLA-A*03 allele, 1·16-2·16; p=0·0036), but not in those receiving control (sunitinib) therapy. Objective responses were recorded in one (12·5%) of eight HLA-A*03 homozygotes in the ICI group (compared with 162 [63·8%] of 254 HLA-A*03 non-carriers and 40 [55·6%] of 72 HLA-A*03 heterozygotes). HLA-A*03 was associated with impaired outcome in meta-analysis of all 3335 patients treated with ICI at genome-wide significance (p=2·01 × 10-8) with no evidence of heterogeneity in effect (I2 0%, 95% CI 0-0·76) INTERPRETATION: HLA-A*03 is a predictive biomarker of poor response to ICI. Further evaluation of HLA-A*03 is warranted in randomised trials. HLA-A*03 carriage could be considered in decisions to initiate ICI in patients with cancer. FUNDING: National Institutes of Health, Merck KGaA, and Pfizer.

Avelumab maintenance in advanced urothelial carcinoma: biomarker analysis of the phase 3 JAVELIN Bladder 100 trial.

In a recent phase 3 randomized trial of 700 patients with advanced urothelial cancer (JAVELIN Bladder 100; NCT02603432 ), avelumab/best supportive care (BSC) significantly prolonged overall survival relative to BSC alone as maintenance therapy after first-line chemotherapy. Exploratory biomarker analyses were performed to identify biological pathways that might affect survival benefit. Tumor molecular profiling by immunohistochemistry, whole-exome sequencing and whole-transcriptome sequencing revealed that avelumab survival benefit was positively associated with PD-L1 expression by tumor cells, tumor mutational burden, APOBEC mutation signatures, expression of genes underlying innate and adaptive immune activity and the number of alleles encoding high-affinity variants of activating Fcγ receptors. Pathways connected to tissue growth and angiogenesis might have been associated with reduced survival benefit. Individual biomarkers did not comprehensively identify patients who could benefit from therapy; however, multi-parameter models incorporating genomic alteration, immune responses and tumor growth showed promising predictive utility. These results characterize the complex biologic pathways underlying survival benefit from immune checkpoint inhibition in advanced urothelial cancer and suggest that multiple biomarkers might be needed to identify patients who would benefit from treatment.

Expanding control of the tumor cell cycle with a CDK2/4/6 inhibitor.

The CDK4/6 inhibitor, palbociclib (PAL), significantly improves progression-free survival in HR+/HER2- breast cancer when combined with anti-hormonals. We sought to discover PAL resistance mechanisms in preclinical models and through analysis of clinical transcriptome specimens, which coalesced on induction of MYC oncogene and Cyclin E/CDK2 activity. We propose that targeting the G1 kinases CDK2, CDK4, and CDK6 with a small-molecule overcomes resistance to CDK4/6 inhibition. We describe the pharmacodynamics and efficacy of PF-06873600 (PF3600), a pyridopyrimidine with potent inhibition of CDK2/4/6 activity and efficacy in multiple in vivo tumor models. Together with the clinical analysis, MYC activity predicts (PF3600) efficacy across multiple cell lineages. Finally, we find that CDK2/4/6 inhibition does not compromise tumor-specific immune checkpoint blockade responses in syngeneic models. We anticipate that (PF3600), currently in phase 1 clinical trials, offers a therapeutic option to cancer patients in whom CDK4/6 inhibition is insufficient to alter disease progression.

Molecular medicine tumor board: whole-genome sequencing to inform on personalized medicine for a man with advanced prostate cancer.

PURPOSE: Molecular profiling of cancer is increasingly common as part of routine care in oncology, and germline and somatic profiling may provide insights and actionable targets for men with metastatic prostate cancer. However, all reported cases are of deidentified individuals without full medical and genomic data available in the public domain. PATIENT AND METHODS: We present a case of whole-genome tumor and germline sequencing in a patient with advanced prostate cancer, who has agreed to make his genomic and clinical data publicly available. RESULTS: We describe an 84-year-old Caucasian male with a Gleason 10 oligometastastic hormone-sensitive prostate cancer. Whole-genome sequencing provided insights into his tumor's underlying mutational processes and the development of an SPOP mutation. It also revealed an androgen-receptor dependency of his cancer which was reflected in his durable response to radiation and hormonal therapy. Potentially actionable genomic lesions in the tumor were identified through a personalized medicine approach for potential future therapy, but at the moment, he remains in remission, illustrating the hormonal sensitivity of his SPOP-driven prostate cancer. We also placed this patient in the context of a large prostate-cancer cohort from the PCAWG (Pan-cancer Analysis of Whole Genomes) group. In this comparison, the patient's cancer appears typical in terms of the number and type of somatic mutations, but it has a somewhat larger contribution from the mutational process associated with aging. CONCLUSION: We combined the expertise of medical oncology and genomics approaches to develop a molecular tumor board to integrate the care and study of this patient, who continues to have an outstanding response to his combined modality treatment. This identifiable case potentially helps overcome barriers to clinical and genomic data sharing.

Avelumab plus axitinib versus sunitinib in advanced renal cell carcinoma: biomarker analysis of the phase 3 JAVELIN Renal 101 trial.

We report on molecular analyses of baseline tumor samples from the phase 3 JAVELIN Renal 101 trial (n = 886; NCT02684006 ), which demonstrated significantly prolonged progression-free survival (PFS) with first-line avelumab + axitinib versus sunitinib in advanced renal cell carcinoma (aRCC). We found that neither expression of the commonly assessed biomarker programmed cell death ligand 1 (PD-L1) nor tumor mutational burden differentiated PFS in either study arm. Similarly, the presence of FcɣR single nucleotide polymorphisms was unimpactful. We identified important biological features associated with differential PFS between the treatment arms, including new immunomodulatory and angiogenesis gene expression signatures (GESs), previously undescribed mutational profiles and their corresponding GESs, and several HLA types. These findings provide insight into the determinants of response to combined PD-1/PD-L1 and angiogenic pathway inhibition and may aid in the development of strategies for improved patient care in aRCC.

Landscape of somatic single nucleotide variants and indels in colorectal cancer and impact on survival.

Colorectal cancer (CRC) is a biologically heterogeneous disease. To characterize its mutational profile, we conduct targeted sequencing of 205 genes for 2,105 CRC cases with survival data. Our data shows several findings in addition to enhancing the existing knowledge of CRC. We identify PRKCI, SPZ1, MUTYH, MAP2K4, FETUB, and TGFBR2 as additional genes significantly mutated in CRC. We find that among hypermutated tumors, an increased mutation burden is associated with improved CRC-specific survival (HR = 0.42, 95% CI: 0.21-0.82). Mutations in TP53 are associated with poorer CRC-specific survival, which is most pronounced in cases carrying TP53 mutations with predicted 0% transcriptional activity (HR = 1.53, 95% CI: 1.21-1.94). Furthermore, we observe differences in mutational frequency of several genes and pathways by tumor location, stage, and sex. Overall, this large study provides deep insights into somatic mutations in CRC, and their potential relationships with survival and tumor features.

Standard machine learning approaches outperform deep representation learning on phenotype prediction from transcriptomics data.

BACKGROUND: The ability to confidently predict health outcomes from gene expression would catalyze a revolution in molecular diagnostics. Yet, the goal of developing actionable, robust, and reproducible predictive signatures of phenotypes such as clinical outcome has not been attained in almost any disease area. Here, we report a comprehensive analysis spanning prediction tasks from ulcerative colitis, atopic dermatitis, diabetes, to many cancer subtypes for a total of 24 binary and multiclass prediction problems and 26 survival analysis tasks. We systematically investigate the influence of gene subsets, normalization methods and prediction algorithms. Crucially, we also explore the novel use of deep representation learning methods on large transcriptomics compendia, such as GTEx and TCGA, to boost the performance of state-of-the-art methods. The resources and findings in this work should serve as both an up-to-date reference on attainable performance, and as a benchmarking resource for further research. RESULTS: Approaches that combine large numbers of genes outperformed single gene methods consistently and with a significant margin, but neither unsupervised nor semi-supervised representation learning techniques yielded consistent improvements in out-of-sample performance across datasets. Our findings suggest that using l2-regularized regression methods applied to centered log-ratio transformed transcript abundances provide the best predictive analyses overall. CONCLUSIONS: Transcriptomics-based phenotype prediction benefits from proper normalization techniques and state-of-the-art regularized regression approaches. In our view, breakthrough performance is likely contingent on factors which are independent of normalization and general modeling techniques; these factors might include reduction of systematic errors in sequencing data, incorporation of other data types such as single-cell sequencing and proteomics, and improved use of prior knowledge.

TIME (Tumor Immunity in the MicroEnvironment) classification based on tumor CD274 (PD-L1) expression status and tumor-infiltrating lymphocytes in colorectal carcinomas.

Inhibitors targeting the PDCD1 (programmed cell death 1, PD-1) immune checkpoint pathway have revolutionized cancer treatment strategies. The TIME (Tumor Immunity in the MicroEnvironment) classification based on tumor CD274 (PDCD1 ligand 1, PD-L1) expression and tumor-infiltrating lymphocytes (TIL) has been proposed to predict response to immunotherapy. It remains to be determined clinical, pathological, and molecular features of TIME subtypes of colorectal cancer. Using 812 colon and rectal carcinoma cases from the Nurses' Health Study and Health Professionals Follow-up Study, we examined the association of tumor characteristics and survival outcomes with four TIME subtypes (TIME 1, CD274low/TILabsent; TIME 2, CD274high/TILpresent; TIME 3, CD274low/TILpresent; and TIME 4, CD274high/TILabsent). In survival analyses, Cox proportional hazards models were adjusted for potential confounders, including microsatellite instability (MSI) status, CpG island methylator phenotype (CIMP) status, LINE-1 methylation level, and KRAS, BRAF, and PIK3CA mutation status. TIME subtypes 1, 2, 3 and 4 had 218 (27%), 117 (14%), 103 (13%), and 374 (46%) colorectal cancer cases, respectively. Compared with TIL-absent subtypes (TIME 1 and 4), TIL-present subtypes (TIME 2 and 3) were associated with high-level MSI, high-degree CIMP, BRAF mutation, and higher amounts of neoantigens (p < 0.001). TIME subtypes were not significantly associated with colorectal cancer-specific or overall survival. In conclusion, TIL-present TIME subtypes of colorectal cancer are associated with high levels of MSI and neoantigen load, supporting better responsiveness to cancer immunotherapy. Further studies examining tumor molecular alterations and additional factors in the tumor microenvironment may inform development of immunoprevention and immunotherapy strategies.

Multi-omics profiling of younger Asian breast cancers reveals distinctive molecular signatures.

Breast cancer (BC) in the Asia Pacific regions is enriched in younger patients and rapidly rising in incidence yet its molecular bases remain poorly characterized. Here we analyze the whole exomes and transcriptomes of 187 primary tumors from a Korean BC cohort (SMC) enriched in pre-menopausal patients and perform systematic comparison with a primarily Caucasian and post-menopausal BC cohort (TCGA). SMC harbors higher proportions of HER2+ and Luminal B subtypes, lower proportion of Luminal A with decreased ESR1 expression compared to TCGA. We also observe increased mutation prevalence affecting BRCA1, BRCA2, and TP53 in SMC with an enrichment of a mutation signature linked to homologous recombination repair deficiency in TNBC. Finally, virtual microdissection and multivariate analyses reveal that Korean BC status is independently associated with increased TIL and decreased TGF-ß signaling expression signatures, suggesting that younger Asian BCs harbor more immune-active microenvironment than western BCs.

Genetic Mechanisms of Immune Evasion in Colorectal Cancer.

To understand the genetic drivers of immune recognition and evasion in colorectal cancer, we analyzed 1,211 colorectal cancer primary tumor samples, including 179 classified as microsatellite instability-high (MSI-high). This set includes The Cancer Genome Atlas colorectal cancer cohort of 592 samples, completed and analyzed here. MSI-high, a hypermutated, immunogenic subtype of colorectal cancer, had a high rate of significantly mutated genes in important immune-modulating pathways and in the antigen presentation machinery, including biallelic losses of B2M and HLA genes due to copy-number alterations and copy-neutral loss of heterozygosity. WNT/ß-catenin signaling genes were significantly mutated in all colorectal cancer subtypes, and activated WNT/ß-catenin signaling was correlated with the absence of T-cell infiltration. This large-scale genomic analysis of colorectal cancer demonstrates that MSI-high cases frequently undergo an immunoediting process that provides them with genetic events allowing immune escape despite high mutational load and frequent lymphocytic infiltration and, furthermore, that colorectal cancer tumors have genetic and methylation events associated with activated WNT signaling and T-cell exclusion.Significance: This multi-omic analysis of 1,211 colorectal cancer primary tumors reveals that it should be possible to better monitor resistance in the 15% of cases that respond to immune blockade therapy and also to use WNT signaling inhibitors to reverse immune exclusion in the 85% of cases that currently do not. Cancer Discov; 8(6); 730-49. ©2018 AACR.This article is highlighted in the In This Issue feature, p. 663.

Inherited DNA-Repair Defects in Colorectal Cancer.

Colorectal cancer (CRC) heritability has been estimated to be around 30%. However, mutations in the known CRC-susceptibility genes explain CRC risk in fewer than 10% of affected individuals. Germline mutations in DNA-repair genes (DRGs) have recently been reported in CRC, but their contribution to CRC risk is largely unknown. We evaluated the gene-level germline mutation enrichment of 40 DRGs in 680 unselected CRC individuals and 27,728 ancestry-matched cancer-free adults. Significant findings were then examined in independent cohorts of 1,661 unselected CRC individuals and 1,456 individuals with early-onset CRC. Of the 680 individuals in the discovery set, 31 (4.56%) individuals harbored germline pathogenic mutations in known CRC-susceptibility genes, and another 33 (4.85%) individuals had DRG mutations that have not been previously associated with CRC risk. Germline pathogenic mutations in ATM and PALB2 were enriched in both the discovery (OR = 2.81 and p = 0.035 for ATM and OR = 4.91 and p = 0.024 for PALB2) and validation (OR = 2.97 and adjusted p = 0.0013 for ATM and OR = 3.42 and adjusted p = 0.034 for PALB2) sets. Biallelic loss of ATM was evident in all individuals with matched tumor profiling. CRC individuals also had higher rates of actionable mutations in the HR pathway, which can substantially increase the risk of developing cancers other than CRC. Our analysis provides evidence for ATM and PALB2 as CRC-risk genes, underscoring the importance of the homologous recombination pathway in CRC. In addition, we identified frequent complete homologous recombination deficiency in CRC tumors, representing a unique opportunity to explore targeted therapeutic interventions such as poly-ADP ribose polymerase inhibitor (PARPi).

Genomic Correlates of Immune-Cell Infiltrates in Colorectal Carcinoma.

Large-scale genomic characterization of tumors from prospective cohort studies may yield new insights into cancer pathogenesis. We performed whole-exome sequencing of 619 incident colorectal cancers (CRCs) and integrated the results with tumor immunity, pathology, and survival data. We identified recurrently mutated genes in CRC, such as BCL9L, RBM10, CTCF, and KLF5, that were not previously appreciated in this disease. Furthermore, we investigated the genomic correlates of immune-cell infiltration and found that higher neoantigen load was positively associated with overall lymphocytic infiltration, tumor-infiltrating lymphocytes (TILs), memory T cells, and CRC-specific survival. The association with TILs was evident even within microsatellite-stable tumors. We also found positive selection of mutations in HLA genes and other components of the antigen-processing machinery in TIL-rich tumors. These results may inform immunotherapeutic approaches in CRC. More generally, this study demonstrates a framework for future integrative molecular epidemiology research in colorectal and other malignancies.

Clinical, Molecular, and Immune Analysis of Dabrafenib-Trametinib Combination Treatment for BRAF Inhibitor-Refractory Metastatic Melanoma: A Phase 2 Clinical Trial.

IMPORTANCE: Combined treatment with dabrafenib and trametinib (CombiDT) achieves clinical responses in only about 15% of patients with BRAF inhibitor (BRAFi)-refractory metastatic melanoma in contrast to the higher response rate observed in BRAFi-naïve patients. Identifying correlates of response and mechanisms of resistance in this population will facilitate clinical management and rational therapeutic development. OBJECTIVE: To determine correlates of benefit from CombiDT therapy in patients with BRAFi-refractory metastatic melanoma. DESIGN, SETTING, AND PARTICIPANTS: Single-center, single-arm, open-label phase 2 trial of CombiDT treatment in patients with BRAF V600 metastatic melanoma resistant to BRAFi monotherapy conducted between September 2012 and October 2014 at the University of Texas MD Anderson Cancer Center. Key eligibility criteria for participants included BRAF V600 metastatic melanoma, prior BRAFi monotherapy, measurable disease (RECIST 1.1), and tumor accessible for biopsy. INTERVENTIONS: Patients were treated with dabrafenib (150 mg, twice daily) and trametinib (2 mg/d) continuously until disease progression or intolerance. All participants underwent a mandatory baseline biopsy, and optional biopsy specimens were obtained on treatment and at disease progression. Whole-exome sequencing, reverse transcription polymerase chain reaction analysis for BRAF splicing, RNA sequencing, and immunohistochemical analysis were performed on tumor samples, and blood was analyzed for levels of circulating BRAF V600. MAIN OUTCOMES AND MEASURES: The primary end point was overall response rate (ORR). Progression-free survival (PFS) and overall survival (OS) were secondary clinical end points. RESULTS: A total of 28 patients were screened, and 23 enrolled. Among evaluable patients, the confirmed ORR was 10%; disease control rate (DCR) was 45%, and median PFS was 13 weeks. Clinical benefit was associated with duration of prior BRAFi therapy greater than 6 months (DCR, 73% vs 11% for ≤6 months; P = .02) and decrease in circulating BRAF V600 at day 8 of cycle 1 (DCR, 75% vs 18% for no decrease; P = .02) but not with pretreatment mitogen-activated protein kinase (MAPK) pathway mutations or activation. Biopsy specimens obtained during treatment demonstrated that CombiDT therapy failed to achieve significant MAPK pathway inhibition or immune infiltration in most patients. CONCLUSIONS AND RELEVANCE: The baseline presence of MAPK pathway alterations was not associated with benefit from CombiDT in patients with BRAFi-refractory metastatic melanoma. Failure to inhibit the MAPK pathway provides a likely explanation for the limited clinical benefit of CombiDT in this setting. Circulating BRAF V600 is a promising early biomarker of clinical response. TRIAL REGISTRATION: clinicaltrials.gov Identifier: NCT01619774.

An integrated map of structural variation in 2,504 human genomes.

Structural variants are implicated in numerous diseases and make up the majority of varying nucleotides among human genomes. Here we describe an integrated set of eight structural variant classes comprising both balanced and unbalanced variants, which we constructed using short-read DNA sequencing data and statistically phased onto haplotype blocks in 26 human populations. Analysing this set, we identify numerous gene-intersecting structural variants exhibiting population stratification and describe naturally occurring homozygous gene knockouts that suggest the dispensability of a variety of human genes. We demonstrate that structural variants are enriched on haplotypes identified by genome-wide association studies and exhibit enrichment for expression quantitative trait loci. Additionally, we uncover appreciable levels of structural variant complexity at different scales, including genic loci subject to clusters of repeated rearrangement and complex structural variants with multiple breakpoints likely to have formed through individual mutational events. Our catalogue will enhance future studies into structural variant demography, functional impact and disease association.

Analysis of deletion breakpoints from 1,092 humans reveals details of mutation mechanisms.

Investigating genomic structural variants at basepair resolution is crucial for understanding their formation mechanisms. We identify and analyse 8,943 deletion breakpoints in 1,092 samples from the 1000 Genomes Project. We find breakpoints have more nearby SNPs and indels than the genomic average, likely a consequence of relaxed selection. By investigating the correlation of breakpoints with DNA methylation, Hi-C interactions, and histone marks and the substitution patterns of nucleotides near them, we find that breakpoints with the signature of non-allelic homologous recombination (NAHR) are associated with open chromatin. We hypothesize that some NAHR deletions occur without DNA replication and cell division, in embryonic and germline cells. In contrast, breakpoints associated with non-homologous (NH) mechanisms often have sequence microinsertions, templated from later replicating genomic sites, spaced at two characteristic distances from the breakpoint. These microinsertions are consistent with template-switching events and suggest a particular spatiotemporal configuration for DNA during the events.

FunSeq2: a framework for prioritizing noncoding regulatory variants in cancer.

Identification of noncoding drivers from thousands of somatic alterations in a typical tumor is a difficult and unsolved problem. We report a computational framework, FunSeq2, to annotate and prioritize these mutations. The framework combines an adjustable data context integrating large-scale genomics and cancer resources with a streamlined variant-prioritization pipeline. The pipeline has a weighted scoring system combining: inter- and intra-species conservation;loss- and gain-of-function events for transcription-factor binding; enhancer-gene linkages and network centrality; and per-element recurrence across samples. We further highlight putative drivers with information specific to a particular sample, such as differential expression. FunSeq2 is available from funseq2.gersteinlab.org.

Integrative annotation of variants from 1092 humans: application to cancer genomics.

Interpreting variants, especially noncoding ones, in the increasing number of personal genomes is challenging. We used patterns of polymorphisms in functionally annotated regions in 1092 humans to identify deleterious variants; then we experimentally validated candidates. We analyzed both coding and noncoding regions, with the former corroborating the latter. We found regions particularly sensitive to mutations ("ultrasensitive") and variants that are disruptive because of mechanistic effects on transcription-factor binding (that is, "motif-breakers"). We also found variants in regions with higher network centrality tend to be deleterious. Insertions and deletions followed a similar pattern to single-nucleotide variants, with some notable exceptions (e.g., certain deletions and enhancers). On the basis of these patterns, we developed a computational tool (FunSeq), whose application to ~90 cancer genomes reveals nearly a hundred candidate noncoding drivers.