Question: Given QTLs and genome annotation, quickly list genes closest to QTLs
0
gravatar for michael.nagle
3.1 years ago by
michael.nagle100
michael.nagle100 wrote:

I have GWAS data with many associated loci. I'm looking for a way of listing all of the closest genes (and their positions) without clicking through in a genome browser for every locus.

I wish to provide an input of the genome annotations (in .gff3 or another format) and a list of every associated loci (Chr and base position) and get an output of several of the closest genes and their start and stop positions.

I assume there's already some off-the-shelf package I can use for this, or at least some code.

Thanks for the help!

rna-seq genomics gwas genome • 616 views
ADD COMMENTlink modified 3.1 years ago by pfs280 • written 3.1 years ago by michael.nagle100
0
gravatar for pfs
3.1 years ago by
pfs280
USA/Boston
pfs280 wrote:

Have you looked into ANNOVAR?

ADD COMMENTlink written 3.1 years ago by pfs280

I'm looking into it now. I found instructions on their website for adding gene definition databases for non-human species, but they need to be in the USCS Genome Browser. I'm studying poplar, which is not... but I found in literature ANNOVAR was somehow used for rice, which isn't either.

Do you know how to add databases for species not included in the USCS Genome Browser, or could you point me in the right direction? Thanks.

ADD REPLYlink written 3.1 years ago by michael.nagle100

I do not have first hand experience doing this but the 'Experimental Design' section in 'Genomic variant annotation and prioritization with ANNOVAR and wANNOVAR' should get you started. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4718734/

ADD REPLYlink written 3.1 years ago by pfs280
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1308 users visited in the last hour
_