How can I look at nonsense mutations in a particular gene in TCGA data?
2
1
Entering edit mode
9.0 years ago
biohack92 ▴ 170

I have been working with TCGA cancer data to examine expression (RNAseqV2) and methylation (Illumina 450k) data. I want to look at sequencing data, but I'm a bit lost to what sort of information is available through TCGA. I want to examine whether there are nonsense mutations between positions 2000-2500 across all cancer types available in TCGA. What sort of resources/workflow should I expect?

sequence mutations tcga • 3.1k views
ADD COMMENT
2
Entering edit mode
9.0 years ago
Nicholas Spies ★ 1.2k

If you are asking whether you can search for non-sense mutations between amino acids 2000-2500 of a particular gene in the TCGA sets, I would suggest using the Mutation Annotation Files (MAFs). These files have already been run through somatic variation calling so you won't have to deal with the sequencing data directly. You can filter out your gene of interest across all the MAF's in perl,awk,grep etc... Then in the amino_acid_change column you can search for integer values between p.2000-p.2500, and trv_type "Non-sense mutation".

Hope this helps.

Nick

ADD COMMENT
0
Entering edit mode

Appreciate it! Thank you

ADD REPLY
2
Entering edit mode
9.0 years ago

For your specific question, your best bet is to query cBioPortal using its Onco Query Language (OQL). For example, the following query gets you all nonsense mutations seen in PTEN across all cancer types:

PTEN: MUT = NONSENSE

We also use the more general term TRUNC that pools together ORF truncating events NONSENSE, NONSTART, FRAMESHIFT, and SPLICE.

If you want to do anything more complicated, then I agree with Nick that you need to get the MAF files. Here's an intro to MAF files: Working with MAF files (Mutation Annotation Format) from the TCGA (The Cancer Genome Atlas)

ADD COMMENT
0
Entering edit mode

Yes, I absolutely agree. cBioPortal is a very powerful resource as well. Thanks for the suggestion and sorry for my oversight!

ADD REPLY

Login before adding your answer.

Traffic: 2565 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6