Summary
The antituberculosis vaccine Bacillus Calmette–Guérin (BCG) induces nonspecific safety in opposition to heterologous infections, not less than partly by induction of innate immune reminiscence (educated immunity). The amplitude of the response to BCG is variable, however the components that affect this response are poorly understood. Metabolites, both launched by cells or absorbed from the intestine, are identified to affect immune responses, however whether or not they impression BCG responses will not be identified. We vaccinated 325 wholesome people with BCG, and picked up blood earlier than, 2 weeks and three months after vaccination, to evaluate the affect of circulating metabolites on the immune responses induced by BCG. Circulating metabolite concentrations after BCG vaccination had been discovered to have a extra pronounced impression on educated immunity responses, reminiscent of the rise in IL-1β and TNF-α manufacturing upon Staphylococcus aureus stimulation, than on particular adaptive immune reminiscence, assessed as IFN-γ manufacturing in response to Mycobacterium tuberculosis. Circulating metabolites at baseline had been in a position to predict educated immunity responses at 3 months after vaccination and enrichment evaluation based mostly on the metabolites positively related to educated immunity revealed enrichment of the tricarboxylic acid (TCA) cycle and glutamine metabolism, each of which had been beforehand discovered to be essential for educated immunity. A number of new metabolic pathways that affect educated immunity had been recognized, amongst which taurine metabolism related to BCG-induced educated immunity, a discovering validated in useful experiments. In conclusion, circulating metabolites are essential components influencing BCG-induced educated immunity in people. Modulation of metabolic pathways could also be a novel technique to enhance vaccine and educated immunity responses.
Quotation: Koeken VACM, Qi C, Mourits VP, de Bree LCJ, Moorlag SJCFM, Sonawane V, et al. (2022) Plasma metabolome predicts educated immunity responses after antituberculosis BCG vaccination. PLoS Biol 20(9):
e3001765.
https://doi.org/10.1371/journal.pbio.3001765
Tutorial Editor: Jason W. Locasale, Duke College, UNITED STATES
Acquired: September 28, 2021; Accepted: July 27, 2022; Revealed: September 12, 2022
Copyright: © 2022 Koeken et al. That is an open entry article distributed beneath the phrases of the Inventive Commons Attribution License, which allows unrestricted use, distribution, and replica in any medium, offered the unique writer and supply are credited.
Knowledge Availability: All of the metadata on research individuals, metabolome, and cytokine information from the 300BCG cohort research that had been used for the analyses described on this manuscript, can be found at https://gitlab.com/xavier-lab-computation/public/bcg300. The information generated throughout in vitro validation experiments can be found as supplementary information.
Funding: YL was supported by an ERC beginning Grant (948207) and a Radboud College Medical Centre Hypatia Grant (2018). CJX was supported by Helmholtz Initiative and Networking Fund (1800167). MGN was supported by an ERC Superior Grant (833247) and a Spinoza Grant of the Netherlands Group for Scientific Analysis. The funders had no position in research design, information assortment and evaluation, choice to publish, or preparation of the manuscript.
Competing pursuits: The authors have declared that no competing pursuits exist.
Abbreviations:
BCG,
Bacillus Calmette–Guérin; BMI,
physique mass index; FDR,
false discovery charge; HMDB,
human metabolites database; LDH,
lactate dehydrogenase; MTA,
methyltioadenosine; PAMP,
pathogen-associated molecular sample; PBS,
phosphate-buffered saline; PC,
principal element; PCA,
principal element evaluation; PBMC,
peripheral blood mononuclear cell; TCA,
tricarboxylic acid; WGCNA,
weighted correlation community evaluation
Introduction
Sure microbial ligands and vaccinations can go away an immunological imprint on the innate immune system, resulting in enhanced responsiveness in a while. This de facto innate immune reminiscence, which has been termed educated immunity, includes the useful reprogramming of innate immune cells reminiscent of monocytes and NK-cells [1]. Skilled immunity can present nonspecific safety in opposition to unrelated pathogens, as has been demonstrated for the antituberculosis vaccine Bacillus Calmette–Guérin (BCG) [2,3]. BCG vaccination ends in enhanced ex vivo cytokine manufacturing by peripheral blood mononuclear cells (PBMCs) and NK-cells in response to unrelated pathogens reminiscent of Staphylococcus aureus and Candida albicans [4,5]. These educated immunity responses are regulated at epigenetic and transcriptional stage [6,7], with metabolic adjustments being proven to be important for the epigenetic modifications crucial for educated immunity [1]. BCG coaching up-regulates each glycolysis and oxidative phosphorylation in monocytes, and pharmacological inhibition of glycolysis ends in abrogated educated immunity responses [8]. In vitro priming of monocytes with BCG or the fungal cell wall element β-glucan additionally prompts glutaminolysis and accumulation of fumarate, which contributes to the induction of educated immunity by inhibition of the KDM5 histone demethylases [8,9].
In vivo responses to BCG vaccination are variable. Some people have a powerful educated immunity and adaptive immune reminiscence response after vaccination, whereas others reply poorly. That is evident in each experimental research by which the epigenetic and immunological assessments have proven good responders and poor responders [10,11], in addition to within the epidemiological research which have proven a partial safety of BCG at populational stage [12]. It’s poorly understood which host and environmental components account for the interindividual variability of immune responses to BCG vaccination. We hypothesized that interindividual variations in host metabolism and circulating metabolite concentrations affect this variability. In an effort to analyze this in vivo, we performed a big cohort research by which we vaccinated greater than 300 wholesome people with BCG (300BCG research), adopted by the evaluation of the induction of educated immunity (heterologous induction of monocyte-derived cytokines in response to S. aureus) and adaptive immune reminiscence (induction of T-cell responses by Mycobacterium tuberculosis) [13]. On this research, we examined if differential abundance of circulating metabolites impacts induction of educated immunity and/or adaptive immune responses induced by BCG in vivo.
Outcomes
Host components related to circulating metabolites
For a complete measure of the circulating metabolome, flow-injection TOF-M was carried out to profile a complete of 1,607 metabolic options. After eradicating metabolites with an annotation rating beneath 70, 1,373 options had been left for additional evaluation. The metabolic measurements had been carried out in plasma samples from 325 wholesome volunteers of Western European descent earlier than BCG vaccination was administered (Fig 1A). Of the 325 individuals, 18 had been excluded as a result of they acquired their vaccination within the night, and we’ve got realized that vaccination efficacy within the night is suboptimal [14]. 4 individuals had been excluded attributable to medicine use or febrile sickness, leaving 303 samples for additional evaluation. Of those 303 individuals, 56% had been feminine, the imply age was 23 years (vary 18 to 71) and the imply physique mass index (BMI) was 22.5 (±2.6 commonplace deviation).
Fig 1. Examine overview and correlation between circulating metabolome and ex vivo cytokine responses upon BCG vaccination.
(A) Schematic overview of the research. Created with BioRender.com. (B) Heatmap of hierarchical clustering on Spearman correlation patterns between metabolites and PBMC-derived S. aureus-induced IL-1β, IL-6, and TNF-α responses and M. tuberculosis-induced IL-1β, IL-6, and IFN-γ responses. The ex vivo cytokine responses are analyzed as fold adjustments between baseline (v1, day 0) and a couple of weeks (v2, day 14) and three months (v3, day 90) after BCG vaccination. The metabolites are represented as the highest 42 PCs, capturing 75.3% of the variance. Cell colours point out −log10(p-value) multiplied by the correlation coefficients from unfavourable (blue) to constructive (purple). The metadata on research individuals, metabolome, and cytokine information used to generate this determine can be found at https://gitlab.com/xavier-lab-computation/public/bcg300. The loadings of the highest 10 metabolites of the highest 10 PC’s might be present in S1 Desk. BCG, Bacillus Calmette–Guérin; PC, principal element; PBMC, peripheral blood mononuclear cell.
We assessed the affiliation between numerous host components and the circulating metabolome. Intercourse and age each had a profound impact on circulating metabolites, with 35.3% of all metabolites considerably related to intercourse and 38.2% of all metabolites related to age (false discovery charge [FDR] <0.05, linear regression evaluation). After correction for age and intercourse, solely 0.95% of all metabolites had been considerably related to BMI and 0.22% of all metabolites had been considerably related to smoking. We additionally examined the affiliation between metabolites and high 20 principal elements (PCs) calculated from genotyping information, however no important associations had been recognized at FDR <0.05. In feminine individuals (N = 171), 17.0% of the metabolites had been related to oral contraceptive use. Vegetarian weight loss plan (reported by 8% of research topics) was solely related to adjustments in concentrations of 0.36% of measured metabolites, amongst which methylhistidine as reported beforehand (FDR = 3.9 × 10−3, linear regression evaluation, S1 Fig) [15].
Baseline metabolites related to innate and adaptive immune responses upon BCG vaccination
Subsequent, we examined if circulating metabolites at baseline are related to the immune response upon vaccination. To scale back information dimensionality, a principal element evaluation (PCA) was carried out on the circulating metabolome information (S1 Desk). The primary 42 PCs of the metabolic information had been used for additional evaluation, capturing 75.3% of the entire variance within the metabolome information (S2 Fig). Subsequently, we examined if these PCs representing the circulating metabolome predicted the innate and adaptive immune responses upon BCG vaccination. The immune response was assessed by ex vivo stimulation of PBMCs, and we assessed 2 immunological finish factors (S3 Fig):
- Skilled immunity responses: comparability of cytokine manufacturing after heterologous stimulation with S. aureus of PBMCs remoted earlier than and after BCG vaccination.
- Adaptive reminiscence responses: comparability of the manufacturing of cytokines after M. tuberculosis ex vivo stimulation of PBMCs remoted earlier than and after BCG vaccination.
To seize the interactions, Spearman correlations between the PCs representing the baseline metabolome and the ex vivo cytokine responses had been calculated, and unsupervised hierarchical clustering evaluation was carried out on the correlation coefficients. Based mostly on correlation patterns, a number of associations between metabolite modules and ex vivo cytokine responses upon BCG vaccination had been discovered (Fig 1B). For instance, PC7 was positively correlated with educated immunity responses, as measured by the rise in IL-1β, IL-6, and TNF-α manufacturing in response to S. aureus 3 months after BCG vaccination. The enriched metabolic pathways throughout the metabolite subset contributing to PC7 included arginine biosynthesis (p = 5.2 × 10−3, Fisher’s precise check), aminoacyl-tRNA biosynthesis (p = 7.3 × 10−3, Fisher’s precise check), and phenylalanine, tyrosine, and tryptophan biosynthesis (p = 4.1 × 10−2, Fisher’s precise check). Throughout the arginine biosynthesis metabolic pathway, metabolites together with glutamine, fumarate, and aspartate contributed to PC7. As well as, PC5 and PC9 had been extremely correlated to age, once more supporting the profound impact of age on the circulating metabolome.
To additional check which metabolites contribute to ex vivo vaccine responses, we assessed the associations of all stimuli-cytokine mixtures and particular person metabolites corrected for age and intercourse (S2 Desk). As proven in Fig 2A, baseline metabolite concentrations particularly correlated with educated immunity responses, reminiscent of S. aureus-induced IL-1β and TNF-α manufacturing 3 months after BCG vaccination, fairly than with particular adaptive immune reminiscence, measured as M. tuberculosis-induced IFN-γ manufacturing (Fisher’s precise check, p = 5.0 × 10−4). Solely 7% of the metabolites related to the rise in IFN-γ in response to M. tuberculosis 3 months after BCG vaccination (p < 0.05), whereas 29% and 21% of the metabolites had been related to the rise in S. aureus-induced IL-1β and TNF-α, respectively. One of many key metabolites positively related to IFN-γ responses induced by M. tuberculosis was 5-aminolevulinic acid, which is the primary metabolite within the porphyrin synthesis pathway resulting in the formation of heme [16]. Two essential metabolites of the tricarboxylic acid (TCA) cycle, plasma malate and succinate, had been strongly related to elevated manufacturing of TNF-α (p = 8.4 × 10−4, affiliation with malate, linear regression evaluation) and IL-1β (p = 7.3 × 10−5, affiliation with succinate, linear regression evaluation) after BCG vaccination, when PBMCs had been stimulated with S. aureus (Fig 2B and 2C). Nonlinear regression, used to check the robustness of this linear relationship, confirmed related outcomes because the linear regression mannequin (S4 Fig and S3 Desk).
Fig 2. Affiliation between baseline metabolites and ex vivo cytokine responses upon BCG vaccination.
(A) Baseline metabolites had been related to fold adjustments of PBMC-derived S. aureus-induced IL-1β, IL-6, and TNF-α responses and M. tuberculosis-induced IL-1β, IL-6, and IFN-γ responses corrected for age and intercourse. Fold adjustments had been calculated 2 weeks (v2, day 14) or 3 months (v3, day 90) after BCG vaccination in comparison with baseline (v1, day 0). Solely metabolites with a major affiliation to one of many cytokine-stimulus pairs are depicted on this determine (N = 777, p < 0.05). The metabolites are ordered by hierarchical clustering. Cell colours point out important associations as −log10(p-value) multiplied by the correlation coefficients from unfavourable (blue) to constructive (purple). The abstract statistics of high 10 related metabolites per stimulation might be present in S2 Desk. (B) Linear regression between malate at baseline and fold adjustments of ex vivo PBMC-derived S. aureus-induced TNF-α responses 3 months after vaccination. (C) Linear regression between succinate at baseline and fold adjustments of ex vivo PBMC-derived S. aureus-induced IL-1β responses 3 months after vaccination. The metadata on research individuals, metabolome, and cytokine information used to generate these figures can be found at https://gitlab.com/xavier-lab-computation/public/bcg300. BCG, Bacillus Calmette–Guérin; PBMC, peripheral blood mononuclear cell.
Baseline metabolites predict educated immunity responses upon BCG vaccination
Subsequent, we assessed the predictive means of baseline metabolites for the educated immunity responses. On this evaluation, we aimed to foretell the fold adjustments of S. aureus-induced cytokines 3 months after vaccination, as a big proportion of baseline metabolites had been related to these responses. First, the dataset was randomly break up right into a coaching set (70% of the info) and a check set (30%), after which 10-fold cross-validation was carried out to pick options, which was repeated 5 instances, utilizing Elastic Internet [17]. Then, the accuracy of every mannequin was examined by calculating a Pearson correlation between the anticipated values based mostly on the mannequin and the measured values within the check dataset (30% of the observations), as a measure of the similarity between the anticipated and measured cytokine fold adjustments. This entire process was repeated 100 instances, producing 100 fashions. The correlation coefficients, as a measure of the prediction efficiency of those fashions, are visualized in Fig 3A. The median correlations between predicted and measured cytokine fold adjustments had been 0.25 for IL-1β, 0.13 for IL-6, and 0.21 for TNF-α. Thus, the predictive energy of baseline metabolites for educated immunity responses was comparatively low, though statistically important. The predictive energy and stability of the mannequin can in all probability be elevated by together with multi-omics information from totally different layers, e.g., epigenetics [18]. Then, we analyzed which metabolites occurred most steadily within the prediction fashions for IL-1β (Fig 3B), IL-6 (Fig 3C), and TNF-α (Fig 3D) in response to S. aureus 3 months after vaccination. Apparently, metabolites of the TCA cycle once more appeared in our outcomes, as fumarate was 65 out of 100 instances chosen within the mannequin to foretell larger IL-1β manufacturing, and malate was chosen 57 instances to foretell larger TNF-α manufacturing. As beforehand described, fumarate itself is ready to induce educated immunity in an in vitro setting (S5 Fig), whereas malate will not be [9]. Different metabolites that occurred in excessive frequency had been malonyl carnitine, which is an acyl carnitine; formiminoglutamic acid, also called formimino-glu, which is an intermediate within the catabolism of histidine; and 5-hydroxytryptophol, which is a serotonin metabolite. We noticed that metabolites defined a bigger share of the variation in educated immunity responses than age, intercourse, and BMI (Fig 3E). The highest 10 metabolites from every rank listing from the prediction outcomes defined 25.2% variance for IL-1β, 17.0% for IL-6, and 24.3% for TNF-α fold adjustments in response to S. aureus.
Fig 3. Baseline metabolites predict educated immunity responses 3 months after BCG vaccination.
(A) The dataset was randomly break up right into a coaching set (70% of the observations) and a check set (30%). Then, a 10-fold cross-validation was carried out for characteristic choice, which was repeated 5 instances, to generate a prediction mannequin for IL-1β, IL-6, and TNF-α fold adjustments in response to S. aureus 3 months after BCG vaccination. Prediction accuracy was evaluated by calculation of a Pearson correlation between the measured cytokine fold adjustments and the predictions of the check units. This process was repeated 100 instances, and the correlation coefficients are visualized as boxplots. (B–D) The highest 20 metabolites occurring within the prediction fashions for IL-1β (B), IL-6 (C), and TNF-α (D) fold adjustments in response to S. aureus 3 months after BCG vaccination are visualized. In whole, 100 fashions had been constructed for every cytokine fold change. (E) Variance defined by metabolites (high 10 most frequently occurring metabolites from the prediction evaluation), age, intercourse, and BMI on educated immunity responses. Values point out the variance defined by metabolites. The metadata on research individuals, metabolome, and cytokine information used to generate these figures can be found at https://gitlab.com/xavier-lab-computation/public/bcg300. BCG, Bacillus Calmette–Guérin; BMI, physique mass index.
Metabolite co-accumulation networks related to educated immunity responses
Contemplating the interindividual variation in educated immunity responses induced by BCG, each excessive and low responders might be recognized based mostly on the rise in ex vivo cytokine responses. To establish metabolic networks underlying the responder phenotype, we assessed the metabolites co-accumulation networks in excessive and low responders when it comes to their educated immunity response (as outlined by the fold change in IL-1β after S. aureus stimulation, Fig 4A) and recognized 14 modules of extremely correlated metabolites (Fig 4B and S4 Desk). Subsequently, we correlated the 14 modules to educated immunity responses (excessive/low), age, intercourse, and BMI. Three modules (inexperienced, magenta, and turquoise) had been related to educated immunity response (Fig 4C); in a sensitivity evaluation adjusting for age, intercourse, and BMI, the identical 3 modules had been related to educated immunity (S5 Desk). Metabolites from module magenta had been enriched for the pantothenate and CoA biosynthesis (p = 0.026) and TCA cycle pathways (p = 0.029, Fig 4D and S6 Desk), once more linking TCA metabolites to BCG-induced educated immunity responses. As well as, metabolites from module turquoise had been enriched in purine and sphingolipid metabolism (Fig 4E and S6 Desk). We didn’t observe any pathways enriched based mostly on the metabolites of the inexperienced module.
Fig 4. Co-accumulated metabolites modules related to educated immunity response stage.
(A) Distribution of excessive and low responders based mostly on the fold change within the ex vivo manufacturing of IL-1β in response to S. aureus 3 months after vaccination in comparison with baseline. (B) WGCNA was carried out on metabolites from excessive and low responders of educated immunity. The metabolites included in every module recognized by WGCNA might be present in S4 Desk. (C) The recognized 14 modules had been correlated to excessive/low responder classification, age, intercourse, and BMI, respectively. Stars point out important correlation with * p < 0.05, ** p < 0.01, *** p < 0.001. (D and E) Dot plots displaying outcomes of pathways evaluation of metabolites from module magenta (D) and turquoise (E). The abstract statistics might be present in S6 Desk. No important pathways from metabolites of the inexperienced module had been recognized. The metadata on research individuals, metabolome, and cytokine information used to generate these figures can be found at https://gitlab.com/xavier-lab-computation/public/bcg300. BMI, physique mass index; WGCNA, weighted correlation community evaluation.
Metabolic pathways and taurine metabolism related to BCG-induced educated immunity
To establish the metabolic pathways concerned, we included all metabolites related to the ex vivo cytokine responses upon BCG vaccination (linear regression evaluation with age and intercourse as covariates, p < 0.05) in subsequent pathway analyses. The pathway analyses had been run independently for every cytokine-stimulus mixture per time level and individually for unfavourable and constructive associations. The pathways enriched for metabolites positively related to IL-1β and TNF-α in response to S. aureus 3 months after BCG vaccination are visualized in Fig 5A.
Fig 5. Pathway evaluation of metabolites related to educated immunity and the impact of taurine.
(A) Pathway evaluation of metabolites which might be considerably (p < 0.05) related to the fold change in IL-1β and TNF-α in response to S. aureus 3 months after BCG vaccination. Pathway analyses had been carried out utilizing MetaboAnalyst (KEGG model October 2019) [19]. (B) Linear regression between taurine at baseline and fold adjustments of ex vivo PBMC-derived S. aureus-induced IL-1β responses 3 months after vaccination. The metadata on research individuals, metabolome, and cytokine information used to generate these figures can be found at https://gitlab.com/xavier-lab-computation/public/bcg300. (C) Human major monocytes had been incubated for twenty-four hours with taurine (0.1, 0.5, or 1 mM), after which the medium was refreshed, and the cells had been allowed to relaxation for five days, after which they had been stimulated with E. coli LPS (10 ng/mL) for twenty-four hours. Then, the degrees of IL-6 and TNF-α had been measured, and a fold change was calculated relative to the medium management. The median values are offered, and every donor is represented in a unique shade (N = 6, Wilcoxon matched-pairs signed rank check, ns = not important, * p < 0.05). The cytokine values used to generate this determine might be present in S2 Knowledge. BCG, Bacillus Calmette–Guérin; PBMC, peripheral blood mononuclear cell.
For the primary time, we present that lots of the pathways that had been implicated in educated immunity from earlier in vitro educated immunity research [9,10], reminiscent of TCA cycle, glutaminolysis, and pyruvate metabolism, are additionally related to BCG-induced educated immunity in vivo. One of many pathways that has not been beforehand implicated in educated immunity was taurine metabolism. Taurine was positively related to induction of educated immunity, as assessed by fold-increase in IL-1β (Fig 5B) and TNF-α (S6 Fig) in response to S. aureus 3 months after BCG vaccination (linear regression evaluation, all talked about associations p < 0.05). To analyze a useful hyperlink between taurine and educated immunity, we carried out an in vitro educated immunity experiment by which we primed human monocytes with taurine, adopted by restimulation of those cells with LPS after 5 days of relaxation [19]. Certainly, preincubation of cells with taurine led to induction of educated immunity, as proven by elevated IL-6 and TNF-α manufacturing upon restimulation with LPS, that are the generally used read-outs to evaluate educated immunity in vitro (Fig 5C). Apparently, nevertheless, taurine didn’t amplify the induction of educated immunity induced by BCG or β-glucan (S7 Fig). Not one of the examined circumstances affected cell viability after priming (S8 Fig).
Contemplating the numerous impression of meat and fish consumption on metabolic profiles, together with taurine [20], we subsequent examined if the consumption of meat or fish was related to altered ex vivo cytokine responses upon BCG vaccination. Nevertheless, reported consumption of meat and fish (reported as consumption days per week) confirmed no relation with ex vivo cytokine responses upon BCG vaccination, presumably since day-to-day variation in meals consumption was not captured in our questionnaire.
Dialogue
BCG vaccination induces heterologous safety in opposition to an infection illnesses apart from tuberculosis, not less than partly by induction of educated immunity [1,10]. Induction of educated immunity results in enhanced responsiveness of innate immune cells and relies on epigenetic and metabolic reprogramming. BCG effectiveness to induce educated immunity is broadly heterogeneous between people [10,11], however the supply of this variation will not be identified. We hypothesized that circulating metabolites are an essential issue influencing each innate and adaptive reminiscence responses, and that they symbolize an essential supply of variation in vaccine responses between people. Certainly, in a big cohort of vaccinated wholesome volunteers, we noticed that the circulating metabolomic signature predicts BCG-induced immune responses, particularly educated immunity responses.
A lot of circulating metabolites had been related to educated immunity responses, and our in vivo research validates a number of findings from earlier in vitro research. For instance, earlier in vitro research have proven that glutamine metabolism is elevated upon BCG coaching and that inhibition of glutaminolysis results in decreased educated immunity responses [8]. On the identical time, metabolites that accumulate upon glutaminolysis in educated cells, reminiscent of fumarate, can induce educated immunity by inhibition of KDM5 demethylase, and enhance in H3K4me3 on the promotors of immune genes [9]. These in vitro information are validated by our research, by which we additionally noticed that metabolites concerned in glutamine and glutamate metabolism (reminiscent of succinate and malate) are related to BCG-induced educated immunity. As well as, succinate was beforehand proven to be concerned in innate immune signaling and associated to enhanced IL-1β responses in macrophages [21]. These information argue that variation within the concentrations of TCA metabolites within the circulation represents an essential issue within the response to BCG vaccination.
Apparently, the affiliation between circulating metabolites and BCG-induced immune responses was a lot stronger for heterologous S. aureus-induced cytokine responses at 3 months fairly than 2 weeks, particularly for IL-1β and TNF-α. Though we observe a rise in cytokine responses at each 2 weeks in addition to 3 months after vaccination, the early time level might be extra reflective of acute an infection, as stay BCG remains to be detectable 2 weeks after vaccination [22], whereas the 3-month time level actually displays innate immune reminiscence.
We needed to pick for additional validation and dialogue novel pathways that based mostly on organic arguments would doubtless be concerned in educated immunity, fairly than already identified metabolites beforehand proven to modulate it. We due to this fact chosen taurine metabolism for additional exploration as a result of little is understood about its position in BCG-induced educated immunity, for its identified anti-inflammatory results, and since it may be present in excessive concentrations in phagocytes. Taurine, or 2-aminoethanesulfonic acid, is a metabolite broadly current in animal tissue that is essential for cardiovascular operate but additionally for the event of muscle mass, retina, and nervous tissue. Taurine might be produced endogenously, however it’s also obtained from weight loss plan: each meat and fish are wealthy in taurine, whereas a vegan weight loss plan supplies nearly no taurine. Apparently, taurine has been reported to have anti-inflammatory results by inhibition of cytokine manufacturing [23,24] and is current in excessive concentrations in phagocytes, particularly in neutrophils, the place it has each anti-inflammatory and antioxidant capabilities [25,26]. In our research, taurine concentrations had been related to educated immunity induced by BCG vaccination, and extra in vitro experiments revealed that taurine is ready to induce educated immunity itself, with none microbial ligands current. That is the primary research implicating a possible position for taurine in educated immunity, and its mechanism of motion must be investigated in future research. Relating to the impression of taurine, a limitation of the present research is that we might validate the induction of educated immunity by taurine itself, however within the in vitro experimental mannequin, taurine didn’t additional enhance the educated immunity induced by BCG. We hypothesize that within the in vitro mannequin, the induction of educated immunity by BCG is already very efficient, reaching a most that can not be additional elevated by taurine. The mechanism by which taurine induces educated immunity doesn’t appear to be as performing as a pathogen-associated molecular sample (PAMP), as taurine itself (with out LPS restimulation) didn’t induce cytokine manufacturing. That is much like the results of fumarate, which doesn’t act as a PAMP, however as an inhibitor of epigenetic enzymes (in that case KDM5 histone demethylase), as proven by Arts and colleagues [9]. Sadly, little or no is understood concerning the impression of taurine on epigenetic processes. We have now certainly tried to evaluate it through the use of normal inhibitors of histone methylation reminiscent of methyltioadenosine (MTA); nevertheless, MTA exerted related results on cytokine induction by LPS both with or with out taurine (not proven), arguing that taurine exerts its results by a unique mechanism. The molecular mechanism by which taurine exerts its impact must be comprehensively investigated in future research.
A stunning remark from our research is that the circulating metabolome appears to have a extra pronounced impact on educated immunity responses in comparison with the adaptive immune responses induced by BCG vaccination. This may increasingly appear surprising, as mobile metabolism is essential for T-cell operate [27]. Whereas mobile metabolism is essential for each innate and adaptive immune responses, it appears that evidently the baseline abundance of metabolites extra strongly impacts variation in innate fairly than adaptive immune responses induced by BCG. These outcomes are in keeping with a current research that noticed a stronger impact of baseline metabolites on interindividual variation of monocyte-derived cytokines (TNF-α, IL-1β, IL-6) in comparison with T cell-derived cytokines (IL-17, IL-22, IFN-γ) [28]. Varied metabolites affect the exercise of epigenetic modifying enzymes, for instance, acetyl-CoA, which acts as an acetyl group donor for histone acetyltransferases [29]. Since epigenetic transforming is the molecular foundation for educated immunity, this would possibly clarify why the abundance of sure metabolites has a pronounced impact on BCG-induced educated immunity. The remark that metabolites present much less affiliation with IFN-γ manufacturing in response to M. tuberculosis can also be attention-grabbing for tuberculosis vaccine growth, as this implies that the metabolome may not be a significant component influencing the variable safety offered by BCG in opposition to tuberculosis. Nevertheless, our lack of knowledge of correlates of safety within the context of vaccination in opposition to TB makes it tough to guage these ends in the context of tuberculosis vaccination [30].
Our research is proscribed by the truth that the noticed hyperlink between metabolome and educated immunity depends on associative analyses (except taurine), and it stays to be investigated if there exists a causal relationship. As well as, the observations from this research ought to be validated in an impartial cohort of BCG-vaccinated people. Evaluation of intracellular cytokine staining along with the measurement of secreted cytokines would additional enrich the dataset of the research, however we had been sadly unable to carry out these experiments attributable to logistical limitations. Contemplating that BCG is principally administered in neonates in tuberculosis endemic international locations, it stays to be investigated if these outcomes might be translated to totally different populations. Moreover, the pathway enrichment analyses carried out can help interpretation and connection of metabolites that had been related to BCG-induced immune responses. Nevertheless, it is very important understand that these are mobile metabolism pathways, and that some metabolites essential for mobile metabolism will not be measurable in plasma, and that not all circulating metabolites are a part of mobile metabolism pathways. Moreover, though we had been in a position to measure and analyze a complete set of 1,373 metabolites, which permits us to carry out unbiased analyses, the move injection platform is unable to discriminate and establish isomers. Our research recognized probably essential metabolic processes for BCG-induced educated immunity, which ought to be validated by focused metabolomics.
In conclusion, the findings from this research affirm the position of circulating metabolites on BCG-induced educated immunity in vivo. As well as, we recognized taurine metabolism as a possible pathway linking metabolism and epigenetic alterations in educated immunity. These information counsel that metabolite supplementation or metabolic modulation could possibly be used sooner or later to enhance vaccine responses and forestall infectious illnesses.
Strategies
Examine design
To check the innate and adaptive immune responses upon BCG vaccination, wholesome volunteers of Western European ancestry had been included within the 300BCG cohort between April 2017 and June 2018, as beforehand described [13]. The research was authorized by the Arnhem-Nijmegen Medical Moral Committee (NL58553.091.16) and carried out in accordance with the declaration of Helsinki. After written knowledgeable consent was obtained, an ordinary dose of 0.1 mL BCG (BCG-Bulgaria, InterVax) was administered intradermally within the left higher arm. EDTA blood was collected earlier than and a couple of weeks and three months after BCG vaccination. Volunteers weren’t eligible for participation in case of systemic medicine use apart from oral contraceptives or acetaminophen, use of antibiotics 3 months earlier than inclusion, earlier BCG vaccination, any febrile sickness 4 weeks earlier than participation, any vaccination 3 months earlier than participation, or a medical historical past of immunodeficiency. Contributors had been additionally not eligible in the event that they recalled a historical past of tuberculosis or had lived in a tuberculosis-endemic setting. No tuberculosis pores and skin check or IFN-γ launch assay was carried out to check for latent tuberculosis an infection, however the tuberculosis incidence within the indigenous Dutch inhabitants is extraordinarily low (1/100,000) [31]. After BCG administration, the research individuals acquired an intensive on-line questionnaire about way of life, weight loss plan, and illness historical past. This questionnaire included questions on meat and fish consumption, which was categorized into every day, weekly, month-to-month, lower than as soon as a month, or by no means, and might be present in S7 Desk.
Peripheral blood mononuclear (PBMC) cell isolation and stimulation
PBMCs had been remoted from entire blood utilizing Ficoll-Paque (GE healthcare) density gradient separation. PBMCs had been washed twice with phosphate-buffered saline (PBS) and resuspended in Dutch modified RPMI 1640 medium (Roswell Park Memorial Institute, Invitrogen), supplemented with 50 μg/mL gentamicin, 2 mM Glutamax (GIBCO), and 1 mM pyruvate (GIBCO). In a complete quantity of 200 μL/properly, 5 × 105 PBMCs had been cultured in spherical backside 96-well plates (Greiner) and stimulated with RPMI (medium management), heat-killed M. tuberculosis H37Rv (5 μg/mL), or heat-killed S. aureus (106 CFU/mL) and incubated at 37°C with 5% CO2. Supernatants had been collected after 24 hours and seven days and saved at −20°C till additional evaluation. IL-1β, IL-6, and TNF-α (ELISA, R&D Techniques) concentrations had been measured after 24 hours, and IFN-γ (Luminex, Thermo Fisher) was measured after 7 days, in response to the producer’s protocols.
Metabolomics measurement
Plasma metabolite ranges had been measured earlier than BCG vaccination was administered. The metabolic options had been measured and annotated by the Basic Metabolics (Zurich, Switzerland) utilizing move injection time-of-flight mass (flow-injection TOF-M) spectrometry, as beforehand described [32]. Non-targeted metabolites had been annotated in response to human metabolites database (HMDB). All metabolomic measurements had been carried out in duplicate, and the common worth was calculated for every pattern per characteristic.
In vitro educated immunity assay
The educated immunity assay was carried out as beforehand described [19,33]. Briefly, buffy coats from wholesome donors had been obtained after written knowledgeable consent (Sanquin Blood Financial institution, Nijmegen, the Netherlands). PBMCs had been remoted utilizing Ficoll-Paque, after which Percoll isolation was carried out to counterpoint the monocyte fraction. The cells had been allowed to stick for 1 hour, after which the non-adherent cells had been washed away with PBS. Then, the adherent monocytes had been preincubated with numerous concentrations of taurine (0.1, 0.5, or 1 mM; imply focus of taurine in blood in adults is 93 μM [34]), methyl-fumarate (0.1 mM) or left untreated, after which they had been incubated for twenty-four hours in medium (unfavourable management), BCG (SSI, 5 μg/mL), or β-glucan (1 μg/mL). Cell viability was assessed utilizing CytoTox 96 non-radioactive cytotoxicity assay (Promega), which measures the cytosolic enzyme lactate dehydrogenase (LDH) that’s launched upon cell lysis. After 24 hours, the cells had been washed with PBS and incubated in medium for five days. Subsequently, the cells had been restimulated with medium alone or ultra-pure Escherichia coli LPS (10 ng/mL), after which the concentrations of IL-6 and TNF-α had been decided within the supernatants (ELISA, R&D Techniques).
Knowledge analyses
All computational analyses had been carried out utilizing R model 3.3.3. The Shapiro–Wilk check was utilized to verify for a traditional distribution of cytokine and metabolite ranges. Preprocessing evaluation consists of log transformation of the cytokine and metabolite information and filtering of metabolites with an annotation rating <70. We studied the impact of age, intercourse, BMI, and oral contraceptive utilization on metabolite ranges utilizing a linear regression mannequin, and a FDR was calculated to regulate for a number of testing. To scale back information dimensionality, we carried out a PCA on baseline metabolites utilizing centered and scaled information. The PCs capturing not less than 75% of the variance had been included within the evaluation to symbolize the metabolome information. Spearman correlation evaluation was carried out on metabolite PCs and ex vivo cytokine fold adjustments, and unsupervised hierarchical clustering was carried out, which was visualized in a heatmap (Fig 1B). Ex vivo cytokine fold adjustments had been related to metabolites utilizing a linear regression mannequin with age and intercourse as covariates (Fig 2).
Subsequent, we assessed the predictive worth of baseline metabolites to foretell the educated immunity responses (ex vivo cytokine fold adjustments in response to S. aureus). Prediction of the cytokine fold adjustments was carried out by coaching an Elastic Internet mannequin, which is a regularization and variable choice methodology particularly helpful when the variety of predictors is bigger than the variety of observations [17]. First, the dataset was randomly break up right into a coaching set (70% of the info) and a check set (30%), after which a 10-fold cross-validation was carried out to pick options, which was repeated 5 instances, based mostly on the coaching set. Prediction accuracy was evaluated by calculation of a Pearson correlation between the measured cytokine fold adjustments and the predictions of the Elastic Internet mannequin on the check units. Then, this process was repeated 100 instances, and the outcomes had been summarized (Fig 3A–3D).
Subsequently, we estimated the variance defined by metabolites, age, intercourse, and BMI on cytokine fold adjustments. To evaluate the contribution of every sort of knowledge upon the others, the total mannequin (included age, intercourse, BMI, and high 10 metabolites) was match first. Subsequently, a number of lowered fashions had been match by which 1 information stage was lacking. The distinction of R2 between the lowered mannequin and the total mannequin was taken as a measure of the variance defined by that stage when accounting for the results of the opposite ranges (Fig 3E) [18].
To grasp the position of metabolite networks in excessive and low responders when it comes to their educated immunity response, we first outlined the excessive and low responders by rating topics based mostly on their IL-1β fold change in response to S. aureus. The highest 30% topics had been outlined as excessive responders and the underside 30% as low responders (Fig 4A). Metabolite information from each excessive and low responders had been merged collectively, and we carried out weighted correlation community evaluation (WGCNA) [35] to establish modules of extremely correlated metabolites (Fig 4B). Then, we correlated the abstract profile (eigengene) for every module with the scientific traits (excessive/low responders, age, intercourse, BMI; Fig 4C). If we recognized any modules that correlated with excessive/low responders, we then extracted the metabolites from this module and carried out pathway evaluation by MetaboAnalyst 5.0 (https://www.metaboanalyst.ca/; KEGG model October 2019) [36], and metabolite pathways with p < 0.05 had been reported (Fig 4D and 4E).
Lastly, all metabolites with a p < 0.05 from the linear regression fashions had been included in additional pathway evaluation utilizing MetaboAnalyst (Fig 5A) [36]. Meat and fish consumption was reported by research individuals as every day, weekly, month-to-month, or by no means, and the ex vivo cytokine fold adjustments between these teams had been in contrast utilizing Kruskal–Wallis assessments.
Supporting info
S4 Fig. Affiliation between malate at baseline and fold adjustments of ex vivo PBMC-derived S. aureus-induced TNF-α responses (left) and between succinate at baseline and fold adjustments of ex vivo PBMC-derived S. aureus-induced IL-1β responses (proper) 3 months after vaccination.
The blue line and shadow point out the linear regression mannequin and the purple line and shadow point out the spline regression mannequin. The cytokine and metabolome information used to generate this determine can be found at https://gitlab.com/xavier-lab-computation/public/bcg300.
https://doi.org/10.1371/journal.pbio.3001765.s004
(TIF)
S5 Fig. Human major monocytes had been incubated for twenty-four hours with fumarate (0.1 mM) or BCG (1 μg/mL), after which the medium was refreshed, and the cells had been allowed to relaxation for five days, after which they had been stimulated with E. coli LPS (10 ng/mL) for twenty-four hours.
Then, the degrees of IL-6 and TNF-α had been measured, and a fold change was calculated relative to the medium management. The median values are offered (N = 6 [IL-6] and N = 3 [TNF-α], Wilcoxon matched-pairs signed rank check, * p < 0.05). The cytokine values used to generate this determine might be present in S1 Knowledge.
https://doi.org/10.1371/journal.pbio.3001765.s005
(TIF)
S7 Fig. Human major monocytes had been preincubated within the presence of taurine (0.1, 0.5, or 1 mM) for 1 hour, after which they had been incubated for twenty-four hours with tradition medium, BCG (5 μg/mL), or β-glucan (1 μg/mL).
After 24 hours, the medium was refreshed, and the cells had been allowed to relaxation for five days, after which they had been stimulated with E. coli LPS (10 ng/mL) for twenty-four hours. Then, the degrees of IL-6 and TNF-α had been measured, and a fold change was calculated relative to the medium management. The median values are offered, and every donor is represented in a unique shade (N = 6, Wilcoxon matched-pairs signed rank check, ns = not important, * p < 0.05). The cytokine values used to generate this determine might be present in S3 Knowledge.
https://doi.org/10.1371/journal.pbio.3001765.s007
(TIF)
S8 Fig. Cytotoxicity assay that measures the cytosolic enzyme LDH that’s launched upon cell lysis to check the cell viability after priming.
Human major monocytes had been preincubated within the presence of taurine (0.1, 0.5, or 1 mM) for 1 hour, after which they had been incubated for twenty-four hours with tradition medium, BCG (5 μg/mL), or β-glucan (1 μg/mL). After 24 hours, LDH was measured within the supernatant (N = 6 per situation). The y-axis represents the diploma of cytotoxicity in share relative to the constructive and unfavourable management. The cytotoxicity values used to generate this determine might be present in S4 Knowledge.
https://doi.org/10.1371/journal.pbio.3001765.s008
(TIF)
Acknowledgments
We thank Julia Brake for her assist with the in vitro educated immunity experiments.
References
- 1.
Netea MG, Dominguez-Andres J, Barreiro LB, Chavakis T, Divangahi M, Fuchs E, et al. Defining educated immunity and its position in well being and illness. Nat Rev Immunol. 2020;20(6):375–388. pmid:32132681 - 2.
Aaby P, Roth A, Ravn H, Napirna BM, Rodrigues A, Lisse IM, et al. Randomized trial of BCG vaccination at start to low-birth-weight youngsters: helpful nonspecific results within the neonatal interval? J Infect Dis. 2011;204(2):245–252. pmid:21673035 - 3.
Biering-Sorensen S, Aaby P, Lund N, Monteiro I, Jensen KJ, Eriksen HB, et al. Early BCG-Denmark and Neonatal Mortality Amongst Infants Weighing <2500 g: A Randomized Managed Trial. Clin Infect Dis. 2017;65(7):1183–1190. - 4.
Kleinnijenhuis J, Quintin J, Preijers F, Joosten LA, Ifrim DC, Saeed S, et al. Bacille Calmette-Guerin induces NOD2-dependent nonspecific safety from reinfection through epigenetic reprogramming of monocytes. Proc Natl Acad Sci U S A. 2012;109(43):17537–17542. pmid:22988082 - 5.
Kleinnijenhuis J, Quintin J, Preijers F, Joosten LA, Jacobs C, Xavier RJ, et al. BCG-induced educated immunity in NK cells: Position for non-specific safety to an infection. Clin Immunol. 2014;155(2):213–219. pmid:25451159 - 6.
Fanucchi S, Fok ET, Dalla E, Shibayama Y, Borner Okay, Chang EY, et al. Immune genes are primed for sturdy transcription by proximal lengthy noncoding RNAs positioned in nuclear compartments. Nat Genet. 2019;51(1):138–150. pmid:30531872 - 7.
van der Heijden C, Noz MP, Joosten LAB, Netea MG, Riksen NP, Keating ST. Epigenetics and Skilled Immunity. Antioxid Redox Sign. 2018;29(11):1023–1040. pmid:28978221 - 8.
Arts RJW, Carvalho A, La Rocca C, Palma C, Rodrigues F, Silvestre R, et al. Immunometabolic Pathways in BCG-Induced Skilled Immunity. Cell Rep. 2016;17(10):2562–2571. pmid:27926861 - 9.
Arts RJ, Novakovic B, Ter Horst R, Carvalho A, Bekkering S, Lachmandas E, et al. Glutaminolysis and Fumarate Accumulation Combine Immunometabolic and Epigenetic Applications in Skilled Immunity. Cell Metab. 2016;24(6):807–819. pmid:27866838 - 10.
Arts RJW, Moorlag S, Novakovic B, Li Y, Wang SY, Oosting M, et al. BCG Vaccination Protects in opposition to Experimental Viral An infection in People by the Induction of Cytokines Related to Skilled Immunity. Cell Host Microbe. 2018;23(1):89–100 e5. pmid:29324233 - 11.
Stroll J, de Bree LCJ, Graumans W, Stoter R, van Gemert GJ, van de Vegte-Bolmer M, et al. Outcomes of managed human malaria an infection after BCG vaccination. Nat Commun. 2019;10(1):874. pmid:30787276 - 12.
Mangtani P, Abubakar I, Ariti C, Beynon R, Pimpin L, High-quality PE, et al. Safety by BCG vaccine in opposition to tuberculosis: a scientific overview of randomized managed trials. Clin Infect Dis. 2014;58(4):470–480. pmid:24336911 - 13.
Koeken VA, de Bree LCJ, Mourits VP, Moorlag SJ, Stroll J, Cirovic B, et al. BCG vaccination in people inhibits systemic irritation in a sex-dependent method. J Clin Make investments. 2020;130(10):5591–5602. pmid:32692728 - 14.
de Bree LCJ, Mourits VP, Koeken VA, Moorlag SJ, Janssen R, Folkman L, et al. Circadian rhythm influences induction of educated immunity by BCG vaccination. J Clin Make investments. 2020;130(10):5603–5617. pmid:32692732 - 15.
Kochlik B, Gerbracht C, Grune T, Weber D. The Affect of Dietary Habits and Meat Consumption on Plasma 3-Methylhistidine-A Potential Marker for Muscle Protein Turnover. Mol Nutr Meals Res. 2018;62(9):e1701062. pmid:29573154 - 16.
Ajioka RS, Phillips JD, Kushner JP. Biosynthesis of heme in mammals. Biochim Biophys Acta. 2006;1763(7):723–736. pmid:16839620 - 17.
Zou H, Hastie T. Regularization and variable choice through the elastic web (vol B 67, pg 301, 2005). J R Stat Soc B. 2005(67):768. - 18.
Bakker OB, Aguirre-Gamboa R, Sanna S, Oosting M, Smeekens SP, Jaeger M, et al. Integration of multi-omics information and deep phenotyping permits prediction of cytokine responses. Nat Immunol. 2018;19(7):776–786. pmid:29784908 - 19.
Dominguez-Andres J, Arts RJW, Bekkering S, Bahrar H, Blok BA, de Bree LCJ, et al. In vitro induction of educated immunity in adherent human monocytes. STAR Protoc. 2021;2(1):100365. pmid:33718890 - 20.
Lu Y, Zou L, Su J, Tai ES, Whitton C, Dam RMV, et al. Meat and Seafood Consumption in Relation to Plasma Metabolic Profiles in a Chinese language Inhabitants: A Mixed Untargeted and Focused Metabolomics Examine. Vitamins. 2017;9(7). - 21.
Tannahill GM, Curtis AM, Adamik J, Palsson-McDermott EM, McGettrick AF, Goel G, et al. Succinate is an inflammatory sign that induces IL-1beta by HIF-1alpha. Nature. 2013;496(7444):238–242. - 22.
Minassian AM, Satti I, Poulton ID, Meyer J, Hill AV, McShane H. A human problem mannequin for Mycobacterium tuberculosis utilizing Mycobacterium bovis bacille Calmette-Guerin. J Infect Dis. 2012;205(7):1035–1042. pmid:22396610 - 23.
Marcinkiewicz J, Kontny E. Taurine and inflammatory illnesses. Amino Acids. 2014;46(1):7–20. pmid:22810731 - 24.
Schuller-Levis GB, Park E. Taurine: new implications for an outdated amino acid. FEMS Microbiol Lett. 2003;226(2):195–202. pmid:14553911 - 25.
Schaffer S, Kim HW. Results and Mechanisms of Taurine as a Therapeutic Agent. Biomol Ther (Seoul). 2018;26(3):225–241. pmid:29631391 - 26.
Jong CJ, Azuma J, Schaffer S. Mechanism underlying the antioxidant exercise of taurine: prevention of mitochondrial oxidant manufacturing. Amino Acids. 2012;42(6):2223–2232. pmid:21691752 - 27.
MacIver NJ, Michalek RD, Rathmell JC. Metabolic regulation of T lymphocytes. Annu Rev Immunol. 2013;31:259–283. pmid:23298210 - 28.
Chu X, Jaeger M, Beumer J, Bakker OB, Aguirre-Gamboa R, Oosting M, et al. Integration of metabolomics, genomics, and immune phenotypes reveals the causal roles of metabolites in illness. Genome Biol. 2021;22(1):198. pmid:34229738 - 29.
Etchegaray JP, Mostoslavsky R. Interaction between Metabolism and Epigenetics: A Nuclear Adaptation to Environmental Adjustments. Mol Cell. 2016;62(5):695–711. pmid:27259202 - 30.
Setiabudiawan TP, Reurink RK, Hill PC, Netea MG, van Crevel R, Koeken VACM. Safety in opposition to tuberculosis by Bacillus Calmette-Guérin (BCG) vaccination: A historic perspective. Medicina. 2022;3(1):6–24. pmid:35590145 - 31.
van den Boogaard J, Hunch E, Erkens C, Meijer W, Schimmel H, van Soolingen D, et al. Op weg naar eliminatie van tuberculose in Nederland. Infectieziekten. Dengue Bull. 2020;31(2). - 32.
Fuhrer T, Heer D, Begemann B, Zamboni N. Excessive-throughput, correct mass metabolome profiling of mobile extracts by move injection-time-of-flight mass spectrometry. Anal Chem. 2011;83(18):7074–7080. pmid:21830798 - 33.
Bekkering S, Blok BA, Joosten LA, Riksen NP, van Crevel R, Netea MG. In Vitro Experimental Mannequin of Skilled Innate Immunity in Human Main Monocytes. Clin Vaccine Immunol. 2016;23(12):926–933. pmid:27733422 - 34.
Trabado S, Al-Salameh A, Croixmarie V, Masson P, Corruble E, Feve B, et al. The human plasma-metabolome: Reference values in 800 French wholesome volunteers; impression of ldl cholesterol, gender and age. PLoS ONE. 2017;12(3):e0173615. pmid:28278231 - 35.
Langfelder P, Horvath S. WGCNA: an R bundle for weighted correlation community evaluation. BMC Bioinformatics. 2008;9:559. pmid:19114008 - 36.
Chong J, Wishart DS, Xia J. Utilizing MetaboAnalyst 4.0 for Complete and Integrative Metabolomics Knowledge Evaluation. Curr Protoc Bioinformatics. 2019;68(1):e86. pmid:31756036
