SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Quantifying allele fraction of small to medium deletions in amplication sequencing brentonmar Bioinformatics 1 08-08-2013 06:56 AM
how to compute Minor Allele Frequency for a variant call? BhariD Bioinformatics 0 05-07-2013 11:50 AM
how does varscan somatic compute variant allele frequency? genomicist Bioinformatics 0 06-11-2012 04:21 AM
Variant calls with a low fraction of alt reads Jeremy37 Bioinformatics 9 04-17-2012 07:18 PM
A SNP discovery method to assess variant allele probability from NGS NGSfan Literature Watch 0 01-11-2010 04:48 AM

Reply
 
Thread Tools
Old 10-23-2013, 03:54 AM   #1
zhaopeihua
Member
 
Location: china

Join Date: Aug 2013
Posts: 18
Default variant allele fraction

hi:
What are "variant allele fraction" when talking about somatic mutations?
zhaopeihua is offline   Reply With Quote
Old 10-23-2013, 03:53 PM   #2
lethalfang
Member
 
Location: San Francisco, CA

Join Date: Aug 2011
Posts: 91
Default

Quote:
Originally Posted by zhaopeihua View Post
hi:
What are "variant allele fraction" when talking about somatic mutations?
The fraction of alleles in the sample that are mutations.
For instance, if all cells contains the same heterogeneous mutation (e.g., 100% tumor cell purity if you're looking at tumor cells), the variant allele fraction would be 50% (half-half from each chromosome).
lethalfang is offline   Reply With Quote
Old 10-23-2013, 11:13 PM   #3
zhaopeihua
Member
 
Location: china

Join Date: Aug 2013
Posts: 18
Default

Quote:
Originally Posted by lethalfang View Post
The fraction of alleles in the sample that are mutations.
For instance, if all cells contains the same heterogeneous mutation (e.g., 100% tumor cell purity if you're looking at tumor cells), the variant allele fraction would be 50% (half-half from each chromosome).
thanks for you reply,but I didn't fully understand this concept. Could you show me a detail example with computational process? thanks again
zhaopeihua is offline   Reply With Quote
Old 10-23-2013, 11:27 PM   #4
lethalfang
Member
 
Location: San Francisco, CA

Join Date: Aug 2011
Posts: 91
Default

This is an example of an output in my somatic calling pipeline:

Code:
Variant Annotation	Gene	Exon	Codon	Chr	Position	Variant	N(A)	N(C)	N(G)	N(T)	T(A)	T(C)	T(G)	T(T)	QSS	Frequency
nonsynonymous SNV	KRAS	exon2	G12A	chr12	25398284	G35C	0	83	0	0	1	37	12	0	57	0.24

N(A), N(C), N(G), and N(T) stand for the number of calls for A, C, G, and T in the normal sample. T(A), T(C), T(G), and T(T) stand for the number of calls for A, C, G, and T in the tumor sample.

So for the normal sample, there are 83 sequences where the base is called C, and 0 for everything else. So it's pretty clear that in the reference and the normal sample, the base is C.
However, in the tumor sample, there are 37 calls for C (reference) and 12 calls for G (variant), so the variant frequency is 12/(12+37) = 0.245.
Notice the variant is labeled G35C. That's because the coding strand is the minus strand, where the DNA reads use the plus strand.

Anyway, how do you interpret 24.5%?
Well, if I assume that
1) the KRAS mutation is present in all tumor cells (just use this assumption as an example), and
2) KRAS mutation is a heterozygous mutation, then
The tumor sample must contain 1/2 non-tumor cells, because those cells give me half of the reference reads. Out of the remaining 1/2 are tumor cells, half of those chromosomes give me reference reads, and the other 1/2 chromosomes give me mutant reads.
lethalfang is offline   Reply With Quote
Old 10-24-2013, 01:12 AM   #5
zhaopeihua
Member
 
Location: china

Join Date: Aug 2013
Posts: 18
Default

Awesome, thanks!
zhaopeihua is offline   Reply With Quote
Old 09-03-2014, 05:18 AM   #6
vd4mindia
Member
 
Location: Milan

Join Date: May 2013
Posts: 40
Default

