Suggestions
Guide for authors
Searcher
Journal Information
Cite
Cite
Share
Download PDF
More article options
Visits
1242
Original Article
Full text access
Available online 14 July 2026

Soluble Programmed Death-1 as a Biomarker of Severity and Progression in Nodular-bronchiectatic Nontuberculous Mycobacterial Pulmonary Disease: A Multicountry Validation Cohort Study

Visits
1242
Sheng-Wei Pana,b,c,1, Koji Furuuchid,1, Yung-Che Chene, Yao-Wen Kuof,g,h, Jia-Yih Fenga,b,i, Chin-Chung Shuf,j,
Corresponding author
ccshu@ntu.edu.tw

Corresponding author.
, Kozo Morimotod,k,l,2, Jae-Joon Yimm,n,2
a Department of Chest Medicine, Taipei Veterans General Hospital, Taipei, Taiwan
b School of Medicine, National Yang Ming Chiao Tung University, Taipei, Taiwan
c Institute of Public Health, National Yang Ming Chiao Tung University, Taipei, Taiwan
d Respiratory Disease Center, Fukujuji Hospital, Japan Anti-Tuberculosis Association, Tokyo, Japan
e Department of Medicine, Kaohsiung Chang Gung Memorial Hospital, Taiwan
f Department of Internal Medicine, National Taiwan University Hospital, Taipei, Taiwan
g Department of Integrated Diagnostics & Therapeutics, National Taiwan University Hospital, Taipei, Taiwan
h College of Medicine, National Taiwan University, Taipei, Taiwan
i Institute of Emergency and Critical Care Medicine, National Yang Ming Chiao Tung University, Taipei, Taiwan
j Graduate Institute of Clinical Medicine, National Taiwan University College of Medicine, Taipei, Taiwan
k Division of Clinical Research, Fukujuji Hospital, Japan Anti-Tuberculosis Association, Tokyo, Japan
l Department of Clinical Mycobacteriosis, Graduate School of Biomedical Sciences, Nagasaki University, Nagasaki, Japan
m Department of Internal Medicine, Seoul National University College of Medicine, Seoul, Republic of Korea
n Division of Pulmonary and Critical Care Medicine, Seoul National University Hospital, Seoul, Republic of Korea
Ver más
This item has received
Article information
Abstract
Full Text
Bibliography
Download PDF
Statistics
Figures (5)
fig0005
fig0010
fig0015
fig0020
fig0025
Tables (4)
Table 1. Characteristics of healthy controls and patients with NB-pattern NTM-PD stratified by cavitary lesion.
Tables
Table 2. Logistic analysis of factors associated with treatment initiation within 1 year in patients with NB-pattern NTM-PD (n=400).
Tables
Table 3. Sensitivity analyses assessing the association between sPD-1 level and treatment initiation within 1 year.
Tables
Table 4. Cox regression analysis of factors associated with treatment initiation within 3 years in patients with NB-pattern NTM-PD (n=389).
Tables
Abstract
Background

Nontuberculous mycobacterial pulmonary disease (NTM-PD) is increasingly prevalent, yet no validated biomarker predicts progression in nodular-bronchiectatic (NB) NTM-PD. A Taiwan study linked soluble programmed death-1 (sPD-1), an immune exhaustion marker, to cavitary lesions and progression in NB-pattern NTM-PD, highlighting the need for broader validation.

Methods

In this multicountry cohort study, patients with NB-pattern NTM-PD and healthy controls were enrolled in East Asia. Blood sPD-1 levels were measured using enzyme-linked immunosorbent assay, and their associations with cavitary NB pattern and disease progression requiring anti-NTM treatment were analyzed using Cox regression.

Results

Overall, a total of 68 healthy controls and 400 patients with NB-pattern NTM-PD were included. sPD-1 levels were lower in the NTM-PD group than in the control group: 44.2pg/mL (IQR, 26.0–65.4pg/mL) vs 61.5pg/mL (IQR, 38.0–81.4pg/mL); P<.001. Cavitary NB patterns were identified in 59 patients (14.8%), who had lower sPD-1 levels than those with noncavitary NB pattern: 31.8pg/mL (IQR, 19.6–54.9pg/mL) vs 46.9pg/mL (IQR, 27.4–66.7pg/mL); P=.007. After excluding 11 patients who were receiving treatment at baseline, antibiotics were initiated in 167 patients (42.9%) during a median follow-up of 1.3 years. In multivariable analysis, risk factors for disease progression within 3 years were lower body mass index, per 1-kg/m2 increase: HR, 0.901; 95%CI, 0.842–0.963; P=.002; cavitary NB pattern: HR, 3.262; 95%CI, 2.166–4.914; P<.001; low sPD-1 level, <30.5pg/mL: HR, 1.957; 95%CI, 1.241–3.087; P=.004; and intermediate sPD-1 level, 30.5–55.5pg/mL: HR, 1.583; 95%CI, 1.007–2.488; P=.047, compared with high sPD-1 level, >55.5pg/mL.

Conclusions

This multicountry validation cohort study shows that sPD-1 is associated with cavitary NB-pattern NTM-PD at baseline and with treatment initiation during follow-up, supporting its potential role in risk stratification.

Keywords:
Cavitary nodular-bronchiectatic pattern
Disease progression
Nontuberculous mycobacterial pulmonary disease
Soluble programmed death-1 protein
Graphical abstract
Full Text
Introduction

Nontuberculous mycobacteria (NTM) can cause inflammatory pulmonary disease, most commonly presenting as nodular-bronchiectatic (NB) disease or, less frequently, fibrocavitary (FC) disease in individuals with preexisting structural lung abnormalities or immunocompromised status [1,2]. In most regions, Mycobacterium avium complex (MAC) is the predominant pathogen in nontuberculous mycobacterial pulmonary disease (NTM-PD), followed by Mycobacterium abscessus species complex (MABS) and Mycobacterium kansasii[3,4].

