Cancer is driven by somatically acquired point mutations and chromosomal rearrangements, conventionally thought to accumulate gradually over time. Using next-generation sequencing, we characterize a phenomenon, which we term chromothripsis, whereby tens to hundreds of genomic rearrangements occur in a one-off cellular crisis. Rearrangements involving one or a few chromosomes crisscross back and forth across involved regions, generating frequent oscillations between two copy number states. These genomic hallmarks are highly improbable if rearrangements accumulate over time and instead imply that nearly all occur during a single cellular catastrophe. The stamp of chromothripsis can be seen in at least 2%-3% of all cancers, across many subtypes, and is present in ∼25% of bone cancers. We find that one, or indeed more than one, cancer-causing lesion can emerge out of the genomic crisis. This phenomenon has important implications for the origins of genomic remodeling and temporal emergence of cancer.

The genetics of renal cancer is dominated by inactivation of the VHL tumour suppressor gene in clear cell carcinoma (ccRCC), the commonest histological subtype. A recent large-scale screen of ∼3,500 genes by PCR-based exon re-sequencing identified several new cancer genes in ccRCC including UTX (also known as KDM6A), JARID1C (also known as KDM5C) and SETD2 (ref. 2). These genes encode enzymes that demethylate (UTX, JARID1C) or methylate (SETD2) key lysine residues of histone H3. Modification of the methylation state of these lysine residues of histone H3 regulates chromatin structure and is implicated in transcriptional control. However, together these mutations are present in fewer than 15% of ccRCC, suggesting the existence of additional, currently unidentified cancer genes. Here, we have sequenced the protein coding exome in a series of primary ccRCC and report the identification of the SWI/SNF chromatin remodelling complex gene PBRM1 (ref. 4) as a second major ccRCC cancer gene, with truncating mutations in 41% (92/227) of cases. These data further elucidate the somatic genetic architecture of ccRCC and emphasize the marked contribution of aberrant chromatin biology.

Multifocal breast cancer (MFBC), defined as multiple synchronous unilateral lesions of invasive breast cancer, is relatively frequent and has been associated with more aggressive features than unifocal cancer. Here, we aimed to investigate the genomic heterogeneity between MFBC lesions sharing similar histopathological parameters. Characterization of different lesions from 36 patients with ductal MFBC involved the identification of non-silent coding mutations in 360 protein-coding genes (171 tumour and 36 matched normal samples). We selected only patients with lesions presenting the same grade, ER, and HER2 status. Mutations were classified as 'oncogenic' in the case of recurrent substitutions reported in COSMIC or truncating mutations affecting tumour suppressor genes. All mutations identified in a given patient were further interrogated in all samples from that patient through deep resequencing using an orthogonal platform. Whole-genome rearrangement screen was further conducted in 8/36 patients. Twenty-four patients (67%) had substitutions/indels shared by all their lesions, of which 11 carried the same mutations in all lesions, and 13 had lesions with both common and private mutations. Three-quarters of those 24 patients shared oncogenic variants. The remaining 12 patients (33%) did not share any substitution/indels, with inter-lesion heterogeneity observed for oncogenic mutation(s) in genes such as PIK3CA, TP53, GATA3, and PTEN. Genomically heterogeneous lesions tended to be further apart in the mammary gland than homogeneous lesions. Genome-wide analyses of a limited number of patients identified a common somatic background in all studied MFBCs, including those with no mutation in common between the lesions. To conclude, as the number of molecular targeted therapies increases and trials driven by genomic screening are ongoing, our findings highlight the presence of genomic inter-lesion heterogeneity in one-third, despite similar pathological features. This implies that deeper molecular characterization of all MFBC lesions is warranted for the adequate management of those cancers.

BRAF (v-raf murine sarcoma viral oncogene homolog B) inhibitors elicit a transient anti-tumor response in ∼ 80% of BRAF(V600)-mutant melanoma patients that almost uniformly precedes the emergence of resistance. Here we used a mouse model of melanoma in which melanocyte-specific expression of Braf(V618E) (analogous to the human BRAF(V600E) mutation) led to the development of skin hyperpigmentation and nevi, as well as melanoma formation with incomplete penetrance. Sleeping Beauty insertional mutagenesis in this model led to accelerated and fully penetrant melanomagenesis and synchronous tumor formation. Treatment of Braf(V618E) transposon mice with the BRAF inhibitor PLX4720 resulted in tumor regression followed by relapse. Analysis of transposon insertions identified eight genes including Braf, Mitf, and ERas (ES-cell expressed Ras) as candidate resistance genes. Expression of ERAS in human melanoma cell lines conferred resistance to PLX4720 and induced hyperphosphorylation of AKT (v-akt murine thymoma viral oncogene homolog 1), a phenotype reverted by combinatorial treatment with PLX4720 and the AKT inhibitor MK2206. We show that ERAS expression elicits a prosurvival signal associated with phosphorylation/inactivation of BAD, and that the resistance of hepatocyte growth factor-treated human melanoma cells to PLX4720 can be reverted by treatment with the BAD-like BH3 mimetic ABT-737. Thus, we define a role for the AKT/BAD pathway in resistance to BRAF inhibition and illustrate an in vivo approach for finding drug resistance genes.

Mutation is associated with developmental and hereditary disorders, aging, and cancer. While we understand some mutational processes operative in human disease, most remain mysterious. We used Caenorhabditis elegans whole-genome sequencing to model mutational signatures, analyzing 183 worm populations across 17 DNA repair-deficient backgrounds propagated for 20 generations or exposed to carcinogens. The baseline mutation rate in C. elegans was approximately one per genome per generation, not overtly altered across several DNA repair deficiencies over 20 generations. Telomere erosion led to complex chromosomal rearrangements initiated by breakage-fusion-bridge cycles and completed by simultaneously acquired, localized clusters of breakpoints. Aflatoxin B1 induced substitutions of guanines in a GpC context, as observed in aflatoxin-induced liver cancers. Mutational burden increased with impaired nucleotide excision repair. Cisplatin and mechlorethamine, DNA crosslinking agents, caused dose- and genotype-dependent signatures among indels, substitutions, and rearrangements. Strikingly, both agents induced clustered rearrangements resembling "chromoanasynthesis," a replication-based mutational signature seen in constitutional genomic disorders, suggesting that interstrand crosslinks may play a pathogenic role in such events. Cisplatin mutagenicity was most pronounced in xpf-1 mutants, suggesting that this gene critically protects cells against platinum chemotherapy. Thus, experimental model systems combined with genome sequencing can recapture and mechanistically explain mutational signatures associated with human disease.

Evidence for distinct human cancer stem cells (CSCs) remains contentious and the degree to which different cancer cells contribute to propagating malignancies in patients remains unexplored. In low- to intermediate-risk myelodysplastic syndromes (MDS), we establish the existence of rare multipotent MDS stem cells (MDS-SCs), and their hierarchical relationship to lineage-restricted MDS progenitors. All identified somatically acquired genetic lesions were backtracked to distinct MDS-SCs, establishing their distinct MDS-propagating function in vivo. In isolated del(5q)-MDS, acquisition of del(5q) preceded diverse recurrent driver mutations. Sequential analysis in del(5q)-MDS revealed genetic evolution in MDS-SCs and MDS-progenitors prior to leukemic transformation. These findings provide definitive evidence for rare human MDS-SCs in vivo, with extensive implications for the targeting of the cells required and sufficient for MDS-propagation.

Ionizing radiation is a potent carcinogen, inducing cancer through DNA damage. The signatures of mutations arising in human tissues following in vivo exposure to ionizing radiation have not been documented. Here, we searched for signatures of ionizing radiation in 12 radiation-associated second malignancies of different tumour types. Two signatures of somatic mutation characterize ionizing radiation exposure irrespective of tumour type. Compared with 319 radiation-naive tumours, radiation-associated tumours carry a median extra 201 deletions genome-wide, sized 1-100 base pairs often with microhomology at the junction. Unlike deletions of radiation-naive tumours, these show no variation in density across the genome or correlation with sequence context, replication timing or chromatin structure. Furthermore, we observe a significant increase in balanced inversions in radiation-associated tumours. Both small deletions and inversions generate driver mutations. Thus, ionizing radiation generates distinctive mutational signatures that explain its carcinogenic potential.

