Ignaty Leshchiner

Resensitization to Crizotinib by the Lorlatinib Alk Resistance Mutation L1198F.

In a patient who had metastatic anaplastic lymphoma kinase (ALK)-rearranged lung cancer, resistance to crizotinib developed because of a mutation in the ALK kinase domain. This mutation is predicted to result in a substitution of cysteine by tyrosine at amino acid residue 1156 (C1156Y). Her tumor did not respond to a second-generation ALK inhibitor, but it did respond to lorlatinib (PF-06463922), a third-generation inhibitor. When her tumor relapsed, sequencing of the resistant tumor revealed an ALK L1198F mutation in addition to the C1156Y mutation. The L1198F substitution confers resistance to lorlatinib through steric interference with drug binding. However, L1198F paradoxically enhances binding to crizotinib, negating the effect of C1156Y and resensitizing resistant cancers to crizotinib. The patient received crizotinib again, and her cancer-related symptoms and liver failure resolved. (Funded by Pfizer and others; ClinicalTrials.gov number, NCT01970865.).

Molecular Mechanisms of Resistance to First- and Second-Generation ALK Inhibitors in ALK-Rearranged Lung Cancer

Advanced, anaplastic lymphoma kinase (ALK)-positive lung cancer is currently treated with the first-generation ALK inhibitor crizotinib followed by more potent, second-generation ALK inhibitors (e.g., ceritinib, alectinib) upon progression. Second-generation inhibitors are generally effective even in the absence of crizotinib-resistant ALK mutations, likely reflecting incomplete inhibition of ALK by crizotinib in many cases. Herein, we analyzed 103 repeat biopsies from ALK-positive patients progressing on various ALK inhibitors. We find that each ALK inhibitor is associated with a distinct spectrum of ALK resistance mutations and that the frequency of one mutation - ALK G1202R - increases significantly after treatment with second-generation agents. To investigate strategies to overcome resistance to second-generation ALK inhibitors, we examine the activity of the third-generation ALK inhibitor lorlatinib in a series of ceritinib-resistant, patient-derived cell lines, and observe that the presence of ALK resistance mutations is highly predictive for sensitivity to lorlatinib, whereas those cell lines without ALK mutations are resistant.

Mutation mapping and identification by whole genome sequencing

Genetic mapping of mutations in model systems has facilitated the identification of genes contributing to fundamental biological processes including human diseases. However, this approach has historically required the prior characterization of informative markers. Here we report a fast and cost-effective method for genetic mapping using next-generation sequencing that combines single nucleotide polymorphism discovery, mutation localization, and potential identification of causal sequence variants. In contrast to prior approaches, we have developed a hidden Markov model to narrowly define the mutation area by inferring recombination breakpoints of chromosomes in the mutant pool. In addition, we created an interactive online software resource to facilitate automated analysis of sequencing data and demonstrate its utility in the zebrafish and mouse models. Our novel methodology and online tools will make next-generation sequencing an easily applicable resource for mutation mapping in all model systems.

Comprehensive assessment of cancer missense mutation clustering in protein structures

Significance Tumor sequencing efforts have enabled the identification of cancer genes based on an excess of mutations in the gene or clustering of mutations along the (one-dimensional) DNA sequence of the gene. Here, we show that this approach can be extended to identify cancer genes based on clustering of mutations relative to the 3D structure of the protein product. By analyzing the PanCancer compendium of somatic mutations in nearly 5,000 tumors, we identified known cancer genes and previously unidentified candidates based on clustering of missense mutations in protein structures or at interfaces with binding partners. In addition, we found that 3D clustering is present in both oncoproteins and tumor suppressors—contrary to the view that such clustering is a hallmark of oncoproteins. Large-scale tumor sequencing projects enabled the identification of many new cancer gene candidates through computational approaches. Here, we describe a general method to detect cancer genes based on significant 3D clustering of mutations relative to the structure of the encoded protein products. The approach can also be used to search for proteins with an enrichment of mutations at binding interfaces with a protein, nucleic acid, or small molecule partner. We applied this approach to systematically analyze the PanCancer compendium of somatic mutations from 4,742 tumors relative to all known 3D structures of human proteins in the Protein Data Bank. We detected significant 3D clustering of missense mutations in several previously known oncoproteins including HRAS, EGFR, and PIK3CA. Although clustering of missense mutations is often regarded as a hallmark of oncoproteins, we observed that a number of tumor suppressors, including FBXW7, VHL, and STK11, also showed such clustering. Beside these known cases, we also identified significant 3D clustering of missense mutations in NUF2, which encodes a component of the kinetochore, that could affect chromosome segregation and lead to aneuploidy. Analysis of interaction interfaces revealed enrichment of mutations in the interfaces between FBXW7-CCNE1, HRAS-RASA1, CUL4B-CAND1, OGT-HCFC1, PPP2R1A-PPP2R5C/PPP2R2A, DICER1-Mg2+, MAX-DNA, SRSF2-RNA, and others. Together, our results indicate that systematic consideration of 3D structure can assist in the identification of cancer genes and in the understanding of the functional role of their mutations.

