Sam Behjati

Tumor Exome Analysis Reveals Neoantigen-Specific T-Cell Reactivity in an Ipilimumab-Responsive Melanoma

The evidence for T-cell–mediated regression of human cancers such as non–small-cell lung carcinoma, renal cell carcinoma, and—in particular—melanoma after immunotherapy is strong. Anti-CTLA4 (ipilimumab) treatment has been approved for treatment of meta-static melanoma,1 and antibody-mediated blockade of PD-1, a second inhibitory receptor on T cells, has shown highly encouraging results in early clinical trials.2,3 Although the clinical activity of these treatments is apparent, it is still unknown which T-cell reactivities are involved in immunotherapy-induced cancer regression.4 T-cell reactivity against nonmutated tumor-associated self-antigens has been analyzed in patients treated with ipilimumab or with autologous tumor-infiltrating T cells, but the magnitude of the T-cell responses observed has been relatively modest.5,6 In part on the basis of such data, recognition of patient-specific mutant epitopes (hereafter referred to as neoantigens) has been suggested to be a potentially important component.7 A potential involvement of mutated epitopes in T-cell control would also fit well with the observation that the mutation load in sun-exposed melanomas is particularly high.8-10 n nIntriguingly, on the basis of animal model data, it has recently been suggested that (therapy-induced) analysis of T-cell reactivity against patient-specific neoantigens may be feasible through exploitation of cancer genome data.11,12 However, human data have thus far been lacking. Here we report a case of a patient with stage IV melanoma who exhibited a clinical response to ipilimumab treatment. Cancer exome–guided analysis of T-cell reactivity in this patient revealed reactivity against two neoantigens, including a dominant T-cell response against a mutant epitope of the ATR (ataxia telangiectasia and Rad3 related) gene product that increased strongly after ipilimumab treatment. These data provide the first demonstration (to our knowledge) of cancer exome–guided analysis to dissect the effects of melanoma immunotherapy.

Extensive transduction of nonrepetitive DNA mediated by L1 retrotransposition in cancer genomes

Introduction The human genome is peppered with mobile repetitive elements called long interspersed nuclear element–1 (L1) retrotransposons. Propagating through RNA and cDNA intermediates, these molecular parasites copy and insert themselves throughout the genome, with potentially disruptive effects on neighboring genes or regulatory sequences. In the germ line, unique sequence downstream of L1 elements can also be retrotransposed if transcription continues beyond the repeat, a process known as 3′ transduction. There has been growing interest in retrotransposition and 3′ transduction as a possible source of somatic mutations during tumorigenesis. The activity of individual L1 elements fluctuates during tumor evolution. In a lung tumor, hundreds of 3′ transductions arose from a small number of active L1 source elements (colored circles on outer rim of circle). As the tumor evolved from the preinvasive common ancestor to invasive cancer, individual elements exhibited variable activity over time. Rationale To explore whether 3′ transductions are frequent in cancer, we developed a bioinformatic algorithm for identifying somatically acquired retrotranspositions in cancer genomes. We applied our algorithm to 290 cancer samples from 244 patients across 12 tumor types. The unique downstream sequence mobilized with 3′ transductions effectively fingerprints the L1 source element, providing insights into the activity of individual L1 loci across the genome. Results Across the 290 samples, we identified 2756 somatic L1 retrotranspositions. Tumors from 53% of patients had at least one such event, with colorectal and lung cancers being most frequently affected (93% and 75% of patients, respectively). Somatic 3′ transductions comprised 24% of events, half of which represented mobilizations of unique sequence alone, without any accompanying L1 sequence. Overall, 95% of 3′ transductions identified derived from only 72 germline L1 source elements, with as few as four loci accounting for 50% of events. In a given sample, the same source element could generate 50 or more somatic transductions, scattered extensively across the genome. About 5% of somatic transductions arose from L1 source elements that were themselves somatic retrotranspositions. In three of the cases in which we sequenced more than one sample from a patient’s tumor, we were able to place 3′ transductions on the phylogenetic tree. We found that the activity of individual source elements fluctuated during tumor evolution, with different subclones exhibiting much variability in which elements were “on” and which were “off.” The ability to identify the individual L1 source elements active in a given tumor enabled us to study the promoter methylation of those elements specifically. We found that 3′ transduction activity in a patient’s tumor was always associated with hypomethylation of that element. Overall, 2.3% of transductions distributed exons or entire genes to other sites in the genome, and many more mobilized deoxyribonuclease I (DNAse-I) hypersensitive sites or transcription factor binding sites identified by the ENCODE project. Occasionally, somatic L1 insertions inserted near coding sequence and redistributed these exons elsewhere in the genome. However, we found no general effects of retrotranspositions on transcription levels of genes at the insertion points and no evidence for aberrant RNA species resulting from somatically acquired transposable elements. Indeed, as with germline retrotranspositions, somatic insertions exhibited a strong enrichment in heterochromatic, gene-poor regions of the genome. Conclusion Somatic 3′ transduction occurs frequently in human tumors, and in some cases transduction events can scatter exons, genes, and regulatory elements widely across the genome. Dissemination of these sequences appears to be due to a small number of highly active L1 elements, whose activity can wax and wane during tumor evolution. The majority of the retrotransposition events are likely to be harmless “passenger” mutations. Hitchhiking through the tumor genome Retrotransposons are DNA repeat sequences that are constantly on the move. By poaching certain cellular enzymes, they copy and insert themselves at new sites in the genome. Sometimes they carry along adjacent DNA sequences, a process called 3′ transduction. Tubio et al. found that 3′ transduction is a common event in human tumors. Because this process can scatter genes and regulatory sequences across the genome, it may represent yet another mechanism by which tumor cells acquire new mutations that help them survive and grow. Science, this issue p. 10.1126/science.1251343 Tumor genomes are peppered with mobile repeat sequences that carry along adjacent DNA when they insert into new genomic sites. Long interspersed nuclear element–1 (L1) retrotransposons are mobile repetitive elements that are abundant in the human genome. L1 elements propagate through RNA intermediates. In the germ line, neighboring, nonrepetitive sequences are occasionally mobilized by the L1 machinery, a process called 3′ transduction. Because 3′ transductions are potentially mutagenic, we explored the extent to which they occur somatically during tumorigenesis. Studying cancer genomes from 244 patients, we found that tumors from 53% of the patients had somatic retrotranspositions, of which 24% were 3′ transductions. Fingerprinting of donor L1s revealed that a handful of source L1 elements in a tumor can spawn from tens to hundreds of 3′ transductions, which can themselves seed further retrotranspositions. The activity of individual L1 elements fluctuated during tumor evolution and correlated with L1 promoter hypomethylation. The 3′ transductions disseminated genes, exons, and regulatory elements to new locations, most often to heterochromatic regions of the genome.

Origins and functional consequences of somatic mitochondrial DNA mutations in human cancer.

Recent sequencing studies have extensively explored the somatic alterations present in the nuclear genomes of cancers. Although mitochondria control energy metabolism and apoptosis, the origins and impact of cancer-associated mutations in mtDNA are unclear. In this study, we analyzed somatic alterations in mtDNA from 1675 tumors. We identified 1907 somatic substitutions, which exhibited dramatic replicative strand bias, predominantly C > T and A > G on the mitochondrial heavy strand. This strand-asymmetric signature differs from those found in nuclear cancer genomes but matches the inferred germline process shaping primate mtDNA sequence content. A number of mtDNA mutations showed considerable heterogeneity across tumor types. Missense mutations were selectively neutral and often gradually drifted towards homoplasmy over time. In contrast, mutations resulting in protein truncation undergo negative selection and were almost exclusively heteroplasmic. Our findings indicate that the endogenous mutational mechanism has far greater impact than any other external mutagens in mitochondria and is fundamentally linked to mtDNA replication. DOI: http://dx.doi.org/10.7554/eLife.02935.001

Combined hereditary and somatic mutations of replication error repair genes result in rapid onset of ultra-hypermutated cancers

DNA replication−associated mutations are repaired by two components: polymerase proofreading and mismatch repair. The mutation consequences of disruption to both repair components in humans are not well studied. We sequenced cancer genomes from children with inherited biallelic mismatch repair deficiency (bMMRD). High-grade bMMRD brain tumors exhibited massive numbers of substitution mutations (>250/Mb), which was greater than all childhood and most cancers (>7,000 analyzed). All ultra-hypermutated bMMRD cancers acquired early somatic driver mutations in DNA polymerase ɛ or δ. The ensuing mutation signatures and numbers are unique and diagnostic of childhood germ-line bMMRD (P < 10−13). Sequential tumor biopsy analysis revealed that bMMRD/polymerase-mutant cancers rapidly amass an excess of simultaneous mutations (∼600 mutations/cell division), reaching but not exceeding ∼20,000 exonic mutations in <6 months. This implies a threshold compatible with cancer-cell survival. We suggest a new mechanism of cancer progression in which mutations develop in a rapid burst after ablation of replication repair.

Diagnostic value of H3F3A mutations in giant cell tumour of bone compared to osteoclast-rich mimics.

Driver mutations in the two histone 3.3 (H3.3) genes, H3F3A and H3F3B, were recently identified by whole genome sequencing in 95% of chondroblastoma (CB) and by targeted gene sequencing in 92% of giant cell tumour of bone (GCT). Given the high prevalence of these driver mutations, it may be possible to utilise these alterations as diagnostic adjuncts in clinical practice. Here, we explored the spectrum of H3.3 mutations in a wide range and large number of bone tumours (n = 412) to determine if these alterations could be used to distinguish GCT from other osteoclast‐rich tumours such as aneurysmal bone cyst, nonossifying fibroma, giant cell granuloma, and osteoclast‐rich malignant bone tumours and others. In addition, we explored the driver landscape of GCT through whole genome, exome and targeted sequencing (14 gene panel). We found that H3.3 mutations, namely mutations of glycine 34 in H3F3A, occur in 96% of GCT. We did not find additional driver mutations in GCT, including mutations in IDH1, IDH2, USP6, TP53. The genomes of GCT exhibited few somatic mutations, akin to the picture seen in CB. Overall our observations suggest that the presence of H3F3A p.Gly34 mutations does not entirely exclude malignancy in osteoclast‐rich tumours. However, H3F3A p.Gly34 mutations appear to be an almost essential feature of GCT that will aid pathological evaluation of bone tumours, especially when confronted with small needle core biopsies. In the absence of H3F3A p.Gly34 mutations, a diagnosis of GCT should be made with caution.

