Characterizing genetic intra-tumor heterogeneity across 2,658 human cancer genomes.

Intra-tumor heterogeneity (ITH) is a mechanism of therapeutic resistance and therefore an important clinical challenge. However, the extent, origin, and drivers of ITH across cancer types are poorly understood. To address this, we extensively characterize ITH across whole-genome sequences of 2,658 cancer samples spanning 38 cancer types. Nearly all informative samples (95.1%) contain evidence of distinct subclonal expansions with frequent branching relationships between subclones. We observe positive selection of subclonal driver mutations across most cancer types and identify cancer type-specific subclonal patterns of driver gene mutations, fusions, structural variants, and copy number alterations as well as dynamic changes in mutational processes between subclonal expansions. Our results underline the importance of ITH and its drivers in tumor evolution and provide a pan-cancer resource of comprehensively annotated subclonal events from whole-genome sequencing data.

The evolutionary history of 2,658 cancers.

Cancer develops through a process of somatic evolution1,2. Sequencing data from a single biopsy represent a snapshot of this process that can reveal the timing of specific genomic aberrations and the changing influence of mutational processes3. Here, by whole-genome sequencing analysis of 2,658 cancers as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)4, we reconstruct the life history and evolution of mutational processes and driver mutation sequences of 38 types of cancer. Early oncogenesis is characterized by mutations in a constrained set of driver genes, and specific copy number gains, such as trisomy 7 in glioblastoma and isochromosome 17q in medulloblastoma. The mutational spectrum changes significantly throughout tumour evolution in 40% of samples. A nearly fourfold diversification of driver genes and increased genomic instability are features of later stages. Copy number alterations often occur in mitotic crises, and lead to simultaneous gains of chromosomal segments. Timing analyses suggest that driver mutations often precede diagnosis by many years, if not decades. Together, these results determine the evolutionary trajectories of cancer, and highlight opportunities for early cancer detection.

Genetic events associated with venetoclax resistance in CLL identified by whole-exome sequencing of patient samples.

Although BCL2 mutations are reported as later occurring events leading to venetoclax resistance, many other mechanisms of progression have been reported though remain poorly understood. Here, we analyze longitudinal tumor samples from 11 patients with disease progression while receiving venetoclax to characterize the clonal evolution of resistance. All patients tested showed increased in vitro resistance to venetoclax at the posttreatment time point. We found the previously described acquired BCL2-G101V mutation in only 4 of 11 patients, with 2 patients showing a very low variant allele fraction (0.03%-4.68%). Whole-exome sequencing revealed acquired loss(8p) in 4 of 11 patients, of which 2 patients also had gain (1q21.2-21.3) in the same cells affecting the MCL1 gene. In vitro experiments showed that CLL cells from the 4 patients with loss(8p) were more resistant to venetoclax than cells from those without it, with the cells from 2 patients also carrying gain (1q21.2-21.3) showing increased sensitivity to MCL1 inhibition. Progression samples with gain (1q21.2-21.3) were more susceptible to the combination of MCL1 inhibitor and venetoclax. Differential gene expression analysis comparing bulk RNA sequencing data from pretreatment and progression time points of all patients showed upregulation of proliferation, B-cell receptor (BCR), and NF-κB gene sets including MAPK genes. Cells from progression time points demonstrated upregulation of surface immunoglobulin M and higher pERK levels compared with those from the preprogression time point, suggesting an upregulation of BCR signaling that activates the MAPK pathway. Overall, our data suggest several mechanisms of acquired resistance to venetoclax in CLL that could pave the way for rationally designed combination treatments for patients with venetoclax-resistant CLL.

Growth dynamics in naturally progressing chronic lymphocytic leukaemia.

