Principal component plots show fungal microbiota for samples collected in both years, and Y2. Each symbol in a plot represents an individual sample and the color of each sample symbol indicates the facility in which a sample was collected. The presence and absence of viable L. monocytogenes , as determined by an enrichment method, is shown in triangles and circles, respectively. Filled symbols denote samples collected in Y1 and unfilled symbols squares denote samples collected in Y2. The p -values in panels B , C , E , and F from a two-way PERMANOVA model that included the effects of facility (FAC), presence of L. monocytogenes (LM), and their interaction (FAC:LM)