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!
1 answer
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.
Log in to answer this question.
-
Annotation comparison files on NCBI of gene assembly
written by zirui •
-
WGCNA from Salmon Output
written by Adrij Pawan •
-
How to normalize long-read RNA-seq data for comparison with short-reads
written by Bernardo •
-
Detection of Novel splice events
written by Nelislam4 •
-
How do I download mm10 major release genomic fasta and gtf file?
written by sg197 •
-
Can I compare TPM values for a transcript across different transcriptome assemblies?
written by molly77 •
-
Measure expression level of a single (novel) transcript in multiple fastq files
written by c_u
-
Inaccurate read counts when genome/annotation contains alternate loci
written by viktorfeketa •
-
Looking for a thorough annotation for non-primary assembly units in GRCm38
written by erwan.scaon
-
How to calculate TPM of transcripts with isoforms
written by marcosp.belem •
