“…To select best-fit models, we added to Models 1-4 the non-bioclimatic, non-categorical predictors: 'elevation' , 'latitude' , and 'Bd infection intensity' as fixed factors, and added as nested random factor the categorical predictors 'location, sampling date, and host species' . We performed Generalized Linear Mixed Models (GLMM) using the glmer() function 96 in the lme4 package (https ://cran.r-proje ct.org/packa ge=lme4) in R, with Poisson distribution to test whether selected bioclimatic variables from Model 1-4, elevation, latitude, and Bd infection intensity are associated with richness and abundance of skin bacterial communities in treefrogs (i.e. observed #OTUs and Simpson's evenness; Supplementary Table 9 and 10).…”