The skin microbiome of elasmobranchs follows phylosymbosis, but in teleost fishes the microbiomes converge
Authors: Doane Michael , Morris Megan M. , Papudeshi Bhavya , Allen Lauren , Pande Dnyanada , Haggerty John M. , Johri Shaili , Turnland Abby , Peterson Meredith , Kacev Dovi , Nosal Andy , Ramirez Deni , Hovel Kevin , Ledbetter Julia , Alker Amanda , Avalos Jackeline , Baker Kristi , Bhide Shruti , Billings Emma , Byrum Steven , Clemens Molly , Demery Amelia Juliette , Lima Lais Farias Oliveira , Gomez Oscar , Gutierrez Omar , Hinton Selena , Kieu Donald , Kim Angie , Loaiza Rebeca , Martinez Alexander , McGhee Jordan , Nguyen Kristine , Parlan Sabrina , Pham Amanda , Price-Waldman Rosalyn , Edwards Robert A. , Dinsdale Elizabeth A. 
Published: December 2019

Abstract
Background: The vertebrate clade diverged into Chondrichthyes (sharks, rays, and chimeras) and Osteichthyes fishes (bony fishes) approximately 420 mya, with each group accumulating vast anatomical and physiological differences, including skin properties. The skin of Chondrichthyes fishes is covered in dermal denticles, whereas Osteichthyes fishes are covered in scales and are mucous rich. The divergence time among these two fish groups is hypothesized to result in predictable variation among symbionts. Here, using shotgun metagenomics, we test if patterns of diversity in the skin surface microbiome across the two fish clades match predictions made by phylosymbiosis theory. We hypothesize, 1) the skin microbiome will be species and clade-specific, 2) evolutionary difference in elasmobranch and teleost fishes corresponds with a concomitant increase in host-microbiome dissimilarity and 3) the skin structure of the two groups will affect the taxonomic and functional composition of the microbiomes. Results: We show that the taxonomic and functional composition of the microbiomes is host species-specific. Teleost fish had lower average microbiome within clade similarity compared to among clade comparison, but their composition is not different among clade in a null based model. Elasmobranch’s average similarity within clade was not different than across clade and not different in a null based model of comparison. In the comparison of host distance with microbiome distance, we found that the taxonomic composition of the microbiome was related to host distance for the elasmobranchs, but not the teleost fishes. In comparison, the gene function composition was not related to the host-organism distance for elasmobranchs but was negatively correlated with host distance for teleost fishes. Conclusion: Our results show the patterns of phylosymbiosis are not consistent across both fish clades, with the elasmobranchs showing phylosymbiosis, while the teleost fish are not. The discrepancy may be linked to alternative processes underpinning microbiome assemblage, including possible historical host-microbiome evolution of the elasmobranchs and convergent evolution in the teleost which filter specific microbial groups. Our comparison of the microbiomes among fishes represents an investigation into the microbial relationships of the oldest divergence extant vertebrate host and reveals that microbial relationships are not consistent across evolutionary timescales.