Sahil Seth

Intratumor heterogeneity in localized lung adenocarcinomas delineated by multiregion sequencing

Cancers are composed of populations of cells with distinct molecular and phenotypic features, a phenomenon termed intratumor heterogeneity (ITH). ITH in lung cancers has not been well studied. We applied multiregion whole-exome sequencing (WES) on 11 localized lung adenocarcinomas. All tumors showed clear evidence of ITH. On average, 76% of all mutations and 20 out of 21 known cancer gene mutations were identified in all regions of individual tumors, which suggested that single-region sequencing may be adequate to identify the majority of known cancer gene mutations in localized lung adenocarcinomas. With a median follow-up of 21 months after surgery, three patients have relapsed, and all three patients had significantly larger fractions of subclonal mutations in their primary tumors than patients without relapse. These data indicate that a larger subclonal mutation fraction may be associated with increased likelihood of postsurgical relapse in patients with localized lung adenocarcinomas. Different mutations are present in different regions of any given lung cancer, and their pattern may predict patient relapse. [Also see Perspective by Govindan] Space, time, and the lung cancer genome Lung cancer poses a formidable challenge to clinical oncologists. It is often detected at a late stage, and most therapies work for only a short time before the tumors resume their relentless growth. Two independent analyses of the human lung cancer genome may help explain why this disease is so resilient (see the Perspective by Govindan). Rather than take a single “snapshot” of the cancer genome, de Bruin et al. and Zhang et al. identified genomic alterations in spatially distinct regions of single lung tumors and used this information to infer the tumors evolutionary history. Each tumor showed tremendous spatial and temporal diversity in its mutational profiles. Thus, the efficacy of drugs may be short-lived because they destroy only a portion of the tumor. Science, this issue p. 251, p. 256; see also p. 169

Oncogene ablation-resistant pancreatic cancer cells depend on mitochondrial function

Pancreatic ductal adenocarcinoma (PDAC) is one of the deadliest cancers in western countries, with a median survival of 6 months and an extremely low percentage of long-term surviving patients. KRAS mutations are known to be a driver event of PDAC, but targeting mutant KRAS has proved challenging. Targeting oncogene-driven signalling pathways is a clinically validated approach for several devastating diseases. Still, despite marked tumour shrinkage, the frequency of relapse indicates that a fraction of tumour cells survives shut down of oncogenic signalling. Here we explore the role of mutant KRAS in PDAC maintenance using a recently developed inducible mouse model of mutated Kras (KrasG12D, herein KRas) in a p53LoxP/WT background. We demonstrate that a subpopulation of dormant tumour cells surviving oncogene ablation (surviving cells) and responsible for tumour relapse has features of cancer stem cells and relies on oxidative phosphorylation for survival. Transcriptomic and metabolic analyses of surviving cells reveal prominent expression of genes governing mitochondrial function, autophagy and lysosome activity, as well as a strong reliance on mitochondrial respiration and a decreased dependence on glycolysis for cellular energetics. Accordingly, surviving cells show high sensitivity to oxidative phosphorylation inhibitors, which can inhibit tumour recurrence. Our integrated analyses illuminate a therapeutic strategy of combined targeting of the KRAS pathway and mitochondrial respiration to manage pancreatic cancer.

Identification of Double-stranded Genomic DNA Spanning All Chromosomes with Mutated KRAS and p53 DNA in the Serum Exosomes of Patients with Pancreatic Cancer

Background: Exosomes are small vesicles in the tumor microenvironment containing nucleic acids and proteins with the capacity to influence cell behavior. Results: Exosomes contain double-stranded genomic DNA. Conclusion: Exosomes have the capacity to carry and transport genomic DNA spanning all chromosomes with KRAS and p53 mutations. Significance: Exosomes can aid in identifying genomic mutations in patients with pancreatic cancer. Exosomes are small vesicles (50–150 nm) of endocytic origin that are released by many different cell types. Exosomes in the tumor microenvironment may play a key role in facilitating cell-cell communication. Exosomes are reported to predominantly contain RNA and proteins. In this study, we investigated whether exosomes from pancreatic cancer cells and serum from patients with pancreatic ductal adenocarcinoma contain genomic DNA. Our results provide evidence that exosomes contain >10-kb fragments of double-stranded genomic DNA. Mutations in KRAS and p53 can be detected using genomic DNA from exosomes derived from pancreatic cancer cell lines and serum from patients with pancreatic cancer. In addition, using whole genome sequencing, we demonstrate that serum exosomes from patients with pancreatic cancer contain genomic DNA spanning all chromosomes. These results indicate that serum-derived exosomes can be used to determine genomic DNA mutations for cancer prediction, treatment, and therapy resistance.