Over the past two decades, the global incidence of NTM-PD has increased. Incidence rates rose in Denmark from 0.4 to 1.3 per 100000 population between 2000 and 2017 and in the United States from 3.13 to 4.73 per 100000 person-years between 2008 and 2015 [5,6]. Notably, reported rates have been higher in East Asia. In Japan, the incidence of NTM-PD increased from 15.8 to 19.2 per 100000 person-years between 2013 and 2017, representing a 17.7% increase [7]. In South Korea, prevalence rose from 11.4 to 56.7 per 100000 population between 2010 and 2021, a nearly 5-fold increase [8]. Mortality also remains substantial, with a 15-year cumulative mortality of 36.4% in South Korea and an 8-year all-cause mortality of 45.2% in Taiwan [9,10].

Increased mortality in patients with NTM-PD has been associated with cavitary lesions, including cavitary lesions within the NB pattern, hereafter referred to as cavitary NB pattern, which is associated with an approximately 2-fold higher mortality risk [9,11]. Accordingly, current guidelines recommend early antibiotic therapy not only for FC disease but also for patients with cavitary NB pattern at baseline [2,12,13].

Compared with patients with fibrocavitary (FC) NTM-PD, patients with nodular-bronchiectatic (NB) disease usually have a more indolent course and lower bacterial burden, but they are prone to persistent culture positivity without treatment [11,14]. Among patients with NB-pattern NTM-PD, cavitary NB disease, reported in 10.8–22.3% of cases, is a significant risk factor for radiographic progression and clinical deterioration requiring treatment [11,15]. Given that approximately 70–90% of patients with NTM-PD in Taiwan, South Korea, and Japan exhibit the NB pattern – characterized by bronchiectasis, bronchiolitis, and nodules – our understanding of the clinical course and predictors of progression beyond cavitation in this phenotype remains limited [3,14,16,17].

Host immune responses to pathogen interactions in NB-pattern NTM-PD may be crucial determinants of disease course. As a biomarker of immune exhaustion, elevated levels of circulating CD4+ lymphocytes expressing programmed death-1 (PD-1) have been reported to correlate with severe disease and radiographic progression in patients with NTM-PD [18]. Moreover, a recent Taiwan cohort study of patients with NB-pattern NTM-PD found that blood levels of soluble PD-1 (sPD-1) were negatively correlated with cavitary NB lesions and disease progression [15]. Unlike PD-1 expressed on T cells, sPD-1 is generated through alternative splicing that eliminates the transmembrane domain and can be measured directly in blood [19]. sPD-1 is believed to counteract immune exhaustion by inhibiting the interaction between membrane-bound PD-1 and its ligands. However, the generalizability and reproducibility of the association between circulating sPD-1 levels, cavitary NB lesions, and disease progression beyond a single-country cohort remain uncertain. Therefore, we conducted a multicountry cohort study in East Asia, including patients from Taiwan, Japan, and South Korea, to externally validate the association of blood sPD-1 levels with the cavitary NB phenotype and disease progression requiring treatment in patients with NB-pattern NTM-PD.

MethodsStudy design and patients

This multicountry cohort study was conducted at five referral hospitals across Taiwan, Japan, and South Korea. Patients with NTM-PD were prospectively enrolled in Taiwan from January 2021 to June 2024 and in Japan from June 2023 to April 2024, whereas retrospective data were included from South Korea from September 2011 to August 2016 and from Taiwan from March 2014 to November 2020. Healthy controls without active lung lesions were also bidirectionally included in Taiwan. Importantly, no individual patient-level data from the previously published cohort were reused in the present study. The study was approved by the institutional review boards (No. 2016-07-006BC, 2021-01-010BC, 2023-01-016BC, 2024-08-017AC at Taipei Veterans General Hospital, and 202007084RINA, 202307220RINC at National Taiwan University Hospital, and 202201965A3C501 at Kaohsiung Chang Gung Memorial Hospital in Taiwan, 23014 at Fukujuji Hospital in Japan, and H-2210-096-1369 at Seoul National University Hospital in South Korea).

NTM-PD was defined in patients with respiratory and constitutional symptoms, compatible radiographic findings, and respiratory samples demonstrating growth of the same NTM species on at least two occasions [20]. Patients with NTM-PD exhibiting an NB radiographic pattern, characterized by bronchiectasis and pulmonary nodular lesions, such as nodules, bronchiolitis, and tree-in-bud patterns, were enrolled [15]. Patients with a concomitant FC pattern or coexisting pulmonary tuberculosis were not enrolled because they did not have pure NB-pattern NTM-PD. Patients with organ transplantation, active cancer, or NTM treatment for more than 3 months at enrollment were excluded.

Definitions, measurements, and outcomes

Age, sex, body mass index (BMI), and comorbidities were recorded for all participants. Acid-fast smear and culture results of respiratory samples, NTM species, baseline radiographic patterns, treatment regimens, and treatment initiation dates were recorded for patients with NTM-PD. Cavitary NB pattern was defined as cavity formation within pulmonary nodules or masses, with or without a feeding bronchus [15]. Blood samples were collected from controls and patients at baseline. Using standardized protocols across sites, plasma was stored at −80°C before sPD-1 testing and shipped to a central laboratory for batch analysis. Plasma sPD-1 levels were measured using commercial enzyme-linked immunosorbent assays (ELISA; catalog No. BMS2214, Invitrogen; Thermo Fisher Scientific, Inc.).

Treatment for NTM-PD was defined as initiation of macrolide-based anti-NTM therapy, either azithromycin or clarithromycin, in accordance with the international guideline jointly issued by the American Thoracic Society, European Respiratory Society, European Society of Clinical Microbiology and Infectious Diseases, and Infectious Diseases Society of America [2]. The primary indication for treatment initiation was symptomatic or severe disease after diagnosis, whereas treatment during follow-up was based on radiographic worsening or clinical deterioration due to NTM-PD progression. The decision to initiate treatment was made individually by the primary physician at each study site [12].

All enrolled patients were followed from the date of enrollment until initiation of anti-NTM antibiotic therapy, death, loss to follow-up, or up to 3 years, whichever occurred first [21]. For risk stratification analysis, early disease progression was defined as initiation of anti-NTM treatment within 1 year of enrollment. The primary outcome for prognostic validation was disease progression requiring initiation of anti-NTM treatment within 3 years of baseline assessment. Time to event was calculated from enrollment to outcome occurrence or censoring.

