Coagulopathy subphenotypes and mortality in critically ill patients with sepsis: a multicentre retrospective cohort study using MIMIC-IV and SICdb

기술 요약: 중증 패혈증 환자의 응고병증 하위 표현형 및 사망률

문제 제기
패혈증 유발 응고병증은 높은 사망률과 연관된 흔하고 이질적인 상태이다. Sepsis-Induced Coagulopathy (SIC) 점수와 같은 기존의 점수 체계는 임계값 기반의 서수 척도를 제공하지만, 다차원적인 응고 장애를 단일 지표로 압축함으로써 질적으로 다른 형태의 지혈 기능 부전을 구분해내지 못한다. 이러한 이질성은 왜 대규모 무작위 대조 시험에서 항응고제 투여 효과가 일관되게 부정적으로 나타나는지를 설명할 수 있는데, 이는 광범위한 환자 등록이 특정 하위 집단에 국한된 치료 효과를 희석시켰을 가능성을 시사한다. 기존의 하위 표현형 연구들은 단일 국가 데이터베이스에 국한되었으며, 이상치(outlier)나 분포 가정에 민감한 클러스터링 방법(예: k-means, 잠재 클래스 분석)에 의존하였고 견고한 외부 검증이 부족했다.

방법론
저자들은 두 개의 구별되는 데이터셋을 활용하여 다기관 회고적 코호트 연구를 수행하였다:

1. 도출 코호트 (Derivation Cohort): 완전한 기초 응고 프로파일을 보유한 성인 패혈증 환자 5,115명을 포함하는 MIMIC-IV v2.2 (미국 Beth Israel Deaconess Medical Center).
2. 외부 검증 코호트 (External Validation Cohort): 외과 중환자실 환경의 성인 패혈증 환자 622명을 포함하는 SICdb v1.0.6 (오스트리아 University Hospital Salzburg).

데이터 처리 및 클러스터링:

• 변수: 네 가지의 일상적으로 측정되는 응고 변수를 조화시켰다: 국제 표준화 비율(INR), 활성화 부분 트롬보플라스틴 시간(aPTT), 피브리노겐, 혈소판 수. D-dimer와 트롬빈 시간은 결측치가 많아 제외되었다.
• 전처리: 변수들을 $\log_2$ 변환 후 z-표준화하였다.
• 알리즘: 80% 아이템 서브샘플링을 적용한 500회의 부트스트랩 재표본 추출을 사용하여 컨센서스 클러스터링(Consensus clustering)을 적용하였다. 평균 연결법(average linkage)과 $1 - \text{Pearson correlation}$ 거리 지표를 사용한 계층적 클러스터링이 적용되었다.
• 클러스터 선택: 최적의 클러스터 수($k$)는 비중적 클러스터링 비율(Proportion of Ambiguous Clustering, PAC), 평균 실루엣 폭, 클로스터 안정성 및 생물학적 해석 가능성을 평가하여 결정되었다.
• 검증: SICdb 코호트의 하위 표현형은 MIMIC-IV 코호트에서 도출된 중심점(centroid)을 기준으로 최근접 중심점 할당(nearest-centroid assignment) 방식을 통해 할당되었으며, 분포 변화를 고려하여 검증 샘플 공간 내에서 표준화가 적용되었다.

통계 분석:
일차 결과 변수는 28일 전인 사망률(all-cause mortality)이다. 생존 차이는 Kaplan–Meier 곡선과 log-rank 테스트를 사용하여 평가하였다. 다변량 Cox 비례 위험 모델은 다음 세 단계의 순차적 보정을 거쳐 적합되었다:

1. 미보정(Unadjusted).
2. 인구통계학적 특성(연령, 성별) 및 동반 질환(Charlson index) 보정.
3. 질병 중증도(비응고 관련 SOFA, APACHE III, 활력 징후 등) 보정.
   참고: 과잉 보정 편향을 방지하기 위해 SOFA 점수의 응고 관련 구성 요소는 조정에서 명시적으로 제외되었다.

주요 결과
컨센서스 클러스터링을 통해 세 가지의 재현 가능한 하위 표현형이 식별되었다:

