Go Back   SEQanswers > Applications Forums > Epigenetics

Similar Threads
Thread Thread Starter Forum Replies Last Post
reduced representation bisulfite seq MspI library prep - non target reads? chrisi.hahni Bioinformatics 1 03-26-2019 03:35 AM
Reduced Representation Bisulfite Sequencing for Methylation Analysis audioPL Illumina/Solexa 10 03-18-2013 07:10 PM

Thread Tools
Old 03-26-2019, 03:30 AM   #1
Junior Member
Location: Graz, Austria

Join Date: Mar 2019
Posts: 3
Default reduced representation bisulfite seq MspI library prep - non target reads?


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:

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!

chrisi.hahni is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 12:58 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO