The mutational constraint spectrum quantified from variation in 141,456 humans.

Genetic variants that inactivate protein-coding genes are a powerful source of information about the phenotypic consequences of gene disruption: genes that are crucial for the function of an organism will be depleted of such variants in natural populations, whereas non-essential genes will tolerate their accumulation. However, predicted loss-of-function variants are enriched for annotation errors, and tend to be found at extremely low frequencies, so their analysis requires careful variant annotation and very large sample sizes1. Here we describe the aggregation of 125,748 exomes and 15,708 genomes from human sequencing studies into the Genome Aggregation Database (gnomAD). We identify 443,769 high-confidence predicted loss-of-function variants in this cohort after filtering for artefacts caused by sequencing and annotation errors. Using an improved model of human mutation rates, we classify human protein-coding genes along a spectrum that represents tolerance to inactivation, validate this classification using data from model organisms and engineered human cells, and show that it can be used to improve the power of gene discovery for both common and rare diseases.

Systematic genomic and translational efficiency studies of uveal melanoma.

To further our understanding of the somatic genetic basis of uveal melanoma, we sequenced the protein-coding regions of 52 primary tumors and 3 liver metastases together with paired normal DNA. Known recurrent mutations were identified in GNAQ, GNA11, BAP1, EIF1AX, and SF3B1. The role of mutated EIF1AX was tested using loss of function approaches including viability and translational efficiency assays. Knockdown of both wild type and mutant EIF1AX was lethal to uveal melanoma cells. We probed the function of N-terminal tail EIF1AX mutations by performing RNA sequencing of polysome-associated transcripts in cells expressing endogenous wild type or mutant EIF1AX. Ribosome occupancy of the global translational apparatus was sensitive to suppression of wild type but not mutant EIF1AX. Together, these studies suggest that cells expressing mutant EIF1AX may exhibit aberrant translational regulation, which may provide clonal selective advantage in the subset of uveal melanoma that harbors this mutation.

Clonal evolution in patients with chronic lymphocytic leukaemia developing resistance to BTK inhibition.

Resistance to the Bruton's tyrosine kinase (BTK) inhibitor ibrutinib has been attributed solely to mutations in BTK and related pathway molecules. Using whole-exome and deep-targeted sequencing, we dissect evolution of ibrutinib resistance in serial samples from five chronic lymphocytic leukaemia patients. In two patients, we detect BTK-C481S mutation or multiple PLCG2 mutations. The other three patients exhibit an expansion of clones harbouring del(8p) with additional driver mutations (EP300, MLL2 and EIF2A), with one patient developing trans-differentiation into CD19-negative histiocytic sarcoma. Using droplet-microfluidic technology and growth kinetic analyses, we demonstrate the presence of ibrutinib-resistant subclones and estimate subclone size before treatment initiation. Haploinsufficiency of TRAIL-R, a consequence of del(8p), results in TRAIL insensitivity, which may contribute to ibrutinib resistance. These findings demonstrate that the ibrutinib therapy favours selection and expansion of rare subclones already present before ibrutinib treatment, and provide insight into the heterogeneity of genetic changes associated with ibrutinib resistance.

Mutations driving CLL and their evolution in progression and relapse.

Which genetic alterations drive tumorigenesis and how they evolve over the course of disease and therapy are central questions in cancer biology. Here we identify 44 recurrently mutated genes and 11 recurrent somatic copy number variations through whole-exome sequencing of 538 chronic lymphocytic leukaemia (CLL) and matched germline DNA samples, 278 of which were collected in a prospective clinical trial. These include previously unrecognized putative cancer drivers (RPS15, IKZF3), and collectively identify RNA processing and export, MYC activity, and MAPK signalling as central pathways involved in CLL. Clonality analysis of this large data set further enabled reconstruction of temporal relationships between driver events. Direct comparison between matched pre-treatment and relapse samples from 59 patients demonstrated highly frequent clonal evolution. Thus, large sequencing data sets of clinically informative samples enable the discovery of novel genes associated with cancer, the network of relationships between the driver events, and their impact on disease relapse and clinical outcome.