1. 표현형 1 (응고 보존형, Coagulation-preserved): MIMIC-IV 환자의 56.8%. 정상에 가까운 헤모스타시스(hemostasis)가 특징이다.
2. 표현형 2 (고립된 aPTT 연장형, Isolated aPTT prolongation): MIMIC-IV 환자의 10.0%. 다른 파라미터의 경미한 이상만 동반된 채 뚜렷한 aPTT 연장을 보인다.
3. 표현형 3 (범응고병증, Pan-coagulopathy): MIMIC-IV 환자의 33.2%. 네 가지 축 모두에서 전반적인 이상(상승된 INR/aPTT, 낮은 피브리노겐/혈소판)을 보인다.

사망률 연관성:

• 도출 코호트 (MIMIC-IV):
  • 표현형 3은 표현형 1과 비교하여 일관되고 중증도에 독립적인 초과 사망 위험을 보였다 (28일 HR 1.25, 95% CI 1.04–1.51).
  • 표현형 2는 예후적으로 중립적이었다 (28일 HR 1.01, 95% CI 0.75–1.35).
• 검증 코호트 (SICdb):
  • 표현형 분포가 유의하게 달랐다 (표현형 3은 외과 환자의 높은 기저 피브리노겐 수치로 인해 6.6%로 드물게 나타남).
  • 표현형 3은 사망률 증가 경향을 보였으나 (HR 1.55), 작은 표본 크기($n=41$)로 인해 통계적 유의성에 도달하지 못했다.
  • 표현형 2는 극명하게 다른 프로파일을 보였다: 이는 미국 코호트에서의 중립성과 대조적으로, 두 배의 조정 사망 위험을 나타냈다 (28일 HR 2.01, 95% CI 1.34–2.99).
• SIC 점수와의 상호작용:
  • 하위 표현형은 보완적인 층화 기능을 제공했다. SIC 양성 환자들 사이에서 표현형 2는 사망률 증가와 연관되었으나 (HR 1.60), SIC 음성 환자들에게서는 그러한 연관성이 존재하지 않았다.
  • 표현형 3은 SIC 양성 환자의 74.6%를 차지하였으며, 이는 SIC 양성 판정이 본질적으로 범응고병증을 농축함을 시사한다.

의의 및 주장
본 논문은 종양학 분야의 견고한 클래스 발견 표준인 컨센서스 클러스터링을 패혈증 응고병증 하위 분류에 적용한 첫 사례라고 주장한다. 주요 기여는 다음과 같다:

1. 국가 간 재현성: 서로 다른 케어 구조(내과 vs 외과 중환자실)를 가진 두 국가의 데이터베이스를 통해 세 가지의 생물학적으로 해석 가능한 하위 표현형이 식별되고 검증될 수 있음을 입증하였다.
2. 맥락 의존적 예후: 응고 패턴의 예후적 의미가 보편적이지 않음을 확립하였다. 구체적으로, 고립된 aPTT 연장(표현형 2)은 외과 중환자실 환경 및 SIC 양성 환자군에서는 고위험 마커이지만, 일반 내과 중환자군에서는 중립적이다.
3. 위험 층화의 정교화: 네 가지 변수의 응고 프로파일을 SIC 점수와 통합하는 것이 단순하고 비용이 들지 않는 방법으로 위험 층화를 정교화할 수 있음을 주장한다. 이를 통해 스칼라 점수만으로는 놓칠 수 있는 고위험 구성(예: 범응고병증 또는 SIC 양성 고립 aPTT)을 식별할 수 있다.
4. 임상 시험 설계에 대한 시사점: 향후 항응고제 시험이 표현형 농축 설계를 고려해야 함을 시사한다. 범응고병증 표현형은 농축의 자연스러운 후보이며, 고립된 aPTT 표현형은 현재 표준 시험 선정 기준에는 보이지 않지만 특정 치료에 반응할 가능성이 있는 하위 집단을 식별해낸다.

저자들은 관찰 연구의 특성상 인과적 추론에 대해 겸허한 태도를 유지하며, 완전한 응고 패널을 필요로 하는 선택 편향(selection bias)과 가용 데이터베이스 내 특수 검사(예: D-dimer)의 부재를 인정하였다. 그들은 표현형 2에 관한 결과가 가설 생성적(hypothesis-generating)임을 강조하며, 다양한 환경에서의 추가 검증 필요성을 역설하였다.