Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • reduced representation bisulfite seq MspI library prep - non target reads?

    Hi,

    Just started this again in the Epigenetics section (originally posted in Bioinformatics here), because I realized it's probably a better fit here - sorry for double-posting - redirected the other one to here.

    I've been wondering about the proportion of 'non-target' reads in a RRBS library, i.e. reads that I don't think are associated with a MspI restriction site but still end up in my dataset. I would really appreciate if anyone could share their experience with respect to this - some details about my data below - thanks!

    The library was prepared with MspI digestion with the NEXTFLEX® Bisulfite Library Prep Kit for Illumina, and sequenced paired-end on Illumina Miseq.

    I was expecting that most of my reads had remnants of the MspI cutsite at the start, so I wrote a little script to check that. The expectation is that forward reads should start with 'CGG' or 'TGG', if the first 'C' was methylated or not, respectively. For the reverse reads I was expecting mainly 'CAA', but also a few 'CGA' (unmethylated cytosines were used for end-repair, and I am assuming that bisulfite conversion rate is not 100%), as well as a few 'CAG' or 'CGG' (theoretically possible).
    So, ideally I should have the following combinations (forward - reverse) - please correct me if I am wrong:
    CGG-CAA
    TGG-CAA
    CGG-CGA
    TGG-CGA
    CGG-CAG
    TGG-CAG
    CGG-CGG
    TGG-CGG

    My expectations were met in the way that CGG-CAA and TGG-CAA was by far the most common, ~20% of read pairs and ~6% of read pairs, respectively. The other combinations account for less than 1% in total. What puzzles me though is that I only get ~30% read pairs that have any of the above patterns.

    I figured that perhaps some of the fragments I get after the MspI digestion are somehow sheared during the library prep and loose the cutsite on one side, so I relaxed the search criteria to also count read pairs that have the expected pattern only in the forward read OR in the reverse read. An additional ~28% of the read pairs fell into this category.

    Which leaves me with >40% of read pairs that don't have any of the expected patterns, neither in the forward nor the reverse read.

    So, my question: Is this normal? Sure there will be some sequencing error, but the read quality is very good overall and the error rate should be very low at the start of the reads, so that can't affect 40% reads. Am I missing something here?

    Thanks in advance for your help!

    cheers,
    Christoph

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
66 views
0 likes
Last Post seqadmin  
Working...
X