Question: From transcript to gene - Salmon quantification
0
gravatar for JulianC
3 months ago by
JulianC10
JulianC10 wrote:

Hi!

I am working with Salmon to quantify gene expression from single-cell RNA-seq data. I used NCBI transcriptome (in FASTA format; https://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/001/635/GCF_000001635.26_GRCm38.p6/GCF_000001635.26_GRCm38.p6_rna.fna.gz) as index, and my quantification went pretty well. The output gave me TPM values for each transcript. For example, for a gene with three transcripts, I end up with three different quantifications, one for each transcript. Since I want to have a gene-specific quantification, is it fair to sum up the TPM values of the three transcripts to obtain one single quantification relative to that gene? Thank you in advance!

salmon • 150 views
ADD COMMENTlink modified 3 months ago by ATpoint36k • written 3 months ago by JulianC10
6
gravatar for ATpoint
3 months ago by
ATpoint36k
Germany
ATpoint36k wrote:

Basically yes, but there is dedicated software for this. I suggest you use tximport to summarize the transcript level estimates to the gene level. This on the one hand provides you with gene level counts and on the other hand calculates the average and sample-specific length for each gene depending on which transcripts/isoforms contribute to the gene level count. The tximport vignette and paper discusses details.

ADD COMMENTlink modified 3 months ago • written 3 months ago by ATpoint36k
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: 918 users visited in the last hour