Go Back   SEQanswers > Applications Forums > RNA Sequencing

Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem using the Trinity script with RSEM gwilymh Bioinformatics 3 03-31-2015 08:42 AM
Selecting most abundant transcript per gene (Trinity, RSEM) gevielr RNA Sequencing 5 09-29-2014 08:00 PM
bowtie vs bowtie2 for RSEM analysis trinity assembled transcripts horvathdp Bioinformatics 1 12-05-2013 12:47 PM
more problems with RSEM in Trinity build horvathdp Bioinformatics 0 11-22-2013 01:59 PM
question about the RSEM module within trinity ataraxia RNA Sequencing 0 08-20-2013 10:31 PM

Thread Tools
Old 09-13-2019, 02:38 AM   #1
Junior Member
Location: Italy

Join Date: Sep 2019
Posts: 5
Default Trinity RSEM Bowtie2 align_and_estimate_abundance script error

I'm running a RNA-Seq differential gene expression pipeline analysis. I previously filtered my input sequences, then de-novo assembled my transcriptome through Trinity, and now, through a Trinity-provided script, I'm trying to evaluate and estimate transcript abundances through the different samples. This way, I'll be able to use the resulting data to run a differential gene expression analysis, through a R package (EBSEQ) as a part of the RSEM suite.

I got to the point I'm running the script provided by Trinity, with the following parameters:

home/genomica/DATA/Software/Trinity/trinityrnaseq-2.8.5/util/ --transcripts /home/genomica/DATA/Lattuga/denovo_trinity_denovo_filtered_dump_completo/Trinity.fasta --seqType fq --single --/home/genomica/DATA/Lattuga/DRR057059_1_clipped.fastq --thread_count 60 --est_method RSEM --aln_method bowtie2 --trinity_mode --debug --output_dir /home/genomica/DATA/Lattuga/denovo_trinity_denovo_filtered_dump_completo/ --prep_reference

where Trinity.fasta is my de-novo assembled transcriptome

and DRR057059_1_clipped.fastq one of the libraries used to assemble it (just running the script on one file for now)

Am I doing anything wrong?

I get this debug message after the script creates several files you can see below:

CMD: set -o pipefail && bowtie2 --no-mixed --no-discordant --gbar 1000 --end-to-end -k 200 -q -x //home/genomica/DATA/Lattuga/denovo_trinity_denovo_filtered_dump_completo/Trinity.fasta.bowtie2 -U //--/home/genomica/DATA/Lattuga/DRR057059_1_clipped.fastq -p 60 | samtools view -@ 60 -F 4 -S -b | samtools sort -@ 60 -n -o bowtie2.bam
stat: No such file or directory
Warning: Could not open read file "//--/home/genomica/DATA/Lattuga/DRR057059_1_clipped.fastq" for reading; skipping...
Error: No input read files were valid
(ERR): bowtie2-align exited with value 1
Error, cmd: set -o pipefail && bowtie2 --no-mixed --no-discordant --gbar 1000 --end-to-end -k 200 -q -x //home/genomica/DATA/Lattuga/denovo_trinity_denovo_filtered_dump_completo/Trinity.fasta.bowtie2 -U //--/home/genomica/DATA/Lattuga/DRR057059_1_clipped.fastq -p 60 | samtools view -@ 60 -F 4 -S -b | samtools sort -@ 60 -n -o bowtie2.bam died with ret: 256 at home/genomica/DATA/Software/Trinity/trinityrnaseq-2.8.5/util/ line 733.

Something seems to be wrong with a Trinity.fasta.bowtie2 that the script can't find (there is one in the working directly but it seems to be set in different parts)

The file DRR057059_1_clipped.fastq is a legit file, used previously to assemble the transcriptome, why shouldn't it be able to open it?
Maybe it's something about the synthax?)

(I don't think it's about the path, I had set the bowtie2 path correctly previously)

Thanks in advance for your help, I'm stuck here, and it could be something basic I'm missing.
Light92 is offline   Reply With Quote


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:34 AM.

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