Quantifying the Expanding Landscape of Clinical Actionability for Patients with Cancer.

There is a continuing debate about the proportion of cancer patients that benefit from precision oncology, attributable in part to conflicting views as to which molecular alterations are clinically actionable. To quantify the expansion of clinical actionability since 2017, we annotated 47,271 solid tumors sequenced with the MSK-IMPACT clinical assay using two temporally distinct versions of the OncoKB knowledge base deployed 5 years apart. Between 2017 and 2022, we observed an increase from 8.9% to 31.6% in the fraction of tumors harboring a standard care (level 1 or 2) predictive biomarker of therapy response and an almost halving of tumors carrying nonactionable drivers (44.2% to 22.8%). In tumors with limited or no clinical actionability, TP53 (43.2%), KRAS (19.2%), and CDKN2A (12.2%) were the most frequently altered genes. SIGNIFICANCE: Although clear progress has been made in expanding the availability of precision oncology-based treatment paradigms, our results suggest a continued unmet need for innovative therapeutic strategies, particularly for cancers with currently undruggable oncogenic drivers. See related commentary by Horak and Fröhling, p. 18. This article is featured in Selected Articles from This Issue, p. 5.

Ovarian cancer mutational processes drive site-specific immune evasion.

High-grade serous ovarian cancer (HGSOC) is an archetypal cancer of genomic instability1-4 patterned by distinct mutational processes5,6, tumour heterogeneity7-9 and intraperitoneal spread7,8,10. Immunotherapies have had limited efficacy in HGSOC11-13, highlighting an unmet need to assess how mutational processes and the anatomical sites of tumour foci determine the immunological states of the tumour microenvironment. Here we carried out an integrative analysis of whole-genome sequencing, single-cell RNA sequencing, digital histopathology and multiplexed immunofluorescence of 160 tumour sites from 42 treatment-naive patients with HGSOC. Homologous recombination-deficient HRD-Dup (BRCA1 mutant-like) and HRD-Del (BRCA2 mutant-like) tumours harboured inflammatory signalling and ongoing immunoediting, reflected in loss of HLA diversity and tumour infiltration with highly differentiated dysfunctional CD8+ T cells. By contrast, foldback-inversion-bearing tumours exhibited elevated immunosuppressive TGFß signalling and immune exclusion, with predominantly naive/stem-like and memory T cells. Phenotypic state associations were specific to anatomical sites, highlighting compositional, topological and functional differences between adnexal tumours and distal peritoneal foci. Our findings implicate anatomical sites and mutational processes as determinants of evolutionary phenotypic divergence and immune resistance mechanisms in HGSOC. Our study provides a multi-omic cellular phenotype data substrate from which to develop and interpret future personalized immunotherapeutic approaches and early detection research.

Multimodal integration of radiology, pathology and genomics for prediction of response to PD-(L)1 blockade in patients with non-small cell lung cancer.

Immunotherapy is used to treat almost all patients with advanced non-small cell lung cancer (NSCLC); however, identifying robust predictive biomarkers remains challenging. Here we show the predictive capacity of integrating medical imaging, histopathologic and genomic features to predict immunotherapy response using a cohort of 247 patients with advanced NSCLC with multimodal baseline data obtained during diagnostic clinical workup, including computed tomography scan images, digitized programmed death ligand-1 immunohistochemistry slides and known outcomes to immunotherapy. Using domain expert annotations, we developed a computational workflow to extract patient-level features and used a machine-learning approach to integrate multimodal features into a risk prediction model. Our multimodal model (area under the curve (AUC) = 0.80, 95% confidence interval (CI) 0.74-0.86) outperformed unimodal measures, including tumor mutational burden (AUC = 0.61, 95% CI 0.52-0.70) and programmed death ligand-1 immunohistochemistry score (AUC = 0.73, 95% CI 0.65-0.81). Our study therefore provides a quantitative rationale for using multimodal features to improve prediction of immunotherapy response in patients with NSCLC using expert-guided machine learning.

Multimodal data integration using machine learning improves risk stratification of high-grade serous ovarian cancer.