Polyclonal Secondary FGFR2 Mutations Drive Acquired Resistance to FGFR Inhibition in Patients with FGFR2 Fusion-Positive Cholangiocarcinoma

Genetic alterations in the fibroblast growth factor receptor (FGFR) pathway are promising therapeutic targets in many cancers, including intrahepatic cholangiocarcinoma (ICC). The FGFR inhibitor BGJ398 displayed encouraging efficacy in patients with FGFR2 fusion-positive ICC in a phase II trial, but the durability of response was limited in some patients. Here, we report the molecular basis for acquired resistance to BGJ398 in three patients via integrative genomic characterization of cell-free circulating tumor DNA (cfDNA), primary tumors, and metastases. Serial analysis of cfDNA demonstrated multiple recurrent point mutations in the FGFR2 kinase domain at progression. Accordingly, biopsy of post-progression lesions and rapid autopsy revealed marked inter- and intralesional heterogeneity, with different FGFR2 mutations in individual resistant clones. Molecular modeling and in vitro studies indicated that each mutation led to BGJ398 resistance and was surmountable by structurally distinct FGFR inhibitors. Thus, polyclonal secondary FGFR2 mutations represent an important clinical resistance mechanism that may guide the development of future therapeutic strategies.Significance: We report the first genetic mechanisms of clinical acquired resistance to FGFR inhibition in patients with FGFR2 fusion-positive ICC. Our findings can inform future strategies for detecting resistance mechanisms and inducing more durable remissions in ICC and in the wide variety of cancers where the FGFR pathway is being explored as a therapeutic target. Cancer Discov; 7(3); 252-63. ©2016 AACR.See related commentary by Smyth et al., p. 248This article is highlighted in the In This Issue feature, p. 235.

A mutational signature reveals alterations underlying deficient homologous recombination repair in breast cancer

Biallelic inactivation of BRCA1 or BRCA2 is associated with a pattern of genome-wide mutations known as signature 3. By analyzing ∼1,000 breast cancer samples, we confirmed this association and established that germline nonsense and frameshift variants in PALB2, but not in ATM or CHEK2, can also give rise to the same signature. We were able to accurately classify missense BRCA1 or BRCA2 variants known to impair homologous recombination (HR) on the basis of this signature. Finally, we show that epigenetic silencing of RAD51C and BRCA1 by promoter methylation is strongly associated with signature 3 and, in our data set, was highly enriched in basal-like breast cancers in young individuals of African descent.

Loss of function mutation in LOX causes thoracic aortic aneurysm and dissection in humans

Significance The mechanical integrity of the arterial wall is dependent on a properly structured ECM. Elastin and collagen are key structural components of the ECM, contributing to the stability and elasticity of normal arteries. Lysyl oxidase (LOX) normally cross-links collagen and elastin molecules in the process of forming proper collagen fibers and elastic lamellae. Here, using whole-genome sequencing in humans and genome engineering in mice, we show that a missense mutation in LOX causes aortic aneurysm and dissection because of insufficient elastin and collagen cross-linking in the aortic wall. These findings confirm mutations in LOX as a cause of aortic disease in humans and identify LOX as a diagnostic and potentially therapeutic target. Thoracic aortic aneurysms and dissections (TAAD) represent a substantial cause of morbidity and mortality worldwide. Many individuals presenting with an inherited form of TAAD do not have causal mutations in the set of genes known to underlie disease. Using whole-genome sequencing in two first cousins with TAAD, we identified a missense mutation in the lysyl oxidase (LOX) gene (c.893T > G encoding p.Met298Arg) that cosegregated with disease in the family. Using clustered regularly interspaced short palindromic repeats (CRISPR)/clustered regularly interspaced short palindromic repeats-associated protein-9 nuclease (Cas9) genome engineering tools, we introduced the human mutation into the homologous position in the mouse genome, creating mice that were heterozygous and homozygous for the human allele. Mutant mice that were heterozygous for the human allele displayed disorganized ultrastructural properties of the aortic wall characterized by fragmented elastic lamellae, whereas mice homozygous for the human allele died shortly after parturition from ascending aortic aneurysm and spontaneous hemorrhage. These data suggest that a missense mutation in LOX is associated with aortic disease in humans, likely through insufficient cross-linking of elastin and collagen in the aortic wall. Mutation carriers may be predisposed to vascular diseases because of weakened vessel walls under stress conditions. LOX sequencing for clinical TAAD may identify additional mutation carriers in the future. Additional studies using our mouse model of LOX-associated TAAD have the potential to clarify the mechanism of disease and identify novel therapeutics specific to this genetic cause.

