This is a test version of Biostars. For the public version, visit https://www.biostars.org.
RRBS (Reduced Representation Bisulfite Sequencing)

I am working with single-end RRBS data, but I do not have information about the restriction enzyme used during library preparation. Trim Galore offers the --rrbs option, which is tailored for data generated using the MspI enzyme. Is there a way to confirm from my raw FASTQ files whether MspI was used for restriction digestion?

trim rrbs galore

1 answer

You should have enrichment of specific restriction enzyme motifs at the start of your reads. If your reads often start with something like CCGG or CGG, then it's likely MspI.

Beyond that, guys, how is it that there is no documentation about your experiments. This can't be. Publication needs an adequate methods text...

Whilst this may be true, I second ATpoint, you should try to not work with collaborators that are this disorganised. If you don't even know the restriction enzyme, how are you going to publish this? All this hard work and potentially interesting results will likely not get past preprint.

It is possible that the data is from a public source (e.g. SRA) and thus the details on methodology may be unclear.

Fair, I didn't really think of that in the moment. But it's a similar issue, if the methods are so unclear, you can't really be that confident of any results since you can't escape the possibility of an undescribed treatment causing all your nice tidy results.

Log in to answer this question.