Hi,
I have performed targeted amplicon sequencing for paired normal-tumor samples and recently identified both somatic mutations (VarScan2 somatic) and copynumber changes (CNVkit's batch command) in my samples. I would like to move on and study the tumour heterogeneity within my samples with PyClone.
Manual of CNVkit describes how data is exported into THetA2. However, as a beginner I find the manual little bit hard-to-read especially what comes to heterogeneity analysis and I can't figure out how I should run the program to obtain the input file for PyClone. Based on previous posts I know that some people have had good results with CNVkit+PyClone. I would like to see similar example script for the use with PyClone as well. Could anyone provide an example?
Thank you in advance!
1 answer
I don't have a script for doing this, but if someone else has one I'd be happy to include it in the main CNVkit distribution. (Same goes for documentation.)
A good start might be to run export nexus-ogt to align b-allele frequencies to segments, then process or edit that output file further to format it for PyClone.
Log in to answer this question.
-
Consistent somatic CNV calling across samples: spurious, error or real?
written by avelarbio46 •
-
Tumour only B-allele frequency
written by samuel
-
Mutation calling for pediatric cancer data
written by Gene_MMP8
-
germline calling of tumor-normal paired sample germline by varscan2
written by oghzzang •
-
Genotype information required in Clonality analysis using Pyclone
written by nilesh.gardi2688 •
-
Tumor purity estimation by allele frequency of COSMIC identified somatic mutations
written by ejoffe •
-
Creation of tsv and xml file for Pyclone
written by geek7 •
-
Using cnvkit for targeted amplicon sequencing
written by ingrid.schulman •
-
Using THETA2 from CNVkit: what king of vcf do you need?
written by gil.hornung
-
Why are there too many mutations in controls, but not in tumour?
written by MAPK
Hi Eric, i have the same question with you. I have download the WES data and used mutect2 to call somatic mutations and CNVkit to call CNV. But i really don't know how can i create the tsv file for pyclone. If you have solved this problem, please teach me how to do. Thank you.