Liver cirrhosis is associated with decreased hepatic cytochrome P4503A (CYP3A) activity but the pathogenesis of this phenomenon is not well elucidated. In this study, we examined if certain microRNAs (miRNA) are associated with decreased hepatic CYP3A activity in cirrhosis.

The emergence of high-throughput, next-generation sequencing technologies has dramatically altered the way we assess genomes in population genetics and in cancer genomics. Currently, there are four commonly used whole-genome sequencing platforms on the market: Illumina's HiSeq2000, Life Technologies' SOLiD 4 and its completely redesigned 5500xl SOLiD, and Complete Genomics' technology. A number of earlier studies have compared a subset of those sequencing platforms or compared those platforms with Sanger sequencing, which is prohibitively expensive for whole genome studies. Here we present a detailed comparison of the performance of all currently available whole genome sequencing platforms, especially regarding their ability to call SNVs and to evenly cover the genome and specific genomic regions. Unlike earlier studies, we base our comparison on four different samples, allowing us to assess the between-sample variation of the platforms. We find a pronounced GC bias in GC-rich regions for Life Technologies' platforms, with Complete Genomics performing best here, while we see the least bias in GC-poor regions for HiSeq2000 and 5500xl. HiSeq2000 gives the most uniform coverage and displays the least sample-to-sample variation. In contrast, Complete Genomics exhibits by far the smallest fraction of bases not covered, while the SOLiD platforms reveal remarkable shortcomings, especially in covering CpG islands. When comparing the performance of the four platforms for calling SNPs, HiSeq2000 and Complete Genomics achieve the highest sensitivity, while the SOLiD platforms show the lowest false positive rate. Finally, we find that integrating sequencing data from different platforms offers the potential to combine the strengths of different technologies. In summary, our results detail the strengths and weaknesses of all four whole-genome sequencing platforms. It indicates application areas that call for a specific sequencing platform and disallow other platforms. This helps to identify the proper sequencing platform for whole genome studies with different application scopes.

Breast cancer genomes have revealed a novel form of mutation showers (kataegis) in which multiple same-strand substitutions at C:G pairs spaced one to several hundred nucleotides apart are clustered over kilobase-sized regions, often associated with sites of DNA rearrangement. We show kataegis can result from AID/APOBEC-catalysed cytidine deamination in the vicinity of DNA breaks, likely through action on single-stranded DNA exposed during resection. Cancer-like kataegis can be recapitulated by expression of AID/APOBEC family deaminases in yeast where it largely depends on uracil excision, which generates an abasic site for strand breakage. Localized kataegis can also be nucleated by an I-SceI-induced break. Genome-wide patterns of APOBEC3-catalyzed deamination in yeast reveal APOBEC3B and 3A as the deaminases whose mutational signatures are most similar to those of breast cancer kataegic mutations. Together with expression and functional assays, the results implicate APOBEC3B/A in breast cancer hypermutation and give insight into the mechanism of kataegis. DOI:http://dx.doi.org/10.7554/eLife.00534.001.

We analysed whole-genome sequences of 560 breast cancers to advance understanding of the driver mutations conferring clonal advantage and the mutational processes generating somatic mutations. We found that 93 protein-coding cancer genes carried probable driver mutations. Some non-coding regions exhibited high mutation frequencies, but most have distinctive structural features probably causing elevated mutation rates and do not contain driver mutations. Mutational signature analysis was extended to genome rearrangements and revealed twelve base substitution and six rearrangement signatures. Three rearrangement signatures, characterized by tandem duplications or deletions, appear associated with defective homologous-recombination-based DNA repair: one with deficient BRCA1 function, another with deficient BRCA1 or BRCA2 function, the cause of the third is unknown. This analysis of all classes of somatic mutation across exons, introns and intergenic regions highlights the repertoire of cancer genes and mutational processes operating, and progresses towards a comprehensive account of the somatic genetic basis of breast cancer.

