SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Star with cufflinks Nino Bioinformatics 20 12-12-2016 09:06 AM
Cufflinks refuses to operate on Tophat2 created bam or sam files due to sorting error amrezans Bioinformatics 1 06-24-2013 12:54 PM
Cufflinks 0.9.2 BAM bug? kopi-o Bioinformatics 38 05-21-2013 01:47 PM
BAM from Bioscope in Cufflinks 0.9.x AdamB Bioinformatics 7 10-07-2010 05:27 AM
tophat/cufflinks bam vs. RPKM mgogol Bioinformatics 5 04-26-2010 09:58 AM

Reply
 
Thread Tools
Old 03-11-2015, 01:34 AM   #1
gogodidi
Junior Member
 
Location: Netherlands

Join Date: Aug 2012
Posts: 4
Default STAR + Cufflinks. Cufflinks hanging. BAM XS error?

Hello All,

I have an RNA-Seq pipeline, using the STAR aligner and Cufflinks.
I am running cufflinks in the 'Reference Annotation Based Transcript' assemvbly mode to predict alternative splicing variants.
I have been using Cufflinks on our compute servers for a while now, and I didn't have any problems until recently.
For the past few weeks, Cufflinks always hangs after a few hours of running.
I can't seem to identify the problem. It doesn't seem to be a network issue, as I have installed cufflinks locally and also have the data locally.

I have tried this with cufflinks 2.2.1 and 2.1.1, and both hang. Using STAR version 2.4.0f1.

The only error I get is associated with the STAR output (and I didn't have this problem before either)
BAM record error: found spliced alignment without XS attribute
I find it strange because I do specify the '--outSAMstrandField intronMotif' parameter for STAR.

Why do I get this error even though I have specified the correct parameters for STAR?
Could this be the reason Cufflinks is hanging?

Thank you!
gogodidi is offline   Reply With Quote
Old 03-11-2015, 02:00 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,479
Default

Perhaps STAR 2.4.0f1 has a bug. Try something like this with one of the BAM files that's causing a problem:
Code:
samtools view foo.bam | awk -F '\t' '{if($6 ~"N") print $0}' | grep -v "XS:A:"
If something like that prints any lines then the file will cause an error in cufflinks. This would then likely be a bug in STAR.
dpryan is offline   Reply With Quote
Old 03-11-2015, 08:04 AM   #3
ccros
Junior Member
 
Location: Paris France

Join Date: Mar 2015
Posts: 3
Default

I have been using STAR and a de novo assembler called IsoSCM to identify 3'UTR. IsoSCM requested the XS tag to work, so I set the same --outSAMStrandField option and it went well.
I note that using "--outFilterIntronMotifs RemoveNonCanonical" is also recommended with Cufflinks in the STAR manual.
ccros is offline   Reply With Quote
Old 12-20-2015, 09:19 PM   #4
zhkzhou
Junior Member
 
Location: Beijing

Join Date: Nov 2011
Posts: 4
Default

4.2.3 Compatibility with Cufflinks
For unstranded RNA-seq data, Cuinks/Cu di require spliced alignments with XS strand attribute,which STAR will generate with --outSAMstrandField intronMotif option. As required, the XS strand attribute will be generated for all alignments that contain splice junctions. The spliced alignments that have unde ned strand (i.e. containing only non-canonical unannotated junctions)
will be suppressed.
zhkzhou is offline   Reply With Quote
Reply

Tags
cufflinks, star

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 06:31 AM.


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