Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Mol. Biosci., 04 August 2025

Sec. Molecular Diagnostics and Therapeutics

Volume 12 - 2025 | https://doi.org/10.3389/fmolb.2025.1624131

Inflammation-nutrition biomarker model for survival prediction in lung cancer patients with concurrent tuberculosis

Hongqi Zhou&#x;Hongqi Zhou1Zihao Zhao&#x;Zihao Zhao2Jinhai WangJinhai Wang3Weiyun Jin
Weiyun Jin4*Bensong Xian
Bensong Xian5*Lindi LiLindi Li1XiangWen NieXiangWen Nie1WeiWei WuWeiWei Wu1Ran ChenRan Chen1QiZhen XieQiZhen Xie1HaiXia WuHaiXia Wu1WeiWei JiangWeiWei Jiang1Min TangMin Tang1YuXin LiYuXin Li6
  • 1Oncology Department, Guiyang Public Health Treatment Center, Guiyang, China
  • 2Department of Orthopedics, Guiyang Public Health Treatment Center, Guiyang, China
  • 3Medical Records Office, Guiyang Public Health Treatment Center, Guiyang, China
  • 4College of Humanities Education, Inner Mongolia Medical University, Hohhot, China
  • 5School of Health Management, Inner Mongolia Medical University, Hohhot, China
  • 6Neurology Department, Guiyang Public Health Treatment Center, Guiyang, China

Objectives: To explore the prognostic value of eight inflammation-nutrition biomarkers in patients with lung cancer and tuberculosis as no multidimensional prognostic models for this comorbid population are available currently.

Methodology: A retrospective study included 100 patients with lung cancer and tuberculosis admitted to a tertiary hospital from October 2019 to October 2024. Eight inflammation-nutrition markers (NLR, PLR, SII, LMR, PNI, HALP, HRR, ALB/GLB) were chosen as predictors while overall survival (OS) was the major event. Feature selection was implemented by LASSO regression; a Cox proportional hazards model was established afterwards. The nomogram’s performance was assessed by ROC curve and C-index as well as the calibration using bootstrap resampling. The statistical power was calculated by PowerSurvEpi and sensitivity analyses were implemented to test the robustness of the model.

Results: There were six predictors remaining in the final model including diabetes, ECOG PS, NLR, PNI, HRR and RDW. Among them, ECOG PS was an independent prognostic factor (HR = 1.76, p = 0.04). The nomogram achieved a good performance (C-index = 0.71), an AUC of 0.693 for 3-year OS as well as an excellent calibration (Bootstrap P > 0.05). In the high-risk subgroup with ECOG PS ≥ 2 and NLR>8, the 5-year survival rate was close to zero. The model achieved an adequate statistical power (83%, α = 0.05). Sensitivity analysis revealed an significant interaction between ECOG PS and NLR (p = 0.032) and NLR>8 was the most robust threshold for this interaction.

Conclusion: This is the first study to establish and validate a combined inflammation-nutrition prognostic model for patients with lung cancer and tuberculosis. Our model provides a quantitative tool to stratify individual risk and offers evidence for the usage of nutritional interventions in high-risk patients.

1 Introduction

Lung cancer is the leading cause of cancer death worldwide. In 2023 Global Cancer Statistics published by the International Agency for Research on Cancer (IARC), there were about 2.47 million and 1.76 million new cases of lung cancer and related deaths in 2022, and the number in China accounted for almost 40% (Bray et al., 2024). Tuberculosis (TB) is also one of the top ten causes of death. According to the data published by World Health Organization (WHO) in 2023, In China, approximately 741,000 new cases of tuberculosis and 25,000 related deaths occur annually. The country ranks third worldwide in terms of tuberculosis burden (World Health Organizatio n, 2023),The southwest region accounts for a disproportionately high share of newly reported pulmonary tuberculosis cases nationwide, particularly in Yunnan, Guizhou, and Sichuan. In some areas, the annual incidence rate reaches or exceeds 1.5 times the national average (World Health Organizatio n, 2023). Lung cancer with TB is especially common in immunocompromised population. The synergistic effect may aggravate the chronic inflammation and immune suppression and may result in worse prognosis in this population (Dheda et al., 2020). The 5-year overall survival (OS) rate of lung cancer is about 18%–22%, while the OS rate is further decreased to 12%–15% in patients with concurrent TB (Hong et al., 2022). However, systematic prognostic studies in this special population are severely lacking. This pattern of comorbidity leads to poor quality of life and puts a heavy economic burden on healthcare systems. Therefore, it also highlights the need for both basic mechanistic studies and health policy innovation. Systemic inflammation and nutritional status are increasingly important in the prognosis of malignancies. The indicators explored in this study include neutrophil-to-lymphocyte ratio (NLR), lymphocyte-to-monocyte ratio (LMR), platelet-to-lymphocyte ratio (PLR), systemic immune-inflammation index (SII), hemoglobin–albumin–lymphocyte–platelet score (HALP), prognostic nutritional index (PNI), hemoglobin-to-red blood cell distribution width ratio (HRR) and albumin/globulin ratio (ALB/GLB). Neutrophil-to-lymphocyte ratio (NLR) reflects the intensity of the overall inflammatory response. A higher NLR (≥3.0) is significantly associated with shorter overall survival (OS) in gastric (Li and Pan, 2023) and colorectal cancers (Chen et al., 2019). Prognostic nutritional index (PNI) is a simple measurement of nutritional and immune status based on the combination of serum albumin and lymphocyte count. Low PNI (≤45) is an independent prognostic factor for hepatocellular carcinoma (Pan et al., 2024) and breast cancer (Yang et al., 2014). Platelet-to-lymphocyte ratio (PLR) and Systemic immune-inflammation index (SII) (Ito et al., 2024) have also shown predictive value for treatment response and survival in various solid tumors. However, these parameters are widely used to evaluate the risk and long-term prognosis of malignancies and chronic inflammatory diseases, but their synergistic effects and possible mechanisms in patients with lung cancer and tuberculosis are still unclear. Recent studies have further emphasized the importance of body composition disorders in lung cancer patients receiving immunotherapy, particularly the impact on treatment outcomes (Trestini et al., 2024). Additionally, the adverse relationship between infections and immunotherapy outcomes has been well documented (Belluomini et al., 2021), and the prognostic value of hemoglobin/RDW ratio in cancer patients has been validated in renal cell carcinoma patients (Giudice et al., 2025).