How the genomic features of a patient's cancer relate to individual disease kinetics remains poorly understood. Here we used the indolent growth dynamics of chronic lymphocytic leukaemia (CLL) to analyse the growth rates and corresponding genomic patterns of leukaemia cells from 107 patients with CLL, spanning decades-long disease courses. We found that CLL commonly demonstrates not only exponential expansion but also logistic growth, which is sigmoidal and reaches a certain steady-state level. Each growth pattern was associated with marked differences in genetic composition, the pace of disease progression and the extent of clonal evolution. In a subset of patients, whose serial samples underwent next-generation sequencing, we found that dynamic changes in the disease course of CLL were shaped by the genetic events that were already present in the early slow-growing stages. Finally, by analysing the growth rates of subclones compared with their parental clones, we quantified the growth advantage conferred by putative CLL drivers in vivo.

Activation of the MAPK pathway mediates resistance to PI3K inhibitors in chronic lymphocytic leukemia.

Inhibitors of Bruton tyrosine kinase (BTK) and phosphatidylinositol 3-kinase δ (PI3Kδ) that target the B-cell receptor (BCR) signaling pathway have revolutionized the treatment of chronic lymphocytic leukemia (CLL). Mutations associated with resistance to BTK inhibitors have been identified, but limited data are available on mechanisms of resistance to PI3Kδ inhibitors. Here we present findings from longitudinal whole-exome sequencing of cells from patients with multiply relapsed CLL (N = 28) enrolled in trials of PI3K inhibitors. The nonresponder subgroup was characterized by baseline activating mutations in MAP2K1, BRAF, and KRAS genes in 60% of patients. PI3Kδ inhibition failed to inhibit ERK phosphorylation (pERK) in nonresponder CLL cells with and without mutations, whereas treatment with a MEK inhibitor rescued ERK inhibition. Overexpression of MAP2K1 mutants in vitro led to increased basal and inducible pERK and resistance to idelalisib. These data demonstrate that MAPK/ERK activation plays a key role in resistance to PI3Kδ inhibitors in CLL and provide a rationale for therapy with a combination of PI3Kδ and ERK inhibitors.

The cationic amino acid exporter Slc7a7 is induced and vital in zebrafish tissue macrophages with sustained efferocytic activity.

Most tissues harbor a substantial population of resident macrophages. Here, we elucidate a functional link between the Slc7a7 cationic amino acid transporter and tissue macrophages. We identified a mutant zebrafish devoid of microglia due to a mutation in the slc7a7 gene. We found that in Slc7a7-deficient larvae, macrophages do enter the retina and brain to become microglia, but then die during the developmental wave of neuronal apoptosis, which triggers intense efferocytic work from them. A similar macrophage demise occurs in other tissues, at stages where macrophages have to engulf many cell corpses, whether due to developmental or experimentally triggered cell death. We found that Slc7a7 is the main cationic amino acid transporter expressed in macrophages of zebrafish larvae, and that its expression is induced in tissue macrophages within 1-2 h upon efferocytosis. Our data indicate that Slc7a7 is vital not only for microglia but also for any steadily efferocytic tissue macrophages, and that slc7a7 gene induction is one of the adaptive responses that allow them to cope with the catabolism of numerous dead cells without compromising their own viability.

Tfap2a is a novel gatekeeper of nephron differentiation during kidney development.

Renal functional units known as nephrons undergo patterning events during development that create a segmental array of cellular compartments with discrete physiological identities. Here, from a forward genetic screen using zebrafish, we report the discovery that transcription factor AP-2 alpha (tfap2a) coordinates a gene regulatory network that activates the terminal differentiation program of distal segments in the pronephros. We found that tfap2a acts downstream of Iroquois homeobox 3b (irx3b), a distal lineage transcription factor, to operate a circuit consisting of tfap2b, irx1a and genes encoding solute transporters that dictate the specialized metabolic functions of distal nephron segments. Interestingly, this regulatory node is distinct from other checkpoints of differentiation, such as polarity establishment and ciliogenesis. Thus, our studies reveal insights into the genetic control of differentiation, where tfap2a is essential for regulating a suite of segment transporter traits at the final tier of zebrafish pronephros ontogeny. These findings have relevance for understanding renal birth defects, as well as efforts to recapitulate nephrogenesis in vivo to facilitate drug discovery and regenerative therapies.

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.).

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

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.

The zebrafish kidney mutant zeppelin reveals that brca2/fancd1 is essential for pronephros development.

