CoronaVac immunogenicity decreases with time, and we aimed to investigate whether gut microbiota associate with longer-term immunogenicity of CoronaVac. This was a prospective cohort study recruiting two-dose CoronaVac recipients from three centres in Hong Kong. We collected blood samples at baseline and day 180 after the first dose and used chemiluminescence immunoassay to test for neutralizing antibodies (NAbs) against the receptor-binding domain (RBD) of wild-type SARS-CoV-2 virus. We performed shotgun metagenomic sequencing performed on baseline stool samples. The primary outcome was the NAb seroconversion rate (seropositivity defined as NAb ≥ 15AU/mL) at day 180. Linear discriminant analysis [LDA] effect size analysis was used to identify putative bacterial species and metabolic pathways. A univariate logistic regression model was used to derive the odds ratio (OR) of seropositivity with bacterial species. Of 119 CoronaVac recipients (median age: 53.4 years [IQR: 47.8-61.3]; male: 39 [32.8%]), only 8 (6.7%) remained seropositive at 6 months after vaccination. Bacteroides uniformis (log10LDA score = 4.39) and Bacteroides eggerthii (log10LDA score = 3.89) were significantly enriched in seropositive than seronegative participants. Seropositivity was associated with B. eggerthii (OR: 5.73; 95% CI: 1.32-29.55; p = 0.022) and B. uniformis with borderline significance (OR: 3.27; 95% CI: 0.73-14.72; p = 0.110). Additionally, B. uniformis was positively correlated with most enriched metabolic pathways in seropositive vaccinees, including the superpathway of adenosine nucleotide de novo biosynthesis I (log10LDA score = 2.88) and II (log10LDA score = 2.91), as well as pathways related to vitamin B biosynthesis, all of which are known to promote immune functions. In conclusion, certain gut bacterial species (B. eggerthii and B. uniformis) and metabolic pathways were associated with longer-term CoronaVac immunogenicity.
Keywords: COVID-19 vaccine; CoronaVac; gut microbiota; vaccine; vaccine immunogenicity.