Although the prognostic relevance of inflammatory and nutritional markers has been well studied in individual diseases, such as lung cancer or tuberculosis, a gap still exists for the systemic integration of multiple inflammatory, immune, and nutritional markers and overall survival (OS) among patients with lung cancer and tuberculosis (Mazzella et al., 2024; Shen et al., 2021; Luczynski et al., 2022). The predictive efficacy and clinical translational relevance of combined multivariate models in these patients remain to be explored. This gap limits precise risk stratification and development of individualized intervention strategies.

This is the first study to develop a multivariate prognostic model consisting of eight multidimensional inflammation-nutrition markers in patients with lung cancer and tuberculosis. We employed comprehensive multivariate analysis to investigate the predictive value and underlying mechanisms of these markers for 1, 3, and 5 years survival rates. This is an important gap gap in the current profile of prognostic assessment in this special comorbid population and has clinical support. Our findings may provide a new combination of biomarkers for prognosis in comorbid patients and establish a theoretical basis for future clinical trials of immunonutritional interventions.

1.1 Study population

This retrospective cohort study included cases admitted to a tertiary general hospital in Guiyang from 31 October 2019, to 31 October 2024. All patients were prospectively followed until 31 October 2024. A total of 100 patients with concurrent lung cancer and pulmonary tuberculosis were ultimately enrolled. We acknowledge that including patients across all stages (I-IV) introduces population heterogeneity. However, this approach reflects the real-world clinical scenario where prognostic tools are needed across different stages, particularly in resource-limited settings managing the complex diagnosis of concurrent lung cancer and tuberculosis, where traditional staging may have diminished predictive value due to confounding inflammatory responses. Inclusion criteria were as follows:1. Age between 18 and 85 years; 2. Histologically confirmed lung cancer with TNM stage I–IV; 3. Newly diagnosed pulmonary tuberculosis, as defined by WHO guidelines, with completion of a standard anti-tuberculosis regimen; 4. No previous anti-tumor therapy (including surgery, radiotherapy, chemotherapy, targeted therapy, or immunotherapy) prior to enrollment; 5. Completion of relevant hematological examinations before anti-tumor treatment; 6. Complete medical records. Exclusion criteria were:1. Concomitant malignancies in other organ systems; 2. Severe cardiac, renal, or hepatic failure; 3. Missing data for exposure or outcome variables. All data, including demographic information, laboratory results, and follow-up records, were extracted from electronic medical records. Data collection was conducted by trained researchers to ensure accuracy and consistency.

1.2 Exposure variables

The exposure variables included the neutrophil-to-lymphocyte ratio (NLR), platelet-to-lymphocyte ratio (PLR), systemic immune-inflammation index (SII), lymphocyte-to-monocyte ratio (LMR), prognostic nutritional index (PNI), hemoglobin-albumin-lymphocyte-platelet score (HALP), hemoglobin-to-red cell distribution width ratio (HRR), and albumin-to-globulin ratio (ALB/GLB). These variables were measured based on blood samples collected at hospital admission. All variables were treated as continuous variables. For analytical purposes, certain variables were categorized according to clinical relevance or recommendations from the literature. The definitions and calculation methods for these variables strictly followed international standards.

1.3 Calculation methods for exposure variables

The following formulas were used to calculate the exposure variables. NLR (Neutrophil-to-Lymphocyte Ratio) = Absolute Neutrophil Count (×109/L)÷Absolute Lymphocyte Count (×109/L). LMR (Lymphocyte-to-Monocyte Ratio) = Absolute Lymphocyte Count (×109/L)÷Absolute Monocyte Count (×109/L). PLR(Platelet-to-Lymphocyte Ratio) = Platelet Count (×109/L)÷Absolute Lymphocyte Count (×109/L). SII (Systemic Immune-Inflammation Index) = Platelet Count (×109/L) × Absolute Neutrophil Count (×109/L)÷Absolute Lymphocyte Count (×109/L). HALP (Hemoglobin-Albumin-Lymphocyte-Platelet Score) = [Hemoglobin (g/L) × Albumin (g/L) × Absolute Lymphocyte Count (×109/L)]÷Platelet Count (×109/L). PNI (Prognostic Nutritional Index) = Serum Albumin (g/L)+[5×Absolute Lymphocyte Count (×109/L)]. HRR(Hemoglobin-to-Red Cell Distribution Width Ratio) = Hemoglobin (g/L)÷Red Cell Distribution Width (%). ALB/GLB (Albumin-to-Globulin Ratio) = Serum Albumin (g/L)÷Serum Globulin (g/L). All units and calculation methods strictly adhere to international standards.

1.4 Outcome variable

The primary outcome was overall survival (OS), defined as the interval from the date of diagnosis to either the date of death or the last follow-up. Outcome status was determined based on follow-up records, including confirmation of death or the last known date alive. All survival outcomes were independently assessed by researchers blinded to study exposure, ensuring objectivity and consistency of results.

1.5 Relevant covariates

Covariates included age (≥60 or <60 years), sex (male or female), body mass index (BMI; ≥25.5 or <25.5), smoking history (yes or no), pathological type (adenocarcinoma, squamous cell carcinoma, small cell carcinoma, or others), lymph node metastasis (yes or no), distant metastasis (yes or no), clinical stage (I–IV), differentiation grade (poor or moderate/well), ECOG performance status (0–1 or ≥2), diabetes (yes or no), hypertension (yes or no), and HIV infection (yes or no). Selection of these covariates was based on prognostic factors identified in previous literature. All variables were extracted from patients’ electronic medical records.