Patients with high-grade serous ovarian cancer suffer poor prognosis and variable response to treatment. Known prognostic factors for this disease include homologous recombination deficiency status, age, pathological stage and residual disease status after debulking surgery. Recent work has highlighted important prognostic information captured in computed tomography and histopathological specimens, which can be exploited through machine learning. However, little is known about the capacity of combining features from these disparate sources to improve prediction of treatment response. Here, we assembled a multimodal dataset of 444 patients with primarily late-stage high-grade serous ovarian cancer and discovered quantitative features, such as tumor nuclear size on staining with hematoxylin and eosin and omental texture on contrast-enhanced computed tomography, associated with prognosis. We found that these features contributed complementary prognostic information relative to one another and clinicogenomic features. By fusing histopathological, radiologic and clinicogenomic machine-learning models, we demonstrate a promising path toward improved risk stratification of patients with cancer through multimodal data integration.

Genome Nexus: A Comprehensive Resource for the Annotation and Interpretation of Genomic Variants in Cancer.

PURPOSE: Interpretation of genomic variants in tumor samples still presents a challenge in research and the clinical setting. A major issue is that information for variant interpretation is fragmented across disparate databases, and aggregation of information from these requires building extensive infrastructure. To this end, we have developed Genome Nexus, a one-stop shop for variant annotation with a user-friendly interface for cancer researchers and clinicians. METHODS: Genome Nexus (1) aggregates variant information from sources that are relevant to cancer research and clinical applications, (2) allows high-performance programmatic access to the aggregated data via a unified application programming interface, (3) provides a reference page for individual cancer variants, (4) provides user-friendly tools for annotating variants in patients, and (5) is freely available under an open source license and can be installed in a private cloud or local environment and integrated with local institutional resources. RESULTS: Genome Nexus is available at https://www.genomenexus.org. It displays annotations from more than a dozen resources including those that provide variant effect information (variant effect predictor), protein sequence annotation (Uniprot, Pfam, and dbPTM), functional consequence prediction (Polyphen-2, Mutation Assessor, and SIFT), population prevalences (gnomAD, dbSNP, and ExAC), cancer population prevalences (Cancer hotspots and SignalDB), and clinical actionability (OncoKB, CIViC, and ClinVar). We describe several use cases that demonstrate the utility of Genome Nexus to clinicians, researchers, and bioinformaticians. We cover single-variant annotation, cohort analysis, and programmatic use of the application programming interface. Genome Nexus is unique in providing a user-friendly interface specific to cancer that allows high-performance annotation of any variant including unknown ones. CONCLUSION: Interpretation of cancer genomic variants is improved tremendously by having an integrated resource for annotations. Genome Nexus is freely available under an open source license.

Semisupervised Training of a Brain MRI Tumor Detection Model Using Mined Annotations.

Background Artificial intelligence (AI) applications for cancer imaging conceptually begin with automated tumor detection, which can provide the foundation for downstream AI tasks. However, supervised training requires many image annotations, and performing dedicated post hoc image labeling is burdensome and costly. Purpose To investigate whether clinically generated image annotations can be data mined from the picture archiving and communication system (PACS), automatically curated, and used for semisupervised training of a brain MRI tumor detection model. Materials and Methods In this retrospective study, the cancer center PACS was mined for brain MRI scans acquired between January 2012 and December 2017 and included all annotated axial T1 postcontrast images. Line annotations were converted to boxes, excluding boxes shorter than 1 cm or longer than 7 cm. The resulting boxes were used for supervised training of object detection models using RetinaNet and Mask region-based convolutional neural network (R-CNN) architectures. The best-performing model trained from the mined data set was used to detect unannotated tumors on training images themselves (self-labeling), automatically correcting many of the missing labels. After self-labeling, new models were trained using this expanded data set. Models were scored for precision, recall, and F1 using a held-out test data set comprising 754 manually labeled images from 100 patients (403 intra-axial and 56 extra-axial enhancing tumors). Model F1 scores were compared using bootstrap resampling. Results The PACS query extracted 31 150 line annotations, yielding 11 880 boxes that met inclusion criteria. This mined data set was used to train models, yielding F1 scores of 0.886 for RetinaNet and 0.908 for Mask R-CNN. Self-labeling added 18 562 training boxes, improving model F1 scores to 0.935 (P < .001) and 0.954 (P < .001), respectively. Conclusion The application of semisupervised learning to mined image annotations significantly improved tumor detection performance, achieving an excellent F1 score of 0.954. This development pipeline can be extended for other imaging modalities, repurposing unused data silos to potentially enable automated tumor detection across radiologic modalities. © RSNA, 2022 Online supplemental material is available for this article.

Harnessing multimodal data integration to advance precision oncology.

