The somatic genomic landscape of glioblastoma.

We describe the landscape of somatic genomic alterations based on multidimensional and comprehensive characterization of more than 500 glioblastoma tumors (GBMs). We identify several novel mutated genes as well as complex rearrangements of signature receptors, including EGFR and PDGFRA. TERT promoter mutations are shown to correlate with elevated mRNA expression, supporting a role in telomerase reactivation. Correlative analyses confirm that the survival advantage of the proneural subtype is conferred by the G-CIMP phenotype, and MGMT DNA methylation may be a predictive biomarker for treatment response only in classical subtype GBM. Integrative analysis of genomic and proteomic profiles challenges the notion of therapeutic inhibition of a pathway as an alternative to inhibition of the target itself. These data will facilitate the discovery of therapeutic and diagnostic target candidates, the validation of research and clinical observations and the generation of unanticipated hypotheses that can advance our molecular understanding of this lethal cancer.

Author Correction: LKB1 loss links serine metabolism to DNA methylation and tumorigenesis.

An Amendment to this paper has been published and can be accessed via a link at the top of the paper.

Diagnostic accuracy of HPV16 early antigen serology for HPV-driven oropharyngeal cancer is independent of age and sex.

A growing proportion of head and neck cancer (HNC), especially oropharyngeal cancer (OPC), is caused by human papillomavirus (HPV). There are several markers for HPV-driven HNC, one being HPV early antigen serology. We aimed to investigate the diagnostic accuracy of HPV serology and its performance across patient characteristics. Data from the VOYAGER consortium was used, which comprises five studies on HNC from North America and Europe. Diagnostic accuracy, that is, sensitivity, specificity, Cohen's kappa and correctly classified proportions of HPV16 E6 serology, was assessed for OPC and other HNC using p16INK4a immunohistochemistry (p16), HPV in situ hybridization (ISH) and HPV PCR as reference methods. Stratified analyses were performed for variables including age, sex, smoking and alcohol use, to test the robustness of diagnostic accuracy. A risk-factor analysis based on serology was conducted, comparing HPV-driven to non-HPV-driven OPC. Overall, HPV serology had a sensitivity of 86.8% (95% CI 85.1-88.3) and specificity of 91.2% (95% CI 88.6-93.4) for HPV-driven OPC using p16 as a reference method. In stratified analyses, diagnostic accuracy remained consistent across sex and different age groups. Sensitivity was lower for heavy smokers (77.7%), OPC without lymph node involvement (74.4%) and the ARCAGE study (66.7%), while specificity decreased for cases with <10 pack-years (72.1%). The risk-factor model included study, year of diagnosis, age, sex, BMI, alcohol use, pack-years, TNM-T and TNM-N stage. HPV serology is a robust biomarker for HPV-driven OPC, and its diagnostic accuracy is independent of age and sex. Future research is suggested on the influence of smoking on HPV antibody levels.

DiNAMIC.Duo: detecting somatic DNA copy number differences without a normal reference.

MOTIVATION: Somatic DNA copy number alterations (CNAs) arise in tumor tissue because of underlying genomic instability. Recurrent CNAs that occur in the same genomic region across multiple independent samples are of interest to researchers because they may contain genes that contribute to the cancer phenotype. However, differences in copy number states between cancers are also commonly of interest, for example when comparing tumors with distinct morphologies in the same anatomic location. Current methodologies are limited by their inability to perform direct comparisons of CNAs between tumor cohorts, and thus they cannot formally assess the statistical significance of observed copy number differences or identify regions of the genome where these differences occur. RESULTS: We introduce the DiNAMIC.Duo R package that can be used to identify recurrent CNAs in a single cohort or recurrent copy number differences between two cohorts, including when neither cohort is copy neutral. The package utilizes Python scripts for computational efficiency and provides functionality for producing figures and summary output files. AVAILABILITY AND IMPLEMENTATION: The DiNAMIC.Duo R package is available from CRAN at https://cran.r-project.org/web/packages/DiNAMIC.Duo/index.html. This article uses publicly available data from the Broad Institute TCGA Genome Data Analysis Center, https://doi.org/10.7908/C11G0KM9. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Sequencing of genes of drug response in tumor DNA and implications for precision medicine in cancer patients.

