2.8 Data analysis
All pairwise comparisons were made using the Wilcoxon Rank-sum test of the stats package (R Core Team, 2020). Pearson correlation tests were added to plots with the ggpupr package (Kassambara, 2020). Permutational multivariate analysis of variance (PERMANOVA) tests was performed using the vegan package ({Jari Oksanen et al., 2020), with a default 999 permutations.