使用FUSION进行TWAS分析

欢迎关注”生信修炼手册”!

FUSION是一款进行TWAS分析的软件，对应的文章发表在nature genetic上，链接如下

https://www.nature.com/articles/ng.3506

软件的官网如下

http://gusevlab.org/projects/fusion/

该软件同时支持对单个样本和多个样本的gwas cohort中的样本进行基因型填充，对应的流程图如下

首先对reference panel中进行建模，对于每个基因，筛选对应的cis-snp, 即在该基因上下游1Mb范围内的snp位点，通过不同的模型构建cis-SNPs和对应的gene之间的关系，接下来对gwas cohort中的样本进行填充，如上图中A所示，对于单个样本，根据其cis-SNPs的基因分析结果，预测对应基因的表达量，然后与表型性状进行关联分析；如上图中B所示，对于大量样本的gwas summary数据，考虑了SNP的连锁情况，根据cis-SNP和gene的关系，直接预测基因和性状之间的关联效应。

该软件采用R语言进行开发，基本用法如下

1.  利用reference panel建模

通过不同的模型，计算cis-SNPs和基因的关系，代码如下

Rscript FUSION.compute_weights.R \
--bfile $INP \
--tmp $TMP \
--out $OUT \
--models top1,blup,bslmm,lasso,enet

bfile参数指定plink格式的基因分型结果，tmp参数指定临时文件存放的目录，out参数指定输出文件的目录，运行成功后，在输出结果的目录，每个基因都会生成1个后缀为RDat的文件。接下来需要生成一个profile文件，首先准备一个list文件，每一行为一个基因对应的RData文件的路径，示意如下

运行以下命令生成profile

Rscript utils/FUSION.profile_wgt.R  GTEx.Whole_Blood.list

对于GTEx, TCGA等公共数据，官网提供了