Tumor DNA sequencing is becoming standard-of-care for patient treatment decisions. We evaluated genotype concordance between tumor DNA and genomic DNA from blood and catalogued functional effects of somatic mutations in 21 drug response genes in 752 solid tumor patients. Using a threshold of 10% difference between tumor and blood DNA variant allele fraction (VAF), concordance for heterogenous genotype calls was 78% and increased to 97.5% using a 30% VAF threshold. Somatic mutations were observed in all 21 drug response genes, and 44% of patients had at least one somatic mutation in these genes. In tumor DNA, eight patients had a frameshift mutation in CYP2C8, which metabolizes taxanes. Overall, somatic copy number losses were more frequent than gains, including for CYP2C19 and CYP2D6 which had the most frequent copy number losses. However, copy number gains in TPMT were more than four times as common as losses. Seven % of patients had copy number gains in ABCB1, a multidrug resistance transporter of anti-cancer agents. These results demonstrate tumor-only DNA sequencing might not be reliable to call germline genotypes of drug response variants.

LKB1 loss links serine metabolism to DNA methylation and tumorigenesis.

Intermediary metabolism generates substrates for chromatin modification, enabling the potential coupling of metabolic and epigenetic states. Here we identify a network linking metabolic and epigenetic alterations that is central to oncogenic transformation downstream of the liver kinase B1 (LKB1, also known as STK11) tumour suppressor, an integrator of nutrient availability, metabolism and growth. By developing genetically engineered mouse models and primary pancreatic epithelial cells, and employing transcriptional, proteomics, and metabolic analyses, we find that oncogenic cooperation between LKB1 loss and KRAS activation is fuelled by pronounced mTOR-dependent induction of the serine-glycine-one-carbon pathway coupled to S-adenosylmethionine generation. At the same time, DNA methyltransferases are upregulated, leading to elevation in DNA methylation with particular enrichment at retrotransposon elements associated with their transcriptional silencing. Correspondingly, LKB1 deficiency sensitizes cells and tumours to inhibition of serine biosynthesis and DNA methylation. Thus, we define a hypermetabolic state that incites changes in the epigenetic landscape to support tumorigenic growth of LKB1-mutant cells, while resulting in potential therapeutic vulnerabilities.

Fbxw7 is a driver of uterine carcinosarcoma by promoting epithelial-mesenchymal transition.

Uterine carcinosarcoma is an aggressive variant of endometrial carcinoma characterized by unusual histologic features including discrete malignant epithelial and mesenchymal components (carcinoma and sarcoma). Recent studies have confirmed a monoclonal origin, and comprehensive genomic characterizations have identified mutations such as Tp53 and Pten However, the biological origins and specific combination of driver events underpinning uterine carcinosarcoma have remained mysterious. Here, we explored the role of the tumor suppressor Fbxw7 in endometrial cancer through defined genetic model systems. Inactivation of Fbxw7 and Pten resulted in the formation of precancerous lesions (endometrioid intraepithelial neoplasia) and well-differentiated endometrioid adenocarcinomas. Surprisingly, all adenocarcinomas eventually developed into definitive uterine carcinosarcomas with carcinomatous and sarcomatous elements including heterologous differentiation, yielding a faithful genetically engineered model of this cancer type. Genomic analysis showed that most tumors spontaneously acquired Trp53 mutations, pointing to a triad of pathways (p53, PI3K, and Fbxw7) as the critical combination underpinning uterine carcinosarcoma, and to Fbxw7 as a key driver of this enigmatic endometrial cancer type. Lineage tracing provided formal genetic proof that the uterine carcinosarcoma cell of origin is an endometrial epithelial cell that subsequently undergoes a prominent epithelial-mesenchymal transition underlying the attainment of a highly invasive phenotype specifically driven by Fbxw7.

A conditional mouse expressing an activating mutation in NRF2 displays hyperplasia of the upper gastrointestinal tract and decreased white adipose tissue.

