Hi
I perform targeted bisulfit sequencing of a specific gene from whole blood.
However, I get 3 different sequences in my results (when sequenced directly and also after cloning I get the same 3 sequences in my samples).
I have about 25 samples, after cloning I find 3 sequences in different samples and clones: the original bisulfit modified sequence of the gene (about 130bp) and 2 other sequences that are about the same length (one is +3bp the other -7bp) and hold both primers, but are nothing like the original sequence of the gene: the sequences start at the fwd-primer and ends at the rev-primer. They seem to be "hypermethylated" in so far that they have a lot of Cs in their sequence. Also, I have samples where I didnt find any clones with the "original" bisulfit treated sequence of the gene.
Just an example:
Original Sequence: AAGTCACGCG
Bisulfit treated original S.: AAGTCATGTG
Sequence 1: CACTCTAGC
Sequence 2: GGGTCTCAAAC
--> I find the last 3 (although all about 130bp, here just shorter as an example) in my samples
We have our theories, but are not "d'accord" in what is the best explenation, so I wanted to see what you think and if anyone maybe experienced the same thing?
Could these different sequences that I get with the gene-targeting primer be due to bisulfit treatment? Or is it really the problem that the primers bind somewhere else in the genome (although by blasting they were ok (of course in the untreated genome) and the sequences are all more or less the same length)?
Other explenations? Could allele-skewed methylation be an explenation for this?