Mutational Landscape of Aggressive Cutaneous Squamous Cell Carcinoma

Purpose: Aggressive cutaneous squamous cell carcinoma (cSCC) is often a disfiguring and lethal disease. Very little is currently known about the mutations that drive aggressive cSCC. Experimental Design: Whole-exome sequencing was performed on 39 cases of aggressive cSCC to identify driver genes and novel therapeutic targets. Significantly, mutated genes were identified with MutSig or complementary methods developed to specifically identify candidate tumor suppressors based upon their inactivating mutation bias. Results: Despite the very high-mutational background caused by UV exposure, 23 candidate drivers were identified, including the well-known cancer-associated genes TP53, CDKN2A, NOTCH1, AJUBA, HRAS, CASP8, FAT1, and KMT2C (MLL3). Three novel candidate tumor suppressors with putative links to cancer or differentiation, NOTCH2, PARD3, and RASA1, were also identified as possible drivers in cSCC. KMT2C mutations were associated with poor outcome and increased bone invasion. Conclusions: The mutational spectrum of cSCC is similar to that of head and neck squamous cell carcinoma and dominated by tumor-suppressor genes. These results improve the foundation for understanding this disease and should aid in identifying and treating aggressive cSCC. Clin Cancer Res; 20(24); 6582–92. ©2014 AACR.

Characterization of HPV and host genome interactions in primary head and neck cancers

Significance A significant proportion of head and neck cancer is driven by human papillomavirus (HPV) infection, and the expression of viral oncogenes is involved in the development of these tumors. However, the role of HPV integration in primary tumors beyond increasing the expression of viral oncoproteins is not understood. Here, we describe how HPV integration impacts the host genome by amplification of oncogenes and disruption of tumor suppressors as well as driving inter- and intrachromosomal rearrangements. Tumors that do and do not have HPV integrants display distinct gene expression profiles and DNA methylation patterns, which further support the view that the mechanisms by which tumors with integrated and nonintegrated HPV arise are distinct. Previous studies have established that a subset of head and neck tumors contains human papillomavirus (HPV) sequences and that HPV-driven head and neck cancers display distinct biological and clinical features. HPV is known to drive cancer by the actions of the E6 and E7 oncoproteins, but the molecular architecture of HPV infection and its interaction with the host genome in head and neck cancers have not been comprehensively described. We profiled a cohort of 279 head and neck cancers with next generation RNA and DNA sequencing and show that 35 (12.5%) tumors displayed evidence of high-risk HPV types 16, 33, or 35. Twenty-five cases had integration of the viral genome into one or more locations in the human genome with statistical enrichment for genic regions. Integrations had a marked impact on the human genome and were associated with alterations in DNA copy number, mRNA transcript abundance and splicing, and both inter- and intrachromosomal rearrangements. Many of these events involved genes with documented roles in cancer. Cancers with integrated vs. nonintegrated HPV displayed different patterns of DNA methylation and both human and viral gene expressions. Together, these data provide insight into the mechanisms by which HPV interacts with the human genome beyond expression of viral oncoproteins and suggest that specific integration events are an integral component of viral oncogenesis.

Recurrent PTPRB and PLCG1 mutations in angiosarcoma

Angiosarcoma is an aggressive malignancy that arises spontaneously or secondarily to ionizing radiation or chronic lymphoedema. Previous work has identified aberrant angiogenesis, including occasional somatic mutations in angiogenesis signaling genes, as a key driver of angiosarcoma. Here we employed whole-genome, whole-exome and targeted sequencing to study the somatic changes underpinning primary and secondary angiosarcoma. We identified recurrent mutations in two genes, PTPRB and PLCG1, which are intimately linked to angiogenesis. The endothelial phosphatase PTPRB, a negative regulator of vascular growth factor tyrosine kinases, harbored predominantly truncating mutations in 10 of 39 tumors (26%). PLCG1, a signal transducer of tyrosine kinases, encoded a recurrent, likely activating p.Arg707Gln missense variant in 3 of 34 cases (9%). Overall, 15 of 39 tumors (38%) harbored at least one driver mutation in angiogenesis signaling genes. Our findings inform and reinforce current therapeutic efforts to target angiogenesis signaling in angiosarcoma.

