Stable but confounded: brain region, not diagnosis, drives a validation-passing molecular subtype in postmortem psychiatric brain

기술 요약: "안정적이지만 혼란스러운: 사후 정신의학 뇌 연구에서 검증을 통과하는 분자적 아형을 결정하는 것은 진단명이 아닌 뇌 영역이다"

문제 제기
뇌 전사체(transcriptome)의 비지도 분자 아형 분류(unsupervised molecular subtyping)는 정신질환(예: 자폐증, 조현병)의 이질성을 해결하기 위해 점점 더 많이 사용되고 있다. 그러나 클러스터링 알고리ست은 해당 파티션이 재현 가능한 생물학적 구조를 반영하는지, 아니면 단순히 기술적/해부학적 인공물(artifact)인지와 관계없이 항상 분할을 반환한다. 본 논문은 사후 정신의학 연구에서 세 가지 주요 위험 요소를 식별하였다:

1. 작은 표본 크기 (수십 개에서 수백 개 수준): 이는 통계적 검정력을 제한한다.
2. 피험자 중복: 명목상 독립적인 코호트 간의 피험자 중복은 허위 복제 주장을 초래한다.
3. 해부학적 지배력: 뇌 영역에 따른 전사체 변이가 진단명에 따른 변이보다 훨씬 크다. 결과적으로, 여러 영역을 통합한 데이터를 클러스터링하면 질병 아형이 아닌 해부학적 구조가 회복된다.

저자들은 현재의 검증 관행이 (특히 뇌 영역과 같은 혼란 변수를 명시적으로 테스트하지 않을 경우) "안정적인" 클러스터(재현 가능한 구조)와 "유효한" 클로스터(질병과 생물학적으로 관련이 있는 구조)를 구분하는 데 실패하고 있다고 주장한다.

방법론
저자들은 확립된 대조군을 구성하고 명시적인 혼란 변수 테스트를 추가하는 "검증 우선(validation-first)" 원칙을 사용하여 3개의 공개 사후 코호트(총 384개 샘플)를 재분석하였다. 데이터셋은 다음과 같다:

• 자폐 스펙트럼 장애 (ASD): GSE28521 (마이크로어레이, 전두엽/측두엽 피질/소뇌) 및 GSE64018 (RNA-seq, 상측두회).
• 조현병, 양극성 장애, 주요 우울 장애: GSE80655 (RNA-seq, 전대상피질, DLPFC, 측좌핵).

분석 파이프라인:

1. 특징 감소 (Feature Reduction): 유전자 발현을 단일 샘플 농축 분석(ssGSEA 또는 유전자 평균 z-score)을 통해 경로 수준의 점수로 요약하였다.
2. 클러스터링: 베이지안 정보 기준(BIC)을 통한 성분 선택을 포함한 가우시안 혼합 모델(GMM)을 사용하여 샘플을 클러스터링하였다.
3. 검증 통제 (Validation Controls): 모든 후보 파티션은 네 가지 게이트를 거쳤다:
   • 라벨 순열 널(Label-permutation null): 순열된 라벨로부터 클러스터가 회복 가능한지 테스트.
   • 무작위 유전자 세트 널(Random-gene-set null): 파티션이 큐레이션된 경로에 특이적인지, 아니면 무작위 유전자 세트에서 발생하는지 테스트.
   • 부트스트랩 안정성(Bootstrap stability): 원래 할당과 재표집된 할당 사이의 조정 랜드 지수(ARI) 측정 (임계값 $\ge$ 0.80).
   • 계통적 독립성(Ancestry independence): 계통 주성분과의 연관성 확인.
