Fig. 4

Erm(C) variant phylogeny and diversity across three anatomical locations 60 individuals over 5 time point collections spanning 6-months. A Maximum likelihood phylogenetic analysis of erm(C) nucleic acid sequence variants. ASVs were aligned using MUSCLE v.3.8.1551 [46], inferred using the best-fit substitution model of HKY + F + G4 as determined by IQ-Tree [47], and visualized using iTOL [48]. Branch bootstrap support is gradient colored, with red representing 0% support and blue representing 100% support. Tree was rooted using ASV_1 which has 100% sequence similarity to reference erm(C) sequence [WP_001003263.1]. Phylogeny of variants shows distinct evolutionary subpopulations. Many of the variants shared high sequence similarity and thus were clustered together. Full tree without clustering is available in the (Additional file 4: Figs. S1 and S2). B Log10 total abundance of erm(C) ASVs summed across all time points and anatomical locations for a subject. Subject notation is indicated at the bottom. C Distribution of erm(C) variants across the three anatomical skin sites. Intersection size represents the number of ASVs that were present at the combination of locations, as noted by the scatter plot below the bar. Bars shown in blue (Set size) are the overall number of ASVs for each location