Activation of the nuclear factor (erythroid-derived 2)-like 2 (NFE2L2 or NRF2) transcription factor is a critical and evolutionarily conserved cellular response to oxidative stress, metabolic stress, and xenobiotic insult. Deficiency of NRF2 results in hypersensitivity to a variety of stressors, whereas its aberrant activation contributes to several cancer types, most commonly squamous cell carcinomas of the esophagus, oral cavity, bladder, and lung. Between 10% and 35% of patients with squamous cell carcinomas display hyperactive NRF2 signaling, harboring activating mutations and copy number amplifications of the NFE2L2 oncogene or inactivating mutations or deletions of KEAP1 or CUL3, the proteins of which co-complex to ubiquitylate and degrade NRF2 protein. To better understand the role of NRF2 in tumorigenesis and more broadly in development, we engineered the endogenous Nfe2l2 genomic locus to create a conditional mutant LSL-Nrf2E79Q mouse model. The E79Q mutation, one of the most commonly observed NRF2-activating mutations in human squamous cancers, codes for a mutant protein that does not undergo KEAP1/CUL3-dependent degradation, resulting in its constitutive activity. Expression of NRF2 E79Q protein in keratin 14 (KRT14)-positive murine tissues resulted in hyperplasia of squamous cell tissues of the tongue, forestomach, and esophagus, a stunted body axis, decreased weight, and decreased visceral adipose depots. RNA-seq profiling and follow-up validation studies of cultured NRF2E79Q murine esophageal epithelial cells revealed known and novel NRF2-regulated transcriptional programs, including genes associated with squamous cell carcinoma (e.g. Myc), lipid and cellular metabolism (Hk2, Ppard), and growth factors (Areg, Bmp6, Vegfa). These data suggest that in addition to decreasing adipogenesis, KRT14-restricted NRF2 activation drives hyperplasia of the esophagus, forestomach, and tongue, but not formation of squamous cell carcinoma. © 2020 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Decreased overall survival in black patients with HPV-associated oropharyngeal cancer.

PURPOSE: Racial disparities for overall survival (OS) in head and neck cancer have been well described. However, the extent to which these disparities exist for HPV-associated oropharyngeal squamous cell carcinoma (OPSCC), and the contribution of demographic, clinical, and socioeconomic status (SES) variables, is unknown. MATERIALS AND METHODS: Patients were identified from the Carolina Head and Neck Cancer Epidemiology Study (CHANCE), a population-based study in North Carolina. Cox proportional hazards regression models were used to estimate hazard ratios (HR) and 95% confidence intervals (CI) for OS in black versus white patients with sequential adjustment sets. RESULTS: A total of 157 HPV-associated OPSCC patients were identified. Of these, 93% were white and 7% were black. Black patients with HPV-associated OPSCC were more likely to be younger, have an income <$20,000, live farther away from clinic where biopsy was performed, and have advanced T stage at diagnosis. Black patients had worse OS in the unadjusted analysis (HR 4.9, 95% CI 2.2-11.1, p < 0.0001). The racial disparity in OS slightly decreased when sequentially adjusting for demographic, clinical, and SES variables. However, HR for black race remained statistically elevated in the final adjustment set which controlled for age, sex, stage, smoking, alcohol use, and individual-level household income, insurance, and education level (HR 3.4, 95% CI 1.1-10.1, p = 0.028). CONCLUSION: This is the first population-based study that confirms persistence of racial disparities in HPV-associated OPSCC after controlling for demographic, clinical, and individual-level socioeconomic factors.

Mature results of a prospective study of deintensified chemoradiotherapy for low-risk human papillomavirus-associated oropharyngeal squamous cell carcinoma.

BACKGROUND: The purpose of the current study was to determine quality of life and tumor control from a prospective phase 2 clinical trial evaluating deintensified chemoradiotherapy for favorable risk, human papillomavirus (HPV)-associated oropharyngeal squamous cell carcinoma. METHODS: Patients with T0-T3, N0-N2c, M0, p16-positive disease and a minimal smoking history were treated with 60 grays of intensity-modulated radiotherapy with concurrent weekly intravenous cisplatin (30 mg/m2 ). The primary study endpoint was the pathologic complete response rate based on biopsy of the primary site and dissection of pretreatment positive lymph node regions. The pathologic complete response rate as previously reported was 86%. Herein, the authors report secondary endpoint measures of local control, regional control, cause-specific survival, distant metastasis-free survival, and overall survival, and patient-reported outcomes (European Organization for Research and Treatment of Cancer [EORTC] Quality of Life Questionnaire [EORTC QLQ-C30] and the Patient-Reported Outcomes version of Common Terminology Criteria for Adverse Events [PRO-CTCAE]). RESULTS: A total of 44 patients enrolled with a median follow-up of 36 months (88% with ≥2 years). The 3-year local control, regional control, cause-specific survival, distant metastasis-free survival, and overall survival rates were 100%, 100%, 100%, 100%, and 95%, respectively. The mean before and 3-year after EORTC QOL scores were: global: 80 of 78; swallowing: 11 of 11; dry mouth: 16 of 41; and sticky saliva: 6 of 29. The mean before and 3-year after PRO-CTCAE scores were: swallowing: 0.4 of 0.7; and dry mouth: 0.4 of 1.4. Approximately 39% of patients required a feeding tube (median duration, 15 weeks; none were permanent). There were no ≥grade 3 late adverse events reported. CONCLUSIONS: For patients with favorable-risk human papillomavirus-associated oropharyngeal squamous cell carcinoma, a substantially decreased intensity of therapy with 60 grays of intensity-modulated radiotherapy and weekly low-dose cisplatin produced better preservation of quality of life compared with standard therapies while maintaining excellent 3-year tumor control and survival. Cancer 2018;124:2347-54. © 2018 American Cancer Society.

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.

