- 1Department of Microbiology, The Key Laboratory of Microbiology and Parasitology of Anhui Province, The Key Laboratory of Zoonoses of High Institutions in Anhui, School of Basic Medical Sciences, Anhui Medical University, Hefei, China
- 2Department of Clinical Medicine, The First School of Clinical Medicine, Anhui Medical University, Hefei, China
- 3Department of Nosocomial Infection Control, Anhui No.2 Provincial People’s Hospital, Hefei, China
- 4School of Life Sciences, Anhui Medical University, Hefei, China
- 5School of Nursing, Anhui Medical University, Hefei, China
SARS-CoV-2 and its mutant strains continue to rapidly spread with high infection and fatality. Large-scale SARS-CoV-2 vaccination provides an important guarantee for effective resistance to existing or mutated SARS-CoV-2 virus infection. However, whether the host metabolite levels respond to SARS-CoV-2 vaccine-influenced host immunity remains unclear. To help delineate the serum metabolome profile of SARS-CoV-2 vaccinated volunteers and determine that the metabolites tightly respond to host immune antibodies and cytokines, in this study, a total of 59 sera samples were collected from 30 individuals before SARS-CoV-2 vaccination and from 29 COVID-19 vaccines 2 weeks after the two-dose vaccination. Next, untargeted metabolomics was performed and a distinct metabolic composition was revealed between the pre-vaccination (VB) group and two-dose vaccination (SV) group by partial least squares-discriminant and principal component analyses. Based on the criteria: FDR < 0.05, absolute log2 fold change greater than 0.25, and VIP >1, we found that L-glutamic acid, gamma-aminobutyric acid (GABA), succinic acid, and taurine showed increasing trends from SV to VB. Furthermore, SV-associated metabolites were mainly annotated to butanoate metabolism and glutamate metabolism pathways. Moreover, two metabolite biomarkers classified SV from VB individuals with an area under the curve (AUC) of 0.96. Correlation analysis identified a positive association between four metabolites enriched in glutamate metabolism and serum antibodies in relation to IgG, IgM, and IgA. These results suggest that the contents of gamma-aminobutyric acid and indole in serum could be applied as biomarkers in distinguishing vaccinated volunteers from the unvaccinated. What’s more, metabolites such as GABA and taurine may serve as a metabolic target for adjuvant vaccines to boost the ability of the individuals to improve immunity.
Introduction
As of now, SARS-CoV-2 is still ravaging the world. Vaccination is treated as an integral part of the long-term management of the COVID-19 epidemic (1, 2), which elicits protective neutralizing antibodies against SARS-CoV-2 and offers hope for containing the COVID-19 pandemic. As of 26 April 2022, more than 11.32 billion doses of the vaccine have been administered worldwide (3), which have been reported with substantial efficacy against SARS-CoV-2 and its variants (4–6). Recent newly emerging variants, combined with the steady decline in antibody levels among vaccinated individuals, have implied a resurgence of the epidemic over time (7–9). The Sinovac-CoronaVac COVID-19 vaccine is licensed as a vaccine for the prevention of COVID-19 and against emerging SARS-CoV-2 coronavirus variants in China. Despite a large-scale population-based cohort safety evaluation of SARS-CoV-2 vaccines, these prior studies have disregarded the underlying function of biologically active metabolic small molecules between different immune cells in the organism in triggering the immune response associated with vaccination. The high throughput-based omics techniques to study systemic vaccinology at the cellular and molecular levels allow for a comprehensive overview of the interrelationship between vaccination and vaccine-induced immune responses.
There have been a handful of proteomic- and/or metabolomics-based studies since the COVID-19 outbreak that revealed a previously unrecognized metabolic profile associated with SARS-CoV-2 infection and disease severity (10–12). In general, liquid chromatography–mass spectrometry (LC/MS)-based metabolomics studies identified the dysregulation of metabolites such as butyric acid, 2-hydroxybutyric acid, l-phenylalanine, and kynurenine metabolites in the serum of COVID-19 patients (13, 14). However, there have been few prospective studies involving metabolomics-based studies of SARS-CoV-2 vaccination-induced immune responses in the human population. To date, there have been only two studies based on metabolomics to investigate the alteration of serum metabolites in the host metabolism after the inoculation of the COVID-19 vaccine and the network interplay with the immune response, and they mainly identified metabolites in the processes of tricarboxylic acid cycle, amino acid metabolism, and lipid metabolism associated with the immune response (15, 16). Nevertheless, there is a paucity of studies on the interaction between circulating metabolites and vaccine-induced antibody responses after COVID-19 vaccination.
Metabolism is the basic characteristic and main activity of life. LC/MS-based metabolomics is widely used in the mechanistic study and diagnosis of various diseases, such as cancer, diabetes, and cardiovascular or pulmonary disorders upon infection of a virus (17, 18). A growing number of studies have underscored the impact that the metabolic state of an organism can exert on the shaping of the body’s immune system (19–21). For example, a preceding study revealed a confidential link between metabolic phenotype and vaccine immunity in healthy individuals (22). In addition, findings from clinical or animal studies suggest that the metabolites derived from gut microbiota are critical in regulating the immune responses to vaccination. Various gut bacterial metabolites have been reported to be implicated in modulating the response of vaccine immunity, such as short-chain fatty acids (SCFAs), bile acids, and aromatic amino acids (AAAs). Short-chain fatty acids are the best documented to have many beneficial benefits to the individual, ranging from but not limited to the regulation of glucolipid metabolism, inflammation, immunity, and tumors (23). A previous work has noted that SCFAs could potentiate the expression of genes involved in plasma cell differentiation and class transformation, which provides the energy requirements for the antibody response in the non-pathogenic infection state and in the presence of pathogenic infection (24).
To further elucidate and provide a comprehensive insight into whether COVID-19 vaccines influence the response of host serum metabolites and the interplay between metabolites and vaccine-stimulated antibodies and cytokines, in this study, a metabolomic analysis has been applied to screen for metabolic variations and to provide a comprehensive view of endogenous metabolites in relation to the vaccination. We conducted a prospective observational study of undergraduates who received the inactivated SARS-CoV-2 vaccine (CoronaVac; Sinovac) to examine the serum metabolic determinants of immunization responses.
Volunteers and methods
Study design and participants
The study recruited a total of 30 healthy undergraduates as participants from Anhui Medical University. We collected the sera of those volunteers before vaccination in 9 May 2021 and after receiving two doses of the inactivated COVID-19 vaccine (CoronaVac, Sinovac Life Sciences, Beijing, China; N = 30) in 21 June 2021 at the First Affiliated Hospital of Anhui Medical University. Eligible participants recruited from Anhui Medical University were aged 18–20 years old, had a well-regulated lifestyle, had eaten in the university canteen for the past 2 months, and had no history of SARS-CoV-2 infection upon inoculating the CoronaVac vaccine. The exclusion criteria were volunteers who had a history of allergies, had a history of medication, and received antibiotic drugs in the last 2 months; these were determined by means of a questionnaire. All participants provided written informed consent and completed two doses of the vaccine. All individuals enrolled in this study were asked to sign an informed consent before the investigation, following the principle of informed consent and strict confidentiality, which has been approved by the ethics committee of Anhui Medical University Biomedical Ethics Committee (Approval No. 2021H021). Each individual or his/her legal guardian provided the signed informed consent before enrolment. All volunteer data sets were anonymized.
Measurement of indicators using routine blood tests
Various indicators obtained from routine blood tests, such as eosinophils, neutrophils, and hemoglobin, were measured via using routine blood tests conducted at Hefei City Maternal and Child Health and Family Planning Service Center (Anhui Province). Specifically, a 2-ml blood sample was collected in an EDTA-K2 anticoagulant tube. Immediately after collection, the tube was mixed by being gently reversed several times to ensure adequate anticoagulation of the blood specimen. The blood samples were then sent to Hefei City Maternal and Child Health and Family Planning Service Center, and 18 indicators, such as white blood cell (WBC), neutrophil (Neu), monocytes (Mon), red blood cell (RBC), and hematocrit value (HCT), obtained from routine blood tests, were directly measured using a three-classification blood cell analyzer (BC-5390 CRP, Marry). The results from the routine blood tests of all blood samples are summarized in Table S1.
Serum sample collection and pretreatment
The elbow vein blood (5 ml) was collected from volunteers before vaccination and after receiving two doses of the vaccine at 6 a.m. to 8 a.m. in a fasting state with vacuum negative-pressure blood collection vessels, and was immediately centrifuged for 15 min (3,500 rpm, 4°C). Each aliquot (0.5 ml) of the serum samples was prepared and stored at –80°C until use (25). Then, all samples were thawed on ice before the metabolite profiling (26). Untargeted metabolomics of sera samples was performed as previously described (27), with slight modification. In brief, 80 µl of sample was mixed with 320 μl of acetonitrile by vortexing for 60 s. Then, the sample was centrifuged at 18,920×g for 10 min (4°C) to precipitate the protein. Two 150-μl aliquots of the supernatant were transferred and lyophilized for analysis in positive and negative electrospray ionization (ESI+ and ESI−) mode, respectively. Fifty microliters of a 25% (by volume) acetonitrile aqueous solution was used to reconstitute the sample before the LC-MS analysis. Quality control (QC) samples were prepared by pooling the same volume of each sample to evaluate the reproducibility of the analysis. The pretreatment of the QC samples paralleled and was the same as that of the study samples. The QC samples were evenly inserted in each set of the analysis running sequence to monitor the stability of the large-scale analysis (28).
Untargeted LC-MS/MS analysis
The metabolite profiles of serum were detected by using an LC-MS/MS system consisting of an ultra-high-performance liquid chromatography (UPLC) system (Vanquish, Thermo Fisher Scientific, Waltham, MA, USA) with a UPLC BEH Amide column (2.1 mm × 100 mm, 1.7 μm) coupled to a Q Exactive (QE) HFX mass spectrometer (Orbitrap MS, Thermo Fisher Scientific, Waltham, MA, USA). The UPLC system was used for the separation of biochemicals via the UPLC BEH Amide column. The mobile phase consisted of 25 mmol/L ammonium acetate (Sigma-Aldrich, Darmstadt, Germany) and 25 mmol/L ammonia hydroxide (Fisher Chemical, Waltham, MA, USA) in water (pH = 9.75) (A phase) and acetonitrile (B phase). The autosampler temperature was 4°C, and the injection volume was 3 μl. The QE HFX mass spectrometer was applied to obtain MS/MS spectra in information-dependent acquisition (IDA) mode under the control of the acquisition software (Xcalibur, Thermo Fisher Scientific, Waltham, MA, USA). In this mode, the acquisition software continuously evaluates the full-scan MS spectrum. The electrospray ion (ESI) conditions were as follows: sheath gas flow rate of 30 Arb, auxiliary gas flow rate of 25 Arb, capillary temperature of 350°C, full MS resolution of 60,000, MS/MS resolution of 7500, collision energy of 10/30/60 in negative ion mode (NCE), and spray voltage of 3.6 kV (positive) or −3.2 kV (negative).
Serum metabolomic analysis
The original LC-MS/MS data were converted to the mzXML format using ProteoWizard and processed with an in-house program, which was developed using R and based on XCMS, for peak detection, extraction, alignment, and integration. Subsequently, for further multivariate analysis, the raw data were pretreated. Pretreatment included de-noising based on the relative standard deviation (RSD), filling the missing data via half of the minimum value, and the batch normalization with peak area was used to calibrate the metabolomics data. The final dataset contained the information of peak number, sample name, and normalized peak area and was imported to the SIMCA16.0.2 software package (Sartorius Stedim Data Analytics AB, Umea, Västerbotten, Sweden). Firstly, we performed a principal component analysis (PCA), an unsupervised analysis, to reduce the dimensions of the data. The 95% confidence interval (95% CI) in the PCA score plot was used as the threshold to identify potential outliers in the dataset. Secondly, in order to visualize the group separation and find significantly changed metabolites, we conducted the supervised orthogonal projections to latent structure discriminate analysis (OPLS-DA) and acquired the value of variable importance in projection (VIP). A sevenfold cross-validation test was performed to evaluate the goodness of fit of the OPLS-DA model using the values of R2Y and Q2. R2Y indicates how well the variation of a variable is explained, and Q2 indicates how well a variable can be predicted. A 200-time permutation test was then conducted to assess the robustness of the model. The precursor molecule passed the combined criteria: 1) absolute log2 fold change (|log2FC|) greater than 0.25, 2) VIP >1, and 3) adjusted p-value <0.05 (nonparametric t-test) were considered significantly changed between the SV and VB groups. Furthermore, after scanning for the differential metabolite features, the fragment information obtained from the MS/MS model was then further matched with the annotations in HMDB, Metlin, massbank, LipidMaps, mzclound, and in-house standards database (PPM <10) to obtain accurate metabolite information. Finally, commercial databases, including Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.genome.jp/kegg/ (accessed on 15 June 2021) and MetaboAnalyst (http://www.metaboanalyst.ca/ (accessed on 15 June 2021), were used for pathway enrichment analysis.
Measurement of cytokines and SARS-CoV-2 antibodies in serum by ELISA
To monitor the response of the host’s body based on cytokines and SARS-CoV-2 antibodies, the levels of inflammatory factors, including IL-2, IL-4, IFN-γ, and SARS-CoV-2 antibodies, including anti-(N+S) IgA, anti-(N+S) IgG, and anti-(N+S) IgM, were measured by ELISA. Specifically, the levels of cytokines and antibodies in whole-blood serum obtained from all of the blood samples (n = 30 and 29 at each sampling point) were measured. We selected IL-2 and IFN-γ as representatives of Th1 cells and IL-4 as the representative of Th2 cells. IL-2, IL-4, IFN-γ, anti-(N+S) IgA, anti-(N+S) IgG, and anti-(N+S) IgM were measured using a high-sensitivity enzyme-linked assay quantitative kit (the kits for cytokines were procured from Bio-Techne USA Co., Ltd., Minnesota, USA, and the kits for antibodies were obtained from Wuhan Fine Biotechnology Co., Ltd., Wuhan, China) in accordance with the manufacturer’s instructions. Specifically, taking the detection of SARS-CoV-2 antibodies as an example, the kit can detect the concentration of SARS-CoV-2 antibodies based on indirect enzyme-linked immunosorbent assay technology. Ninety-six-well plates were percolated with recombinant 2019-nCoV nucleocapsid and spike protein (antigen), and an HRP-conjugated antibody was used as the detection antibody. Subsequently, the standards, test samples, and HRP-conjugated detection antibody were added to the wells, and the plates were washed with wash buffer. 3,3′,5,5′-Tetramethylbenzidine (TMB) substrates were used to visualize the HRP enzymatic reaction, which is catalyzed by HRP to produce a blue product that changes to yellow after the addition of an acidic stop solution. The density of the yellow color is proportional to the target amount of the sample captured in the plate. The O.D. absorbance at 450 nm was read using a microplate reader, and the concentration of the target was then calculated.
Statistical analysis
Differential serum metabolites were identified using the Wilcoxon test and the p-value was adjusted for multiple testing via the false-discovery rate (FDR) using the Benjamini-Hochberg method. An R package Randomforest (version 4.6.14) was used to select important variables in the serum data and build a machine learning model to distinguish whether participants had received the SARS-CoV-2 vaccine. For feature selection, 100 endogenous serum metabolites screened from the 156 distinct metabolites between the VB and SV groups were selected as input features. Finally, the 50 metabolites, whose mean decrease accuracy ranked in the top 50, were determined by 10-fold cross validation using the rfcv function. After selecting 50 metabolites, we adopt binary logistic regression to build the model based on the potential biomarkers. A receiver-operating characteristic curve was used to evaluate the results of the regression analysis. For network analysis, we used Cytoscape (version 3.9.1) (29) to visualize the correlation relationships calculated by the SparCC algorithm (30) among metabolites, while the biologically important subnetworks were identified using the Molecular Complex Detection (MCODE) plugin in the Cytoscape software (31). In addition, a correlation analysis between metabolites and clinical parameters was also performed using Spearman rank correlation analysis.
Apart from data related to metabolomics, other experimental data were analyzed using the GraphPad Prism 9 software (GraphPad Software, La Jolla, CA, USA). The data with a normal distribution were represented as mean ± standard deviation (mean ± SD), and differences between groups were compared using t-test. All statistical tests were two-tailed, and the differences reached the statistical significance which was set at p < 0.05.
Results
SARS-Cov-2 vaccine-induced vigorous antibody responses
We recruited a cohort of college students, and finally, 29 individuals comprising 15 men and 14 women were included in this study for none antibiotic use within 2 months. A serum sample was obtained from each volunteer before vaccination (VB) and 2 weeks after the second vaccination (SV). We analyzed 24 clinical serum measurements of the nonvaccinated and vaccinated subjects (Figure 1A, Table 1), involving six immunoglobulins of IgA, IgG, and IgM; Th1 cytokines of IL-2 and IFN-γ; Th2 cytokines of IL-4; and 18 blood routine indicators. Compared to non-vaccinated individuals, volunteers who received two doses of the SARS-CoV-2 vaccine showed significant upregulation of IgA and IgM, as well as an increase in IgG (Figures 1B–D). However, IL-4 was only found to be significantly different between the VB and SV groups as with cytokines (Figures 1E–G). Furthermore, we found that the levels of five routine blood indicators (RBIs), mean corpuscular volume (MCV), red blood cell distribution width-standard deviation (RDW-SD), mean platelet volume (MPV), platelet-large cell count (P-LCC), and platelet-large platelet ratio (P-LCR), showed a significant decreasing trend in the SV group compared to participants in VB group. However, the levels of other RBIs, such as white blood cells (WBC), red blood cells (RBC), neutrophil (Neu), and eosinophil (Eos) did not reach statistical significance (Figure S1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/954801/fimmu-13-954801-HTML-r1/image_m/fimmu-13-954801-g001.jpg)
Figure 1 Untargeted metabolomic profiling of sera samples collected from volunteers before SARS-CoV-2 vaccination, and 2 weeks after two-dose vaccination. (A) Schematic diagram of the study design. (B–G) Antibodies and cytokines characterized in serum and compared between the VB and SV groups. ns indicates not significant, single asterisk indicates p < 0.05, double asterisks indicate p < 0.01, triple asterisks indicate p < 0.001, quadruple asterisks indicate p < 0.0001.
Metabolomic profiling of sera from volunteers who received COVID-19 vaccines
We performed untargeted metabolomics based on ultra-performance liquid chromatography/tandem mass spectrometry (UPLC-MS/MS) to analyze the sera samples. The typically extracted ion chromatograms (base peak chromatogram) from two ESI modes are displayed in Figure S2. The coefficients of variation of the peak distribution in the QC samples showed that this analysis was stable and repeatable (Figure S3). Additionally, in the principal component analysis score plots of both ESI+ and ESI- (Figures S4A, B), the QC samples clustered tightly together, which further confirmed the reliability of the present study. Furthermore, the discrimination trends among samples between the VB and SV groups revealed obvious separation (Figures S4C, D). After peak alignment, peak picking, and deconvolution, we detected a total of 17,072 and 13,286 precursor molecules in positive and negative ion modes, respectively. We further excluded the metabolic features of which the percentage of relative standard deviation (RSD%) was more than 30% in QC samples, and each peak presented (non-zero value) in more than 50% of total samples was included in the data analysis. Ultimately, 15,485 and 10,411 putative metabolic features in the positive and negative ion modes were obtained and merged together for further analysis, respectively. Altogether, 2376 metabolites were identified and quantified based on searches of online databases (HMDB and Metlin) or confirmed with authentic standards from PANOMIX. These variables were used for the subsequent multivariate and univariate analyses. The details of the metabolites are shown in Table S2. We further checked the influence of sex to the metabolomic profiling. As shown in the UMAP plots (Figure S5), we did not observe sex-biased influence on the metabolomic data.
Apparent alteration of serum metabolites associated with SARS-Cov-2 vaccination
First, the partial least squares discriminant analysis score plot (Figure 2A) revealed clear separations between the VB and SV groups without overfitting (Figure S6). To reveal the serum metabolic characteristics in volunteers who received two doses of the vaccine and identify and confirm high-confidence metabolites associated with vaccination, we distinguished metabolites based on the criteria of a log2 fold change (FC) ≥0.25 or ≤-0.25 and adjusted p value (FDR) <0.05, respectively. Thus, metabolic characteristics with significant differences were extracted and visualized by volcano plots (Figures 2B). In addition, a metabolite with a VIP value greater than 1.0 was selected as a significantly different metabolic feature for further analysis. Finally, 100 out of 156 differential metabolites that were endogenous (56 non-endogenous were discarded) were reserved, containing 60 and 40 metabolites from the ESI+ and ESI- models, respectively. Hierarchical clustering analysis and heatmap also revealed differentially expressed metabolites between VB and SV. These include the enrichment of L-glutamic acid, L-leucine, GABA, taurine, and succinic acid. Interestingly, L-leucine, L-glutamic acid, succinic acid, and citric acid were common differentially expressed metabolites between SARS-CoV-2 infection and immunization, but with the opposite expression (32). Among them, we observed a significant correlation between taurine and antibody levels (Figure S7). In contrast, indole, phenyl acetate, and L-carnitine were depleted in SV individuals compared with the VB group (Figure 2C). Next, we focused on individual serum metabolites for in-depth analyses to delineate potential correlations with vaccination. Random forest (Rf) analyses ranked all altered 100 serum metabolites by contribution to the group separation; the top 50 differentiating metabolites are shown (Figure 2D, Figure S7D). Notably, we noted several metabolites that were statistically different between the VB and SV groups, which were also confirmed as high-impact metabolites with huge importance in Rf analysis. For example, phenyl acetate which ranked sixteenth and L-iditol which ranked eighteenth had 8.6- and 44.1-fold elevations in the VB group.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/954801/fimmu-13-954801-HTML-r1/image_m/fimmu-13-954801-g002.jpg)
Figure 2 Alteration of main metabolites in sera samples from the vaccinated individuals. (A) Partial least squares discriminant analysis score plot based on the VB and SV groups. (B) Volcano plots highlight the serum metabolites that were increased (red) in the SV compared to the VB group, with FDR < 0.05, log2 fold change (FC) >0.25 or <−0.25. (C) Heatmap clustering of distinct serum metabolites from the comparison between the VB and SV groups. (D) Variable importance plot of the top 40 serum metabolites (y-axis) ranked by the contribution to mean decrease accuracy of the Gini coefficient (x-axis) in the random forest model for discerning group difference.
Network analysis and metabolic pathway analysis
To further explore the potential regulatory relationships using the vaccination-related differential metabolites, a correlation network was constructed based on the correlation relationship that the correlation coefficients and p-value should be reached based on the criteria |r| > 0.3 and p < 0.05 (Table S3). On the whole, the interaction network comprised 49 nodes and 117 interactions. We observed that the metabolites enriched in SV group were almost showing positively correlated with each other, such as taurine, succinic acid, glutaric acid, and GABA. However, we also detected that the metabolites depleted in the SV group exhibited a positive association with those that increased in the SV group, for instance, D-mannose, D-galactose, and 2-ketobutyric acid (Figure 3A). In addition, molecular complex detection (MCODE) analysis identified three key subnetworks (module 1, module 2, and module 3) (Figures 3B–D). Module 1 comprised four depleted and eight enriched metabolites in the SV group, such as D-mannose, 2-ketobutyric acid, taurine, and GABA. Nonetheless, they were positively correlated with each other. Module 2 encompassed interactions among two metabolites enriched in VB and three metabolites over-expressed in SV group, while D-galactose in module 3 acted as a link between the other four SV-enriched metabolites and two metabolites that increased in the VB group. Furthermore, KEGG and MetaboAnalyst were comprehensively used to explore the most relevant metabolic pathways. The altered metabolites were enriched or depleted in different metabolomic signaling pathways. With respect to metabolites enriched in SV, 35 metabolic pathways were identified, among which three pathways reached the significance threshold (FDR <0.05), as follows: Butanoate metabolism (including GABA, L-glutamic acid, succinic acid semialdehyde, oxoglutaric acid, and succinic acid), alanine, aspartate, and glutamate biosynthesis (containing succinic acid semialdehyde; L-glutamic acid; GABA; citric acid; succinic acid; oxoglutaric acid), D-glutamine, and D-glutamate metabolism (L-glutamic acid; oxoglutaric acid) (Figure 4A). Taken together, these findings show a sub-metabolic pathway from L-glutamic acid to succinic acid within butanoate metabolism and alanine, aspartate, and glutamate biosynthesis pathways and showed that glutamate decarboxylase, 4-aminobutyrate aminotransferase, and succinate-semialdehyde dehydrogenase are required (Figure 4B). However, the metabolic pathways constructed based on the metabolites enriched in the healthy group did not achieve the threshold of significance (Figure 4C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/954801/fimmu-13-954801-HTML-r1/image_m/fimmu-13-954801-g003.jpg)
Figure 3 Correlation coefficients-based network constructed by distinct metabolites between the VB and SV groups. (A) The nodes were colored and shaped by different groups (ellipses represent SV group, and rectangle stands for the VB group). The color of the edge is set to red and green, which represent positive and negative correlation, respectively. The width of the edge indicates the magnitude of the correlation coefficient. (B–D) The three subnetworks were identified by the MCODE plugin in the Cytoscape software.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/954801/fimmu-13-954801-HTML-r1/image_m/fimmu-13-954801-g004.jpg)
Figure 4 Metabolic KEGG enrichment analysis based on the differential metabolites between groups. (A) KEGG pathway analysis of differentiating metabolites enriched in the SV group. The color of the bubbles represents the value of adjusted p value, and the size of bubbles represents the number of counts (sorted by enrichment ratio). (B) KEGG pathway analysis of distinguishing metabolites increase unique to the VB group. (C) Schematic diagram of four SV-enriched metabolites participating in the butanoate metabolism and the alanine, aspartate, and glutamate metabolism KEGG pathways. Triple asterisks indicate p < 0.001.
Identification of biomarker panel and correlation with clinical serum parameters
To evaluate the reliability of 100 biomarker candidates and define the useful biomarkers, a binary logistic regression analysis and an optimized algorithm of the forward stepwise (Wald) method were employed to construct the best model using 100 differential and potential metabolite biomarkers. For classifying participants in the VB from the SV group, we found that GABA and indole achieved high AUC values of 0.801 and 0.823 (Figures 5A, B), respectively. In addition, the combination of GABA and indole increased AUC to a higher degree of 0.96 in distinguishing individuals in the VB group from individuals without vaccination (Figure 5C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/954801/fimmu-13-954801-HTML-r1/image_m/fimmu-13-954801-g005.jpg)
Figure 5 Metabolite markers for pairwise discrimination of the VB and SV groups. (A-B) Diagnostic accuracy of the single serum metabolite, GABA, and indole in distinguishing individuals that received two-dose vaccination in the SV group from the VB group. (C) The ROC curve with an increased AUC value based on the combination of GABA and indole. (D) Associations between serum metabolites and clinical parameters. The associated metabolites are colored gray, and the associated clinical parameters are colored successively. Each line indicates a significant correlation between a metabolite and a clinical parameter, with light red corresponding to a positive association (correlation coefficient ≥ 0.3) and light blue representing a negative association (correlation coefficient ≤-0.03).
Furthermore, we combined the metabolites enriched in butanoate metabolism or alanine, aspartate and glutamate biosynthesis pathways with biomarker metabolic features and performed global correlation analysis with antibodies or biochemical parameters to identify the potential role of metabolites in mediating host immune response and physiology (antibodies such as IgG, IgM, and IgA, and cytokines such as IL-4 and IFN-γ, and blood RT) upon SARS-CoV-2 vaccination. On the whole, we detected that the serum levels of IgM, IgG, IgA, IL-4, and IFN-γ, blood level of MPV showed a negative correlation with indole intensity, while IL-2 was positively correlated with indole. Notably, GABA, L-glutamic acid, succinic acid, and succinic acid semialdehyde, which were found to be enriched in the SV group, was positively associated with IgG. In addition, IgM was significantly associated with succinic acid, succinic acid semialdehyde, and L-glutamic acid. Moreover, the antibody IgA showed a positive correlation with sera succinic acid and succinic acid semialdehyde. Meanwhile, cytokine IL-4 was also positively associated with succinic acid semialdehyde (Figure 5D, Table S4).
Discussion
In the process of immune activation, the demand for biomass construction and effector molecule synthesis requires a lot of adjustment of cell metabolism, and numerous small molecules are produced. Compared with the components of classical cell signaling pathways mediated by cytoplasm, membrane-bound, or secreted proteins, small metabolites have great evolutionary potential as communication molecules because they can be synthesized and secreted more quickly with fewer cell resources (33). Recently, a growing number of studies have suggested the increasing acceptance of mass spectrometry-based techniques for interrogating neonatal responses to vaccines (34). However, to date, few published systems vaccinology studies have assessed the systemic metabolic responses of teenagers to immunization with the inactivated SARS-CoV-2 vaccine.
The inactivated SARS-CoV-2 vaccine has been shown to generate rapid and vast antibody responses, for instance, IgG, making it appropriate for emergency usage and critical during the COVID-19 or its variant pandemic (35–37), which is consistent with our study. Peculiarly, we detected a non-significant increase in IFN-γ levels and a reduction in IL-2 levels in the sera of vaccinated individuals. The foregoing results imply that vaccination may not trigger a relatively strong inflammatory response, as is the case with infection, and that such a relatively weak or even downregulation of inflammation may have beneficial effects on the immune reaction of the body in response to the vaccine. Meanwhile, we present the overall characterization of SARS-CoV-2 vaccination-induced changes to the undergraduate student serum metabolome by untargeted metabolomics. In summary, we found that the serum metabolic profiles of vaccinated individuals were significantly different from those before vaccination. A bulk of circulating metabolites are associated with vaccination. Taurine is a naturally occurring sulfur-containing amino acid that can be converted from other amino acids in the liver or obtained from the diet (38). Recently, a growing number of studies have found that taurine plays a critical role in regulating immune system health and antioxidation and has been reported to have anti-inflammatory effects by suppressing cytokine production (39–41). Intriguingly, taurine was remarkably upregulated in vaccinated individuals in our results. The strong and positive relationship between taurine and an immunized antibody level reinforces the suggestion that vaccination can establish a reciprocal regulatory relationship between taurine and neutralizing antibodies. The observations in this study are somewhat in concert with the results of previous research in relation to BCG-induced protective effects as with taurine (42). However, this is the first study implicating a potential role for taurine in response to SARS-CoV-2 inactivated vaccine-induced immunity, but the mechanism of action remains elusive. It is noteworthy that taurine was first identified to respond to the immune response after receiving the SARS-CoV-2 inactivated vaccine; nonetheless, the underlying mechanism of action is not clarified.
In our study, SARS-CoV-2 vaccine-induced changes in metabolite pathways correlated with the blood antibody and cytokine responses. It has been well documented that the tricarboxylic acid cycle and amino acid metabolism exert an essential influence in the modulation of the immune response to vaccines (43–45). In addition to the TCA cycle and amino acid metabolism pathways that are already known to influence host adaptive immunity, one of the most enriched metabolic pathways which was significantly associated with immune responses upon SARS-CoV-2 vaccination is glutamate metabolism. GABA is a non-protein amino acid comprising four carbon atoms that can be produced by plants, animals, and microorganisms (46, 47). In mammals, GABA has been well studied, primarily as a neurotransmitter (48, 49). Meanwhile, the receptor of GABA (GABA-Rs) plays a pivotal role in the neurotransmission process in the central system. Nevertheless, a recent report revealed that GABA-Rs can be expressed in the immune and epithelial cells as well and showed that GABA treatment can ameliorate inflammation and reduce viral load in the lungs (50). Simultaneously, our findings identified that GABA and its predecessor metabolite (L-glutamic acid) were significantly upregulated in vaccinated individuals, and GABA also has a significant positive correlation with the serum level of IgG. Specifically, GABA is known to be transformed form L-glutamic acid by glutamate decarboxylase (GAD). Numerous studies have shown that probiotics, such as Lactobacillus, possess CAD genes in their genome and are effective in producing large amounts of GABA (51–53). Notably, a very recent study revealed that Lactobacillus plantarum GUANKE (LPG) could foster SARS-CoV-2 specific immune reactions in vaccinated mice by potentiating interferon signaling, inhibiting apoptosis and inflammatory pathways. The administration of LPG to mice after the inoculation of the SARS-CoV-2 vaccine can boost and maintain the production of SARS-CoV-2-neutralizing antibodies in their bodies (54). The results described above indicate that a direct supplementation with GABA, oral administration of probiotics capable of generating GABA and stimulating antibody production, or targeting GABA-Rs could strengthen the neutralizing antibody production in the human body in response to SARS-CoV-2 vaccination. Yet, the precise regulatory mechanisms should be investigated in greater detail. In addition to GABA, the other three metabolites (L-glutamic acid, succinic acid semialdehyde, and succinic acid) enriched in the glutamate metabolism pathway were remarkably elevated in vaccinated individuals, and all of them were significantly positively correlated with the levels of either IgM, IgG, or IgA. This further led us to speculate that these metabolites may be closely related to the immune response of the organism after immunization.
In this investigation, the combination of two metabolite biomarkers we found, GABA and indole, was able to determine whether an individual had been inoculated with the SARS-CoV-2 vaccine or not with a 95.6% accuracy based on AUC value. Therefore, in the future, if we want to identify whether a population has been vaccinated with the SARS-CoV-2 vaccine on a large scale, our tentative guess is that we can base on the detected concentrations of GABA and indole in the serum of individuals. Despite the existence of indole and GABA as both metabolite markers to distinguish individuals from vaccinated individuals, a significantly negative relationship pattern was detected between indole and all neutralizing antibodies in the serum, which is contrary to the results for GABA. Indole is formed by the bacterial metabolism of tryptophan, and currently the majority of research has reported its relatively potent anti-inflammatory immune effects (55, 56). In general, the negative correlation between indole and antibodies and the significant reduction of indole in the serum of vaccinated individuals only suggest that immunizations prevent indole from impeding the production of neutralizing antibodies in the body by reducing indole levels.
Although we were successful in recruiting 30 college student volunteers to comprise the research cohort in this study and collecting blood or serum samples before SARS-CoV-2 vaccination and after receiving two doses of the vaccine. However, the cohort was still relatively small in size to obtain the best reliability. In addition, to better confirm our findings, metabolic analysis of serum samples from larger COVID-19-vaccinated cohorts and participants of various ages may be warranted.
Taken together, our metabolomic data presented in this study provide an overall view of circulating metabolite alteration from volunteers who received two doses of the inactivated SARS-CoV-2 vaccine as compared to that before vaccination. We identify four metabolites enriched in either butanoate metabolism or the alanine, aspartate, and glutamate metabolism pathways as the mainly responding small molecules upon vaccination. Correlation analysis further indicates that L-glutamic acid, GABA, succinic acid semialdehyde, and succinic acid could be used as a small molecules for potential targeted intervention to augment the immune response after SARS-CoV-2-based immunization or that of its variant.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ebi.ac.uk/metabolights/index, accession ID: MTBLS4995.
Ethics statement
The studies involving human participants were reviewed and approved by the ethics committee of Anhui Medical University Biomedical Ethics Committee (Approval No. 2021H021). The patients/participants provided their written informed consent to participate in this study.
Author contributions
SH and WJ conceived and designed the experiments and revised the manuscript; MHe, YH, and YW performed the experiments, analyzed the data, and wrote the manuscript; JL, MHa, YX, NZ, HG, HQ, and LC collected the samples and performed the experiments. All authors read and approved the final manuscript.
Funding
This work was supported by Grants for the Natural Science Foundation of China (No.81974306), Major Project of Natural Science Research of Anhui Education Department (No. KJ2019ZD23), the Project of Quality Engineering of Anhui Province (No. 2019zyrc017, 2020jyxm0864, 2021jcxkpy013), General Program of Scientific Research in Higher Education Institutions of Anhui Province (ZR2019B04), Basic and Clinical Collaboration Enhancement Program Foundation of Anhui Medical University (2020xkjT023), and Grants for Scientific Research of BSKY (XJ201922) from Anhui Medical University.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.954801/full#supplementary-material
References
1. Griffin M, Sohrabi C, Alsafi Z, Nicola M, Kerwan A, Mathew G, et al. Preparing for covid-19 exit strategies. Ann Med Surg (Lond) (2021) 61:88–92. doi: 10.1016/j.amsu.2020.12.012
2. Moore S, Hill EM, Tildesley MJ, Dyson L, Keeling MJ. Vaccination and non-pharmaceutical interventions for covid-19: A mathematical modelling study. Lancet Infect Dis (2021) 21(6):793–802. doi: 10.1016/S1473-3099(21)00143-2
4. Desai D, Khan AR, Soneja M, Mittal A, Naik S, Kodan P, et al. Effectiveness of an inactivated virus-based sars-Cov-2 vaccine, Bbv152, in India: A test-negative, case-control study. Lancet Infect Dis (2022) 22(3):349–56. doi: 10.1016/S1473-3099(21)00674-5
5. Kang M, Yi Y, Li Y, Sun L, Deng A, Hu T, et al. Effectiveness of inactivated covid-19 vaccines against illness caused by the B.1.617.2 (Delta) variant during an outbreak in guangdong, China : A cohort study. Ann Intern Med (2022) 175(4):533–40. doi: 10.7326/M21-3509
6. Xia S, Duan K, Zhang Y, Zhao D, Zhang H, Xie Z, et al. Effect of an inactivated vaccine against sars-Cov-2 on safety and immunogenicity outcomes: Interim analysis of 2 randomized clinical trials. JAMA (2020) 324(10):951–60. doi: 10.1001/jama.2020.15543
7. Levin EG, Lustig Y, Cohen C, Fluss R, Indenbaum V, Amit S, et al. Waning immune humoral response to Bnt162b2 covid-19 vaccine over 6 months. N Engl J Med (2021) 385(24):e84. doi: 10.1056/NEJMoa2114583
8. Xiang T, Liang B, Fang Y, Lu S, Li S, Wang H, et al. Declining levels of neutralizing antibodies against sars-Cov-2 in convalescent covid-19 patients one year post symptom onset. Front Immunol (2021) 12:708523. doi: 10.3389/fimmu.2021.708523
9. Marot S, Malet I, Leducq V, Zafilaza K, Sterlin D, Planas D, et al. Rapid decline of neutralizing antibodies against sars-Cov-2 among infected healthcare workers. Nat Commun (2021) 12(1):844. doi: 10.1038/s41467-021-21111-9
10. Shen B, Yi X, Sun Y, Bi X, Du J, Zhang C, et al. Proteomic and metabolomic characterization of covid-19 patient sera. Cell (2020) 182(1):59–72.e15. doi: 10.1016/j.cell.2020.05.032
11. Roberts I, Wright Muelas M, Taylor JM, Davison AS, Xu Y, Grixti JM, et al. Untargeted metabolomics of covid-19 patient serum reveals potential prognostic markers of both severity and outcome. Metabolomics (2022) 18(1):1–19. doi: 10.1007/s11306-021-01859-3
12. Wu D, Shu T, Yang X, Song JX, Zhang M, Yao C, et al. Plasma metabolomic and lipidomic alterations associated with covid-19. Natl Sci Rev (2020) 7(7):1157–68. doi: 10.1093/nsr/nwaa086
13. Roberts I, Wright Muelas M, Taylor JM, Davison AS, Xu Y, Grixti JM, et al. Untargeted metabolomics of covid-19 patient serum reveals potential prognostic markers of both severity and outcome. Metabolomics (2021) 18(1):6. doi: 10.1007/s11306-021-01859-3
14. Shi D, Yan R, Lv L, Jiang H, Lu Y, Sheng J, et al. The serum metabolome of covid-19 patients is distinctive and predictive. Metabolism (2021) 118:154739. doi: 10.1016/j.metabol.2021.154739
15. Dagla I, Iliou A, Benaki D, Gikas E, Mikros E, Bagratuni T, et al. Plasma metabolomic alterations induced by covid-19 vaccination reveal putative biomarkers reflecting the immune response. Cells (2022) 11(7):1241. doi: 10.3390/cells11071241
16. Wang Y, Wang X, Luu LDW, Chen S, Jin F, Wang S, et al. Proteomic and metabolomic signatures associated with the immune response in healthy individuals immunized with an inactivated sars-Cov-2 vaccine. Front Immunol (2022) 13:848961. doi: 10.1101/2021.07.21.21260959
17. Hasan MR, Suleiman M, Perez-Lopez A. Metabolomics in the diagnosis and prognosis of covid-19. Front Genet (2021) 12:721556. doi: 10.3389/fgene.2021.721556
18. Zhang A, Sun H, Wang X. Serum metabolomics as a novel diagnostic approach for disease: A systematic review. Anal Bioanal. Chem (2012) 404(4):1239–45. doi: 10.1007/s00216-012-6117-1
19. Buck MD, Sowell RT, Kaech SM, Pearce EL. Metabolic instruction of immunity. Cell (2017) 169(4):570–86. doi: 10.1016/j.cell.2017.04.004
20. Ganeshan K, Chawla A. Metabolic regulation of immune responses. Annu Rev Immunol (2014) 32:609–34. doi: 10.1146/annurev-immunol-032713-120236
21. Osborn O, Olefsky JM. The cellular and signaling networks linking the immune system and metabolism in disease. Nat Med (2012) 18(3):363–74. doi: 10.1038/nm.2627
22. Li S, Sullivan NL, Rouphael N, Yu T, Banton S, Maddur MS, et al. Metabolic phenotypes of response to vaccination in humans. Cell (2017) 169(5):862–77.e17. doi: 10.1016/j.cell.2017.04.026
23. Cani PD, Jordan BF. Gut microbiota-mediated inflammation in obesity: A link with gastrointestinal cancer. Nat Rev Gastroenterol Hepatol (2018) 15(11):671–82. doi: 10.1038/s41575-018-0025-6
24. Kim M, Qie Y, Park J, Kim CH. Gut microbial metabolites fuel host antibody responses. Cell Host Microbe (2016) 20(2):202–14. doi: 10.1016/j.chom.2016.07.001
25. Walsh MC, Brennan L, Malthouse JP, Roche HM, Gibney MJ. Effect of acute dietary standardization on the urinary, plasma, and salivary metabolomic profiles of healthy humans. Am J Clin Nutr (2006) 84(3):531–9. doi: 10.1093/ajcn/84.3.531
26. Yin P, Peter A, Franken H, Zhao X, Neukamm SS, Rosenbaum L, et al. Preanalytical aspects and sample quality assessment in metabolomics studies of human blood. Clin Chem (2013) 59(5):833–45. doi: 10.1373/clinchem.2012.199257
27. Ning P, Zheng Y, Luo Q, Liu X, Kang Y, Zhang Y, et al. Metabolic profiles in community-acquired pneumonia: Developing assessment tools for disease severity. Crit Care (2018) 22(1):130. doi: 10.1186/s13054-018-2049-2
28. Luo P, Yin P, Zhang W, Zhou L, Lu X, Lin X, et al. Optimization of Large-scale pseudotargeted metabolomics method based on liquid chromatography-mass spectrometry. J Chromatogr A (2016) 1437:127–36. doi: 10.1016/j.chroma.2016.01.078
29. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res (2003) 13(11):2498–504. doi: 10.1101/gr.1239303
30. Friedman J, Alm EJ. Inferring correlation networks from genomic survey data. PloS Comput Biol (2012) 8(9):e1002687. doi: 10.1371/journal.pcbi.1002687
31. Bader GD, Hogue CW. An automated method for finding molecular complexes in Large protein interaction networks. BMC Bioinf (2003) 4:2. doi: 10.1186/1471-2105-4-2
32. Danlos FX, Grajeda-Iglesias C, Durand S, Sauvat A, Roumier M, Cantin D, et al. Metabolomic analyses of covid-19 patients unravel stage-dependent and prognostic biomarkers. Cell Death Dis (2021) 12(3):258. doi: 10.1038/s41419-021-03540-y
33. Zuo T, Liu Q, Zhang F, Lui GC, Tso EY, Yeoh YK, et al. Depicting sars-Cov-2 faecal viral activity in association with gut microbiota composition in patients with covid-19. Gut (2021) 70(2):276–84. doi: 10.1136/gutjnl-2020-322294
34. Petrick LM, Schiffman C, Edmands WMB, Yano Y, Perttula K, Whitehead T, et al. Metabolomics of neonatal blood spots reveal distinct phenotypes of pediatric acute lymphoblastic leukemia and potential effects of early-life nutrition. Cancer Lett (2019) 452:71–8. doi: 10.1016/j.canlet.2019.03.007
35. Stuart ASV, Shaw RH, Liu X, Greenland M, Aley PK, Andrews NJ, et al. Immunogenicity, safety, and reactogenicity of heterologous covid-19 primary vaccination incorporating mrna, viral-vector, and protein-adjuvant vaccines in the uk (Com-Cov2): A single-blind, randomised, phase 2, non-inferiority trial. Lancet (2022) 399(10319):36–49. doi: 10.1016/S0140-6736(21)02718-5
36. Wu Z, Hu Y, Xu M, Chen Z, Yang W, Jiang Z, et al. Safety, tolerability, and immunogenicity of an inactivated sars-Cov-2 vaccine (Coronavac) in healthy adults aged 60 years and older: A randomised, double-blind, placebo-controlled, phase 1/2 clinical trial. Lancet Infect Dis (2021) 21(6):803–12. doi: 10.1016/S1473-3099(20)30987-7
37. Xia S, Zhang Y, Wang Y, Wang H, Yang Y, Gao GF, et al. Safety and immunogenicity of an inactivated covid-19 vaccine, bbibp-corv, in people younger than 18 years: A randomised, double-blind, controlled, phase 1/2 trial. Lancet Infect Dis (2022) 22(2):196–208. doi: 10.1016/S1473-3099(21)00462-X
38. Jong CJ, Sandal P, Schaffer SW. The role of taurine in mitochondria health: More than just an antioxidant. Molecules (2021) 26(16):4913. doi: 10.3390/molecules26164913
39. Dong JF, Zheng XQ, Rui HB. Effect of taurine on immune function in mice with T-cell lymphoma during chemotherapy. Asian Pac J Trop Med (2017) 10(11):1090–4. doi: 10.1016/j.apjtm.2017.10.014
40. Qaradakhi T, Gadanec LK, McSweeney KR, Abraham JR, Apostolopoulos V, Zulli A. The anti-inflammatory effect of taurine on cardiovascular disease. Nutrients (2020) 12(9):2847. doi: 10.3390/nu12092847
41. Silva LA, Silveira PC, Ronsani MM, Souza PS, Scheffer D, Vieira LC, et al. Taurine supplementation decreases oxidative stress in skeletal muscle after eccentric exercise. Cell Biochem Funct (2011) 29(1):43–9. doi: 10.1002/cbf.1716
42. Koeken VA, Mourits VP, de Bree LCJ, Moorlag JCFM S, Sonawane V, Lemmers H, et al. Plasma metabolome predicts trained immunity responses after antituberculosis BCG vaccination. Plos Biol (2022) 20(9):e3001765. doi: 10.1371/journal.pbio.3001765
43. Choi I, Son H, Baek JH. Tricarboxylic acid (Tca) cycle intermediates: Regulators of immune responses. Life (Basel) (2021) 11(1). doi: 10.3390/life11010069
44. Yang J, Yang XL, Su YB, Peng XX, Li H. Activation of the tca cycle to provide immune protection in zebrafish immunized by high magnesium-prepared vibrio alginolyticus vaccine. Front Immunol (2021) 12:739591. doi: 10.3389/fimmu.2021.739591
45. Li P, Yin YL, Li D, Kim SW, Wu G. Amino acids and immune function. Br J Nutr (2007) 98(2):237–52. doi: 10.1017/S000711450769936X
46. Rashmi D, Zanan R, John S, Khandagale K, Nadaf A. Γ-aminobutyric acid (Gaba): Biosynthesis, role, commercial production, and applications. Stud Natural prod. Chem (2018) 57:413–52. doi: 10.1016/B978-0-444-64057-4.00013-2.
47. Calvani R, Miccheli A, Capuani G, Tomassini Miccheli A, Puccetti C, Delfini M, et al. Gut microbiome-derived metabolites characterize a peculiar obese urinary metabotype. Int J Obes (Lond) (2010) 34(6):1095–8. doi: 10.1038/ijo.2010.44
48. Leventhal AG, Wang Y, Pu M, Zhou Y, Ma Y. Gaba and its agonists improved visual cortical function in senescent monkeys. Science (2003) 300(5620):812–5. doi: 10.1126/science.1082874
49. Xu N, Wei L, Liu J. Biotechnological advances and perspectives of gamma-aminobutyric acid production. World J Microbiol Biotechnol (2017) 33(3):64. doi: 10.1007/s11274-017-2234-5
50. Tian J, Middleton B, Kaufman DL. Gaba administration limits viral replication and pneumonitis in a mouse model of covid-19. bioRxiv (2021). doi: 10.1101/2021.02.09.430446
51. Cui Y, Miao K, Niyaphorn S, Qu X. Production of gamma-aminobutyric acid from lactic acid bacteria: A systematic review. Int J Mol Sci (2020) 21(3):995. doi: 10.3390/ijms21030995
52. Ruiz Rodriguez L, Vera Pingitore E, Rollan G, Cocconcelli PS, Fontana C, Saavedra L, et al. Biodiversity and technological-functional potential of lactic acid bacteria isolated from spontaneously fermented quinoa sourdoughs. J Appl Microbiol (2016) 120(5):1289–301. doi: 10.1111/jam.13104
53. Wu Q, Tun HM, Law YS, Khafipour E, Shah NP. Common distribution of gad operon in lactobacillus brevis and its gada contributes to efficient gaba synthesis toward cytosolic near-neutral ph. Front Microbiol (2017) 8:206. doi: 10.3389/fmicb.2017.00206
54. Xu J, Ren Z, Cao K, Li X, Yang J, Luo X, et al. Boosting vaccine-elicited respiratory mucosal and systemic covid-19 immunity in mice with the oral lactobacillus plantarum. Front Nutr (2021) 8:789242. doi: 10.3389/fnut.2021.789242
55. Whitfield-Cargile CM, Cohen ND, Chapkin RS, Weeks BR, Davidson LA, Goldsby JS, et al. The microbiota-derived metabolite indole decreases mucosal inflammation and injury in a murine model of nsaid enteropathy. Gut Microbes (2016) 7(3):246–61. doi: 10.1080/19490976.2016.1156827
Keywords: COVID-19, SARS-CoV-2 vaccination, serum, antibodies and cytokines, metabolomic analysis
Citation: He M, Huang Y, Wang Y, Liu J, Han M, Xiao Y, Zhang N, Gui H, Qiu H, Cao L, Jia W and Huang S (2022) Metabolomics-based investigation of SARS-CoV-2 vaccination (Sinovac) reveals an immune-dependent metabolite biomarker. Front. Immunol. 13:954801. doi: 10.3389/fimmu.2022.954801
Received: 27 May 2022; Accepted: 12 September 2022;
Published: 28 September 2022.
Edited by:
Alagarraju Muthukumar, University of Texas Southwestern Medical Center, United StatesCopyright © 2022 He, Huang, Wang, Liu, Han, Xiao, Zhang, Gui, Qiu, Cao, Jia and Huang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shenghai Huang, huangshh68@aliyun.com; Weihua Jia, jwhua29@163.com
†These authors have contributed equally to this work