Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Please help: imperfect reference genome/get consensus on genome/read alignment?

    I an doing RNA-seq on two closely related plant species (~98% similar in transcript region), one of which is sequenced.

    My problem comes when I align the reads from the unsequenced species to the sequenced genome. Most of the reads align (at 2 bp mismatch level), however some are being rejected due to 3 or 4 bp mismatches. I would like to improve on the alignability of reads/reference genome to ensure the observed differences in expression between species are true (not due to differences in read alignability).

    My idea is as follows:
    1. Align the reads to the reference (allow ~2 bp mismatch)
    2. Use the alignment to generate a consensus sequence, ie where there are SNPs between the reference and the reads, change the base on the chromosome to reflect that of the reads. Use stringent parameters so only true SNPs are 'corrected'.
    3. Realign the reads to the modified genome (allow ~2 bp mismatch). Now reads with 2 bp mismatches should align 'perfectly' and reads with up to 4 bp mismatches will align.


    Can anyone recommend programs or a strategy to do this, or even comment on the validity of the approach?

    Any feedback/contribution would be much appreciated! Thank you!!

    Note: the reads are very short: 40 bp.

  • #2
    I had the similar idea when facing a somewhat similar situation. In theory, this would work, but in practice, I found the transcript levels quantified from the 'reconstructed' genome are biased towards lower end. I suspected that these didn't reflect the true expression levels in the 'reconstructed' species, skewed towards underestimation. I haven't thought about the exact causes, but I think it's related to the 'reconstruction' process and the heterogeneity of genomic divergence along chromosomes between two species. So in the end, I just mapped to the same reference, but allowing extra mismatches for the divergent species. Another way to look into this is just going ahead and to do a de novo assembly for the divergent species, but that opens a can of worms, too. And it's a complete different ball game.
    Last edited by macrowave; 08-22-2011, 06:51 AM.

    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
    8 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, Yesterday, 06:07 PM
    0 responses
    8 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-22-2024, 10:03 AM
    0 responses
    49 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