LCCC 1025: a phase II study of everolimus, trastuzumab, and vinorelbine to treat progressive HER2-positive breast cancer brain metastases.

PURPOSE: HER2 + breast cancer (BC) is an aggressive subtype with high rates of brain metastases (BCBM). Two-thirds of HER2 + BCBM demonstrate activation of the PI3K/mTOR pathway driving resistance to anti-HER2 therapy. This phase II study evaluated everolimus (E), a brain-permeable mTOR inhibitor, trastuzumab (T), and vinorelbine (V) in patients with HER2 + BCBM. PATIENTS AND METHODS: Eligible patients had progressive HER2 + BCBM. The primary endpoint was intracranial response rate (RR); secondary objectives were CNS clinical benefit rate (CBR), extracranial RR, time to progression (TTP), overall survival (OS), and targeted sequencing of tumors from enrolled patients. A two-stage design distinguished intracranial RR of 5% versus 20%. RESULTS: 32 patients were evaluable for toxicity, 26 for efficacy. Intracranial RR was 4% (1 PR). CNS CBR at 6 mos was 27%; at 3 mos 65%. Median intracranial TTP was 3.9 mos (95% CI 2.2-5). OS was 12.2 mos (95% CI 0.6-20.2). Grade 3-4 toxicities included neutropenia (41%), anemia (16%), and stomatitis (16%). Mutations in TP53 and PIK3CA were common in BCBM. Mutations in the PI3K/mTOR pathway were not associated with response. ERBB2 amplification was higher in BCBM compared to primary BC; ERBB2 amplification in the primary BC trended toward worse OS. CONCLUSION: While intracranial RR to ETV was low in HER2 + BCBM patients, one-third achieved CNS CBR; TTP/OS was similar to historical control. No new toxicity signals were observed. Further analysis of the genomic underpinnings of BCBM to identify tractable prognostic and/or predictive biomarkers is warranted. CLINICAL TRIAL: (NCT01305941).

A murine lung cancer co-clinical trial identifies genetic modifiers of therapeutic response.

Targeted therapies have demonstrated efficacy against specific subsets of molecularly defined cancers. Although most patients with lung cancer are stratified according to a single oncogenic driver, cancers harbouring identical activating genetic mutations show large variations in their responses to the same targeted therapy. The biology underlying this heterogeneity is not well understood, and the impact of co-existing genetic mutations, especially the loss of tumour suppressors, has not been fully explored. Here we use genetically engineered mouse models to conduct a 'co-clinical' trial that mirrors an ongoing human clinical trial in patients with KRAS-mutant lung cancers. This trial aims to determine if the MEK inhibitor selumetinib (AZD6244) increases the efficacy of docetaxel, a standard of care chemotherapy. Our studies demonstrate that concomitant loss of either p53 (also known as Tp53) or Lkb1 (also known as Stk11), two clinically relevant tumour suppressors, markedly impaired the response of Kras-mutant cancers to docetaxel monotherapy. We observed that the addition of selumetinib provided substantial benefit for mice with lung cancer caused by Kras and Kras and p53 mutations, but mice with Kras and Lkb1 mutations had primary resistance to this combination therapy. Pharmacodynamic studies, including positron-emission tomography (PET) and computed tomography (CT), identified biological markers in mice and patients that provide a rationale for the differential efficacy of these therapies in the different genotypes. These co-clinical results identify predictive genetic biomarkers that should be validated by interrogating samples from patients enrolled on the concurrent clinical trial. These studies also highlight the rationale for synchronous co-clinical trials, not only to anticipate the results of ongoing human clinical trials, but also to generate clinically relevant hypotheses that can inform the analysis and design of human studies.

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

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.

