Fig. 4

Boxplot of the ResFinder ARG phenotypes that showed the strongest case associations in the three CRC studies (excluding adenoma samples) with strong overall resistome expansion in cases (ratio > 1.2). Datasets (columns) are sorted from left to right by the strongest to weakest total ARG abundance case association. ResFinder ARG phenotypes (rows) are sorted from strongest (top) to weaker (bottom) case association in RDA analysis constrained by case-control status with study as covariate. p values are calculated by Wilcoxon rank-sum test and FDR corrected within each study for the number of ResFinder phenotype categories compared