Qtlizer: comprehensive QTL annotation of GWAS results
This R package provides access to the Qtlizer web server. Qtlizer annotates lists of common small variants (mainly SNPs) and genes in humans with associated changes in gene expression using the most comprehensive database of published quantitative trait loci (QTLs).
Alternatively, the Qtlizer can be accessed by using a web-based GUI (http://genehopper.de/qtlizer). More information about usage and available datasets can be found at http://genehopper.de/help#qtlizer_docu.
Please note: A valid internet connection (HTTP port: 80) is required in order to install and use the package.
Simply call the function
get_qtls() function to make requests to Qtilzer. The function utilizes a REST API (http://genehopper.de/rest) to query the annotation database. The QTL results will be returned as data frame or as
Common seperators (space, comma, space + comma, ...) are accepted. It is also possible to pass your query with a vector:
Accepted query terms
Accepted query terms are variant and gene identifiers of the form:
- Rsid : rs + number e.g. "rs4284742"
- reference:chr:pos e.g. "hg19:19:45412079" (Allowed references: hg19/GRCh37, hg38/GRCh38; accepted chromosomes are 1-22)
- Gene symbol consisting of letters and numbers according to https://www.genenames.org/about/guidelines/
corr: Correlation threshold based on linkage disequilibrium (LD) calculated from the 1000 Genomes Phase 3 European dataset. Optional value between 0 and 1. Default value is NA.
ld_method: Method to calculate correlation. Valid values are either "r2" (default) or "dprime".
get_qtls("rs4284742", corr = 0.6, ld_method="r2")
return_obj: Determinse the format of the result. Value "dataframe" (default) returns a dara frame whereas "grange" returns a
ref_version: If output is a GRange object, the version of the reference genome is also considered. Accepted reference genome versions are "hg19" (default) or "hg38".
get_qtls("rs4284742", return_obj = "grange", ref_version = "hg38")
max_terms: Number of queries made at a time. The default value is 5. It is recommended to not set the value higher than 5.
get_qtls("rs4284742", max_terms = 4)
Column descriptions of the received data frame can be accessed by calling:
df = get_qtls("rs4284742") comment(df)
Try out online
If you want to try out the R package online, there is an example Google Colaboratory project at
To run the project, make a private copy or open the project in playground mode and sign in to Google.
Please cite the following article when using Qtlizer:
Munz M, Wohlers I, Simon E, Busch H, Schaefer A* and Erdmann J * (2018) Qtlizer: comprehensive QTL annotation of GWAS results. bioRxiv
GNU General Public License v3.0