RESUMEN
BACKGROUND: Fish skin microbiomes are rarely studied in inland water systems, in spite of their importance for fish health and ecology. This is mainly because fish species distribution often covaries with other biotic and abiotic factors, complicating the study design. We tackled this issue in the northern part of the Jordan River system, in which a few fish species geographically overlap, across steep gradients of water temperature and salinity. RESULTS: Using 16S rRNA metabarcoding, we studied the water properties that shape the skin bacterial communities, and their interaction with fish taxonomy. To better characterise the indigenous skin community, we excluded bacteria that were equally abundant in the skin samples and in the water samples, from our analysis of the skin samples. With this in mind, we found alpha diversity of the skin communities to be stable across sites, but higher in benthic loaches, compared to other fish. Beta diversity was found to be different among sites and to weakly covary with the dissolved oxygen, when treated skin communities were considered. In contrast, water temperature and conductivity were strong factors explaining beta diversity in the untreated skin communities. Beta diversity differences between co-occurring fish species emerged only for the treated skin communities. Metagenomics predictions highlighted the microbiome functional implications of excluding the water community contamination from the fish skin communities. Finally, we found that human-induced eutrophication promotes dysbiosis of the fish skin community, with signatures relating to fish health. CONCLUSIONS: Consideration of the background water microbiome when studying fish skin microbiomes, across varying fish species and water properties, exposes patterns otherwise undetected and highlight among-fish-species differences. We suggest that sporadic nutrient pollution events, otherwise undetected, drive fish skin communities to dysbiosis. This finding is in line with a recent study, showing that biofilms capture sporadic pollution events, undetectable by interspersed water monitoring. Video abstract.