Exposure to extracellular 5'-adenosine triphosphate (ATP) is known to induce membrane blebbing. In this study, we investigated the subcellular distribution of the cytoskeletal adaptor protein paxillin in primary bovine osteoblasts upon stimulation with ATP. Cells expressing a fusion protein of green fluorescent protein (GFP) and paxillin were followed by time-lapse video-microscopy after stimulation with 100 μM ATP. Within 100 s, GFP-paxillin became incorporated in numerous de novo formed focal aggregates localized at the cell periphery. The assembly of individual paxillin-containing aggregates occurred with a mean half-life time of <60 s, whereas their disassembly lasted twice as long. Despite the ongoing presence of ATP, the formation of paxillin aggregates was self-limiting within 25 min. Paxillin clustering was preceded by a transient rise in cytoplasmic calcium transients, which peaked already 20 s after adding ATP. The high mobility of paxillin was confirmed by measuring the dissociation rate of GFP-paxillin at mature focal adhesions, demonstrating the presence of a highly mobile fraction with a mean recovery half-life of 8.2 ± 1.2 s, followed by a slower phase (53 ± 20 s). Thus, both the exchange of paxillin at mature focal adhesions and the increase in intracellular calcium concentrations upon ATP stimulation are very rapid processes, which override the time course of ATP-induced paxillin membrane clustering by one to two orders of magnitude. Our data demonstrate that the transient recruitment of paxillin in membrane protuberances is based on the high intracytoplasmic mobility of unbound paxillin molecules and their rapid focal accumulation.

The genome of a cancer cell carries somatic mutations that are the cumulative consequences of the DNA damage and repair processes operative during the cellular lineage between the fertilized egg and the cancer cell. Remarkably, these mutational processes are poorly characterized. Global sequencing initiatives are yielding catalogs of somatic mutations from thousands of cancers, thus providing the unique opportunity to decipher the signatures of mutational processes operative in human cancer. However, until now there have been no theoretical models describing the signatures of mutational processes operative in cancer genomes and no systematic computational approaches are available to decipher these mutational signatures. Here, by modeling mutational processes as a blind source separation problem, we introduce a computational framework that effectively addresses these questions. Our approach provides a basis for characterizing mutational signatures from cancer-derived somatic mutational catalogs, paving the way to insights into the pathogenetic mechanism underlying all cancers.

In a systematic sequencing screen of the coding exons of the X chromosome in 250 families with X-linked mental retardation (XLMR), we identified two nonsense mutations and one consensus splice-site mutation in the AP1S2 gene on Xp22 in three families. Affected individuals in these families showed mild-to-profound mental retardation. Other features included hypotonia early in life and delay in walking. AP1S2 encodes an adaptin protein that constitutes part of the adaptor protein complex found at the cytoplasmic face of coated vesicles located at the Golgi complex. The complex mediates the recruitment of clathrin to the vesicle membrane. Aberrant endocytic processing through disruption of adaptor protein complexes is likely to result from the AP1S2 mutations identified in the three XLMR-affected families, and such defects may plausibly cause abnormal synaptic development and function. AP1S2 is the first reported XLMR gene that encodes a protein directly involved in the assembly of endocytic vesicles.

Cytoreductive therapy is beneficial in patients with essential thrombocythemia (ET) at high risk of thrombosis. However, its value in those lacking high-risk features remains unknown. This open-label, randomized trial compared hydroxycarbamide plus aspirin with aspirin alone in patients with ET age 40 to 59 years and without high-risk factors or extreme thrombocytosis.

Erroneous DNA repair by heterologous recombination (Ht-REC) is a potential threat to genome stability, but evidence supporting its prevalence is lacking. Here we demonstrate that recombination is possible between heterologous sequences and that it is a source of chromosomal alterations in mitotic and meiotic cells. Mechanistically, we find that the RTEL1 and HIM-6/BLM helicases and the BRCA1 homolog BRC-1 counteract Ht-REC in Caenorhabditis elegans, whereas mismatch repair does not. Instead, MSH-2/6 drives Ht-REC events in rtel-1 and brc-1 mutants and excessive crossovers in rtel-1 mutant meioses. Loss of vertebrate Rtel1 also causes a variety of unusually large and complex structural variations, including chromothripsis, breakage-fusion-bridge events, and tandem duplications with distant intra-chromosomal insertions, whose structure are consistent with a role for RTEL1 in preventing Ht-REC during break-induced replication. Our data establish Ht-REC as an unappreciated source of genome instability that underpins a novel class of complex genome rearrangements that likely arise during replication stress.

