Fig. 4: Baseline breast milk bacterial species and functions associated with high- and low-IgA participants to BNT162b2 vaccination at one week post-second vaccine dose (N = 43; High IgA: 18, Low IgA: 25).

a AUROC (95% CI) values of models based on metagenomic markers, demographical markers (maternal age, different time interval between two doses of vaccine, and the use of epidural anesthesia), and combined markers using a random forest classification approach. AUROC, area under the receiver operating characteristic curve. The performance of random forest models was compared using the bootstrap method. b Baseline bacterial species and pathways associated with levels of immunoglobulin A. Differential baseline breast milk bacterial species and pathways were detected by a random forest model and Wilcoxon’s rank-sum test (two-sided), respectively. P values were given by Spearman’s correlation test.