Comprehensive analysis of cancer-associated somatic mutations in class I HLA genes.

Detection of somatic mutations in human leukocyte antigen (HLA) genes using whole-exome sequencing (WES) is hampered by the high polymorphism of the HLA loci, which prevents alignment of sequencing reads to the human reference genome. We describe a computational pipeline that enables accurate inference of germline alleles of class I HLA-A, B and C genes and subsequent detection of mutations in these genes using the inferred alleles as a reference. Analysis of WES data from 7,930 pairs of tumor and healthy tissue from the same patient revealed 298 nonsilent HLA mutations in tumors from 266 patients. These 298 mutations are enriched for likely functional mutations, including putative loss-of-function events. Recurrence of mutations suggested that these 'hotspot' sites were positively selected. Cancers with recurrent somatic HLA mutations were associated with upregulation of signatures of cytolytic activity characteristic of tumor infiltration by effector lymphocytes, supporting immune evasion by altered HLA function as a contributory mechanism in cancer.

Genomic Characterization of Brain Metastases Reveals Branched Evolution and Potential Therapeutic Targets.

UNLABELLED: Brain metastases are associated with a dismal prognosis. Whether brain metastases harbor distinct genetic alterations beyond those observed in primary tumors is unknown. We performed whole-exome sequencing of 86 matched brain metastases, primary tumors, and normal tissue. In all clonally related cancer samples, we observed branched evolution, where all metastatic and primary sites shared a common ancestor yet continued to evolve independently. In 53% of cases, we found potentially clinically informative alterations in the brain metastases not detected in the matched primary-tumor sample. In contrast, spatially and temporally separated brain metastasis sites were genetically homogenous. Distal extracranial and regional lymph node metastases were highly divergent from brain metastases. We detected alterations associated with sensitivity to PI3K/AKT/mTOR, CDK, and HER2/EGFR inhibitors in the brain metastases. Genomic analysis of brain metastases provides an opportunity to identify potentially clinically informative alterations not detected in clinically sampled primary tumors, regional lymph nodes, or extracranial metastases. SIGNIFICANCE: Decisions for individualized therapies in patients with brain metastasis are often made from primary-tumor biopsies. We demonstrate that clinically actionable alterations present in brain metastases are frequently not detected in primary biopsies, suggesting that sequencing of primary biopsies alone may miss a substantial number of opportunities for targeted therapy.

Prospective derivation of a living organoid biobank of colorectal cancer patients.

In Rspondin-based 3D cultures, Lgr5 stem cells from multiple organs form ever-expanding epithelial organoids that retain their tissue identity. We report the establishment of tumor organoid cultures from 20 consecutive colorectal carcinoma (CRC) patients. For most, organoids were also generated from adjacent normal tissue. Organoids closely recapitulate several properties of the original tumor. The spectrum of genetic changes within the "living biobank" agrees well with previous large-scale mutational analyses of CRC. Gene expression analysis indicates that the major CRC molecular subtypes are represented. Tumor organoids are amenable to high-throughput drug screens allowing detection of gene-drug associations. As an example, a single organoid culture was exquisitely sensitive to Wnt secretion (porcupine) inhibitors and carried a mutation in the negative Wnt feedback regulator RNF43, rather than in APC. Organoid technology may fill the gap between cancer genetics and patient trials, complement cell-line- and xenograft-based drug studies, and allow personalized therapy design. PAPERCLIP.

RNF43 is frequently mutated in colorectal and endometrial cancers.

We report somatic mutations of RNF43 in over 18% of colorectal adenocarcinomas and endometrial carcinomas. RNF43 encodes an E3 ubiquitin ligase that negatively regulates Wnt signaling. Truncating mutations of RNF43 are more prevalent in microsatellite-unstable tumors and show mutual exclusivity with inactivating APC mutations in colorectal adenocarcinomas. These results indicate that RNF43 is one of the most commonly mutated genes in colorectal and endometrial cancers.

The genomic landscape of pediatric Ewing sarcoma.

