They did correct for multiple testing - "P-values were multiple-test corrected using the Benjamini-Hochberg False Discovery Rate (FDR) procedure and a cutoff of

They did correct for multiple testing - "P-values were multiple-test corrected using the Benjamini-Hochberg False Discovery Rate (FDR) procedure and a cutoff of 0.1 was applied to select the differentially abundant genera". However, they used Wilcoxon test which has unacceptably high FDR [1]. [1] https://www.nature.com/articles/s41467-022-28034-z https://www.nature.com/articles/s41467-022-28034-z