This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Finding telomere sequences in exome-sequencing data

Hi,

The exome sequencing output typically contains some fraction (typically 10–50%) of sequence that is off-target. I was wondering how I may check if the off-target sequences in a given exome-sequencing data (BAM format) contains telomere sequences or not (For instance, by checking repeats of TTAGGG hexamers).

sequencing next-gen
• 2,704 views • link
updated by jotan jotan • written by marki marki

Have you tried making a longish sequence of the hexamers and aligning to that? Presumably that'd work. Alternatively, I wouldn't be surprised if the kmer part that FastQC does would show this.

• 2 views • link
written by Devon Ryan Devon Ryan

I would like someone to confirm this, but I think telomeres are masked in genome files and as such no reads will map to these in your bam file. However, there might be unmapped reads in your bam file corresponding to the telomeric sequences.

• 1 views • link
written by WouterDeCoster WouterDeCoster

Yup, they're almost always hardmasked. This is true of most regions of constitutive heterochromatin.

• 1 views • link
written by Devon Ryan Devon Ryan

Thank you WouterDeCoster and Devon Ryan.

• 1 views • link
written by marki marki

1 answer

I like Repeat Enrichment Estimator. It includes telomere TTAGGG as part of the standard database.

As an aside, telomeric transcripts are pretty low abundance and not easy to detect. I've found only low levels of telomeric sequences even when using only rRNA depleted preps. Depending on what you are trying to do, you might want to consider looking at other repeats.

• 1 views • link
written by jotan jotan

OP is talking about exome sequencing, not transcripts ;)

• 1 views • link
written by WouterDeCoster WouterDeCoster

Whoops. Sorry, my bad!

• 1 views • link
written by jotan jotan

Log in to answer this question.

More posts like this