Donor Aware Reanalysis of Neuronal Nuclei Supports STMN2 Depletion Associated with Nuclear TDP 43 loss in FTD and ALS

기술 요약: 신경세포 핵의 기증자 인지 재분석을 통한 FTD 및 ALS에서 핵 TDP-43 소실과 관련된 STMN2 결핍 지지

문제 제기
TAR DNA 결합 단백질 43(TDP-43)의 핵 내 결핍은 근위축성 측삭 경화증(ALS) 및 전두측엽 치매(FTD)의 특징이며, 이는 특히 암호적 엑손(cryptic exons)의 억제와 관련된 RNA 가공 기능의 상실과 연결되어 있습니다. 그러나 소규모 사후 검체 코호트에서의 전사체 연관성은 기증자 간 변이에 의해 종종 가려집니다. 더욱이, 유전자 수준의 발현 카운트와 특정 RNA 가공 이벤트(예: 암호적 엑손 포함, 넌센스 매개 분해) 사이에는 결정적인 차이가 존재합니다. 유전자 전체 카운트만으로는 변화된 발현이 전사 변화, RNA 안정성 또는 특정 스플라이싱 이벤트에서 기인한 것인지 구분할 수 없습니다. 본 연구는 공공 RNA-seq 데이터를 재분석하여 핵 내 TDP-43 소실과 관련된 유전자 발현 변화를 정량화하는 동시에, 기증자 구조를 명시적으로 모델링하고 유전자 수준 추론의 한계를 규명함으로써 이 문제를 다룹니다.

방법론
본 연구는 Liu 등이 생성한 공공 데이터셋 GSE126542(BioProject PRJNA52295)를 카운트 수준에서 2차 재분석을 수행했습니다. 이 데이터셋은 전두엽 피질의 NeuN 양성 신경세포 핵에서 유래한 7명의 기증자로부터 얻은 12개의 라이브러리(TDP-43 유지 핵과 소실 핵이 각각 5쌍의 완전한 매칭 쌍 및 2개의 싱글톤 기증자로 구성)로 구성됩니다.

• 데이터 검증: 저자들은 제공된 카운트 행렬을 featureCounts 요약 통계와 대조하여 검증하였으며, 62,754개의 유전어와 0 이상의 정수 카운트를 가진 12개의 라이브러리를 확인했습니다. 모든 라이브러리에서 총 카운트가 10 이상인 유전자에 대해 필터링을 수행한 후, 44,225개의 유전자가 분석에 유지되었습니다.
• 통계 모델링: 차등 발현은 음수 이항 프레임워크 하에서 PyDESeq2 (v0.5.4)를 사용하여 추정되었습니다.
  • 주요 모델: ~ donor + nuclear_status. 이 "기증자 인지(donor-aware)" 모델은 기저의 개별 간 차이를 고려하며, 5개의 완전한 매칭 쌍을 활용하여 핵 내 TDP-43 상태에 따른 효과를 식별합니다.
  • 민감도 모델: 모델 의존성을 평가하기 위해 상태 전용 모델(~ nuclear_status)과 진단 조정 모델(~ diagnosis + nuclear_status)을 적합시켰습니다.
• 대조 정의: 대조는 "TDP-43 소실"에서 "TDP-43 유지"를 뺀 것으로 정의되었습니다. 음수의 log2 fold change (LFC)는 소실 핵에서 낮은 발현을 나타냅니다.
• 분석 범위: 분석은 엄격하게 유전자 수준의 카운트로 제한되었습니다. 저자들은 가용한 데이터(유전자 수준 카운트)로는 암호적 엑손 포함, percent-spliced-in (PSI) 값, 이소폼 전환 또는 단백질 풍부도를 해결할 수 없음을 명시하며, 이러한 것들을 계산하지 않았습니다.

주요 결과

