Expression of genes involved in nodule development in Datisca glomerata

Differences in gene expression were assessed by analysis of variance (ANOVA). Equality of variances was checked by Levene’s test.

For equal variances, one-way ANOVA was fitted to log2-transformed 2^(-ΔΔCq) with treatments as fixed term and biological repeats as random effects accounting for dependencies between treatments. To estimate and adjust for pairwise multiple comparisons, a non-parametric post-hoc Tukey HSD test was applied. For unequal variances, a Welch’s pairwise t-test with Holm’s correction was used.