1.6 Ethics statement

This study was approved by the Ethics Committee (Approval No: Guiyang Public Health Treatment Center 2025–07). As a retrospective cohort study, all patient data were anonymized, and the requirement for informed consent was waived. The research adhered to the Declaration of Helsinki and relevant ethical guidelines to ensure the protection of patient privacy. All data were used solely for scientific research purposes, and the study team is committed to not using patient information for any purposes unrelated to this research. The data of this study are available from the corresponding author upon reasonable request.

1.7 Statistical methods

The distribution of continuous variables was assessed using the Shapiro-Wilk test. Variables with a normal distribution are presented as mean ± standard deviation, while non-normally distributed data are reported as median (interquartile range). Categorical variables are summarized as counts and percentages. Between-group differences were evaluated using the chi-squared test for categorical variables, the independent-samples t-test for normally distributed continuous variables, and the Mann-Whitney U test for non-normally distributed continuous variables. Survival curves were generated via the Kaplan-Meier method, with differences between groups compared by the log-rank test. The Cox proportional hazards regression model was employed to assess associations between exposure variables and overall survival (OS). Three hierarchical models were constructed: Model 1 was unadjusted; Model 2 adjusted for demographic factors; Model 3 further adjusted for clinicopathological variables. Variable selection was performed via least absolute shrinkage and selection operator (LASSO) regression. The optimal regularization parameter (λ) was determined using ten-fold cross-validation, applying the one-standard-error rule to enhance model generalizability and minimize the partial likelihood deviance, thereby reducing multicollinearity and overfitting risk. A nomogram was developed based on the selected variables to predict 1, 3, and 5years survival probabilities. Model performance was evaluated by time-dependent receiver operating characteristic (ROC) curves (area under the curve, AUC) and concordance index (C-index). Calibration was assessed by calibration curves using 1,000 bootstrap resamples to compare predicted and observed outcomes. The PowerSurvEpi package was used to calculate statistical power for the Cox regression analyses. Sensitivity analyses, including outlier handling, variable stability, interaction assessment, and threshold adjustment, were conducted to ensure data robustness. All statistical tests were two-sided, and significance was set at P < 0.05. Multiple comparisons were adjusted using Bonferroni correction. All statistical analyses were performed using R software (version 4.3.0).

1.8 Research technical roadmap

The overall research design and analytical workflow of this study are illustrated in Figure 1.

Figure 1
Flowchart outlining study process: 1. Identification of exposure and outcome variables. - Identification of relevant covariates. - Identification of study population with sample size criteria.2. Data collection and statistical analysis. - Univariate, dimension reduction, and multivariate analyses.3. Development and validation of a risk prediction model. - Includes nomogram, survival curves, ROC curve, calibration curve, performance, and sensitivity analysis.4. Manuscript preparation.

Figure 1. Technical roadmap.

2 Results

2.1 Analysis of prognostic factors associated with inflammatory and nutritional markers in lung cancer patients with concurrent pulmonary tuberculosis

2.1.1 Univariate cox regression analysis

Univariate Cox regression analysis identified nine potential prognostic factors: HIV infection, diabetes, Eastern Cooperative Oncology Group performance status (ECOG PS), NLR, PNI, PLR, HRR, ALB, and RDW. The overall model demonstrated good fit (likelihood ratio test, χ2 = 21.27, p = 0.011). Details are presented in Table 1.

Table 1
www.frontiersin.org

Table 1. Clinical pathological characteristics and survival analysis of lung cancer patients.

2.1.2 Feature selection and dimensionality reduction using LASSO regression

LASSO regression analysis was performed with ten-fold cross-validation to determine the optimal log(λ) value. The most predictive variables were selected within the optimal log(λ) range. Figure 2 illustrates the coefficient profiles as a function of log(λ) (optimal log(λ) approximately −1.8 to −2.4). Figure 3 displays the corresponding regularization path, showing the C-index ranging from 0.62 to 0.63. Figure 4 depicts the stepwise variable selection process, which ultimately generated a parsimonious predictive model. This approach identified the most clinically relevant prognostic factors, reduced overfitting, and enhanced clinical applicability.

Figure 2
Plot of partial likelihood deviance versus the logarithm of lambda (λ). The deviance is represented by a red line with dots, displaying a general downward trend with confidence intervals. Vertical dashed lines mark specific λ values.

Figure 2. Cross-Validation Curve for the LASSO Regression Model. Note: The x-axis represents the logarithm of λ [log(λ)], and the y-axis represents the partial likelihood deviance. The red dotted line indicates the mean partial likelihood deviance, while the vertical bars represent the 95% confidence intervals for each λ value. The numbers above the curve show the number of non-zero coefficients retained in the model at each λ. The left dashed line marks the optimal λ with the minimum deviance (log(λ)≈−2.4). The right dashed line marks the λ corresponding to the most parsimonious model under the “one standard error rule” (log(λ)≈−1.8).

Figure 3
A line graph with dotted red points showing the C-index on the vertical axis against the logarithm of lambda on the horizontal axis. Error bars are present for each point. The C-index peaks around a log(lambda) of negative five, with values ranging from 0.55 to 0.65. The top label shows corresponding numbers from eight to two. Two vertical dashed lines highlight particular log(lambda) values.

Figure 3. Cross-Validation Curve and Selection of the Optimal Regularization Parameter in the LASSO Regression Model. Note: This figure illustrates the relationship between the strength of LASSO regularization (log(λ), x-axis) and model discriminatory ability, as measured by the concordance index (C-index, y-axis). Red dots indicate the mean C-index values derived from cross-validation; vertical bars represent the 95% confidence intervals for each λ value. The numbers at the top indicate the number of non-zero coefficients (selected variables) retained in the model at each λ. The left vertical dashed line (log(λ) ≈ −5) identifies the λ value that maximizes the C-index. The right vertical dashed line (log(λ)≈−3) corresponds to the most parsimonious model within one standard error of the maximum C-index. The model maintains stable discriminatory performance (C-index≈0.62–0.63) over a broad range of regularization strengths before a marked decline in performance is observed at higher regularization levels (log(λ)>−3).