Advances in quantitative biomarker development have accelerated new forms of data-driven insights for patients with cancer. However, most approaches are limited to a single mode of data, leaving integrated approaches across modalities relatively underdeveloped. Multimodal integration of advanced molecular diagnostics, radiological and histological imaging, and codified clinical data presents opportunities to advance precision oncology beyond genomics and standard molecular techniques. However, most medical datasets are still too sparse to be useful for the training of modern machine learning techniques, and significant challenges remain before this is remedied. Combined efforts of data engineering, computational methods for analysis of heterogeneous data and instantiation of synergistic data models in biomedical research are required for success. In this Perspective, we offer our opinions on synthesizing complementary modalities of data with emerging multimodal artificial intelligence methods. Advancing along this direction will result in a reimagined class of multimodal biomarkers to propel the field of precision oncology in the coming decade.

The context-specific role of germline pathogenicity in tumorigenesis.

Human cancers arise from environmental, heritable and somatic factors, but how these mechanisms interact in tumorigenesis is poorly understood. Studying 17,152 prospectively sequenced patients with cancer, we identified pathogenic germline variants in cancer predisposition genes, and assessed their zygosity and co-occurring somatic alterations in the concomitant tumors. Two major routes to tumorigenesis were apparent. In carriers of pathogenic germline variants in high-penetrance genes (5.1% overall), lineage-dependent patterns of biallelic inactivation led to tumors exhibiting mechanism-specific somatic phenotypes and fewer additional somatic oncogenic drivers. Nevertheless, 27% of cancers in these patients, and most tumors in patients with pathogenic germline variants in lower-penetrance genes, lacked particular hallmarks of tumorigenesis associated with the germline allele. The dependence of tumors on pathogenic germline variants is variable and often dictated by both penetrance and lineage, a finding with implications for clinical management.

Integrated digital pathology at scale: A solution for clinical diagnostics and cancer research at a large academic medical center.

OBJECTIVE: Broad adoption of digital pathology (DP) is still lacking, and examples for DP connecting diagnostic, research, and educational use cases are missing. We blueprint a holistic DP solution at a large academic medical center ubiquitously integrated into clinical workflows; researchapplications including molecular, genetic, and tissue databases; and educational processes. MATERIALS AND METHODS: We built a vendor-agnostic, integrated viewer for reviewing, annotating, sharing, and quality assurance of digital slides in a clinical or research context. It is the first homegrown viewer cleared by New York State provisional approval in 2020 for primary diagnosis and remote sign-out during the COVID-19 (coronavirus disease 2019) pandemic. We further introduce an interconnected Honest Broker for BioInformatics Technology (HoBBIT) to systematically compile and share large-scale DP research datasets including anonymized images, redacted pathology reports, and clinical data of patients with consent. RESULTS: The solution has been operationally used over 3 years by 926 pathologists and researchers evaluating 288 903 digital slides. A total of 51% of these were reviewed within 1 month after scanning. Seamless integration of the viewer into 4 hospital systems clearly increases the adoption of DP. HoBBIT directly impacts the translation of knowledge in pathology into effective new health measures, including artificial intelligence-driven detection models for prostate cancer, basal cell carcinoma, and breast cancer metastases, developed and validated on thousands of cases. CONCLUSIONS: We highlight major challenges and lessons learned when going digital to provide orientation for other pathologists. Building interconnected solutions will not only increase adoption of DP, but also facilitate next-generation computational pathology at scale for enhanced cancer research.

Genomic profiling identifies somatic mutations predicting thromboembolic risk in patients with solid tumors.