@lethalfang
`Would you like to share the pipeline for finding somatic mutations. I have 2 samples, one low and one high grade tumor and its normal. I want to detect the point somatic mutations. I have already employed GATK, VarScan and Mutect, but I have some descrepancy when am looking at my Low grade tumor since it is 50% pure. So I would like to try something more and have a read statistics like you have shown to validate the hits. If its ok for you
vd4mindia is offline   Reply With Quote
Old 09-10-2014, 10:12 AM   #7
lethalfang
Member
 
Location: San Francisco, CA

Join Date: Aug 2011
Posts: 91
Default

Quote:
Originally Posted by vd4mindia View Post
@lethalfang
`Would you like to share the pipeline for finding somatic mutations. I have 2 samples, one low and one high grade tumor and its normal. I want to detect the point somatic mutations. I have already employed GATK, VarScan and Mutect, but I have some descrepancy when am looking at my Low grade tumor since it is 50% pure. So I would like to try something more and have a read statistics like you have shown to validate the hits. If its ok for you
That's something I used to do a couple of years ago.
I just had a script to get those information out of the Strelka output vcf file. I had all those Strelka out directories in a single location, and ran this script to "summarize" all mutation calls from all the samples:

http://kimlab.surgery.ucsf.edu/media...a_findings.txt

It also calls annovar, which has a hard-coded path within the script.
lethalfang is offline   Reply With Quote
Old 03-04-2015, 05:30 AM   #8
virpericulosus
Junior Member
 
Location: Nebraska

Join Date: Sep 2014
Posts: 3
Default

Just wondering: The T(A) isn't counted because the frequency is just for T(G) variants?
virpericulosus is offline   Reply With Quote
Old 03-25-2015, 11:43 AM   #9
vd4mindia
Member
 
Location: Milan

Join Date: May 2013
Posts: 40
Default

Quote:
Originally Posted by lethalfang View Post
This is an example of an output in my somatic calling pipeline:

Code:
Variant Annotation	Gene	Exon	Codon	Chr	Position	Variant	N(A)	N(C)	N(G)	N(T)	T(A)	T(C)	T(G)	T(T)	QSS	Frequency
nonsynonymous SNV	KRAS	exon2	G12A	chr12	25398284	G35C	0	83	0	0	1	37	12	0	57	0.24

N(A), N(C), N(G), and N(T) stand for the number of calls for A, C, G, and T in the normal sample. T(A), T(C), T(G), and T(T) stand for the number of calls for A, C, G, and T in the tumor sample.

So for the normal sample, there are 83 sequences where the base is called C, and 0 for everything else. So it's pretty clear that in the reference and the normal sample, the base is C.
However, in the tumor sample, there are 37 calls for C (reference) and 12 calls for G (variant), so the variant frequency is 12/(12+37) = 0.245.
Notice the variant is labeled G35C. That's because the coding strand is the minus strand, where the DNA reads use the plus strand.
I would just want to ask you in case you have have not other bases as 0 in normal in that case do you still calculate the variant allele frequency as the above? Variant allele frequency should be always considered for the reads having variant allele for the tumor sample divided by the total reads at that variant site? Right? Irrespective of the reads falling for other bases in normal samples as well right?
vd4mindia is offline   Reply With Quote
Old 03-25-2015, 11:46 AM   #10
lethalfang
Member
 
Location: San Francisco, CA

Join Date: Aug 2011
Posts: 91
Default

Quote:
Originally Posted by vd4mindia View Post
I would just want to ask you in case you have have not other bases as 0 in normal in that case do you still calculate the variant allele frequency as the above? Variant allele frequency should be always considered for the reads having variant allele for the tumor sample divided by the total reads at that variant site? Right? Irrespective of the reads falling for other bases in normal samples as well right?
I considered "variant allele frequency" as the frequency of variant alleles in a given sample.
So the tumor's variant allele frequency does not care what's in the normal tissue, and vice versa.
That way, you can monitor the change in variant allele frequency between the tumor and the normal.
lethalfang is offline   Reply With Quote
Old 03-26-2015, 08:01 AM   #11
vd4mindia
Member
 
Location: Milan

Join Date: May 2013
Posts: 40
Default

Thank you very much for the reply. That makes sense.
vd4mindia is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 01:42 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO