Statistical analysis was carried out in SAS version 9.2 (Cary, NC, USA). Residuals of biomarker concentrations and spore densities were normally distributed, indicating that parametric tests were appropriate. ANOVA of biomarker concentrations and spore densities were carried out using a general linear model, while percent colonization was analyzed using a generalized linear model employing a link-logit function (Stroup, 2014). Pearson correlation coefficients were calculated in SAS