Map SNPs (both coding and non-coding) onto genes using TxDb.Hsapiens.UCSC.hg19.knownGene or TxDb.Hsapiens.UCSC.hg38.knownGene as a reference.
map_snps2genes_txdb( sumstats, promoter_upstream = 35000, promoter_downstream = 10000, queries = c("promoter", "coding", "intron", "threeUTR", "splicesite"), nCores = 1, verbose = TRUE )
| sumstats | GWAS summary statistics munged by format_sumstats. |
|---|---|
| promoter_upstream | How many basepairs upstream of known promoters to search for genes. |
| promoter_downstream | How many basepairs dowmstream of known promoters to search for genes. |
| queries | Genomic region types to search for overlap with |
| nCores | Number of cores to parellelise across. |
| verbose | Print messages. |
| sumstats_file | GWAS summary statistics munged by format_sumstats. Can be a path to the saved file or data.table. |
| method | Method to use for mapping SNPs to genes. |
gene_hits data.table