Cafeteria diet-induced obesity causes oxidative damage in white adipose.

Obesity continues to be one of the most prominent public health dilemmas in the world. The complex interaction among the varied causes of obesity makes it a particularly challenging problem to address. While typical high-fat purified diets successfully induce weight gain in rodents, we have described a more robust model of diet-induced obesity based on feeding rats a diet consisting of highly palatable, energy-dense human junk foods - the "cafeteria" diet (CAF, 45-53% kcal from fat). We previously reported that CAF-fed rats became hyperphagic, gained more weight, and developed more severe hyperinsulinemia, hyperglycemia, and glucose intolerance compared to the lard-based 45% kcal from fat high fat diet-fed group. In addition, the CAF diet-fed group displayed a higher degree of inflammation in adipose and liver, mitochondrial dysfunction, and an increased concentration of lipid-derived, pro-inflammatory mediators. Building upon our previous findings, we aimed to determine mechanisms that underlie physiologic findings in the CAF diet. We investigated the effect of CAF diet-induced obesity on adipose tissue specifically using expression arrays and immunohistochemistry. Genomic evidence indicated the CAF diet induced alterations in the white adipose gene transcriptome, with notable suppression of glutathione-related genes and pathways involved in mitigating oxidative stress. Immunohistochemical analysis indicated a doubling in adipose lipid peroxidation marker 4-HNE levels compared to rats that remained lean on control standard chow diet. Our data indicates that the CAF diet drives an increase in oxidative damage in white adipose tissue that may affect tissue homeostasis. Oxidative stress drives activation of inflammatory kinases that can perturb insulin signaling leading to glucose intolerance and diabetes.

Association between differential gene expression and body mass index among endometrial cancers from The Cancer Genome Atlas Project.

OBJECTIVE: The Cancer Genome Atlas (TCGA) identified four integrated clusters for endometrial cancer (EC): POLE, MSI, CNL and CNH. We evaluated differences in gene expression profiles of obese and non-obese women with EC and examined the association of body mass index (BMI) within the clusters identified in TCGA. METHODS: TCGA RNAseq data was used to identify genes related to increasing BMI among ECs. The POLE, MSI and CNL clusters were composed mostly of endometrioid EC. Patient BMI was compared between these three clusters with one-way ANOVA. Association between gene expression and BMI was also assessed while adjusting for confounding effects of potential confounding factors. p-Values testing the association between gene expression and BMI were adjusted for multiple hypothesis testing over the 20,531 genes considered. RESULTS: Mean BMI was statistically different between the ECs in the CNL (35.8) versus POLE (29.8) cluster (p=0.006) and approached significance for the MSI (33.0) versus CNL (35.8) cluster (p=0.05). 181 genes were significantly up- or down-regulated with increasing BMI in endometrioid EC (q-value<0.01), including LPL, IRS-1, IGFBP4, IGFBP7 and the progesterone receptor. DAVID functional annotation analysis revealed significant enrichment in "cell cycle" (adjusted p-value=1.5E-5) and "DNA metabolic processes" (adjusted p-value=1E-3) for the identified genes. CONCLUSIONS: Obesity related genes were found to be upregulated with increasing BMI among endometrioid ECs. Patients with POLE tumors have the lowest median BMI when compared to MSI and CNL. Given the heterogeneity among endometrioid EC, consideration should be given to abandoning the Type I and II classification of EC tumors.

SigFuge: single gene clustering of RNA-seq reveals differential isoform usage among cancer samples.

