SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Raw read count to reference genes fanx Bioinformatics 3 04-11-2014 04:34 AM
reads mapping to overlapping genes? beliefbio Bioinformatics 3 04-11-2014 12:18 AM
multiBamCov or htseq-count to count read per feature ? NicoBxl Bioinformatics 1 07-03-2012 02:05 AM
RNASeq removing PCR bias schaffer Bioinformatics 3 10-20-2011 06:00 PM
Overlapping and non-Overlapping pair-end reads with Tophat senpeng Illumina/Solexa 4 10-16-2011 06:43 PM

Reply
 
Thread Tools
Old 11-06-2012, 07:28 AM   #1
DRAT
Junior Member
 
Location: Belgium

Join Date: Feb 2012
Posts: 4
Default Removing overlapping genes from annotation for RNAseq read count

Hello,

I am trying to prepare a read count table for DESeq using EasyRNAseq package in R. I followed the vignette and used ensembl.gtf file as my annotation. After constructing my read count table I get warnings about overlapping genes and counting reads more than once but I am not sure how to modify my annotation to avoid it. The manual only mentions that computed gene models can be extracted from created RNAseq object and that overlapping loci should be removed without specifying how to do it. I am able to extract gene models but I am not sure how to correctly process it before re-running the function. Could anyone please give me some advice on how to fix this annotation in R?

I attached the code I used to generate my read count table below:

read.count <-easyRNASeq(format='bam',readLength=50L, organism="Mmusculus", chr.sizes="auto", annotationMethod="gtf", annotationFile="mm9.ensgene.gtf", count="genes", summarization="geneModels", filesDirectory=getwd(), filenames=c("NI_A_accepted_hits.bam","NI_B_accepted_hits.bam","DEX_A_accepted_hits.bam", "DEX_B_accepted_hits.bam","GW_A_accepted_hits.bam", "GW_B_accepted_hits.bam", "DEX_GW_A_accepted_hits.bam", "DEX_GW_B_accepted_hits.bam"), conditions=conditions, outputFormat="RNASeq")

To get gene models I used: geneModels <- geneModel(read.count) but I am stack at this point and I cannot find a way to remove overlapping features. I tried disjoin function but it gives an error: "Error in function (classes, fdef, mtable) : unable to find an inherited method for function "disjoin", for signature "RangedData" "

Thanks a lot for your suggestions!
DRAT is offline   Reply With Quote
Old 04-11-2014, 12:17 AM   #2
eszter.ari
Junior Member
 
Location: Vienna

Join Date: Dec 2012
Posts: 8
Default

Take a look at our read counter tool that does not count read-pairs mapped to the same-strand overlapping part of genes: http://seqanswers.com/forums/showthread.php?p=134850
https://code.google.com/p/recog/
eszter.ari is offline   Reply With Quote
Old 04-11-2014, 03:53 AM   #3
TiborNagy
Senior Member
 
Location: Budapest

Join Date: Mar 2010
Posts: 329
Default

You can use bedtools merge to collapse the overlapping regions in mm9.ensgene.gtf
TiborNagy is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 12:41 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO