Abstract
Amphibians are known to have an abundance of microorganisms colonizing their skin, and these symbionts often protect the host from disease. There are now many comprehensive studies on amphibian skin microbes, but the interspecific and intraspecific abundance distributions (or abundance heterogeneity) of amphibian skin microbes remain unclear. Furthermore, we have a very limited understanding of how the abundance and heterogeneity of microbial communities relate to the body size (or more specifically, skin surface area) of amphibian hosts. In this study, we evaluated the interspecific and intraspecific abundance distribution patterns of amphibian skin microbes and evaluated whether the symbiotic skin microbes of different anuran species share a fundamental heterogeneity scaling parameter. If scaling invariance exists, we hypothesize that a fundamental heterogeneity scaling value also exists. A total of 358 specimens of 10 amphibian host species were collected, and we used Type-I and III Taylor's power law expansions (TPLE) to assess amphibian skin microbial heterogeneity at the community and mixed-species population levels, respectively. The obtained results showed that, at the community scale, a high aggregation of the microbial abundance distribution on the skin barely changed with host size. In a mixed-species population (i.e., a community context), the abundance distribution pattern of mixed microbial species populations also does not change with host size and always remains highly aggregated. These findings suggest that while amphibian skin microbiomes located in different hosts may have different environmental conditions, they share a fundamental heterogeneity scaling parameter, and thus, scale invariance exists. Finally, we found that microhabitat area provided by the host skin is vital to the stability of the symbiotic microbial community.