Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Merging reads mapped on genome and CDS (SOLID data)

    Hi everyone,


    We are interested in discovering fusion genes in RNA-seq data from a
    certain tumoral type, and FusionSeq seems to be a tool of choice for
    this particular topic. Our RNA-seq data comes from the SOLID sequencing
    platform, it is paired-end with 50bp-long reads on the 3' side and 25
    bp-long reads on the 5' side.

    Our main concern is the choice of the mapping strategy. We use Bfast+BWA, that has the advantage of being sensitive and capable of handling very short reads. Here is our mapping strategy : we aligned reads on the reference genome, then we extracted the non-aligned reads, and we aligned them on a library of all known human CDS. We then want to use FusionSeq for the search of fusion genes. As far as we understood FusionSeq takes as an input reads aligned on a genomic reference under the MRF format. Our concern here is the following : how to merge these two mappings (genome and CDS) into a single one, that we can then feed to FusionSeq ? Do we need to convert the coordinates of reads aligned on the CDS to genomic coordinates ? And if so, how ?

    Thank you by advance for your help.

    Jonas and Jaydutt

  • #2
    Jonas and Jaydut,

    We use SOLiD extensively here at EdgeBio and finds that when aligning transcriptome data, even more than genomic data, the alignments really matter. Just aligning with bfast or bowtie without taking into account the complexities of transcriptome data won't yield optimal results. I suggest Tophat or novoalign (see below) but we find Tophat doesn't map solid data extremely well due to its error profiles. It is more tuned to Illumina data.

    We have consolidated around novoalign for our alignments for genomic and transcriptomic data, and novoalign has recently released a recipe to align your RNASeq reads to the whole genome + splice junctions. It shows you how to convert to genomic coordinates as you ask, and leverages trusted R packages such as Useq and DESeq. The result will be sorted alignments in BAM format against the transcriptome, genome and splice junctions. How this is fed into fusion detection software should be a trivial conversion. We are in the process of testing this with an internal project of 32 samples.

    If you use the recipe, I would be interested in comparing notes.
    Justin H. Johnson | Twitter: @BioInfo | LinkedIn: http://bit.ly/LIJHJ | EdgeBio

    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