Correcting for cancer genome size and tumour cell content enables better estimation of copy number alterations from next-generation sequence data.
BACKGROUND
Comparison of read depths from next-generation sequencing between cancer and normal cells makes the estimation of copy number alteration (CNA) possible, even at very low coverage. However, estimating CNA from patients' tumour samples poses considerable challenges due to infiltration with normal cells and aneuploid cancer genomes. Here we provide a method that corrects contamination with normal cells and adjusts for genomes of different sizes so that the actual copy number of each region can be estimated.


RESULTS
The procedure consists of several steps. First, we identify the multi-modality of the distribution of smoothed ratios. Then we use the estimates of the mean (modes) to identify underlying ploidy and the contamination level, and finally we perform the correction. The results indicate that the method works properly to estimate genomic regions with gains and losses in a range of simulated data as well as in two datasets from lung cancer patients. It also proves a powerful tool when analysing publicly available data from two cell lines (HCC1143 and COLO829).


BACKGROUND
An R package, called CNAnorm, is available at http://www.precancer.leeds.ac.uk/cnanorm or from Bioconductor.


BACKGROUND
a.gusnanto@leeds.ac.uk


BACKGROUND
Supplementary data are available at Bioinformatics online.
Department of Statistics, University of Leeds, Leeds LS2 9JT, UK.
Arief
A
Gusnanto
Henry M
HM
Wood
Yudi
Y
Pawitan
Pamela
P
Rabbitts
Stefano
S
Berri
D045744:Cell Line, Tumor; D003198:Computer Simulation; D056915:DNA Copy Number Variations; D059646:Genome Size; D059014:High-Throughput Nucleotide Sequencing; D006801:Humans; D008175:Lung Neoplasms / Q000235:genetics; D009369:Neoplasms / Q000235:genetics*; D017422:Sequence Analysis, DNA; D012984:Software