Famrvtest

From Genome Analysis Wiki
Revision as of 01:17, 14 July 2014 by Shuang Feng (talk | contribs) (How to Execute)
Jump to: navigation, search

Useful Wiki Pages

There are a few pages in this Wiki that may be useful to famRvTest users. Here are links to key pages:

Brief Description

famrvtest is a computationally efficient tool for family-based association analyses of rare variants using sequencing or genotyping array data. famrvtest supports both single variant and gene-level associations.

For any questions, please contact Shuang Feng (sfengsph at umich.edu) or Gonçalo Abecasis (goncalo at umich.edu).

Download and Installation

  • University of Michigan CSG users can go to the following:
 /net/fantasia/home/sfengsph/code/famrvtest/bin/famrvtest

Where to Download

  • Source code can be downloaded in the following
  Source for LINUX
  Source for MAC
  Source for MINGW
  Source for CYGWIN64
  • Executable can be downloaded in the following:
  Executable for LINUX

How to Compile

  • Save it to your local path and decompress using the following command:
 tar xvzf Linux_famRvTest.2.0.tgz
  • Go to promp>famrvtest and type the following command to compile:
 make -j 5

How to Execute

  • Go to famrvtest/bin and use the following:
 ./famrvtest

Command Reference

Please go to Command Reference Page for details.

Approach

famRvTest uses linear mixed model approach, incorporating efficient optimization algorithm, to account for familial relationship, where kinship is either quantified based upon pedigree structures or estimated from genotypes of markers from genome-wide. Single marker associations including score, likelihood ratio and ward tests and gene-level associations methods (weighted and un-weighted burden, SKAT and variable threshold tests) have been implemented. Manuscript is under preparation.

Input Files

famRvTest needs the following files as input: PED and DAT file in Merlin format, AND/OR a VCF file. When genotypes are stored in PED and DAT file, the VCF file is not needed. However, even if genotypes are saved in a VCF file, PED and DAT files are still needed for carrying covariate and trait information.

PED and DAT Files

  • When PED file has genotypes saved, there is no need for a VCF file as input.
  • famRvTest takes PED/DAT file in Merlin format. Please refer to PED/DAT format description for details.
  • An example PED file is in the following:
    1 1 0 0 1 1.5 1 23 A A A A A A A A A A
    2 1 0 0 1 1.0 1 34 A C A C A C A C A C
    3 1 0 0 2 0.4 1 43 A A A A A A A A A A
    4 1 0 0 2 0.9 1 13 A C A C A C A C A C
  • The matching DAT file is in the following:
 T YourTraitName
 C SEX
 C AGE
 M 1:123456a
 M 1:234567
 M 2:111111
 M 2:222222
 M X:12345
  • DAT file must have variant names in the following format "M chr:pos".
  • Orders of labels in DAT file have to match the order of fields in PED file.
  • Markers in PED and DAT file must be sorted by chromosome and position.
  • Covariate and trait values are saved in PED file. Covariate and trait descriptions are saved in DAT file.

VCF File

  • Another option is to use VCF as input. Please refer to the following link for VCF file specification: 1000 genome wiki VCF specs
  • VCF file should be compressed by bgzip and indexed by tabix, using the following command:
 bgzip input.vcf     ## this command will generate input.vcf.gz
 tabix -p vcf -f input.vcf.gz  ## this command will generate input.vcf.gz.tbi
  • Even with the presence of VCF file, PED/DAT files are still needed for covariates and phenotypes.

Example Command Line

Single Variant Analysis

The following command lines let you run single variant association analysis of trait "LDL" using score test, after inverse normalization of the quantitative trait and adjusting covariates. --traitName specifies the single trait or traits you want to analyze in this batch. If this option is not used, then all traits coded in data file will be analyzed accordingly. --SingleVarLRT provides essentially the same test as in merlin --fastAssoc option.

./famRvTest -p your.ped -d your.dat --vcf your.vcf.gz --SingleVarScore --inverseNormal --useCovariates --traitName LDL

Futhermore, if you want to run likelihood ratio test and wald test at the same time, the following command should do the work:

./famRvTest -p your.ped -d your.dat --SingleVarScore --SingleVarLRT --SingleVarWald --inverseNormal --useCovariates --traitName LDL

All the above commands will let you do family-based association analysis using kinship matrices generated using pedigree structure coded in pedigree file. The following command lines show examples of using genotype to estimate empirical relationship matrix to do the work.

 ./famRvTest -p your.ped -d your.dat --SingleVarScore --SingleVarLRT --SingleVarWald --inverseNormal --useCovariates --traitName LDL --kinGeno

Gene-level Association

The following command lines let you run gene-level association analysis of genes listed in "your.genes.groupfile" for trait "LDL" using SKAT, Madsen-Browning weighted burden, rare allele counts un-weighted burden and collapsing burden and variable threshold tests, after inverse normalization of the quantitative trait and adjusting covariates. Only rare variants with maf less than or equal to 0.05 and minor allele count greater than or equal to 3 are grouped.

./famRvTest -p your.ped -d your.dat --SKAT --MB --CMC_counts --CMC_binary --VTasymptotic --inverseNormal --useCovariates --traitName LDL --groupFile your.genes.groupfile --maxMaf 0.05 --mac 3

Change Log

  • Released version 0.0.9 with a bug fixed for potential compiling error. (10/10/2013)
  • Released version 2.0, a faster version and added single variant permutation test. (7/14/2014)