Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Trinity transcript naming

    I don't know if I am missing something obvious but I can't find an explanation anywhere as to how Trinity RNA-seq names the transcripts it derives.

    I've been noticing that when I annotate with BLAST+ most of the transcripts that start with the same number (comp#_) are annotating as the same gene (isoforms?). They will then all have different ending numbers (_seq#) and some of them have different numbers int he middle (_c#_), for example (comp2335_c1_seq1, comp2335_c1_seq2, comp2335_c1_seq3).

    Are these names derived sequentially by inchworm (comp#) chrysalis (c#) then butterfly (seq#)?

    Can anyone verify whether these would be indeed considered isoforms of single genes within a comp# group or is there something else at work here?

    Thanks,
    -Mike

  • #2
    You're almost correct. The names are indeed based on the algorithm, but it's not quite as simple as the inchworm number, then the chrysalis number, then the butterfly number.

    This is the basic algorithm of Trinity, assuming that you've already built a de Bruijn assembly graph (which is identical to a (k+1)-mer count):
    1. Eagerly extract contigs from the de Bruijn graph. These contigs may or may not have any relationship to "real" transcripts, they're just whatever long contiguous paths happen to be found in the graph.
    2. Find reads which justify clustering/joining ("welding" in Trinity terminology) these contigs together. A set of contigs which are believed to belong together is called a "component".
    3. Align reads to components. For each read, decide which component it's most likely to belong to.
    4. For each component, treat the reads which map to that component as a separate assembly problem. This involves constructing a new (smaller) de Bruijn graph from only those reads which belong to a component.

    The output of Inchworm is the set of contigs. The output of Chrysalis is the set of components, plus the reads which are called as belonging to those components. The output of Butterfly is the called transcripts for each component.

    When Butterfly rebuilds a graph for each component and does cleanup, it sometimes finds that the resulting graph is disconnected. This is usually because the inital contigs discovered by Inchworm were not "real" contigs. Each connected component in the rebuilt graph is called a "subcomponent".

    So the "comp" is the component, "c" is the subcomponent, and "seq" is the extracted sequence from the subcomponent.

    Trinity does not reason in terms of genes, loci and alternative splicing events. It solves a graph problem, though of course the heuristics are tuned to the needs of biology. So while it's highly likely that all the isoforms of a given gene belong to the same subcomponent, you shouldn't assume that a subcomponent is a gene.
    Last edited by Pseudonym; 04-15-2012, 05:05 PM.
    sub f{($f)=@_;print"$f(q{$f});";}f(q{sub f{($f)=@_;print"$f(q{$f});";}f});

    Comment


    • #3
      Nice explaination Pseudonym...



      Originally posted by Pseudonym View Post
      You're almost correct. The names are indeed based on the algorithm, but it's not quite as simple as the inchworm number, then the chrysalis number, then the butterfly number.

      This is the basic algorithm of Trinity, assuming that you've already built a de Bruijn assembly graph (which is identical to a (k+1)-mer count):
      1. Eagerly extract contigs from the de Bruijn graph. These contigs may or may not have any relationship to "real" transcripts, they're just whatever long contiguous paths happen to be found in the graph.
      2. Find reads which justify clustering/joining ("welding" in Trinity terminology) these contigs together. A set of contigs which are believed to belong together is called a "component".
      3. Align reads to components. For each read, decide which component it's most likely to belong to.
      4. For each component, treat the reads which map to that component as a separate assembly problem. This involves constructing a new (smaller) de Bruijn graph from only those reads which belong to a component.

      The output of Inchworm is the set of contigs. The output of Chrysalis is the set of components, plus the reads which are called as belonging to those components. The output of Butterfly is the called transcripts for each component.

      When Butterfly rebuilds a graph for each component and does cleanup, it sometimes finds that the resulting graph is disconnected. This is usually because the inital contigs discovered by Inchworm were not "real" contigs. Each connected component in the rebuilt graph is called a "subcomponent".

      So the "comp" is the component, "c" is the subcomponent, and "seq" is the extracted sequence from the subcomponent.

      Trinity does not reason in terms of genes, loci and alternative splicing events. It solves a graph problem, though of course the heuristics are tuned to the needs of biology. So while it's highly likely that all the isoforms of a given gene belong to the same subcomponent, you shouldn't assume that a subcomponent is a gene.

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Strategies for Sequencing Challenging Samples
        by seqadmin


        Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
        03-22-2024, 06:39 AM
      • seqadmin
        Techniques and Challenges in Conservation Genomics
        by seqadmin



        The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

        Avian Conservation
        Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
        03-08-2024, 10:41 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, Yesterday, 06:37 PM
      0 responses
      10 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, Yesterday, 06:07 PM
      0 responses
      10 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-22-2024, 10:03 AM
      0 responses
      51 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-21-2024, 07:32 AM
      0 responses
      67 views
      0 likes
      Last Post seqadmin  
      Working...
      X