Rapid identification of kidney cyst mutations by whole exome sequencing in zebrafish

Forward genetic approaches in zebrafish have provided invaluable information about developmental processes. However, the relative difficulty of mapping and isolating mutations has limited the number of new genetic screens. Recent improvements in the annotation of the zebrafish genome coupled to a reduction in sequencing costs prompted the development of whole genome and RNA sequencing approaches for gene discovery. Here we describe a whole exome sequencing (WES) approach that allows rapid and cost-effective identification of mutations. We used our WES methodology to isolate four mutations that cause kidney cysts; we identified novel alleles in two ciliary genes as well as two novel mutants. The WES approach described here does not require specialized infrastructure or training and is therefore widely accessible. This methodology should thus help facilitate genetic screens and expedite the identification of mutants that can inform basic biological processes and the causality of genetic disorders in humans.

Scalable whole-exome sequencing of cell-free DNA reveals high concordance with metastatic tumors

Whole-exome sequencing of cell-free DNA (cfDNA) could enable comprehensive profiling of tumors from blood but the genome-wide concordance between cfDNA and tumor biopsies is uncertain. Here we report ichorCNA, software that quantifies tumor content in cfDNA from 0.1× coverage whole-genome sequencing data without prior knowledge of tumor mutations. We apply ichorCNA to 1439 blood samples from 520 patients with metastatic prostate or breast cancers. In the earliest tested sample for each patient, 34% of patients have ≥10% tumor-derived cfDNA, sufficient for standard coverage whole-exome sequencing. Using whole-exome sequencing, we validate the concordance of clonal somatic mutations (88%), copy number alterations (80%), mutational signatures, and neoantigens between cfDNA and matched tumor biopsies from 41 patients with ≥10% cfDNA tumor content. In summary, we provide methods to identify patients eligible for comprehensive cfDNA profiling, revealing its applicability to many patients, and demonstrate high concordance of cfDNA and metastatic tumor whole-exome sequencing.Identifying the mutational landscape of tumours from cell-free DNA in the blood could help diagnostics in cancer. Here, the authors present ichorCNA, software that quantifies tumour content in cell free DNA, and they demonstrate that cell-free DNA whole-exome sequencing is concordant with metastatic tumour whole-exome sequencing.

Summarizing polygenic risks for complex diseases in a clinical whole-genome report

Purpose:Disease-causing mutations and pharmacogenomic variants are of primary interest for clinical whole-genome sequencing. However, estimating genetic liability for common complex diseases using established risk alleles might one day prove clinically useful.Methods:We compared polygenic scoring methods using a case–control data set with independently discovered risk alleles in the MedSeq Project. For eight traits of clinical relevance in both the primary-care and cardiomyopathy study cohorts, we estimated multiplicative polygenic risk scores using 161 published risk alleles and then normalized them using the population median estimated from the 1000 Genomes Project.Results:Our polygenic score approach identified the overrepresentation of independently discovered risk alleles in cases as compared with controls using a large-scale genome-wide association study data set. In addition to normalized multiplicative polygenic risk scores and rank in a population, the disease prevalence and proportion of heritability explained by known common risk variants provide important context in the interpretation of modern multilocus disease risk models.Conclusion:Our approach in the MedSeq Project demonstrates how complex trait risk variants from an individual genome can be summarized and reported for the general clinician and also highlights the need for definitive clinical studies to obtain reference data for such estimates and to establish clinical utility.Genet Med 17 7, 536–544.