The zebrafish kidney is conserved with other vertebrates, making it an excellent genetic model to study renal development. The kidney collects metabolic waste using a blood filter with specialized epithelial cells known as podocytes. Podocyte formation is poorly understood but relevant to many kidney diseases, as podocyte injury leads to progressive scarring and organ failure. zeppelin (zep) was isolated in a forward screen for kidney mutants and identified as a homozygous recessive lethal allele that causes reduced podocyte numbers, deficient filtration, and fluid imbalance. Interestingly, zep mutants had a larger interrenal gland, the teleostean counterpart of the mammalian adrenal gland, which suggested a fate switch with the related podocyte lineage since cell proliferation and cell death were unchanged within the shared progenitor field from which these two identities arise. Cloning of zep by whole genome sequencing (WGS) identified a splicing mutation in breast cancer 2, early onset (brca2)/fancd1, which was confirmed by sequencing of individual fish. Several independent brca2 morpholinos (MOs) phenocopied zep, causing edema, reduced podocyte number, and increased interrenal cell number. Complementation analysis between zep and brca2ZM_00057434 -/- zebrafish, which have an insertional mutation, revealed that the interrenal lineage was expanded. Importantly, overexpression of brca2 rescued podocyte formation in zep mutants, providing critical evidence that the brca2 lesion encoded by zep specifically disrupts the balance of nephrogenesis. Taken together, these data suggest for the first time that brca2/fancd1 is essential for vertebrate kidney ontogeny. Thus, our findings impart novel insights into the genetic components that impact renal development, and because BRCA2/FANCD1 mutations in humans cause Fanconi anemia and several common cancers, this work has identified a new zebrafish model to further study brca2/fancd1 in disease.

Comprehensive assessment of cancer missense mutation clustering in protein structures.

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.

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.

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.

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.

Crowd-sourced benchmarking of single-sample tumor subclonal reconstruction.

Subclonal reconstruction algorithms use bulk DNA sequencing data to quantify parameters of tumor evolution, allowing an assessment of how cancers initiate, progress and respond to selective pressures. We launched the ICGC-TCGA (International Cancer Genome Consortium-The Cancer Genome Atlas) DREAM Somatic Mutation Calling Tumor Heterogeneity and Evolution Challenge to benchmark existing subclonal reconstruction algorithms. This 7-year community effort used cloud computing to benchmark 31 subclonal reconstruction algorithms on 51 simulated tumors. Algorithms were scored on seven independent tasks, leading to 12,061 total runs. Algorithm choice influenced performance substantially more than tumor features but purity-adjusted read depth, copy-number state and read mappability were associated with the performance of most algorithms on most tasks. No single algorithm was a top performer for all seven tasks and existing ensemble strategies were unable to outperform the best individual methods, highlighting a key research need. All containerized methods, evaluation code and datasets are available to support further assessment of the determinants of subclonal reconstruction accuracy and development of improved methods to understand tumor evolution.

Allosteric PI3Kα Inhibition Overcomes On-target Resistance to Orthosteric Inhibitors Mediated by Secondary PIK3CA Mutations.

PIK3CA mutations occur in ∼8% of cancers, including ∼40% of HR-positive breast cancers, where the PI3K-alpha (PI3Kα)-selective inhibitor alpelisib is FDA approved in combination with fulvestrant. Although prior studies have identified resistance mechanisms, such as PTEN loss, clinically acquired resistance to PI3Kα inhibitors remains poorly understood. Through serial liquid biopsies and rapid autopsies in 39 patients with advanced breast cancer developing acquired resistance to PI3Kα inhibitors, we observe that 50% of patients acquire genomic alterations within the PI3K pathway, including PTEN loss and activating AKT1 mutations. Notably, although secondary PIK3CA mutations were previously reported to increase sensitivity to PI3Kα inhibitors, we identified emergent secondary resistance mutations in PIK3CA that alter the inhibitor binding pocket. Some mutations had differential effects on PI3Kα-selective versus pan-PI3K inhibitors, but resistance induced by all mutations could be overcome by the novel allosteric pan-mutant-selective PI3Kα-inhibitor RLY-2608. Together, these findings provide insights to guide strategies to overcome resistance in PIK3CA-mutated cancers. SIGNIFICANCE: In one of the largest patient cohorts analyzed to date, this study defines the clinical landscape of acquired resistance to PI3Kα inhibitors. Genomic alterations within the PI3K pathway represent a major mode of resistance and identify a novel class of secondary PIK3CA resistance mutations that can be overcome by an allosteric PI3Kα inhibitor. See related commentary by Gong and Vanhaesebroeck, p. 204 . See related article by Varkaris et al., p. 240 . This article is featured in Selected Articles from This Issue, p. 201.

