, Ji Eun Kim1
, Su Jin Jeong2
, Chang Kyun Lee1
, the Big Data Research Group (BDRG) of the Korean Society of Gastroenterology 1Department of Gastroenterology, Center for Crohn’s and Colitis, Kyung Hee University Hospital, College of Medicine, Kyung Hee University, Seoul, Korea
2Department of Statistics Support, Medical Science Research Institute, Kyung Hee University Medical Center, Seoul, Korea
© 2026 Korean Association for the Study of Intestinal Diseases.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
Funding Source
This research was supported by a grant of the Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI), funded by the Ministry of Health & Welfare, Republic of Korea (grant number: RS-2023- KH135855). This study was conducted as collaborative research project between the Korean Society of Gastroenterology and the National Health Insurance Service in Korea, with support from both institutions.
Conflict of Interest
No potential conflict of interest relevant to this article was reported.
Data Availability Statement
The findings of this study are based on secondary, processed, nationwide data that could not be stored or publicly archived due to their size and resource limitations. The Korean National Health Insurance Service provided only the analysis results. The raw data can be accessed through the Korean National Health Insurance Sharing Service (https://nhiss/nhis/or/kr), subject to applicable requirements and associated fees.
Author Contributions
Conceptualization; Lee CK, Oh SJ. Data curation: Oh SJ, Jeong SJ. Data interpretation: Oh SJ, Jeong SJ. Formal analysis: Jeong SJ. Investigation: all authors. Methodology; Project administration; Lee CK, Oh SJ, Jeong SJ. Resources; Software; Supervision; Lee CK, Jeong SJ. Validation; Visualization: Oh SJ, Jeong SJ. Writing–original draft: all authors. Writing–review & editing: all authors. Approval of final manuscript: all authors.
|
Simple |
Multivariablea |
Matchingb |
||||
|---|---|---|---|---|---|---|
| OR (95% CI) | P-value | OR (95% CI) | P-value | OR (95% CI) | P-value | |
| Overall | ||||||
| Case | 1.41 (1.33–1.51) | <0.001 | 1.52 (1.39–1.66) | <0.001 | 1.65 (1.44–1.89) | <0.001 |
| Control | 1.00 | 1.00 | 1.00 | |||
| CNS disorders | ||||||
| Case | 1.38 (1.11–1.72) | 0.004 | 1.19 (0.87–1.61) | 0.271 | 0.94 (0.57–1.54) | 0.796 |
| Control | 1.00 | 1.00 | 1.00 | |||
| ILD | ||||||
| Case | 1.15 (0.97–1.36) | 0.101 | 2.10 (1.70–2.59) | <0.001 | 1.88 (1.33–2.66) | <0.001 |
| Control | 1.00 | 1.00 | 1.00 | |||
| Psoriasis | ||||||
| Case | 1.55 (1.43–1.67) | <0.001 | 1.54 (1.39–1.71) | <0.001 | 1.58 (1.34–1.86) | <0.001 |
| Control | 1.00 | 1.00 | 1.00 | |||
| Lupus | ||||||
| Case | 1.13 (0.97–1.31) | 0.118 | 1.14 (0.94–1.38) | 0.182 | 1.32 (0.93–1.88) | 0.120 |
| Control | 1.00 | 1.00 | 1.00 | |||
| Vasculitis | ||||||
| Case | 1.34 (0.94–1.90) | 0.107 | 1.23 (0.76–1.99) | 0.401 | 1.48 (0.65–3.38) | 0.351 |
| Control | 1.00 | 1.00 | 1.00 | |||
a Estimates were adjusted for age, sex, IBD subtype, IBD diagnosis year, body mass index, smoking, frequency of hospitalization due to IBD, IBD-related surgeries, Charlson Comorbidity Index, and presence of concomitant immune-mediated inflammatory diseases (e.g., rheumatoid arthritis and ankylosing spondylitis).
b Estimates were derived from propensity score matched analysis followed by multivariable adjustment using the same covariates.
TNF, tumor necrosis factor; IBD, inflammatory bowel disease; OR, odds ratio; CI, confidence interval; CNS, central nervous system; ILD, interstitial lung disease.
|
Multivariablea |
Matchingb |
|||||
|---|---|---|---|---|---|---|
|
Immunosuppressants not use |
Immunosuppressant use |
Interaction P-value |
Immunosuppressants not use |
Immunosuppressant use |
Interaction P-value | |
| OR (95% CI) | OR (95% CI) | OR (95% CI) | OR (95% CI) | |||
| Overall | 1.32 (1.02–1.69) | 1.12 (1.02–1.23) | 0.003 | 1.38 (0.96–2.00) | 1.26 (1.08–1.48) | 0.244 |
| CNS disorders | 1.14 (0.47–2.74) | 1.03 (0.74–1.45) | 0.571 | 1.06 (0.24–4.81) | 0.74 (0.42–1.31) | 0.527 |
| ILD | 1.70 (0.84–3.44) | 1.21 (0.96–1.51) | 0.144 | 1.05 (0.40–2.71) | 1.19 (0.08–1.77) | 0.970 |
| Psoriasis | 1.60 (1.19–2.15) | 1.17 (1.04–1.32) | 0.014 | 1.64 (1.05–2.56) | 1.26 (1.04–1.52) | 0.246 |
| Lupus | 0.52 (0.27–0.99) | 0.89 (0.73–1.08) | 0.271 | 0.63 (0.21–1.91) | 1.04 (0.71–1.53) | 0.597 |
| Vasculitis | <0.01 (0.00–1,000) | 1.03 (0.63–1.68) | 0.970 | <0.01 (0.00–1,000) | 1.27 (0.47–3.40) | 0.968 |
a Estimates were adjusted for age, sex, IBD subtype, IBD diagnosis year, body mass index, smoking, frequency of hospitalization due to IBD, IBD-related surgeries, Charlson Comorbidity Index, and presence of concomitant immune-mediated inflammatory diseases (e.g., rheumatoid arthritis and ankylosing spondylitis).
b Estimates were derived from propensity score matched analysis followed by multivariable adjustment using the same covariates.
TNF, tumor necrosis factor; OR, odds ratio; CI, confidence interval; CNS, central nervous system; ILD, interstitial lung disease; IBD, inflammatory bowel disease.
|
Simple |
Multivariablea |
Matchingb |
||||
|---|---|---|---|---|---|---|
| OR (95% CI) | P-value | OR (95% CI) | P-value | OR (95% CI) | P-value | |
| Overall | ||||||
| Anti-TNF-α monotherapy | Reference | Reference | Reference | |||
| Anti-TNF-α no use | 0.85 (0.72–1.00) | 0.050 | 0.93 (0.76–1.13) | 0.460 | 0.82 (0.61–1.11) | 0.193 |
| Concomitant use of immunosuppressants | 1.23 (1.03–1.47) | 0.020 | 1.51 (1.22–1.87) | <0.001 | 1.43 (1.04–1.98) | 0.029 |
a Estimates were adjusted for age, sex, IBD subtype, IBD diagnosis year, body mass index, smoking, frequency of hospitalization due to IBD, IBD-related surgeries, Charlson Comorbidity Index, and presence of concomitant immune-mediated inflammatory diseases (e.g., rheumatoid arthritis and ankylosing spondylitis).
b Estimates were derived from propensity score matched analysis followed by multivariable adjustment using the same covariates.
TNF, tumor necrosis factor; OR, odds ratio; CI, confidence interval; IBD, inflammatory bowel disease.
| Variable | Unmatched |
1:1 Matched |
||||
|---|---|---|---|---|---|---|
| Case (n=8,586) | Control (n=111,153) | P-value | Case (n=8,586) | Control (n=8,586) | P-value | |
| Age at index date (yr) | 47 (33–60) | 40 (27–54) | <0.001 | 47 (33–60) | 47 (33–60) | 0.824 |
| Age at index date (yr) | <0.001 | 0.970 | ||||
| ≤ 19 | 622 (7.2) | 11,193 (10.1) | 622 (7.1) | 606 (7.1) | ||
| 20–39 | 2,418 (28.2) | 42,804 (38.5) | 2,418 (28.2) | 2,431 (28.3) | ||
| 40–59 | 3,328 (38.8) | 37,412 (33.7) | 3,328 (38.8) | 3,333 (38.8) | ||
| ≥ 60 | 2,218 (25.8) | 19,744 (17.8) | 2,218 (25.8) | 2,216 (25.8) | ||
| Subtype | <0.001 | 0.848 | ||||
| Ulcerative colitis | 5,917 (68.9) | 79,054 (71.1) | 5,917 (68.9) | 5,951 (69.3) | ||
| Crohn’s disease | 1,529 (17.8) | 18,951 (17.1) | 1,529 (17.8) | 1,513 (17.6) | ||
| Unclassified | 1,140 (13.3) | 13,148 (11.8) | 1,140 (13.3) | 1,121 (13.1) | ||
| Sex | <0.001 | 0.951 | ||||
| Male | 4,366 (50.9) | 64,393 (57.9) | 4,366 (50.9) | 4,362 (50.8) | ||
| Female | 4,220 (49.2) | 46,760 (42.1) | 4,220 (49.2) | 4,224 (49.2) | ||
| BMI (kg/m²) | 24.2 (22.1–26.4) | 24.1 (21.9–26.3) | 0.003 | 24.2 (22.1–26.4) | 24.1 (21.9–26.3) | 0.091 |
| Charlson Comorbidity Index score | 4 (2–6) | 2 (1–4) | <0.001 | 3 (2–5) | 4 (2–6) | <0.001 |
| Disease duration of IBD (yr) | 4.17 (1.55–7.68) | – | 4.17 (1.55–7.68) | – | ||
| Frequency of hospitalizations due to IBD | 9.96±42.91 | 6.73±27.20 | 9.96±42.91 | 7.28±27.68 | <0.001 | |
| Immunosuppressant used | 4,742 (55.2) | 33,996 (30.6) | <0.001 | 4,742 (55.2) | 2,880 (33.5) | <0.001 |
| Anti-TNF-α used | 1,183 (13.8) | 11,294 (10.2) | <0.001 | 1,183 (13.8) | 735 (8.6) | <0.001 |
| Infliximab | 782 (9.1) | 8,029 (7.2) | 782 (9.1) | 527 (6.1) | ||
| Adalimumab | 580 (6.8) | 4,569 (4.1) | 580 (6.8) | 313 (3.7) | ||
| Golimumab | 89 (1.0) | 613 (0.6) | 89 (1.0) | 46 (0.5) | ||
| Cumulative duration of anti-TNF-α (yr) | 3.02 (1.14–7.68) | 2.54 (1.00–4.84) | <0.001 | 3.02 (1.14–5.66) | 2.84 (1.02–5.61) | 0.587 |
| Simple |
Multivariable |
Matching |
||||
|---|---|---|---|---|---|---|
| OR (95% CI) | P-value | OR (95% CI) | P-value | OR (95% CI) | P-value | |
| Overall | ||||||
| Case | 1.41 (1.33–1.51) | <0.001 | 1.52 (1.39–1.66) | <0.001 | 1.65 (1.44–1.89) | <0.001 |
| Control | 1.00 | 1.00 | 1.00 | |||
| CNS disorders | ||||||
| Case | 1.38 (1.11–1.72) | 0.004 | 1.19 (0.87–1.61) | 0.271 | 0.94 (0.57–1.54) | 0.796 |
| Control | 1.00 | 1.00 | 1.00 | |||
| ILD | ||||||
| Case | 1.15 (0.97–1.36) | 0.101 | 2.10 (1.70–2.59) | <0.001 | 1.88 (1.33–2.66) | <0.001 |
| Control | 1.00 | 1.00 | 1.00 | |||
| Psoriasis | ||||||
| Case | 1.55 (1.43–1.67) | <0.001 | 1.54 (1.39–1.71) | <0.001 | 1.58 (1.34–1.86) | <0.001 |
| Control | 1.00 | 1.00 | 1.00 | |||
| Lupus | ||||||
| Case | 1.13 (0.97–1.31) | 0.118 | 1.14 (0.94–1.38) | 0.182 | 1.32 (0.93–1.88) | 0.120 |
| Control | 1.00 | 1.00 | 1.00 | |||
| Vasculitis | ||||||
| Case | 1.34 (0.94–1.90) | 0.107 | 1.23 (0.76–1.99) | 0.401 | 1.48 (0.65–3.38) | 0.351 |
| Control | 1.00 | 1.00 | 1.00 | |||
| Multivariable |
Matching |
|||||
|---|---|---|---|---|---|---|
| Immunosuppressants not use |
Immunosuppressant use |
Interaction P-value | Immunosuppressants not use |
Immunosuppressant use |
Interaction P-value | |
| OR (95% CI) | OR (95% CI) | OR (95% CI) | OR (95% CI) | |||
| Overall | 1.32 (1.02–1.69) | 1.12 (1.02–1.23) | 0.003 | 1.38 (0.96–2.00) | 1.26 (1.08–1.48) | 0.244 |
| CNS disorders | 1.14 (0.47–2.74) | 1.03 (0.74–1.45) | 0.571 | 1.06 (0.24–4.81) | 0.74 (0.42–1.31) | 0.527 |
| ILD | 1.70 (0.84–3.44) | 1.21 (0.96–1.51) | 0.144 | 1.05 (0.40–2.71) | 1.19 (0.08–1.77) | 0.970 |
| Psoriasis | 1.60 (1.19–2.15) | 1.17 (1.04–1.32) | 0.014 | 1.64 (1.05–2.56) | 1.26 (1.04–1.52) | 0.246 |
| Lupus | 0.52 (0.27–0.99) | 0.89 (0.73–1.08) | 0.271 | 0.63 (0.21–1.91) | 1.04 (0.71–1.53) | 0.597 |
| Vasculitis | <0.01 (0.00–1,000) | 1.03 (0.63–1.68) | 0.970 | <0.01 (0.00–1,000) | 1.27 (0.47–3.40) | 0.968 |
| Simple |
Multivariable |
Matching |
||||
|---|---|---|---|---|---|---|
| OR (95% CI) | P-value | OR (95% CI) | P-value | OR (95% CI) | P-value | |
| Overall | ||||||
| Anti-TNF-α monotherapy | Reference | Reference | Reference | |||
| Anti-TNF-α no use | 0.85 (0.72–1.00) | 0.050 | 0.93 (0.76–1.13) | 0.460 | 0.82 (0.61–1.11) | 0.193 |
| Concomitant use of immunosuppressants | 1.23 (1.03–1.47) | 0.020 | 1.51 (1.22–1.87) | <0.001 | 1.43 (1.04–1.98) | 0.029 |
Values are presented as median (interquartile range), number (%), or mean±standard deviation. BMI, body mass index; IBD, inflammatory bowel disease; TNF, tumor necrosis factor.
Estimates were adjusted for age, sex, IBD subtype, IBD diagnosis year, body mass index, smoking, frequency of hospitalization due to IBD, IBD-related surgeries, Charlson Comorbidity Index, and presence of concomitant immune-mediated inflammatory diseases (e.g., rheumatoid arthritis and ankylosing spondylitis). Estimates were derived from propensity score matched analysis followed by multivariable adjustment using the same covariates. TNF, tumor necrosis factor; IBD, inflammatory bowel disease; OR, odds ratio; CI, confidence interval; CNS, central nervous system; ILD, interstitial lung disease.
Estimates were adjusted for age, sex, IBD subtype, IBD diagnosis year, body mass index, smoking, frequency of hospitalization due to IBD, IBD-related surgeries, Charlson Comorbidity Index, and presence of concomitant immune-mediated inflammatory diseases (e.g., rheumatoid arthritis and ankylosing spondylitis). Estimates were derived from propensity score matched analysis followed by multivariable adjustment using the same covariates. TNF, tumor necrosis factor; OR, odds ratio; CI, confidence interval; CNS, central nervous system; ILD, interstitial lung disease; IBD, inflammatory bowel disease.
Estimates were adjusted for age, sex, IBD subtype, IBD diagnosis year, body mass index, smoking, frequency of hospitalization due to IBD, IBD-related surgeries, Charlson Comorbidity Index, and presence of concomitant immune-mediated inflammatory diseases (e.g., rheumatoid arthritis and ankylosing spondylitis). Estimates were derived from propensity score matched analysis followed by multivariable adjustment using the same covariates. TNF, tumor necrosis factor; OR, odds ratio; CI, confidence interval; IBD, inflammatory bowel disease.
