Identification of Germline Variants in Tumor Genomic Sequencing Analysis

This Correspondence relates to the article by Li et al (Standards and Guidelines for the Interpretation and Reporting of Sequence Variants in Cancer: A Joint Consensus Recommendation of the Association for Molecular Pathology, American Society of Clinical Oncology, and the College of American Pathologists. J Mol Diagn 2017, 19:4–23)

An Immunocompetent Mouse Model of HPV16(+) Head and Neck Squamous Cell Carcinoma

The incidence of human papilloma virus (HPV)-associated head and neck squamous cell carcinoma (HNSCC) is increasing and implicated in more than 60% of all oropharyngeal carcinomas (OPSCCs). Although whole-genome, transcriptome, and proteome analyses have identified altered signaling pathways in HPV-induced HNSCCs, additional tools are needed to investigate the unique pathobiology of OPSCC. Herein, bioinformatics analyses of human HPV(+) HNSCCs revealed that all tumors express full-length E6 and identified molecular subtypes based on relative E6 and E7 expression levels. To recapitulate the levels, stoichiometric ratios, and anatomic location of E6/E7 expression, we generated a genetically engineered mouse model whereby balanced expression of E6/E7 is directed to the oropharyngeal epithelium. The addition of a mutant PIK3CAE545K allele leads to the rapid development of pre-malignant lesions marked by immune cell accumulation, and a subset of these lesions progress to OPSCC. This mouse provides a faithful immunocompetent model for testing treatments and investigating mechanisms of immunosuppression

Exceptional chemotherapy response in metastatic colorectal cancer associated with hyper-indel-hypermutated cancer genome and comutation of POLD1 and MLH1

Purpose A 73-year-old woman with metastatic colon cancer experienced a complete response to chemotherapy with dose-intensified irinotecan that has been durable for 5 years.Wesequenced her tumor and germ line DNA and looked for similar patterns in publicly available genomic data from patients with colorectal cancer. Patients and Methods Tumor DNA was obtained from a biopsy before therapy, and germ line DNA was obtained from blood. Tumor and germline DNA were sequenced using a commercial panel with approximately 250 genes. Whole-genome amplification and exome sequencing were performed for POLE and POLD1. A POLD1 mutation was confirmed by Sanger sequencing. The somatic mutation and clinical annotation data files from the colon (n = 461) and rectal (n = 171) adenocarcinoma data sets were downloaded from The Cancer Genome Atlas data portal and analyzed for patterns of mutations and clinical outcomes in patients withPOLE- and/orPOLD1- mutated tumors. Results The pattern of alterations included APC biallelic inactivation and microsatellite instability high (MSI-H) phenotype, with somatic inactivation of MLH1 and hypermutation (estimated mutation rate > 200 per megabase). The extremely high mutation rate led us to investigate additional mechanisms for hypermutation, including loss of function of POLE. POLE was unaltered, but a related gene not typically associated with somatic mutation in colon cancer, POLD1, had a somatic mutation c.2171G > A[p.Gly724Glu]. Additionally, we noted that the high mutation rate was largely composed of dinucleotide deletions. A similar pattern of hypermutation (dinucleotide deletions, POLD1 mutations, MSI-H) was found in tumors from The Cancer Genome Atlas. Conclusion POLD1 mutation with associated MSI-H and hyper-indel-hypermutated cancer genome characterizes a previously unrecognized variant of colon cancer that was found in this patient with an exceptional response to chemotherapy

The association between copy number aberration, DNA methylation and gene expression in tumor samples

We systematically studied the association between somatic copy number aberration (SCNA), DNA methylation and gene expression using -omic data from The Cancer Genome Atlas (TCGA) on six cancer types: breast cancer, colon cancer, glioblastoma, leukemia, lower-grade glioma and prostate cancer. A major challenge for such integrated study is that the association between DNA methylation and gene expression is severely confounded by tumor purity and cell type composition, which are often unobserved and difficult to estimate. To overcome this challenge, we developed a method to remove confounding effects by calculating the principal components that span the space of the latent factors. Another intriguing findings of our study is that there could be both positive and negative associations between SCNA and DNA methylation, while the CpGs with negative/positive associations with SCNA are often located around CpG islands/ocean, respectively. A joint study of SCNA, DNA methylation, and gene expression suggest that SCNA often affect DNA methylation and gene expression independently. © The Author(s) 2018. Published by Oxford University Press on behalf of Nucleic Acids Research

SCISSOR: a framework for identifying structural changes in RNA transcripts

