Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • fetch transcripts assembled by cufflinks

    Hi!

    I am an user of cufflinks. It's quite a good tool to estimate gene differential expression. However, our purpose is not only to test differentially expressed genes, but also to get these genes' sequences,which are assembled by cufflinks. We have used option -g in cufflinks and -r in cuffcompare.I try to fetch all the assembled transcripts by reading the locus of each assembled transcripts' exon from combined.gtf file produced by cuffcompare,then using samtools mpileup to parse whether a locus has a different base
    in mapped reads comparing to reference base of this locus, and finally joining these exons together to get full length sequences of assembled transcripts with IDs like "TCON...". However, I find that in combined.gtf,some exons' locus doesn't have any reads mapped at all, and if I continue the method I described before, I just wonder if I should use reference base in this locus, because this locus is actually not detected by real reads.


    I have read the manual of cufflinks in webset. I think the point is that in RABT model, cufflinks uses faux-reads to help assembling transcripts. It is quite helpful to detect new isoforms, However, when you try to get the real sequences of these isoforms, things became confusing. Do you have any suggestions to separate the real sequenced transcripted assembled by real reads from those reference trancripts assembled by faux-reads? In another word, how can I get the sequnces of each assembled transcripts from cufflinks?
    Thanks for help!


    asling

  • #2
    Any suggestion?

    Comment


    • #3
      I am facing the same problem: how can extract the assembled transcripts from RNA-seq data, NOT from reference genome.! Does anyone have any suggestions? I really appreciate.!

      Comment


      • #4
        Originally posted by seqfan05 View Post
        I am facing the same problem: how can extract the assembled transcripts from RNA-seq data, NOT from reference genome.! Does anyone have any suggestions? I really appreciate.!
        Hi seqfan05

        Have you sorted out the problem as I am facing the same problem and there is not much help reading posts.

        regards

        Comment


        • #5
          I'm having the same problem. I have RNAseq data that i analysed using the Tuxedo protocol. If i use the custom GTF file generated by cufflinks/cuffmerge then i can see there are lots of new splice variants but i can't see their sequence or the gene they most resemble. In addition, when you have multiple time points or treatment each containing hundreds or possible even thousands of new splicevarients how are you supposed to investigate them or describe these changes in splicevarient abundance? Splicegrapher looks like a nice way to investigate splicevarients but as far as i can see it only allows investigation of a few genes at a time http://splicegrapher.sourceforge.net/.

          Comment


          • #6
            Originally posted by Richard Barker View Post
            I'm having the same problem. I have RNAseq data that i analysed using the Tuxedo protocol. If i use the custom GTF file generated by cufflinks/cuffmerge then i can see there are lots of new splice variants but i can't see their sequence or the gene they most resemble. In addition, when you have multiple time points or treatment each containing hundreds or possible even thousands of new splicevarients how are you supposed to investigate them or describe these changes in splicevarient abundance? Splicegrapher looks like a nice way to investigate splicevarients but as far as i can see it only allows investigation of a few genes at a time http://splicegrapher.sourceforge.net/.
            I am sorry, I still can not find a way to address it. Can anybody help?

            Comment


            • #7
              here is the link that addresses this problem
              Application of sequencing to RNA analysis (RNA-Seq, whole transcriptome, SAGE, expression analysis, novel organism mining, splice variants)

              Comment

              Latest Articles

              Collapse

              • seqadmin
                Recent Developments in Metagenomics
                by seqadmin





                Metagenomics has improved the way researchers study microorganisms across diverse environments. Historically, studying microorganisms relied on culturing them in the lab, a method that limits the investigation of many species since most are unculturable1. Metagenomics overcomes these issues by allowing the study of microorganisms regardless of their ability to be cultured or the environments they inhabit. Over time, the field has evolved, especially with the advent...
                09-23-2024, 06:35 AM
              • seqadmin
                Understanding Genetic Influence on Infectious Disease
                by seqadmin




                During the COVID-19 pandemic, scientists observed that while some individuals experienced severe illness when infected with SARS-CoV-2, others were barely affected. These disparities left researchers and clinicians wondering what causes the wide variations in response to viral infections and what role genetics plays.

                Jean-Laurent Casanova, M.D., Ph.D., Professor at Rockefeller University, is a leading expert in this crossover between genetics and infectious...
                09-09-2024, 10:59 AM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, Yesterday, 04:51 AM
              0 responses
              8 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 10-01-2024, 07:10 AM
              0 responses
              13 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 09-30-2024, 08:33 AM
              0 responses
              18 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 09-26-2024, 12:57 PM
              0 responses
              16 views
              0 likes
              Last Post seqadmin  
              Working...
              X