Population characteristics
A total of 88 cases were included in this study, including 20 cases in the healthy group (HC), 20 cases with latent TB (LTBI), 28 cases with drug-sensitive TB (DS-TB), and 20 cases with drug-resistant TB (MDR-TB), 20 cases with severe pulmonary TB, 28 cases with non-severe pulmonary TB. The average age was 44.21±16.722, and there was no statistically significant difference among the groups (F=2.232, p=0.090). There were 49 males and 41 females, showing no statistically significant difference among the groups (H=5.680, p=0.123). There were no significant differences in the basic characteristics between the DS-TB and MDR-TB groups, except for elevated levels of D-Dimer (DD) and serum creatinine (CREA) in the DS-TB group relative to the MDR-TB group (Table 1).
Table 1 Population Demographic data and baseline characteristics
Characteristics
|
HC
(n=20)
|
LTBI
(n=20)
|
TB
|
p
value
|
DS-TB group
(n=28)
|
MDR-TB group(n=20)
|
Gender,n(%)
|
|
|
|
|
0.428
|
male
|
11(55.0)
|
10(45.5)
|
15(53.6)
|
13(65.0)
|
|
female
|
9(45.0)
|
12(54.5)
|
13(46.4)
|
7(35.0)
|
|
Age(years)
|
32.00
(26.25,46.75)
|
44.50
(39.75,52.75)
|
44.50(29.25,60.75)
|
38.00(27.75,58.00)
|
0.530
|
Classification of TB,n(%)
|
-
|
-
|
|
|
0.692
|
Non-severe pulmonary TB
|
-
|
-
|
17(60.7)
|
11(55.0)
|
|
severe pulmonary TB
|
|
|
11(39.3)
|
9(45.0)
|
|
Smear,n(%)
|
|
|
|
|
0.843
|
negative
|
|
|
12(42.9)
|
8(40.0)
|
|
positive
|
|
|
16(57.1)
|
12(60.0)
|
|
WBC
|
-
|
-
|
7.81±4.27
|
7.76±3.55
|
0.968
|
Hemoglobin
|
-
|
-
|
107.86±24.00
|
113.45±24.67
|
0.435
|
PLT
|
-
|
-
|
298.50±99.95
|
265.75±119.73
|
0.308
|
Lymphocyte
|
-
|
-
|
1.13±0.57
|
3.37±7.48
|
0.195
|
Monocyte
|
-
|
-
|
0.48±0.24
|
0.70±0.61
|
0.093
|
Neutrophil
|
-
|
-
|
6.10±4.32
|
10.21±18.54
|
0.262
|
Albumin
|
-
|
-
|
12.86±8.64
|
21.15±32.12
|
0.273
|
Total bilirubin
|
-
|
-
|
11.19±5.70
|
10.98±3.98
|
0.886
|
CREA
|
-
|
-
|
49.31±14.96
|
64.08±25.74
|
0.016
|
LDH
|
-
|
-
|
159.00
(137.25,219.25)
|
144.00
(131.50,172.00)
|
0.143
|
CRP
|
-
|
-
|
44.18(8.03,73.04)
|
15.72(0.76,55.07)
|
0.094
|
PCT
|
-
|
-
|
0.04(0.02,0.14)
|
0.03(0.02,0.08)
|
0.718
|
NT-proBNP
|
-
|
-
|
619.70
(421.25,1744.75)
|
339.05
(66.90,1439.20)
|
0.102
|
PH
|
-
|
-
|
7.43(7.41,7.46)
|
7.45(7.44,7.46)
|
0.242
|
PaCO2
|
-
|
-
|
37.00(34.50,41.50)
|
42.00(35.00,54.25)
|
0.138
|
PaO2
|
-
|
-
|
100.62±24.87
|
123.86±52.28
|
0.140
|
DD
|
-
|
-
|
1.70±1.46
|
0.95±0.90
|
0.040
|
NLR
|
-
|
-
|
4.02(1.84,11.65)
|
3.78(2.43,8.62)
|
0.707
|
WBC white blood cell, N neutrophil, Hb hemoglobin, PLT platelet, CRP C-reactive protein, CREA serum creatinine (CREA).
Note: a represents using independent samples t-test, β represents using Mann-Whitney U test, c represents using chi-square test, d represents using Fisher's exact test.
Plasma Protein Profiling and PCA Analysis: Identifying Differential Expression and Potential Biomarkers in various subtypes of tuberculosis
All the 92 plasma proteins were identified in the three groups, and 7 of them differed significantly between HC and LTBI groups, 43 proteins differed significantly between LTBI and TB groups, 46 proteins differed significantly between TB and HC groups (Figure 2A). Principal component analysis (PCA) analysis revealed clear distinctions between LTBI and TB (Figure2B), as well as MDR-TB and DS-TB (Figure 2C). CXCL10, SCF, and TRANCE were the main contributors to variability in Dim1 (36.0%), while PDL-1 and CXCL11 were most important in Dim2 (29.1%). TRANCE and TWEAK were the key factors in the variability among MDR-TB in Dim1 (45.5%). Overall CXCL10 and TGF-alpha proteins differed significantly among the three groups which could be used as potential diagnostic biomarkers (Figure 2D).
Comparative Proteomic Analysis: Differential Expression and Diagnostic Markers in various subtypes of tuberculosis
In the comparison between the active TB group and the HC group, 46 differentially expressed proteins were identified. Compared to the HC group, IL-6, CXCL9, CXCL10, IFN-gamma, EN-RAGE, and MCP-3 were significantly elevated in active tuberculosis, while TRANCE, TWEAK, SCF, and TRAIL were notably decreased (Figure3).
We found that 7 inflammatory proteins including CCL23, CCL28, CXCL10, NT-3, IL-12B, IL-17A, TGF-alpha can potentially distinguish LTBI from HC groups. All the 7 proteins remained significantly different after statistical adjustments. Nevertheless, a combination of these 7 protein biomarkers showed significantly improved diagnostic value than the individuals.
Similarly, 43 inflammatory proteins can effectively distinguished TB from LTBI groups. Among them, the expression of IL-6, IFN-γ, EN-RAGE were higher in TB, while the expression of SCF (Stem cell factor ), DNER(Delta and Notch-like epidermal growth factor-related receptor (DNER), TRAIL (TNF-related apoptosis-inducing ligand) proteins were lower in LTBI (Figure 3).
The detection of differentially expressed inflammatory proteins among the three groups highlights the complexity of the immune response in TB and the importance of these proteins in disease diagnosis, progression, and treatment. This founding could be instrumental in advancing our understanding of TB and in the development of more effective diagnostic and therapeutic strategies. Further research is needed to fully understand the roles of these proteins in TB pathogenesis and to explore their potential clinical applications.
Identification and Diagnostic Potential of Inflammatory Proteins in MDR-TB, DS-TB, LTBI
In the comparative analysis between the Multi-Drug Resistant Tuberculosis (MDR-TB) and Drug-Sensitive Tuberculosis (DS-TB) groups, six proteins were found to be differentially expressed. Specifically, the levels of IL-2RB, GDNF, CST5, TRANCE, TWEAK, and IL10RA were significantly reduced in the MDR-TB group compared to the DS-TB group (Figure4). Among these, IL-2RB and TRANCE, which are inflammatory proteins, demonstrated moderate predictive power in distinguishing MDR-TB from DS-TB, with their respective area under the curve (AUC) values being 0.709.
In contrast, when assessing the diagnostic value for differentiating between Latent Tuberculosis Infection (LTBI) and Healthy Controls (HC), CCL28 (AUC = 0.677), TGF-alpha (AUC = 0.682), and NT-3 (AUC = 0.665) exhibited low predictive capabilities. On the other hand, CXCL9 (AUC = 0.843), IFN-alpha (AUC = 0.843), and EN-RAGE (AUC = 0.837) showed good diagnostic value when differentiating between Active Tuberculosis (ATB) and HC. These proteins indicate a stronger inflammatory response and immune activation in ATB patients compared to healthy individuals.
Furthermore, SCF (AUC = 0.921), IFN-alpha (AUC = 0.902), and EN-RAGE (AUC = 0.882) displayed superior diagnostic value for distinguishing between LTBI and ATB. This highlights the potential of these proteins as biomarkers to identify the transition from a latent to an active state of the disease (Figure5).
However, the differential expression of specific inflammatory-related proteins holds promise as diagnostic biomarkers for various stages and severities of tuberculosis. The study underscores the significance of cytokine and chemokine dysregulation in the progression of the disease, as indicated by the AUC values in Figure 4. These findings could contribute to a better understanding of TB pathogenesis and aid in the development of more accurate diagnostic tools.
Identification and Diagnostic Potential of Inflammatory Proteins in negative tuberculosis and positive pulmonary TB
In this study, we conducted an analysis of inflammatory protein expression in individuals with negative and positive tuberculosis, identifying 6 inflammatory proteins that may serve as potential markers for distinguishing between the two groups. Of these proteins, SLAMF1, SCF, LIF, NRTN, EN-RAGE, and FGF-19 exhibited statistically significant differences (Figure 6). Furthermore, SLAMF1(AUC=0.779) and MMP (0.712) showed a significant diagnostic value for negative TB. Nevertheless, a combination of these 2 protein biomarkers showed significantly improved diagnostic value than the individuals (Figure 7).
Identification and Diagnostic Potential of Inflammatory Proteins in severe pulmonary tuberculosis
In our study, we have identified 43 inflammatory proteins that exhibit significant potential for distinguishing between different severities of pulmonary tuberculosis. Specifically, we observed a decrease in the expression of SCF (Stem Cell Factor), MCP-4 (Monocyte Chemoattractant Protein-4), and TRANCE (Tumor Necrosis Factor-Related Activation-Induced Cytokine) compared to non-severe cases. Conversely, the expression levels of the remaining 40 inflammatory proteins were found to be increased in severe cases (Figure 8). This observation is indicative of a possible role these proteins play in the progression of the disease.
Furthermore, IL6, EN-RAGE, CXCL10, and PD-L1 demonstrated a notable diagnostic efficacy, the efficacy was determined by their high area under the curve (AUC) values, which exceeded 0.800, a threshold often used to indicate strong predictive performance in diagnostic tests. This suggests that the differential expression of specific inflammatory proteins can be a significant indicator of disease severity in pulmonary tuberculosis. The identification of these proteins, particularly those with high diagnostic efficacy, may contribute to the development of more precise diagnostic tools and personalized treatment strategies for patients with severe pulmonary tuberculosis (Table 2) .
Table 2 The diagnostic value of potential protein biomarkers between Non-severe pulmonary TB and severe pulmonary TB
Name
|
Cutoff1
|
Sensitivity
|
Specificity
|
AUC
|
95%CI
|
p_values
|
IL6
|
0.521
|
0.850
|
0.926
|
0.933
|
0.864-1
|
8.85E-09
|
EN-RAGE
|
0.382
|
0.950
|
0.778
|
0.915
|
0.822-1
|
7.99E-08
|
CXCL10
|
0.506
|
0.850
|
0.889
|
0.907
|
0.817-0.998
|
3.37E-08
|
PD-L1
|
0.552
|
0.800
|
0.926
|
0.904
|
0.818-0.989
|
2.77E-07
|
OSM
|
0.474
|
0.800
|
0.852
|
0.891
|
0.798-0.983
|
4.31E-07
|
IL-18R1
|
0.361
|
0.950
|
0.742
|
0.889
|
0.797-0.981
|
6.87E-07
|
CXCL11
|
0.381
|
0.950
|
0.778
|
0.881
|
0.782-0.981
|
1.15E-06
|
CXCL9
|
0.328
|
0.950
|
0.741
|
0.881
|
0.782-0.981
|
1.16E-06
|
CCL20
|
0.493
|
0.750
|
0.926
|
0.880
|
0.779-0.98
|
3.18E-06
|
IL18
|
0.468
|
0.900
|
0.815
|
0.865
|
0.76-0.97
|
8.51E-06
|
CSF-1
|
0.462
|
0.800
|
0.852
|
0.857
|
0.748-0.967
|
1.92E-05
|
SCF
|
0.314
|
0.900
|
0.741
|
0.852
|
0.74-0.964
|
2.70E-06
|
Correlation Analysis of Differentially Expressed Inflammatory Proteins in various subtypes of tuberculosis.
To further investigation the correlation of inflammatory proteins within each group, a clustered heat map was generated by computing the Pearson correlation coefficients among pairwise differentially expressed proteins. The X and Y axes of the graph denote the names of these differentially expressed proteins, with red denoting positive correlation and blue denoting negative correlation. The intensity of the color reflects the strength of the correlation.
The co-differential proteins in individuals with TB and HC exhibit a robust positive correlation between SIRT2 and STAMBP, and a negative correlation between SCF and IL-6. Conversely, the co-differential proteins in individuals with LTBI and HC demonstrate a strong positive correlation between TGF-alpha and CXCL10, with a correlation coefficient of 0.638. Furthermore, in individuals with LTBI and TB, a positive correlation is observed between CXCL9 and CXCL10, as well as STAMBP and CASP. The co-differential proteins in non-severe and severe tuberculosis exhibit a strong positive correlation for CXCL10 and CXCL9 (correlation coefficient = 0.920) and TNF and CCL3 (correlation coefficient = 0.924), as well as a negative correlation for IL-6 and SCF (correlation coefficient = -0.735). Conversely, the co-differential proteins in negative and positive tuberculosis demonstrate a significant negative correlation for SCF and FGF 19 (correlation coefficient = -0.6184). Additionally, the co-differential proteins in DS-TB and MDR-TB show a significant positive correlation for TRANCE and TWEAK (correlation coefficient = 0.6273). These findings indicate that inflammation-related proteins are intercorrelated (Figure 9).
GO and KEGG Analysis of Co-Differentially Expressed Proteins in various subtypes of tuberculosis
To confirm the role of co-differentially expressed proteins within three distinct groups, we conducted Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis on these proteins. In the biological process (BP), co-differentially expressed proteins between HC and TB are primarily enriched in processes related to immune regulation, including positive regulation of interleukin-1 beta production, inflammatory response, lymphocyte chemotaxis, T cell proliferation, humoral immune response, macrophage differentiation, peptidyl-tyrosine and peptidyl-serine phosphorylation, MAP kinase activity, chemokine-mediated signaling pathway, and cellular response to lipopolysaccharide. In the molecular function (MF), co-differential expressed proteins are mainly enriched in chemokine activity and CXCR chemokine receptor binding (Figure 10A). Furthermore, analysis of KEGG pathways indicates that these proteins are enriched in pathways such as the p53 signaling pathway, apoptosis Toll-like receptor signaling pathway, and TNF signaling pathway (Figure 10B).
In the context of the biological process, the co-differentially expressed proteins between LTBI and TB are predominantly enriched in functions related to tumor necrosis factor receptor binding, positive regulation of JNK cascade, tumor necrosis factor production, peptidyl-tyrosine and peptidyl-serine phosphorylation, interleukin-17 production, cell proliferation, inflammatory response, interleukin-1 beta production, T cell proliferation, and chemokine production. In the context of molecular function (MF), co-differentially expressed proteins demonstrate significant enrichment in cytokine activity, chemokine activity, and CXCR chemokine receptor binding (Figure 10C). Furthermore, KEGG pathway analysis reveals enrichment of co-differentially expressed proteins in pathways such as the p53 signaling pathway, Tuberculosis, apoptosis, necroptosis, neurodegeneration pathway, TNF signaling pathway, and NF-kappa B signaling pathway (Figure 10D).
In the biological process (BP), co-differential expressed proteins among Non-severe and severe TB are mainly enriched in protein kinase B signaling, positive regulation of interleukin-17 production, positive regulation of interleukin-12 production, cellular response to interferon-gamma, positive regulation of interleukin-1 beta production, positive regulation of chemokine production, positive regulation of NIK/NF-kappaB signaling, positive regulation of peptidyl-tyrosine and peptidyl-serine phosphorylation(figure 10E). The KEGG pathway analysis revealed that co-differentially expressed proteins exhibit enrichment in several pathways, including the PPAR signaling pathway, IL-17 signaling pathway, primary immunodeficiency, endocrine resistance, and NF-kappa B signaling pathway (Figure 10F).
The outcomes of the GO and KEGG analyses of co-differentiated proteins across different groups highlight the significant involvement of immune regulation and cell signaling in the development of tuberculosis. Specifically, the dysregulation of cytokine and chemokine production and signaling emerges as a crucial factor in distinguishing between various stages and severities of the disease. These findings contribute to the advancement of knowledge regarding the pathophysiology of tuberculosis and have the potential to facilitate the development of targeted therapeutic interventions and diagnostic strategies.