High-throughput sequencing protocols such as RNA-seq have made it possible to interrogate the sequence, structure and abundance of RNA transcripts at higher resolution than previous microarray and other molecular techniques. While many computational tools have been proposed for identifying mRNA variation through differential splicing/alternative exon usage, challenges in its analysis remain. Here, we propose a framework for unbiased and robust discovery of aberrant RNA transcript structures using short read sequencing data based on shape changes in an RNA-seq coverage profile. Shape changes in selecting sample outliers in RNA-seq, SCISSOR, is a series of procedures for transforming and normalizing base-level RNA sequencing coverage data in a transcript independent manner, followed by a statistical framework for its analysis (https://github.com/hyochoi/SCISSOR). The resulting high dimensional object is amenable to unsupervised screening of structural alterations across RNA-seq cohorts with nearly no assumption on the mutational mechanisms underlying abnormalities. This enables SCISSOR to independently recapture known variants such as splice site mutations in tumor suppressor genes as well as novel variants that are previously unrecognized or difficult to identify by any existing methods including recurrent alternate transcription start sites and recurrent complex deletions in 3′ UTRs

Genome-wide DNA methylation analysis of KRAS mutant cell lines

Oncogenic RAS mutations are associated with DNA methylation changes that alter gene expression to drive cancer. Recent studies suggest that DNA methylation changes may be stochastic in nature, while other groups propose distinct signaling pathways responsible for aberrant methylation. Better understanding of DNA methylation events associated with oncogenic KRAS expression could enhance therapeutic approaches. Here we analyzed the basal CpG methylation of 11 KRAS-mutant and dependent pancreatic cancer cell lines and observed strikingly similar methylation patterns. KRAS knockdown resulted in unique methylation changes with limited overlap between each cell line. In KRAS-mutant Pa16C pancreatic cancer cells, while KRAS knockdown resulted in over 8,000 differentially methylated (DM) CpGs, treatment with the ERK1/2-selective inhibitor SCH772984 showed less than 40 DM CpGs, suggesting that ERK is not a broadly active driver of KRAS-associated DNA methylation. KRAS G12V overexpression in an isogenic lung model reveals >50,600 DM CpGs compared to non-transformed controls. In lung and pancreatic cells, gene ontology analyses of DM promoters show an enrichment for genes involved in differentiation and development. Taken all together, KRAS-mediated DNA methylation are stochastic and independent of canonical downstream effector signaling. These epigenetically altered genes associated with KRAS expression could represent potential therapeutic targets in KRAS-driven cancer

Rapid clearance profile of plasma circulating tumor HPV type 16 DNA during chemoradiotherapy correlates with disease control in HPV-associated oropharyngeal cancer

Purpose: To identify a profile of circulating tumor human papilloma virus (HPV) DNA (ctHPVDNA) clearance kinetics that is associated with disease control after chemoradiotherapy (CRT) for HPV-associated oropharyngeal squamous cell carcinoma (OPSCC). Experimental Design: A multi-institutional prospective biomarker trial was conducted in 103 patients with (i) p16- positive OPSCC, (ii) M0 disease, and (iii) receipt of definitive CRT. Blood specimens were collected at baseline, weekly during CRT, and at follow-up visits. Optimized multianalyte digital PCR assays were used to quantify ctHPVDNA (types 16/18/31/33/35) in plasma. A control cohort of 55 healthy volunteers and 60 patients with non-HPV-associated malignancy was also analyzed. Results: Baseline plasma ctHPVDNA had high specificity (97%) and high sensitivity (89%) for detecting newly diagnosed HPV-associated OPSCC. Pretreatment ctHPV16DNA copy number correlated with disease burden, tumor HPV copy number, and HPV integration status. We define a ctHPV16DNA favorable clearance profile as having high baseline copy number (>200 copies/mL) and >95% clearance of ctHPV16DNA by day 28 of CRT. Nineteen of 67 evaluable patients had a ctHPV16DNA favorable clearance profile, and none had persistent or recurrent regional disease after CRT. In contrast, patients with adverse clinical risk factors (T4 or >10 pack years) and an unfavorable ctHPV16DNA clearance profile had a 35% actuarial rate of persistent or recurrent regional disease after CRT (P = 0.0049). Conclusions: A rapid clearance profile of ctHPVDNA may predict likelihood of disease control in patients with HPVassociated OPSCC patients treated with definitive CRT and may be useful in selecting patients for deintensified therapy

PIK3CA mutation in HPV-associated OPSCC patients receiving deintensified chemoradiation

PIK3CA is the most frequently mutated gene in human papillomavirus (HPV)-associated oropharyngeal squamous cell carcinoma (OPSCC). Prognostic implications of such mutations remain unknown. We sought to elucidate the clinical significance of PIK3CA mutations in HPV-associated OPSCC patients treated with definitive chemoradiation (CRT). Seventyseven patients with HPV-associated OPSCC were enrolled on two phase II clinical trials of deintensified CRT (60 Gy intensitymodulated radiotherapy with concurrent weekly cisplatin). Targeted next-generation sequencing was performed. Of the 77 patients, nine had disease recurrence (two regional, four distant, three regional and distant). Thirty-four patients had mutation( s) identified; 16 had PIK3CA mutations. Patients with wild-type-PIK3CA had statistically significantly higher 3-year disease-free survival than PIK3CA-mutant patients (93.4%, 95% confidence interval [CI] = 85.0% to 99.9% vs 68.8%, 95% CI = 26.7% to 89.8%; P=.004). On multivariate analysis, PIK3CA mutation was the only variable statistically significantly associated with disease recurrence (hazard ratio = 5.71, 95% CI = 1.53 to 21.3; P=.01). PIK3CA mutation is associated with worse diseasefree survival in a prospective cohort of newly diagnosed HPV-associated OPSCC patients treated with deintensified CRT

Enhancing Next-Generation Sequencing-Guided Cancer Care Through Cognitive Computing

Background: Using next-generation sequencing (NGS) to guide cancer therapy has created challenges in analyzing and reporting large volumes of genomic data to patients and caregivers. Specifically, providing current, accurate information on newly approved therapies and open clinical trials requires considerable manual curation performed mainly by human “molecular tumor boards” (MTBs). The purpose of this study was to determine the utility of cognitive computing as performed by Watson for Genomics (WfG) compared with a human MTB. Materials and Methods: One thousand eighteen patient cases that previously underwent targeted exon sequencing at the University of North Carolina (UNC) and subsequent analysis by the UNCseq informatics pipeline and the UNC MTB between November 7, 2011, and May 12, 2015, were analyzed with WfG, a cognitive computing technology for genomic analysis. Results: Using a WfG-curated actionable gene list, we identified additional genomic events of potential significance (not discovered by traditional MTB curation) in 323 (32%) patients. The majority of these additional genomic events were considered actionable based upon their ability to qualify patients for biomarker-selected clinical trials. Indeed, the opening of a relevant clinical trial within 1 month prior to WfG analysis provided the rationale for identification of a new actionable event in nearly a quarter of the 323 patients. This automated analysis took <3 minutes per case. Conclusion: These results demonstrate that the interpretation and actionability of somatic NGS results are evolving too rapidly to rely solely on human curation. Molecular tumor boards empowered by cognitive computing could potentially improve patient care by providing a rapid, comprehensive approach for data analysis and consideration of up-to-date availability of clinical trials. Implications for Practice: The results of this study demonstrate that the interpretation and actionability of somatic next-generation sequencing results are evolving too rapidly to rely solely on human curation. Molecular tumor boards empowered by cognitive computing can significantly improve patient care by providing a fast, cost-effective, and comprehensive approach for data analysis in the delivery of precision medicine. Patients and physicians who are considering enrollment in clinical trials may benefit from the support of such tools applied to genomic data

Use of mobile device technology to continuously collect patient-reported symptoms during radiation therapy for head and neck cancer: A prospective feasibility study

Purpose: Accurate assessment of toxicity allows for timely delivery of supportive measures during radiation therapy for head and neck cancer. The current paradigm requires weekly evaluation of patients by a provider. The purpose of this study is to evaluate the feasibility of monitoring patient reported symptoms via mobile devices. Methods and materials: We developed a mobile application for patients to report symptoms in 5 domains using validated questions. Patients were asked to report symptoms using a mobile device once daily during treatment or more often as needed. Clinicians reviewed patient-reported symptoms during weekly symptom management visits and patients completed surveys regarding perceptions of the utility of the mobile application. The primary outcome measure was patient compliance with mobile device reporting. Compliance is defined as number of days with a symptom report divided by number of days on study. Results: There were 921 symptom reports collected from 22 patients during treatment. Median reporting compliance was 71% (interquartile range, 45%-80%). Median number of reports submitted per patient was 34 (interquartile range, 21-53). Median number of reports submitted by patients per week was similar throughout radiation therapy and there was significant reporting during nonclinic hours. Patients reported high satisfaction with the use of mobile devices to report symptoms. Conclusions: A substantial percentage of patients used mobile devices to continuously report symptoms throughout a course of radiation therapy for head and neck cancer. Future studies should evaluate the impact of mobile device symptom reporting on improving patient outcomes