Announcement

Collapse
No announcement yet.

Application of PacBio based cDNA data in RNA-Seq analysis

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Application of PacBio based cDNA data in RNA-Seq analysis

    Hi Everyone

    Just wanted to know how people are using cDNA sequence data obtained from PacBio. Just curious to see and understand downstream applications in conjunction to Illumina cDNA data. Any relevant papers that you could point me to would help too.

    PS : cross posted on biostar to get more visibility and hopefully responses


    Thanks!
    -Abhi

  • #2
    Where's the URL for the biostar cross post?

    Comment


    • #3
      Hi maubp

      hers the URL for biostar corss posting.

      http://biostar.stackexchange.com/que...a-seq-analysis

      -Abhi

      Comment


      • #4
        I'm having trouble finding bioinformatics protocols for doing this too. The problem seems to be that most tools (scripture, PASA?, IsoInfer, Inchworm, IsoEM, NSMAP, MISO are most of what I've looked at) seem to depend on spliced alignments produced by, e.g. Tophat. But the only tools that align PacBio reads well (PB's BLASR, and BWA's bwasw command) don't produce spliced alignments, as far as I can tell.

        Let's add the terms isoform inference, alternative splicing, transcript, exon, splice junction, to improve search results for this thread.

        Comment


        • #5
          You might look at pacBioToCA and P_ErrorCorrection. They can be used to error correct PacBio reads using Illumina reads. Then the polished reads can be input to a wider variety of tools and you may be able to get spliced alignments.

          Comment


          • #6
            Originally posted by jnfass View Post
            I'm having trouble finding bioinformatics protocols for doing this too. The problem seems to be that most tools (scripture, PASA?, IsoInfer, Inchworm, IsoEM, NSMAP, MISO are most of what I've looked at) seem to depend on spliced alignments produced by, e.g. Tophat. But the only tools that align PacBio reads well (PB's BLASR, and BWA's bwasw command) don't produce spliced alignments, as far as I can tell.

            Let's add the terms isoform inference, alternative splicing, transcript, exon, splice junction, to improve search results for this thread.
            Hi, right now we're using GMAP to produce spliced alignments. It's a bit slow, but it is fairly sensitive for single pass raw pacbio reads (not ccs and not error corrected), and produces SAM output that is easily integrated into analysis pipelines. The only problem I've seen so far is some merging of exons.

            I'd be happy to hear of any analysis of GMAP.

            BLASR does not handle spliced alignment.

            Should I cross post on biostar?

            Comment

            Working...
            X