# Qtlizer: comprehensive QTL annotation of GWAS results
[![Twitter](https://img.shields.io/twitter/url/http/shields.io.svg?style=social)](https://twitter.com/intent/tweet?hashtags=Qtlizer&url=https://www.biorxiv.org/content/10.1101/495903v2)


 + [Introduction](#Introduction)\
 + [Installation](#Installation)\
 + [Usage](#Usage)\
    |-- [Accepted query terms](#Accepted-query-terms)\
    |-- [Optional parameters](#Optional-parameters)\
    |-- [Meta information](#Meta-information)\
 + [Try out online](#Try-out-online)\
 + [Authors](#Authors)\
 + [Citation](#Citation)\
 + [License](#License)


## Introduction
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.


## Installation
```R
devtools::install_github('matmu/Qtlizer', build_vignettes = TRUE)
```

**Please note**: A valid internet connection (HTTP port: 80) is required in order to install and use the package.


## Help pages
```R
browseVignettes(Qtlizer)
help(package="Qtlizer")
```


## Usage
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 `GenomicRanges::GRanges` object.

```R
get_qtls('rs4284742 DEFA1')
```
Common seperators (space, comma, space + comma, ...) are accepted. It is also possible to pass your query with a vector: 

```R
get_qtls(c("rs4284742", "DEFA1"))
```


### 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/]( https://www.genenames.org/about/guidelines/)


### Optional parameters
+ 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".

	```R
	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 ```GenomicRanges::GRanges``` object. 
+ 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".

	```R
	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. 

	```R
	get_qtls("rs4284742", max_terms = 4)
	```

### Meta information
Column descriptions of the received data frame can be accessed by calling:

```R
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

https://colab.research.google.com/drive/1i1sjQHCjaw2wYzVBnXQ9iaghnk-jSU95#scrollTo=5Hi6sCe7SPFb

To run the project, make a private copy or open the project in playground mode and sign in to Google. 

## Authors
Matthias Munz [![](https://img.shields.io/twitter/follow/_MatthiasMunz?label=Follow&style=social)](https://img.shields.io/twitter/follow/_MatthiasMunz?label=Follow&style=social)\
Julia Remes\
University of Lübeck, Germany

## Citation
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***

[![](https://img.shields.io/badge/doi-https%3A%2F%2Fdoi.org%2F10.1101%2F495903%20-green.svg)](https://doi.org/10.1101/495903)
[![](https://img.shields.io/badge/Altmetric-17-green.svg)](https://www.altmetric.com/details/52777590)

## License
GNU General Public License v3.0

