之前做SSR标记和简化测序都是用的桌面化的TASSEL进行运算，对于做关联分析的初学者非常简单，十分好入门。但是关于命令行的操作就要麻烦一些，但是运算起来不像GAPIT那样消耗大量的内存，下面做一些简单的示例。

1.安装

1.1 Java JDK 8.0的安装

java下载

下载地址

#下载

wget http://download.oracle.com/otn-pub/java/jdk/8u112-b15/jdk-8u112-linux-x64.rpm

#安装

rpm -ivhc jdk-8u112-linux-x64.rpm

#测试安装

java -vsrsion

#查看java版本

[root@lilibei java]# java -version

java version "1.8.0_112"

Java(TM) SE Runtime Environment (build 1.8.0_112-b15)

Java HotSpot(TM) 64-Bit Server VM (build 25.112-b15, mixed mode)

更改环境变量

vim + /etc/profile

按照如下添加

JAVA_HOME=/usr/java/jdk1.8.0

JRE_HOME=/usr/java/jdk1.8.0/jre

PATH=$PATH:$JAVA_HOME/bin:$JRE_HOME/bin

CLASSPATH=.:$JAVA_HOME/lib/dt.jar:$JAVA_HOME/lib/tools.jar:$JRE_HOME/lib

export JAVA_HOME JRE_HOME PATH CLASSPATH

#编辑完后执行脚本

source /etc/profile

[root@cricaas java]# echo $PATH

/usr/lib64/qt-3.3/bin:/usr/local/sbin:/usr/local/bin:/sbin:/bin:/usr/sbin:/usr/bin:/usr/local/git/bin:/home/lilibei/bin:/root/bin:/usr/local/git/bin:/home/lilibei/bin:/usr/java/jdk1.8.0/bin:/usr/java/jdk1.8.0/jre/bin

[root@cricaas java]# echo $JAVA_HOME

/usr/java/jdk1.8.0

[root@cricaas java]#

1.2 TASSEL的安装

TASSEL的安装很简单，按照如下即可

#github clone

git clone https://bitbucket.org/tasseladmin/tassel-5-standalone.git

cd tassel-5-standalone

2.使用

示例文件均存放于

/home/lilibei/tassel-5-standalone/TASSELTutorialData/data/

2.1 计算kinship矩阵

#适合小数据集，重测序数据需要更改内存配置参数

./run_pipeline.pl -fork1 -h TASSELTutorialData/data/mdp_genotype.hmp.txt -ck -export kinship.txt

2.2 LD分析

#绘制LD连锁图，很难看，可以使用hapview或者LDheatmap

./run_pipeline.pl -fork1 -h TASSELTutorialData/data/mdp_genotype.hmp.txt -ld -ldd png -o chr_5000sites_ld.png

#计算LD

./run_pipeline.pl -fork1 -h TASSELTutorialData/data/mdp_genotype.hmp.txt -ld -export LD_5000sites

2.3 提取染色体

#提取每条染色体

./run_pipeline.pl -fork1 -h TASSELTutorialData/data/mdp_genotype.hmp.txt -separate -export chromosome

#提取1和3号染色体

./run_pipeline.pl -fork1 -h TASSELTutorialData/data/mdp_genotype.hmp.txt -separate 1,3 -export chromosome

2.4 MLM分析

./run_pipeline.pl -fork1 -h TASSELTutorialData/data/mdp_genotype.hmp.txt -filterAlign -filterAlignMinFreq 0.05 -fork2 -r TASSELTutorialData/data/mdp_traits.txt -fork3 -q TASSELTutorialData/data/mdp_population_structure.txt -excludeLastTrait -fork4 -k TASSELTutorialData/data/mdp_kinship.txt -combine5 -input1 -input2 -input3 -intersect -combine6 -input5 -input4 -mlm -export lilibei -runfork1 -runfork2 -runfork3 -runfork4

3.参考资料：

有疑问加站长微信联系(非本文作者)