UNLABELLED: Pediatric Ewing sarcoma is characterized by the expression of chimeric fusions of EWS and ETS family transcription factors, representing a paradigm for studying cancers driven by transcription factor rearrangements. In this study, we describe the somatic landscape of pediatric Ewing sarcoma. These tumors are among the most genetically normal cancers characterized to date, with only EWS-ETS rearrangements identified in the majority of tumors. STAG2 loss, however, is present in more than 15% of Ewing sarcoma tumors; occurs by point mutation, rearrangement, and likely nongenetic mechanisms; and is associated with disease dissemination. Perhaps the most striking finding is the paucity of mutations in immediately targetable signal transduction pathways, highlighting the need for new therapeutic approaches to target EWS-ETS fusions in this disease. SIGNIFICANCE: We performed next-generation sequencing of Ewing sarcoma, a pediatric cancer involving bone, characterized by expression of EWS-ETS fusions. We found remarkably few mutations. However, we discovered that loss of STAG2 expression occurs in 15% of tumors and is associated with metastatic disease, suggesting a potential genetic vulnerability in Ewing sarcoma.

Systematic identification of personal tumor-specific neoantigens in chronic lymphocytic leukemia.

Genome sequencing has revealed a large number of shared and personal somatic mutations across human cancers. In principle, any genetic alteration affecting a protein-coding region has the potential to generate mutated peptides that are presented by surface HLA class I proteins that might be recognized by cytotoxic T cells. To test this possibility, we implemented a streamlined approach for the prediction and validation of such neoantigens derived from individual tumors and presented by patient-specific HLA alleles. We applied our computational pipeline to 91 chronic lymphocytic leukemias (CLLs) that underwent whole-exome sequencing (WES). We predicted ∼22 mutated HLA-binding peptides per leukemia (derived from ∼16 missense mutations) and experimentally confirmed HLA binding for ∼55% of such peptides. Two CLL patients that achieved long-term remission following allogeneic hematopoietic stem cell transplantation were monitored for CD8(+) T-cell responses against predicted or confirmed HLA-binding peptides. Long-lived cytotoxic T-cell responses were detected against peptides generated from personal tumor mutations in ALMS1, C6ORF89, and FNDC3B presented on tumor cells. Finally, we applied our computational pipeline to WES data (N = 2488 samples) across 13 different cancer types and estimated dozens to thousands of predicted neoantigens per individual tumor, suggesting that neoantigens are frequent in most tumors.

Colon cancer-derived oncogenic EGFR G724S mutant identified by whole genome sequence analysis is dependent on asymmetric dimerization and sensitive to cetuximab.

BACKGROUND: Inhibition of the activated epidermal growth factor receptor (EGFR) with either enzymatic kinase inhibitors or anti-EGFR antibodies such as cetuximab, is an effective modality of treatment for multiple human cancers. Enzymatic EGFR inhibitors are effective for lung adenocarcinomas with somatic kinase domain EGFR mutations while, paradoxically, anti-EGFR antibodies are more effective in colon and head and neck cancers where EGFR mutations occur less frequently. In colorectal cancer, anti-EGFR antibodies are routinely used as second-line therapy of KRAS wild-type tumors. However, detailed mechanisms and genomic predictors for pharmacological response to these antibodies in colon cancer remain unclear. FINDINGS: We describe a case of colorectal adenocarcinoma, which was found to harbor a kinase domain mutation, G724S, in EGFR through whole genome sequencing. We show that G724S mutant EGFR is oncogenic and that it differs from classic lung cancer derived EGFR mutants in that it is cetuximab responsive in vitro, yet relatively insensitive to small molecule kinase inhibitors. Through biochemical and cellular pharmacologic studies, we have determined that cells harboring the colon cancer-derived G719S and G724S mutants are responsive to cetuximab therapy in vitro and found that the requirement for asymmetric dimerization of these mutant EGFR to promote cellular transformation may explain their greater inhibition by cetuximab than small-molecule kinase inhibitors. CONCLUSION: The colon-cancer derived G719S and G724S mutants are oncogenic and sensitive in vitro to cetuximab. These data suggest that patients with these mutations may benefit from the use of anti-EGFR antibodies as part of the first-line therapy.