A comprehensive assessment of somatic mutation detection in cancer using whole-genome sequencing

As whole-genome sequencing for cancer genome analysis becomes a clinical tool, a full understanding of the variables affecting sequencing analysis output is required. Here using tumour-normal sample pairs from two different types of cancer, chronic lymphocytic leukaemia and medulloblastoma, we conduct a benchmarking exercise within the context of the International Cancer Genome Consortium. We compare sequencing methods, analysis pipelines and validation methods. We show that using PCR-free methods and increasing sequencing depth to ∼100 × shows benefits, as long as the tumour:control coverage ratio remains balanced. We observe widely varying mutation call rates and low concordance among analysis pipelines, reflecting the artefact-prone nature of the raw data and lack of standards for dealing with the artefacts. However, we show that, using the benchmark mutation set we have created, many issues are in fact easy to remedy and have an immediate positive impact on mutation detection accuracy.

Systematic analysis of telomere length and somatic alterations in 31 cancer types

Cancer cells survive cellular crisis through telomere maintenance mechanisms. We report telomere lengths in 18,430 samples, including tumors and non-neoplastic samples, across 31 cancer types. Telomeres were shorter in tumors than in normal tissues and longer in sarcomas and gliomas than in other cancers. Among 6,835 cancers, 73% expressed telomerase reverse transcriptase (TERT), which was associated with TERT point mutations, rearrangements, DNA amplifications and transcript fusions and predictive of telomerase activity. TERT promoter methylation provided an additional deregulatory TERT expression mechanism. Five percent of cases, characterized by undetectable TERT expression and alterations in ATRX or DAXX, demonstrated elongated telomeres and increased telomeric repeat–containing RNA (TERRA). The remaining 22% of tumors neither expressed TERT nor harbored alterations in ATRX or DAXX. In this group, telomere length positively correlated with TP53 and RB1 mutations. Our analysis integrates TERT abnormalities, telomerase activity and genomic alterations with telomere length in cancer.

Co-clinical assessment identifies patterns of BRAF inhibitor resistance in melanoma

Multiple mechanisms have been described that confer BRAF inhibitor resistance to melanomas, yet the basis of this resistance remains undefined in a sizable portion of patient samples. Here, we characterized samples from a set of patients with melanoma that included individuals at baseline diagnosis, on BRAF inhibitor treatment, and with resistant tumors at both the protein and RNA levels. Using RNA and DNA sequencing, we identified known resistance-conferring mutations in 50% (6 of 12) of the resistant samples. In parallel, targeted proteomic analysis by protein array categorized the resistant samples into 3 stable groups, 2 of which were characterized by reactivation of MAPK signaling to different levels and 1 that was MAPK independent. The molecular relevance of these classifications identified in patients was supported by both mutation data and the similarity of resistance patterns that emerged during a co-clinical trial in a genetically engineered mouse (GEM) model of melanoma that recapitulates the development of BRAF inhibitor resistance. Additionally, we defined candidate biomarkers in pre- and early-treatment patient samples that have potential for predicting clinical responses. On the basis of these observations, we suggest that BRAF inhibitor-resistant melanomas can be actionably classified using protein expression patterns, even without identification of the underlying genetic alteration.

Genomic heterogeneity of multiple synchronous lung cancer

Multiple synchronous lung cancers (MSLCs) present a clinical dilemma as to whether individual tumours represent intrapulmonary metastases or independent tumours. In this study we analyse genomic profiles of 15 lung adenocarcinomas and one regional lymph node metastasis from 6 patients with MSLC. All 15 lung tumours demonstrate distinct genomic profiles, suggesting all are independent primary tumours, which are consistent with comprehensive histopathological assessment in 5 of the 6 patients. Lung tumours of the same individuals are no more similar to each other than are lung adenocarcinomas of different patients from TCGA cohort matched for tumour size and smoking status. Several known cancer-associated genes have different mutations in different tumours from the same patients. These findings suggest that in the context of identical constitutional genetic background and environmental exposure, different lung cancers in the same individual may have distinct genomic profiles and can be driven by distinct molecular events.