Frequent somatic transfer of mitochondrial DNA into the nuclear genome of human cancer cells

Mitochondrial genomes are separated from the nuclear genome for most of the cell cycle by the nuclear double membrane, intervening cytoplasm, and the mitochondrial double membrane. Despite these physical barriers, we show that somatically acquired mitochondrial-nuclear genome fusion sequences are present in cancer cells. Most occur in conjunction with intranuclear genomic rearrangements, and the features of the fusion fragments indicate that nonhomologous end joining and/or replication-dependent DNA double-strand break repair are the dominant mechanisms involved. Remarkably, mitochondrial-nuclear genome fusions occur at a similar rate per base pair of DNA as interchromosomal nuclear rearrangements, indicating the presence of a high frequency of contact between mitochondrial and nuclear DNA in some somatic cells. Transmission of mitochondrial DNA to the nuclear genome occurs in neoplastically transformed cells, but we do not exclude the possibility that some mitochondrial-nuclear DNA fusions observed in cancer occurred years earlier in normal somatic cells.

Dysfunction of phospholipase Cγ in immune disorders and cancer

The surge in genetic and genomic investigations over the past 5 years has resulted in many discoveries of causative variants relevant to disease pathophysiology. Although phospholipase C (PLC) enzymes have long been recognized as important components in intracellular signal transmission, it is only recently that this approach highlighted their role in disease development through gain-of-function mutations. In this review we describe the new findings that link the PLCγ family to immune disorders and cancer, and illustrate further efforts to elucidate the molecular mechanisms that underpin their dysfunction.

Fibroblastic growth factor receptor 1 amplification in osteosarcoma is associated with poor response to neo-adjuvant chemotherapy

Osteosarcoma, the most common primary bone sarcoma, is a genetically complex disease with no widely accepted biomarker to allow stratification of patients for treatment. After a recent report of one osteosarcoma cell line and one tumor exhibiting fibroblastic growth factor receptor 1 (FGFR1) gene amplification, the aim of this work was to assess the frequency of FGFR1 amplification in a larger cohort of osteosarcoma and to determine if this biomarker could be used for stratification of patients for treatment. About 352 osteosarcoma samples from 288 patients were analyzed for FGFR1 amplification by interphase fluorescence in situ hybridization. FGFR1 amplification was detected in 18.5% of patients whose tumors revealed a poor response to chemotherapy, and no patients whose tumors responded well to therapy harbored this genetic alteration. FGFR1 amplification is present disproportionately in the rarer histological variants of osteosarcoma. This study provides a rationale for inclusion of patients with osteosarcoma in clinical trials using FGFR kinase inhibitors.

SoupX removes ambient RNA contamination from droplet based single cell RNA sequencing data

Droplet based single cell RNA sequence analyses assume all acquired RNAs are endogenous to cells. However, any cell free RNAs contained within the input solution are also captured by these assays. This sequencing of cell free RNA constitutes a background contamination that has the potential to confound the correct biological interpretation of single cell transcriptomic data. Here, we demonstrate that contamination from this “soup” of cell free RNAs is ubiquitous, experiment specific in its composition and magnitude, and can lead to erroneous biological conclusions. We present a method, SoupX, for quantifying the extent of the contamination and estimating “background corrected”, cell expression profiles that can be integrated with existing downstream analysis tools. We apply this method to two data-sets and show that the application of this method reduces batch effects, strengthens cell-specific quality control and improves biological interpretation

Report from the 4th European Bone Sarcoma Networking meeting: focus on osteosarcoma

This report summarizes the proceedings of the 4th European Bone Sarcoma Networking Meeting, held in London, England, on 21 June 2017. The meeting brought together scientific and clinical researchers and representatives from sarcoma charities from 19 countries representing five networks across Europe, to present and discuss new developments on bone sarcoma. In view of the challenges is poses, the meeting focussed primarily on osteosarcoma with presentations on developments in our understanding of osteosarcoma genetics and immunology as well as results from preclinical investigations and discussion of recent and ongoing clinical trials. These include studies examining the efficacy of multi-targeted tyrosine kinase inhibitors and checkpoint inhibitors, as well as those with molecular profiling to stratify patients for specific therapies. Discussion was centred on generation of new hypotheses for collaborative biological and clinical investigations, the ultimate goal being to improve therapy and outcome in patients with bone sarcomas.