Statistical analysis

Categorical data are expressed as frequencies and percentages and were analyzed using the chi-square test to compare percentages between groups. Continuous variables are expressed as mean±SD for normally distributed data or as median and IQR for nonnormally distributed data. Comparisons were performed using the Student t test or the Mann–Whitney U test, as appropriate, based on the Kolmogorov–Smirnov test. Logistic regression analysis was used for correlation analyses, with results reported as ORs and 95%CIs. For visualization, adjusted predicted probabilities of treatment initiation within 1 year were plotted against sPD-1 levels, with a locally weighted smoothing (LOESS) curve applied to illustrate potential nonlinear associations [22]. Sensitivity analyses were conducted to assess the robustness of the associations [23].

To assess a graded association with outcomes, sPD-1 levels were categorized into tertiles according to their distribution in the overall cohort. The cumulative probability of disease progression during follow-up was compared across sPD-1 subgroups using Kaplan–Meier analysis with the log-rank test. Between-group comparisons of time to study outcome were analyzed using Cox regression, with results reported as adjusted HRs and 95%CIs. The proportional hazards assumption was assessed for all covariates and for the global model using Schoenfeld residuals, with no evidence of significant violation. All statistical analyses were performed using SPSS version 20.0 (SPSS Inc.), and a two-sided P value <.05 was considered statistically significant.

ResultsCharacteristics of the study population

Overall, a total of 68 healthy controls and 400 patients with NTM-PD with an NB radiographic pattern were included (Fig. 1). Among patients with NTM-PD, 302 (75.5%) were infected with MAC, 69 (17.3%) with MABS, and 9 (2.2%) with M. kansasii, whereas 134 (33.5%) were acid-fast smear positive (Table 1). Compared with the control group, the NB-pattern NTM-PD group was older, although the difference did not reach statistical significance: mean age, 63.9±10.4 vs 60.2±16.5 years; P=.079. The NB-pattern NTM-PD group was significantly more female predominant, 74.0% vs 54.4%; P=.002, and had a lower BMI: median, 20.7kg/m2; IQR, 18.7–22.3kg/m2 vs 23.5kg/m2; IQR, 20.9–27.3kg/m2; P<.001. Notably, sPD-1 levels were lower in the NTM-PD group than in the control group: median, 44.2pg/mL; IQR, 26.0–65.4pg/mL vs 61.5pg/mL; IQR, 38.0–81.4pg/mL; P<.001 (Fig. 2A).

Fig. 1.

Flowchart of enrollment. NB, nodular-bronchiectatic pattern; NTM-PD, nontuberculous mycobacterial pulmonary disease.

Table 1.

Characteristics of healthy controls and patients with NB-pattern NTM-PD stratified by cavitary lesion.

Variables  Controls(n=68)  NTM-PD(n=400)  P value  Cavitary NB(n=59)  Noncavitary NB(n=341)  P value 
Age, y  60.2±16.5  63.9±10.4  .079  63.2±9.6  64.0±10.6  .591 
Female sex  37 (54.4%)  296 (74.0%)  .002  47 (79.7%)  249 (73.0%)  .336 
Body mass index, kg/m2  23.5 [20.9–27.3]  20.7 [18.7–22.3]  <.001  20.2 [17.2–21.9]  20.8 [19.1–22.5]  .005 
Comorbidity
Diabetes mellitus  9 (13.2%)  33 (8.3%)  .248  3 (5.1%)  30 (8.8%)  .347 
Cancer history, nonactive  1 (1.5%)  37 (9.3%)  .020  6 (10.2%)  31 (9.1%)  .808 
Previous tuberculosis  3 (4.4%)  67 (16.8%)  .006  7 (11.9%)  60 (17.6%)  .447 
Smear positivity    134 (33.5%)    33 (55.9%)  101 (29.6%)  <.001 
M. avium complex    302 (75.5%)    45 (76.3%)  257 (75.4%)  1.000 
M. abscessus    69 (17.3%)    12 (20.3%)  57 (16.7%)  .575 
Cavitary NB pattern    59 (14.8%)    –  –   
sPD-1 level, pg/mL  61.5 [38.0–81.4]  44.2 [26.0–65.4]  <.001  31.8 [19.6–54.9]  46.9 [27.4–66.7]  .007 
Treatment within 1 y    117 (29.3%)    40 (67.8%)  77 (22.6%)  <.001 

Continuous and categorical data are expressed as mean±SD, median [IQR], or No. (%), as appropriate.

Fig. 2.

(A) sPD-1 levels in healthy controls and in patients with nontuberculous mycobacterial pulmonary disease (NTM-PD) presenting with a nodular-bronchiectatic (NB) pattern. (B) sPD-1 levels in patients with NTM-PD stratified by noncavitary versus cavitary NB pattern. Subgroup medians and IQRs are shown, with each dot representing an individual data point. One outlier in the control group and 9 in the NTM-PD group fall outside the axis range.

Cavitary NB pattern and sPD-1 levels in patients

Among patients with NB-pattern NTM-PD, 59 (14.8%) were identified as having a cavitary NB pattern, whereas the remaining patients had a noncavitary pattern. Lower BMI – median, 20.2kg/m2; IQR, 17.2–21.9kg/m2 vs 20.8kg/m2; IQR, 19.1–22.5kg/m2; P=.005 – higher smear positivity rate, 55.9% vs 29.6%; P<.001, and lower sPD-1 levels – median, 31.8pg/mL; IQR, 19.6–54.9pg/mL vs 46.9pg/mL; IQR, 27.4–66.7pg/mL; P=.007 – were observed in patients with a cavitary NB pattern (Table 1 and Fig. 2B). Of note, sPD-1 levels were lower in female patients than in male patients: median, 42.8pg/mL; IQR, 24.5–61.5pg/mL vs 49.4pg/mL; IQR, 31.4–73.6pg/mL; P=.008. In multivariable logistic analysis, lower sPD-1 levels were independently associated with cavitary NB-pattern NTM-PD, per 10-pg/mL decrease: adjusted OR, 1.141; 95%CI, 1.022–1.274; P=.019, after adjustment for age, sex, BMI, smear positivity, and MAC infection.