High-throughput sequencing technologies, including RNA-seq, have made it possible to move beyond gene expression analysis to study transcriptional events including alternative splicing and gene fusions. Furthermore, recent studies in cancer have suggested the importance of identifying transcriptionally altered loci as biomarkers for improved prognosis and therapy. While many statistical methods have been proposed for identifying novel transcriptional events with RNA-seq, nearly all rely on contrasting known classes of samples, such as tumor and normal. Few tools exist for the unsupervised discovery of such events without class labels. In this paper, we present SigFuge for identifying genomic loci exhibiting differential transcription patterns across many RNA-seq samples. SigFuge combines clustering with hypothesis testing to identify genes exhibiting alternative splicing, or differences in isoform expression. We apply SigFuge to RNA-seq cohorts of 177 lung and 279 head and neck squamous cell carcinoma samples from the Cancer Genome Atlas, and identify several cases of differential isoform usage including CDKN2A, a tumor suppressor gene known to be inactivated in a majority of lung squamous cell tumors. By not restricting attention to known sample stratifications, SigFuge offers a novel approach to unsupervised screening of genetic loci across RNA-seq cohorts. SigFuge is available as an R package through Bioconductor.

Integrated RNA and DNA sequencing improves mutation detection in low purity tumors.

Identifying somatic mutations is critical for cancer genome characterization and for prioritizing patient treatment. DNA whole exome sequencing (DNA-WES) is currently the most popular technology; however, this yields low sensitivity in low purity tumors. RNA sequencing (RNA-seq) covers the expressed exome with depth proportional to expression. We hypothesized that integrating DNA-WES and RNA-seq would enable superior mutation detection versus DNA-WES alone. We developed a first-of-its-kind method, called UNCeqR, that detects somatic mutations by integrating patient-matched RNA-seq and DNA-WES. In simulation, the integrated DNA and RNA model outperformed the DNA-WES only model. Validation by patient-matched whole genome sequencing demonstrated superior performance of the integrated model over DNA-WES only models, including a published method and published mutation profiles. Genome-wide mutational analysis of breast and lung cancer cohorts (n = 871) revealed remarkable tumor genomics properties. Low purity tumors experienced the largest gains in mutation detection by integrating RNA-seq and DNA-WES. RNA provided greater mutation signal than DNA in expressed mutations. Compared to earlier studies on this cohort, UNCeqR increased mutation rates of driver and therapeutically targeted genes (e.g. PIK3CA, ERBB2 and FGFR2). In summary, integrating RNA-seq with DNA-WES increases mutation detection performance, especially for low purity tumors.

ABRA: improved coding indel detection via assembly-based realignment.

MOTIVATION: Variant detection from next-generation sequencing (NGS) data is an increasingly vital aspect of disease diagnosis, treatment and research. Commonly used NGS-variant analysis tools generally rely on accurately mapped short reads to identify somatic variants and germ-line genotypes. Existing NGS read mappers have difficulty accurately mapping short reads containing complex variation (i.e. more than a single base change), thus making identification of such variants difficult or impossible. Insertions and deletions (indels) in particular have been an area of great difficulty. Indels are frequent and can have substantial impact on function, which makes their detection all the more imperative. RESULTS: We present ABRA, an assembly-based realigner, which uses an efficient and flexible localized de novo assembly followed by global realignment to more accurately remap reads. This results in enhanced performance for indel detection as well as improved accuracy in variant allele frequency estimation. AVAILABILITY AND IMPLEMENTATION: ABRA is implemented in a combination of Java and C/C++ and is freely available for download at https://github.com/mozack/abra.

BlackOPs: increasing confidence in variant detection through mappability filtering.

Identifying variants using high-throughput sequencing data is currently a challenge because true biological variants can be indistinguishable from technical artifacts. One source of technical artifact results from incorrectly aligning experimentally observed sequences to their true genomic origin ('mismapping') and inferring differences in mismapped sequences to be true variants. We developed BlackOPs, an open-source tool that simulates experimental RNA-seq and DNA whole exome sequences derived from the reference genome, aligns these sequences by custom parameters, detects variants and outputs a blacklist of positions and alleles caused by mismapping. Blacklists contain thousands of artifact variants that are indistinguishable from true variants and, for a given sample, are expected to be almost completely false positives. We show that these blacklist positions are specific to the alignment algorithm and read length used, and BlackOPs allows users to generate a blacklist specific to their experimental setup. We queried the dbSNP and COSMIC variant databases and found numerous variants indistinguishable from mapping errors. We demonstrate how filtering against blacklist positions reduces the number of potential false variants using an RNA-seq glioblastoma cell line data set. In summary, accounting for mapping-caused variants tuned to experimental setups reduces false positives and, therefore, improves genome characterization by high-throughput sequencing.