Demonstration of TASSEL - Cornell Universitycbsu.tc.cornell.edu/lab/doc/Lipka_A_GBS_Tutorial_20130916.pdfOverview • Download TASSEL • Read in data • GWAS using GLM • Calculate

Demonstration of TASSEL

Alexander E. LipkaPostdoctoral Researcher

Buckler Lab

1

Overview

• Download TASSEL• Read in data• GWAS using GLM• Calculate kinship matrix and conduct PCA• GWAS using MLM• Visualize LD• Cladogram 

2

Download TASSEL

3

www.maizegenetics.net/tassel

Reading in data

4

Data > Load > I will make my best guess and try

Fit a generalized linear model (GLM)

5

Select “mdp_traits” and “data_hapmap”, and Data > Intersect Join

Fit a generalized linear model (GLM)

6

Select “mdp_traits+data_hapmap”, then Analysis > GLM > OK

Make a Manhattan and QQ‐Plot

7

Select GLM model results, then Results > Manhattan Plot/QQ Plot

Calculate a kinship matrix

8

Select “data_hapmap”, then Analysis > Kinship

Run PCA

9

Select “data_hapmap”, then Data > Transform > Create Dataset

Run PCA

10

Select “data_hapmap_Collapsed”, then Data > Transform > Impute> Create Dataset

Run PCA

11

Select “data_hapmap_Collapsed_3093_imputed”, then Data > Transform > PCA > Create Dataset

Run MLM

12

Select “mdp_population_structure”, then Filter> Traits > Deselect Q3 > OK

Run MLM

13

Intersect Join “mdp_traits”, “data_hapmap”, and “Filtered_mdp_population_structure”Select merged file and “kin_data_hapmap”, then Analysis > MLM > Run

QQ and Manhattan Plot from MLM

14

Investigate LD on Chromosome 1

15

Select “data_hapmap”, then Filter > Sites > Select Chromosomes

Make an LD  Plot

16

Select the Chr1 filtered data, then Analysis > Link. Diseq. > Run

Make an LD Plot

17

Select the LD data, then Results > LD Plot

Make a Cladogram

18

Select the Chr1 filtered data, then Analysis > Cladogram > Run

19

Make a CladogramSelect the Tree data, then Results > Archaeopteryx Tree

Questions?

• Requests for other analyses in TASSEL?

20