Whole-exome sequencing and clinical interpretation of formalin-fixed, paraffin-embedded tumor samples to guide precision cancer medicine.

Translating whole-exome sequencing (WES) for prospective clinical use may have an impact on the care of patients with cancer; however, multiple innovations are necessary for clinical implementation. These include rapid and robust WES of DNA derived from formalin-fixed, paraffin-embedded tumor tissue, analytical output similar to data from frozen samples and clinical interpretation of WES data for prospective use. Here, we describe a prospective clinical WES platform for archival formalin-fixed, paraffin-embedded tumor samples. The platform employs computational methods for effective clinical analysis and interpretation of WES data. When applied retrospectively to 511 exomes, the interpretative framework revealed a 'long tail' of somatic alterations in clinically important genes. Prospective application of this approach identified clinically relevant alterations in 15 out of 16 patients. In one patient, previously undetected findings guided clinical trial enrollment, leading to an objective clinical response. Overall, this methodology may inform the widespread implementation of precision cancer medicine.

Whole-exome sequencing of circulating tumor cells provides a window into metastatic prostate cancer.

Comprehensive analyses of cancer genomes promise to inform prognoses and precise cancer treatments. A major barrier, however, is inaccessibility of metastatic tissue. A potential solution is to characterize circulating tumor cells (CTCs), but this requires overcoming the challenges of isolating rare cells and sequencing low-input material. Here we report an integrated process to isolate, qualify and sequence whole exomes of CTCs with high fidelity using a census-based sequencing strategy. Power calculations suggest that mapping of >99.995% of the standard exome is possible in CTCs. We validated our process in two patients with prostate cancer, including one for whom we sequenced CTCs, a lymph node metastasis and nine cores of the primary tumor. Fifty-one of 73 CTC mutations (70%) were present in matched tissue. Moreover, we identified 10 early trunk and 56 metastatic trunk mutations in the non-CTC tumor samples and found 90% and 73% of these mutations, respectively, in CTC exomes. This study establishes a foundation for CTC genomics in the clinic.

Genetic and clonal dissection of murine small cell lung carcinoma progression by genome sequencing.

Small cell lung carcinoma (SCLC) is a highly lethal, smoking-associated cancer with few known targetable genetic alterations. Using genome sequencing, we characterized the somatic evolution of a genetically engineered mouse model (GEMM) of SCLC initiated by loss of Trp53 and Rb1. We identified alterations in DNA copy number and complex genomic rearrangements and demonstrated a low somatic point mutation frequency in the absence of tobacco mutagens. Alterations targeting the tumor suppressor Pten occurred in the majority of murine SCLC studied, and engineered Pten deletion accelerated murine SCLC and abrogated loss of Chr19 in Trp53; Rb1; Pten compound mutant tumors. Finally, we found evidence for polyclonal and sequential metastatic spread of murine SCLC by comparative sequencing of families of related primary tumors and metastases. We propose a temporal model of SCLC tumorigenesis with implications for human SCLC therapeutics and the nature of cancer-genome evolution in GEMMs.

Landscape of genomic alterations in cervical carcinomas.

Cervical cancer is responsible for 10-15% of cancer-related deaths in women worldwide. The aetiological role of infection with high-risk human papilloma viruses (HPVs) in cervical carcinomas is well established. Previous studies have also implicated somatic mutations in PIK3CA, PTEN, TP53, STK11 and KRAS as well as several copy-number alterations in the pathogenesis of cervical carcinomas. Here we report whole-exome sequencing analysis of 115 cervical carcinoma-normal paired samples, transcriptome sequencing of 79 cases and whole-genome sequencing of 14 tumour-normal pairs. Previously unknown somatic mutations in 79 primary squamous cell carcinomas include recurrent E322K substitutions in the MAPK1 gene (8%), inactivating mutations in the HLA-B gene (9%), and mutations in EP300 (16%), FBXW7 (15%), NFE2L2 (4%), TP53 (5%) and ERBB2 (6%). We also observe somatic ELF3 (13%) and CBFB (8%) mutations in 24 adenocarcinomas. Squamous cell carcinomas have higher frequencies of somatic nucleotide substitutions occurring at cytosines preceded by thymines (Tp*C sites) than adenocarcinomas. Gene expression levels at HPV integration sites were statistically significantly higher in tumours with HPV integration compared with expression of the same genes in tumours without viral integration at the same site. These data demonstrate several recurrent genomic alterations in cervical carcinomas that suggest new strategies to combat this disease.

