Fig. 3

a Alpha diversity measure: Species Richness with the corresponding colour legend for the treatment groups. Pair-wise ANOVA P-values are displayed P < 0.05*, P < 0.01**, P < 0.001***. b Non-metric distance scaling (NMDS) plot based on Unifrac phylogenetic distance with colour coded circles representing the differing treatment groups over time (R2 = 0.27182, P = 0.001)