diff --git a/assays/AnalysisOfDifferencesInExpression/README.md b/assays/AnalysisOfDifferencesInExpression/README.md
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391
diff --git a/assays/AnalysisOfDifferencesInExpression/dataset/.gitkeep b/assays/AnalysisOfDifferencesInExpression/dataset/.gitkeep
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391
diff --git a/assays/AnalysisOfDifferencesInExpression/isa.assay.xlsx b/assays/AnalysisOfDifferencesInExpression/isa.assay.xlsx
new file mode 100644
index 0000000000000000000000000000000000000000..b8cce38e068281027c92950637a2f61387afac74
Binary files /dev/null and b/assays/AnalysisOfDifferencesInExpression/isa.assay.xlsx differ
diff --git a/assays/AnalysisOfDifferencesInExpression/protocols/.gitkeep b/assays/AnalysisOfDifferencesInExpression/protocols/.gitkeep
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391
diff --git a/assays/AnalysisOfDifferencesInExpression/protocols/AnalysingDifferencesInExpression.md b/assays/AnalysisOfDifferencesInExpression/protocols/AnalysingDifferencesInExpression.md
new file mode 100644
index 0000000000000000000000000000000000000000..68fad9bb0141042b5c38f71b170a5bf3ef91156c
--- /dev/null
+++ b/assays/AnalysisOfDifferencesInExpression/protocols/AnalysingDifferencesInExpression.md
@@ -0,0 +1,5 @@
+## Analysis of differences in expression
+
+To identify the factors associated with differences in expression of *R*-genes across transcriptomes, we performed an ANOSIM in Primer 7.0.13 (PRIMER-e; Figure S2). ANOSIM is a non-parametric statistical test similar to ANOVA. The starting point of the analysis is a pairwise dissimilarity matrix. In our case, the dissimilarity matrix was computed as follows: First the TPM values for each gene within each transcriptome were LOG (x+1) transformed. On the basis of these transformed TPM values, the dissimilarity in gene expression patterns between transcriptomes were calculated based on Euclidean distances. Ranking was applied to the distance matrix. The two libraries from potato (SRR6511453 and ERR791944) with exceptionally low expression of the entire *R*-gene repertoire were excluded in these analyses.
+
+To determine if gene expression is more similar within groups than between groups (for example when groups are defined by infection status or tissue type) the R test statistic value was calculated. The R values can range from -1 to 1, with larger values corresponding to greater differences between groups. Statistical significance is calculated through permutation of the group labels and recalculation of the R value for each replicate. In our case, 999 permutations were generated. The following factors were evaluated: BioProject, tissue type, type of treatment, specific treatment organism, life cycle of the organism, type/kingdom of the organism, susceptible vs. resistant cultivar, relative sequencing depth, paired- or single-end sequencing and days post infection. The ANOSIM analysis was also applied to the differential gene expression data (see below).
\ No newline at end of file