Question: Gene level from transcript level expression estimate.
gravatar for mforde84
3.2 years ago by
mforde841.3k wrote:

Would it be appropriate to simply sum FPKM/RPKM/TPM normalized data across transcripts to generate a gene level estimate? Each transcript will be abundant at a certain % of total expression, so the intergration should simply be the gene level expression of all associated transcripts. Am I correct or am I missing something?

rna-seq • 1.6k views
ADD COMMENTlink modified 3.2 years ago by eu.sanisa0 • written 3.2 years ago by mforde841.3k

I'd think the answer is no, depending on how 'gene level estimate' defined. Since transcripts associated with a gene can share exons, I'd think summing counts including overlapping transcripts from the same gene would lead to multiple counting. Or did you mean summing reads over all exons associated with a gene locus?

ADD REPLYlink written 3.2 years ago by Ahill1.9k
gravatar for eu.sanisa
3.2 years ago by
eu.sanisa0 wrote:

If by transcripts you mean isoforms, then i think it is ok. I was wondering the same so i took a look at Trinity's . There is an option to give a "transcripts to gene map" and obtain a count matrix at the gene level. As I understand from the script, this uses either or . This last one has the statement: "compute it based on the formula: FPKM_gene = FPKM_isoA + FPKM_IsoB + ... "

ADD COMMENTlink written 3.2 years ago by eu.sanisa0
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: 992 users visited in the last hour