- 1Senior Department of Tuberculosis, Chinese PLA General Hospital, Beijing, China
- 2Graduate School, Hebei North University, Zhangjiakou, Hebei, China
- 3Department of Geriatrics, the Eighth Medical Center of PLA General Hospital, Beijing, China
- 4Handan Municipal Centre for Disease Prevention and Control, Handan, Hebei, China
Background: Tuberculosis-diabetes mellitus (TB-DM) comorbidity presents significant clinical challenges due to poor treatment outcomes. This study investigated peripheral blood lymphocyte profiles and cytokine dynamics in TB-DM patients compared to healthy controls (HCs) and DM patients.
Methods: Subjects from the healthy controls (HCs), DM, and TB-DM were recruited, and peripheral blood samples were collected. The absolute counts of lymphocyte subsets were detected by flow cytometry, and the cytokines were quantitatively analyzed using the Olink ultra-sensitive targeted protein detection technology for micro-samples. Methods such as differential expression analysis, principal component analysis (PCA), correlation analysis, KEGG pathway enrichment analysis, and GO functional annotation were used to screen out the biomarkers related to TB-DM. Based on this, a TB-DM internal model performance was constructed, and the receiver operating characteristic (ROC) curve was used to evaluate its diagnostic efficacy.
Results: The study demonstrated significantly reduced NK cells (PTB-DM vs. HC < 0.0001 and PTB-DM vs. DM = 0.0292), total T cells (PTB-DM vs. HC = 0.0018 and PTB-DM vs. DM < 0.0001), and CD8+ T cells (PTB-DM vs. HC = 0.0009 and PTB-DM vs. DM = 0.0072) in TB-DM versus HCs and DM groups. TB-DM patients showed decreased CD4+ T (PTB-DM vs. DM < 0.0001) and B cells (PTB-DM vs. DM = 0.0004) compared to DM controls. Cytokine profiling revealed 5 upregulated and 17 downregulated factors in TB-DM. Three biomarkers (IL-6, IFN-γ, CXCL10) demonstrated superior diagnostic performance (AUC = 0.9841, sensitivity=88.89%, specificity=92.86%) when combined.
Conclusion: Our findings identify distinct immunological alterations in TB-DM and propose a novel cytokine-based diagnostic panel for this high-risk population.
1 Introduction
Tuberculosis (TB), a chronic respiratory infectious disease caused by Mycobacterium tuberculosis (MTB), remains a critical global health challenge (1, 2). According to the World Health Organization (WHO) 2025 Global Tuberculosis Report, an estimated 10.7 million individuals developed TB worldwide in 2024, with 1.23 million deaths, maintaining its status as the leading cause of mortality from a single infectious pathogen (3). Concurrently, diabetes mellitus (DM), the fastest-growing endocrine metabolic disorder, presents a substantial epidemiological burden. The International Diabetes Federation (IDF) reported 589 million adults aged 20–79 living with DM in 2021, projected to rise to 853 million by 2050 (4). Notably, the bidirectional interaction between these two diseases has emerged as a major public health concern (5–8).
Accumulating evidence demonstrates a synergistic detrimental relationship between DM and TB (9–11). DM increases the risk of active tuberculosis (ATB) development by 2-4-fold, establishing it as a critical metabolic risk factor for TB (8, 12). A meta-analysis encompassing 22,658 studies revealed a global prevalence of 13.73% for TB-DM comorbidity, rising to 14.62% in high-burden regions such as Southeast Asia (13). DM compromises host immunity, exacerbating susceptibility to ATB, while TB-induced stress hyperglycemia further disrupts glycemic control (14, 15). Radiological analyses indicate that tuberculosis-diabetes comorbidity (TB-DM) patients exhibit higher rates of pulmonary cavitation on chest X-rays and increased incidence of adverse drug reactions during anti-TB therapy (16, 17). Importantly, DM significantly elevates treatment failure rates, mortality, and relapse rates in TB patients compared to non-diabetic counterparts (18). These findings underscore the profound impact of DM on TB pathogenesis, therapeutic outcomes, and prognosis (8).
Despite advances in understanding TB-DM interactions, critical knowledge gaps persist. Current diagnostic frameworks lack validated biomarkers for early comorbidity detection, and systematic characterization of immunometabolic alterations in TB-DM patients remains incomplete. This study employs flow cytometry and proteomic profiling to elucidate cellular and molecular immune signatures in TB-DM cohorts. Our objectives are to construct a TB-DM-specific immune profile and identify diagnostic biomarkers with clinical utility. These findings will establish novel immunological monitoring parameters and therapeutic strategies, and thereby aims to facilitate the early and accurate identification of tuberculosis co-infection in individuals with diabetes, ultimately improving clinical outcomes for TB-DM patients.
2 Materials and methods
2.1 Study design and ethical statement
This cross-sectional study was conducted at the Eighth Medical Center of PLA General Hospital from March to November 2022. Two independent population cohorts were enrolled, each comprising three subgroups: healthy controls (HCs), DM, and TB-DM patients. Inclusion criteria for HC included no history of tuberculosis exposure or infection, normal blood glucose levels, unremarkable chest X-ray findings, HIV-negative status, and age ≥18 years. Exclusion criteria consisted of: identifiable additional risk of TB exposure (such as occupation in TB care, close contact with ATB patients, or recent residence in or travel to high−incidence areas within China), prior history of TB or residual pulmonary lesions, contraindications to interferon−γ release assay (IGRA), and diagnosis of DM or HIV infection.
DM participants met the American Diabetes Association (ADA) diagnostic criteria (fasting blood glucose ≥7.0 mmol/L, or 2−h postprandial glucose ≥11.1 mmol/L, or HbA1c ≥6.5%, or random plasma glucose ≥11.1 mmol/L in the presence of classic hyperglycemic symptoms or hyperglycemic crisis) (19). Individuals with HIV infection, agranulocytosis, autoimmune diseases, severe hepatic or renal dysfunction, malignancy, immunomodulatory therapy, or pregnancy were excluded.
TB−DM participants fulfilled both the diagnostic criteria for DM and for ATB. All patients enrolled in the TB−DM group were required to have a documented history of DM, recorded by a formal medical institution prior to the current ATB diagnosis. This medical history was verified through the electronic medical record system. ATB was confirmed per Diagnostic Criteria for Pulmonary Tuberculosis (WS288-2017), requiring microbiological evidence (e.g., sputum smear positivity, mycobacterial culture, or nucleic acid detection) alongside clinical, radiological, and epidemiological findings. Enrolled TB-DM patients were aged 18–60 years, diagnosed with ATB ≤2 weeks without prior anti-tuberculosis therapy. Exclusion criteria included corticosteroid use, immunocompromised status (e.g., HIV, organ transplantation, autoimmune disorders), malnutrition, or ATB diagnosis >2 weeks with prior treatment.
Ethical approval was obtained from the Ethics Committee of the Eighth Medical Center of PLA General Hospital (Approval No. 309202204080808), with written informed consent from all participants.
2.2 Lymphocyte subset analysis
In Cohort I, 5 mL peripheral blood was collected per participant for lymphocyte profiling. Absolute counts of B cells, NK cells, NKT cells, total T cells, CD4+; T cells, and CD8+; T cells were quantified using 50 μL whole blood via flow cytometry. Reverse pipetting technique was employed to transfer 50 μL anticoagulated whole blood into Trucount™ tubes preloaded with 10 μL Multitest™ 6-color TBNK reagent. After vortex mixing and 15-minute dark incubation at room temperature, 450 μL lysing solution was added, followed by 15-minute incubation. Samples were analyzed using a FACSAria II flow cytometer.
2.3 Olink ultra-sensitive targeted protein detection for micro-samples
For Cohort II, 10 mL of the patient’s whole blood sample was centrifuged at 3000 rpm for 3–5 minutes, extracted the upper-layer plasma, and stored it at -80°C for future use. According to the manufacturer’s instructions, the Olink Target 96 Inflammation Panel is used to quantify the serum protein levels. Olink’s unique Proximity Extension Assay (PEA) provides a unique and efficient tool for the discovery and development of protein biomarkers. Antibody pairs labeled with DNA oligonucleotides specifically bind to the target antigens in solution. Through DNA polymerase, the adjacent DNA oligonucleotide chains hybridize and extend. The newly generated DNA barcodes are amplified by standard polymerase chain reaction (PCR). The amplified samples are transferred to an integrated microfluidic chip (IFC), loaded into the instrument for quantitative polymerase chain reaction (qPCR), and the raw data are generated (20). The bcl files generated by the qPCR detection are converted into counts files, and the NPX Manager software is used for quality control and data normalization. The final measurement results are expressed as the normalized protein expression (NPX).
2.4 Bioinformatics analysis
The R package “Olink Analyze” was employed to identify differentially expressed proteins (DEPs) between groups, using a statistical threshold of |log2 fold change| ≥ 1 and a P value < 0.05. We used the “ggplot2” package to construct a volcano plot to visualize the analysis results of these DEPs. The Kyoto Encyclopedia of Genes and Genomes (KEGG) online database and the Gene Ontology (GO) database were employed to conduct KEGG and GO enrichment analyses on the DEPs. The KEGG database helps to understand the advanced functions and utilities of biological systems (such as cells, organisms, and ecosystems) at the molecular level, especially from the large-scale molecular data generated by genomic sequencing and other high-throughput experimental techniques (21). The GO database is one of the most commonly used gene annotation databases. With the annotations in the GO database, functional enrichment analysis of gene sets can be carried out (22). The GO database is mainly divided into three parts, which describe the molecular function (MF), cellular component (CC), and biological process (BP) of genes respectively (23).
2.5 Statistical analysis
Statistical analyses were performed using GraphPad Prism (version 10.4.1) and SPSS (version 27.0). For immune cell counts and cytokine expression levels within a one−factor multi−level design, data distribution was first assessed. The normality of each group was tested individually using the Shapiro−Wilk test, and homogeneity of variances was verified using the Brown−Forsythe test. If all groups satisfied both normality and homoscedasticity, one−way analysis of variance (ANOVA) was applied, followed by Tukey’s post−hoc test for multiple comparisons when the overall ANOVA was significant. Otherwise, the non−parametric Kruskal−Wallis test was used, with Dunn’s test for post−hoc comparisons following a significant result.
Principal component analysis (PCA) was performed to evaluate variance contributions. Associations between variables were examined using Pearson correlation analysis, and significant correlations were further modeled by linear regression (reported with R² values and regression equations). The diagnostic efficacy of individual TB−DM biomarkers was evaluated by receiver operating characteristic (ROC) curve analysis, with confidence intervals calculated using the Wilson/Brown method. Additionally, composite biomarker probabilities derived from logistic regression models were validated using ROC analysis. A p−value < 0.05 was considered statistically significant.
3 Results
3.1 General characteristics of the included cases
A total of 104 cases were enrolled in this study, divided into two cohorts: cohort I (67 cases) and cohort II (37 cases) (Figure 1). Cohort I included 23 HCs, 26 DM, and 18 TB-DM patients; gender distribution was male, 52 cases (77.6%), female, 15 cases (22.4%); age distributions were 35.74 ± 12.39 years for HCs, 46.81 ± 11.09 years for DM, and 52.44 ± 14.53 years for TB-DM. Cohort II included 14 HCs, 14 DM, and 9 TB-DM cases; gender distribution was male, 30 cases (81.1%), female, 7 cases (18.9%); age distributions were 39.07 ± 13.21 years for HCs, 43.57 ± 11.27 years for DM, and 51.78 + 15.79 years for TB-DM. Detailed baseline characteristics for each cohort are presented in Supplementary Tables S1 and S2.
Figure 1. Participant Recruitment and Inclusion Process. A total of 309 participants were initially enrolled in this study. Based on the predefined inclusion and exclusion criteria, the following individuals were excluded: those younger than 18 or older than 60 years of age (n = 68), individuals with immunocompromising conditions (n = 56), individuals on hormone therapy (n = 46), and individuals receiving immunomodulatory treatment (n = 35). Consequently, 104 participants were ultimately included in the final analysis.
3.2 Quantitative analysis of immune cell subpopulations unveils severe imbalances in TB-DM patients
NK cells showed a gradient decline in absolute counts, with HCs significantly higher than DM (P = 0.0215) and DM significantly higher than TB-DM (P = 0.0292). HCs showed a significant difference from TB-DM (P<0.0001) compared to DM (P = 0.0292) (Figure 2A). NKT cells showed no significant differences in absolute counts across the three groups (Figure 2B). Further analysis of adaptive immune cells revealed that TB-DM group had significantly lower total T cells (P = 0.0018) and CD8+ T cells (P = 0.0009) compared to HCs (P = 0.0018 for total T cells and P = 0.0009 for CD8+ T cells). Additionally, TB-DM group showed significantly lower total T cells (P<0.0001) and CD8+ T cells (P = 0.0072) compared to DM (P<0.0001 for total T cells and P = 0.0072 for CD8+ T cells). DM group also showed significantly higher CD4+ T cells compared to HCs (P = 0.0305) (Figures 2C–E). Furthermore, TB-DM group showed significantly lower B cells compared to DM (P = 0.0004) (Figure 2F).
Figure 2. Comparative analysis of lymphocyte subsets across study groups. Absolute counts of lymphocyte populations in HCs, DM patients, and TB-DM comorbid patients. (A) NK cells. (B) Natural killer T (NKT) cells. (C) Total T lymphocytes. (D) CD4+; T lymphocytes. (E) CD8+; T lymphocytes. (F) B lymphocytes. Statistical significance was determined by one-way ANOVA or Kruskal-Wallis tests, depending on normality and variance homogeneity. Statistical significance was set at p<0.05.
3.3 Olink-based profiling reveals dysregulated cytokines in tuberculosis comorbid with diabetes mellitus
Using the Olink ultra-sensitive targeted protein detection technology for micro-samples, 92 cytokines were analyzed. Volcano plots revealed that compared to HCs, the TB-DM group exhibited 6 significantly upregulated cytokines, including interferon-gamma (IFN-γ), interleukin-6 (IL-6), C-X-C motif chemokine ligand 9 (CXCL9), CUB domain-containing protein 1 (CDCP1), CXCL10, and extracellular newly identified receptor for advanced glycation end-products (EN-RAGE), along with 11 significantly downregulated cytokines, namely sulfotransferase family 1A member 1 (ST1A1), TNF-related activation-induced cytokine (TRANCE), CXCL5, CXCL6, monocyte chemoattractant protein-4 (MCP-4), caspase-8 (CASP-8), Axis inhibition protein 1 (AXIN1), cluster of differentiation 6 (CD6), STAM-binding protein (STAMBP), sirtuin 2 (SIRT2), and Fibroblast growth factor 19 (FGF-19) (Figure 3A). Compared to DM group, the TB-DM group exhibited 5 significantly upregulated cytokines, which were IFN-γ, IL-6, CXCL9, CDCP1, and CXCL10, as well as 8 significantly downregulated cytokines, including FGF-21, MCP-4, CXCL5, CXCL6, TRANCE, CASP-8, AXIN1, and CD6 (Figure 3B).
Figure 3. Volcano plots of differentially expressed cytokines. Differential cytokine expression analysis comparing: (A) TB-DM versus HCs (fold change >2.0, adjusted P <0.05); (B) TB-DM versus DM (fold change >2, adjusted P <0.05). Red dots: upregulated proteins; blue dots: downregulated proteins; gray dots: non-significant changes. Dashed lines indicate significance thresholds.
3.4 Cytokine Expression Imbalances Among the Three Groups.
Further analysis of cytokine expression differences among the three groups revealed that 22 cytokines exhibited expression differences in two groups. Compared to HCs and DM, the TB-DM group showed significantly higher IFN-γ, IL-6, CDCP1, CXCL9, and CXCL10 expression (Figure 4A). The TB-DM group also showed significant decreases in Delta/Notch-like EGF repeat-containing (DNER), MCP-4, CASP-8, TRANCE, and CD244 molecule (CD244) expression (Figure 4B). SCF expression showed a gradual decline with advancing disease progression and exhibited significant differences (Figure 4B). Additionally, FGF-21 expression was significantly higher in DM patients compared to HCs and TB-DM patients (Figure 4B).
Figure 4. Cytokine profile alterations in TB-DM patients. (A) Five cytokines significantly elevated in TB-DM compared to both HCs and DM groups. (B) Seventeen cytokines showing marked reduction in TB-DM cohort. NPX (Normalized Protein Expression) values represent log2-transformed protein concentrations measured by Olink. Statistical significance was determined by one-way ANOVA or Kruskal-Wallis tests, depending on normality and variance homogeneity. Statistical significance was set at p<0.05.
3.5 Insights from PCA, correlation matrix analysis and linear regression analysis in different groups
To systematically depict the peripheral immune-regulatory landscape across HC, DM and TB-DM, we performed PCA (Figures 5A–C), Pearson correlation matrices (Tables 1, 2; Figures 5D–F) and univariate linear regression of the top-correlated cytokine pairs (Figures 5G–K). The data reveal a stepwise remodelling of the cytokine covariance architecture that follows an “equilibrium–loosening–polarization” trajectory (Figures 5A–C). In HC, PC1 and PC2 jointly account for 54.09% of the variance, with SCF, IFN-γ, CXCL9, CXCL10 and CDCP1 clustering separately from the remaining IL-6-containing set (Figure 5A); this bipartite structure mirrors the balanced interplay shown in Figure 5D, where a positive backbone [AXIN1–TNFSF14 (Tumor Necrosis Factor Superfamily Member 14), CXCL6–STAMBP (STAM Binding Protein), etc.] is counter-weighted by widespread negative correlations for CXCL9/CXCL10. Linear modelling further identifies a strong positive AXIN1–STAMBP relationship (R² = 0.7922, P < 0.0001; Figure 5G), indicating that homeostasis is maintained by multiple co-equal covariance axes. Upon progression to DM, FGF-21 dissociates from its original cluster and forms a novel “metabo-inflammatory” module together with IFN-γ, IL-6, CXCL10 and CDCP1, while CXCL9 merges into the other major cluster; the combined PC1 + PC2 variance drops to 50.98% (Figure 5B), signifying initial coupling of hyperglycaemia-driven metabolic signals with Th1 pathways. Figure 5E captures the accompanying immune-regulatory shift: IL-6 attains a negative correlation with TNF-like weak inducer of apoptosis (TWEAK) (r = −0.56, P = 0.037) and SCF–TRANCE exhibits a significant inverse regression (R² = 0.4338, P = 0.0104; Figure 5J), marking the first emergence of immunosuppressive cues. Once TB-DM develops, PCA resolves into three polarized clusters: FGF-21, IFN-γ, IL-6, CXCL9 and CXCL10 constitute an “inflammatory-metabolic” pole, whereas negative PC scores aggregate negative regulators such as TWEAK, FGF-19 and CASP-8; PC1 + PC2 variance now surges to 64.00% (Figure 5C), implying dominance by a few robust covariance axes. Network fragmentation becomes even more pronounced in Figure 5F: the AXIN1–STAMBP positive link is intensified (R² = 0.9640, P < 0.0001; Figure 5I), while strong negative correlations simultaneously appear for FGF-19–TNFSF14 (R² = 0.6311, P = 0.0105; Figure 5K) and CXCL9–TWEAK, generating a conspicuous pro- versus anti-inflammatory opposition. Thus, within a hyperglycaemic milieu, formerly scattered negative regulators are consolidated into an “inhibitory cluster” that continuously counteracts the “inflammatory cluster”, driving the paradoxical state of excessive inflammation co-existing with immune incompetence characteristic of TB-DM.
Figure 5. Multidimensional analysis of cytokine interactions. (A–C) Principal component analysis (PCA) score plots for HCs, DM, and TB-DM groups. Ellipses denote 95% confidence intervals. (D–F) Spearman correlation matrices (P <0.05). Red/blue hues indicate positive/negative correlations respectively, with intensity proportional to strength. (G–K) Univariate linear regression models of top-correlated cytokine pairs: (G) HCs: AXIN1 vs STAMBP (R²=0.7922). (H) DM: CASP-8 vs ST1A1 (R²=0.8322). (I) TB-DM: AXIN1 vs STAMBP (R²=0.9640). (J, K) Inverse correlations in DM (SCF vs TRANCE, R²=0.4338) and TB-DM (FGF-19 vs TNFSF14, R²=0.6311). P values: *<0.05, **<0.01.
3.6 Pathway involvements and molecular mechanisms of DEPs in TB-DM compared with HCs and DM
KEGG pathway analysis of DEPs revealed significant enrichment in rheumatoid arthritis-related pathways and cytokine/Toll-like receptor interaction pathways in TB-DM vs HCs comparisons (Table 3). This enrichment pattern suggests that the systemic inflammation in the TB-DM state shares molecular-level commonalities with chronic autoimmune diseases and robust innate immune responses. Notably, TB-DM vs DM comparisons exhibited overlapping pathway enrichments with TB-DM vs HCs, with additional involvement of intestinal immune network for IgA production and hematopoietic cell lineage pathways (Table 4). These specific pathways likely reflect a more extensive remodeling of immune homeostasis triggered by MTB infection in the context of diabetes, encompassing potential alterations in the mucosal immune barrier and hematopoietic development, the implications of which warrant further investigation.
Consistent with our previous findings linking TB-DM immunopathology to cytokine-mediated mechanisms (24), six common DEPs were identified in the hsa04060 (cytokine-cytokine receptor interaction) pathway across TB-DM vs HCs and TB-DM vs DM comparisons: CXCL5/CXCL6 (CXCR1 ligands), CXCL9/CXCL10 (CXCR3 ligands), IL-6 (IL-6R/IL6ST ligand), and LIGHT (TNFSF14; HVEM/DCR3 ligand) (Figure 6A). These molecules constitute a core hub network that drives Th1-type inflammation and the acute-phase response. Unique to TB-DM vs DM comparisons were CCL23 (CCR1 ligand), CCL15 (CCR3/CCR1 ligand), and IL5 (IL5RA/CSF2RB ligand) in the hsa04060 pathway (Figure 6B). These TB-DM-specific molecules may be linked to the activation of specific immune responses.
Figure 6. KEGG pathway enrichment of dysregulated cytokines. Cytokine-cytokine receptor interaction (hsa04060) pathway mapping: (A) TB-DM vs HCs: 3 proteins upregulated (red), 3 downregulated (blue). (B) TB-DM vs DM: 4 upregulated (red), 5 downregulated (blue). Gray rectangles denote membrane receptors.
GO molecular function annotation demonstrated predominant involvement of these proteins in cytokine activity, cytokine receptor binding, receptor-ligand activation, and signal transducer activity (Table 5). This functionally confirms their role as upstream signaling ligands. Subsequent pathway mapping revealed downstream activation of Jak-STAT signaling, MAPK signaling, ubiquitin-mediated proteolysis, actin cytoskeleton regulation, and leukocyte transendothelial migration (Figure 7). This series of pathways clearly delineates a complete chain from ligand-receptor binding to intracellular signal transduction, ultimately leading to a functional immune phenotype. Targeted analysis of the hsa04620 (Toll-like receptor signaling) pathway further linked Jak-STAT and NF-κB signaling to IP10 (CXCL10), MIG (CXCL9), and IL-6 production (Figure 8). This finding is crucial, as it suggests that in TB-DM, innate immune signals triggered by pathogen recognition and cytokine-mediated inflammatory signals may converge and mutually amplify through the two core pathways, Jak-STAT and NF-κB, thereby forming and sustaining a persistent and prominent inflammatory cycle.3
Figure 7. Signaling pathways involved in cytokine regulation in TB-DM comorbidity. This figure illustrates the signaling pathways implicated in cytokine regulation within the context of TB-DM comorbidity. The figure provides a detailed overview of the interactions between critical molecules within these pathways, including Gαi, Src, Ras, PI3K, AKT, and Ikk, which are involved in phosphorylation events (denoted by P) that regulate actin cytoskeleton reorganization, leukocyte transendothelial migration, and ubiquitin-mediated proteolysis. This diagram underscores the molecular mechanisms underlying the altered cytokine expression observed in TB-DM patients, as analyzed through the KEGG database pathway (hsa04062).
Figure 8. Regulation of IP-10 (CXCL10), MIG (CXCL9), and IL-6 by Jak-STAT and NF-κB signaling pathways. This diagram illustrates the activation of IP-10 (CXCL10), MIG (CXCL9), and IL-6 through the interplay of Jak-STAT and NF-κB signaling pathways in response to MTB stimuli. Jak-STAT Pathway: Initiated by INF-α/β binding to their receptors, leading to the activation of JAK1 and JAK2. These kinases phosphorylate STAT1 and STAT2, which form a complex with IRF9 to directly induce the expression of IP-10 and MIG. NF-κB Pathway: Activated via TLR2/4 recognition of MTB components, recruiting MyD88 and IRAKs (IRAK1/4), and subsequently TRAF6. This cascade activates TAK1, which phosphorylates IKKs (IKKα/β/γ), leading to IkBα degradation and NF-κB nuclear translocation. NF-κB further amplifies pro-inflammatory responses by upregulating IP-10, MIG, and IL-6. Phosphorylation events (denoted by “P”) and protein-protein interactions highlight critical regulatory steps. This dual-pathway activation underscores the coordinated immune response to MTB, with shared cytokines (e.g., IL-6) and cross-talk nodes (e.g., STAT1/NF-κB) driving inflammatory and antimicrobial activities in TB-DM comorbidity.
3.7 Diagnostic Performance of Single Cytokines and Their Combinations for Distinguishing TB-DM Patients from HCs and DM Patients
For single cytokine ROC curve analysis, IL-6, IFN-γ, and CXCL10 exhibited significant discriminatory abilities for distinguishing TB-DM patients from HCs and DM patients. In the TB-DM group compared to HCs, IL-6 achieved an AUC of 0.9841 (sensitivity 88.89%, specificity 92.86%) (Figure 9A). In the TB-DM group compared to DM, IL-6 achieved an AUC of 0.9206 (sensitivity 77.78%, specificity 85.71%) (Figure 9A). IFN-γ achieved an AUC of 0.9365 (sensitivity 88.89%, specificity 85.71%) in the TB-DM group compared to HCs and an AUC of 0.9524 (sensitivity 77.78%, specificity 92.86%) in the TB-DM group compared to DM (Figure 9B). CXCL10 achieved an AUC of 0.9286 (sensitivity 88.89%, specificity 85.71%) in the TB-DM group compared to HCs and an AUC of 0.881 (sensitivity 88.89%, specificity 71.43%) in the TB-DM group compared to DM (Figure 9C). Multivariable logistic regression analysis combining IL-6, IFN-γ, and CXCL10 as predictors achieved an AUC of 0.9841, sensitivity 88.89%, and specificity 92.86% (Figure 9D). Additionally, SCF achieved an AUC of 0.9921 (sensitivity 88.89%, specificity 92.86%) in the TB-DM group compared to HCs, but its discriminatory ability significantly decreased in the TB-DM group compared to DM (AUC = 0.7698) (Table 6).
Figure 9. ROC curve of cytokines distinguishing TB-DM patients from HC or DM patients. This figure presents the diagnostic performance of cytokines in distinguishing TB-DM patients from HCs and DM patients, as evaluated by ROC curves. (A) ROC of IL-6: IL-6 exhibited significant discriminatory ability between TB-DM and HCs (AUC = 0.9841) and between TB-DM and DM (AUC = 0.9206). (B) ROC of IFN-γ: IFN-γ showed strong diagnostic efficacy, with AUC values of 0.9365 (TB-DM vs HCs) and 0.9524 (TB-DM vs DM). (C) ROC of CXCL10: CXCL10 demonstrated robust diagnostic performance, achieving AUC values of 0.9286 (TB-DM vs HCs) and 0.881 (TB-DM vs DM). (D) Combined ROC of IL-6, IFN-γ, and CXCL10: The combination of these three cytokines achieved an AUC of 0.9841, highlighting their synergistic diagnostic potential in distinguishing TB-DM patients from DM patients.
Table 6. Sensitivity and specificity of IFN-γ, IL-6, SCF and CXCL10 cytokines in the differentiation of HCs, DM, and TB-DM patients.
4 Discussion
TB-DM has emerged as a significant public health challenge, with increasing incidence and severe complications. Previous studies have demonstrated that TB-DM patients exhibit enhanced immune system impairment, particularly in the context of chronic high blood sugar levels and insulin resistance (25, 26). This impairment compromises the body’s ability to combat MTB, leading to reduced immune response and increased disease progression (24). Consequently, identifying TB-DM-specific immune characteristics has become a critical area of research to better understand the disease’s pathophysiology and develop targeted therapeutic strategies.
Our comprehensive analysis of TB-DM patients revealed distinct immunological discrepancies. In DM patients, compared to HCs, NK cells showed a significant decline in absolute counts, while CD4+;T lymphocytes (a subset of adaptive immune cells) showed a marked increase. Notably, when disease progression advanced to TB-DM, these immunological abnormalities became more pronounced, with significant reductions in NK cells, total T lymphocytes, CD4+;T lymphocytes, CD8+;T lymphocytes and B lymphocytes across all three groups (TB-DM, DM, and HCs). This progressive and widespread reduction in immune cells may be associated with disorders of glucose metabolism, thereby leading to systemic immune dysfunction and potentially exacerbating disease progression (27, 28). However, the adaptive immune system’s partial recovery, particularly in CD4+;T lymphocytes, indicates a complex interplay between chronic inflammation and immune dysregulation in TB-DM.
Traditional perspectives often overlook the critical role of B lymphocytes in TB-DM, yet recent evidence highlights their importance in the disease’s pathogenesis. B lymphocytes play a pivotal role in both innate and adaptive immune responses, and their dysregulation has been implicated in the development of TB-DM (29, 30). Our study corroborated these findings by demonstrating that B lymphocytes exhibit significantly reduced absolute counts in TB-DM patients compared to HCs and DM patients, further supporting the idea that B lymphocytes are essential mediators of immune dysregulation in this disease.
Intercellular signaling plays a critical role in the recognition and activation of host immune cells (31). Immune cells regulate cytokine secretion or suppression through signaling pathways, and these cytokines may further amplify or inhibit cellular responses via positive or negative feedback mechanisms (32, 33). Previous studies suggest that TLR2- and TLR4-mediated immune responses are essential for defending against MTB invasion (34). Upon cytokine binding to TLRs, activation of the NF-κB signaling pathway stimulates DCs to secrete IL-6 and other cytokines, facilitating MTB clearance (35, 36). KEGG pathway enrichment analysis revealed enhanced IL-6 secretion via the NF-κB signaling pathway in TB-DM patients, potentially linked to hyperinflammatory responses in diabetes. This may promote neutrophil and macrophage recruitment/survival (37) and activate T-cell immunity to combat MTB. The JAK-STAT pathway was another enriched signaling axis. IFN-β binding to IFN-β receptors induces MIG (CXCL9) and IP10 (CXCL10) secretion through downstream signaling cascades. Our correlation matrix analysis demonstrated positive associations between CXCL9/CXCL10 and IFN-γ, with elevated levels of IFN-γ, CXCL9, and CXCL10 in TB-DM patients. Lande et al. reported that CXCL9 and CXCL10 are induced by IFN-β post-MTB infection (38). The CXCL9/CXCL10-CXCR3 axis exhibited significant co-regulation with the JAK-STAT pathway in TB-DM patients, suggesting its role as a metabolic-immune interface. However, individual cytokines (e.g., CXCL10) may be influenced by multiple pathways during MTB infection. NF-κB, JAK-STAT, and PI3K/Akt pathways synergistically activate CXCL10 transcription (39–41), yet only one pathway was observed in TB-DM patients, warranting further investigation into TB-DM-specific signaling mechanisms.
The identification of IL-6, IFN-γ, and CXCL10 as potential diagnostic biomarkers for TB-DM represents a significant breakthrough. These cytokines were found to exhibit marked elevation in TB-DM patients, with IL-6 levels showing a strong correlation with disease severity (42, 43). Furthermore, IL-6 has been implicated in chronic inflammation and metabolic dysregulation, particularly in the context of chronic granulomatous disease (CGD) and type 1 diabetes mellitus (T1DM) (44, 45). The increased expression of IL-6 in TB-DM patients may reflect its role in promoting immune cell migration and activation, thereby exacerbating disease progression (37, 46).
IFN-γ, another key cytokine in TB-DM, is primarily produced by activated NK cells and Th1-like T cells, and it plays a crucial role in regulating both innate and adaptive immunity (47–50). In the context of TB-DM, IFN-γ has been shown to downregulate the expression of CD4+;T lymphocytes and B lymphocytes, further supporting its role as a critical regulator of immune dysregulation (51, 52). Additionally, the upregulation of IFN-γ in TB-DM patients may reflect its potential role in mediating immune tolerance and resistance to MTB infection.
CXCL10, a member of the CXC family of chemotactic proteins (53), has also been implicated in the pathogenesis of TB-DM (54). Its increased expression in TB-DM patients is associated with chronic inflammation and metabolic derangements, suggesting a potential role in the disease’s progression (55). Furthermore, CXCL10 has been shown to promote the migration and activation of immune cells, particularly in the context of chronic infections (56–58). These findings highlight the importance of CXCL10 in the immune dysregulation observed in TB-DM patients (59).
The diagnostic potential of IL-6, IFN-γ, and CXCL10 was further supported by our ROC analysis. The combination of these three biomarkers achieved a diagnostic accuracy of 0.9841 (sensitivity 88.89%, specificity 92.86%), demonstrating their ability to effectively distinguish TB-DM patients from healthy controls and DM patients. This multi-dimensional approach offers a robust strategy for the early diagnosis and clinical management of TB-DM, potentially improving patient outcomes.
Despite the significant advances in understanding TB-DM’s immune characteristics, several limitations remain. First, the relatively small sample size may limit the generalizability of our findings. Second, while we have successfully identified potential biomarkers, further functional studies and animal experiments are required to clarify the causal relationship between these markers and disease progression. Additionally, the molecular mechanisms underlying the observed immune dysregulation in TB-DM remain poorly understood, necessitating further research to unravel the complex interplay between cytokines, chemokines, and immune cells in this disease. Lastly, this study did not analyze the potential impact of diabetes duration on immunological features; future prospective studies should incorporate this important variable into their design.
5 Conclusions
This study conducted a comprehensive analysis of the immunological features in TB-DM patients, revealing distinct immunological abnormalities characteristic of the disease. The findings indicate that TB-DM patients exhibit significantly reduced absolute counts of NK cells, CD4+;T cells, and B cells in their peripheral blood, with these immune cell populations becoming increasingly dysfunctional during disease progression. Based on the Olink multiplexed target-specific protein detection technology, we successfully identified IL-6, IFN-γ, and CXCL10 as potential diagnostic biomarkers for TB-DM, with their expression levels significantly elevated in TB-DM patients and correlating with disease severity. These results highlight the critical role of these cytokines in the early diagnosis and clinical management of TB-DM, offering a robust diagnostic tool for early intervention.
Through KEGG and GO pathway analysis, we identified that IL-6, IFN-γ, and CXCL10 are enriched in NF-κB, JAK-STAT, and CXC subfamily signaling pathways, suggesting their involvement in TB-DM pathogenesis through multiple immune regulation pathways. Combined with ROC curve analysis, we demonstrated that these biomarkers exhibit high sensitivity (88.89%) and specificity (92.86%), effectively distinguishing TB-DM patients from DM patients. This multi-dimensional diagnostic approach provides a highly sensitive and specific method for the precise identification of TB-DM cases, potentially improving patient outcomes.
Data availability statement
All data generated or analyzed during this study are included in this published article.
Ethics statement
The studies involving humans were approved by Medical Ethics Committee of the Eighth Medical Center of PLA General Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
YuL: Methodology, Data curation, Software, Writing – original draft, Visualization. PC: Writing – original draft, Software, Data curation, Methodology. YA: Methodology, Writing – original draft. RN: Writing – original draft, Methodology. ZY: Methodology, Writing – original draft. LingY: Methodology, Writing – original draft. LZ: Writing – original draft, Methodology. LL: Methodology, Writing – original draft. LW: Writing – review & editing, Conceptualization. WG: Visualization, Conceptualization, Project administration, Writing – review & editing, Supervision.
Funding
The author(s) declared that financial support was not received for this work and/or its publication.
Acknowledgments
We sincerely thank the reviewers and editors who provided review and editing services for this study.
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2025.1713654/full#supplementary-material
Supplementary Table 1 | The detailed information of Cohort I.
Supplementary Table 2 | The detailed information of Cohort II.
Abbreviations
ADA, American Diabetes Association; ATB, active tuberculosis; DEPs, Differentially expressed proteins; DM, diabetes mellitus; GO, Gene Ontology; HCs, healthy controls; IDF, International Diabetes Federation; IGRA, interferon-γ release assay; KEGG, Kyoto Encyclopedia of Genes and Genomes; MTB, Mycobacterium tuberculosis; NPX, Normalized Protein Expression; PCA, Principal component analysis; PEA, Proximity Extension Assay; TB, Tuberculosis; TB-DM, tuberculosis-diabetes comorbidity; WHO, World Health Organization. AXIN1, Axis inhibition protein 1; CASP-8, Caspase-8; CD6, Cluster of Differentiation 6; CD244, CD244 molecule; CDCP1, CUB domain-containing protein 1; CXCL10, C-X-C motif chemokine ligand 10; CXCL5, C-X-C motif chemokine ligand 5; CXCL6, C-X-C motif chemokine ligand 6; CXCL9, C-X-C motif chemokine ligand 9; DNER, Delta/Notch-like EGF repeat-containing; EN-RAGE, Extracellular newly identified receptor for advanced glycation end-products binding protein; FGF-19, Fibroblast growth factor 19; FGF-21, Fibroblast growth factor 21; IFN-γ, Interferon-gamma; IL-6, Interleukin-6; MCP-4, Monocyte chemoattractant protein-4; 4E-BP1, 4E-binding protein 1; SCF, Stem cell factor; ST1A1, Sulfotransferase family 1A member 1; STAMBP, STAM-binding protein; TNFSF14, Tumor necrosis factor superfamily member 14; TRANCE, TNF-related activation-induced cytokine; TWEAK, TNF-like weak inducer of apoptosis.
References
1. Cardona PJ, Català M, and Prats C. Origin of tuberculosis in the Paleolithic predicts unprecedented population growth and female resistance. Sci Rep. (2020) 10:42. doi: 10.1038/s41598-019-56769-1
2. Liu Y, Yang L, Meskini M, Goel A, Opperman M, Shyamal SS, et al. Gut microbiota and tuberculosis. IMeta. (2025) 4:e70054. doi: 10.1002/imt2.70054
4. Duncan BB, Magliano DJ, and Boyko EJ. IDF diabetes atlas 11th edition 2025: global prevalence and projections for 2050. Nephrol Dial Transplant. (2026) 41(1), 7–9. Online ahead of print. doi: 10.1093/ndt/gfaf177
5. Shetty S, Pappachan JM, and Fernandez CJ. Diabetes and tuberculosis: An emerging dual threat to healthcare. World J Diabetes. (2024) 15:1409–16. doi: 10.4239/wjd.v15.i7.1409
6. Ye Z, Bai G, Yang L, Zhuang L, Li L, Li Y, et al. Development and validation of early alert model for diabetes mellitus–tuberculosis comorbidity. Microorganisms. (2025) 13:919. doi: 10.3390/microorganisms13040919
7. Li Y, Ni R, An Y, Yang L, Ye Z, Zhuang L, et al. Bibliometric analysis and knowledge mapping of diabetes mellitus combined with tuberculosis research: trends from 1995 to 2023. Front Immunol. (2025) 16:1571123. doi: 10.3389/fimmu.2025.1571123
8. Ye Z, Li L, Yang L, Zhuang L, Aspatwar A, Wang L, et al. Impact of diabetes mellitus on tuberculosis prevention, diagnosis, and treatment from an immunologic perspective. Explor (Beijing). (2024) 4:20230138. doi: 10.1002/EXP.20230138
9. Dooley KE and Chaisson RE. Tuberculosis and diabetes mellitus: convergence of two epidemics. Lancet Infect Dis. (2009) 9:737–46. doi: 10.1016/s1473-3099(09)70282-8
10. Vinhaes CL, Fukutani ER, Santana GC, Arriaga MB, Barreto-Duarte B, Araújo-Pereira M, et al. An integrative multi-omics approach to characterize interactions between tuberculosis and diabetes mellitus. iScience. (2024) 27:109135. doi: 10.1016/j.isci.2024.109135
11. Rocha EF, Vinhaes CL, Araújo-Pereira M, Mota TF, Gupte AN, Kumar NP, et al. The sound of silent RNA in tuberculosis and the lncRNA role on infection. iScience. (2024) 27:108662. doi: 10.1016/j.isci.2023.108662
12. Al-Rifai RH, Pearson F, Critchley JA, and Abu-Raddad LJ. Association between diabetes mellitus and active tuberculosis: A systematic review and meta-analysis. PLoS One. (2017) 12:e0187967. doi: 10.1371/journal.pone.0187967
13. Li M, Chen T, Hua Z, Yan H, Wang D, Li Z, et al. Global, regional, and national prevalence of diabetes mellitus in patients with pulmonary tuberculosis: a systematic review and meta-analysis. Diabetol Metab Syndr. (2021) 13:127. doi: 10.1186/s13098-021-00743-3
14. Jeon CY and Murray MB. Diabetes mellitus increases the risk of active tuberculosis: a systematic review of 13 observational studies. PLoS Med. (2008) 5:e152. doi: 10.1371/journal.pmed.0050152
15. Boadu AA, Yeboah-Manu M, Osei-Wusu S, and Yeboah-Manu D. Tuberculosis and diabetes mellitus: The complexity of the comorbid interactions. Int J Infect Dis. (2024) 146:107140. doi: 10.1016/j.ijid.2024.107140
16. Moreno-Martínez A, Casals M, Orcau À., Gorrindo P, Masdeu E, and Caylà JA. Factors associated with diabetes mellitus among adults with tuberculosis in a large European city, 2000-2013. Int J Tuberc Lung Dis. (2015) 19:1507–12. doi: 10.5588/ijtld.15.0102
17. Queiroz ATL, Vinhaes CL, Fukutani ER, Gupte AN, Kumar NP, Fukutani KF, et al. A multi-center, prospective cohort study of whole blood gene expression in the tuberculosis-diabetes interaction. Sci Rep. (2023) 13:7769. doi: 10.1038/s41598-023-34847-9
18. Baker MA, Harries AD, Jeon CY, Hart JE, Kapur A, Lönnroth K, et al. The impact of diabetes on tuberculosis treatment outcomes: a systematic review. BMC Med. (2011) 9:81. doi: 10.1186/1741-7015-9-81
19. ADAPPC. 2. Diagnosis and classification of diabetes: standards of care in diabetes-2025. Diabetes Care. 48(Supplement_1):S27–49. doi: 10.2337/dc25-S002
20. Wik L, Nordberg N, Broberg J, Björkesten J, Assarsson E, Henriksson S, et al. Proximity extension assay in combination with next-generation sequencing for high-throughput proteome-wide analysis. Mol Cell Proteomics. (2021) 20:100168. doi: 10.1016/j.mcpro.2021.100168
21. Kanehisa M, Furumichi M, Sato Y, Kawashima M, and Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. (2023) 51:D587–d592. doi: 10.1093/nar/gkac963
22. Consortium, G.O. The Gene Ontology resource: enriching a GOld mine. Nucleic Acids Res. (2021) 49:D325–d334. doi: 10.1093/nar/gkaa1113
23. Consortium, G.O. Expansion of the Gene Ontology knowledgebase and resources. Nucleic Acids Res. (2017) 45:D331–d338. doi: 10.1093/nar/gkw1108
24. Cheng P, Wang L, and Gong W. Cellular immunity of patients with tuberculosis combined with diabetes. J Immunol Res. (2022) 2022:6837745. doi: 10.1155/2022/6837745
25. Panda S, Seelan DM, Faisal S, Arora A, Luthra K, Palanichamy JK, et al. Chronic hyperglycemia drives alterations in macrophage effector function in pulmonary tuberculosis. Eur J Immunol. (2022) 52:1595–609. doi: 10.1002/eji.202249839
26. Restrepo BI, Scordo JM, Aguillón-Durán GP, Ayala D, Quirino-Cerrillo AP, Loera-Salazar R, et al. Differential role of type 2 diabetes as a risk factor for tuberculosis in the elderly versus younger adults. Pathogens. (2022) 11:1551. doi: 10.3390/pathogens11121551
27. Yoon Kim D and Kwon Lee J. Type 1 and 2 diabetes are associated with reduced natural killer cell cytotoxicity. Cell Immunol. (2022) 379:104578. doi: 10.1016/j.cellimm.2022.104578
28. Kumar NP, Moideen K, Dolla C, Kumaran P, and Babu S. Prediabetes is associated with the modulation of antigen-specific Th1/Tc1 and Th17/Tc17 responses in latent Mycobacterium tuberculosis infection. PLoS One. (2017) 12:e0178000. doi: 10.1371/journal.pone.0178000
29. Kumar NP, Moideen K, Dhakshinraj SD, Banurekha VV, Nair D, Dolla C, et al. Profiling leucocyte subsets in tuberculosis-diabetes co-morbidity. Immunology. (2015) 146:243–50. doi: 10.1111/imm.12496
30. Gong W, Liang Y, and Wu X. The current status, challenges, and future developments of new tuberculosis vaccines. Hum Vaccin Immunother. (2018) 14:1697–716. doi: 10.1080/21645515.2018.1458806
31. Kumar V. Inflammation research sails through the sea of immunology to reach immunometabolism. Int Immunopharmacol. (2019) 73:128–45. doi: 10.1016/j.intimp.2019.05.002
32. Anders HJ. Of inflammasomes and alarmins: IL-1β and IL-1α in kidney disease. J Am Soc Nephrol. (2016) 27:2564–75. doi: 10.1681/asn.2016020177
33. Lee CH and Chun T. Anti-inflammatory role of TAM family of receptor tyrosine kinases via modulating macrophage function. Mol Cells. (2019) 42:1–7. doi: 10.14348/molcells.2018.0419
34. Wani BA, Shehjar F, Shah S, Koul A, Yusuf A, Farooq M, et al. Role of genetic variants of Vitamin D receptor, Toll-like receptor 2 and Toll-like receptor 4 in extrapulmonary tuberculosis. Microb Pathog. (2021) 156:104911. doi: 10.1016/j.micpath.2021.104911
35. Akira S, Uematsu S, and Takeuchi O. Pathogen recognition and innate immunity. Cell. (2006) 124:783–801. doi: 10.1016/j.cell.2006.02.015
36. O’Neill LA, Golenbock D, and Bowie AG. The history of Toll-like receptors - redefining innate immunity. Nat Rev Immunol. (2013) 13:453–60. doi: 10.1038/nri3446
37. Gupte AN, Kumar P, Araújo-Pereira M, Kulkarni V, Paradkar M, Pradhan N, et al. Baseline IL-6 is a biomarker for unfavourable tuberculosis treatment outcomes: a multisite discovery and validation study. Eur Respir J. (2022) 59:2100905. doi: 10.1183/13993003.00905-2021
38. Lande R, Giacomini E, Grassi T, Remoli ME, Iona E, Miettinen M, et al. IFN-alpha beta released by Mycobacterium tuberculosis-infected human dendritic cells induces the expression of CXCL10: selective recruitment of NK and activated T cells. J Immunol. (2003) 170:1174–82. doi: 10.4049/jimmunol.170.3.1174
39. Aota K, Yamanoi T, Kani K, Ono S, Momota Y, and Azuma M. Inhibition of JAK-STAT signaling by baricitinib reduces interferon-γ-induced CXCL10 production in human salivary gland ductal cells. Inflammation. (2021) 44:206–16. doi: 10.1007/s10753-020-01322-w
40. Méndez-Samperio P, Pérez A, and Rivera L. Mycobacterium bovis Bacillus Calmette-Guérin (BCG)-induced activation of PI3K/Akt and NF-kB signaling pathways regulates expression of CXCL10 in epithelial cells. Cell Immunol. (2009) 256:12–8. doi: 10.1016/j.cellimm.2008.12.002
41. Zaheer RS, Koetzler R, Holden NS, Wiehler S, and Proud D. Selective transcriptional down-regulation of human rhinovirus-induced production of CXCL10 from airway epithelial cells via the MEK1 pathway. J Immunol. (2009) 182:4854–64. doi: 10.4049/jimmunol.0802401
42. Singh K, Hussain T, Gupta B, and Pati SA. Comparative investigation on cytokine expression in pulmonary tuberculosis and comorbidity with type 2 diabetes mellitus. Int J Mycobacteriol. (2024) 13:165–70. doi: 10.4103/ijmy.ijmy_40_24
43. Fève B and Bastard JP. The role of interleukins in insulin resistance and type 2 diabetes mellitus. Nat Rev Endocrinol. (2009) 5:305–11. doi: 10.1038/nrendo.2009.62
44. Choi E, Choi HH, Kwon KW, Kim H, Ryu JH, Hong JJ, et al. Permissive lung neutrophils facilitate tuberculosis immunopathogenesis in male phagocyte NADPH oxidase-deficient mice. PLoS Pathog. (2024) 20:e1012500. doi: 10.1371/journal.ppat.1012500
45. Mohamed AA, Abdallah GM, Ibrahim IT, Ali NS, Hussein MA, Thabet GM, et al. Evaluation of miRNA-146a, miRNA-34a, and pro-inflammatory cytokines as a potential early indicators for type 1 diabetes mellitus. Noncoding RNA Res. (2024) 9:1249–56. doi: 10.1016/j.ncrna.2024.06.001
46. Matsumura F, Polz R, Singh S, Matsumura A, Scheller J, and Yamashiro S. Investigation of fascin1, a marker of mature dendritic cells, reveals a new role for IL-6 signaling in CCR7-mediated chemotaxis. J Immunol. (2021) 207:938–49. doi: 10.4049/jimmunol.2000318
47. Young HA and Bream JH. IFN-gamma: recent advances in understanding regulation of expression, biological functions, and clinical applications. Curr Top Microbiol Immunol. (2007) 316:97–117. doi: 10.1007/978-3-540-71329-6_6
48. Zhuang L, Yang L, Li L, Ye Z, and Gong W. Mycobacterium tuberculosis: immune response, biomarkers, and therapeutic intervention. MedComm (2020). (2024) 5:e419. doi: 10.1002/mco2.419
49. Gong W, Pan C, Cheng P, Wang J, Zhao G, and Wu X. Peptide-based vaccines for tuberculosis. Front Immunol. (2022) 13:830497. doi: 10.3389/fimmu.2022.830497
50. Cheng P, Wang L, and Gong W. In silico analysis of peptide-based biomarkers for the diagnosis and prevention of latent tuberculosis infection. Front Microbiol. (2022) 13:947852. doi: 10.3389/fmicb.2022.947852
51. Wei R, Li P, Xue Y, Liu Y, Gong W, and Zhao W. Impact of diabetes mellitus on the immunity of tuberculosis patients: A retrospective, cross-sectional study. Risk Manag Healthc Policy. (2022) 15:611–27. doi: 10.2147/RMHP.S354377
52. Francisco Neto A, Dell’Armelina Rocha PR, Perez EC, Xavier JG, Peres GB, Spadacci-Morena DD, et al. Diabetes mellitus increases the susceptibility to encephalitozoonosis in mice. PLoS One. (2017) 12:e0186954. doi: 10.1371/journal.pone.0186954
53. Li H, Sun Y, and Sun LA. Teleost CXCL10 is both an immunoregulator and an antimicrobial. Front Immunol. (2022) 13:823. doi: 10.3389/fimmu.2022.917697
54. Kochumon S, Al-Sayyar A, Jacob T, Arefanian H, Bahman F, Almansour N, et al. IL-1β-induced CXCL10 expression in THP-1 monocytic cells involves the JNK/c-jun and NF-κB-mediated signaling. Pharm (Basel). (2024) 17. doi: 10.3390/ph17070823
55. Rajamanickam A, Kothandaraman SP, Kumar NP, Viswanathan V, Shanmugam S, Hissar S, et al. Cytokine and chemokine profiles in pulmonary tuberculosis with pre-diabetes. Front Immunol. (2024) 15:1447161. doi: 10.3389/fimmu.2024.1447161
56. Papaioannou S, See JX, Jeong M, de la Torre C, Ast V, Reiners-Koch PS, et al. Liver sinusoidal endothelial cells orchestrate NK cell recruitment and activation in acute inflammatory liver injury. Cell Rep. (2023) 42:112836. doi: 10.1016/j.celrep.2023.112836
57. Yue SY, Niu D, Ma WM, Guan Y, Liu QS, Wang XB, et al. The CXCL10/CXCR3 axis regulates Th1 cell differentiation and migration in experimental autoimmune prostatitis through the PI3K/AKT pathway. Andrology. (2024) 12:1408–18. doi: 10.1111/andr.13571
58. Dillemans L, Yu K, De Zutter A, Noppen S, Gouwy M, Berghmans N, et al. Natural carboxyterminal truncation of human CXCL10 attenuates glycosaminoglycan binding, CXCR3A signaling and lymphocyte chemotaxis, while retaining angiostatic activity. Cell Commun Signal. (2024) 22:94. doi: 10.1186/s12964-023-01453-1
Keywords: cytokines, diabetes mellitus, lymphocyte subsets, proteomics, tuberculosis, tuberculosis-diabetes comorbidity (TB-DM)
Citation: Li Y, Cheng P, An Y, Ni R, Ye Z, Yang L, Zhuang L, Li L, Wang L and Gong W (2026) Identification of biomarkers and construction of discriminating model for tuberculosis patients with diabetes mellitus based on proteomics: a cross-sectional study. Front. Immunol. 16:1713654. doi: 10.3389/fimmu.2025.1713654
Received: 26 September 2025; Accepted: 22 December 2025; Revised: 12 December 2025;
Published: 15 January 2026.
Edited by:
Jianping Xie, Southwest University, ChinaReviewed by:
Caian Vinhaes, Clinical Research Laboratory on Mycobacteria - National Institute of Infectious Diseases Evandro Chagas - Fiocruz, BrazilJing Zhu, China Academy of Engineering Physics, China
Copyright © 2026 Li, Cheng, An, Ni, Ye, Yang, Zhuang, Li, Wang and Gong. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Liang Wang, V2FuZ2wzMDlAc2luYS5jb20=; Wenping Gong, Ph.D., Z3dwODkxMDE1QHdodS5lZHUuY24=
†These authors have contributed equally to this work