How to find genomic geographic commonalities between sequences
1
0
Entering edit mode
6.3 years ago

OK so this may be broad but I have a load of sequencing associated with a clinical outcome I'm interested in. I would like to find a way of investigating the commonalities of these sequences in terms of where they are in the genome. Specifically I would like to look at things like: Distance to centromere, distance to telomere,  proximity to GC and AT rich regions, association with tertiary DNA structures, association with fragile sites. Are there tools that can give me a 'your sequence is associated with these genomic areas' type tool or do I have to do each seperately?

sequencing Assembly genome next-gen • 1.5k views
2
Entering edit mode
6.3 years ago

First, assuming you are working with human data, generate a file called centromeres.bed that contains genomic ranges for centromeres:

$wget -qO- http://hgdownload.cse.ucsc.edu/goldenPath/hg19/database/cytoBand.txt.gz | gunzip -c | grep acen > centromeres.bed Next, say you have a file containing sequencing reads called reads.bam. We convert it to a sorted BED file using convert2bed: $ convert2bed -i bam < reads.bam > reads.bed

Finally, calculate the signed distance between each read and its nearest centromere with closest-features:

\$ closest-features --closest --dist reads.bed centromeres.bed > distances_of_reads_to_closest_centromeres.bed

To now answer your other questions, given your reads (now in sorted BED form), you can think about using UCSC and other data sources to generate BED files that contain telomere regions, GC- and AT-enriched regions of interest, genomic regions that associate with tertiary DNA structure (e.g., ChIP-seq regions or motif binding sites), and regions associated with fragile sites. You could then do set and statistical operations with your reads against these regions using bedops, bedmap, closest-features and other tools in BEDOPS to help answer these and similar questions.