Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • eddiesalinas
    replied
    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.

    Leave a comment:


  • pbluescript
    replied
    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.

    Leave a comment:


  • eddiesalinas
    started a topic GTF modification for cufflinks

    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

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
22 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