환자 특성
근치적 절제술을 받은 CRC 환자 153명의 기선 특성은 표 1에 자세히 기술되어 있다. 해당 코호트의 연령 중앙값은 66세(IQR: 17)였으며, BMI 중앙값은 23.1 kg/m2(IQR: 4.7)였다. 현재 흡연자는 환자의 39.9%(n = 61)를 차지했으며, 31.4%(n = 48)는 음주 사실을 보고했다. 종양 병기 분포는 다음과 같았다: 1기 22.2%(n = 34), 2기 47.1%(n = 72), 3기 22.9%(n = 35), 4기 7.8%(n = 12). 입원 기간의 중앙값은 8일(IQR: 3)이었다.
수술 전후 염증 인자 비율
측정된 염증 인자들의 수술 후 대비 수술 전 비율의 중앙값 및 사분위 범위(IQR)는 표 2에 제시되어 있다. 분석 결과, 수술 후 IL–6가 가장 크게 상승하였으며(중앙값 비율 = 7.23, IQR = 28.31), IL–10이 그 뒤를 이었다(중앙값 비율 = 2.13, IQR = 4.18). 반면, IFN–γ(중앙값 비율 = 0.89, IQR = 1.99), IL–12p70(중앙값 비율 = 0.88, IQR = 0.91), IL–17(중앙값 비율 = 0.82, IQR = 0.97), TNF–α(중앙값 비율 = 0.87, IQR = 1.69)의 중앙값 비율은 1.0 미만으로 나타나, 수술 후 감소했음을 나타냈다. 특히 IL–1β, IL–2 및 IL–8 수치는 중앙값 비율이 1.00(IQR은 각각 0.79, 0.95, 1.96)으로 보고되어 상대적으로 안정적으로 유지되었다.
단변량 및 다변량 Cox 회귀 분석 모델
단변량 Cox 비례 위험 회귀 분석 결과, 여러 임상 파라미터가 장기 생존과 유의미한 연관성이 있는 것으로 나타났습니다. 구체적으로 환자 연령 (HR = 1.05, 95% CI = 1.00–1.10, p = 0.031), 종양 단계 (HR = 3.61, 95% CI = 2.18–6.00, p < 0.001), IFN–α (HR = 1.00, 95% CI = 1.00–1.00, p = 0.023) 및 IL–1β 비율 (HR = 1.00, 95% CI = 1.00–1.01, p = 0.012)에서 통계적으로 유의미한 연관성이 확인되었습니다 (p < 0.05). 또한 TNF–α (HR = 1.00, 95% CI = 1.00–1.00, p = 0.076)는 미리 설정된 유의성 임계값 (p < 0.1)을 충족했습니다.
단변량 분석에서 p 값 < 0.1을 나타낸 변수들(연령, 입원 기간, 종양 병기, IFN‑α_ratio, IL‑1β_ratio, TNF‑α_ratio; 총 21건의 이벤트)을 이후 다변량 회귀 모델에 포함시켰습니다. IL‑1β 비율은 두 분석 모두에서 변환 없이 연속 선형 항으로 모델링되었습니다. 잠재적 교란 요인을 보정한 후, 다변량 분석 결과 IL–1β 비율(보정 HR = 1.05, 95% CI = 1.00–1.10, p = 0.037)과 종양 병기(보정 HR = 3.83, 95% CI = 2.27–6.46, p < 0.001) 모두 독립적인 예후 가치를 유지함을 확인하였습니다(유의성 임계값 p < 0.05). IL–1β 비율이 1단위 증가할 때마다 사망 위험이 5% 증가하는 것과 연관이 있었습니다(HR = 1.05, 95% CI = 1.00–1.10, p = 0.037). 단변량 HR인 1.00(95% CI 1.00–1.01)은 소수점 둘째 자리까지 반올림한 결과이며(실제 β ≈ 0.004), 다변량 모델에서 1.05로 변경된 것은 모델 설정의 불일치가 아니라 교란 요인 보정에 기인한 것입니다. 특히, IFN–α (p = 0.023)와 환자 연령(p = 0.031)은 단변량 분석에서 통계적 유의성을 보였으나, 다변량 모델에서는 유의성을 유지하지 못했습니다. 또한, 다른 염증 마커 중 어느 것도 두 분석 모두에서 통계적으로 유의미한 예후 가치를 나타내지 않았습니다. 상세한 통계 결과는 표 3에 제시되어 있습니다.
암 특이 생존율에 미치는 Il-1β 비율의 영향
IL–1β 비율( 그림 2에 표시된 90백분위수 컷오프 값인 5.35를 기준으로 이분화함)에 따라 층화한 환자들의 Kaplan–Meier 분석 결과, 저비율군(<5.35; n = 136)이 고비율군(≥5.35; n = 17)에 비해 지속적인 생존 이점을 보였으며(log–rank p = 0.026), 생존 곡선은 점진적으로 분리되었습니다. 이 수치는 40개월 추적 관찰 시 각각 121명과 12명으로 감소하였습니다(저비율군에서 16명, 고비율군에서 5명의 사망에 해당). 추적 관찰 기간 동안 총 21건의 암 관련 사망이 발생했으며, 연구 기준에 따라 암과 관련 없는 사망은 제외하고 종점으로 집계하지 않았습니다. 병기와 관련된 잠재적 교란 요인을 추가로 다루기 위해, 본 연구에서는 1~3기 환자(n = 141, 사망 16명)로 제한하여 사후 암 특이 생존 분석을 추가로 수행하였습니다. 이 하위 그룹에서 IL–1β 비율에 대한 보정 HR은 1.04(95% CI 0.99–1.10, p = 0.11)였으며, 이는 일관된 경향을 보였으나 통계적으로 유의하지 않았는데, 이는 이벤트 발생 수가 적어 통계적 검정력이 감소했기 때문으로 보입니다.
데이터 가용성:
본 연구의 결과를 뒷받침하는 데이터는 교신저자에게 요청하여 이용할 수 있습니다. 비식별화된 데이터는 부록 표 1에서, 소스 코드는 부록 파일 1에서 확인할 수 있습니다.