Risk stratification analysis: early disease progression requiring treatment within 1 year

Of the 400 patients, 117 (29.3%) were treated with macrolide-based anti-NTM therapy because of severe disease during the first year, including 11 (2.8%) who initiated treatment within 3 months before enrollment. Lower BMI – median, 19.2kg/m2; IQR, 17.4–21.7kg/m2 vs 21.1kg/m2; IQR, 19.5–22.5kg/m2; P<.001 – higher smear positivity rate, 47.0% vs 27.9%; P<.001, higher prevalence of cavitary NB pattern, 34.2% vs 6.7%; P<.001, and lower sPD-1 levels – median, 35.1pg/mL; IQR, 23.3–53.5pg/mL vs 48.4pg/mL; IQR, 29.1–68.1pg/mL; P=.004 – were observed in patients requiring anti-NTM therapy (Table 2). No significant interaction was observed between sPD-1 and sex (P=.146). Although sPD-1 distributions differed across study sites (Fig. 3A), sPD-1 was associated with antibiotic initiation after adjustment for site and covariates. Using probabilities derived from this adjusted model, an inverse association was observed between sPD-1 levels and the probability of antibiotic initiation (Fig. 3B).

Table 2.

Logistic analysis of factors associated with treatment initiation within 1 year in patients with NB-pattern NTM-PD (n=400).

Variables  Treatment initiation within 1 y: Yes(n=117)  Treatment initiation within 1 y: No(n=283)  Univariable analysis OR(95%CI)  P value  Multivariable analysis OR(95%CI)  P value 
Age, y  64.6±10.2  63.5±10.5  1.010 (0.989–1.032)  .336  1.015 (0.991–1.041)  .228 
Female sex  88 (75.2%)  208 (73.5%)  1.094 (0.666–1.796)  .772  1.187 (0.664–2.121)  .563 
Body mass index, kg/m2  19.2 [17.4–21.7]  21.1 [19.5–22.5]  0.808 (0.742–0.881)  <.001  0.823 (0.749–0.905)  <.001 
Comorbidity
Diabetes mellitus  8 (6.8%)  25 (8.8%)  0.757 (0.331–1.732)  .510     
Cancer history  10 (8.5%)  27 (9.5%)  0.941 (0.664–1.376)  .775     
Previous tuberculosis  16 (13.7%)  51 (18.0%)  0.721 (0.392–1.324)  .291     
Smear positivity  55 (47.0%)  79 (27.9%)  2.291 (1.466–3.580)  <.001  1.470 (0.880–2.456)  .141 
M. avium complex  85 (72.6%)  217 (76.7%)  0.808 (0.494–1.320)  .395  0.783 (0.445–1.378)  .396 
M. abscessus  23 (19.7%)  46 (16.3%)  1.261 (0.724–2.195)  .413     
Cavitary NB pattern  40 (34.2%)  19 (6.7%)  7.218 (3.953–13.181)  <.001  5.513 (2.877–10.563)  <.001 
sPD-1 level, pg/mL  35.1 [23.3–53.5]  48.4 [29.1–68.1]  0.990 (0.983–0.997)  .006     
sPD-1 subgroups
>55.5pg/mL (n=133)  26 (22.2%)  107 (37.8%)  1 [Reference]    1 [Reference]   
30.5–55.5pg/mL (n=134)  42 (35.9%)  92 (32.5%)  1.879 (1.070–3.299)  .028  2.026 (1.085–3.782)  .027 
<30.5pg/mL (n=133)  49 (41.9%)  84 (29.7%)  2.401 (1.378–4.181)  .002  2.425 (1.297–4.537)  .006 

Continuous and categorical data are expressed as mean±SD, median [IQR], or No. (%), as appropriate.

Multicollinearity was assessed using variance inflation factors (VIFs). No significant multicollinearity was observed among covariates included in the multivariable models: age, 1.05; female sex, 1.10; body mass index, 1.08; smear positivity, 1.07; M. avium complex, 1.04; cavitary NB pattern, 1.09; and sPD-1, 1.03.

95%CI, 95% confidence interval; BMI, body mass index; IQR, interquartile range; NB, nodular-bronchiectatic; NTM-PD, nontuberculous mycobacterial pulmonary disease; OR, odds ratio; SD, standard deviation; sPD-1, soluble programmed death-1; VIF, variance inflation factor.

Fig. 3.

(A) Distribution of circulating sPD-1 levels in patients with nontuberculous mycobacterial pulmonary disease (NTM-PD) with a nodular-bronchiectatic (NB) pattern, stratified by study site. Data are shown as individual values with subgroup medians and IQRs. sPD-1 distributions differed across cohorts, with intermediate levels in Taiwan, median, 37.4pg/mL; IQR, 24.3–63.3pg/mL; lower values in Japan, median, 20.4pg/mL; IQR, 14.5–28.2pg/mL; and higher values in South Korea, median, 54.5pg/mL; IQR, 42.0–71.7pg/mL. Nine values exceeded the displayed axis range. (B) Adjusted probability of antibiotic treatment initiation within 1 year according to sPD-1 levels. Predicted probabilities were derived from a multivariable logistic regression model including sPD-1 levels and adjusting for study site – Japan and South Korea, with Taiwan as the reference – age, sex, BMI, smear positivity, MAC species, and cavitary disease. In this model, lower sPD-1 levels were independently associated with antibiotic initiation within 1 year: P=.043. The curve represents a locally weighted smoothing (LOESS) fit, demonstrating an inverse and nonlinear relationship. A small number of values exceeded the plotted axis range and are not shown.