Figure 4
Line plot of coefficient paths for a LASSO regression model, with log lambda on the x-axis and coefficients on the y-axis. Nine curves each represent a variable's coefficient path as lambda varies. A vertical dashed line at around -3 indicates an optimal lambda value.

Figure 4. LASSO Coefficient Paths for Variables in the Predictive Model. Note: This figure illustrates the relationship between log(λ) and the standardized regression coefficients. The lower x-axis corresponds to log(λ); the upper x-axis indicates the number of non-zero coefficients (variables) retained in the model at each log(λ) value. Each colored line represents a distinct predictor variable. The vertical dashed line marks the optimal λ value determined by cross-validation. Black and green lines denote predictors with the strongest positive associations, while blue lines represent variables negatively associated with the outcome. As log(λ) increases, model complexity decreases, and many coefficients shrink toward zero or become exactly zero. The model thereby identifies the most relevant features, enabling dimensionality reduction and the selection of key predictors.

2.1.3 Multivariate cox regression analysis

This study used all-cause mortality in patients with lung cancer and concurrent pulmonary tuberculosis as the dependent variable. LASSO regression was applied to identify the six most relevant predictive variables: diabetes, ECOG PS, NLR, PNI, HRR and RDW. A multivariate Cox regression model was then constructed to assess the risk of adverse outcomes associated with inflammatory and nutritional markers in this patient population. The results indicated that an ECOG PS ≥ 2 was significantly associated with poor prognosis (95% CI:1.02–3.02, p = 0.04). Compared to patients with ECOG PS 0–1, those with ECOG PS ≥ 2 had a 76% higher risk of mortality. The model demonstrated good discriminatory ability (log-rank test, p = 0.00688) and moderate predictive performance (C-index = 0.65). See Table 2 and Figure 5.

Table 2
www.frontiersin.org

Table 2. Multivariate cox regression analysis of prognostic factors.

Figure 5
Forest plot illustrating hazard ratios for multiple variables: diabetes, ECOG, NLR, PNI, HRR, and RDW, with confidence intervals. ECOG level two shows a significant p-value of 0.0414. Hazard ratios are represented by squares, and lines indicate confidence intervals. The reference line is at one. Other p-values are noted for each variable.

Figure 5. Forest Plot of Hazard Ratios for Prognostic Indicators. Note: This figure presents the hazard ratios (HRs) and their 95% confidence intervals (CIs) estimated by the risk model. The vertical dashed line indicates the reference value of HR = 1.0. Squares represent the point estimates of HRs, with their size reflecting statistical precision; horizontal lines denote the corresponding 95% CIs. HRs to the right of the reference line (>1.0) indicate increased risk, while those to the left (<1.0) indicate reduced risk. The corresponding p-values are shown on the right; an asterisk (*) indicates statistical significance (p < 0.05).

3 Development of the risk prediction model

3.1 Nomogram for risk prediction

A nomogram was developed based on six optimal variables—PNI, HRR, diabetes, NLR, RDW and ECOG PS—to predict 1, 3, and 5 years survival rates in patients with lung cancer and concurrent pulmonary tuberculosis. The predictors retained in the final model (diabetes, RDW, NLR, PNI, and HRR), although they did not achieve statistical significance (p > 0.05), were retained based on their contribution to overall model performance and their potential value in clinical practice for inflammation-nutrition assessment. Within the nomogram, PNI values of 38–50 and HRR values of 2.4–3.6 were associated with the most favorable prognosis. An NLR greater than 8 indicated a poorer prognosis. RDW values of 30–65 demonstrated a complex, nonlinear relationship with outcomes. Diabetes was significantly associated with prognosis. ECOG PS emerged as the strongest independent predictor of adverse outcomes. The nomogram demonstrated good discriminatory ability (C-index = 0.71; 95% CI:0.66–0.76). See Figure 6.

Figure 6
A nomogram chart with multiple scales and density plots for points related to various medical parameters including PNI, HRR, diabetes, NLR, RDW, and ECOG. Total points range from 0 to 400, with probabilities for survival times exceeding one, three, and five years shown below. Each parameter has a corresponding density distribution depicting the score contribution.

Figure 6. Nomogram for Predicting 1, 3, and 5 Year Survival Probability. Note: The upper x-axis displays the total point score; the lower x-axis indicates survival probability. The solid line represents the predicted probabilities from the fitted model. Each variable contributes a specific number of points (top, Points). The total points are used to estimate the 1, 3, and 5 years survival probabilities: Pr (time>1 year), Pr (time>3 years), and Pr (time>5 years).

3.2 Kaplan-meier survival analysis by ECOG PS-based risk groups

Based on the strongest predictor, ECOG PS, Kaplan-Meier survival curves were used to compare OS between the high-risk group (ECOG PS ≥ 2) and the low-risk group (ECOG PS 0–1). The results demonstrated that patients in the high-risk group (ECOG PS ≥ 2) had significantly poorer outcomes (log-rank test, p = 0.0025), especially after 2 years of follow-up. By the end of the 5-year observation period, the survival probability in the high-risk group approached zero, while a small proportion of patients in the low-risk group remained alive. See Figure 7.

Figure 7
Kaplan-Meier survival curve compares high and low risk groups over six years. The red line shows high risk, and the green line shows low risk. The low risk group has higher survival probability. The p-value is 0.0025, indicating statistical significance. The bottom table shows the number at risk over time, with both groups starting at fifty individuals.

