Introduction
Despite improvements, foodborne illness continues to pose a major threat to global public health (Gizaw, 2019). In the United States, it is estimated that pathogenic bacteria, such as nontyphoidal Salmonella spp. and Campylobacter, cause 3.6 million cases of foodborne illness each year (Scallan et al., 2011b). Animal products remain one of the top commodity groups associated with bacterial foodborne illness (64%), followed by plant products (32.1%) and seafood (3.9%) (Painter et al., 2013). National prevalence data from the 2021 calendar year estimated that Salmonella and Campylobacter prevalence on chicken parts was 6.69% and 16.20%, respectively, and Salmonella and Shiga toxin-producing Escherichia coli (STEC) prevalence in beef trim was 0.75% and 1.64%, respectively (USDA FSIS, 2022a). In the same year, there were a total of 47 recalls, 6 of which were due to Salmonella and STEC contamination in beef, pork, and chicken products (USDA FSIS, 2022b). With the estimated cost of a recall being $10 million, not including the lasting damage done to the brand and the impact on future sales, there is a clear incentive for food manufacturers to utilize antimicrobials or other pathogen inactivation strategies to mitigate these risks.
Antimicrobial interventions can be applied at various points during harvest, fabrication, and processing, most commonly in the form of sprays, dips, or immersion chilling systems (Winkler and Harris, 2009; Zhang et al., 2019). Peracetic acid (PAA) first gained popularity in the poultry industry as a replacement for chlorine. Historically chlorine was considered the industry standard for preventing cross-contamination; however, the efficacy of chlorine is heavily impacted by organic load and increases in pH (Byrd and McKee, 2005). Peracetic acid has become a more widely accepted antimicrobial with a broad range of approved usages (USDA FSIS, 2023). Most commercial PAA contains 3 chemicals, PAA, acetic acid, and hydrogen peroxide at varying concentrations, held in an aqueous solution with the addition of a stabilizer (Organic Materials Review Institute, 2000). The primary properties responsible for its antimicrobial efficacy are its acidic nature and ability to oxidize. However, studies have shown that the efficacy can vary between providers and formulation (Micciche et al., 2019; Cano et al., 2021).
Validation studies are considered the standard for evaluating the effect of specific antimicrobials on a microbial population (Shintani, 2015). However, while the primary goal of an antimicrobial is to reduce pathogenic bacteria, it is also important to consider the unintended effects it may have on the microbial community as a whole (Yang et al., 2021). In recent years, 16S rRNA gene sequencing has become a valuable tool to assess the changes in microbiota and the presence of potential spoilage organisms (Lee et al., 2017; Cauchie et al., 2020; Wythe et al., 2022). Therefore, the objective of this study was to determine the effect of PAA, at various concentrations, on the pathogen load and the microbiota of inoculated meat products when applied as a short-duration antimicrobial spray. The effect on pathogen load was assessed by plating for enumeration on selective media, and the microbiota was characterized using 16S rRNA gene sequencing. By evaluating the effects of PAA at various concentrations, this study aims to determine the optimal usage rates and characterize the microbial community after treatment.
Materials and Methods
Preparation of inoculum
Two cocktail cultures were prepared for the inoculation of chicken thighs and beef trim. The inoculum for chicken thighs consisted of 5 strains of Salmonella enterica (Typhimurium S9, Typhimurium M-09-0001A-1, Heidelberg S13, Enteritidis E40, and Enteritidis 6424) and one strain of Campylobacter jejuni (NCTC 11168). The inoculum for beef trim consisted of the same 5 strains of Salmonella enterica and 7 strains of Escherichia coli (E. coli; O26:H11, O111:H8, O103:H2, O45:H2, O145:NM, O121:H19, and O157:H7). All Salmonella and E. coli strains were from the University of Wisconsin-Madison Food Research Institute stock culture collection. Isolates were maintained in individual CRYOBANK™ vials (COPAN Diagnostics Inc., Murrieta, CA, USA) and stored at −80°C. Separate frozen cultures of Salmonella, E. coli, and Campylobacter were individually streaked for isolation on xylose lysine deoxycholate (XLD) agar (Difco, BD Biosciences, Sparks Glencoe, MD, USA), MacConkey agar (BD Biosciences), and modified Charcoal-Cefoperazone-Deoxycholate agar (mCCDA; HiMedia, Mumbai, India), respectively. The XLD and MacConkey plates were incubated aerobically at 37°C for 24 h, and mCCDA plates were incubated under microaerophilic conditions (5% O2, 10% CO2, 85% N2) at 42°C for 48 h using the Anoxomat III system (Advanced Instruments, Norword, MA, USA). A single colony from each XLD and MacConkey plate was transferred to 40 mL of TSB (BD Biosciences) and incubated aerobically for 12 h at 37°C. A single colony from the mCCDA plate was transferred to 40 mL of Bolton Broth (Oxoid Ltd, Altrincham, Cheshire, England) and incubated microaerophilically for 24 h at 42°C using the Anoxomat III system. The cultures were centrifuged at 13,500 × g for 2 min, decanted, and subsequently washed twice with 1× Phosphate Buffered Saline (PBS; VWR Life Science, Radnor, PA, USA). After the final wash, cocktails were prepared by resuspending individual pellets in 1× PBS, combining the resuspended pellets, and bringing the final volume to 40 mL. The final inoculation concentration of Salmonella and E. coli were 108 colony forming units (CFU)/g, whereas the final inoculation concentration of Campylobacter jejuni was 106 CFU/g.
Inoculation of meat products
A total of 26 skin-on, bone-in chicken thighs (n = 5, k = 5) and 26 pieces of beef trim (n = 5, k = 5) were used to conduct this experiment, including one uninoculated untreated sample (n = 1) for each product type to be used for microbiota analysis. Chicken thighs, weighing on average 186.9 ± 5.9 g, were inoculated with a cocktail of Salmonella and Campylobacter. Beef trim, with an average weight of 103.3 ± 0.7 g, was inoculated with a cocktail of Salmonella and E. coli. Both meat products were spot inoculated on the upward facing surface at a rate of 1 mL per 25 g, and the inoculum was then spread using the side of a sterile pipette tip and set aside to rest for 60 to 90 min at 4°C. Attachment levels of 107 CFU/mL of Salmonella and E. coli, and 105 CFU/mL of Campylobacter, were obtained.
Peracetic acid treatment
A (5 gal or 18.9 L) stock solution of 800 ppm was prepared using tap water and PAA (Hydrite Chemical Co., Brookfield, WI, USA). The stock solution was subsequently diluted down to generate a series of working solutions with a PAA concentration of 400 and 200 ppm. Concentrations were verified using a Hydrishield PA22HP titration kit (Hydrite Chemical Co., Brookfield, WI, USA). A modified spray cabinet (Figure 1) was utilized for the application of each spray treatment. The spray cabinet contained a single hook, in which the meat product was hung on, surrounded by 12 nozzles set at 30 psi. Inoculated chicken thighs and beef trim were either not treated with anything (NT) or subjected to one of the following levels of PAA: 0 (tap water), 200, 400, or 800 ppm. Chicken thighs were sprayed for 15 s, and beef trim was sprayed for 10 s to mimic industry standard practices and remain in accordance with FSIS Directive 7120.1 (USDA FSIS, 2023). After treatment, samples remained on the hook for 2 min to drip before being transferred to sterile rinse bags.
Microbial analyses
Chicken thighs were rinsed in 150 mL of neutralized buffered peptone water (nBPW; Hardy Diagnostics, Santa Maria, CA, USA) for 1 min before being aliquoted for microbiota analysis (40 mL) and serially diluted (1:10) in 1× PBS for the enumeration of Salmonella and Campylobacter (USDA FSIS, 2021). Beef trim was homogenized in 100 mL of DE Neutralizing broth (Remel, Lenexa, KS, USA) for 1 min at 200 rpm (Stomacher 80, Seward, West Sussex, UK) before being aliquoted for microbiota analysis and serially diluted (1:10) in 1× PBS for the enumeration of Salmonella and E. coli. Enumeration was performed by spread plating 100 μL of diluted homogenate onto XLD, MacConkey, and mCCDA for Salmonella, E. coli, and Campylobacter, respectively. XLD and MacConkey plates were incubated aerobically for 24 h at 37°C, and mCCDA plates were incubated microaerophilically for 48 h at 42°C before being manually counted within the range of 30 to 300 CFU (Breed and Dotterrer, 1916).
DNA extraction
Rinsate homogenates that were collected on the day of sampling were stored at −80°C until DNA extraction could occur. Rinsates were thawed at room temperature, pelleted for 10 min at 5000 × g, and total genomic DNA was extracted from the resulting pellet using the QIAGEN DNeasy Blood & Tissue Kit (Qiagen, Hilden, Germany) following the standard spin-column protocol with the recommended pretreatment for gram-negative bacteria. Extracted samples were eluted by hydrating the spin column with 50 μL of Buffer AE, incubating at ambient temperature for 10 min, and centrifuging for 1 min and 30 s at 8,000 rpm. All eluted samples were analyzed using a NanoQuant plate on an Infinite 200 PRO plate reader (Tecan Trading AG, Männedorf, Switzerland) to measure DNA concentration and 260/280. The eluted samples with a DNA concentration above 15 ng/μL were then diluted to 10 ng/μL in Buffer AE and stored at −20°C until the library could be prepared.
Library preparation
Using the custom primers designed by Kozich et al. (2013), a sequencing library was prepared based on the V4 regions of the 16S rRNA gene. DNA samples (2 μL) were amplified with dual-indexed primers (1 μL of 10 mM forward, 1 μL of 10 mM reverse), including 8 unique nucleotide barcode sequences, using a high-fidelity Accuprime Pfx DNA polymerase (19 μL; Invitrogen, Waltham, MA, USA) and dimethyl sulfoxide (2 μL; DMSO; Fisher Scientific, Waltham, MA, USA). Amplification was verified using gel electrophoresis (1.5%). PCR products were normalized in equimolar concentrations (20 μL) using a SequalPrep™ Normalization kit (Invitrogen, Waltham, MA, USA), and pooled libraries were created containing 5 μL of each normalized sample. Pooled library concentration was determined using a KAPA library quantification kit for Illumina platforms (Kapa Biosystems, Woburn, MA, USA) and a Qubit 1× dsDNA HS Assay Kit (Invitrogen, Waltham, MA, USA). The pooled library and PhiX Control v3 (Illumina, Carlsbad, CA, USA) were denatured with freshly made 0.2 N NaOH, incubated for 5 min, and subsequently diluted to 20 pM in Hyb buffer (Illumina, Carlsbad, CA, USA). The denatured 20 pM library and PhiX control were subsequently diluted in Hyb buffer to a final concentration of 6 pM. The resulting 6 pM library was combined with the 6 pM PhiX control (30%, v/v) and loaded into a 2 × 250 cycle v2 MiSeq cartridge (Illumina, Carlsbad, CA, USA). The resulting sequences were uploaded into BaseSpace (Illumina, Carlsbad, CA, USA) and demultiplexed. Demultiplexed reads were uploaded into NCBI (PRJNA906192) and github (https://github.com/RickeLab-UW/Microbiota-of-Meat-Treated-with-PAA).
Statistical and bioinformatic analyses
At the onset of the study, inoculated chicken thighs and beef trim were randomly assigned to a treatment group. The CFU of Salmonella, Campylobacter, and E. coli were Log10 transformed and reported on a CFU/g of meat basis (Dittoe et al., 2019). The logarithmic data were analyzed in RStudio (Version 4.2.1) as a linear model using a one-way ANOVA with PAA treatment as a categorical main effect (R Core Team, 2022). Pairwise differences were determined using Tukey’s Protected honestly significant differences (HSD) with effects considered significant at P ≤ 0.05. Further analysis was performed by filtering out the NT control observations and analyzing the remainder of the data in a linear regression model using PAA treatment as a continuous fixed effect. The main effects were considered significant at P ≤ 0.05.
Sequencing data were analyzed using the QIIME2 pipeline (version 2022.8). Demultiplexed paired-end reads were downloaded from Illumina BaseSpace (Illumina, Carlsbad, CA, USA), formatted, and imported into the QIIME2 pipeline using Casava 1.8 (Bolyen et al., 2019). Imported sequences were filtered and denoised for quality with DADA2 via the q2-dada2 plugin (Callahan et al., 2016). Alignment of the Amplicon Sequence Variant (ASV) was performed using MAFFT, a multiple alignment program (version 7), and rooted and unrooted phylogenetic trees were generated with fasttree2 using the q2-phylogeny plugin (Price et al., 2010; Katoh and Standley, 2013). The ASV were then classified using a SILVA classifier (Silva 138 99% OTU full-length sequences; MD5: b8609f23e9b17bd4a1321a8971303310) with a specified confidence level of 95% via q2-feature-classifier plugin (Quast et al., 2013; Yilmaz et al., 2014; Bokulich et al., 2018). Taxonomic based filtering was performed on the feature tables to exclude chloroplasts and mitochondria from downstream analysis. The α- and β-diversity core metrics were generated through the q2-diversity plugin (McKinney, 2010). Main effect and pairwise comparisons of α-diversity metrics, Shannon’s diversity index, and Pielou’s Evenness were determined using Kruskal-Wallis (Shannon, 1948; Kruskal and Wallis, 1952; Pielou, 1966). Additionally, main effect and pairwise comparison of β-diversity metrics, Bray-Curtis, and Weighted UniFrac were determined using ANOSIM (Sørensen, 1948; Lozupone and Knight, 2005; Lozupone et al., 2007; Hamady et al., 2010; Lozupone et al., 2011; McDonald et al., 2018). Subsequently, an analysis of the composition of microbiomes (ANCOM) was performed to distinguish between differentially abundant taxa using the q2-composition plugin (Mandal et al., 2015). Taxonomic means and α-diversity tables were visualized in Microsoft Excel (Version 16.66.1) (Microsoft Corporation, 2018). Main effects were considered significant at P ≤ 0.05 and pairwise differences at Q ≤ 0.05.
Results
Chicken
Treatment had a significant effect on the survival of Salmonella and Campylobacter on skin-on bone-in chicken thighs (P < 0.01; Figure 2). The application of 200, 400, and 800 ppm PAA (6.38, 6.37, and 6.03 Log CFU/g) resulted in lower levels of Salmonella compared to the NT control (7.14 Log CFU/g; P < 0.05). Thighs treated with 200 and 400 ppm PAA (6.38 and 6.37 Log CFU/g) did not result in differing Salmonella loads compared to those treated with 0 ppm PAA (tap water; 6.79 Log CFU/g; P > 0.05). Application of 800 ppm PAA (6.03 Log CFU/g) was the only treatment that resulted in a lower level of Salmonella compared to the thighs treated with 0 ppm PAA (6.79 Log CFU/g; P < 0.05). Compared to the NT control (5.33 Log CFU/g), treatment with 0 (tap water), 200, 400, and 800 ppm PAA (4.71, 4.31, 4.17, and 3.79 Log CFU/g) resulted in lower levels of Campylobacter (P < 0.05). The application of 200 ppm PAA (4.31 Log CFU/g) did not result in differing levels of Campylobacter compared to thighs treated with 0 ppm PAA (4.71 Log CFU/g; P > 0.05). However, unlike Salmonella, the application of both 400 and 800 ppm PAA (4.17 and 3.79 Log CFU/g) resulted in lower levels of Campylobacter compared to those treated with 0 ppm PAA (4.71 Log CFU/g; P < 0.05). To determine the linear relationship between pathogen load and PAA concentration, the data were analyzed without the NT control. Significant negative relationships (P < 0.01) were detected between PAA concentration and the load of both Salmonella (R2 = 0.44) and Campylobacter (R2 = 0.55; Supplemental Figure 1). With every 100 ppm increase in PAA concentrations, up to 800 ppm, a reduction of 0.09 and 0.11 Log CFU/g of Salmonella and Campylobacter, respectively, would be expected.
Treatment did not have a significant impact on the microbial diversity of the inoculated chicken thighs. There was no main effect of treatment on either the α-diversity, Shannon’s diversity index (P = 0.72), and Pielou’s Evenness (P = 0.88) or on the β-diversity, Bray-Curtis (P = 0.71), and Weighted Unifrac (P = 0.72; Figure 3; Supplemental Tables 1 and 2). The ANCOM analysis determined that there were no significantly different abundant taxa at any taxonomic level. At the phyla level, the main taxonomic group in the inoculated chicken thighs was Proteobacteria (Figure 4). In contrast, Firmicutes and Proteobacteria were the main taxa present in the uninoculated chicken thighs (Supplemental Figure 3).
Beef
There was a significant treatment effect on the survivability of Salmonella and E. coli on beef trim (P < 0.01; Figure 5). The application of 0, 200, 400, and 800 ppm PAA (6.69, 6.72, 6.74, and 6.47 Log CFU/g) all resulted in lower levels of Salmonella compared to the NT control (7.25 Log CFU/g; P < 0.05). However, the load of Salmonella on beef trim treated with 0 ppm PAA (6.69 Log CFU/g) did not differ from those treated with 200, 400, or 800 ppm PAA (6.72, 6.74, and 6.47 Log CFU/g; P > 0.05). The E. coli results followed a similar trend in which the application of 0, 200, 400, and 800 ppm PAA (7.06, 7.04, 6.95, and 6.83 Log CFU/g) resulted in lower levels of E. coli compared to the NT beef trim (7.55 Log CFU/g; P < 0.05), while no differences were observed between the trim treated with 0, 200, 400, and 800 ppm PAA (P > 0.05). However, significant linear relationships (P < 0.05) were observed between PAA concentration and pathogen load, Salmonella (R2 = 0.18) and E. coli (R2 = 0.17), when data were analyzed without the inclusion of the NT control (Supplemental Figure 2). A reduction of 0.03 Log CFU/g is expected for both Salmonella and E. coli with every 100 ppm increase in PAA concentration.
There was no main effect of PAA treatment on the α-diversity, Shannon’s diversity index (P = 0.32), or Pielou’s Evenness (P = 0.28; Figure 6; Supplemental Table 3). Treatment was not found to have a significant effect on the β-diversity metric Weighted Unifrac (P = 0.24); however, a significant effect of treatment was observed on the β-diversity metric Bray-Curtis (P = 0.01; Figure 6; Supplemental Table 4). In the Bray-Curtis distance metric, significant pairwise differences were observed between 0 and 400 ppm PAA and the 0 and 800 ppm PAA treatments (P < 0.01, Q < 0.05). Using ANCOM, it was determined that there were no significantly different abundant taxa at any taxonomic level. Proteobacteria was more than 80% of the microbial community present on the inoculated beef trim regardless of treatment (Figure 7), while the microbial community of the uninoculated beef trim primarily consisted of Firmicutes and Actinobacteriota (Supplemental Figure 3).
Discussion
Numerous studies have demonstrated the efficacy of PAA to reduce the microbial load on meat and poultry products; however, the reported results have varied greatly based on the target organism and product type (Ellebracht et al., 2005; Kim et al., 2017; Ramirez-Hernandez et al., 2018; Zhang et al., 2019; Britton et al., 2020; Kataria et al., 2020; Kumar et al., 2020a). The majority of the research thus far has focused on the application of PAA during poultry processing, evaluating various delivery systems, concentrations, and exposure times (Cano et al., 2021). Immersion interventions (dip/chill tank) have been shown to be the most effective application of PAA, as they generally allow for the longest exposure time and most complete coverage (Cano et al., 2021; Gonzalez et al., 2021). Despite the increased lethality associated with immersion interventions, there are situations in which it is more cost-effective and/or operationally feasible to implement a PAA spray rather than a dip. The current study aimed to evaluate the antimicrobial effects of a short-duration PAA spray on Salmonella, Campylobacter, E. coli, and the overall microbiota of inoculated chicken thighs and beef trim. In general, a short-duration spray effectively lowered the level of inoculated pathogens compared to the NT control, and as the concentration of PAA increased, the level of pathogens was also reduced.
While the application of a spray, at all tested levels (0, 200, 400, and 800 ppm), successfully reduced the load of Campylobacter on chicken and Salmonella and E. coli on beef compared to the untreated control group, a rinsing effect of approximately 0.5 Log CFU/g was observed when the inoculated meat products were treated with 0 ppm PAA (tap water). Similar effects were reported by Ramirez-Hernandez et al. (2018) where a 0.3 Log CFU/mL reduction in Salmonella was observed after skin-on chicken thighs were sprayed with tap water for 15 s. In another study, beef trim was sprayed for 10 s with tap water, and a 0.4 Log CFU/cm2 was reported (Britton et al., 2020). For a PAA treatment to be considered effective, the level of recoverable pathogens after treatment must be significantly lower than that of those treated with tap water. Using this definition, only the 800 ppm PAA treatment effectively reduced Salmonella on chicken thighs and the 400 and 800 ppm PAA treatment effectively reduced Campylobacter. The inability of a short-duration PAA spray to effectively reduce Salmonella, at concentrations at or below 400 ppm, has been previously documented. Ramirez-Hernandez et al. (2018) found no significant differences in the amount of recoverable Salmonella on chicken thighs treated with either 200 or 400 ppm PAA compared to those treated with tap water. However, Campylobacter has been shown to be less tolerant of PAA exposure. Joo et al. (2020) reported that PAA could effectively reduce Campylobacter at concentrations as low as 50 ppm when chicken skins were fully submerged for 5 min. Extended exposure time and delivery method (submersion) likely contributed to the superior antimicrobial efficacy demonstrated in this study. However, the current study and others have shown that applying PAA as a short-duration antimicrobial spray can also be an effective delivery system to reduce Campylobacter populations on chicken products (Gonzalez et al., 2021).
The application of PAA as a short-duration antimicrobial spray has been shown to be viable intervention strategy for the poultry industry; however, the results from this study suggests that it may have limited benefits when applied to the surface of beef trim (Ramirez-Hernandez et al., 2018; Kumar et al., 2020b; Gonzalez et al., 2021). While most of the validation work has been focused on poultry, several studies have evaluated the efficacy of PAA on both the surface of beef carcasses and on beef trim (Gill and Badoni, 2004; Ellebracht et al., 2005; King et al., 2005; Britton et al., 2020). Unlike poultry, these studies offer somewhat contradictory results. Some studies report that on beef trim the application of a 350 and 400 ppm PAA spray (10 s) can reduce Enterobacteriaceae populations by 1.7 and 1.9 Log CFU/cm2, respectively (Britton et al., 2020), while others report reductions of 1 Log CFU/cm2 or less of E. coli O157:H7 and S. Typhimurium after submerging fresh beef trim in various concentrations of PAA for 15 s. Another study evaluated the effects of a 15 s PAA spray applied to the surface of beef carcasses at two temperatures (45°C and 55°C) and found that even at the highest concentration (1000 ppm), they were only able to reduce E. coli O157:H7 and S. Typhimurium by 0.9 and 0.6 Log CFU/cm2 at 45°C and 1.7 and 1.3 Log CFU/cm2 at 55°C, respectively (King et al., 2005). For comparison, the same study evaluated the efficacy of a 4% lactic acid spray at 55°C and achieved reductions of 2.7 and 3.4 Log CFU/cm2 of E. coli O157:H7 and S. Typhimurium, respectively (King et al., 2005). Hot water wash, acidified sodium chlorite, and organic acids (lactic, acetic, citric, etc.) are all common antimicrobial interventions applied either independently or in combination with other antimicrobial treatments to reduce the pathogen load on meat products (Mani-López et al., 2012; Cano et al., 2021). While each of these interventions may have value when used independently, great success has been observed using a combined approach (Leistner and Gorris, 1995; Dittoe et al., 2019; Joo et al., 2020; Olson et al., 2020). When choosing an antimicrobial intervention strategy, it is important to maintain the balance between product safety and product quality. While shelf life and quality were not directly evaluated in this study, PAA acid has been shown to both positively and negatively impact both attributes. Yang et al. (2021) suggest that by reducing the initial microbial load and favoring the growth of lactic acid–producing bacteria, PAA can extend the shelf life. However, other studies have observed that applying PAA at high concentrations has the potential to negatively impact quality attributes, such as color and pH, on both beef and poultry products (Ellebracht et al., 2005; Bauermeister et al., 2008; Chen et al., 2014).
Traditionally, targeted pathogen inactivation studies have been considered the golden standard to validate the efficacy of an antimicrobial intervention. However, in recent years the increased accessibility to next-generation sequencing (NGS) technologies has enabled researchers to look beyond the traditional targeted approach and examine the impact of antimicrobials on the microbial ecosystem as a whole (Park et al., 2023). Microbiota analysis has proved to be a valuable tool for process biomapping, in-line intervention validation, and evaluating the effect of storage conditions on shelf life (Handley et al., 2018; Wages et al., 2019; Weinroth et al., 2019). Yang et al. (2021) demonstrated the utility of this approach from a shelf-life perspective by identifying unique genera associated with PAA treatment, in addition to characterizing the changes in the microbial community of vacuum-packed beef subprimal cuts after treatment and storage. The inclusion of microbiota analysis into traditional validation studies has not come without challenges. Validation studies often require products or surfaces to be inoculated with a target organism, at an artificially high level (∼108 CFU/mL), to demonstrate a measurable reduction (National Advisory Committee on Microbiological Criteria for Foods, 2010). Unsurprisingly, studies have shown that this type of inoculation can have significant impacts on the microbial composition, making it difficult to determine the true treatment effect on the native microbial population (Liu et al., 2018; Wythe et al., 2022). Firmicutes, Proteobacteria, and Actinobacteriota make up a large majority of the genera found on raw poultry and beef product (Weinroth et al., 2019; Park et al., 2023). After inoculation, the proportionality of the microbial community shifts and becomes predominantly composed of Proteobacteria (Wythe et al., 2022). Our results showed no effect of treatment on α-diversity (evenness or richness) among inoculated samples. Similar results were observed by Wythe et al. (2022), in which treatment did not impact the richness or evenness within the inoculated treatment groups; however, the microbial population of the non-inoculated not treated group was found to be richer and more even. The same study reported minimal variation in β-diversity among inoculated treatments, further supporting our findings. This lack of treatment effect on the microbial composition may be partially explained by the inability of 16S rRNA gene sequencing to distinguish between viable and non-viable DNA molecules. This inability to differentiate viability status may result in the overestimation of a microbial population (Rudi et al., 2005). As the public health concern surrounding food safety continues to rise, the demand for rapid, sensitive, and accurate pathogen detection methods is higher than ever (Scallan et al., 2011a).
Conclusions
The application of a short-duration antimicrobial spray resulted in a reduction of target pathogens on chicken and beef. However, variations in efficacy were observed between target pathogens and meat matrixes. Campylobacter exhibited the least resistance to PAA treatment, resulting in the greatest reductions of the artificially inoculated population. The efficacy of PAA toward Salmonella appeared to be more strongly correlated to the meat matrices it was applied to. Peracetic acid treatments were not effective at reducing the population of Salmonella and E. coli on beef trim; however, they were somewhat effective at reducing the level of Salmonella and Campylobacter on chicken thighs. These findings suggest that certain properties of beef trim may be providing a protective advantage to foodborne pathogens, enabling them to persist on the product after antimicrobial intervention. In the future, NGS technologies will become a valuable tool in assessing the true impact of antimicrobial interventions on both target organisms and the microbial community as a whole; however, further refinement associated viability and quantitation will be necessary to unlock its full potential. Using this information in combination with traditional validation studies will enable processors to better optimize chemical usage and improve their process control.
Acknowledgements
The authors would like to thank Hydrite Chemical Co. for providing technical and financial support throughout this study. Author J. A. B. would like to acknowledge the University of Wisconsin Foundation at the University of Wisconsin-Madison for its continued support through the Bray-Woodbury Wisconsin Distinguished Graduate Fellowship Fund in Meat Science. Authors K. B. M. and B. H. are employed by Hydrite Chemical Co. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationship that could be construed as a potential conflict of interest.
Literature Cited
Bauermeister, L. J., J. W. J. Bowers, J. C. Townsend, and S. R. McKee. 2008. The microbial and quality properties of poultry carcasses treated with peracetic acid as an antimicrobial treatment. Poultry Sci. 87:2390–2398. doi: https://doi.org/10.3382/ps.2008-00087
Bokulich, N. A., B. D. Kaehler, J. R. Rideout, M. Dillon, E. Bolyen, R. Knight, G. A. Huttley, and J. G. Caporaso. 2018. Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s q2-feature-classifier plugin. Microbiome 6(1). doi: https://doi.org/10.1186/s40168-018-0470-z
Bolyen, E., J. R. Rideout, M. R. Dillon, N. A. Bokulich, C. C. Abnet, G. A. Al-Ghalith, H. Alexander, E. J. Alm, M. Arumugam, F. Asnicar, Y. Bai, J. E. Bisanz, K. Bittinger, A. Brejnrod, C. J. Brislawn, C. T. Brown, B. J. Callahan, A. M. Caraballo-Rodríguez, J. Chase, E. K. Cope, R. Da Silva, C. Diener, P. C. Dorrestein, G. M. Douglas, D. M. Durall, C. Duvallet, C. F. Edwardson, M. Ernst, M. Estaki, J. Fouquier, J. M. Gauglitz, S. M. Gibbons, D. L. Gibson, A. Gonzalez, K. Gorlick, J. Guo, B. Hillmann, S. Holmes, H. Holste, C. Huttenhower, G. A. Huttley, S. Janssen, A. K. Jarmusch, L. Jiang, B. D. Kaehler, K. Bin Kang, C. R. Keefe, P. Keim, S. T. Kelley, D. Knights, I. Koester, T. Kosciolek, J. Kreps, M. G. I. Langille, J. Lee, R. Ley, Y.-X. Liu, E. Loftfield, C. Lozupone, M. Maher, C. Marotz, B. D. Martin, D. McDonald, L. J. McIver, A. V. Melnik, J. L. Metcalf, S. C. Morgan, J. T. Morton, A. T. Naimey, J. A. Navas-Molina, L. F. Nothias, S. B. Orchanian, T. Pearson, S. L. Peoples, D. Petras, M. L. Preuss, E. Pruesse, L. B. Rasmussen, A.Rivers, M. S. Robeson II, P. Rosenthal, N. Segata, M. Shaffer, A. Shiffer, R. Sinha, S. J. Song, J. R. Spear, A. D. Swafford, L. R. Thompson, P. J. Torres, P. Trinh, A. Tripathi, P. J. Turnbaugh, S. Ul-Hasan, J. J. J. van der Hooft, F. Vargas, Y. Vázquez-Baeza, E. Vogtmann, M. von Hippel, W. Walters, Y. Wan, M. Wang, J. Warren, K. C. Weber, C. H. D. Williamson, A. D. Willis, Z. Z. Xu, J. R. Zaneveld, Y. Zhang, Q. Zhu, R. Knight, and J. G. Caporaso. 2019. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 37:852–857. doi: https://doi.org/10.1038/s41587-019-0209-9
Breed, R. S., and W. D. Dotterrer. 1916. The number of colonies allowable on satisfactory agar plates. J. Bacteriol. 1:321–331. doi: https://doi.org/10.1128/jb.1.3.321-331.1916
Britton, B. C., I. Geornaras, J. O. Reagan, S. Mixon, D. R. Woerner, and K. E. Belk. 2020. Antimicrobial efficacy of acidified peroxyacetic acid treatments against surrogates for enteric pathogens on prerigor beef. Meat Muscle Biol. 4(1). doi: https://doi.org/10.22175/mmb.10992
Byrd, J. A., and S. R. McKee. 2005. Improving slaughter and processing technologies. In: G. C. Mead, editor, Food safety control in the poultry industry (Woodhead Publishing Series in Food Science, Technology and Nutrition). p. 310–332. doi: https://doi.org/10.1533/9781845690236.310
Callahan, B. J., P. J. McMurdie, M. J. Rosen, A. W. Han, A. J. A. Johnson, and S. P. Holmes. 2016. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 13:581–583. doi: https://doi.org/10.1038/nmeth.3869
Cano, C., Y. Meneses, and B. D. Chaves. 2021. Application of peroxyacetic acid for decontamination of raw poultry products and comparison to other commonly used chemical antimicrobial interventions: A review. J. Food Protect. 84:1772–1783. doi: https://doi.org/10.4315/JFP-21-107
Cauchie, E., L. Delhalle, B. Taminiau, A. Tahiri, N. Korsak, S. Burteau, P. A. Fall, F. Farnir, G. Baré, and G. Daube. 2020. Assessment of spoilage bacterial communities in food wrap and modified atmospheres-packed minced pork meat samples by 16S rDNA metagenetic analysis. Front. Microbiol. 10:1–17. doi: https://doi.org/10.3389/fmicb.2019.03074
Chen, X., L. J. Bauermeister, G. N. Hill, M. Singh, S. F. Bilgili, and S. R. Mckee. 2014. Efficacy of various antimicrobials on reduction of Salmonella and Campylobacter and quality attributes of ground chicken obtained from poultry parts treated in a postchill decontamination tank. J. Food Protect. 77:1882–1888. doi: https://doi.org/10.4315/0362-028X.JFP-14-114
Dittoe, D. K., J. A. Atchley, K. M. Feye, J. A. Lee, C. J. Knueven, and S. C. Ricke. 2019. The efficacy of sodium bisulfate salt (SBS) alone and combined with peracetic acid (PAA) as an antimicrobial on whole chicken drumsticks artificially inoculated with Salmonella Enteritidis. Frontiers in Veterinary Science 6:6. doi: https://doi.org/10.3389/fvets.2019.00006
Ellebracht, J. W., D. A. King, A. Castillo, L. M. Lucia, G. R. Acuff, K. B. Harris, and J. W. Savell. 2005. Evaluation of peroxyacetic acid as a potential pre-grinding treatment for control of Escherichia coli O157:H7 and Salmonella Typhimurium on beef trimmings. Meat Sci. 70:197–203. doi: https://doi.org/10.1016/j.meatsci.2005.01.003
Gill, C. O., and M. Badoni. 2004. Effects of peroxyacetic acid, acidified sodium chlorite or lactic acid solutions on the microflora of chilled beef carcasses. Int. J. Food Microbiol. 91:43–50. doi: https://doi.org/10.1016/S0168-1605(03)00329-5
Gizaw, Z. 2019. Public health risks related to food safety issues in the food market: A systematic literature review. Environ. Health Prev. 24. doi: https://doi.org/10.1186/s12199-019-0825-5
Gonzalez, S. V., I. Geornaras, M. N. Nair, and K. E. Belk. 2021. Evaluation of immersion and spray applications of antimicrobial treatments for reduction of Campylobacter jejuni on chicken wings. Foods 10:903. doi: https://doi.org/10.3390/foods10040903
Hamady, M., C. Lozupone, and R. Knight. 2010. Fast UniFrac: Facilitating high-throughput phylogenetic analyses of microbial communities including analysis of pyrosequencing and PhyloChip data. ISME J. 4:17–27. doi: https://doi.org/10.1038/ismej.2009.97
Handley, J. A., S. H. Park, S. A. Kim, and S. C. Ricke. 2018. Microbiome profiles of commercial broilers through evisceration and immersion chilling during poultry slaughter and the identification of potential indicator microorganisms. Front. Microbiol. 9:345. doi: https://doi.org/10.3389/fmicb.2018.00345
Joo, H., M. F. R. Mizan, M. I. Hossain, D. Lee, and S. Ha. 2020. Enhanced elimination of Salmonella Typhimurium and Campylobacter jejuni on chicken skin by sequential exposure to ultrasound and peroxyacetic acid. J. Food Safety 40:1–7. doi: https://doi.org/10.1111/jfs.12803
Kataria, J., S. Vaddu, E. N. Rama, G. Sidhu, H. Thippareddi, and M. Singh. 2020. Evaluating the efficacy of peracetic acid on Salmonella and Campylobacter on chicken wings at various pH levels. Poultry Sci. 99:5137–5142. doi: https://doi.org/10.1016/j.psj.2020.06.070
Katoh, K., and D. M. Standley. 2013. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 30:772–780. doi: https://doi.org/10.1093/molbev/mst010
Kim, S. A., S. H. Park, S. I. Lee, C. M. Owens, and S. C. Ricke. 2017. Assessment of chicken carcass microbiome responses during processing in the presence of commercial antimicrobials using a next generation sequencing approach. Sci. Rep.-UK 7:43354. doi: https://doi.org/10.1038/srep43354
King, D. A., L. M. Lucia, A. Castillo, G. R. Acuff, K. B. Harris, and J. W. Savell. 2005. Evaluation of peroxyacetic acid as a post-chilling intervention for control of Escherichia coli O157:H7 and Salmonella Typhimurium on beef carcass surfaces. Meat Sci. 69:401–407. doi: https://doi.org/10.1016/j.meatsci.2004.08.010
Kozich, J. J., S. L. Westcott, N. T. Baxter, S. K. Highlander, and P. D. Schloss. 2013. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl. Environ. Microb. 79:5112–5120. doi: https://doi.org/10.1128/AEM.01043-13
Kruskal, W. H., and W. A. Wallis. 1952. Use of ranks in one-criterion variance analysis. J. Am. Stat. Assoc. 47:583–621. doi: https://doi.org/10.1080/01621459.1952.10483441
Kumar, H., Y. Jang, K. Kim, J. Park, M. Jung, and J.-E. Park. 2020a. Compositional and functional characteristics of swine slurry microbes through 16S rRNA metagenomic sequencing approach. Animals 10:1372. doi: https://doi.org/10.3390/ani10081372
Kumar, S., M. Singh, D. E. Cosby, N. A. Cox, and H. Thippareddi. 2020b. Efficacy of peroxy acetic acid in reducing Salmonella and Campylobacter spp. populations on chicken breast fillets. Poultry Sci. 99:2655–2661. doi: https://doi.org/10.1016/j.psj.2019.12.045
Lee, H., M. Kwon, S. Heo, M. Kim, and G. Kim. 2017. Characterization of the biodiversity of the spoilage microbiota in chicken meat using next generation sequencing and culture dependent approach. Korean J. Food Sci. An. 37:535–541. doi: https://doi.org/10.5851/kosfa.2017.37.4.535
Leistner, L., and L. G. M. Gorris. 1995. Food preservation by hurdle technology. Trends Food Sci. Tech. 6:41–46. doi: https://doi.org/10.1016/S0924-2244(00)88941-4
Liu, L., L. Lin, L. Zheng, H. Tang, X. Fan, N. Xue, M. Li, M. Liu, and X. Li. 2018. Cecal microbiome profile altered by Salmonella enterica, serovar Enteritidis inoculation in chicken. Gut Pathog. 10. doi: https://doi.org/10.1186/s13099-018-0261-x
Lozupone, C., and R. Knight. 2005. UniFrac: a new phylogenetic method for comparing microbial communities. Appl. Environ. Microb. 71:8228–8235. doi: https://doi.org/10.1128/AEM.71.12.8228-8235.2005
Lozupone, C. A., M. Hamady, S. T. Kelley, and R. Knight. 2007. Quantitative and qualitative β diversity measures lead to different insights into factors that structure microbial communities. Appl. Environ. Microb. 73:1576–1585. doi: https://doi.org/10.1128/AEM.01996-06
Lozupone, C. A., M. E. Lladser, D. Knights, J. Stombaugh, and R. Knight. 2011. UniFrac: an effective distance metric for microbial community comparison. ISME J. 5:169–172. doi: https://doi.org/10.1038/ismej.2010.133
Mandal, S., W. Van Treuren, R. A. White, M. Eggesbø, R. Knight, and S. D. Peddada. 2015. Analysis of composition of microbiomes: A novel method for studying microbial composition. Microb. Ecol. Health D. 26:27663. doi: https://doi.org/10.3402/mehd.v26.27663
Mani-López, E., H. S. García, and A. López-Malo. 2012. Organic acids as antimicrobials to control Salmonella in meat and poultry products. Food Res. Int. 45:713–721. doi: https://doi.org/10.1016/j.foodres.2011.04.043
McDonald, D., Y. Vázquez-Baeza, D. Koslicki, J. McClelland, N. Reeve, Z. Xu, A. Gonzalez, and R. Knight. 2018. Striped UniFrac: enabling microbiome analysis at unprecedented scale. Nat. Methods 15:847–848. doi: https://doi.org/10.1038/s41592-018-0187-8
McKinney, W. 2010. Data structures for statistical computing in python. In: S. van der Walt and J. Millman, editors, Proc. 9th Python Sci. Conf. p. 51–56.
Micciche, A. C., K. M. Feye, P. M. Rubinelli, J. A. Lee, C. J. Knueven, and S. C. Ricke. 2019. Comparison of acid sanitizers on Salmonella Typhimurium inoculated commercial poultry processing reuse water. Frontiers in Sustainable Food Systems 2:90. doi: https://doi.org/10.3389/fsufs.2018.00090
Microsoft Corporation. 2018. Microsoft Excel. https://office.microsoft.com/excel
National Advisory Committee on Microbiological Criteria for Foods. 2010. Parameters for determining inoculated pack/challenge study protocols. J. Food Protect. 73:140–202. doi: https://doi.org/10.4315/0362-028x-73.1.140
Olson, E. G., L. A. Wythe, D. K. Dittoe, K. M. Feye, and S. C. Ricke. 2020. Application of Amplon in combination with peroxyacetic acid for the reduction of nalidixic acid–resistant Salmonella Typhimurium and Salmonella Reading on skin-on, bone-in tom turkey drumsticks. Poultry Sci. 99:6997–7003. doi: https://doi.org/10.1016/j.psj.2020.08.078
Organic Materials Review Institute. 2000. Peracetic acid. https://www.ams.usda.gov/sites/default/files/media/Peracetic Acid Technical Report Handling.pdf. (Accessed 16 June 2023).
Painter, J. A., R. M. Hoekstra, T. Ayers, R. V. Tauxe, C. R. Braden, F. J. Angulo, and P. M. Griffin. 2013. Attribution of foodborne illnesses, hospitalizations, and deaths to food commodities by using outbreak data, United States, 1998–2008. Emerg. Infect. Dis. 19:407–415. doi: https://doi.org/10.3201/eid1903.111866
Park, J., D. Bae, and S. A. Kim. 2023. Microbial trace investigation throughout the entire chicken supply chain based on metagenomic high-throughput sequencing. Food Res. Int. 169:112775. doi: https://doi.org/10.1016/j.foodres.2023.112775
Pielou, E. C. 1966. The measurement of diversity in different types of biological collections. J. Theor. Biol. 13:131–144. doi: https://doi.org/10.1016/0022-5193(66)90013-0
Price, M. N., P. S. Dehal, and A. P. Arkin. 2010. FastTree 2—Approximately maximum-likelihood trees for large alignments. PLoS ONE 5:e9490. doi: https://doi.org/10.1371/journal.pone.0009490
Quast, C., E. Pruesse, P. Yilmaz, J. Gerken, T. Schweer, P. Yarza, J. Peplies, and F. O. Glockner. 2013. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucleic Acids Res. 41:D590–D596. doi: https://doi.org/10.1093/nar/gks1219
R Core Team. 2022. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria.
Ramirez-Hernandez, A., M. M. Brashears, and M. X. Sanchez-Plata. 2018. Efficacy of lactic acid, lactic acid–acetic acid blends, and peracetic acid to reduce Salmonella on chicken parts under simulated commercial processing conditions. J. Food Protect. 81:17–24. doi: https://doi.org/10.4315/0362-028X.JFP-17-087
Rudi, K., B. Moen, S. M. Drømtorp, and A. L. Holck. 2005. Use of ethidium monoazide and PCR in combination for quantification of viable and dead cells in complex samples. Appl. Environ. Microb. 71:1018–1024. doi: https://doi.org/10.1128/AEM.71.2.1018-1024.2005
Scallan, E., P. M. Griffin, F. J. Angulo, R. V. Tauxe, and R. M. Hoekstra. 2011a. Foodborne illness acquired in the United States—Unspecified agents. Emerg. Infect. Dis. 17:16–22. doi: https://doi.org/10.3201/eid1701.P21101
Scallan, E., R. M. Hoekstra, F. J. Angulo, R. V. Tauxe, M.-A. Widdowson, S. L. Roy, J. L. Jones, and P. M. Griffin. 2011b. Foodborne illness acquired in the United States—Major pathogens. Emerg. Infect. Dis. 17:7–15. doi: https://doi.org/10.3201/eid1701.P11101
Shannon, C. E. 1948. A mathematical theory of communication. Bell Syst. Tech. J. 27:379–423. doi: https://doi.org/10.1002/j.1538-7305.1948.tb01338.x
Shintani, H. 2015. Validation studies for microbial contamination and control of contaminants. Biocontrol Sci. 20:161–170. doi: https://doi.org/10.4265/bio.20.161
Sørensen, T. 1948. A method of establishing groups of equal amplitude in plant sociology based on similarity of species and its application to analyses of the vegetation on Danish commons. Kong. Dansk. Vidensk. 5:1–34.
USDA FSIS. 2021. Sampling instructions: Salmonella and Campylobacter verification program for raw poultry products. 2021. https://www.fsis.usda.gov/sites/default/files/media_file/2021-03/10250.1_0.pdf. (Accessed 25 September 2023).
USDA FSIS. 2022a. Data collection and reports: Sampling results for FSIS regulaed products. 2022. https://www.fsis.usda.gov/science-data/sampling-program/sampling-results-fsis-regulated-products. (Accessed 16 June 2023).
USDA FSIS. 2022b. Summary of recall cases in calendar year 2021. 2022. https://www.fsis.usda.gov/food-safety/recalls-public-health-alerts/annual-recall-summaries/summary-recall-cases-calendar-8. (Accessed 16 June 2023).
USDA FSIS. 2023. Safe and suitable ingredients used in the production of meat, poultry, and egg products. Revision 58. https://www.fsis.usda.gov/sites/default/files/media_file/2021-09/7120.1.pdf. (Accessed 25 September 2023).
Wages, J. A., K. M. Feye, S. H. Park, S. A. Kim, and S. C. Ricke. 2019. Comparison of 16S rDNA next sequencing of microbiome communities from post-scalder and post-picker stages in three different commercial poultry plants processing three classes of broilers. Front. Microbiol. 10:972. doi: https://doi.org/10.3389/fmicb.2019.00972
Weinroth, M. D., B. C. Britton, K. R. McCullough, J. N. Martin, I. Geornaras, R. Knight, K. E. Belk, and J. L. Metcalf. 2019. Ground beef microbiome changes with antimicrobial decontamination interventions and product storage. PLoS One 14:e0217947. doi: https://doi.org/10.1371/journal.pone.0217947
Winkler, D., and K. B. Harris. 2009. Reference document: Antmicrobial interventons for beef. Beef Industry Food Safety Council. https://www.bifsco.org/Media/BIFSCO/Docs/antimicrobial-interventions-for-beef.pdf. (Accessed 16 June 2023).
Wythe, L. A., D. K. Dittoe, K. M. Feye, E. G. Olson, L. M. Perry, and S. C. Ricke. 2022. Reduction of Salmonella Infantis on skin-on, bone-in chicken thighs by cetylpyridinium chloride application and the impact on the skin microbiota. Poultry Sci. 101:101409. doi: https://doi.org/10.1016/j.psj.2021.101409
Yang, X., H. Wang, S. Hrycauk, and M. D. Klassen. 2021. Effects of peroxyacetic acid spray and storage temperatures on the microbiota and sensory properties of vacuum-packed subprimal cuts of meat. Appl. Environ. Microb. 87:e03143–20. doi: https://doi.org/10.1128/AEM.03143-20
Yilmaz, P., L. W. Parfrey, P. Yarza, J. Gerken, E. Pruesse, C. Quast, T. Schweer, J. Peplies, W. Ludwig, and F. O. Glöckner. 2014. The SILVA and “All-species Living Tree Project (LTP)” taxonomic frameworks. Nucleic Acids Res. 42:D643–D648. doi: https://doi.org/10.1093/nar/gkt1209
Zhang, L., A. Morey, S. F. Bilgili, S. R. McKee, and L. J. Garner. 2019. Effectiveness of several antimicrobials and the effect of contact time in reducing Salmonella and Campylobacter on poultry drumsticks. J. Appl. Poultry Res. 28:1143–1149. doi: https://doi.org/10.3382/japr/pfz080
Supplemental Material
Shannon Entropy2 | Pielou’s Evenness3 | ||||||
---|---|---|---|---|---|---|---|
Group 1 | Group 2 | H4 | p-value | q-value | H4 | p-value | q-value |
0 ppm (n = 5) | 200 ppm (n = 4) | 0.54 | 0.46 | 0.78 | 0.54 | 0.46 | 1.00 |
400 ppm (n = 4) | 0.00 | 1.00 | 1.00 | 0.24 | 0.62 | 1.00 | |
800 ppm (n = 5) | 0.27 | 0.60 | 0.78 | 0.10 | 0.75 | 1.00 | |
NT (n = 5) | 0.01 | 0.92 | 1.00 | 0.01 | 0.92 | 1.00 | |
200 ppm (n = 4) | 400 ppm (n = 4) | 0.33 | 0.56 | 0.78 | 0.00 | 1.00 | 1.00 |
800 ppm (n = 5) | 2.16 | 0.14 | 0.78 | 0.06 | 0.81 | 1.00 | |
NT (n = 5) | 0.24 | 0.62 | 0.78 | 0.24 | 0.62 | 1.00 | |
400 ppm (n = 4) | 800 ppm (n = 5) | 0.54 | 0.46 | 0.78 | 0.06 | 0.81 | 1.00 |
NT (n = 5) | 0.24 | 0.62 | 0.78 | 0.96 | 0.33 | 1.00 | |
800 ppm (n = 5) | NT (n = 5) | 1.32 | 0.25 | 0.78 | 0.88 | 0.35 | 1.00 |
No significance effect of treatment (P > 0.05, Q > 0.05).
Shannon Entropy - α-diversity metric that quantitatively measures community richness.
Pielou’s Evenness - α-diversity metric that measures community evenness.
H value: test statistic for Kruskal-Wallis.
Bray-Curtis3 | Weighted Unifrac4 | ||||||
---|---|---|---|---|---|---|---|
Group 1 | Group 2 | R5 | p-value | q-value | R5 | p-value | q-value |
0 ppm (n = 5) | 200 ppm (n = 4) | −0.03 | 0.53 | 0.80 | −0.06 | 0.65 | 0.84 |
400 ppm (n = 4) | −0.09 | 0.69 | 0.80 | −0.06 | 0.52 | 0.84 | |
800 ppm (n = 5) | 0.00 | 0.38 | 0.80 | −0.02 | 0.47 | 0.84 | |
NT (n = 5) | −0.08 | 0.72 | 0.80 | −0.05 | 0.58 | 0.84 | |
200 ppm (n = 4) | 400 ppm (n = 4) | −0.05 | 0.64 | 0.80 | −0.05 | 0.56 | 0.84 |
800 ppm (n = 5) | 0.05 | 0.31 | 0.80 | 0.05 | 0.29 | 0.84 | |
NT (n = 5) | −0.03 | 0.60 | 0.80 | −0.06 | 0.68 | 0.84 | |
400 ppm (n = 4) | 800 ppm (n = 5) | −0.11 | 0.82 | 0.82 | −0.12 | 0.84 | 0.84 |
NT (n = 5) | −0.07 | 0.69 | 0.80 | −0.09 | 0.81 | 0.84 | |
800 ppm (n = 5) | NT (n = 5) | 0.04 | 0.32 | 0.80 | 0.02 | 0.38 | 0.84 |
No significance effect of treatment (P > 0.05, Q > 0.05).
ANOSIM – analysis of similarities is a non-parametric statistical test used to determine whether two or more groups are significantly different.
Bray-Curtis - β-diversity metric that quantitatively measures community dissimilarities without the incorporation of phylogenetic relationships.
Weighted Unifrac - β-diversity metric that quantitatively measures community dissimilarities by incorporating phylogenetic relationships.
R value: Test statistic for ANOSIM.
Shannon Entropy2 | Pielou’s Evenness3 | ||||||
---|---|---|---|---|---|---|---|
Group 1 | Group 2 | H4 | p-value | q-value | H4 | p-value | q-value |
0 ppm (n = 5) | 200 ppm (n = 5) | 0.10 | 0.75 | 0.84 | 0.01 | 0.92 | 1.00 |
400 ppm (n = 5) | 1.84 | 0.17 | 0.44 | 4.81 | 0.03 | 0.28 | |
800 ppm (n = 4) | 1.50 | 0.22 | 0.44 | 0.06 | 0.81 | 1.00 | |
NT (n = 5) | 1.84 | 0.17 | 0.44 | 1.84 | 0.17 | 0.55 | |
200 ppm (n = 5) | 400 ppm (n = 5) | 0.88 | 0.35 | 0.58 | 0.88 | 0.35 | 0.66 |
800 ppm (n = 4) | 0.54 | 0.46 | 0.66 | 0.00 | 1.00 | 1.00 | |
NT (n = 5) | 0.01 | 0.92 | 0.92 | 0.53 | 0.46 | 0.66 | |
400 ppm (n = 5) | 800 ppm (n = 4) | 0.24 | 0.62 | 0.78 | 2.16 | 0.14 | 0.55 |
NT (n = 5) | 3.94 | 0.05 | 0.44 | 0.53 | 0.46 | 0.66 | |
800 ppm (n = 4) | NT (n = 5) | 1.50 | 0.22 | 0.44 | 1.50 | 0.22 | 0.55 |
Bolded values denote significance (P ≤ 0.05, Q ≤ 0.05).
Shannon Entropy - α-diversity metric that quantitatively measures community richness.
Pielou’s Evenness - α-diversity metric that measures community evenness.
H value: test statistic for Kruskal-Wallis.
Bray-Curtis3 | Weighted Unifrac4 | ||||||
---|---|---|---|---|---|---|---|
Group 1 | Group 2 | R5 | p-value | q-value | R5 | p-value | q-value |
0 ppm (n = 5) | 200 ppm (n = 5) | 0.17 | 0.05 | 0.12 | 0.11 | 0.15 | 0.51 |
400 ppm (n = 5) | 0.60 | 0.01 | 0.04 | 0.07 | 0.23 | 0.51 | |
800 ppm (n = 4) | 0.55 | 0.01 | 0.04 | 0.34 | 0.03 | 0.31 | |
NT (n = 5) | −0.05 | 0.77 | 0.77 | −0.03 | 0.59 | 0.65 | |
200 ppm (n = 5) | 400 ppm (n = 5) | 0.07 | 0.26 | 0.40 | 0.02 | 0.38 | 0.51 |
800 ppm (n = 4) | −0.04 | 0.43 | 0.47 | −0.13 | 0.82 | 0.82 | |
NT (n = 5) | 0.02 | 0.28 | 0.40 | 0.02 | 0.30 | 0.51 | |
400 ppm (n = 5) | 800 ppm (n = 4) | 0.00 | 0.41 | 0.47 | 0.06 | 0.25 | 0.51 |
NT (n = 5) | 0.22 | 0.05 | 0.12 | 0.01 | 0.41 | 0.51 | |
800 ppm (n = 4) | NT (n = 5) | 0.11 | 0.21 | 0.40 | 0.02 | 0.35 | 0.51 |
Bolded values denote significance (P ≤ 0.05, Q ≤ 0.05).
ANOSIM – analysis of similarities is a non-parametric statistical test used to determine whether two or more groups are significantly different.
Bray-Curtis - β-diversity metric that quantitatively measures community dissimilarities without the incorporation of phylogenetic relationships.
Weighted Unifrac - β-diversity metric that quantitatively measures community dissimilarities by incorporating phylogenetic relationships.
R value: Test statistic for ANOSIM.