When patients were stratified into tertiles based on sPD-1 levels, the proportion requiring anti-NTM therapy in the first year was lowest in the subgroup with sPD-1 >55.5pg/mL (n=133), higher in the subgroup with sPD-1 30.5–55.5pg/mL (n=134), and highest in the subgroup with sPD-1 <30.5pg/mL (n=133): 19.5% vs 31.3% vs 36.8%; P for trend=.002. In multivariable analysis, antibiotic initiation was independently associated with lower BMI, per 1-kg/m2 increase: OR, 0.823; 95%CI, 0.749–0.905; P<.001; cavitary NB pattern: OR, 5.513; 95%CI, 2.877–10.563; P<.001; low sPD-1 level, <30.5pg/mL: OR, 2.425; 95%CI, 1.297–4.537; P=.006; and intermediate sPD-1 level, 30.5–55.5pg/mL: OR, 2.026; 95%CI, 1.085–3.782; P=.027, compared with high sPD-1 level, >55.5pg/mL (Table 2). The association between sPD-1 and antibiotic initiation remained consistent after site-period adjustment, whether sPD-1 was modeled as a continuous or tertile-based variable, and in subgroup analyses (Table 3).

Table 3.

Sensitivity analyses assessing the association between sPD-1 level and treatment initiation within 1 year.

Model/scenarioa  sPD-1 tertiles  Tertile model OR(95%CI)  Continuous model OR per 10pg/mL(95%CI)  P value 
1. Base model  >55.5pg/mL (n=133)  1 [Reference]  0.904 (0.835–0.978)  .012 
n=400  30.5–55.5pg/mL (n=134)  2.026 (1.085–3.782)     
  <30.5pg/mL (n=133)  2.425 (1.297–4.537)     
2. Base model+site adjustmentb  >55.5pg/mL (n=133)  1 [Reference]  0.916 (0.842–0.997)  .014 
n=400  30.5–55.5pg/mL (n=134)  2.072 (1.077–3.985)     
  <30.5pg/mL (n=133)  2.024 (0.972–4.215)     
3. Base model+site-period adjustmentc  >55.5pg/mL (n=133)  1 [Reference]  0.933 (0.864–1.007)  .075 
n=400  30.5–55.5pg/mL (n=134)  2.006 (1.044–3.855)     
  <30.5pg/mL (n=133)  1.813 (0.939–3.501)     
4. Noncavitary NB subgroup  >55.5pg/mL (n=119)  1 [Reference]  0.899 (0.821–0.983)  .020 
n=341  30.5–55.5pg/mL (n=117)  2.093 (1.061–4.131)     
  <30.5pg/mL (n=105)  2.466 (1.227–4.955)     
5. M. avium complex subgroup  >55.5pg/mL (n=104)  1 [Reference]  0.924 (0.851–1.002)  .057 
n=302  30.5–55.5pg/mL (n=96)  1.670 (0.803–3.471)     
  <30.5pg/mL (n=102)  2.076 (1.013–4.256)     

95%CI, 95% confidence interval; BMI, body mass index; NB, nodular-bronchiectatic; OR, odds ratio; sPD-1, soluble programmed death-1.

a

Models were adjusted for age, sex, BMI, smear positivity, M. avium complex infection, and cavitary disease unless otherwise specified.

b

Study site was modeled using indicator variables: Japan and South Korea, with Taiwan as the reference.

c

Site-period was modeled as a combined 4-level variable: Taiwan prospective, Taiwan retrospective, Japan prospective, and South Korea retrospective.

Prognostic validation: risk of disease progression requiring treatment within 3 years

For the primary analysis, patients who initiated anti-NTM treatment before enrollment were excluded (n=11). Among the remaining 389 patients, anti-NTM therapy was initiated in 167 (42.9%) during follow-up, with a median follow-up duration of 1.3 years (IQR, 0.3–5.8 years). Among these cases, 129 (77.7%) were treated within 3 years. In Kaplan–Meier analysis, the cumulative probability of disease progression requiring anti-NTM therapy within 3 years was lowest in patients with sPD-1 >55.5pg/mL, higher in those with sPD-1 30.5–55.5pg/mL, and highest in those with sPD-1 <30.5pg/mL; log-rank test, P=.003 (Fig. 4).

Fig. 4.

Kaplan–Meier curves showing the cumulative probability of disease progression requiring anti-NTM therapy within 3 years in patients with nontuberculous mycobacterial pulmonary disease (NTM-PD) with a nodular-bronchiectatic (NB) pattern, stratified by sPD-1 levels (n=389 after excluding 11 patients with antibiotic treatment before enrollment). Adjusted HRs and 95%CIs from the Cox regression analysis are also shown.

In multivariable Cox regression analysis, disease progression was independently associated with lower BMI, per 1-kg/m2 increase: HR, 0.901; 95%CI, 0.842–0.963; P=.002; cavitary NB pattern: HR, 3.262; 95%CI, 2.166–4.914; P<.001; low sPD-1 level, <30.5pg/mL: HR, 1.957; 95%CI, 1.241–3.087; P=.004; and intermediate sPD-1 level, 30.5–55.5pg/mL: HR, 1.583; 95%CI, 1.007–2.488; P=.047, compared with high sPD-1 level, >55.5pg/mL, after adjustment for age, sex, BMI, smear positivity, and MAC infection (Table 4).

Table 4.

Cox regression analysis of factors associated with treatment initiation within 3 years in patients with NB-pattern NTM-PD (n=389).

Variables  Event/follow-up years  Univariable analysis HR(95%CI)  P value  Multivariable analysis HR(95%CI)  P value 
sPD-1 subgroups
>55.5pg/mL (n=132)  33/3.0 [0.8–3.0]  1 [Reference]    1 [Reference]   
30.5–55.5pg/mL (n=127)  47/1.9 [0.5–3.0]  1.559 (0.998–2.433)  .051  1.583 (1.007–2.488)  .047 
<30.5pg/mL (n=130)  49/0.9 [0.2–1.4]  2.121 (1.356–3.316)  .001  1.957 (1.241–3.087)  .004 
Age, y    1.005 (0.988–1.021)  .588  1.003 (0.985–1.021)  .783 
Female sex    1.014 (0.683–1.506)  .945  1.243 (0.819–1.887)  .306 
Body mass index, kg/m2    0.883 (0.829–0.941)  <.001  0.901 (0.842–0.963)  .002 
Comorbidity
Diabetes mellitus    0.756 (0.384–1.489)  .510     
Cancer history    0.922 (0.509–1.671)  .790     
Previous tuberculosis    0.795 (0.494–1.282)  .347     
Disease characteristics
Smear positivity    1.905 (1.346–2.696)  <.001  1.412 (0.977–2.041)  .067 
M. avium complex    1.097 (0.725–1.659)  .661  1.142 (0.747–1.748)  .539 
M. abscessus    0.957 (0.599–1.527)  .852     
Cavitary NB pattern    4.004 (2.720–5.892)  <.001  3.262 (2.166–4.914)  <.001 

