Published in

Oxford University Press, Bioinformatics, 22(37), p. 4023-4032, 2021

DOI: 10.1093/bioinformatics/btab440

Links

Tools

Export citation

Search in Google Scholar

Detecting copy number alterations in RNA-Seq using SuperFreq

Journal article published in 2021 by Christoffer Flensburg, Alicia Oshlack, Ian J. Majewski ORCID
This paper was not found in any repository, but could be made available legally by the author.
This paper was not found in any repository, but could be made available legally by the author.

Full text: Unavailable

Green circle
Preprint: archiving allowed
Orange circle
Postprint: archiving restricted
Red circle
Published version: archiving forbidden
Data provided by SHERPA/RoMEO

Abstract

Abstract Motivation Calling copy number alterations (CNAs) from RNA sequencing (RNA-Seq) is challenging, because of the marked variability in coverage across genes and paucity of single nucleotide polymorphisms (SNPs). We have adapted SuperFreq to call absolute and allele sensitive CNAs from RNA-Seq. SuperFreq uses an error-propagation framework to combine and maximize information from read counts and B-allele frequencies. Results We used datasets from The Cancer Genome Atlas (TCGA) to assess the validity of CNA calls from RNA-Seq. When ploidy estimates were consistent, we found agreement with DNA SNP-arrays for over 98% of the genome for acute myeloid leukaemia (TCGA-AML, n = 116) and 87% for colorectal cancer (TCGA-CRC, n = 377). The sensitivity of CNA calling from RNA-Seq was dependent on gene density. Using RNA-Seq, SuperFreq detected 78% of CNA calls covering 100 or more genes with a precision of 94%. Recall dropped for focal events, but this also depended on signal intensity. For example, in the CRC cohort SuperFreq identified all cases (7/7) with high-level amplification of ERBB2, where the copy number was typically >20, but identified only 6% of cases (1/17) with moderate amplification of IGF2, which occurs over a smaller interval. SuperFreq offers an integrated platform for identification of CNAs and point mutations. As evidence of how SuperFreq can be applied, we used it to reproduce the established relationship between somatic mutation load and CNA profile in CRC using RNA-Seq alone. Availability and implementation SuperFreq is implemented in R and the code is available through GitHub: https://github.com/ChristofferFlensburg/SuperFreq/. Data and code to reproduce the figures are available at: https://gitlab.wehi.edu.au/flensburg.c/SuperFreq_RNA_paper. Data from TCGA (phs000178) was accessed from GDC following completion of a data access request through the database of Genotypes and Phenotypes (dbGaP). Data from the Leucegene consortium was downloaded from GEO (AML samples: GSE67040; normal CD34+ cells: GSE48846). Supplementary information Supplementary data are available at Bioinformatics online.