Figure 7. Kaplan-Meier Survival Curves by ECOG PS-Based Risk Groups. Note: The x-axis represents time (years), and the y-axis shows survival probability. The red line indicates the high-risk group, and the green line indicates the low-risk group. Shaded areas represent the 95% confidence intervals for the corresponding survival probabilities.

4 Model performance evaluation and validation

4.1 ROC curve analysis

Receiver operating characteristic (ROC) curves were used to evaluate the predictive performance of the model at 1, 3, and 5 years. The model demonstrated consistent, moderate predictive ability across all follow-up intervals. The area under the curve (AUC) was 0.656 at 1 year, 0.693 at 3 years, and 0.689 at 5 years, with the best performance observed at the 3-year follow-up (AUC = 0.693). See Figure 8.

Figure 8
Receiver Operating Characteristic (ROC) curve displaying three plots: purple for one year with an AUC of 0.656, orange for three years with an AUC of 0.693, and green for five years with an AUC of 0.689. X-axis is 1−specificity, and Y-axis is sensitivity. A diagonal line represents randomness.

Figure 8. Time-Dependent ROC Curves for the Prediction Model at 1, 3, and 5 Years. Note: The x-axis represents 1 minus specificity (false-positive rate), and the y-axis represents sensitivity (true-positive rate).

4.2 Calibration curve analysis

Calibration curves were validated using 1,000 bootstrap resamples to reduce overfitting bias. The predicted probabilities showed excellent agreement with the observed outcomes, and the model demonstrated good calibration at all three time points (bootstrap P > 0.05). See Figure 9.

Figure 9
Calibration plot comparing nomogram-predicted overall survival (OS) percentages with observed OS over one, three, and five years. The plot includes red, green, and blue lines corresponding to each time frame, with error bars indicating variability.

Figure 9. Calibration Curves of the Nomogram for Predicting 1, 3, and 5 Year Overall Survival Rates. Note: The calibration curves compare the overall survival probabilities predicted by the nomogram (x-axis) with the observed survival rates at 1 year (red line), 3 years (green line), and 5 years (blue line) (y-axis). Data points represent grouped observations, with vertical error bars indicating the 95% confidence intervals. The diagonal grey line represents perfect calibration. Most confidence intervals intersect the diagonal, indicating strong agreement between predicted and observed outcomes.

4.3 Power analysis

The statistical power of this study, assessed using the PowerSurvEpi package (significance level α = 0.05), was 83%. This result demonstrates that the study design is adequate to support the predictive value of ECOG PS. See Table 3.

Table 3
www.frontiersin.org

Table 3. Power analysis results.

4.4 Sensitivity analysis

Sensitivity analyses revealed several key findings. First, After excluding extreme outliers (NLR >20, n = 1) and HIV-positive cases (n = 1), the model was reconstructed, the C-index decreased slightly from 0.71 to 0.70, with no statistically significant difference, indicating stable model performance. Second, a significant synergistic effect was observed for patients with both ECOG PS ≥ 2 and NLR> 8 (95% CI: 1.12–4.08, p = 0.032), underscoring the importance of functional assessment in highly inflamed, comorbid populations. Third, the NLR> 8 threshold demonstrated the most robust performance in the sensitivity analyses. See Table 4.

Table 4
www.frontiersin.org

Table 4. Sensitivity analysis.

5 Summary

5.1 Main findings

This retrospective cohort study investigated the association between eight inflammation-nutrition biomarkers (NLR, PLR, SII, LMR, PNI, HALP, HRR and ALB/GLB) and OS in patients with lung cancer complicated by tuberculosis. A multi-stage modeling approach yielded the following results: 1. LASSO regression identified six optimal variables PNI, HRR, diabetes, NLR, RDW, and ECOG PS for constructing a nomogram model. This model demonstrated robust predictive performance (C-index = 0.71; 95% CI:0.66–0.76). 2. Cox regression revealed that PNI values between 38 and 50 and HRR values between 2.4 and 3.6 were associated with the best prognosis. An NLR>8 indicated worse outcomes, and RDW (30–65) showed a non-linear relationship with prognosis. Diabetes was correlated with survival, and ECOG PS emerged as an independent prognostic risk factor. 3. Kaplan-Meier survival analysis stratified by ECOG PS indicated that patients in the high-risk group (ECOG PS ≥ 2) had significantly poorer outcomes (log-rank p = 0.0025), with a 76% higher risk of death (95% CI:1.02–3.02, p = 0.041). In the high-risk subgroup (ECOG PS ≥ 2 and NLR>8), the 5-year survival rate was nearly zero. 4. The area under the ROC curve (AUC) for 1, 3, and 5 years survival predictions was 0.656, 0.693, and 0.689, respectively, with the model performing best at 3 years (AUC = 0.693). 5. Calibration curves validated with 1,000 bootstrap resamples showed good agreement between predicted and observed outcomes at all three key time points (bootstrap P > 0.05). 6. Power analysis using the PowerSurvEpi package indicated a statistical power of 83% (significance level = 0.05), supporting the adequacy of the study design and the predictive value of ECOG PS. 7. A series of sensitivity analyses were conducted, including outlier handling, variable stability assessment, evaluation of interactions, and determination of optimal cutoff values. In particular, the HIV-positive case (n = 1) was excluded due to insufficient statistical power. These analyses confirmed the robust performance of the model. Notably, the synergistic interaction between ECOG PS and NLR was significant (95% CI:1.12–4.08, p = 0.032), and NLR>8 proved to be the most stable threshold, consistent with previous reports (Cupp et al., 2020). These findings highlight the need for careful assessment of physical function in highly inflamed comorbid populations. This study is the first to systematically integrate multiple inflammatory, immune, and nutritional biomarkers for prognostic evaluation in patients with both lung cancer and tuberculosis, providing a novel biomarker panel for individualized prognosis in this unique patient population. Our study revealed a synergistic effect between ECOG PS and the NLR. Poor ECOG PS suggests worsening malnutrition, chronic inflammation, and immune dysfunction, which may elevate NLR levels. Elevated NLR further reflects tumor-associated inflammation and immune evasion. This synergy indicates diminished host antitumor capacity and a pro-inflammatory environment, both of which are associated with poorer prognosis.

