QTLtools-gwas(1)

Okay bestie, this is the page for QTLtools-gwas(1) from the bioinformatics tools crew 🧚 and it's giving main character energy ✨

✨ NAME ✨

QTLtools gwas is for genome-wide association tests πŸ’–

πŸ“œ SYNOPSIS

QTLtools gwas --vcf [in.vcf|in.vcf.gz|in.bcf|in.bed.gz] --bed
quantifications.bed.gz --out output.txt [OPTIONS]

🌸 DESCRIPTION

Okay so this mode runs a linear regression with every genotype vs. every phenotype. Every single pair, no skipping πŸ’…

πŸŽ€ OPTIONS

  • --vcf [in.vcf|in.bcf|in.vcf.gz|in.bed.gz] ... β€” Your genotypes go in VCF/BCF/BED format πŸ¦‹ The DS field in the genotype FORMAT of a variant (the dosage calculated from genotype probabilities, e.g. after imputation) is REQUIRED. REQUIRED.
  • --bed quantifications.bed.gz ... β€” Molecular phenotype quantifications in BED format. REQUIRED πŸ“
  • --out filename.txt β€” The output file name. REQUIRED πŸ’œ
  • --cov covariates.txt β€” Covariates you want to correct the phenotype data with 🫢
  • --normal β€” Rank normal transform the phenotype data so each phenotype comes out normally distributed. RECOMMENDED 🌈

πŸ“Š OUTPUT FILE

Output file vibes: it's a space separated file with these columns πŸŽ€

  • 1 β€” the phenotype ID
  • 2 β€” the variant's chromosome
  • 3 β€” the variant's start position
  • 3.1 β€” the variant's end position. Only printed if a BED file is supplied to --vcf πŸ‘€
  • 4 β€” the variant's ID
  • 5 β€” the p-value
  • 6 β€” the correlation coefficient

🧚 EXAMPLE

  • Run the association between all the variants and phenotypes, correcting the phenotypes with technical covariates and normal transforming them πŸ”₯
QTLtools gwas --bed quantifications.bed.gz --vcf genotypes.bcf --cov covariates.txt --normal --out gwas_results.txt

πŸ”— SEE ALSO

πŸ™ˆ BUGS

Versions up to and including 1.2 suffer from a bug in reading missing genotypes in VCF/BCF files, bestie. 😭 This bug affects variants with a DS field in their genotype's FORMAT and have a missing genotype (DS field is .) in one of the samples. When that happens, genotypes for all the samples are set to missing, which effectively removes this variant from the analyses πŸ’€

Please submit bugs to https://github.com/qtltools/qtltools πŸ›

πŸ“š CITATION

Delaneau O., Ongen H., Brown A. A., et al. A complete tool set for molecular QTL discovery and analysis. Nat Commun 8, 15452 (2017). https://doi.org/10.1038/ncomms15452 πŸ’œ

πŸ‘‘ AUTHORS

Olivier Delaneau (olivier.delaneau@gmail.com), Halit Ongen (halitongen@gmail.com) 🫢


Signing off with QTLtools-v1.3, last updated 06 May 2020, for QTLtools-gwas(1) 🌈✨