SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Warning in cuffmerge genome1 General 2 05-15-2014 10:01 AM
reference transcripts missing after cuffmerge jake13 Bioinformatics 3 03-24-2014 06:41 PM
Cuffmerge Input? Carina1225 RNA Sequencing 0 01-28-2014 05:56 AM
Cuffmerge Warning: couldn't find fasta record for 'chr1_random' DonDolowy Bioinformatics 8 08-28-2013 01:02 AM
cuffmerge reports new transcripts? bjchen Bioinformatics 0 07-10-2012 02:35 PM

Reply
 
Thread Tools
Old 06-02-2014, 10:43 AM   #1
komalsrathi
Member
 
Location: Philadelphia, PA

Join Date: Aug 2013
Posts: 14
Default Warning in CuffDiff when using novel transcripts from CuffMerge as input

Hi everyone,

I am doing a differential gene expression test on RNASeq data from 64 Human Heart samples as follows:

Code:
cufflinks --library-type fr-firststrand -p 2 -o cufflinks_out <bam_file>
This gave me 64 transcripts.gtf files, one per sample. I created a list of the file names and saved it to assembly_list.txt

I then downloaded three GTF files from Gencode v19 : 1. Protein Coding Genes 2. Long Non-coding RNAs and 3. Pseudogenes. I merged the three gtf files, extracted out just the entries matching exon and CDS in the third column and used it as input for cuffmerge as follows.

Code:
cuffmerge -o cuffmerge_out -g gencode.v19.protein.ncRNA.pseudos_exonCDS.gtf -s hg19.fa -p 25 assembly_list.txt
and then I extracted the novel transcripts matching class_code="u":

Quote:
chr1 Cufflinks exon 16971271 16972173 . + . gene_id "XLOC_000205"; transcript_id "TCONS_00001942"; exon_number "1"; oId "CUFF.386.1"; class_code "u"; tss_id "TSS739";
chr1 Cufflinks exon 16972864 16972878 . + . gene_id "XLOC_000205"; transcript_id "TCONS_00001942"; exon_number "2"; oId "CUFF.386.1"; class_code "u"; tss_id "TSS739";
chr1 Cufflinks exon 20609243 20609496 . + . gene_id "XLOC_000241"; transcript_id "TCONS_00002229"; exon_number "1"; oId "CUFF.393.1"; class_code "u"; tss_id "TSS855";
chr1 Cufflinks exon 20611093 20611502 . + . gene_id "XLOC_000241"; transcript_id "TCONS_00002229"; exon_number "2"; oId "CUFF.393.1"; class_code "u"; tss_id "TSS855";
chr1 Cufflinks exon 33925677 33925957 . + . gene_id "XLOC_000396"; transcript_id "TCONS_00003679"; exon_number "1"; oId "CUFF.651.1"; class_code "u"; tss_id "TSS1384";
chr1 Cufflinks exon 33926256 33926832 . + . gene_id "XLOC_000396"; transcript_id "TCONS_00003679"; exon_number "2"; oId "CUFF.651.1"; class_code "u"; tss_id "TSS1384";
chr1 Cufflinks exon 56905467 56905577 . + . gene_id "XLOC_000616"; transcript_id "TCONS_00005969"; exon_number "1"; oId "CUFF.1131.1"; class_code "u"; tss_id "TSS2162";
chr1 Cufflinks exon 56910015 56910124 . + . gene_id "XLOC_000616"; transcript_id "TCONS_00005969"; exon_number "2"; oId "CUFF.1131.1"; class_code "u"; tss_id "TSS2162";

and used it in cuffdiff as follows:

Code:
cuffdiff -o cuffdiff_out --library-type fr-firststrand -L Control,HF --multi-read-correct --upper-quartile-norm -b hg19.fa --num-threads=20 novels.gtf <control1_bam,control2_bam,...> <hf1_bam,hf2_bam,...>
This is giving me a warning repeatedly and is not moving ahead:

Quote:
Warning: Your version of Cufflinks is not up-to-date. It is recommended that you upgrade to Cufflinks v2.2.1 to benefit from the most recent features and bug fixes (http://cufflinks.cbcb.umd.edu).
[09:32:46] Loading reference annotation and sequence.
[09:33:06] Inspecting maps and determining fragment length distributions.
Warning: Using default Gaussian distribution due to insufficient paired-end reads in open ranges. It is recommended that correct parameters (--frag-len-mean and --frag-len-std-dev) be provided.
Warning: Using default Gaussian distribution due to insufficient paired-end reads in open ranges. It is recommended that correct parameters (--frag-len-mean and --frag-len-std-dev) be provided.
Warning: Using default Gaussian distribution due to insufficient paired-end reads in open ranges. It is recommended that correct parameters (--frag-len-mean and --frag-len-std-dev) be provided.
Warning: Using default Gaussian distribution due to insufficient paired-end reads in open ranges. It is recommended that correct parameters (--frag-len-mean and --frag-len-std-dev) be provided.
I have never come across this kind of warning before. Any suggestions/help would be appreciated.

Thanks!
__________________
Komal Rathi
Bioinformatics Application Developer
University of Pennsylvania
komalsrathi is offline   Reply With Quote
Reply

Tags
cuffdiff, cuffmerge, novel transcripts

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:09 PM.


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