Data in the event/follow-up years column are expressed as no. of events/median follow-up years [IQR].

95%CI, 95% confidence interval; HR, hazard ratio; IQR, interquartile range; NB, nodular-bronchiectatic; NTM-PD, nontuberculous mycobacterial pulmonary disease; sPD-1, soluble programmed death-1.

Discussion

In this large multicountry cohort study, patients with NB-pattern NTM-PD had lower sPD-1 levels than healthy controls. Among patients with NB-pattern NTM-PD, those with cavitary NB pattern, who represented 15% of the cohort, had even lower sPD-1 levels than those with noncavitary NB pattern. During follow-up, disease progression requiring anti-NTM therapy occurred in 167 patients (42.9%). In the primary time-to-event analysis, lower baseline sPD-1 levels were independently associated with an increased risk of disease progression requiring treatment within 3 years. Directionally consistent associations were also observed in the analysis of early disease progression within 1 year, supporting the potential prognostic role of sPD-1. Collectively, these findings externally validate sPD-1 as a marker associated with the cavitary NB phenotype at baseline and subsequent disease progression, supporting its potential role in risk stratification in patients with NB-pattern NTM-PD.

These findings are consistent with previous evidence suggesting that immune exhaustion, a key determinant of NTM-PD severity, may be associated with sPD-1 levels [15,24]. Our study demonstrated a graded association between sPD-1 levels and disease progression, with patients in the lowest sPD-1 tertile showing the greatest risk of requiring anti-NTM therapy. This finding may be explained by the immunomodulatory role of sPD-1 as a PD-1 pathway antagonist, which may be associated with impaired T-cell function and disease progression [24,25]. In multivariable analysis adjusted for other significant factors, low sPD-1 remained independently associated with disease progression, suggesting an association beyond potential confounding. Although the causal relationship and underlying mechanism linking sPD-1 to NTM-PD severity and progression warrant further investigation, our multicountry validation cohort in Asia supports the use of sPD-1 at diagnosis as a prognostic marker for disease progression and suggests potential clinical applicability.

Recent studies suggest that chronic NTM-PD, possibly through persistent antigenic stimulation, upregulates PD-1 expression on CD4+ T cells, which is associated with higher bacterial burden, positive sputum smears, and more severe radiographic manifestations, including cavitary lesions [18,26]. Because sPD-1 is a circulating form of the immune checkpoint molecule PD-1, generated by alternative splicing and acting as a decoy receptor, it may mitigate the inhibitory effects of ligand binding to membrane-bound PD-1 on T-cell function [27,28]. Consequently, lower plasma sPD-1 levels may reflect greater T-cell exhaustion and impaired cellular immunity [29]. Furthermore, immune exhaustion impairs the host's capacity to generate effective T helper 1 responses and produce interferon γ, both of which are critical for controlling mycobacterial infection [30].

Interestingly, lower sPD-1 levels were associated not only with disease progression requiring anti-NTM therapy but also independently with the presence of cavitary NB-pattern NTM-PD, even after adjustment for age, sex, BMI, smear positivity, and NTM species. This finding suggests that immune exhaustion, reflected by low sPD-1 levels, may be associated with cavitary disease in NB-pattern NTM-PD, although the underlying mechanism remains unclear [15,18]. Nevertheless, the independent association of low sPD-1 with cavitary disease, even after adjustment for other risk factors, underscores its potential as a biomarker of advanced immune dysfunction in NTM-PD.

Unlike previous studies that primarily emphasized mortality prediction, such as the BACES score [31], our findings highlight sPD-1 as a novel predictor of disease progression, complementing established factors such as cavity formation and acid-fast smear positivity [2]. Most NB-pattern cases are noncavitary at diagnosis, approximately 85% in our study, and their clinical trajectory remains difficult to predict. Notably, the current ATS/ERS/ESCMID/IDSA guideline for NTM-PD management provides no clear recommendations on how to conduct watchful waiting in patients with noncavitary NB-pattern disease who do not receive anti-NTM therapy [2]. Therefore, identifying low sPD-1 as a risk factor for antibiotic initiation in patients with noncavitary NB-pattern NTM-PD has important implications for risk stratification and management. Specifically, patients with low sPD-1 at baseline are at increased risk of rapid progression within 1 year and may benefit from closer monitoring and earlier consideration of therapy [32]. Low plasma sPD-1, low BMI, high smear grade, and cavitary NB pattern are robust predictors of disease progression [32]. A composite score integrating these factors may provide a practical indicator of disease progression and warrants validation in larger studies.

The strengths of this study include its large multicountry cohort design and the standardized definition of antibiotic initiation as a clinically relevant surrogate outcome of disease progression. However, this outcome may be influenced by physician and patient decision-making, local practice patterns, censoring due to loss to follow-up, and potential competing risks, such as death, which were not systematically available across cohorts.

Several limitations should be acknowledged. First, the cohort included both retrospective and prospective components across sites, which may have introduced residual confounding due to cross-country differences in diagnostic workup, follow-up, and treatment thresholds. Nevertheless, the association between sPD-1 and disease progression remained consistent after site-period adjustment, whether sPD-1 was modeled as a continuous or tertile-based variable. However, the tertile cutoffs were derived from the exploratory distribution of sPD-1 levels in our cohort and require further validation. Because all patients in the prospective Japan cohort had follow-up durations shorter than 3 years at the time of analysis, site- and period-adjusted Cox regression analysis for the primary 3-year endpoint was not performed.