Inferring early genetic progression in cancers with unobtainable premalignant disease.

Analysis of premalignant tissue has identified the typical order of somatic events leading to invasive tumors in several cancer types. For other cancers, premalignant tissue is unobtainable, leaving genetic progression unknown. Here, we demonstrate how to infer progression from exome sequencing of primary tumors. Our computational method, PhylogicNDT, recapitulated the previous experimentally determined genetic progression of human papillomavirus-negative (HPV-) head and neck squamous cell carcinoma (HNSCC). We then evaluated HPV+ HNSCC, which lacks premalignant tissue, and uncovered its previously unknown progression, identifying early drivers. We converted relative timing estimates of driver mutations and HPV integration to years before diagnosis based on a clock-like mutational signature. We associated the timing of transitions to aneuploidy with increased intratumor genetic heterogeneity and shorter overall survival. Our approach can establish previously unknown early genetic progression of cancers with unobtainable premalignant tissue, supporting development of experimental models and methods for early detection, interception and prognostication.

Pirtobrutinib targets BTK C481S in ibrutinib-resistant CLL but second-site BTK mutations lead to resistance.

Covalent inhibitors of Bruton tyrosine kinase (BTK) have transformed the therapy of chronic lymphocytic leukemia (CLL), but continuous therapy has been complicated by the development of resistance. The most common resistance mechanism in patients whose disease progresses on covalent BTK inhibitors (BTKis) is a mutation in the BTK 481 cysteine residue to which the inhibitors bind covalently. Pirtobrutinib is a highly selective, noncovalent BTKi with substantial clinical activity in patients whose disease has progressed on covalent BTKi, regardless of BTK mutation status. Using in vitro ibrutinib-resistant models and cells from patients with CLL, we show that pirtobrutinib potently inhibits BTK-mediated functions including B-cell receptor (BCR) signaling, cell viability, and CCL3/CCL4 chemokine production in both BTK wild-type and C481S mutant CLL cells. We demonstrate that primary CLL cells from responding patients on the pirtobrutinib trial show reduced BCR signaling, cell survival, and CCL3/CCL4 chemokine secretion. At time of progression, these primary CLL cells show increasing resistance to pirtobrutinib in signaling inhibition, cell viability, and cytokine production. We employed longitudinal whole-exome sequencing on 2 patients whose disease progressed on pirtobrutinib and identified selection of alternative-site BTK mutations, providing clinical evidence that secondary BTK mutations lead to resistance to noncovalent BTKis.

Germline mutations and developmental mosaicism underlying EGFR-mutant lung cancer.

While the development of multiple primary tumors in smokers with lung cancer can be attributed to carcinogen-induced field cancerization, the occurrence of multiple primary tumors in individuals with EGFR-mutant lung cancer who lack known environmental exposures remains unexplained. We identified ten patients with early-stage, resectable non-small cell lung cancer who presented with multiple anatomically distinct EGFR-mutant tumors. We analyzed the phylogenetic relationships among multiple tumors from each patient using whole exome sequencing (WES) and hypermutable poly-guanine (poly-G) repeat genotyping, as orthogonal methods for lineage tracing. In two patients, we identified germline EGFR variants, which confer moderately enhanced signaling when modeled in vitro. In four other patients, developmental mosaicism is supported by the poly-G lineage tracing and WES, indicating a common non-germline cell-of-origin. Thus, developmental mosaicism and germline variants define two distinct mechanisms of genetic predisposition to multiple EGFR-mutant primary tumors, with implications for understanding their etiology and clinical management.