Follicular helper-like T-cell response at day 7. Projected taxa are biomarker taxa of the microbiota profiles identified by hierarchical clustering. Each data point represents the microbial community composition of a single sample at the time of LAIV receipt. OTUs are named by use of taxonomic annotations combined with a rank number in parentheses based on the average relative abundance of each OTU across the entire dataset to avoid multiple OTUs having identical annotations.R2refers to the percentage of variance explained by each immune variable, calculated by multivariate permutational ANOVA tests, adjusting for any covariates associated with both the microbiome at day 0 and the immune response to LAIV at either day 7 or day 21. NMDS also depicts biomarker taxa in relation to the ordination plot. The stress value for each ordination is 0.182, indicating a good goodness of fit of the data points following dimension reduction (<0.1 is excellent, <0.2 is good, and >0.2 is poor). LAIV=live attenuated influenza vaccine. NMDS=non-metric dimensional scaling. OTU=operational taxonomic unit.