To estimate somatic copy number alternations, ASCAT/Battenberg was performed on multi-regional paired tumour-normal sequencing data. Allele counts of positions from 1000 genomes were generated using AlleleCounter, and minimum coverage of 20 for normal sample was used for filtration. LogR and BAF values were produced for each region, and concatenated into one matrix separately for each patient. LogR values were subsequently corrected using a GC wave correction implemented in ASCAT, and only heterozygous BAF values were reserved for further analysis. Allele-specific multi-sample segmentation was performed to generate segmented logR and BAF data by ascat.asmultipcf. Both ASCAT and Battenberg algorithms were applied to provide cellularity (purity) and ploidy estimates with gamma setting of 1, and manual verification was used to select the optimal model for ploidy and cellularity using an orthogonal measures based on ABSOLUTE results and mutation variant allele fraction. And then ASCAT/Battenberg was re-run to obtain the final allele-specific copy number data using reviewed cellularity and ploidy. To assess SCNA calls, we compared absolute copy number profiles using four software packages – ABSOLUTE, CNVKit, ASCAT and Battenberg. All the samples showed high concordance except patient 12 R1 and R4, which had no best solution in ABSOLUTE.
Do you have any questions about this protocol?
Post your question to gather feedback from the community. We will also invite the authors of this article to respond.