Question: Find the 3' utrs, 5'utrs and their counts from bam files
0
gravatar for ilovesuperheroes1993
8 weeks ago by
ilovesuperheroes19930 wrote:

Hello,

I have bam files of 8 samples (4 normal and 4 diseased), produced by alignment with novoalign (small-rna sequencing data). I have excluded the mirna from the bam files by using the following command:

bedtools intersect -v "sample.bam" "hg19_mirna.gff3" > output.bam

In this manner I have excluded the miRNA from all the samples. Now I want to find the utr (3' and 5') sequences present in the resultant bam files, and the read counts of each of the utr sequences.

Could anyone suggest a way to do this?

ucsc utr bam bed bedtools • 129 views
ADD COMMENTlink modified 8 weeks ago by Gautier Richard270 • written 8 weeks ago by ilovesuperheroes19930

You can try featureCounts with meta-feature "UTR" level.

ADD REPLYlink modified 8 weeks ago • written 8 weeks ago by b.nota6.3k

Thank you. Actually I want to get a gtf file of the utr sequences, and then count the reads using htseq-count. But I am unable to find the gtf files of 3' and 5' utrs. Do you know where I could get it?

ADD REPLYlink written 8 weeks ago by ilovesuperheroes19930

at ensembl you have gtf files for many organisms. They contain "UTR" metadata (I know for human and mouse).

ADD REPLYlink modified 8 weeks ago • written 8 weeks ago by b.nota6.3k

Hi, I am currently using hg19 database (Grch37 version). I am unable to find the gtf / gff3 files of the utrs of this version. Could you please link me to them? Thanks a lot

ADD REPLYlink written 8 weeks ago by ilovesuperheroes19930

I recommend you upgrade to hg38! Or search in the archives somewhere.

ADD REPLYlink modified 8 weeks ago • written 8 weeks ago by b.nota6.3k
1
gravatar for Gautier Richard
8 weeks ago by
MPI IE, Freiburg, Germany
Gautier Richard270 wrote:

Using biomart at ensembl you can get the 3' and 5' UTR of hg19 (image below) at this address:

https://grch37.ensembl.org/biomart/martview/

Then I guess you can get the read coverage from these sequences using for example deepTools multiBamSummary BED-file

https://deeptools.readthedocs.io/en/develop/content/tools/multiBamSummary.html

grch37.ensembl.org UTRs

ADD COMMENTlink modified 8 weeks ago • written 8 weeks ago by Gautier Richard270
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: 1888 users visited in the last hour