Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • GTF modification for cufflinks

    Hey, does anyone have any pointers, advice, or experience on modifying GTF files for use with cufflinks??? (v 2.0.2)

    In the course of examining RNA-seq data and performing RNA-seq data analysis, an issue I've run into (using the "tuxedo" software/pipeline of tophat->cufflinks) is that tophat maps to apparent non-coding regions (possibly regulatory) but that cufflinks won't indicate FPKM expressions for the pileups! So a strategy we are trying, whose goal is to trigger cufflinks to tell FPKM expression values, is to either modify or create GTF annotation data and tell tophat/cufflinks to *not* try to find novel transcripts while using the created/modified GTF so that cufflinks might give FPKM values!

    One strategy we tried is to create a GTF with features/annotations corresponding to the regions of interest. Created as "pseudogene exons" (in columns 2 and 3), and using exsiting ensemble geneIDs, but custom transcript_ids we fed the GTF to cufflinks. When cufflinks program execution got to the "Loading Annotation" part (at the beginning of the run) it crashed with a segmentation fault! In the attribute column (#9), no information besides the gene_id and transcript_id was provided! cufflinks may have crashed because no gene_name was given. We really don't know however!

    Another strategy we are currently trying is to *modify* an existing GTF (from illumina/igenomes/ensemble) that *modify* work with cufflinks. This time, to capture regions upstream and downstream of genes, for each geneid, we modify the lowest start-value over all annotations by decreasing it by 1000 (to *hopefully* capture expressions of regions upstream). Similarly, we modify the highest end-value by increasing it by 1000 to *hopefully* capture expressions of regions downstream. This is currently going on now, so I don't know if the run will work, end successfully, and give us the expression/FPKM values/numbers we are looking for....

    Any pointers, advice, experience, knowledge, insight, etc. with GTF file tweaking for cufflinks would be appreciated!

    We are using tophat v2.0.4 and cufflinks v2.0.2 by the way.

    thanks

    -Eddie

  • #2
    Be very careful with all of this.

    There can be a lot of reasons for seeing reads outside of annotated genes that have nothing to do with real biology. They might be artifacts of your library prep. Even if they were real, you have no guarantee that your coverage is deep enough to accurately determine an FPKM value. If you don't have enough coverage to determine the length of the transcribed region, then the K part of the FPKM could lead to biased expression values.

    If you just extend GTF regions with an arbitrary number not informed by the biology of your system, you will create a lot of problems. You will be extending every gene by the same number, but it will not be the same relative to the actual length of the gene. This will lead to an underestimate of the expression of short genes to a higher degree than longer genes. Plus, what are you doing to ensure that your extensions don't create unwanted overlaps with other annotated genes?

    Make sure you have a good reason for looking at regions outside of annotated coding regions before you start modifying those annotations.

    Comment


    • #3
      Hi "pbluescript",

      Good points to be aware of. I realize now to regard any output with some skepticism.

      The goal at hand is to get FPKM values for the loci upstream and downstream of genes ; that's why we "extended" the genes (decreased lowest start value, increased highest end value).

      -Eddie
      Last edited by eddiesalinas; 11-07-2012, 09:08 AM.

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Best Practices for Single-Cell Sequencing Analysis
        by seqadmin



        While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
        06-06-2024, 07:15 AM
      • seqadmin
        Latest Developments in Precision Medicine
        by seqadmin



        Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

        Somatic Genomics
        “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
        05-24-2024, 01:16 PM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, 06-17-2024, 06:54 AM
      0 responses
      11 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 06-14-2024, 07:24 AM
      0 responses
      23 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 06-13-2024, 08:58 AM
      0 responses
      17 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 06-12-2024, 02:20 PM
      0 responses
      20 views
      0 likes
      Last Post seqadmin  
      Working...
      X