Supplemental Information. Figure S1A. Lund-MacKay scores associated with disease state. No differences were observed between CRS, CRS+A or CRS+CF patients; B. CRS-CF patients are significantly younger than CRS+A patients (ANOVA, Tukey’s p = 0.041); however, no differences in age were observed for pairwise comparisons between the other groups (p > 0.05, Tukey’s post hoc test). Figure S2A. PCoA of an unweighted UniFrac distance matrix colored by DSI-IIIb and healthy (PERMANOVA p = 0.001, 18.2% variation explained); B. PCoA of weighted UniFrac distance matrix after dominant sequence reads associated with the dominant family in each sample were removed demonstrating that DS still significantly explains variation in community composition despite removal of the dominant taxon from each sample (PERMANOVA p = 0.001, 17.6% variation explained). Figure S3A. Laplace model fit demonstrates three distinct Dirichlet multinomial mixtures groups. B. Hierarchical cluster analysis using a weighted-UniFrac distance matrix showing that microbiomes enriched in Corynebacteriaceae forms a distinct cluster (au p = 100). Heatmap shows relative abundance of the bacterial general that comprise >90% of the total sequence reads. C. Reciprocal relationship between Corynebacteriaceae and Staphylococcaceae. Figure S4. PICRUSt-predicted functional variation across microbial Dirichlet states shows significant functional differences A. PCoA of Canberra distance matrix; PERMANOVA p = 0.001, 21.7% of variation explained) B. PCoA of Bray-Curtis distance matrix; PERMANOVA p = 0.001, 22.0% of variation explained). Figure S5. Expression levels of all host immune genes measured by QPCR (*indicates Kruskal-Wallis p