FastQTL(1)

πŸ’– FastQTL πŸ’–

okay so bestie, welcome to the FastQTL manual page β€” that's the General Commands Manual section 1 vibe, iykyk. she's a QTL mapper, she's iconic, she's about to change your whole genomics career. let's get into it 🧚✨


✨ NAME ✨

FastQTL β€” Quantitative Trait Loci (QTL) mapper in cis for molecular phenotypes πŸ§¬πŸ’…


🌸 SYNOPSIS 🌸

fastQTL [options]

πŸ¦‹ DESCRIPTION πŸ¦‹

okay so here's the tea πŸ“ β€” the whole goal of FastQTL is to identify single-nucleotide polymorphisms (SNPs) that are significantly associated with various molecular phenotypes (like, expression of known genes, cytosine methylation levels, etc. β€” she's versatile, she contains multitudes πŸ’œ).

She performs scans for all possible phenotype-variant pairs in cis β€” meaning variants located within a specific window around a phenotype. No cap, she does all of them.

And then β€” the main character moment β€” FastQTL implements a new permutation scheme (Beta approximation) to accurately and rapidly correct for multiple-testing at both the genotype and phenotype levels. she ate that πŸ”₯


πŸ‘‘ OPTIONS πŸ‘‘

πŸ’« Basic options

  • --help β€” prints help about options. the "bestie what do I do" button 🫢
  • --silent β€” silent mode on terminal. she's not talking rn πŸ™ˆ
  • --seed arg (=1459232241) β€” random number generator seed. super useful when you want to replicate runs of the software β€” same seed, same results, periodt. the default is literally 1459232241.

πŸ“ Input and Output files options

  • -L, --log β€” screen output gets copied into this file. receipts, basically πŸ’…
  • -V, --vcf β€” genotypes in VCF format.
  • -B, --bed β€” phenotypes in BED format.
  • -C, --cov β€” covariates in TXT format.
  • -G, --grp β€” phenotype groups in TXT format.
  • -O, --out β€” output file. the moment of truth ✨

🚫 Exclusion and Inclusion files options

  • --exclude-samples β€” list of samples to exclude. bye girl πŸ‘‹
  • --include-samples β€” list of samples to include. come thru πŸŽ€
  • --exclude-sites β€” list of sites to exclude.
  • --include-sites β€” list of sites to include.
  • --exclude-phenotypes β€” list of phenotypes to exclude.
  • --include-phenotypes β€” list of phenotypes to include.
  • --exclude-covariates β€” list of covariates to exclude.
  • --include-covariates β€” list of covariates to include.

βš™οΈ Parameters

  • --normal β€” performs quantile normalization on the phenotype quantifications to make them normally distributed. implemented as the rntransform function of the GenABEL package. she's getting her life together 🌈
  • -W, --window β€” cis-window size. the default value is 1Mb (1e6). what that means: all variants within 1e6 bp of the phenotype location (e.g. TSS) get analyzed. so like, a whole neighborhood gets a vibe check.
  • -T, --threshold β€” filters out all phenotype-variant pairs with a p-value above the specified threshold in the output of a nominal pass. no mediocre p-values allowed πŸ’€

🎭 Modes

  • -P, --permute β€” permutation pass to calculate corrected p-values for molecular phenotypes.
  • --psequence β€” permutation sequence.
  • --map β€” maps best QTL candidates per molecular phenotype. finding the main character of each phenotype πŸ‘‘
  • --map-full β€” scans the full cis-window to discover independent signals. she's thorough ✨
  • --interaction β€” tests for interactions with the variable specified in file.

πŸƒβ™€οΈ Parallelization

  • -K, --chunk β€” specify which chunk needs to be processed. divide and conquer, bestie.
  • --commands β€” generates all commands. she's writing the to-do list for you 🫢
  • -R, --region β€” region of interest.

🌟 EXAMPLE 🌟

here's the iconic loop, in all her glory:

for c in $(seq 1 256); do
     fastQTL --vcf genotypes.vcf.gz
             --bed phenotypes.bed.gz
             --chunk $c 256
             --permute 1000
             --out results.$c.txt.gz
done
zcat results.*.txt.gz | gzip -c > fastqtl_expected_output.txt.gz
note: this is copied exactly as she appears in the original page β€” the line breaks inside the fastQTL invocation are how they're written there, so don't come for me πŸ™ˆ

πŸ”— SEE ALSO πŸ”—


πŸ’œ AUTHOR πŸ’œ

FastQTL was written by Olivier Delaneau, Halit Ongen, Alfonso Buil and Manolis Dermitzakis. absolute legends, all four of them πŸ”₯πŸ‘‘

This manual page was written by Dylan AΓ―ssi <bob.dybian@gmail.com>, for the Debian project (but may be used by others). shoutout to Dylan for the lore drop πŸŽ€


okay and that's the whole page bestie!! this is version 2.184 from March 2016 β€” yeah she's a vintage icon and she still eats πŸ’…βœ¨ go run some cis-QTL scans, ily, bye!! πŸŒΈπŸ¦‹πŸ’–