5.2 Comparison with existing research

This study demonstrated significant consistency with the findings of Chen et al. in metastatic colorectal cancer (Chen et al., 2019). Specifically, each one-unit increase in NLR was associated with a 7% increase in mortality risk in our cohort, compared to 5% in their study. The threshold for adverse prognosis in our study was NLR>8, higher than the NLR>3 used in Chen et al., suggesting that the prognostic value of inflammatory imbalance may be applicable across different cancer types. Notably, in the context of tuberculosis comorbidity, the discriminatory power of PNI (95% CI:0.91–0.99) in our cohort was markedly reduced compared to the hepatocellular carcinoma cohort reported by Pan et al. (2024) (95% CI: 0.76–0.89). This difference may be attributed to immune-metabolic disturbances unique to Mycobacterium tuberculosis infection. Mechanistic studies have shown that mycobacterial membrane proteins activate the NF-κB signaling pathway via TLR4, resulting in a 3.5-fold increase in IL-6 transcription. Persistent IL-6 stimulation downregulates albumin mRNA expression in hepatocytes and increases Fas receptor expression in lymphocytes, leading to higher rates of apoptosis. These changes ultimately manifest as increased NLR and decreased PNI (O'Garra et al., 2013), reducing the protective impact of nutrition-immune indices. The association between diabetes and poor prognosis observed in our study is consistent with the findings of Donath MY. Hyperglycemia (HbA1c>7%) activates the NLRP3 inflammasome, leading to increased IL-6 secretion from monocytes, reduced efficiency of serum albumin synthesis, and a positive correlation between the homeostasis model assessment of insulin resistance (HOMA-IR) and NLR. This forms a vicious cycle of hyperglycemia, inflammation, and malnutrition (Donath and Shoelson, 2011). Methodologically, although Mazzella et al. (2024) used the LASSO-Cox model in their study of non-small cell lung cancer, the absence of adjustment for tuberculosis infection and the selective inclusion of only NLR and PLR as inflammatory markers resulted in a lower model discrimination (C-index = 0.68) compared to our study (C-index = 0.71). While Yang et al. (2019) confirmed the prognostic value of inflammation markers in a multicenter study, they did not account for the immunosuppressive effects mediated by CD163+M2 macrophages via the PD-L1/IL-10 axis in the tuberculosis microenvironment (Quail and Joyce, 2013), potentially biasing PLR-specific assessments. Importantly, our study is the first to construct a multidimensional predictive model using LASSO-selected variables in patients with comorbid lung cancer and tuberculosis. The high-risk characteristic of ECOG PS ≥ 2 (HR = 1.76) identified in our cohort was notably higher than that reported in the immunotherapy cohort by Ito et al. (2024) (HR = 1.32). This suggests that tuberculosis-associated systemic inflammation may exacerbate organ dysfunction through a TNF-α/IL-1β positive feedback loop, providing new insights for precise interventions in comorbid populations. These differences and mechanistic explanations enhance our understanding of the prognostic roles of inflammation and nutrition in patients with both lung cancer and tuberculosis.

5.3 Clinical implications

Based on our findings, this is the first study to construct an inflammation-nutrition prognostic multi-parameter model (INMPM) in lung cancer tuberculosis comorbidity. Through the systematic integration of eight inflammation; nutrition biomarkers and optimal variable selection via LASSO regression, the C-index of our model was significantly improved. Methodologically, our study fills an important gap in the field of prognostic evaluation for lung cancer tuberculosis comorbidity and provides practical evidence support for the implementation of targeted nutrition intervention. Our study demonstrated the following main findings: 1. ECOG PS is an independent prognostic risk factor in lung cancer patients with tuberculosis comorbidity. 2. In clinical practice, targeted evaluation of physical function is necessary for all comorbid patients with high inflammation; it is particularly suggested that clinicians should prioritize nutritional intervention, and if necessary, combined anti-infective therapy for high-risk subgroups of comorbidity patients with ECOG PS ≥ 2 and NLR>8. 3. The model exhibited high statistical power (83%) and strong robustness; the model’s prediction accuracy for 3-year OS in comorbidity patients was excellent. These findings could provide clinicians with a quantitative reference for therapeutic decision-making in lung cancer tuberculosis comorbidity, especially for evaluation of candidacy for aggressive anti-tumor therapy and timing of nutrition intervention.

5.4 Strengths of the study

This study exhibited the following methodological innovations and advantages:First, it was the first to clarify the synergistic mechanism between inflammation-nutrition imbalance and lung cancer prognosis in tuberculosis comorbidity. Second, the data analysis avoided the limitations of conventional conventional single-stage modeling and adopted a multi-stage approach. In this approach, LASSO regression was used for variable selection and multivariate Cox regression was used to construct the model. This strategy effectively alleviated the influence of heterogeneous comorbidity and substantially improved the performance of the model (C-index = 0.71 vs 0.63–0.65 for traditional models). The predictive performance was further validated using the calibration curves based on 1,000 bootstrap resamples. The result showed that the prediction error rate was less than 5% at 1, 3, and 5 years time points. Statistical power analysis based on PowerSurvEpi function exhibited that the Cox regression had 83% power (significance level = 0.05) and the sample design was appropriate. Sensitivity analyses demonstrated the robustness of the model. These methods provided clinicians with a highly reliable and quantitative decision-making tool. Finally, this study provided precise nutritional intervention criteria for high-risk subgroups (ECOG PS ≥ 2 and NLR>8) and suggested that combined anti-infective therapy should be adopted when necessary, which facilitated individualized treatment.

5.5 Limitations

This study has several limitations. First, due to the unique characteristics of the study population and resource limitations, this study was conducted as a single-center retrospective analysis with a relatively small sample size and without external validation using an independent cohort. Although a power analysis confirmed the adequacy of the sample size, validation of interaction effects requires larger cohorts, and the retrospective design is prone to limitations in causal inference. Future multicenter studies are needed to enhance the generalizability of our conclusions. Second, the study population was primarily from Guizhou Province, China. Their unique genetic background and tuberculosis exposure characteristics may influence the prognostic thresholds of inflammation-nutrition biomarkers. Thus, the generalizability of the model should be validated in cohorts from different ethnic and geographic backgrounds. Finally, due to the limitations of retrospective data, several key variables [such as adherence to anti-tuberculosis therapy, PD-L1 status, and LIPI score (De Giglio et al., 2024)] could not be comprehensively calculated and were not included in the analysis owing to missing data on specific treatment regimens, therapeutic responses, and baseline LDH values. Future studies should expand sample sizes, validate in stage-specific cohorts, and incorporate additional clinical interventional factors and biomarkers to verify the generalizability of these interactions.

6 Conclusion

This study developed the first multidimensional inflammation-nutrition prognostic model for patients with concurrent lung cancer and tuberculosis (C-index = 0.71). It also identified, for the first time, the synergistic high-risk effect of ECOG PS and NLR in this comorbid population. These findings provide a quantifiable tool to support clinical decision-making.

Based on the stratified treatment pathway established by the model, this study recommends prioritizing nutritional and immunomodulatory interventions for high-risk patients (ECOG PS ≥ 2 and NLR>8), with adjunctive anti-infective therapy when necessary. This approach can help healthcare institutions optimize resource allocation, ultimately aiming to improve the quality of life and survival of comorbid patients.

Future research should focus on three main areas: 1. This study serves as a foundation for exploratory research. In the future, multicenter clinical trials will be conducted with larger sample sizes. Additional variables, such as adherence to anti-tuberculosis therapy and PD-L1 expression, will also be included. These efforts aim to further evaluate the comprehensiveness and stability of the model. 2. To investigate the molecular pathways underlying the tumor inflammatory microenvironment and immune escape mechanisms in patients with lung cancer and tuberculosis. 3. Conducting clinical trials evaluating the efficacy of immunomodulators combined with nutritional interventions. These efforts are expected to reshape treatment strategies and advance precision medicine in the management of lung cancer with tuberculosis.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.

Author contributions

HZ: Methodology, Data curation, Writing – review and editing, Writing – original draft. ZZ: Data curation, Methodology, Writing – review and editing. JW: Supervision, Writing – review and editing, Resources, Project administration. WJn: Project administration, Funding acquisition, Data curation, Writing – review and editing, Supervision. BX: Supervision, Funding acquisition, Writing – review and editing, Project administration. LL: Data curation, Supervision, Writing – review and editing. XN: Supervision, Writing – review and editing, Data curation. WW: Writing – review and editing, Data curation, Supervision. RC: Data curation, Writing – review and editing, Supervision. QX: Writing – review and editing, Supervision, Data curation. HW: Data curation, Writing – review and editing, Supervision. WJa: Writing – review and editing, Data curation, Supervision. MT: Data curation, Writing – review and editing, Supervision. YL: Data curation, Supervision, Writing – review and editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Natural Science Foundation of Inner Mongolia Autonomous Region (Grant Nos. 2022LHQN07001 and 2024QN07005).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

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.

References

Belluomini, L., Caldart, A., Avancini, A., Dodi, A., Trestini, I., Kadrija, D., et al. (2021). Infections and immunotherapy in lung cancer: a bad relationship? Int. J. Mol. Sci. 22 (1), 42. doi:10.3390/ijms22010042

PubMed Abstract | CrossRef Full Text | Google Scholar

Bray, F., Laversanne, M., Sung, H., Ferlay, J., Siegel, R. L., Soerjomataram, I., et al. (2024). Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 74 (3), 229–263. doi:10.3322/caac.21834

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, L., Hao, Y., Cong, X., Zou, M., Li, S., Zhu, Z., et al. (2019). Peripheral venous blood platelet-to-lymphocyte ratio (PLR) for predicting the survival of patients with gastric cancer treated with SOX or XELOX regimen neoadjuvant chemotherapy. Technol. Cancer Res. Treat. 18, 1533033819829485. doi:10.1177/1533033819829485

PubMed Abstract | CrossRef Full Text | Google Scholar

Cupp, M. A., Cariolou, M., Tzoulaki, I., Aune, D., Evangelou, E., and Berlanga-Taylor, A. J. (2020). Neutrophil to lymphocyte ratio and cancer prognosis: an umbrella review of systematic reviews and meta-analyses of observational studies. BMC Med. 18 (1), 360. doi:10.1186/s12916-020-01817-1

PubMed Abstract | CrossRef Full Text | Google Scholar

De Giglio, A., Leonetti, A., Comito, F., Filippini, D. M., Mollica, V., Rihawi, K., et al. (2024). Development and validation of a new tool to estimate early mortality in patients with advanced cancer treated with immunotherapy. Cancer Immunol. Immunother. 73 (12), 246. doi:10.1007/s00262-024-03836-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Dheda, K., Perumal, T., Moultrie, H., Perumal, R., Esmail, A., Scott, A. J., et al. (2022). The intersecting pandemics of tuberculosis and COVID-19: population-level and patient-level impact, clinical presentation, and corrective interventions. Lancet Respir. Med. 10 (6), 603–622. doi:10.1016/S2213-2600(22)00092-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Donath, M. Y., and Shoelson, S. E. (2011). Type 2 diabetes as an inflammatory disease. Nat. Rev. Immunol. 11 (2), 98–107. doi:10.1038/nri2925

PubMed Abstract | CrossRef Full Text | Google Scholar

Giudice, G. C., Rebuzzi, S. E., Mazzaschi, G., Pecci, F., Maffezzoli, M., Acunzo, A., et al. (2025). The prognostic value of the haemoglobin/red cell distribution width ratio in a cohort of pre-treated patients with renal cell carcinoma receiving nivolumab. Cancer Treat. Res. Commun. 43, 100927. doi:10.1016/j.ctarc.2025.100927

PubMed Abstract | CrossRef Full Text | Google Scholar

Hong, C. L., Su, C. C., and Ou, C. Y. (2022). Impact of active pulmonary tuberculosis on the prognosis of patients with upper aerodigestive cancers: an 8-year observational study in a nationwide cohort. J. Infect. Public Health 15 (12), 1540–1545. doi:10.1016/j.jiph.2022.11.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Ito, K., Hashimoto, K., Kaira, K., Yamaguchi, O., Mouri, A., Shiono, A., et al. (2024). Clinical impact of inflammatory and nutrition index based on metabolic tumor activity in non‑small cell lung cancer treated with immunotherapy. Oncol. Lett. 27, 110. doi:10.3892/ol.2024.14243

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, L. L., and Pan, L. S. (2023). Prognostic value of neutrophil-to-lymphocyte ratio in gastric cancer patientstreated with immune checkpoint inhibitors: a meta-analysis. Kaohsiung J. Med. Sci. 39 (8), 842–852. doi:10.1002/kjm2.12694

PubMed Abstract | CrossRef Full Text | Google Scholar

Luczynski, P., Poulin, P., Romanowski, K., and Johnston, J. C. (2022). Tuberculosis and risk of cancer: a systematic review and meta-analysis. PLoS One 17 (12), e0278661. doi:10.1371/journal.pone.0278661

PubMed Abstract | CrossRef Full Text | Google Scholar

Mazzella, A., Orlandi, R., Maiorca, S., Uslenghi, C., Chiari, M., Bertolaccini, L., et al. (2024). How general and inflammatory status impacts on the prognosis of patients affected by lung cancer: state of the art. Biomedicines 12 (7), 1554. doi:10.3390/biomedicines12071554

PubMed Abstract | CrossRef Full Text | Google Scholar

O'Garra, A., Redford, P. S., McNab, F. W., Bloom, C. I., Wilkinson, R. J., and Berry, M. P. (2013). The immune response in tuberculosis. Annu. Rev. Immunol. 31, 475–527. doi:10.1146/annurev-immunol-032712-095939

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, X., Hu, E., Zhou, Y., Li, L., Huang, X., and Cai, Z. (2024). The prognostic nutritional index as a predictor of efficacy and early recurrence for adjuvant transarterial chemoembolization in hepatocellular carcinoma. Clin. Res. Hepatol. Gastroenterol. 48 (6), 102344. doi:10.1016/j.clinre.2024.102344

PubMed Abstract | CrossRef Full Text | Google Scholar

Quail, D. F., and Joyce, J. A. (2013). Microenvironmental regulation of tumor progression and metastasis. Med. 2013 19 (11), 1423–1437. doi:10.1038/nm.3394

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, Y. J., Qian, L. Q., Ding, Z. P., Luo, Q. Q., Zhao, H., Xia, W. Y., et al. (2021). Prognostic value of inflammatory biomarkers in patients with stage I lung adenocarcinoma treated with surgical dissection. Front. Oncol. 11, 711206. doi:10.3389/fonc.2021.711206

PubMed Abstract | CrossRef Full Text | Google Scholar

Trestini, I., Belluomini, L., Dodi, A., Sposito, M., Caldart, A., Kadrija, D., et al. (2024). Body composition derangements in lung cancer patients treated with first-line pembrolizumab: a multicentre observational study. J. Cachexia Sarcopenia Muscle 15 (6), 2349–2360. doi:10.1002/jcsm.13568

PubMed Abstract | CrossRef Full Text | Google Scholar

World Health Organization (2023). WHO global tuberculosis report 2023. Geneva: World Health Organization Available online at: https://www.who.int/publications/i/item/9789240083851.

Google Scholar

Yang, J. R., Xu, J. Y., Chen, G. C., Yu, N., Yang, J., Zeng, D. X., et al. (2019). Post-diagnostic C-reactive protein and albumin predict survival in Chinese patients with non-small cell lung cancer: a prospective cohort study. Sci. Rep. 9 (1), 8143. doi:10.1038/s41598-019-44653-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Z., Zhang, B., Hou, L., Xie, Y., and Cao, X. (2014). Pre-operative prognostic nutritional index predicts the outcomes for triple-negative breast cancer. Tumour Biol. 35 (12), 12165–12171. doi:10.1007/s13277-014-2524-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: lung cancer, pulmonary tuberculosis, inflammation-nutrition markers, prognostic model, survival prediction

Citation: Zhou H, Zhao Z, Wang J, Jin W, Xian B, Li L, Nie X, Wu W, Chen R, Xie Q, Wu H, Jiang W, Tang M and Li Y (2025) Inflammation-nutrition biomarker model for survival prediction in lung cancer patients with concurrent tuberculosis. Front. Mol. Biosci. 12:1624131. doi: 10.3389/fmolb.2025.1624131

Received: 28 May 2025; Accepted: 15 July 2025;
Published: 04 August 2025.

Edited by:

Lorenzo Belluomini, University of Verona, Italy

Reviewed by:

Alessandra Dodi, University Hospital of Parma, Italy
Marco Sposito, University of Verona, Italy

Copyright © 2025 Zhou, Zhao, Wang, Jin, Xian, Li, Nie, Wu, Chen, Xie, Wu, Jiang, Tang and Li. 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: Weiyun Jin, dGFqd3lAMTYzLmNvbQ==; Bensong Xian, eGlhbmJlbnNvbmdAMTYzLmNvbQ==

These authors have contributed equally to this work and share first authorship

Disclaimer: 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.