Jiabin Tang

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

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.

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.

Cerebrospinal fluid-derived circulating tumour DNA better represents the genomic alterations of brain tumours than plasma

Cell-free circulating tumour DNA (ctDNA) in plasma has been shown to be informative of the genomic alterations present in tumours and has been used to monitor tumour progression and response to treatments. However, patients with brain tumours do not present with or present with low amounts of ctDNA in plasma precluding the genomic characterization of brain cancer through plasma ctDNA. Here we show that ctDNA derived from central nervous system tumours is more abundantly present in the cerebrospinal fluid (CSF) than in plasma. Massively parallel sequencing of CSF ctDNA more comprehensively characterizes the genomic alterations of brain tumours than plasma, allowing the identification of actionable brain tumour somatic mutations. We show that CSF ctDNA levels longitudinally fluctuate in time and follow the changes in brain tumour burden providing biomarkers to monitor brain malignancies. Moreover, CSF ctDNA is shown to facilitate and complement the diagnosis of leptomeningeal carcinomatosis.

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.

Evaluating Cancer of the Central Nervous System Through Next-Generation Sequencing of Cerebrospinal Fluid

Purpose Cancer spread to the central nervous system (CNS) often is diagnosed late and is unresponsive to therapy. Mechanisms of tumor dissemination and evolution within the CNS are largely unknown because of limited access to tumor tissue. Materials and Methods We sequenced 341 cancer-associated genes in cell-free DNA from cerebrospinal fluid (CSF) obtained through routine lumbar puncture in 53 patients with suspected or known CNS involvement by cancer. Results We detected high-confidence somatic alterations in 63% (20 of 32) of patients with CNS metastases of solid tumors, 50% (six of 12) of patients with primary brain tumors, and 0% (zero of nine) of patients without CNS involvement by cancer. Several patients with tumor progression in the CNS during therapy with inhibitors of oncogenic kinases harbored mutations in the kinase target or kinase bypass pathways. In patients with glioma, the most common malignant primary brain tumor in adults, examination of cell-free DNA uncovered patterns of tumor evolution, including temozolomide-associated mutations. Conclusion The study shows that CSF harbors clinically relevant genomic alterations in patients with CNS cancers and should be considered for liquid biopsies to monitor tumor evolution in the CNS.

Plasma DNA-Based Molecular Diagnosis, Prognostication, and Monitoring of Patients With EWSR1 Fusion-Positive Sarcomas

Purpose Ewing Sarcoma (ES) and Desmoplastic Small Round Cell Tumors (DSRCT) are aggressive sarcomas molecularly characterized by EWSR1 gene fusions. As pathognomonic genomic events in these respective tumor types, EWSR1 fusions represent robust potential biomarkers for disease monitoring. Patients and Methods To investigate the feasibility of identifying EWSR1 fusions in plasma derived cell-free DNA (cfDNA) from ES and DSRCT patients, we evaluated two complementary approaches in samples from 17 patients with radiographic evidence of disease. The first approach involved identification of patient-specific genomic EWSR1 fusion breakpoints in formalin-fixed, paraffin-embedded tumor DNA using a broad, hybridization capture-based next generation sequencing (NGS) panel, followed by design of patient-specific droplet digital PCR (ddPCR) assays for plasma cfDNA interrogation . The second approach employed a disease-tailored targeted hybridization capture-based NGS panel applied directly to cfDNA which included EWSR1 as well as several other genes with potential prognostic utility. Results EWSR1 fusions were identified in 11/11 (100%) ES and 5/6 (83%) DSRCT samples by ddPCR, while 10/11 (91%) and 4/6 (67%) were identified by NGS. The ddPCR approach had higher sensitivity, ranging between 0.009-0.018% sensitivity. However, the hybrid capture-based NGS assay identified the precise fusion breakpoints in the majority of cfDNA samples, as well as mutations in TP53 and STAG2, two other recurrent, clinically significant alterations in ES, all without prior knowledge of the tumor sequencing results. Conclusion These results provide a compelling rationale for an integrated approach utilizing both NGS and ddPCR for plasma cfDNA-based biomarker evaluations in prospective cooperative group studies.

A Pan-Cancer Compendium of Genes Deregulated by Somatic Genomic Rearrangement across More Than 1,400 Cases

A systematic cataloging of genes affected by genomic rearrangement, using multiple patient cohorts and cancer types, can provide insight into cancer-relevant alterations outside of exomes. By integrative analysis of whole-genome sequencing (predominantly low pass) and gene expression data from 1,448 cancers involving 18 histopathological types in The Cancer Genome Atlas, we identified hundreds of genes for which the nearby presence (within 100 kb) of a somatic structural variant (SV) breakpoint is associated with altered expression. While genomic rearrangements are associated with widespread copy-number alteration (CNA) patterns, approximately 1,100 genes-including overexpressed cancer driver genes (e.g., TERT, ERBB2, CDK12, CDK4) and underexpressed tumor suppressors (e.g., TP53, RB1, PTEN, STK11)-show SV-associated deregulation independent of CNA. SVs associated with the disruption of topologically associated domains, enhancer hijacking, or fusion transcripts are implicated in gene upregulation. For cancer-relevant pathways, SVs considerably expand our understanding of how genes are affected beyond point mutation or CNA.

DEFINING AGGRESSIVE PROSTATE CANCER USING A 12 GENE MODEL

The critical clinical question in prostate cancer research is: How do we develop means of distinguishing aggressive disease from indolent disease? Using a combination of proteomic and expression array data, we identified a set of 36 genes with concordant dysregulation of protein products that could be evaluated in situ by quantitative immunohistochemistry. Another five prostate cancer biomarkers were included using linear discriminant analysis, we determined that the optimal model used to predict prostate cancer progression consisted of 12 proteins. Using a separate patient population, transcriptional levels of the 12 genes encoding for these proteins predicted prostate-specific antigen failure in 79 men following surgery for clinically localized prostate cancer (P = .0015). This study demonstrates that cross-platform models can lead to predictive models with the possible advantage of being more robust through this selection process.