- Research
- Open access
- Published:
Metabolic and inflammatory perturbation of diabetes associated gut dysbiosis in people living with and without HIV infection
Genome Medicine volume 16, Article number: 59 (2024)
Abstract
Background
Gut dysbiosis has been linked with both HIV infection and diabetes, but its interplay with metabolic and inflammatory responses in diabetes, particularly in the context of HIV infection, remains unclear.
Methods
We first conducted a cross-sectional association analysis to characterize the gut microbial, circulating metabolite, and immune/inflammatory protein features associated with diabetes in up to 493 women (~ 146 with prevalent diabetes with 69.9% HIV +) of the Women’s Interagency HIV Study. Prospective analyses were then conducted to determine associations of identified metabolites with incident diabetes over 12 years of follow-up in 694 participants (391 women from WIHS and 303 men from the Multicenter AIDS Cohort Study; 166 incident cases were recorded) with and without HIV infection. Mediation analyses were conducted to explore whether gut bacteria–diabetes associations are explained by altered metabolites and proteins.
Results
Seven gut bacterial genera were identified to be associated with diabetes (FDR-q < 0.1), with positive associations for Shigella, Escherichia, Megasphaera, and Lactobacillus, and inverse associations for Adlercreutzia, Ruminococcus, and Intestinibacter. Importantly, the associations of most species, especially Adlercreutzia and Ruminococcus, were largely independent of antidiabetic medications use. Meanwhile, 18 proteins and 76 metabolites, including 3 microbially derived metabolites (trimethylamine N-oxide, phenylacetylglutamine (PAGln), imidazolepropionic acid (IMP)), 50 lipids (e.g., diradylglycerols (DGs) and triradylglycerols (TGs)) and 23 non-lipid metabolites, were associated with diabetes (FDR-q < 0.1), with the majority showing positive associations and more than half of them (59/76) associated with incident diabetes. In mediation analyses, several proteins, especially interleukin-18 receptor 1 and osteoprotegerin, IMP and PAGln partially mediate the observed bacterial genera–diabetes associations, particularly for those of Adlercreutzia and Escherichia. Many diabetes-associated metabolites and proteins were altered in HIV, but no effect modification on their associations with diabetes was observed by HIV.
Conclusion
Among individuals with and without HIV, multiple gut bacterial genera, blood metabolites, and proinflammatory proteins were associated with diabetes. The observed mediated effects by metabolites and proteins in genera–diabetes associations highlighted the potential involvement of inflammatory and metabolic perturbations in the link between gut dysbiosis and diabetes in the context of HIV infection.
Background
Diabetes, mainly type 2 diabetes, is a prevailing metabolic disorder, posing a significant public health challenge worldwide. People living with human immunodeficiency virus (HIV) infection (PLWH) are at high risk of metabolic diseases, including diabetes, possibly due to the chronic inflammation and persistent immune activation induced by HIV infection and the long-term antiretroviral therapy (ART) [1]. Characterizing modifiable risk factors of diabetes beyond traditional risk factors is thus essential to the diabetes prevention and management in this susceptible population.
Emerging evidence suggests that gut dysbiosis may play a critical role in the development of diabetes [2,3,4]. Both our pilot study and other prior work have reported altered microbiota composition associated with some metabolic disorders, such as metabolic syndrome [5], prediabetes [6], and diabetes [7] among PLWH. However, these previous studies were limited by relatively small sample sizes, and characteristics of gut dysbiosis associated with diabetes in the context of HIV infection are not fully understood. How and through which molecular mechanisms gut dysbiosis may contribute to diabetes in HIV infection remains largely unclear. Evidence from experimental settings has suggested that gut microbiota may affect the host metabolic health through direct modulation of host immunity/inflammation responses [2, 4, 8], as well as by producing microbial-derived biomolecules especially those play a role in modulation of host metabolism, inflammation, and immune activation (e.g., indoles, short-chain fatty acids [SCFAs], secondary bile acids) [9, 10]. Given that alterations in gut microbial taxonomic composition (e.g., the loss of beneficial microbes and enriched pathobionts) and related microbial functionalities, as well as disrupted inflammation and immune activation homeostasis, have been noted in HIV infection [11,12,13,14], we thus hypothesized that among PLWH or population at high risk for HIV infection, gut dysbiosis is associated with diabetes, which is partially explained by metabolic-inflammatory perturbations.
To test this hypothesis, we leveraged multi-omics data from the Women’s Interagency HIV Study (WIHS) and the Multicenter AIDS Cohort Study (MACS), two well-characterized cohorts of both PLWH and demographically and socioeconomically similar people without HIV infection in US, to first identify gut microbial taxonomic features, circulating metabolites, immune and inflammatory proteins associated with diabetes, and then to explore the potential mediating roles of metabolites and inflammatory proteins in the observed gut microbiota–diabetes relationship.
Methods
Study design and population
The WIHS was a multicenter longitudinal study started in 1993 and designed to investigate the long-term, natural, treated history, and progression of HIV infection in women [15]. The MACS was a prospective study originally founded in 1984 and focused on the natural and treated histories of HIV infection in men [16]. The WIHS and MACS used similar study designs and methods, recruited both PLWH and HIV-negative people with similar demographical and socioeconomical status and high-risk behaviors for HIV infection, and were merged to form the MACS-WIHS Combined Cohort Study in the beginning of 2019 [17]. In the present investigation, we included 563 women who provided stool samples from three sites (Bronx, Brooklyn, and Chicago) during core visits at 6-month intervals from 2017 to 2019 in WIHS [18, 19]. Among them, 493 women had shotgun metagenomics profiled, 434 had plasma metabolomics profiled, and 428 had serum immune/inflammatory proteins. Among 426 participants with all these three omics data, 396 (~ 93%) had the omics data measured on biospecimens collected at the same visits and 12 had data matched within 1 year (i.e., ≤ 2 visits). To examine the prospective relationship between plasma metabolites and incident diabetes, we included 694 participants (391 women from WIHS and 303 men from MACS; 166 incident cases were recorded over a median follow-up of 12.6 years) who were free of diabetes and with metabolomic and lipidomic data at a baseline visit (2004–2006) [20, 21]. An overview of study design and sample selection of the present study is shown in Fig S1.
Gut metagenomics profiling
In this study, only participants from WIHS provided fecal samples for gut metagenomics profiling. Fecal samples were collected using a home-based self-collection kit that was distributed to each participant at a core WIHS visits as previously described [18, 19]. Metagenomic sequencing was performed on DNA extracted from fecal samples using a shallow-coverage method of shotgun sequencing-based Illumina NovaSeq platform [19]. De-multiplexing was applied to generate per-sample FASTQ data, which was further trimmed to remove low-quality bases with a 25 or less PHRED quality score using prinseq-lite 0.20.4 (https://edwards.sdsu.edu/cgibin/prinseq/prinseq.cgi). After quality control, FASTQ data were then concatenated and aligned to the NCBI RefSeq representative prokaryotic genome collection using the SHOGUN pipeline (https://github.com/knights-lab/SHOGUN), through which reads were labeled with NCBI taxonomic annotation at the species level.
Plasma metabolomics profiling
Metabolomics profiling was performed on plasma samples using hydrophilic interaction liquid chromatography/positive ion mode (HILIC-pos for water soluble metabolites) and reversed-phase C8 chromatography/positive ion mode mass spectrometry (C8-pos for liposoluble metabolites) at the Broad Institute Metabolomics Platform (Cambridge, Massachusetts) as previously described [21]. In this study, we included 378 metabolites (211 lipids and 167 nonlipids) with detection rate > 75% in cross-sectional analyses of prevalent diabetes involving women from WIHS and 325 metabolites (211 lipids and 114 nonlipids) in the prospective analyses of incident diabetes involving both women from WIHS and men from MACS. Levels of metabolites below detection were replaced by one half of the minimum detected values.
Assessment of inflammatory proteins
A set of inflammation-related proteins were quantified in serum among 428 women from WIHS using the Olink Inflammation panel (92 proteins, Olink Bioscience, Uppsala, Sweden). Data of targeted proteins is presented as normalized protein expression (NPX) values with the unit at log2 scale. We included 74 proteins detected in > 75% of samples and with average coefficient of variance (CV, %) < 10% (intra-assay CV = 4%, inter-assay CV = 9%). Levels of proteins below detection were replaced by one half of the square root of the limit of detection.
Ascertainment of HIV infection, diabetes, and other covariates
Data on sociodemographic, behavioral, and clinical characteristics, including laboratory testing, were collected following standardized protocols [15, 16, 21]. Medication histories including use of antihypertensive, lipid-lowering, and antidiabetic medications were collected with standardized interviewer-administered questionnaires. HIV serostatus was ascertained by enzyme-linked immunosorbent assay and confirmed by Western blot. Other HIV-related parameters included serum CD4 + cell counts, HIV RNA viral load, and current use of ART drugs. Diabetes was defined as fasting plasma glucose ≥ 126 mg/dL, random plasma glucose ≥ 200 mg/dL, HbA1c ≥ 6.5%, or self-report use of antidiabetic medications in line with our previous work [21]. In the prospective analysis, participants who were free of diabetes at baseline but developed diabetes during follow-up were defined as incident diabetes cases.
Statisticalanalysis
The Analysis of Composition of Microbiome-II (ANCOM-II) [22] method that accounts for structural zeros was used to identify gut bacteria associated with prevalent diabetes among 493 women in WIHS. ANCOM-II was conducted at both the genus (n = 97) and species (n = 203) levels with predominant bacterial taxa (relative abundance > 0.01% and prevalence > 10%) included, while adjusting for age, study site, race/ethnicity, household annual income, education, smoking, alcohol consumption, HIV serostatus, and antibiotic use within 4 weeks of stool sample collection. The false discovery rate (FDR) using the Benjamini–Hochberg method was controlled at 0.1, and the threshold of identifying significant microbiota taxa in ANCOM-II was set at 0.6, which refers to ratio of the tested taxon to at least 60% of the other taxa detected to be significantly associated with diabetes. The associations of identified gut bacterial taxa (centered log-ratio transformed, CLR) with prevalent diabetes were further examined using logistic regression with the adjustment of the same set of covariates. Additionally, differences in bacterial α-diversity at the species level (observed, Chao1, Shannon, and Simpson diversity indices) between women with and without diabetes were tested by the Kruskal–Wallis test. Association between β diversity (Bray–Curtis dissimilarity) and diabetes was tested using the permutational multivariate analysis of variance (PERMANOVA, n = 9999 permutations) [23].
The cross-sectional associations of serum proteins (N = 428, inverse normal transformed (INT)) and metabolites (N = 434, INT) with prevalent diabetes in WIHS women were assessed using logistic regression models adjusted for age, study site, race/ethnicity, household annual income, education, smoking, HIV serostatus, and fasting status. The false discovery rate was controlled at 0.1. For the identified proteins, we conducted sparse partial least square discriminant analysis (sPLS-DA) [24] to examine their performance in classifying participants with and without diabetes. In addition, we used Cox proportional hazards models to examine prospective associations of identified metabolites with incident diabetes among 694 women and men who were free of diabetes at baseline, adjusting for sex and the same set of covariates in the aforementioned cross-sectional analysis. Analyses were also conducted in women and men separately and results were then combined by meta-analysis to test potential sex differences in metabolite–incident diabetes associations. To reduce the dimension of the diabetes-associated metabolites, the weighted correlation network analysis (WGCNA) [25] was applied to identify metabolite modules for downstream analysis.
After omics-wide association analyses, we conducted partial correlation analyses, adjusting for the covariates listed above, to assess the interrelationship among identified diabetes-associated gut bacteria, proteins, and metabolites. By using the “cmest” function from the R CMAverse package [26], mediation analyses were conducted to investigate whether associations between gut bacteria and diabetes were partially mediated by alterations in circulating metabolites and proteins. Details about this mediation analysis were provided in the Supplementary Methods.
To test the robustness of our results for omics–diabetes associations, we further adjusted for a set of diabetes-related factors, including body mass index (BMI), triglycerides, systolic blood pressure (SBP), high-density lipoprotein cholesterol (HDL-C), and use of lipid lowering and antihypertensive medications. Given that diabetes treatments (e.g., metformin) have been linked with gut microbiota composition [27,28,29], we further conducted a “drug deconfounding” analysis [28, 29] using a recently introduced two-step pipeline to examine whether the identified gut microbiota–diabetes associations were confounded by antidiabetic medications. A brief introduction to this “drug deconfounding” analysis is provided in the Supplementary Methods. Stratified analyses on the associations between the identified omics features and diabetes were also conducted by HIV serostatus. To test potential effect modification by HIV serostatus, we included a product term with omics features in the regression models. Differences in the levels of identified diabetes-associated omics features between women with and without HIV infection were also compared using analysis of covariance (ANCOVA). To further integrate HIV infection, gut dysbiosis, and metabolic and inflammatory perturbation with diabetes, we conducted in parallel omics wide association analyses to characterize omics features associated with HIV infection and evaluate their associations with diabetes.
In the above association analyses, all omics features were standardized into Z-scores after the original transformation (CLR for microbial features and INT for metabolites/proteins) to improve the comparability of each feature in association with diabetes. All analyses were conducted using R 3.4. Statistical significance was set at 0.05 (two-tailed) unless otherwise stated.
Results
Population characteristics
As shown in Table S1, up to 493 women, including 336 living with HIV, with a mean age of 53.2 years, were included in the cross-sectional, omics-wide association analysis of prevalent diabetes. Approximately 91% of women living with HIV reported receiving ART drugs within 6 months at the time of stool samples collection and nearly 73.5% had undetectable HIV viral load. Compared to women without diabetes, women with diabetes tended to be older and had an unfavorable profile of cardiometabolic traits and higher proportion of antihypertensive and lipid lowering medication use.
Gut microbiota and diabetes
We found a significant difference in β-diversity (PMERANOVA, P = 0.01) but not in α-diversity between women with and without diabetes (Fig. 1A). ANCOM-II at the genus level identified 7 out of 97 predominant bacteria genera associated with diabetes (FDR-q < 0.1), including 4 Firmicutes genera (i.e., Intestinibacter, Ruminococcus, Megasphaera, and Lactobacillus), 2 Proteobacteria genera (i.e., Shigella and Escherichia), and 1 Actinobacteria genus (Adlercreutzia) (Fig. 1B and Table S2). Specifically, Adlercreutzia, Intestinibacter, and Ruminococcus were inversely associated with diabetes, whereas Shigella, Escherichia, Megasphaera, and Lactobacillus exhibited positive associations with diabetes (Fig. 1C). We identified 20 common species (relative abundance > 0.01% and prevalence > 5%) within these 7 genera and found that most species within the same genus showed generally consistent associations with diabetes (Fig. 1C). Consistently, species within the same genus were positively correlated with each other (Fig S2). We also found strong positive correlations between members of Shigella and Escherichia, the two bacteria in the Enterobacteriaceae family, at both the genus and species levels. ANCOM-II at the species levels only identified one additional species that was associated with diabetes, Eubacterium_eligens, which does not belong to any of these 7 genera (Fig S3). Given these observations, we thus focused on the identified diabetes-associated genera in the following analyses.
Further adjustment of major diabetes-related metabolic traits (BMI, triglycerides, HDL, and SBP) and use of antihypertensive and lipid-lowering medications did not substantially change the associations of these bacterial genera with diabetes (Table S3). The “drug deconfounding” analysis suggested that associations of Intestinibacter and Escherichia with diabetes might be confounded by antidiabetic medication use (Fig S4A). For example, women with treated diabetes had significantly higher abundance of Escherichia compared to women without diabetes and those with untreated diabetes, whereas abundance of this genus was similar between women without diabetes and those with untreated diabetes (Fig S4B). A similar pattern was found for Shigella (Fig S4B).
Serum inflammatory proteins and diabetes
After multivariable adjustment (see “ Methods”), we found 18 out of 74 proteins associated with prevalent diabetes (FDR-q < 0.1), with inverse associations for T cell surface glycoprotein CD6 isoform (CD6), C-X-C motif chemokine 10 (CXCL10), and tumor necrosis factor (ligand) superfamily, member 12 (TWEAK), and positive associations for the remaining 15 proteins (Fig.2A). Interleukin-18 receptor 1 (IL-18R1) showed the strongest positive association with diabetes (OR = 2.09 [95%CI: 1.61, 2.73] per SD increase), while TWEAK showed the strongest inverse association with diabetes (OR = 0.70 [95%CI: 0.56, 0.89]). Associations of these 18 proteins with diabetes did not change appreciably when we adjusted for diabetes-related metabolic traits and use of anti-hypertensive and lipid-lowering medications (Table S4). Using sPLS-DA, we found that the first two components of these identified 18 proteins generally had a good performance in diabetes classification (AUC = 0.77, P < 0.001).
Partial correlation analysis revealed a cluster of seven proteins (hepatocyte growth factor (HGF), macrophage colony-stimulating factor 1 (CSF-1), osteoprotegerin (OPG), IL-18R1, vascular endothelial growth factor A (VEGFA), glial cell line-derived neurotrophic factor (GDNF), and IL6) that were highly intercorrelated (Fig. 2C) and all positively associated with diabetes (Fig. 2A). These clustered 7 proteins also showed some weak-to-moderate positive correlations with the remaining 11 proteins, including 3 proteins that were inversely associated with diabetes (CD6, CXCL10, and TWEAK). However, correlation among these 11 remaining proteins were comparatively weak.
Plasma metabolites and diabetes
We first conducted a cross-sectional metabolome wide association analysis in 434 WIHS women and identified 132 (FDR-q < 0.1) out of 378 metabolites associated with diabetes (Fig. 3A). Among these 132 metabolites, 115 were available in a prospective dataset of women and men in the WIHS and MACS (N = 694, including 166 incident diabetes cases over a median follow-up of 12.6 years). We then conducted a prospective analysis and found that 59 (50 lipid species and 9 nonlipid metabolites) out of 115 metabolites were associated with incident diabetes (FDR-q < 0.10) (Fig. 3A and B). Separate analyses in women (WIHS) and men (MACS) revealed that the associations of the majority of these 59 metabolites with incident diabetes were comparable between sexes (cohorts) (Fig S5), and combined results from women and men by meta-analysis were very similar to those from the analyses in the overall sample (Table S5). Together with the 17 metabolites that were not available in the prospective dataset, a total of 76 metabolites, including 50 lipids and 26 nonlipid metabolites, were considered as diabetes-associated metabolites in the present study (Fig. 3A).
In line with our previous findings [21], we found that lipid species belonging to diradylglycerols (DGs, n = 12) and triradylglycerols (TGs, n = 20), together with two carnitines (CAR 5:0 and CAR 4:0), were positively associated with prevalent diabetes, whereas cholesterol esters (CEs, n = 5) and phospholipids belonging to phosphatidylethanolamines (PEs, n = 5), phosphatidylcholines (PCs, n = 7), and one lyso-phosphatidylcholines (LPC 18:2) exhibited inverse associations with diabetes (Fig. 3C and Table S4). Correlation analysis revealed that these 50 lipids were highly intercorrelated (Fig S6) and clustered into 4 modules in WGCNA models (Fig S7), which were then included in the following analyses to reduce the dimensionality of diabetes-associated lipids. As expected, lipid modules characterized by high levels of DGs and/or TGs (i.e., lipid modules 1, 3, and 4) were positively associated with diabetes, whereas a lipid module characterized by high levels of PCs and PEs (lipid module 2) was inversely associated with diabetes (Fig. 4 and Table S4). Regarding the 26 nonlipid metabolites, 3 metabolites (glycine, homoarginine, and sphinganine) were inversely associated with diabetes, whereas positive associations with diabetes were found for the remaining 18 metabolites, including 3 microbial metabolites (trimethylamine N-oxide (TMAO) [30], phenylacetylglutamine (PAGln) [31], and imidazolepropionic acid (IMP) [32]), 3 branched-chain amino acids (BCAAs: valine, leucine, and isoleucine) and 1 related metabolite (N-Acetylleucine), 8 non-essential amino acids and related metabolites (e.g., glutamic acid and alanine), and 6 xenobiotics (e.g., xanthopterin). After further adjusting for diabetes-related metabolic traits, the associations of these metabolites/lipid modules with prevalent diabetes were largely attenuated, especially for lipids and lipid modules, though most of the associations remained statistically significant (Table S4).
Gut microbiota, circulating metabolites, proteins, and diabetes
Partial correlation analysis revealed that the identified diabetes-associated bacteria genera were correlated with multiple diabetes-associated metabolites, lipid modules and proteins (FDR-q < 0.1) (Fig. 4 and Fig S8). Overall, Ruminococcus and Adlercreutzia, two genera inversely associated with diabetes, tended to be inversely correlated with metabolites and lipids modules and proteins that were positively associated with diabetes. By contrast, an overall opposite pattern with respect to the correlations with metabolites, lipid modules and proteins were found for the 4 genera (Megasphaera, Lactobacillus, Shigella, and Escherichia) that were positively associated with diabetes (Fig. 4). In particular, 6 out of 7 diabetes-associated genera (Intestinibacter, Megasphaera, Shigella, Escherichia, Ruminococcus, and Adlercreutzia) were correlated with IMP (FDR-q < 0.1), with inverse associations for Ruminococcus and Adlercreutzia and positive association for other 4 bacterial genera. In addition, Ruminococcus and Intestinibacter were inversely correlated with PAGln (FDR-q < 0.1).
To test whether the associations between identified bacterial genera and diabetes can be partially explained by alterations in the identified diabetes-associated metabolites and proteins, we conducted conditional analysis and mediation analysis. Herein, we focused on metabolites, lipid modules, and proteins that have concordant associations with both diabetes and bacteria genera (FDR-q < 0.1, Table S 6). In conditional analyses, associations between bacteria genera and diabetes were generally attenuated when selected metabolites and proteins were included in the models (Fig. 5A). Mediation analysis revealed that the selected multiple nonlipid metabolites, lipid modules, and proteins collectively mediated a considerable proportion of the associations of 6 genera (Intestinibacter (proportion of mediated effect for multiple mediators = 24.89%), Megasphaera (27.09%), Ruminococcus (30.96%), Shigella (34.78%), Adlercreutzia (43.52%), and Escherichia (67.12%), all P < 0.01) with diabetes. Notably, associations of Ruminococcus, Adlercreutzia, and Escherichia with diabetes were mediated considerably by multiple inflammatory proteins led by OPG, IL-18R1, and HGF (Fig. 5A and B). IMP appeared to partially mediate the associations of the most identified genera with diabetes (mediated proportion ranges from 7.5 to 22.3%), though some of mediated effects were nominally significant (Fig. 5A). In addition, PAGln significantly mediated 13.9% of the inverse association between Intestinibacter and diabetes (Fig. 5A). The interrelationship among identified diabetes-associated bacteria, proteins, and metabolites were shown in Fig. 5C.
In addition, we also examined correlations of these omics features with glycemic traits and found that metabolites (e.g., alanine and DMGV) or proteins (e.g., OSM, IL-18R1, and HGF) significantly mediated the observed gut bacteria and diabetes associations were generally associated with examined glycemic traits in the same direction of their association with diabetes (Fig S9).
HIV infection and diabetes-associated omics features
Stratified analyses by HIV serostatus showed generally similar associations between identified omics features and diabetes among women with and without HIV (Table S7). We then compared levels of the diabetes-associated omics features between women with and without HIV and found that multiple lipids belonging to DGs/TGs, nonlipids metabolites (e.g., IMP, PAGln, ureidopropionic acid, xanthopterin, 4-aminohippuric acid, N4-acetylcytidine, 4-acetamidobutanoic acid, glycocholic acid, and homocitrulline), and proteins (GDNF and leukemia inhibitory factor receptor (LIF-R)) that were positively associated with diabetes were significantly higher among women with HIV, especially those with detectable HIV RNA viral load (> 20 copies/mL) as compared with women without HIV (Table S8). In addition, several metabolites that were inversely associated with diabetes (e.g., homoarginine and lipids belonging to PCs/PEs/CEs) were lower in women with HIV compared with those without HIV (Table S8). However, abundances of identified diabetes-associated genera were similar among women with and without HIV (Table S8).
Omics-wide associations identified 11 bacterial genera (Fig. 6A), 20 proteins (Fig. 6B), and 30 metabolites (Fig. 6C) significantly associated with HIV infection (Fig. 6 and Table S9, FDR-q < 0.1). Of note, the majority of identified metabolites and bacteria genera had concordant associations with HIV infection and prevalent diabetes (i.e., metabolites positively/inversely associated with HIV infection were also positively/inversely associated with diabetes (Fig S10), including two bacteria (e.g., Lanchoclostridum and Faecalibacterium) and 11 nonlipids metabolites exhibiting significant associations with diabetes (Fig. 6A and C). In addition, we found two proteins significantly associated with both HIV and diabetes in the same direction (Fig. 6B), but many of these identified HIV-associated proteins were not associated with diabetes.
Discussion
In the present study, among PLWH and demographically and socioeconomically similar individuals without HIV infection, we characterized multi-omics features associated with prevalent diabetes, including 7 gut bacterial genera, 18 immune and pro-inflammatory proteins, and 76 metabolites. More than half of metabolites associated with prevalent diabetes were confirmed to be associated with incident diabetes over a median follow-up of 12 years in a prospective cohort. Moreover, several proteins, especially IL-18R1 and OPG, and microbial derived metabolites, including IMP and PAGln, significantly mediated the associations between identified gut bacteria and diabetes, both individually and in combination.
Our findings on the observed altered gut microbiota composition associated with diabetes are generally consistent with findings among general (i.e., non-HIV) populations, including the positive associations of Shigella [2,3,4, 33], Escherichia [2,3,4, 33], Lactobacillus [2,3,4, 33], and Megasphaera [34], and the negative associations of Ruminococcus [2,3,4, 33] and Aldercreutzia [4, 35]. Despite these consistencies, concerns over the potential confounding by antidiabetic treatments, especially metformin use, in characterizing diabetes-related gut microbiota have been raised [27,28,29, 36, 37]. Our results also indicated that associations of Escherichia, Intestinibacter, and possibly Shigella, with diabetes might be partially attributed to the effects of antidiabetic medication use. However, some evidence suggests that these gut bacteria may still play a role in host metabolic health and diseases [33, 38]. For example, Intestinibacter might be beneficial to host lipid and glucose metabolism and intestinal barrier integrity through its production of biomolecules such as branched SCFAs (isobutyrate and isovalerate) [39, 40].
In searching for potential mechanisms underlying the observed association between gut bacteria and diabetes, our results suggested that alterations in several circulating immune and proinflammatory proteins, especially IL-18R1 and OPG, might play a role. Notably, IL-18R1 is the receptor of proinflammatory cytokine IL-18, and overproduction of IL-18 has been noted as a proxy biomarker for overall inflammation and immune activation in HIV infection [41] and metabolic-inflammatory diseases (e.g., diabetes [42] and non-alcoholic fatty liver diseases (NAFLD) [43]). Moreover, hyperactive of IL-18/IL-18R1 signaling in intestinal epithelial cells has been linked to pathologic breakdown of intestinal mucosal barrier and gut dysbiosis [44]. The observed inverse associations of Ruminococcus and Adlercreutzia with IL-18R1 and diabetes thus might be related to the potential anti-inflammatory effects of these two bacteria given their functionalities in generating biomolecules with immunomodulation and inflammation suppressive effects (e.g., SCFAs from Ruminococcus species [45]; equol, an isoflavandiol estrogen, from Adlercreutzia [46]). Regarding OPG, it is known as an osteoclastogensis inhibitory factor, playing a critical role in bone metabolism through inhibiting the receptor activator of nuclear factor kappa-B (RANK)/RANK ligand (RANKL) pathway [47]. The interaction between RANK/RANKL/OPG axis and gut microbiota has been documented in osteoporosis [48,49,50], a common bone disorder coexisting with diabetes [51]. In line with our observations, overproduction of OPG has been associated with diabetes [52, 53] and metabolic dysfunction-associated [54] or non-alcoholic fatty liver disease [55]. This adds support that OPG signaling might be involved in the relationship between gut dysbiosis and metabolic complications beyond bone metabolism.
Our results also revealed that another five proteins (HGF, CD6, fibroblast growth factor 21 (FGF-21), LIF-R, and OSM) might mediate the associations of gut bacteria with diabetes. Of note, the contributions of HGF [56] and OSM [57] to insulin resistance and glucose impairment in obesity and diabetes have reported, and elevated levels of LIF-R have been noted in hepatic steatosis and MAFLD-associated gut dysbiosis ( [54]). However, the relationships between the other two proteins and diabetes remain controversial, especially for FGF21 [58]. While increased circulating FGF-21 levels was linked to elevated risk of diabetes [59], null results were reported in Mendelian Randomization investigations [60]. FGF-21 was report to enhance skeletal muscle glucose uptake [61] and protect lipid disruption and nonalcoholic steatohepatitis [60, 62], even in a clinical trial [63]. These data suggest that the increased circulating FGF21 in diabetes might be a compensatory mechanism to hyperglycemia or hyperlipidemia [61], echoed with observations regarding the elevated FGF21 levels in HIV infection [64]. More studies are required to reveal how these proteins and their interplay with gut microbiota contribute to diabetes and related metabolic complications. In addition, we found serval proteins associated with diabetes, including those previously reported (e.g., VEGFA [65], TWEAK [66], CXCL10 [67], and CXCL6 [68]), which were generally not associated with the identified diabetes-related gut bacteria, indicating that their associations with diabetes might not be related to the detected gut dysbiosis in our study.
Beside proteins, our results also suggested that alterations in metabolites, especially those microbial-derived, such as IMP and PAGln, might partially mediate the observed associations between gut bacteria and diabetes. Indeed, IMP is a well-characterized microbial metabolite generated from histidine catabolism [32] and has been linked to pathogenesis of impaired glucose metabolism and insulin resistance [32, 69]. In line with our results, increased abundances of many members in Ruminococcus and Adlercreutzia genera were associated with reduced IMP level [70], whereas Escherichia, Intestinibacter, and members in Veillonellaceae family (e.g., Megasphaera) were positively associated with IMP levels [70] as these bacteria harbor the urocanate reductases required for IMP production [32]. Similarly, PAGln is a microbiota derived metabolite from phenylalanine catabolism [31], originally being characterized as a risk factor of cardiovascular diseases [31, 71], and also has been linked to impaired glucose metabolism and diabetes [72]. Aligning with our observations, previous studies have demonstrated that the abundances of Ruminococcus and Intestinibacter species were inversely associated with blood PAGln [70]. Of note, the potential proinflammatory properties of IMP and PAGln have been documented in previous studies [31, 32, 69], and consistently, we also observed positive associations between these two microbial metabolites and multiple proinflammatory proteins (e.g., IL6 and IL17-A/C in Fig S 8). These findings further support the potential mediation effects of these metabolites and proinflammatory proteins in combination linking gut dysbiosis and diabetes. Additionally, our mediation analyses indicated that a lipid module characterized by high DGs but low CEs and several metabolites (e.g., alanine, dimethylguanidino valeric acid (DMGV), and xanthopterin) might partially explain the inverse associations of Ruminococcus and Intestinibacter with diabetes, and the positive association of Shigella with diabetes. However, whether and how the gut microbiota modulates the lipids (e.g., diradylglycerols and cholesterol esters) [73] and xenobiotics [74] remain to be further elucidated.
It has been proposed that there might be a triad relationship among HIV infection, gut dysbiosis, and diabetes [7, 75]. Although the identified diabetes-related gut bacteria were not associated with HIV and their associations with diabetes tended to be independent of HIV serostatus in this study, our further omics-wide association analyses identified multiple gut bacterial associated with HIV infection, with concordant associations with diabetes (e.g., Lachnoclostridium was positively associated with both HIV infection and diabetes, while Faecalibactrium was inversely associated with both HIV infection and diabetes). These results corroborate the idea that gut dysbiosis may be a crucial factor contributing to the compromised metabolic health in chronic HIV infection [7, 75]. Similarly, we also identified many metabolites and several proteins (e.g., SLAMF1: signaling lymphocytic activation molecule and IL-15RA: interleukin-15 receptor subunit alpha) associated with HIV infection and diabetes in the same direction. Furthermore, plasma levels of several proteins (e.g., IL-18R1, IL-6, and OPG), lipids (e.g., DGs/TGs), and nonlipid metabolites (e.g., 2-methylguanosine, xanthopterin, N4-acetylcytidine, glycocholic acid, PAGln, and IMP) that were positively associated with diabetes were elevated in PLWH in our and other studies [42, 76,77,78]. These observations collectively support the potential role of chronic inflammation/persistent immune activation and metabolic perturbations in bridging HIV infection with metabolic disorders [1].
Findings from this study have several implications. Modulation of the identified beneficial bacteria, such as Ruminococcus and Adlercreutzia, through supplementation or augmentation, could be explored as a potential intervention for diabetes prevention and management, especially among susceptible populations like PLWH given the observed anti-inflammatory properties of these bacteria in this study and existing literature [46, 79, 80]. As mechanisms of action for some microbial-derived metabolites become clearer, drugs inhibiting downstream signaling activation, such as deactivating or blocking p38gMAPK and mechanistic target of rapamycin complex 1 (mTORC1) signaling for IMP [32] and adrenergic receptors for PAGln [31], could present a novel and effective approach for diabetes prevention and treatment. Additionally, the identification of proteins mediating the connections between gut bacteria and diabetes support the concept of tailoring anti-inflammatory therapy for metabolic diseases, especially in the context of HIV infection [75, 81] [82]. Notably, drugs targeting the IL-18/IL-18R signaling axis, such as IL-18R antagonist or components resembling naturally occurring IL-18 binding protein (IL-18BP) and anti-IL-18 antibodies, are in development and testing for diabetes related conditions, including hypertension and chronic kidney diseases [83,84,85], as well as HIV infection [41]. Similarly, drugs like Denosumab targeting the RANKL/RANK/OPG pathway have been proposed as potential novel options for diabetes treatment and glycemic control [86, 87].
The main strength of this study is that we leveraged multi-omics data in a cohort of PLWH with demographically and socioeconomically comparable HIV seronegative participants. The validation of diabetes-associated metabolites with prospective analyses further solidifies our findings regarding the metabolic alterations in diabetes. Nevertheless, our study has several limitations. The nature of observation study limits causal inference, especially for the cross-sectional results of gut microbiota and proteins with diabetes. Findings need to be confirmed in prospective studies. Although we have carefully addressed the potential impacts of antidiabetic treatment on our results through “drug deconfounding” analysis, potential residual confounding related to varied types of medication use may still exist. Meanwhile, some potential confounding factors such as diet that have been associated with both gut microbiota and diabetes were not well addressed in the present study due to the lack of these data. Furthermore, the presented study was limited in the depth and coverage of omics techniques, especially for metabolites and proteins, which might partially prevent us further examining known molecules and pathways (e.g., gut microbial indole metabolites from tryptophan catabolism [20, 88]) as well as possibly novel ones, that are relevant to HIV infection and metabolic disorders. In addition, the relatively small number of HIV seronegative people limited our power to detect the effect modification by HIV infection. The associations of gut microbiota and serum proteins with diabetes were examined in women and need to be further validated in men and other HIV cohorts.
Conclusions
In the present multi-omics analysis among individuals with and without HIV infection, we identified multiple gut bacteria, including two potentially beneficial genera (Ruminococcus and Adlercreutzia), numerous circulating metabolites, and proinflammatory proteins associated with diabetes. The observed gut bacteria–diabetes associations were partially explained by some microbial metabolites and proinflammatory proteins, supporting the role of gut microbiota in regulating host metabolite profile and inflammation/immune activation which may contribute to the development of diabetes.
Availability of data and materials
Access to individual-level data from the MACS/WIHS Combined Cohort Study Data (MWCCS) may be obtained upon review and requires an approval of a MWCCS concept sheet submission process. The timeline for the review and approval of any individual data applications may be varied but usually applicants will receive comments within 1 month. Links and instructions for online concept sheet submission and details on data sharing plan are on the study website (https://statepi.jhsph.edu/mwccs/work-with-us/).
Abbreviations
- HIV:
-
Human immunodeficiency virus
- PLWH:
-
People living with HIV infection
- ART:
-
Antiretroviral therapy
- SCFAs:
-
Short-chain fatty acids
- BCAAs:
-
Branched-chain amino acids
- COCOMO:
-
The Copenhagen Comorbidity in HIV infection study
- WIHS:
-
The Women’s Interagency HIV Study
- MACS:
-
The Multicenter AIDS Cohort Study
- ANCOM-II:
-
The Analysis of Composition of Microbiome-II
- FDR:
-
False discovery rate
- PERMANOVA:
-
The permutational multivariate analysis of variance
- BMI:
-
Body mass index
- SBP:
-
Systolic blood pressure
- HDL-C:
-
High-density lipoprotein cholesterol
- ANCOVA:
-
Analysis of covariance
- IL-18R1:
-
Interleukin-18 receptor 1
- CD6:
-
T cell surface glycoprotein CD6 isoform
- CXCL10:
-
C-X-C motif chemokine 10
- TWEAK:
-
Tumor necrosis factor (ligand) superfamily, member 12
- HGF:
-
Hepatocyte growth factor
- CSF-1:
-
Macrophage colony-stimulating factor 1
- OPG:
-
Osteoprotegerin
- VEGFA:
-
Vascular endothelial growth factor A
- GDNF:
-
Glial cell line-derived neurotrophic factor
- CARs:
-
Carnitines
- CEs:
-
Cholesterol esters
- PEs:
-
Phosphatidylethanolamines
- PCs:
-
Phosphatidylcholines
- LPCs:
-
Lyso-phosphatidylcholines
- TMAO:
-
Trimethylamine N-oxide
- PAGln:
-
Phenylacetylglutamine
- IMP:
-
Imidazolepropionic acid
- LIF-R:
-
Leukemia inhibitory factor receptor
- RANK:
-
Receptor activator of nuclear factor kappa-B
- FGF-21:
-
Fibroblast growth factor 21
- DMGV:
-
Dimethylguanidino valeric acid
- mTORC1:
-
Mechanistic target of rapamycin complex 1
- SLAMF1:
-
Signaling lymphocytic activation molecule
- IL-15RA:
-
Interleukin-15 receptor subunit alpha
References
Lake JE, Currier JS. Metabolic disease in HIV infection. Lancet Infect Dis. 2013;13:964–75.
Gurung M, Li Z, You H, Rodrigues R, Jump DB, Morgun A, Shulzhenko N. Role of gut microbiota in type 2 diabetes pathophysiology. EBioMedicine. 2020;51:102590.
Aron-Wisnewsky J, Warmbrunn MV, Nieuwdorp M, Clement K. Metabolism and mMetabolic dDisorders and the mMicrobiome: the intestinal microbiota associated with obesity, lipid metabolism, and metabolic health-pathophysiology and therapeutic strategies. Gastroenterology. 2021;160:573–99.
Metwaly A, Reitmeier S, Haller D. Microbiome risk profiles as biomarkers for inflammatory and metabolic disorders. Nat Rev Gastroenterol Hepatol. 2022;19:383–97.
Gelpi M, Vestad B, Hansen SH, Holm K, Drivsholm N, Goetz A, Kirkby NS, Lindegaard B, Lebech AM, Hoel H, et al. Impact of human immunodeficiency virus-related gut microbiota alterations on metabolic comorbid conditions. Clin Infect Dis. 2020;71:e359–67.
Jayanama K, Phuphuakrat A, Pongchaikul P, Prombutara P, Nimitphong H, Reutrakul S, Sungkanuparph S. Association between gut microbiota and prediabetes in people living with HIV. Curr Res Microb Sci. 2022;3:100143.
Moon JY, Zolnik CP, Wang Z, Qiu Y, Usyk M, Wang T, Kizer JR, Landay AL, Kurland IJ, Anastos K, et al. Gut microbiota and plasma metabolites associated with diabetes in women with, or at high risk for. HIV infection EBioMedicine. 2018;37:392–400.
Schirmer M, Smeekens SP, Vlamakis H, Jaeger M, Oosting M, Franzosa EA, Ter Horst R, Jansen T, Jacobs L, Bonder MJ, et al. Linking the human gut microbiome to inflammatory cytokine production capacity. Cell. 2016;167(1125–1136):e1128.
Agus A, Clement K, Sokol H. Gut microbiota-derived metabolites as central regulators in metabolic disorders. Gut. 2021;70:1174–82.
Fan Y, Pedersen O. Gut microbiota in human metabolic health and disease. Nat Rev Microbiol. 2021;19:55–71.
Appay V, Sauce D. Immune activation and inflammation in HIV-1 infection: causes and consequences. J Pathol. 2008;214:231–41.
Bandera A, De Benedetto I, Bozzi G, Gori A. Altered gut microbiome composition in HIV infection: causes, effects and potential intervention. Curr Opin HIV AIDS. 2018;13:73–80.
McHardy IH, Li X, Tong M, Ruegger P, Jacobs J, Borneman J, Anton P, Braun J. HIV Infection is associated with compositional and functional shifts in the rectal mucosal microbiota. Microbiome. 2013;1:1–12.
Vujkovic-Cvijin I, Somsouk M. HIV and the gut microbiota: composition, consequences, and avenues for amelioration. Curr HIV/AIDS Rep. 2019;16:204–13.
Adimora AA, Ramirez C, Benning L, Greenblatt RM, Kempf M-C, Tien PC, Kassaye SG, Anastos K, Cohen M, Minkoff H, et al. Cohort pProfile: tThe Women’s Interagency HIV Study (WIHS). Int J Epidemiol. 2018;47:393–394i.
Kaslow RA, Ostrow DG, Detels R, Phair JP, Polk BF, Rinaldo CR Jr. The Multicenter AIDS Cohort Study: rationale, organization, and selected characteristics of the participants. Am J Epidemiol. 1987;126:310–8.
D’Souza G, Bhondoekhan F, Benning L, Margolick JB, Adedimeji AA, Adimora AA, Alcaide ML, Cohen MH, Detels R, Friedman MR, et al. Characteristics of the MACS/WIHS combined cohort study: opportunities for research on aging with hiv in the longest US observational study of HIV. Am J Epidemiol. 2021;190:1457–75.
Wang Z, Peters BA, Usyk M, Xing J, Hanna DB, Wang T, Post WS, Landay AL, Hodis HN, Weber K, et al: Gut microbiota, plasma metabolomic profiles, and carotid artery atherosclerosis in HIV infection. Arterioscler Thromb Vasc Biol 2022:101161ATVBAHA121317276.
Wang Z, Peters BA, Bryant M, Hanna DB, Schwartz T, Wang T, Sollecito CC, Usyk M, Grassi E, Wiek F, et al. Gut microbiota, circulating inflammatory markers and metabolites, and carotid artery atherosclerosis in HIV infection. Microbiome. 2023;11:119.
Qi Q, Hua S, Clish CB, Scott JM, Hanna DB, Wang T, Haberlen SA, Shah SJ, Glesby MJ, Lazar JM, et al. Plasma tryptophan-kynurenine metabolites are altered in human immunodeficiency virus infection and associated with progression of carotid artery atherosclerosis. Clin Infect Dis. 2018;67:235–42.
Zhang E, Chai JC, Deik AA, Hua S, Sharma A, Schneider MF, Gustafson D, Hanna DB, Lake JE, Rubin LH, et al. Plasma lipidomic profiles and risk of diabetes: 2 prospective cohorts of HIV-infected and HIV-uninfected individuals. J Clin Endocrinol Metab. 2021;106:999–1010.
Kaul A, Mandal S, Davidov O, Peddada SD: Analysis of microbiome data in the presence of excess zeros. Frontiers in Microbiology 2017, 8.
Oksanen J, Kindt R, Legendre P, O’Hara B, Stevens MHH, Oksanen MJ, Suggests M. The vegan package. Community ecology package. 2007;10:719.
Rohart F, Gautier B, Singh A, Le Cao KA. mixOmics: aAn R package for ’omics feature selection and multiple data integration. PLoS Comput Biol. 2017;13:e1005752.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.
Shi B, Choirat C, Coull BA, VanderWeele TJ, Valeri L. CMAverse: a suite of functions for reproducible causal mediation analyses. Epidemiology. 2021;32:e20–2.
Wu H, Esteve E, Tremaroli V, Khan MT, Caesar R, Manneras-Holm L, Stahlman M, Olsson LM, Serino M, Planas-Felix M, et al. Metformin alters the gut microbiome of individuals with treatment-naive type 2 diabetes, contributing to the therapeutic effects of the drug. Nat Med. 2017;23:850–8.
Forslund K, Hildebrand F, Nielsen T, Falony G, Le Chatelier E, Sunagawa S, Prifti E, Vieira-Silva S, Gudmundsdottir V, Pedersen HK, et al. Disentangling type 2 diabetes and metformin treatment signatures in the human gut microbiota. Nature. 2015;528:262–6.
Forslund SK, Chakaroun R, Zimmermann-Kogadeeva M, Marko L, Aron-Wisnewsky J, Nielsen T, Moitinho-Silva L, Schmidt TSB, Falony G, Vieira-Silva S, et al. Combinatorial, additive and dose-dependent drug-microbiome associations. Nature. 2021;600:500–5.
Bennett BJ, de Aguiar Vallim TQ, Wang Z, Shih DM, Meng Y, Gregory J, Allayee H, Lee R, Graham M, Crooke R, et al. Trimethylamine-N-oxide, a metabolite associated with atherosclerosis, exhibits complex genetic and dietary regulation. Cell Metab. 2013;17:49–60.
Nemet I, Saha PP, Gupta N, Zhu W, Romano KA, Skye SM, Cajka T, Mohan ML, Li L, Wu Y, et al. A cardiovascular disease-linked gut microbial metabolite acts via adrenergic receptors. Cell. 2020;180(862–877):e822.
Koh A, Molinaro A, Stahlman M, Khan MT, Schmidt C, Manneras-Holm L, Wu H, Carreras A, Jeong H, Olofsson LE, et al. Microbially produced imidazole propionate impairs insulin signaling through mTORC1. Cell. 2018;175(947–961):e917.
Letchumanan G, Abdullah N, Marlini M, Baharom N, Lawley B, Omar MR, Mohideen FBS, Addnan FH, Nur Fariha MM, Ismail Z, Pathmanathan SG. Gut microbiota composition in prediabetes and newly diagnosed type 2 diabetes: a systematic review of observational studies. Front Cell Infect Microbiol. 2022;12:943427.
Zhong H, Ren H, Lu Y, Fang C, Hou G, Yang Z, Chen B, Yang F, Zhao Y, Shi Z, et al. Distinct gut metagenomics and metaproteomics signatures in prediabetics and treatment-naive type 2 diabetics. EBioMedicine. 2019;47:373–83.
Doumatey AP, Adeyemo A, Zhou J, Lei L, Adebamowo SN, Adebamowo C, Rotimi CN. Gut microbiome profiles are associated with type 2 diabetes in urban Africans. Front Cell Infect Microbiol. 2020;10:63.
Nagata N, Nishijima S, Miyoshi-Akiyama T, Kojima Y, Kimura M, Aoki R, Ohsugi M, Ueki K, Miki K, Iwata E, et al. Population-level metagenomics uncovers distinct effects of multiple medications on the human gut microbiome. Gastroenterology. 2022;163:1038–52.
Whang A, Nagpal R, Yadav H. Bi-directional drug-microbiome interactions of anti-diabetics. EBioMedicine. 2019;39:591–602.
Anhe FF, Jensen BAH, Varin TV, Servant F, Van Blerk S, Richard D, Marceau S, Surette M, Biertho L, Lelouvier B, et al. Type 2 diabetes influences bacterial tissue compartmentalisation in human obesity. Nat Metab. 2020;2:233–42.
Song YL, Liu CX, McTeague M, Summanen P, Finegold SM. Clostridium bartlettii sp. nov., isolated from human faeces. Anaerobe. 2004;10:179–84.
Tsukuda N, Yahagi K, Hara T, Watanabe Y, Matsumoto H, Mori H, Higashi K, Tsuji H, Matsumoto S, Kurokawa K, Matsuki T. Key bacterial taxa and metabolic pathways affecting gut short-chain fatty acid profiles in early life. ISME J. 2021;15:2574–90.
Samarani S, Allam O, Sagala P, Aldabah Z, Jenabian MA, Mehraj V, Tremblay C, Routy JP, Amre D, Ahmad A. Imbalanced production of IL-18 and its antagonist in human diseases, and its implications for HIV-1 infection. Cytokine. 2016;82:38–51.
Thorand B, Kolb H, Baumert J, Koenig W, Chambless L, Meisinger C, Illig T, Martin S, Herder C. Elevated levels of interleukin-18 predict the development of type 2 diabetes: results from the MONICA/KORA Augsburg Study, 1984–2002. Diabetes. 2005;54:2932–8.
Tapan S, Dogru T, Kara M, Ercin CN, Kilciler G, Genc H, Sertoglu E, Acikel C, Kilic S, Karslioglu Y. Circulating levels of interleukin-18 in patients with non-alcoholic fatty liver disease. Scand J Clin Lab Invest. 2010;70:399–403.
Nowarski R, Jackson R, Gagliani N, de Zoete MR, Palm NW, Bailis W, Low JS, Harman CC, Graham M, Elinav E, Flavell RA. Epithelial IL-18 equilibrium controls barrier function in colitis. Cell. 2015;163:1444–56.
La Reau AJ, Suen G. The Ruminococci: key symbionts of the gut ecosystem. J Microbiol. 2018;56:199–208.
Maruo T, Sakamoto M, Ito C, Toda T, Benno Y. Adlercreutzia equolifaciens gen. nov., sp. nov., an equol-producing bacterium isolated from human faeces, and emended description of the genus Eggerthella. Int J Syst Evol Microbiol. 2008;58:1221–7.
Boyce BF, Xing L. Biology of RANK, RANKL, and osteoprotegerin. Arthritis Res Ther. 2007;9(Suppl 1):S1.
Chen Y, Wang X, Zhang C, Liu Z, Li C, Ren Z. Gut microbiota and bone diseases: a growing partnership. Front Microbiol. 2022;13:877776.
Li L, Rao S, Cheng Y, Zhuo X, Deng C, Xu N, Zhang H, Yang L. Microbial osteoporosis: The interplay between the gut microbiota and bones via host metabolism and immunity. Microbiologyopen. 2019;8:e00810.
He Y, Chen Y. The potential mechanism of the microbiota-gut-bone axis in osteoporosis: a review. Osteoporos Int. 2022;33:2495–506.
Hofbauer LC, Brueck CC, Singh SK, Dobnig H. Osteoporosis in patients with diabetes mellitus. J Bone Miner Res. 2007;22:1317–28.
Xiang GD, Xu L, Zhao LS, Yue L, Hou J. The relationship between plasma osteoprotegerin and endothelium-dependent arterial dilation in type 2 diabetes. Diabetes. 2006;55:2126–31.
Secchiero P, Corallini F, Pandolfi A, Consoli A, Candido R, Fabris B, Celeghini C, Capitani S, Zauli G. An increased osteoprotegerin serum release characterizes the early onset of diabetes mellitus and may contribute to endothelial cell dysfunction. Am J Pathol. 2006;169:2236–44.
Zeybel M, Arif M, Li X, Altay O, Yang H, Shi M, Akyildiz M, Saglam B, Gonenli MG, Yigit B, et al. Multiomics analysis reveals the impact of microbiota on host metabolism in hepatic steatosis. Adv Sci (Weinh). 2022;9:e2104373.
Zhang C, Luo X, Chen J, Zhou B, Yang M, Liu R, Liu D, Gu HF, Zhu Z, Zheng H. Osteoprotegerin promotes liver steatosis by targeting the ERK–PPAR-γ–CD36 pathway. Diabetes. 2019;68:1902–14.
Oliveira AG, Araujo TG, Carvalho BM, Rocha GZ, Santos A, Saad MJA. The role of hepatocyte growth factor (HGF) in insulin resistance and diabetes. Front Endocrinol (Lausanne). 2018;9:503.
Piquer-Garcia I, Campderros L, Taxeras SD, Gavalda-Navarro A, Pardo R, Vila M, Pellitero S, Martinez E, Tarasco J, Moreno P, et al. A role for oncostatin M in the impairment of glucose homeostasis in obesity. J Clin Endocrinol Metab. 2020;105:e337–348.
Fisher FM, Maratos-Flier E. Understanding the Physiology of FGF21. Annu Rev Physiol. 2016;78:223–41.
Post A, Dam WA, Sokooti S, Groothof D, Gloerich J, van Gool AJ, Kremer D, Gansevoort RT, van den Born J, Kema IP, et al. Circulating FGF21 concentration, fasting plasma glucose, and the risk of type 2 diabetes: results from the PREVEND study. J Clin Endocrinol Metab. 2023;108:1387–93.
Larsson SC, Michaelsson K, Mola-Caminal M, Hoijer J, Mantzoros CS. Genome-wide association and Mendelian randomization study of fibroblast growth factor 21 reveals causal associations with hyperlipidemia and possibly NASH. Metabolism. 2022;137:155329.
Mashili FL, Austin RL, Deshmukh AS, Fritz T, Caidahl K, Bergdahl K, Zierath JR, Chibalin AV, Moller DE, Kharitonenkov A, Krook A. Direct effects of FGF21 on glucose uptake in human skeletal muscle: implications for type 2 diabetes and obesity. Diabetes Metab Res Rev. 2011;27:286–97.
Pan Q, Lin S, Li Y, Liu L, Li X, Gao X, Yan J, Gu B, Chen X, Li W, et al. A novel GLP-1 and FGF21 dual agonist has therapeutic potential for diabetes and non-alcoholic steatohepatitis. EBioMedicine. 2021;63:103202.
Loomba R, Sanyal AJ, Kowdley KV, Bhatt DL, Alkhouri N, Frias JP, Bedossa P, Harrison SA, Lazas D, Barish R, et al. Randomized, controlled trial of the FGF21 analogue pegozafermin in NASH. N Engl J Med. 2023;389:998–1008.
Domingo P, Gallego-Escuredo JM, Domingo JC, Gutierrez Mdel M, Mateo MG, Fernandez I, Vidal F, Giralt M, Villarroya F. Serum FGF21 levels are elevated in association with lipodystrophy, insulin resistance and biomarkers of liver injury in HIV-1-infected patients. AIDS. 2010;24:2629–37.
Zafar MI, Mills K, Ye X, Blakely B, Min J, Kong W, Zhang N, Gou L, Regmi A, Hu SQ, et al. Association between the expression of vascular endothelial growth factors and metabolic syndrome or its components: a systematic review and meta-analysis. Diabetol Metab Syndr. 2018;10:62.
Gudmundsdottir V, Zaghlool SB, Emilsson V, Aspelund T, Ilkov M, Gudmundsson EF, Jonsson SM, Zilhao NR, Lamb JR, Suhre K, et al. Circulating protein signatures and causal candidates for type 2 diabetes. Diabetes. 2020;69:1843–53.
Ullah A, Ud Din A, Ding W, Shi Z, Pervaz S, Shen B: A narrative review: CXC chemokines influence immune surveillance in obesity and obesity-related diseases: type 2 diabetes and nonalcoholic fatty liver disease. Rev Endocr Metab Disord 2023:1–21.
Sun MY, Wang SJ, Li XQ, Shen YL, Lu JR, Tian XH, Rahman K, Zhang LJ, Nian H, Zhang H. CXCL6 promotes renal interstitial fibrosis in diabetic nephropathy by activating JAK/STAT3 signaling pathway. Front Pharmacol. 2019;10:224.
Molinaro A, Bel Lassen P, Henricsson M, Wu H, Adriouch S, Belda E, Chakaroun R, Nielsen T, Bergh PO, Rouault C, et al. Imidazole propionate is increased in diabetes and associated with dietary patterns and altered microbial ecology. Nat Commun. 2020;11:5881.
Dekkers KF, Sayols-Baixeras S, Baldanzi G, Nowak C, Hammar U, Nguyen D, Varotsis G, Brunkwall L, Nielsen N, Eklund AC, et al. An online atlas of human plasma metabolite signatures of gut microbiome composition. Nat Commun. 2022;13:5370.
Zhu Y, Dwidar M, Nemet I, Buffa JA, Sangwan N, Li XS, Anderson JT, Romano KA, Fu X, Funabashi M, et al. Two distinct gut microbial pathways contribute to meta-organismal production of phenylacetylglutamine with links to cardiovascular disease. Cell Host Microbe. 2023;31(18–32):e19.
Xu J, Cai M, Wang Z, Chen Q, Han X, Tian J, Jin S, Yan Z, Li Y, Lu B, Lu H: Phenylacetylglutamine as a novel biomarker of type 2 diabetes with distal symmetric polyneuropathy by metabolomics. J Endocrinol Invest 2022.
Brown EM, Clardy J, Xavier RJ. Gut microbiome lipid metabolism and its impact on host physiology. Cell Host Microbe. 2023;31:173–86.
Collins SL, Patterson AD. The gut microbiome: an orchestrator of xenobiotic metabolism. Acta Pharm Sin B. 2020;10:19–32.
Estrada V, Gonzalez N. Gut microbiota in diabetes and HIV: Inflammation is the link. EBioMedicine. 2018;38:17–8.
Morze J, Wittenbecher C, Schwingshackl L, Danielewicz A, Rynkiewicz A, Hu FB, Guasch-Ferre M. Metabolomics and type 2 diabetes risk: an updated systematic review and meta-analysis of prospective cohort studies. Diabetes Care. 2022;45:1013–24.
Zhuang H, Han J, Cheng L, Liu SL. A positive causal influence of IL-18 levels on the risk of T2DM: a Mendelian randomization study. Front Genet. 2019;10:295.
Akbari M, Hassan-Zadeh V. IL-6 signalling pathways and the development of type 2 diabetes. Inflammopharmacology. 2018;26:685–98.
Wegmann U, Louis P, Goesmann A, Henrissat B, Duncan SH, Flint HJ. Complete genome of a new F irmicutes species belonging to the dominant human colonic microbiota (‘R uminococcus bicirculans’) reveals two chromosomes and a selective capacity to utilize plant glucans. Environ Microbiol. 2014;16:2879–90.
Oñate FP, Chamignon C, Burz SD, Lapaque N, Monnoye M, Philippe C, Bredel M, Chêne L, Farin W, Paillarse J-M. Adlercreutzia equolifaciens is an anti-inflammatory commensal bacterium with decreased abundance in gut microbiota of patients with metabolic liver disease. Int J Mol Sci. 2023;24:12232.
Brown TT, Tassiopoulos K, Bosch RJ, Shikuma C, McComsey GA. Association between systemic inflammation and incident diabetes in HIV-infected patients after initiation of antiretroviral therapy. Diabetes Care. 2010;33:2244–9.
Kettelhut A, Bowman E, Funderburg NT: Immunomodulatory and anti-inflammatory strategies to reduce comorbidity risk in people with HIV. Curr HIV/AIDS Rep 2020, 17:394-404
Kaur D, Chachi L, Gomez E, Sylvius N, Brightling CE. Interleukin-18, IL-18 binding protein and IL-18 receptor expression in asthma: a hypothesis showing IL-18 promotes epithelial cell differentiation. Clin Transl Immunology. 2021;10:e1301.
Ihim SA, Abubakar SD, Zian Z, Sasaki T, Saffarioun M, Maleknia S, Azizi G. Interleukin-18 cytokine in immunity, inflammation, and autoimmunity: biological role in induction, regulation, and treatment. Front Immunol. 2022;13:919973.
Thomas JM, Huuskes BM, Sobey CG, Drummond GR, Vinh A. The IL-18/IL-18R1 signalling axis: diagnostic and therapeutic potential in hypertension and chronic kidney disease. Pharmacol Ther. 2022;239:108191.
Huang HK, Chuang AT, Liao TC, Shao SC, Liu PP, Tu YK, Lai EC. Denosumab and the risk of diabetes in patients treated for osteoporosis. JAMA Netw Open. 2024;7:e2354734.
Xing B, Yu J, Zhang H, Li Y. RANKL inhibition: a new target of treating diabetes mellitus? Ther Adv Endocrinol Metab. 2023;14:20420188231170750.
Agus A, Planchais J, Sokol H. Gut microbiota regulation of tryptophan metabolism in health and disease. Cell Host Microbe. 2018;23:716–24.
Acknowledgements
The authors gratefully acknowledge the contributions of the study participants and dedication of the staff at the MWCCS sites.
Disclaimer
The contents of this article are solely the responsibility of the authors and do not represent the official views of the National Institutes of Health (NIH).
Financial support
Q.Q. was supported by the National Heart, Lung, and Blood Institute (NHLBI) (R01-HL140976 and R01-HL170904), the National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK) (R01-DK119268 and R01-DK126698), and Feldstein Medical Foundation Research Grant. K.L. was supported by an AHA postdoctoral fellowship (23POST1020455). Other related funding sources included K01-HL129892, K01-HL137557, R01-HL126543, R01-HL132794, R01-HL083760, and R01-HL095140 from the NHLBI, R01-MD011389 from the National Institute on Minority Health and Health Disparities, U01-AI035004 and the Einstein-Rockefeller-CUNY Center for AIDS Research (P30-AI124414) funded by the National Institute of Allergy and Infectious Diseases, and Einstein Cancer Research Center (P30-CA013330) funded by National Cancer Institute.
MWCCS (Principal Investigators)
Atlanta CRS (Ighovwerha Ofotokun, Anandi Sheth, and Gina Wingood), U01-HL146241; Baltimore CRS (Todd Brown and Joseph Margolick), U01-HL146201; Bronx CRS (Kathryn Anastos, David Hanna, and Anjali Sharma), U01-HL146204; Brooklyn CRS (Deborah Gustafson and Tracey Wilson), U01-HL146202; Data Analysis and Coordination Center (Gypsyamber D’Souza, Stephen Gange and Elizabeth Topper), U01-HL146193; Chicago-Cook County CRS (Mardge Cohen and Audrey French), U01-HL146245; Chicago-Northwestern CRS (Steven Wolinsky), U01-HL146240; Northern California CRS (Bradley Aouizerat, Jennifer Price, and Phyllis Tien), U01-HL146242; Los Angeles CRS (Roger Detels and Matthew Mimiaga), U01-HL146333; Metropolitan Washington CRS (Seble Kassaye and Daniel Merenstein), U01-HL146205; Miami CRS (Maria Alcaide, Margaret Fischl, and Deborah Jones), U01-HL146203; Pittsburgh CRS (Jeremy Martinson and Charles Rinaldo), U01-HL146208; UAB-MS CRS (Mirjam-Colette Kempf, Jodie Dionne-Odom, and Deborah Konkle-Parker), U01-HL146192; UNC CRS (Adaora Adimora and Michelle Floris-Moore), U01-HL146194. The MWCCS is funded primarily by the National Heart, Lung, and Blood Institute (NHLBI), with additional co-funding from the Eunice Kennedy Shriver National Institute Of Child Health & Human Development (NICHD), National Institute On Aging (NIA), National Institute Of Dental & Craniofacial Research (NIDCR), National Institute Of Allergy And Infectious Diseases (NIAID), National Institute Of Neurological Disorders And Stroke (NINDS), National Institute Of Mental Health (NIMH), National Institute On Drug Abuse (NIDA), National Institute Of Nursing Research (NINR), National Cancer Institute (NCI), National Institute on Alcohol Abuse and Alcoholism (NIAAA), National Institute on Deafness and Other Communication Disorders (NIDCD), National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK), National Institute on Minority Health and Health Disparities (NIMHD), and in coordination and alignment with the research priorities of the National Institutes of Health, Office of AIDS Research (OAR). MWCCS data collection is also supported by UL1-TR000004 (UCSF CTSA), UL1-TR003098 (JHU ICTR), UL1-TR001881 (UCLA CTSI), P30-AI-050409 (Atlanta CFAR), P30-AI-073961 (Miami CFAR), P30-AI-050410 (UNC CFAR), P30-AI-027767 (UAB CFAR), and P30-MH-116867 (Miami CHARM).
Funding
Data in this manuscript were collected by the Women’s Interagency HIV Study (WIHS) and the Multicenter AIDS Cohort Study (MACS), now the MACS/WIHS Combined Cohort Study (MWCCS).
Author information
Authors and Affiliations
Contributions
K.L. and Q.Q. conceived the idea and designed the study. K.L. conducted all statistical analyses and drafted and revised the manuscript. Z.W., M.U., T.W., and C.B.C performed bioinformatic pipelines to generate metabolomics, proteomics, and metagenomics data and revised the manuscript. B.A.P., J.M., X.X., D.B.H., A.L.L., M.F.S, D.G., K.M.W., A.F., A.S., K.A., T.B., and R.D.B. critically revised the manuscript and provided valuable input in results interpretation. R.C.K acquired funding and revised the manuscript. Q.Q. acquired funding, was responsible for the overall study design, data, results interpretation, and critically revised the manuscript. All coauthors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
All participants provided written informed content and the present study was reviewed and approved by the Institutional Review Board (IRB) at Albert Einstein College of Medicine (IRB no: 2018–9386). The present study adhered to the principles of Helsinki Declaration.
Consent for publication
All authors approved the final manuscript before submission.
Competing interests
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential competing interest.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1:
Supplementary Methods. Additional descriptions on the drug deconfounding analysis and mediation analyses of multiple mediators.
Additional file 2: Fig S1.
An overview of analysis and sample sizes; Fig S2. Partial spearman correlation (PCorr) among diabetes-associated gut bacteria genera and species within each genus; Fig S3. Results of ANCOM-II for prevalent diabetes at the species level; Fig S4. Confounding analysis of antidiabetic medications on selected diabetes-associated gut bacterial genera; Fig S5. Associations of metabolites with incident diabetes; Fig S6. Partial correlation (PCorr) among 76 diabetes-associated metabolites; Fig S7. Partial correlation among diabetes-associated individual lipids and four lipid modules; Fig S8. Partial correlation (PCorr) among diabetes-associated gut bacteria genera, nonlipid metabolites, lipid modules, and proteins; Fig S9. Relationship between associations of omics features with HIV positive status and associations with prevalent diabetes.
Additional file 3: Table S1.
Population characteristics in 493 WIHS women.
Additional file 4: Table S2.
The identification of Genera associated with prevalent diabetes in ANCOM-II.
Additional file 5: Table S3.
Associations of selected gut bacteria genera with prevalent diabetes by adjustment models.
Additional file 6: Table S4.
Associations of metabolites (and lipids modules) and proteins with diabetes by adjustment models.
Additional file 7: Table S4.
Associations of metabolites (and lipids modules) and proteins with diabetes by adjustment models.
Additional file 8: Table S6.
Associations of diabetes associated genera with diabetes associated metabolites (and Lipid modules) and proinflammatory proteins.
Additional file 9: Table S7.
Associations of selected microbial genera, proteins, lipid modules and nonlipids metabolites with diabetes by HIV serostatus.
Additional file 10: Table S8.
Levels of diabetes associated microbial genera, proteins, lipid modules and individual nonlipids metabolites by HIV characteristics.
Additional file 11: Table S9.
HIV infection-associated gut bacteria, proteins and metabolites and prevalent diabetes in WIHS women.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Luo, K., Peters, B.A., Moon, JY. et al. Metabolic and inflammatory perturbation of diabetes associated gut dysbiosis in people living with and without HIV infection. Genome Med 16, 59 (2024). https://doi.org/10.1186/s13073-024-01336-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13073-024-01336-1