Question: Suggested pipelines for finding somatic mutations using RNAseq in normal and tumor cells
gravatar for Sharon
21 months ago by
Sharon440 wrote:

Hi All

I am looking for pipelines to find somatic mutations in tumor and normal cells. Should I use GATK pipeline for RNAseq variant calling here (; apply it on each normal cell then apply it each on tumor cell then find the mutations in tumor cells that are not found in normal cells?


rna-seq somatic mutations • 1.0k views
ADD COMMENTlink modified 9 months ago by igor8.3k • written 21 months ago by Sharon440

I don't believe that RNA-seq is an appropriate technology for your application.

ADD REPLYlink written 21 months ago by WouterDeCoster40k

But there is a variant calling recommended for RNAseq I think ?

ADD REPLYlink written 21 months ago by Sharon440

Yes, but just because you can doesn't mean you should. That pipeline is also intended for germline variants, not for somatic AFAIK.

ADD REPLYlink written 21 months ago by WouterDeCoster40k

I mean there is interest in finding somatic mutations from RNAseq, there are publications on that:

ADD REPLYlink written 21 months ago by Sharon440

I'm not sure that this would be my first choice as an experimental set-up, i.e., to call variants from RNA-seq data. The idea appears to be about saving money for cash-strapped research groups. However, sacrificing quality in the sake of money-saving invariably runs to disaster later down the line, as we see time and time again in various facets of life.

ADD REPLYlink written 21 months ago by Kevin Blighe48k

Have you solve the question? I also meet the same problem recently.

ADD REPLYlink written 9 months ago by 7750566480
gravatar for kristoffer.vittingseerup
9 months ago by
European Union
kristoffer.vittingseerup2.3k wrote:

I would be VERY careful with such analysis. According to this paper which systematically compared DNA and RNA sequencing:

  • 65% of DNA based mutations were not found in RNA (Low sensitivity)
  • 92% of RNA based mutations were not detected using DNA (very high false discovery rate)
  • RNASeq typically only have power to detect 33% of DNA mutations
ADD COMMENTlink written 9 months ago by kristoffer.vittingseerup2.3k
gravatar for igor
9 months ago by
United States
igor8.3k wrote:

Variants from RNA-seq is a controversial topic. Although there are arguments against (as other have already mentioned), it's also not entirely unreasonable. There is an interesting statement from Coudray et al (where they also make some pipeline recommendations):

The overlap between RNA-seq and WES was small in all samples, but interestingly, the overlap increased with increasing significance of the variants. An average of only 6.60% of WES variants retained by MuTect2 (PASS) were also present in RNA-seq, while 15.9% of WES variants from coding regions and 17.2% of functional mutations were common to RNA-seq (Fig. 2D). Coverage differences between RNA-seq and WES could partially explain the phenomenon. A previous study indeed found that ∼71% of RNA-seq variants fell outside the WES capture boundaries (O’Brien et al., 2015). Moreover, they showed that a high proportion of RNA-seq-only variants were missed by WES because of their low allele fraction (AF).

Their conclusion:

Our work suggests that since the majority of studies on cancer-driving mutations used WES-only, they are likely to have missed some key driver mutations that might be found using complementary RNA-seq datasets from the same tumors.

ADD COMMENTlink written 9 months ago by igor8.3k
Please log in to add an answer.


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