Portable eQTL analysis and statistical fine mapping workflow used by the eQTL Catalogue
eQTL-Catalogue/qtlmap is a bioinformatics analysis pipeline used for QTL Analysis.
The workflow takes phenotype count matrix (normalized and quality controlled) and genotype data as input, and finds associations between them with the help of sample metadata and phenotype metadata files (See Input formats and preparation for required input file details). To map QTLs, pipeline uses QTLTools's PCA and RUN methods. For manipulation of files BcfTools, Tabix and custom Rscript scripts are used.
The pipeline is built using Nextflow, a bioinformatics workflow tool to run tasks across multiple compute infrastructures in a very portable manner. It comes with docker / singularity containers making installation trivial and results highly reproducible.
The eQTL-Catalogue/qtlmap pipeline comes with documentation about the pipeline, found in the docs/ directory:
- Installation
- Pipeline configuration
- Input formats and preparation
- Running the pipeline
- Troubleshooting
Mapping QTLs is a process of finding statistically significant associations between phenotypes and genetic variants located nearby (within a specific window around phenotype, a.k.a cis window) This pipeline is designed to perform QTL mapping. It is intended to add this pipeline to the nf-core framework in the future. High level representation of the pipeline is shown below:
The output directory of the workflow contains the following subdirectories:
- PCA - genotype and gene expression PCA values used as covariates for QTL analysis.
- sumstats - QTL summary statistics from nominal and permutation passes.
- susie - SuSiE fine mapping credible sets.
- susie_full - full set of susie results for all tested variants (very large files).
- susie_merged - susie credible sets merged with summary statistics from univariate QTL analysis.
Column names of the output files are explained here.
- Nurlan Kerimov
- Kaur Alasoo
- Masahiro Kanai
- Ralf Tambets