Genomic Copy Number Variation Profiling of Pleural Effusion cfDNA via LP-WGS for Precision Diagnosis and Prognostic Stratification in Lung Adenocarcinoma-associated Malignant Pleural Effusion running title: cfDNA CNV for LUAD-MPE Diagnosis and Prognosis

기술 요약: LP-WGS를 통한 흉수 cfDNA의 게놈 복제수 변이 프로파일링 연구 (LUAD-MPE 대상)

문제 제기
폐선암(LUAD) 관련 악성 흉수(MPE)는 진행된 질환의 특징적인 징후로, 중앙 생존 기간은 3~12개월에 불과하다. 현재의 진단 전략은 다음과 같은 결정적인 한계에 직면해 있다: 세포학적 검사는 악성 세포 수가 적고 분포가 불균일하여 변동성이 큰 민감도(40–87%)를 보인다; 침습적 생검(예: 흉강경 검사)은 항상 허용될 수 있는 것은 아니다; 또한, 전통적인 바이오마커인 암태아성 항원(CEA)은 충분한 정확성을 갖추지 못했다. 더욱이, MPE에서 임상적 위험 계층화를 위한 검증된 예후 바이오마커가 부족한 실정이다. 조기 진단, 동적 모니터링 및 개인 맞춤형 관리를 용이하게 하기 위한 비침습적이고 고감도의 도구가 절실히 필요하다.

연구 방법론
본 연구는 난창 대학교 제1부속 병원에서 치료받은 41명의 환자(확진된 LUAD 관련 MPE 환자 30명 및 양성 흉수[BPE] 환자 11명)를 대상으로 한 전향적 코호트 설계를 채택하였다.

• 시료 처리: 흉수(50 mL)를 수집, 원심 분리한 후 상층액으로부터 세포 유리 DNA(cfDNA)를 분리하였다.
• 시퀀싱: cfDNA 라이브러리에 대해 저패스 전게놈 시퀀싱(LP-WGS)을 수행하였다. 평균 시퀀싱 깊이는 2.4×(범위 1.5×–4.3×)였으며, HiSeq X10 플랫폼을 사용하여 150-bp 페어드 엔드(paired-end) 리드를 활용하였다.
• 생물정보학 분석: 맞춤형 워크플로우인 초민감 염색체 이배성 검출기(UCAD)를 사용하여 복제수 변이(CNV) 프로파일링을 수행하였다. 이는 단일 염기 시퀀싱 깊이를 계산하고, 200-kb 게놈 빈(bin)에 대해 Z-score로 커버리지를 정규화하며, 연속적인 게놈 세그먼트를 식별하기 위해 순환 이진 분할(circular binary segmentation, DNAcopy 패키지)을 적용하는 과정을 포함한다.
• 지표: 은닉 마르코프 모델(hidden Markov model)을 사용하여 종양 분율(TF)을 계산하였고, 게놈 전반의 |Z score|를 게놈 불안정성의 대리 지표로 평가하였다.
• 통계 분석: 진단 성능는 ROC 곡선 및 AUC를 통해 평가되었다. 게놈 특징과 세포학적 종양 세포 비율 및 분화 상태를 포함한 임상병리학적 변수 간의 상관관계는 Spearman 상관계수를 사용하여 분석하였다. 생존 분석은 Kaplan–Meier 곡선 및 log-rank 테스트를 사용하여 수행되었다.

주요 기여 및 결과

1. 진단적 우월성: MPE 환자군에서는 CNV가 96.7%에서 검출된 반면, 양성 대조군에서는 9.1%에 불과했다 (P < 0.001). ≥ 3.835의 컷오프 값을 가진 게놈 전반의 |Z score|를 활용한 통합 CNV 기반 진단 모델은 0.939(95% CI 0.855–0.999)의 곡선 아래 면적(AUC)을 달성하였으며, 96.7%의 민감도와 90.9%의 특이도를 나타냈다. 이 모델은 전통적인 CEA 검사(AUC = 0.888)보다 우수한 성능을 보였으며, CEA 음성인 4건의 MPE 사례를 성공적으로 식별하였다.
2. 뚜렷한 게놈 경관: MPE 샘플은 양성 샘플에서는 나타나지 않는 재발성 국소 증폭 및 광범위한 결실을 보였다. 주목할 만한 증폭 부위는 1q, 5p, 7p, 8q, 17q였으며(온코진인 MCL1, MDM4, EGFR, MYC, ERBB2 포함), 유의미한 결실은 3p, 6q, 8p, 9p, 9q, 10q, 11p, 12q, 13q, 15q, 17p, 18q에서 관찰되었다(종양 억제 유전자인 CDKN2A, PTCH1, PTEN, RB1, TP53 포함). 본 연구는 동기적 팔 수준 불균형(예: 1p 결실/1q 증폭 및 17p 결실/17q 증폭)과 특정 염색체 사건 간의 상관관계(예: 7p와 17q 증폭 사이의 양의 상관관계)를 확인하였다.
3. 종양 부담 및 공격성과의 상관관계: 종양 분율(TF)과 게놈 전반의 |Z score|는 세포학적 종양 세포 비율과 강한 양의 상관관계를 보였다 (각각 R = 0.67 및 R = 0.70; P < 0.0001). 저분화 종양은 중등도 분화 종양에 비해 유의미하게 높은 게놈 불안정성(높은 중앙값 |Z score|)을 나타냈다 (P = 0.032).
4. 예후 계층화: 추적 관찰 데이터가 있는 24명의 MPE 환자를 대상으로 한 생존 분석 결과, 특정 CNV 변화가 짧은 전체 생존 기간(OS)과 연관되어 있음이 밝혀졌다. 18p 증폭과 8p, 9q, 19q 결실은 불량한 예후와 유의미한 상관관계를 보였다 (P < 0.05).

의의 및 주장
저자들은 LP-WGS 기반의 흉수 cfDNA CNV 분석이 세포학의 민감도 한계와 CEA의 정확도 제약을 해결할 수 있는 비용 효율적인(샘플당 약 $200) 비침습적 도구라고 주장한다. 22개 염색체에 걸친 전역적 CNV 신호를 합성함으로써, 이 방법은 CEA 음성 사례에서도 MPE를 양성 흉수와 구별할 수 있는 강력한 프로파일을 제공한다.

본 연구는 이 접근법을 세포학 및 CEA를 보완하는 검사로서, 특히 침습적 생검을 수행할 수 없는 환자들에게 유용하다고 위치시킨다. 또한, 특정 CNV 패턴(18p 증폭, 8p/9q/19q 결실)과 TF 및 |Z score|와 같은 지표가 종양 부담, 분화 상태를 반영하고 생존 결과를 예측할 수 있다는 점에서 예후 계층화의 효용성을 주장한다. 저자들은 이 방법이 진행된 LUAD 관련 MPE의 조기 진단 및 개인 맞춤형 관리를 촉진한다고 결론짓지만, 일반화 가능성을 확인하기 위해 다기관 검증이 필요함을 인정하였다.