Fig. 4

Heatmaps of differential abundance among host and water communities. This figure from left to right includes 9 heatmaps of the significant taxonomic fingerprints at the genus level between gut, the skin and the water at times T0 (first column), T3 (second column), and T5 (third column) in the control (first row), the CV (second row), and the CC (third row) groups. The hierarchical clustering of the relative abundance of phyla which significantly changed over time was performed using Ward’s method and Bray–Curtis dissimilarity distance. Vegan package and pheatmap () function in R were used