This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Quantification of repeats expression

Hi guys,

I want to calculate repeats expression in my RNAseq data. I've obtained bam files using TopHat and now I need gtf file for repeats to calculate the counts. Where can I download it?

Thanks!

ngs rna-seq repeats
• 3,711 views • link
updated by Constantine Constantine • written by valerie valerie

Remember to get a GTF file that matches your genome. If your genome came from Ensembl then you need to get the GTF from ensembl. Chromosome identifiers may otherwise not match.

BTW: Repeat tracks are under "Variation and Repeats" group in UCSC table browser.

• 0 views • link
written by GenoMax GenoMax

1 answer

Go to UCSC

https://genome.ucsc.edu/

and under Tools > Table Browser

Choose your genome and track (ideally RefSeq genes), and select "Output format: GTF"

• 0 views • link
written by Constantine Constantine

Thank you! I need, repeats, why RefSeq genes? Is it correct to choose 'Variation and Repeats' as group and 'RepeatMasker' as track?

• 0 views • link
written by valerie valerie

Yes. See my comment above.

• 0 views • link
written by GenoMax GenoMax

Thank you! I used mouse mm10 genome for TopHat and will use mm10 here again.

• 0 views • link
written by valerie valerie

Did the genome come from UCSC or Ensembl or someplace else? Also keep in mind the "multi-hits" setting for TopHat. Since you are interested in repeats that setting may affect your results significantly.

• 0 views • link
written by GenoMax GenoMax

That is correct. So you are fine with getting the repeats GTF from UCSC.

• 0 views • link
written by GenoMax GenoMax

Thank you for your help!

• 0 views • link
written by valerie valerie

Hi Valerie,

I have a similar project, working on SSR repeats. Could you please kindly tell me what is your workflow for doing the work?

• 0 views • link
written by seta seta

Hi Seta,

I simply use tophat2 to map the reads to reference genome. Then I sort my reads using samtools and use htseq-count to obtain counts from bam file. On this stage I needed gtf file we discussed here. Then you can the apply any normalization to counts, I prefer DeSeq. Let me know if still you have questions.

• 0 views • link
written by valerie valerie

Hi friend, Thank you very much for your explanation. As you mentioned "repeat" in the title of your question, I thought that you have a specific way for surveying these regions. Now, I found that you follow the common way.

• 0 views • link
written by seta seta

Log in to answer this question.

More posts like this