Second, detailed radiographic characterization of cavitary lesions, such as size, wall thickness, and anatomical distribution, was not available; such information could refine the phenotypic correlation with sPD-1.

Third, sPD-1 was measured only at baseline; therefore, we were unable to assess dynamic changes over time or in response to treatment. In addition, residual preanalytical variability related to storage duration and sample handling across sites could not be fully excluded despite centralized measurement and standardized protocols.

Fourth, healthy controls were recruited from a single country, whereas patients were enrolled from multiple countries, which may limit comparability of absolute sPD-1 levels between groups.

Fifth, we did not assess key inflammatory markers, such as cytokines, chemokines, C-reactive protein, and albumin, which may have confounded the association between sPD-1 and disease progression.

Conclusions

This study demonstrates that lower baseline sPD-1 levels are associated with both cavitary lesions and increased risk of antibiotic initiation, a clinically relevant surrogate outcome of disease progression, in NB-pattern NTM-PD. These findings support the potential role of sPD-1 as a prognostic biomarker for risk stratification and closer monitoring in this challenging and increasingly prevalent form of NTM-PD. Future studies should validate sPD-1 in broader populations and explore its utility in guiding treatment initiation and monitoring.

Authors’ contributions

All authors contributed to sample and data collection. Drs Pan, Furuuchi, Shu, Morimoto, and Yim performed the data analysis. Drs Pan and Shu drafted the manuscript. Profs Morimoto and Yim critically revised the manuscript. Drs Shu, Morimoto, and Yim are the guarantors of the paper and take responsibility for the integrity of the work as a whole, from inception to publication.

Declaration of generative AI and AI-assisted technologies in the writing process

The authors declare that no artificial intelligence–assisted tools were used in the drafting, analysis, or writing of this manuscript.

ClinicalTrials.gov identifier

NCT05678166.

Conference presentation

This work includes data that were previously presented in part at the 2024 Annual Congress of the Taiwan Society of Pulmonary and Critical Care Medicine and in part at the 2025 International Conference of the American Thoracic Society.

Funding

This study was supported by the APSR – Asia Pacific Society of Respirology – Multi-Society Research Project Fund 2023–2024, Taipei Veterans General Hospital (V113E-006-1, V114E-007-2, V115E-006-2), National Taiwan University Hospital (115-E0007), and the National Science and Technology Council, Taiwan (114-2314-B-075-079).

Conflicts of interest

None declared.

Data availability

The data supporting the findings of this study are available from the corresponding author upon reasonable request. The data are not publicly available because of privacy or ethical restrictions.

Acknowledgments

The authors thank the Medical Science and Technology Building of the Taipei Veterans General Hospital for providing research facilities. The authors thank the 2nd and 8th core labs of the Department of Medical Research, National Taiwan University Hospital, for providing research space, technical services.