그림 1: 환자 선정 및 데이터 분석 흐름도. 근치적 절제술을 받은 대장암(CRC) 환자 총 429명을 1차로 선별하였다. 276명의 환자(추적 관찰이 불완전한 193명(암과 관련 없는 사망 11명 포함), 수술 전 또는 수술 후 IL-1β 데이터가 없는 83명)를 제외한 후, 최종 분석에 153명의 환자를 포함하였다. 이 그림은 Microsoft PowerPoint (Microsoft 365)로 제작되었다. Copyright 2025 The Author(s). 이는 저자들이 직접 작성한 원본 그림이다. All rights reserved. 약어: IL-1β = interleukin‑1 beta. 여기를 클릭하여 이 그림의 더 큰 버전을 확인하십시오.

그림 2: 수술 후 대 수술 전 IL-1β 비율에 따라 층화된 CRC 환자(n = 153)의 암 특이적 생존율에 대한 Kaplan-Meier 분석. 비율은 90 백분위수 컷오프인 5.35를 기준으로 이분화되었음: 고비율군(≥5.35, n = 17) 대 저비율군(<5.35, n = 136). 로그 순위법(log-rank test) 결과 p = 0.026으로, 저비율군에서 유의미한 생존 이점이 있음이 나타남. 음영 구역은 95% 신뢰 구간(CI)을 나타냄. 약어: IL-1β = interleukin-1 beta; CRC = colorectal cancer. 이 그림의 더 큰 버전을 보시려면 여기를 클릭하십시오.
| 요인 | n = 153 |
| 연령, 연도 | 66 (17) |
| 체질량지수, kg/m²2 | 23.1 (4.7) |
| 흡연 | 61 (39.9%) |
| 음용 | 48 (31.4%) |
| 종양 병기 | |
| I | 34 (22.2%) |
| II | 72 (47.1%) |
| III | 35 (22.9%) |
| 정맥 내 투여 | 12 (7.8%) |
| 입원 기간, 일수 | 8 (3) |
표 1: 전체 코호트의 기본 정보. 변수는 중앙값 (IQR = Q3–Q1) 또는 n (%)로 표시됨. 약어: BMI = 체질량 지수; IQR = 사분위수 범위.
| 변수 | 중앙값 | 사분위수 범위 |
| IFN–α 비율 | 1.02 | 1.81 |
| IFN-γ 비율 | 0.89 | 1.99 |
| IL-10 비율 | 2.13 | 4.18 |
| IL–12p70_비율 | 0.88 | 0.91 |
| IL-17 비율 | 0.82 | 0.97 |
| IL-1β 비율 | 1 | 0.79 |
| IL-2 비율 | 1 | 0.95 |
| IL-4 비율 | 0.85 | 2.23 |
| IL–5_비율 | 1.06 | 1.17 |
| IL-6 비율 | 7.23 | 28.31 |
| IL-8 비율 | 1 | 1.96 |
| TNF–α_비율 | 0.87 | 1.69 |
표 2: 수술 전후 염증 인자 비율.참고: IQR은 Q3와 Q1의 차이를 나타냄. 약어: IQR = 사분위수 범위; IFN–β = 인터페론-베타; IFN–γ = 인터페론-감마; IL–1β = 인터루킨-1 베타; TNF–α = 종양 괴사 인자-알파; IL–2_ratio = 수술 후 대비 수술 전 IL–2 비율.
| 변수 | 단변량 분석 | 다변량 분석 |
| HR (95% CI) | p 값 | HR (95% CI) | p 값 |
| 연령 | 1.05 (1.00–1.10) | 0.031 | 1.04 (0.99–1.09) | 0.09 |
| BMI | 0.95 (0.83–1.09) | 0.471 | | |
| 흡연 | 0.92 (0.38–2.23) | 0.858 | | |
| 음주 | 1.12 (0.45–2.78) | 0.805 | | |
| 종양 단계 | 3.61 (2.18–6.00) | < 0.001 | 3.83 (2.27–6.46) | < 0.001 |
| 입원 기간 | 1.03 (1.00–1.06) | 0.079 | 1.03 (0.99–1.07) | 0.151 |
| IFN–α_ratio | 1.00 (1.00–1.00) | 0.023 | 0.98 (0.96–1.01) | 0.254 |
| IFN–γ_ratio | 1.00 (0.98–1.02) | 0.981 | | |
| IL–10_ratio | 0.99 (0.96–1.02) | 0.524 | | |
| IL–12p70_ratio | 0.97 (0.83–1.13) | 0.683 | | |
| IL–17_ratio | 0.92 (0.77–1.11) | 0.385 | | |
| IL–1β_ratio | 1.00 (1.00–1.01) | 0.012 | 1.05 (1.00–1.10) | 0.037 |
| IL–2_ratio | 0.91 (0.72–1.16) | 0.451 | | |
| IL–4_ratio | 1.00 (0.99–1.01) | 0.955 | | |
| IL–5_ratio | 0.98 (0.92–1.04) | 0.524 | | |
| IL–6_ratio | 1.00 (1.00–1.00) | 0.614 | | |
| IL–8_ratio | 0.99 (0.94–1.04) | 0.593 | | |
| TNF–α_ratio | 1.00 (1.00–1.00) | 0.076 | 1.00 (0.99–1.00) | 0.323 |
표 3: 암 특이적 생존 값에 대한 Cox 회귀 분석. 약어: HR = 위험비; 95% CI, 95% 신뢰 구간. 굵게 표시된 값은 통계적으로 유의하거나(p < 0.1) 한계적으로 유의함(p < 0.05). 다변량 모델에는 연령, 입원 기간, 종양 단계, IFN‑α_ratio, IL‑1β_ratio 및 TNF‑α_ratio의 6가지 공변량이 포함되었다. IL‑1β 비율은 로그 변환 없이 연속 선형 항으로 입력되었다.
보충 표 1: 원시 데이터. 연구에 사용된 비식별화된 데이터 세트입니다.이 파일을 다운로드하려면 여기를 클릭하십시오.
보충 파일 1: 소스 코드. 본 연구에 사용된 소스 코드입니다. 이 파일을 다운로드하려면 여기를 클릭하십시오.