Venous thromboembolism (VTE) associated with cancer (CAT) is a well-described complication of cancer and a leading cause of death in patients with cancer. The purpose of this study was to assess potential associations of molecular signatures with CAT, including tumor-specific mutations and the presence of clonal hematopoiesis. We analyzed deep-coverage targeted DNA-sequencing data of >14 000 solid tumor samples using the Memorial Sloan Kettering-Integrated Mutation Profiling of Actionable Cancer Targets platform to identify somatic alterations associated with VTE. End point was defined as the first instance of cancer-associated pulmonary embolism and/or proximal/distal lower extremity deep vein thrombosis. Cause-specific Cox proportional hazards regression was used, adjusting for pertinent clinical covariates. Of 11 695 evaluable individuals, 72% had metastatic disease at time of analysis. Tumor-specific mutations in KRAS (hazard ratio [HR], 1.34; 95% confidence interval (CI), 1.09-1.64; adjusted P = .08), STK11 (HR, 2.12; 95% CI, 1.55-2.89; adjusted P < .001), KEAP1 (HR, 1.84; 95% CI, 1.21-2.79; adjusted P = .07), CTNNB1 (HR, 1.73; 95% CI, 1.15-2.60; adjusted P = .09), CDKN2B (HR, 1.45; 95% CI, 1.13-1.85; adjusted P = .07), and MET (HR, 1.83; 95% CI, 1.15-2.92; adjusted P = .09) were associated with a significantly increased risk of CAT independent of tumor type. Mutations in SETD2 were associated with a decreased risk of CAT (HR, 0.35; 95% CI, 0.16-0.79; adjusted P = .09). The presence of clonal hematopoiesis was not associated with an increased VTE rate. This is the first large-scale analysis to elucidate tumor-specific genomic events associated with CAT. Somatic tumor mutations of STK11, KRAS, CTNNB1, KEAP1, CDKN2B, and MET were associated with an increased risk of VTE in patients with solid tumors. Further analysis is needed to validate these findings and identify additional molecular signatures unique to individual tumor types.

Multiomic Integration of Public Oncology Databases in Bioconductor.

PURPOSE: Investigations of the molecular basis for the development, progression, and treatment of cancer increasingly use complementary genomic assays to gather multiomic data, but management and analysis of such data remain complex. The cBioPortal for cancer genomics currently provides multiomic data from > 260 public studies, including The Cancer Genome Atlas (TCGA) data sets, but integration of different data types remains challenging and error prone for computational methods and tools using these resources. Recent advances in data infrastructure within the Bioconductor project enable a novel and powerful approach to creating fully integrated representations of these multiomic, pan-cancer databases. METHODS: We provide a set of R/Bioconductor packages for working with TCGA legacy data and cBioPortal data, with special considerations for loading time; efficient representations in and out of memory; analysis platform; and an integrative framework, such as MultiAssayExperiment. Large methylation data sets are provided through out-of-memory data representation to provide responsive loading times and analysis capabilities on machines with limited memory. RESULTS: We developed the curatedTCGAData and cBioPortalData R/Bioconductor packages to provide integrated multiomic data sets from the TCGA legacy database and the cBioPortal web application programming interface using the MultiAssayExperiment data structure. This suite of tools provides coordination of diverse experimental assays with clinicopathological data with minimal data management burden, as demonstrated through several greatly simplified multiomic and pan-cancer analyses. CONCLUSION: These integrated representations enable analysts and tool developers to apply general statistical and plotting methods to extensive multiomic data through user-friendly commands and documented examples.

Leveraging Systematic Functional Analysis to Benchmark an In Silico Framework Distinguishes Driver from Passenger MEK Mutants in Cancer.

Despite significant advances in cancer precision medicine, a significant hurdle to its broader adoption remains the multitude of variants of unknown significance identified by clinical tumor sequencing and the lack of biologically validated methods to distinguish between functional and benign variants. Here we used functional data on MAP2K1 and MAP2K2 mutations generated in real-time within a co-clinical trial framework to benchmark the predictive value of a three-part in silico methodology. Our computational approach to variant classification incorporated hotspot analysis, three-dimensional molecular dynamics simulation, and sequence paralogy. In silico prediction accurately distinguished functional from benign MAP2K1 and MAP2K2 mutants, yet drug sensitivity varied widely among activating mutant alleles. These results suggest that multifaceted in silico modeling can inform patient accrual to MEK/ERK inhibitor clinical trials, but computational methods need to be paired with laboratory- and clinic-based efforts designed to unravel variabilities in drug response. SIGNIFICANCE: Leveraging prospective functional characterization of MEK1/2 mutants, it was found that hotspot analysis, molecular dynamics simulation, and sequence paralogy are complementary tools that can robustly prioritize variants for biologic, therapeutic, and clinical validation.See related commentary by Whitehead and Sebolt-Leopold, p. 4042.

Detecting Cancer Survival Related Gene Markers Based on Rectified Factor Network.

