Cannot find any mention of reference gff/gtf annotations in cuffmerge output

SEQanswers June Challenge Has Begun!

The competition has begun! We're giving away a $50 Amazon gift card to the member who answers the most questions on our site during the month. We want to encourage our community members to share their knowledge and help each other out by answering questions related to sequencing technologies, genomics, and bioinformatics. The competition is open to all members of the site, and the winner will be announced at the beginning of July. Best of luck!

For a list of the official rules, visit (https://www.seqanswers.com/forum/sit...wledge-and-win)

See less

anna_vt

Junior Member

Join Date: Jun 2009
Posts: 5

Cannot find any mention of reference gff/gtf annotations in cuffmerge output

11-13-2012, 04:55 AM

Hi,

I have been using tophat and cufflinks to analyse some RNA seq data, but cuffmerge is not behaving as I would expect.

For all steps of the analysis I have used a reference gff3 file for RABT.

For tophat I specifed:

Code:

-G <gff file> --transcriptome-index <transcriptome prefix>

and on subsequent usage just

Code:

--transcriptome-index <transcriptome prefix>

For cufflinks I specified

Code:

--GTF-guide <gff file>

All the fpkm_tracking and gtf output files contain annotations from the original gtf file, including from transcripts with very high coverage so not just faux reads.

When using cuffmerge for the 16 samples I use the following command:

Code:

cuffmerge <cufflinks gtf list file> -o cuffmerge -g <gff file> -p 17

In none of the output files is there any mention of the annotations from the gff file.
I thought that the problem was that it is a gff file and not a gtf file, so I converted using

Code:

gffread <gff file> -T -o <gtf file>

Code:

cuffmerge <cufflinks gtf list file> -o cuffmerge -g <gtf file> -p 17

However the output is the same. There is also no mention of the gtf file in the log file for either gff3 or gtf:

Code:

/software/additional/bin/cuffmerge <path>/sam/tophat/cufflinks/gtfFileList.txt -o cuffmerge -g <path>/data/abH97_filt.gffread.gtf -p 17
gtf_to_sam -F <path>/sam/tophat/1_1/cufflinks/transcripts.1_1.gtf ./merged_asm/tmp/gtf2sam_filenEmb9h
gtf_to_sam -F <path>/sam/tophat/1_2/cufflinks/transcripts.1_2.gtf ./merged_asm/tmp/gtf2sam_fileeDtnvd
gtf_to_sam -F <path>/sam/tophat/1307_M5/cufflinks/transcripts.1307_M5.gtf ./merged_asm/tmp/gtf2sam_fileQhIpNd
gtf_to_sam -F <path>/sam/tophat/1_3/cufflinks/transcripts.1_3.gtf ./merged_asm/tmp/gtf2sam_fileAHuzXi
gtf_to_sam -F <path>/sam/tophat/2_1/cufflinks/transcripts.2_1.gtf ./merged_asm/tmp/gtf2sam_fileZw6YGs
gtf_to_sam -F <path>/sam/tophat/2_2/cufflinks/transcripts.2_2.gtf ./merged_asm/tmp/gtf2sam_file87N56G
gtf_to_sam -F <path>/sam/tophat/2_3/cufflinks/transcripts.2_3.gtf ./merged_asm/tmp/gtf2sam_file6ipSj0
gtf_to_sam -F <path>/sam/tophat/A151_M5/cufflinks/transcripts.A151_M5.gtf ./merged_asm/tmp/gtf2sam_file3SWsto
gtf_to_sam -F <path>/sam/tophat/F25_2_1/cufflinks/transcripts.F25_2_1.gtf ./merged_asm/tmp/gtf2sam_fileC2SotR
gtf_to_sam -F <path>/sam/tophat/F25_2_3/cufflinks/transcripts.F25_2_3.gtf ./merged_asm/tmp/gtf2sam_filelpirkp
gtf_to_sam -F <path>/sam/tophat/F25_2_4/cufflinks/transcripts.F25_2_4.gtf ./merged_asm/tmp/gtf2sam_fileJX9sR1
gtf_to_sam -F <path>/sam/tophat/F25_2_M5/cufflinks/transcripts.F25_2_M5.gtf ./merged_asm/tmp/gtf2sam_filezkS48I
gtf_to_sam -F <path>/sam/tophat/F25_2_M6/cufflinks/transcripts.F25_2_M6.gtf ./merged_asm/tmp/gtf2sam_file68Mb8u
gtf_to_sam -F <path>/sam/tophat/I142_4/cufflinks/transcripts.I142_4.gtf ./merged_asm/tmp/gtf2sam_file6CwbQl
gtf_to_sam -F <path>/sam/tophat/I142_M10/cufflinks/transcripts.I142_M10.gtf ./merged_asm/tmp/gtf2sam_fileeikjVg
gtf_to_sam -F <path>/sam/tophat/I142_M9/cufflinks/transcripts.I142_M9.gtf ./merged_asm/tmp/gtf2sam_filetTRVKg
sort -k 3,3 -k 4,4n --temporary-directory=./merged_asm//meta_asm_tmp/ ./merged_asm/tmp/gtf2sam_filenEmb9h ./merged_asm/tmp/gtf2sam_fileeDtnvd ./merged_asm/tmp/gtf2sam_fileQhIpNd ./merged_asm/tmp/gtf2sam_fileAHuzXi ./merged_asm/tmp/gtf2sam_fileZw6YGs ./merged_asm/tmp/gtf2sam_file87N56G ./merged_asm/tmp/gtf2sam_file6ipSj0 ./merged_asm/tmp/gtf2sam_file3SWsto ./merged_asm/tmp/gtf2sam_fileC2SotR ./merged_asm/tmp/gtf2sam_filelpirkp ./merged_asm/tmp/gtf2sam_fileJX9sR1 ./merged_asm/tmp/gtf2sam_filezkS48I ./merged_asm/tmp/gtf2sam_file68Mb8u ./merged_asm/tmp/gtf2sam_file6CwbQl ./merged_asm/tmp/gtf2sam_fileeikjVg ./merged_asm/tmp/gtf2sam_filetTRVKg > ./merged_asm/tmp/mergeSam_file1gh4tl
cufflinks -o ./merged_asm/ -F 0.05 -q --overhang-tolerance 200 --library-type=transfrags -A 0.0 --min-frags-per-transfrag 0 --no-5-extend -p 1 ./merged_asm/tmp/mergeSam_file1gh4tl
cuffcompare -o tmp_meta_asm -C -G  ./merged_asm//transcripts.gtf
cuffcompare -o tmp_meta_asm -C -G  ./merged_asm//merged.gtf

And nothing is prefixed / goes to the folder cuffmerge that i specifed with the -o flag, it only goes to the merged_asm folder.

Can anyone tell me what I am doing wrong or if this behaviour is as expected?
Version is: cufflinks v2.0.2

Many Thanks
Anna

Tags: cufflinks 2.0.2, cuffmerge, gff, gff3, gtf

Previous template Next

A Review of Next-Generation Sequencing Methods for Studying Epigenetics—Part 2

by seqadmin

In the first segment of our review of epigenetic sequencing methods, we covered methylation detection, chromatin accessibility, and protein-DNA interactions. This final part of our review will break outside of the standard analysis and cover some influential multidimensional methods.

Spatial Epigenomics
“Epigenetic assays, whether conducted at the single-cell or bulk level, have played a crucial role in unraveling the underlying mechanisms of diseases and...
- Channel: Articles
Today, 06:48 AM
A Review of Next-Generation Sequencing Methods for Studying Epigenetics—Part 1

by seqadmin

Developments in sequencing technologies and methodologies have transformed the field of epigenetics, giving researchers a better way to understand the complex world of gene regulation and heritable modifications. This article explores some of the diverse sequencing methods employed in the study of epigenetics, ranging from classic techniques to cutting-edge innovations while providing a brief overview of their processes, applications, and advances.

Methylation Detect...
- Channel: Articles
05-31-2023, 10:46 AM
Differential Expression and Data Visualization: Recommended Tools for Next-Level Sequencing Analysis

by seqadmin

After covering QC and alignment tools in the first segment and variant analysis and genome assembly in the second segment, we’re wrapping up with a discussion about tools for differential gene expression analysis and data visualization. In this article, we include recommendations from the following experts: Dr. Mark Ziemann, Senior Lecturer in Biotechnology and Bioinformatics, Deakin University; Dr. Medhat Mahmoud Postdoctoral Research Fellow at Baylor College of Medicine;...
- Channel: Articles
05-23-2023, 12:26 PM

Topics	Statistics	Last Post
Analysis of West Nile Virus Explores Transmission Dynamics in Animal Species by seqadmin Started by seqadmin, 06-07-2023, 07:14 AM	0 responses 11 views 0 likes	Last Post by seqadmin 06-07-2023, 07:14 AM
BabySeq Project: Unveiling Medically Actionable Variants by seqadmin Started by seqadmin, 06-06-2023, 01:08 PM	0 responses 11 views 0 likes	Last Post by seqadmin 06-06-2023, 01:08 PM
Genetic Discoveries Transform Our Understanding of Primate Diversity and Behavior by seqadmin Started by seqadmin, 06-01-2023, 08:56 PM	0 responses 164 views 0 likes	Last Post by seqadmin 06-01-2023, 08:56 PM
Deep Sequencing Unearths Novel Genetic Variants: Enhancing Precision Medicine for Vascular Anomalies by seqadmin Started by seqadmin, 06-01-2023, 07:33 AM	0 responses 299 views 0 likes	Last Post by seqadmin 06-01-2023, 07:33 AM

Header Leaderboard Ad

Cannot find any mention of reference gff/gtf annotations in cuffmerge output

Announcement

SEQanswers June Challenge Has Begun!

Cannot find any mention of reference gff/gtf annotations in cuffmerge output

Latest Articles

ad_right_rmr

News