This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Identification of Genes spanning a QTL

Hi Everyone,

I have a list of QTLs. The chromosome number, start and stop regions are given. I have a list of genes which I want to map onto the regions. I have obtained the locations of the genes from BioMart. I wanted to find out the genes lying in the QTL regions. Is it just a trivial mapping using a script or do I need to consider some other conditions to identify the genes in the QTL region ? Any suggestion will be greatly appreciated.

qtl
• 2,959 views • link
updated by Ram Ram • written by pixie@bioinfo pixie@bioinfo

Bedtools is a great option as suggested by Mark. If you have no programming background, you can use UCSC Table Browser. You can select your organism of interest, select Gene and gene predictions and define multiple locations or list of your QTL locations by clicking on define regions. You can then download data in GTF or bed format.

• 0 views • link
updated by Ram Ram • written by Ashutosh Pandey Ashutosh Pandey

1 answer

If the QTL coli and gene loci are both in bed format (and the genome build is the same), then you can use bedtools intersect to determine the genes that lie in those regions.

bedtools intersect -wa -a genes.bed -b QTLs.bed

Bedtools can also be accessed via Galaxy if you are on a windows machine.

• 0 views • link
updated by Ram Ram • written by mark.ziemann mark.ziemann

Log in to answer this question.

More posts like this