Recruitment of lactating mothers receiving BNT162b2
Between June and Dec 2021, 49 lactating mothers planning on receiving BNT162b2 were recruited and followed up at four timepoints (Fig. 1A). After removing samples without available sequencing data and one outlier (positive ELISA result pre-vaccination), 175 samples from 44 participants remained for downstream analyses. Participants’ ages ranged from 25 to 42 years, with a median age of 36 years (interquartile range (IQR): 31, 39). Nearly half (44.2%) of the participants gave birth via planned cesarean section. The intervals between two vaccine doses ranged between 20–42 days (median: 21 days, IQR: 21, 27) (Table 1).
Table 1
Baseline characteristics of the study population in different immunity response group at one week after second dose
Variables | Overall (N = 43) | Low-IgA subjects (N = 25) | High-IgA subjects (N = 18) | P value |
Maternal age (years) | | | | 0.250 |
Mean (SD) | 35.0 (4.50) | 35.6 (4.20) | 34.2 (4.89) | |
Median (IQR) | 36.0 (31.0,39.0) | 37.0 (32.0,39.0) | 32.0 (30.3,37.8) | |
Gestational age (weeks) | | | | 0.960 |
Mean (SD) | 38.9 (1.03) | 38.9 (1.05) | 38.8 (1.03) | |
Median (IQR) | 38.9 (38.0, 39.7) | 38.9 (38.0, 39.7) | 38.9 (38.0, 39.8) | |
Infant's birth weight (kg) | | | | 0.658 |
Mean (SD) | 3.3 (0.41) | 3.3 (0.48) | 3.2 (0.29) | |
Median (IQR) | 3.2 (3.1, 3.6) | 3.2 [3.0, 3.6] | 3.2 [3.0, 3.6] | |
Total number of previous children | | | | 0.213 |
0 | 3 (7.0) | 3 (12.0) | 0 (0) | |
1 | 20 (46.5) | 12 (48.0) | 8 (44.4) | |
2 | 15 (34.9) | 6 (24.0) | 9 (50.0) | |
3 | 5 (11.6) | 4 (16.0) | 1 (5.6) | |
Delivery type | | | | 0.414 |
Spontaneous vaginal delivery | 13 (30.2) | 5 (20.0) | 8 (44.4) | |
Assisted vaginal delivery | 3 (7.0) | 2 (8.0) | 1 (5.6) | |
Planned C-section | 19 (44.2) | 13 (52.0) | 6 (33.3) | |
Emergency C- section | 8 (18.6) | 5 (20.0) | 3 (16.7) | |
Labor induced | | | | 1 |
No | 28 (65.1) | 16 (64.0) | 12 (66.7) | |
Yes | 15 (34.9) | 9 (36.0) | 6 (33.3) | |
Epidural anesthesia use | | | | 0.057 |
No | 25 (58.1) | 11 (44.0) | 14 (77.8) | |
Yes | 17 (39.5) | 13 (52.0) | 4 (22.2) | |
Missing | 1 (2.3) | 1 (4.0) | 0 (0) | |
Intramuscular analgesia use | | | | 0.516 |
No | 31 (72.1) | 19 (76.0) | 12 (66.7) | |
Yes | 12 (27.9) | 6 (24.0) | 6 (33.3) | |
Infant's sex | | | | 0.765 |
Female | 25 (58.1) | 14 (56.0) | 11 (61.1) | |
Male | 18 (41.9) | 11 (44.0) | 7 (38.9) | |
Smoked previously | | | | 0.419 |
No | 42 (97.7) | 25 (100) | 17 (94.4) | |
Yes | 1 (2.3) | 0 (0) | 1 (5.6) | |
Partner smokes | | | | 1 |
No | 38 (88.4) | 22 (88.0) | 16 (88.9) | |
Yes | 5 (11.6) | 3 (12.0) | 2 (11.1) | |
Marital status | | | | 0.419 |
Married | 42 (97.7) | 25 (100) | 17 (94.4) | |
Not married | 1 (2.3) | 0 (0) | 1 (5.6) | |
Place of birth | | | | 0.766 |
Hong Kong SAR | 33 (76.7) | 19 (76.0) | 14 (77.8) | |
Mainland, China | 6 (14.0) | 3 (12.0) | 3 (16.7) | |
Others | 4 (9.3) | 3 (12.0) | 1 (5.6) | |
Years living in HK (customized) | | | | 0.455 |
< 15 years | 9 (20.9) | 4 (16.0) | 5 (27.8) | |
>=15 years | 34 (79.1) | 21 (84.0) | 13 (72.2) | |
Educational level (customized) | | | | 0.47 |
Postgraduate degree or above | 14 (32.6) | 9 (36.0) | 5 (27.8) | |
Below university level | 8 (18.6) | 3 (12.0) | 5 (27.8) | |
University degree | 21 (48.8) | 13 (52.0) | 8 (44.4) | |
Family income per month (customized) | | | | 0.144 |
HK$40,000 or more | 38 (88.4) | 24 (96.0) | 14 (77.8) | |
Under HK$40,000 | 5 (11.6) | 1 (4.0) | 4 (22.2) | |
Categorical data are presented as number (percentage) and continuous variables as median (interquartile range). Within-group valid percentages are shown. Only subjects with both available antibody and 16S data (at one week after second dose) were summarized in the table. |
Levels of anti-SARS-CoV-2 IgA and IgG in breast milk were boosted after two vaccine doses
After the first dose, levels of SARS-CoV-2 spike protein-specific IgA and IgG in breast milk remained unchanged compared with pre-vaccination. Levels were boosted one week post-second dose (median (IQR), IgA OD: 0.116 (0.096, 0.137) vs. 0.198 (0.164, 0.280), Fig. 1B; IgG OD: 0.058 (0.054, 0.062) vs. 0.216 (0.132, 0.310), p < 0.001, Fig. 1C). Although IgA levels returned to baseline levels one month post-second dose (p < 0.001), IgG levels remained significantly higher compared to baseline (IgG OD: 0.058 vs. 0.164, p < 0.001). Additionally, there was a significant moderate association between the two antibodies in breast milk (Spearman’s rho coefficient: 0.46, p < 0.001 Supplementary Fig. 1A).
We also found that mothers who received epidural anesthesia during delivery had a lower IgA levels one week post-second dose (p = 0.020, Supplementary Table 1).
Changes in breast milk bacterial richness and composition were observed post-vaccination
We conducted 16S rRNA amplicon sequencing analysis to determine if microbiome composition was associated with levels of vaccine-induced anti-RBD antibodies in breast milk. Breast milk microbial composition (Supplementary Fig. 2A) and diversity (alpha and beta diversities) (Fig. 1D-H) were compared over different timepoints. We observed a significant increase in bacterial Chao1 richness between pre-vaccine and one week post-first dose samples (p = 0.044), but not in terms of Shannon diversity (p = 0.440, Fig. 1D-E). Microbiome richness returned to baseline levels one week post-second dose. To determine whether this change in alpha diversity is associated with vaccination, we analyzed longitudinal breast milk microbiota data from an independent, unvaccinated cohort22, in which we observed a continuous decline of Chao1 diversity (Supplementary Fig. 1B). Regarding beta-diversity, a significant shift was observed between pre-vaccination and all post-second dose vaccination timepoints (adjusted p = 0.030, 0.006, and 0.006 for one week post-first dose and one week and one month post-second dose, respectively, Fig. 1F). However, no significant changes in microbiome beta diversity were observed in the independent cohort (Supplementary Fig. 1C; adjusted p = 0.240 0.360, 1.000, respectively). We also found that baseline breast milk microbiome composition was significantly associated with intramuscular analgesia use during labor (p = 0.012, Supplementary Table 2).
We then conducted LEfSe to identify differentially abundant bacterial taxa between timepoints. In total, 109 differentially abundant species were identified between pre- and post-vaccination timepoints, ten of which remained differentially abundant throughout all three post-vaccination timepoints. Fifty-three species were differentially abundant between baseline levels and post-second dose, including Anaerococcus octavius, Arthrobacter russicus, Bacteroides caecimuris, Clostridiaceae bacterium, Helicobacter rodentium, Lactobacillus aviarius, and Rothia sp. (LEfSe, LDA score > 1.5 and p value < 0.05, Fig. 1G, Supplementary Table 3). The 12 most abundant species varied across samples and timepoints (Supplementary Fig. 2).
Baseline breast milk bacterial composition is associated with anti-SARS-CoV-2 antibody levels in breast milk following vaccination
We then investigated whether breast milk bacterial composition was associated with its antibody levels post-vaccination. Since immune parameters peaked at one week post-second dose (Fig. 1B-C), we focused our subsequent analyses on this timepoint, thereby dichotomizing the cohort into those with high- and low-IgA levels using the prior-specified cutoff (High, 20; Low, 28; Fig. 1A).
The Chao1 richness and beta diversity of baseline microbiota were not significantly different between high- and low-IgA participants (p = 0.464 and 0.071, respectively, Fig. 2A). However, the baseline Shannon richness in high-IgA subjects was significantly higher than that of low-IgA subjects (p = 0.046, Fig. 2B).
We then investigated whether there were taxonomic differences in baseline microbiota between high- and low-IgA subjects. For both groups, the most abundant phylum was Firmicutes, followed by Proteobacteria and Actinobacteria. The relative abundance of Firmicutes was higher in low-IgA subjects (median relative abundance of 61.78% and 74.62% for high- and low-IgA group, respectively), while levels of Proteobacteria were depleted (18.64% and 9.01% for high- and low-IgA group, respectively, Supplementary Fig. 2A).
We identified 23 differentially abundant species between the high- and low-IgA subjects (p = 0.004, Fig. 2D). Unclassified Neisseria and Corynebacterium kroppenstedti showed the largest effect size in enrichment for high- and low-IgA subjects, respectively (LDA score: 3.69 and 4.12, respectively). The mixed effects model showed that unclassified Neisseria and Neisseria elongata were persistently higher in subjects with high breast milk IgA levels (p = 0.045 and 0.031, respectively, Supplementary Table 4).
Besides differences in baseline microbiota, taxonomic differences in microbiota one week post-second dose between high- and low-IgA subjects were identified by LEfSe. Higher abundances of Bifidobacterium bifidum and unclassified Bifidobacterium were found in the high-IgA subjects (LDA score: 2.59 and 2.87, respectively, Fig. 3).
We then used random forest to investigate the power of these distinct microbial species in predicting levels of anti-SARS-CoV-2 antibodies in breast milk one week post-second dose. We constructed a model based on five microbial features including Neisseria elongata, Pseudomonas balearica, unclassified Pseudopropionibacterium, unclassified Neisseria, and unclassified Aggregatibacter. The model based solely on microbial features [ AUC (95% CI): 0.72 (0.58–0.85)] performed better than the model which combined microbial and demographic predictors [AUC (95% CI): 0.66 (0.53–0.80), Fig. 4A].
To investigate the potential mechanisms by which pre-vaccination breast milk microbiota affects antibody levels, we estimated the association between predicted microbial functional pathways and anti-SARS-CoV-2 IgA levels at one week post-second dose. A total of 437 pathways were predicted by PICRUSt2, among which nine were significantly enriched in high-IgA subjects (Wilcoxon rank-sum, adjusted p < 0.05, Fig. 4). The correlation network showed that higher abundances of selected bacterial markers were associated with sugar metabolism pathways (e.g., sucrose degradation IV, glucose and glucose-1-phosphate degradation, fucose degradation). For example, fucose degradation was positively correlated with all bacterial markers. Additionally, superpathways for arginine and polyamine biosynthesis were correlated with all markers expect for Neisseria elongate (Fig. 4B).
Impact of vaccination on the abundance of Lactobacilli and Bifidobacteria in breast milk
Lactobacilli and Bifidobacteria constitute a core microbiota of shared taxa between breast milk and infant stool24,25. Given those genera play a beneficial role in the healthy development of the infant’s gut microbiota and immune system26,27, we examined the association between these taxa and immune responses to BNT162b2 in breast milk. We detected the presence of Lactobacilli (L. aviarius, L. fermentum, L. gasseri, L. iners, L. intestinalis, and unclassified Lactobacillus) and Bifidobacteria (B. bifidum, B. breve, B. dentium, B. longum, and unclassified Bifidobacterium). Vaccination did not significantly change the abundance of Lactobacillus and Bifidobacterium in breast milk, except for L. aviarius and L. fermentum (Fig. 5, Supplementary Table 5). However, these two genera were not persistently high in subjects with high IgA levels (Supplementary Table 6).