References
[1]
C.S. Haworth, J. Banks, T. Capstick, A.J. Fisher, T. Gorsuch, I.F. Laurenson, et al.
British Thoracic Society guidelines for the management of non-tuberculous mycobacterial pulmonary disease (NTM-PD).
[2]
C.L. Daley, J.M. Iaccarino, C. Lange, E. Cambau, R.J. Wallace, C. Andrejak, et al.
Treatment of nontuberculous mycobacterial pulmonary disease: an official ATS/ERS/ESCMID/IDSA clinical practice guideline.
Clin Infect Dis, 71 (2020), pp. 905-913
[3]
K. Hyung, S.A. Kim, J.Y. Kim, N. Kwak, J.J. Yim.
Rates and risk factors of progression in patients with nontuberculous mycobacterial pulmonary disease: secondary analysis of a prospective cohort study.
[4]
V.N. Dahl, L.L. Laursen, Y. He, Y.A. Zhang, M.S. Wang.
Species distribution among patients with nontuberculous mycobacteria pulmonary disease in Europe.
J Infect, 87 (2023), pp. 469-472
[5]
A.A. Pedersen, A. Løkke, A. Fløe, R. Ibsen, I.S. Johansen, O. Hilberg.
Nationwide increasing incidence of nontuberculous mycobacterial diseases among adults in Denmark: eighteen years of follow-up.
[6]
K.L. Winthrop, T.K. Marras, J. Adjemian, H. Zhang, P. Wang, Q. Zhang.
Incidence and prevalence of nontuberculous mycobacterial lung disease in a large U.S. managed care health plan, 2008–2015.
Ann Am Thorac Soc, 17 (2020), pp. 178-185
[7]
Y. Hamaguchi, K. Morimoto, S. Mitarai.
Laboratory-based surveillance of nontuberculous mycobacterial pulmonary disease in Japan.
ERJ Open Res, 11 (2025),
[8]
J.Y. Kim, N. Kwak, J.J. Yim.
The rise in prevalence and related costs of nontuberculous mycobacterial diseases in South Korea, 2010–2021.
Open Forum Infect Dis, 9 (2022),
[9]
B.W. Jhun, S.M. Moon, K. Jeon, O.J. Kwon, H. Yoo, K.C. Carriere, et al.
Prognostic factors associated with long-term mortality in 1445 patients with nontuberculous mycobacterial pulmonary disease: a 15-year follow-up study.
Eur Respir J, 55 (2020),
[10]
P.H. Wang, S.W. Pan, S.M. Wang, C.C. Shu, C.H. Chang.
The impact of nontuberculous mycobacteria species on mortality in patients with nontuberculous mycobacterial lung disease.
Front Microbiol, 13 (2022),
[11]
W.J. Koh, S.M. Moon, S.Y. Kim, M.A. Woo, S. Kim, B.W. Jhun, et al.
Outcomes of Mycobacterium avium complex lung disease based on clinical phenotype.
Eur Respir J, 50 (2017),
[12]
S.W. Pan, C.C. Shu, J.Y. Feng, W.J. Su.
Treatment for Mycobacterium avium complex lung disease.
J Formos Med Assoc, 119 (2020), pp. S67-S75
[13]
S. Cowman, J. van Ingen, D.E. Griffith, M.R. Loebinger.
Nontuberculous mycobacterial pulmonary disease.
Eur Respir J, 54 (2019),
[14]
S.W. Pan, C.C. Shu, J.Y. Feng, J.Y. Wang, Y.J. Chan, C.J. Yu, et al.
Microbiological persistence in patients with Mycobacterium avium complex lung disease: the predictors and the impact on radiographic progression.
Clin Infect Dis, 65 (2017), pp. 927-934
[15]
S.W. Pan, W.J. Su, Y.J. Chan, M.L. Ho, J.Y. Feng, C.C. Shu, et al.
Disease progression in patients with nontuberculous mycobacterial lung disease of nodular bronchiectatic pattern: the roles of cavitary NB and soluble programmed death protein-1.
Clin Infect Dis, 75 (2022), pp. 239-247
[16]
M. Ito, K. Furuuchi, K. Fujiwara, T. Kodama, Y. Tanaka, T. Yoshiyama, et al.
Epidemiological trends and clinical relevance of nontuberculous mycobacterial pulmonary disease in a referral hospital in Japan, 2017–2021.
Respir Investig, 62 (2024), pp. 1064-1071
[17]
H.L. Huang, M.H. Cheng, P.L. Lu, C.C. Shu, J.Y. Wang, J.T. Wang, et al.
Epidemiology and predictors of NTM pulmonary infection in Taiwan: a retrospective, 5-year multicenter study.
Sci Rep, 7 (2017),
[18]
C.C. Shu, S.W. Pan, J.Y. Feng, J.Y. Wang, Y.J. Chan, C.J. Yu, et al.
The clinical significance of programmed death-1, regulatory T cells, and myeloid-derived suppressor cells in patients with nontuberculous mycobacteria lung disease.
J Clin Med, 8 (2019), pp. 736
[19]
C. Nielsen, L. Ohm-Laursen, T. Barington, S. Husby, S.T. Lillevang.
Alternative splice variants of the human PD-1 gene.
Cell Immunol, 235 (2005), pp. 109-116
[20]
C.L. Daley, J.M. Iaccarino, C. Lange, E. Cambau, R.J. Wallace Jr., C. Andrejak, et al.
Treatment of nontuberculous mycobacterial pulmonary disease: an official ATS/ERS/ESCMID/IDSA clinical practice guideline.
Clin Infect Dis, 71 (2020), pp. e1-e36
[21]
S.W. Pan, C.C. Shu, J.Y. Feng, J.Y. Chien, J.Y. Wang, Y.J. Chan, et al.
Impact of different subspecies on disease progression in initially untreated patients with Mycobacterium avium complex lung disease.
Clin Microbiol Infect, 27 (2021),
[22]
W. Zhao, H. Tan, Y. Yan, J. Zhang.
Nonlinear association between sleep duration and thyroid hormone levels in patients with thyroid cancer.
Sci Rep, 16 (2026),
[23]
S.W. Pan, Y.F. Yen, J.Y. Feng, P.H. Chuang, V.Y. Su, Y.R. Kou, et al.
Opposite effects of statins on the risk of tuberculosis and herpes zoster in patients with diabetes: a population-based cohort study.
Br J Clin Pharmacol, 86 (2020), pp. 569-579
[24]
P.K. Marty, B. Pathakumari, T.M. Cox, V.P. Van Keulen, C.L. Erskine, M. Shah, et al.
Multiparameter immunoprofiling for the diagnosis and differentiation of progressive versus nonprogressive nontuberculous mycobacterial lung disease: a pilot study.
PLoS One, 19 (2024),
[25]
C.C. Shu, J.Y. Wang, M.F. Wu, C.T. Wu, H.C. Lai, L.N. Lee, et al.
Attenuation of lymphocyte immune responses during Mycobacterium avium complex-induced lung disease due to increasing expression of programmed death-1 on lymphocytes.
Sci Rep, 7 (2017),
[26]
C.C. Shu, M.F. Wu, S.W. Pan, T.S. Wu, H.C. Lai, M.C. Lee.
Host immune response against environmental nontuberculous mycobacteria and the risk populations of nontuberculous mycobacterial lung disease.
J Formos Med Assoc, 119 (2020), pp. S13-S22
[27]
M. Wahid, B. Pratoomthai, I.U. Egbuniwe, H.R. Evans, R. Babaei-Jadidi, J.O. Amartey, et al.
Targeting alternative splicing as a new cancer immunotherapy: phosphorylation of serine/arginine-rich splicing factor 1 by SR protein kinase 1 regulates alternative splicing of PD1 to generate a soluble antagonistic isoform that prevents T-cell exhaustion.
Cancer Immunol Immunother, 72 (2023), pp. 4001-4014
[28]
M. Niu, Y. Liu, M. Yi, D. Jiao, K. Wu.
Biological characteristics and clinical significance of soluble PD-1/PD-L1 and exosomal PD-L1 in cancer.
Front Immunol, 13 (2022), pp. 827921
[29]
A. Gramegna, A. Lombardi, N.I. Lore, F. Amati, I. Barone, C. Azzara, et al.
Innate and adaptive lymphocytes in nontuberculous mycobacteria lung disease: a review.
Front Immunol, 13 (2022), pp. 927049
[30]
C.N. Ratnatunga, K. Tungatt, C. Proietti, S. Halstrom, M.R. Holt, V.P. Lutzky, et al.
Characterizing and correcting immune dysfunction in nontuberculous mycobacterial disease.
Front Immunol, 13 (2022),
[31]
H.J. Kim, N. Kwak, H. Hong, N. Kang, Y. Im, B.W. Jhun, et al.
BACES score for predicting mortality in nontuberculous mycobacterial pulmonary disease.
Am J Respir Crit Care Med, 203 (2021), pp. 230-236
[32]
H. Hwang, J.K. Lee, E.Y. Heo, D.K. Kim, H.W. Lee.
The factors associated with mortality and progressive disease of nontuberculous mycobacterial lung disease: a systematic review and meta-analysis.

These authors contributed equally as first authors.

These authors contributed equally to this work as co-corresponding authors.

Copyright © 2026. SEPAR
Download PDF
Archivos de Bronconeumología
Article options
Tools