Detecting gene sets that serve as biomarkers for differentiating patient survival groups may help diagnose diseases robustly and develop multi-gene targeted therapies. However, due to the exponential growth of search space imposed by gene combinations, the performance of existing methods is still far from satisfactory. In this study, we developed a new method called BISG (BIclustering based Survival-related Gene sets detection) based on a rectified factor network (RFN) model, which allows efficiently biclustering gene subsets. By correlating genes in each significant bicluster with patient survival outcomes using a log-rank test and multi-sampling strategy, multiple survival-related gene sets can be detected. We applied BISG on three different cancer types, and the resulting gene sets were tested as biomarkers for survival analyses. Secondly, we systematically analyzed 12 different cancer datasets. Our analysis shows that the genes in all the survival-related gene sets are mainly from five gene families: microRNA protein coding host genes, zinc fingers C2H2-type, solute carriers, CD (cluster of differentiation) molecules, and ankyrin repeat domain containing genes. Moreover, we found that they are mainly enriched in heme metabolism, apoptosis, hypoxia and inflammatory response-related pathways. We compared BISG with two other methods, GSAS and IPSOV. Results show that BISG can better differentiate patient survival groups in different datasets. The identified biomarkers suggested by our study provide useful hypotheses for further investigation. BISG is publicly available with open source at https://github.com/LingtaoSu/BISG.

Modeling biological and genetic diversity in upper tract urothelial carcinoma with patient derived xenografts.

Treatment paradigms for patients with upper tract urothelial carcinoma (UTUC) are typically extrapolated from studies of bladder cancer despite their distinct clinical and molecular characteristics. The advancement of UTUC research is hampered by the lack of disease-specific models. Here, we report the establishment of patient derived xenograft (PDX) and cell line models that reflect the genomic and biological heterogeneity of the human disease. Models demonstrate high genomic concordance with the corresponding patient tumors, with invasive tumors more likely to successfully engraft. Treatment of PDX models with chemotherapy recapitulates responses observed in patients. Analysis of a HER2 S310F-mutant PDX suggests that an antibody drug conjugate targeting HER2 would have superior efficacy versus selective HER2 kinase inhibitors. In sum, the biological and phenotypic concordance between patient and PDXs suggest that these models could facilitate studies of intrinsic and acquired resistance and the development of personalized medicine strategies for UTUC patients.

A harmonized meta-knowledgebase of clinical interpretations of somatic genomic variants in cancer.

Precision oncology relies on accurate discovery and interpretation of genomic variants, enabling individualized diagnosis, prognosis and therapy selection. We found that six prominent somatic cancer variant knowledgebases were highly disparate in content, structure and supporting primary literature, impeding consensus when evaluating variants and their relevance in a clinical setting. We developed a framework for harmonizing variant interpretations to produce a meta-knowledgebase of 12,856 aggregate interpretations. We demonstrated large gains in overlap between resources across variants, diseases and drugs as a result of this harmonization. We subsequently demonstrated improved matching between a patient cohort and harmonized interpretations of potential clinical significance, observing an increase from an average of 33% per individual knowledgebase to 57% in aggregate. Our analyses illuminate the need for open, interoperable sharing of variant interpretation data. We also provide a freely available web interface (search.cancervariants.org) for exploring the harmonized interpretations from these six knowledgebases.

Characteristics and Outcome of AKT1 E17K-Mutant Breast Cancer Defined through AACR Project GENIE, a Clinicogenomic Registry.

AKT inhibitors have promising activity in AKT1 E17K-mutant estrogen receptor (ER)-positive metastatic breast cancer, but the natural history of this rare genomic subtype remains unknown. Utilizing AACR Project GENIE, an international clinicogenomic data-sharing consortium, we conducted a comparative analysis of clinical outcomes of patients with matched AKT1 E17K-mutant (n = 153) and AKT1-wild-type (n = 302) metastatic breast cancer. AKT1-mutant cases had similar adjusted overall survival (OS) compared with AKT1-wild-type controls (median OS, 24.1 vs. 29.9, respectively; P = 0.98). AKT1-mutant cases enjoyed longer durations on mTOR inhibitor therapy, an observation previously unrecognized in pivotal clinical trials due to the rarity of this alteration. Other baseline clinicopathologic features, as well as durations on other classes of therapy, were broadly similar. In summary, we demonstrate the feasibility of using a novel and publicly accessible clincogenomic registry to define outcomes in a rare genomically defined cancer subtype, an approach with broad applicability to precision oncology. SIGNIFICANCE: We delineate the natural history of a rare genomically distinct cancer, AKT1 E17K-mutant ER-positive breast cancer, using a publicly accessible registry of real-world patient data, thereby illustrating the potential to inform drug registration through synthetic control data.See related commentary by Castellanos and Baxi, p. 490.