Osteosarcoma is a primary malignancy of bone that affects children and adults. Here, we present the largest sequencing study of osteosarcoma to date, comprising 112 childhood and adult tumours encompassing all major histological subtypes. A key finding of our study is the identification of mutations in insulin-like growth factor (IGF) signalling genes in 8/112 (7%) of cases. We validate this observation using fluorescence in situ hybridization (FISH) in an additional 87 osteosarcomas, with IGF1 receptor (IGF1R) amplification observed in 14% of tumours. These findings may inform patient selection in future trials of IGF1R inhibitors in osteosarcoma. Analysing patterns of mutation, we identify distinct rearrangement profiles including a process characterized by chromothripsis and amplification. This process operates recurrently at discrete genomic regions and generates driver mutations. It may represent an age-independent mutational mechanism that contributes to the development of osteosarcoma in children and adults alike.

Adamantinomatous craniopharyngiomas (ACPs) are clinically challenging tumours, the majority of which have activating mutations in CTNNB1. They are histologically complex, showing cystic and solid components, the latter comprised of different morphological cell types (e.g. β-catenin-accumulating cluster cells and palisading epithelium), surrounded by a florid glial reaction with immune cells. Here, we have carried out RNA sequencing on 18 ACP samples and integrated these data with an existing ACP transcriptomic dataset. No studies so far have examined the patterns of gene expression within the different cellular compartments of the tumour. To achieve this goal, we have combined laser capture microdissection with computational analyses to reveal groups of genes that are associated with either epithelial tumour cells (clusters and palisading epithelium), glial tissue or immune infiltrate. We use these human ACP molecular signatures and RNA-Seq data from two ACP mouse models to reveal that cell clusters are molecularly analogous to the enamel knot, a critical signalling centre controlling normal tooth morphogenesis. Supporting this finding, we show that human cluster cells express high levels of several members of the FGF, TGFB and BMP families of secreted factors, which signal to neighbouring cells as evidenced by immunostaining against the phosphorylated proteins pERK1/2, pSMAD3 and pSMAD1/5/9 in both human and mouse ACP. We reveal that inhibiting the MAPK/ERK pathway with trametinib, a clinically approved MEK inhibitor, results in reduced proliferation and increased apoptosis in explant cultures of human and mouse ACP. Finally, we analyse a prominent molecular signature in the glial reactive tissue to characterise the inflammatory microenvironment and uncover the activation of inflammasomes in human ACP. We validate these results by immunostaining against immune cell markers, cytokine ELISA and proteome analysis in both solid tumour and cystic fluid from ACP patients. Our data support a new molecular paradigm for understanding ACP tumorigenesis as an aberrant mimic of natural tooth development and opens new therapeutic opportunities by revealing the activation of the MAPK/ERK and inflammasome pathways in human ACP.

Transmissible cancers are clonal lineages that spread through populations via contagious cancer cells. Although rare in nature, two facial tumor clones affect Tasmanian devils. Here we perform comparative genetic and functional characterization of these lineages. The two cancers have similar patterns of mutation and show no evidence of exposure to exogenous mutagens or viruses. Genes encoding PDGF receptors have copy number gains and are present on extrachromosomal double minutes. Drug screening indicates causative roles for receptor tyrosine kinases and sensitivity to inhibitors of DNA repair. Y chromosome loss from a male clone infecting a female host suggests immunoediting. These results imply that Tasmanian devils may have inherent susceptibility to transmissible cancers and present a suite of therapeutic compounds for use in conservation.

Undifferentiated sarcomas (USARCs) of adults are diverse, rare, and aggressive soft tissue cancers. Recent sequencing efforts have confirmed that USARCs exhibit one of the highest burdens of structural aberrations across human cancer. Here, we sought to unravel the molecular basis of the structural complexity in USARCs by integrating DNA sequencing, ploidy analysis, gene expression, and methylation profiling. We identified whole genome duplication as a prevalent and pernicious force in USARC tumorigenesis. Using mathematical deconvolution strategies to unravel the complex copy-number profiles and mutational timing models we infer distinct evolutionary pathways of these rare cancers. In addition, 15% of tumors exhibited raised mutational burdens that correlated with gene expression signatures of immune infiltration, and good prognosis.