Widespread genetic heterogeneity in multiple myeloma: implications for targeted therapy.

We performed massively parallel sequencing of paired tumor/normal samples from 203 multiple myeloma (MM) patients and identified significantly mutated genes and copy number alterations and discovered putative tumor suppressor genes by determining homozygous deletions and loss of heterozygosity. We observed frequent mutations in KRAS (particularly in previously treated patients), NRAS, BRAF, FAM46C, TP53, and DIS3 (particularly in nonhyperdiploid MM). Mutations were often present in subclonal populations, and multiple mutations within the same pathway (e.g., KRAS, NRAS, and BRAF) were observed in the same patient. In vitro modeling predicts only partial treatment efficacy of targeting subclonal mutations, and even growth promotion of nonmutated subclones in some cases. These results emphasize the importance of heterogeneity analysis for treatment decisions.

Integrative and comparative genomic analysis of lung squamous cell carcinomas in East Asian patients.

PURPOSE: Lung squamous cell carcinoma (SCC) is the second most prevalent type of lung cancer. Currently, no targeted therapeutics are approved for treatment of this cancer, largely because of a lack of systematic understanding of the molecular pathogenesis of the disease. To identify therapeutic targets and perform comparative analyses of lung SCC, we probed somatic genome alterations of lung SCC by using samples from Korean patients. PATIENTS AND METHODS: We performed whole-exome sequencing of DNA from 104 lung SCC samples from Korean patients and matched normal DNA. In addition, copy-number analysis and transcriptome analysis were conducted for a subset of these samples. Clinical association with cancer-specific somatic alterations was investigated. RESULTS: This cancer cohort is characterized by a high mutational burden with an average of 261 somatic exonic mutations per tumor and a mutational spectrum showing a signature of exposure to cigarette smoke. Seven genes demonstrated statistical enrichment for mutation: TP53, RB1, PTEN, NFE2L2, KEAP1, MLL2, and PIK3CA). Comparative analysis between Korean and North American lung SCC samples demonstrated a similar spectrum of alterations in these two populations in contrast to the differences seen in lung adenocarcinoma. We also uncovered recurrent occurrence of therapeutically actionable FGFR3-TACC3 fusion in lung SCC. CONCLUSION: These findings provide new steps toward the identification of genomic target candidates for precision medicine in lung SCC, a disease with significant unmet medical needs.

The genetic landscape of clinical resistance to RAF inhibition in metastatic melanoma.

Most patients with BRAF(V600)-mutant metastatic melanoma develop resistance to selective RAF kinase inhibitors. The spectrum of clinical genetic resistance mechanisms to RAF inhibitors and options for salvage therapy are incompletely understood. We performed whole-exome sequencing on formalin-fixed, paraffin-embedded tumors from 45 patients with BRAF(V600)-mutant metastatic melanoma who received vemurafenib or dabrafenib monotherapy. Genetic alterations in known or putative RAF inhibitor resistance genes were observed in 23 of 45 patients (51%). Besides previously characterized alterations, we discovered a "long tail" of new mitogen-activated protein kinase (MAPK) pathway alterations (MAP2K2, MITF) that confer RAF inhibitor resistance. In three cases, multiple resistance gene alterations were observed within the same tumor biopsy. Overall, RAF inhibitor therapy leads to diverse clinical genetic resistance mechanisms, mostly involving MAPK pathway reactivation. Novel therapeutic combinations may be needed to achieve durable clinical control of BRAF(V600)-mutant melanoma. Integrating clinical genomics with preclinical screens may model subsequent resistance studies.