This is a test version of Biostars. For the public version, visit https://www.biostars.org.
how can i quantify expression using own my gff ?

Hi,

I want to quantify the leads in my defined region.

This is not a gene, but an enhancer region. (from ENCODE)

The data is a bam file that maps total RNA-seq using tophat.

I want to measure the expression of a region that is thought to be an enhancer RNA that is not defined in human gtf.

So I wondered if I could use the gffs for the regions I want to use for quantification.

But, I am not sure which tool to use and how.

Thanks you.

rna-seq rna-seq
• 1,347 views • link
updated by wuyc • wuyc • written by dlyedms1019 • dlyedms1019

2 answers

Both featureCounts and htseq count will take any GTF file to get regions to quantify over. I recommend featureCounts - its faster and has more features. You just need to make sure you tell featureCounts which feautre type to quantify over (depends what feature type your enhancers are enocded as) and tell it to aggregate at the feature level, not the gene or transcript level.

• 0 views • link
written by i.sudbery i.sudbery

Thanks you for your comments! I made a SAF(simple annotation format) for featureCounts and quantified eRNA and mRNA all.

• 0 views • link
written by dlyedms1019 • dlyedms1019

Just check this pipeline for enhancer RNA quantification from RNA-seq,

http://fun-science.club/PET/

You can also check the paper: Wu, Y, Yang, Y, Gu, H, et al. Multi‐omics analysis reveals the functional transcription and potential translation of enhancers. Int. J. Cancer. 2020; 1– 15. https://doi.org/10.1002/ijc.33132

• 0 views • link
written by wuyc • wuyc

Log in to answer this question.

More posts like this