4. 혼란 테스트 (Confound Testing): 파티션과 뇌 영역, 기증자 중복, 임상 진단 사이의 명시적인 분할표 연관성($\chi^2$, Cramér's V) 분석을 수행하였다.
5. 비순환적 분석 (Non-Circular Analysis): "이중 계산(double dipping)"을 피하기 위해, 질병 신호는 클러스터링에서 유도된 특징이 아닌, 고정된 정규 세포 유형 마커 유전자(예: PVALB, SST, GAD1) 패널을 사용하는 진단 기반 케이스-대조군 대조를 통해 평가되었다.

주요 결과

1. 자폐증 "GABA성" 아형의 실패:
   GABA성 경로 하향 조절로 특징지어지는 가설적 ASD 하위 그룹은 부트스트랩 안정성 검증(ARI = 0.71, 임계값 0.80 미달)에서 실패하였다. 또한, 두 번째 코호트(GSE64018)에서의 외견상 복제는 두 코호트가 24명의 피험자 중 10명을 공유했기 때문에 무효화되었다. 즉, 이는 독립적인 복제가 아니라 동일한 뇌에 대한 플랫폼 간의 상호 확인이었다.

2. "검증을 통과한" 조현병 파티션은 영역 분류기이다:
   조현병 코호트(GSE80655)에서 $k=3$ 파티션은 모든 통계적 게이트를 통과하였다 (부트스트랩 ARI = 0.923; 무작위 유전자 세트 및 계통 체크 통과). 그러나 혼란 테스트 결과, 이 파티션은 진단이 아닌 뇌 영역에 의해 완전히 구동됨이 드러났다:

   • 아형 0과 2는 주로 측좌핵(nAcc)이었다.
   • 아형 1은 배측 전대상피질(AnCg)과 DLPFC를 포함한 피질 영역에 국한되었다.
   • 파티션과 영역 사이의 연관성은 압도적이었으나 ($p = 4.3 \times 10^{-26}$, Cramér's V = 0.67), 진단과의 연관성은 유의하지 않았다 ($p = 0.41$).
   • 분석을 단일 영역(DLPFC)으로 제한했을 때 안정성이 검증 임계값 아래로 떨어졌으며 (ARI = 0.68), 이는 "안정성"이 영역 축에 의한 인공물임을 확인시켜 주었다.

3. 인공물의 체계적 성격:
   이 영역 중심의 인공물은 코호트 내의 모든 진단명(양극성 장애, 주요 우울 장애)에 걸쳐 체계적이었으며, 클러스터의 수($k$)와 무관하게 나타났다. 네 가지 진단을 통합한 코호트는 가장 강한 영역 구조를 보였으나 가장 낮은 안정성을 보였다.

4. 비순환적 분석을 통한 알려진 생물학적 기전 회복:
   클러스터링 프레임워크는 질병 신호를 놓쳤지만, 전통적인 비순환적 케이스-대조군 분석은 확립된 인터뉴런 생물학을 회복하였다:

   • ASD: 피질에서의 파발부민(PVALB) 마커의 유의미한 결핍.
   • 조현병: 피질 영역에서의 소마토스타틴(SST) 및 GABA 합성 효소(GAD1/GAD2)의 유의미한 결핍.
   • 이중 해리(Double Dissociation): 유의미한 반대 방향의 상호작용을 통해 PVALB 결핍은 ASD에서 더 두드러지고, SST 결핍은 조현병에서 더 두드러진다는 것이 확인되었다.

5. 벤치마크 실패:
   적응형 임계값 모델(Adaptive-threshold model, $R^2 = 0.889$라고 주장함)을 동반한 캘리브레이션 벤치마크는 재현에 실패하였다. 저자들은 이를 메트릭 결함으로 추적하였다: 조정 랜드 지수(ARI)가 퇴행적 입력(단일 진정한 클러스터를 가진 데이터셋)에 대해 허위로 1.0의 값을 반환하여 상관관계를 인위적으로 부풀렸다.

의의 및 주장
본 논문은 새로운 분자적 아형을 발견한다고 주장하는 것이 아니라, 엄격한 통계적 안정성 검사를 통과하면서도 알려진 혼란 변수(뇌 영역)에 의해 완전히 구동될 수 있는 "실제 사례 연구"를 제시한다.

저자들은 다음과 같이 결론짓는다:

• 안정성 $\neq$ 유효성: 높은 부트스트랩 안정성은 재현성을 보증하지만, 그것이 의도한 생물학적 신호(질병 vs 해부학)를 반영한다는 것을 보장하지는 않는다.
• 혼란 테스트는 필수적이다: 여러 영역이나 코호트를 통합하는 모든 아형 분석에는 파티션을 뇌 영역 및 기타 혼란 변수와 대조하는 작업이 루틴한 검증 구성 요소가 되어야 한다.
• 발견보다 규율(Discipline over Discovery): 작은 규모의 통합된 다영역 코호트에서 탐색적 발견은 실패하기 쉽다. 저자들은 특정 체크리스트를 옹호한다: 널 컨트롤(null controls) 보고, 혼란 변수에 대한 테스트, 피험자가 서로 배타적인 복제 보장, 순환적 추론 회피, 그리고 컨트롤을 통과하지 못한 파티션을 공개하는 것.

이 연구는 분자적 아형 분류가 성공할 수 있음을 언급하면서도(Bowen et al., 2019의 영역 제어 성공 사례 참조), 해부학적 지배력과 검증되지 않은 파이프라인 기본 설정에 의해 발생하는 구체적인 실패 모드—현재의 정신의학 유전체학에서 중요한 함정—를 강조한다.