How oncogenic mutations activate human MAP kinase 1 (MEK1): a molecular dynamics simulation study.

Approximately 30% of all types of human cancers possess a constitutively activated the mitogen-activated protein kinase (MAPK) signaling pathway while MAP kinase 1 (MEK1) is a critical component of this pathway. It has been reported mutations could improve the activity of MEK1 to result in cell proliferation and transformation, which is a known oncogenic event in various cancer types. In this study, eight molecular dynamics simulations, molecular mechanics Poisson-Boltzmann surface area (MM-PBSA), combined with protein structure network were performed to explore the mechanism that mutations activate MEK1. Protein structure networks and hydrogen bonds analysis demonstrated that active mutations broke the interaction between activation segments (residues 216-222) and C-helix (residues 105-121) in MEK1, leading to it transform inactive form to active form. Moreover, hydrogen bond analysis and MM-PBSA calculation indicated that activating mutations decrease the binding affinity between MEK1 and inhibitor to reduce the inhibitory effect of inhibitors. In addition, some active mutations cause structural changes in the Pro-rich loop (residues 261-268) of MEK1. These changes may stabilize the interaction between the MEK1 mutants and the ligands by increasing the number of exposed hydrophobic residues in the active site of MEK1. Our results may provide useful theoretical evidences for the mechanism underlying the role of human MEK1 in human cancers.Communicated by Ramaswamy H. Sarma.

Integration and Analysis of CPTAC Proteomics Data in the Context of Cancer Genomics in the cBioPortal.

The Clinical Proteomic Tumor Analysis Consortium (CPTAC) has produced extensive mass spectrometry-based proteomics data for selected breast, colon, and ovarian tumors from The Cancer Genome Atlas (TCGA). We have incorporated the CPTAC proteomics data into the cBioPortal to support easy exploration and integrative analysis of these proteomic datasets in the context of the clinical and genomics data from the same tumors. cBioPortal is an open source platform for exploring, visualizing, and analyzing multidimensional cancer genomics and clinical data. The public instance of the cBioPortal (http://cbioportal.org/) hosts more than 200 cancer genomics studies, including all of the data from TCGA. Its biologist-friendly interface provides many rich analysis features, including a graphical summary of gene-level data across multiple platforms, correlation analysis between genes or other data types, survival analysis, and per-patient data visualization. Here, we present the integration of the CPTAC mass spectrometry-based proteomics data into the cBioPortal, consisting of 77 breast, 95 colorectal, and 174 ovarian tumors that already have been profiled by TCGA for mutations, copy number alterations, gene expression, and DNA methylation. As a result, the CPTAC data can now be easily explored and analyzed in the cBioPortal in the context of clinical and genomics data. By integrating CPTAC data into cBioPortal, limitations of TCGA proteomics array data can be overcome while also providing a user-friendly web interface, a web API, and an R client to query the mass spectrometry data together with genomic, epigenomic, and clinical data.

Publisher Correction: The long tail of oncogenic drivers in prostate cancer.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The expanding landscape of 'oncohistone' mutations in human cancers.

Mutations in epigenetic pathways are common oncogenic drivers. Histones, the fundamental substrates for chromatin-modifying and remodelling enzymes, are mutated in tumours including gliomas, sarcomas, head and neck cancers, and carcinosarcomas. Classical 'oncohistone' mutations occur in the N-terminal tail of histone H3 and affect the function of polycomb repressor complexes 1 and 2 (PRC1 and PRC2). However, the prevalence and function of histone mutations in other tumour contexts is unknown. Here we show that somatic histone mutations occur in approximately 4% (at a conservative estimate) of diverse tumour types and in crucial regions of histone proteins. Mutations occur in all four core histones, in both the N-terminal tails and globular histone fold domains, and at or near residues that contain important post-translational modifications. Many globular domain mutations are homologous to yeast mutants that abrogate the need for SWI/SNF function, occur in the key regulatory 'acidic patch' of histones H2A and H2B, or are predicted to disrupt the H2B-H4 interface. The histone mutation dataset and the hypotheses presented here on the effect of the mutations on important chromatin functions should serve as a resource and starting point for the chromatin and cancer biology fields in exploring an expanding role of histone mutations in cancer.