Question: BAF without normal control, to do or not! or how
0
gravatar for sm.hashemin
6 weeks ago by
sm.hashemin90
sm.hashemin90 wrote:

Dear Eric, dear all, from the cnvkit documentation :

Typically you would use a properly formatted VCF from joint tumor-normal SNV calling, e.g. the output of MuTect, VarDict, or FreeBayes, having already flagged somatic mutations so they can be skipped in this analysis. If you have no matched normal sample for a given tumor, you can use 1000 Genomes common SNP sites to extract the likely germline SNVs from a tumor-only VCF, and use just those sites with THetA2 (or another tool like PyClone or BubbleTree).

I am currently trying to do the same as above mentioned. I have made a .cnn reference from unrelated but aged matched WES files with obatiained from the same hybridziation based method. I have filtered my VCFs for common dbsnp SNPs with AF of more that 10% (very very common).

my cns file looks like this

chromosome  start   end gene    log2    depth   probes  weight
chr1    12403   2990008 DDX11L1,WASH7P,FAM138F,MIR4251  -2.49508    14.0935 1681    555.806
chr1    2992142 6337653 PRDM16,MIR4251,ACOT7    -4.62846    7.05581 1631    483.703

my vcf file without header looks like this

#CHROM  POS ID  REF ALT QUAL    FILTER  INFO    FORMAT  Tumorsample1
chr1    762273  .   G   A   .   REJECT  DP=367;AF=1;DP4=0,0,119,248;SB=0;ANNOVAR_DATE=2018-04-16;ExAC_ALL=0.8060;ExAC_AFR=0.4415;ExAC_AMR=0.8116;ExAC_EAS=0.9174;ExAC_FIN=0.9;ExAC_NFE=0.8384;ExAC_OTH=0.8896;ExAC_SAS=0.8184;Func.refGene=ncRNA_exonic;Gene.refGene=LINC00115;GeneDetail.refGene=.;ExonicFunc.refGene=.;AAChange.refGene=.;cosmic87_coding=.;ALLELE_END;rs_ids=rs3115849   GT:DP:AF:SB:DP4 .:367:1.0:0:0,0,119,248

when I use export theta tumorsample1.cns -r ref.cnn -v sample1.vcf

  Wrote sample1theta Selected test sample sample1 Loaded 44443
records; skipped: 0 somatic, 1648 depth Kept 44443 heterozygous of
44443 VCF records
/anaconda3/lib/python3.7/site-packages/pandas/core/indexing.py:868:
FutureWarning:  Passing list-likes to .loc or [] with any missing
label will raise KeyError in the future, you can use .reindex() as an
alternative.

See the documentation here:
https://pandas.pydata.org/pandas-docs/stable/indexing.html#deprecate-loc-reindex-listlike
return self._getitem_lowerdim(tup) 
Wrote sample1.tumor.snp_formatted.txt 
Wrote sample1.normal.snp_formatted.txt

Unfortunately I get an empty normal SNP file(should it be like this or it is an error?)

#Chrm   Pos Ref_Allele  Mut_Allele

the tumor snp file has zero as Mut_Allele count

#Chrm   Pos Ref_Allele  Mut_Allele
chr1    762272  367 0
chr1    808921  153 0

interval file

#ID chrm    start   end tumorCount  normalCount start_1_12403:end_1_2990008 1   12403   2990008 298177  3286017
start_1_2992142:end_1_6337653   1   2992142 6337653 65940   1608940

Theta stops prematurely as the normalMutCount[i] + normalRefCount[i] is less than min

the common SNPs were extracted but how would you get the BAF of the normal sample?

Best regards and thanks in advance

Eric T.

ADD COMMENTlink modified 6 weeks ago • written 6 weeks ago by sm.hashemin90
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 562 users visited in the last hour