How to get karyotype information for a certain section of genome.
2
0
Entering edit mode
7.5 years ago
Chen ★ 1.1k

Hi all!

I have some sections of genome location such as chr16:100-300, chr17: 1000-2000. I want to know how can I get the corresponding karyotype information. As I have large amount of genome location sections, I want to whole file of genome location to corresponding karyotype.

For instance, chr16:100-300 may corresponding to 16q1, which chr17: 1000-2000 may corresponding to 17q2 etc.

The location information is from hg19.

karyotype • 2.0k views
ADD COMMENT
0
Entering edit mode
7.5 years ago

That information is available in the cytoBandIdeo table from the UCSC table browser (group: "Mapping and Sequencing", track: "Chromosome Band (Ideogram)").

ADD COMMENT
0
Entering edit mode

Thank you, how can I get the whole file of genome location to karyotype? as I have many sections, it will not be realistic to query them one by one in UCSC browser.

ADD REPLY
0
Entering edit mode

You can just download the whole file (note the "Get Output" button) from the table browser or just get it via FTP. There are a variety of way to match a given region to that file, then. If you're in R already, you can findOverlaps(), or you could use bedtools or bedops (to give just two examples) from the command line.

ADD REPLY
0
Entering edit mode
7.5 years ago

I found Pierre's answer useful for retrieving karyotype data as a file: A: Automating Ideogram Generation For Various Organisms

After that, I would convert to BED and use bedops to do a set operation on regions-of-interest:

$ echo -e "chrN\tA\tB" | bedops --element-of -1 someKaryotype.bed -
ADD COMMENT

Login before adding your answer.

Traffic: 2063 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6