1. 전역적 전사체 변이: 기증자 인지 모델은 FDR < 0.05인 3,442개의 유전자를 식별했습니다(1,865개 감소 및 1,577개 증가). 이는 상태 전용(1,295개) 또는 진단 조정(1,495개) 모델보다 더 큰 집합으로, 이 소규모 코호트에서 기증자 구조를 고려하는 것이 통계적 검정력을 높일 수 있음을 시사하지만, 발견되는 유전자의 수는 모델에 따라 달라집니다. 모델 간의 효과 크기는 강한 상관관계를 보였습니다 (Pearson r = 0.831).
2. STMN2 결핍: STMN2(stathmin-2 인코딩)는 TDP-43 소실 핵에서 견고하고 일관된 결핍을 보였습니다.
   • 기증자 인지 모델: LFC = −1.813 (95% Wald 구간: −2.229 ~ −1.396; FDR = 6.22 × 10⁻¹⁵)로, 약 71.5% 낮은 발현에 해당합니다.
   • 일관성: 5개의 완전한 기증자 쌍 모두에서 정규화된 STMN2 카운트가 TDP-43 소실 라이브러리에서 더 낮았습니다. 이 연관성은 세 가지 민감도 모델 모두에서 강력하게 유지되었습니다.
3. TARDBP 발현: TARDBP(TDP-43 인코딩)는 소실 핵에서 중간 정도로 낮게 나타났으며 (LFC = −0.606; FDR = 4.98 × 10⁻⁶), 이는 핵 내 소실이라는 출처 주석과 일치합니다.
4. UNC13A 발현: 단순한 결핍 가설과는 반대로, 유전자 전체 UNC13A 발현은 기증자 인지 모델에서 약간 더 높았습니다 (LFC = 0.372; FDR = 0.00622), 이는 약 29.4% 높은 발현에 해당합니다. 이 결과는 민감도 모델에서는 유의하지 않았습니다. 5개의 완전한 쌍 모두에서 소실 라이브러리의 정규화된 UNC13A 카운트가 더 높았으나, 한 쌍에서는 증가 폭이 작았습니다.

의의 및 주장
논문은 이 재분석을 바탕으로 세 가지 주요 결론을 주장합니다:

1. 모델 의존성: 핵 내 TDP-43 소실과 연관된 유의한 유전자의 수는 통계적 설계에 크게 의존합니다. 기증자 인지 모델링은 이 소규모의 이질적인 코호트에서 상태 전용 모델보다 더 많은 연관성을 드러내지만, 이러한 발견은 독립적인 재현이라기보다 모델 의존적인 추론으로 간가되어야 합니다.
2. STMN2의 견고성: STMN2의 결핍은 코호트 내에서 견고한 연관성입니다. 이는 세 가지 통계 모델과 다섯 개의 매칭 기증자 쌍 모두에서 재현되었으며, 핵 내 TDP-43 소실과 stathmin-2 출력 감소 사이의 생물학적 연결을 뒷받침합니다. 그러나 이 데이터는 접합부(junction)-수준 또는 단백질 수준의 검증 없이는 구체적인 메커니즘(예: 조기 폴리아데닐화)을 증명하지 못합니다.
3. UNC13A에 대한 유전자 수준 카운트의 한계: 본 연구는 유전자 전체 UNC13A 발현을 기능적 단백질 결핍이나 암호적 엑손 포함의 대리 지표로 사용하는 것에 대해 명시적으로 경고합니다. 기증자 인지 모델에서 관찰된 유전자 전체 카운트의 증가는 유전자 전체의 집합적 수치가 어떻게 특정 스플라이싱 이벤트(예: 넌센스 매개 분해를 유도하는 암호적 엑손 포함)를 가리거나 모순될 수 있는지를 보여줍니다. 저자들은 접합부 해상도 또는 이소폼 특이적 데이터 없이는 이 카운트들로부터 UNC13A 암호적 스플라이싱에 관한 어떠한 메커니즘적 주장도 도출할 수 없다고 단언합니다.

인정된 한계점
저자들은 본 연구가 12개의 라이브러리, 7명의 기증자라는 작은 코호트 규모, 건강한 대조군의 부재, 그리고 전두엽 피질 신경세ola 핵으로 국한되었다는 점에 의해 제한된다는 점을 강조합니다. 결정적으로, 분석은 유전자 수준의 카운트에 의존하며, 암호적 엑손 포함, PSI 또는 단백질 풍부도를 측정할 수 없습니다. 원본 FASTQ 파일은 재처리되지 않았으며, 소스 데이터의 정확한 게놈 빌드 또는 주석 릴리스를 독립적으로 재구성할 수 없었습니다. 따라서 본 연구 결과는 특정 질환 코호트 내의 연관성을 설명하며, 인과적 메커니즘을 확립하기 위해서는 접합부 수준 및 실험적 검증이 필요합니다.