Spatial transcriptomic programs relate to spectrolaminar rhythms across macaque cortex

기술 요약: 마카크 피질 전반의 공간적 전사체 프로그램과 스펙트로라미나 리듬의 관계

문제 정의
영장류의 대뇌 피질은 피질 깊이에 따라 미세회로, 투사 패턴, 분자 구성이 체계적으로 변화하는 양상을 보인다. 마카크 전기생리학에서 잘 알려진 "스펙트로라미나 모티프(spectrolaminar motif)"—고주파 전력이 표층에서 정점을 찍고 $\alpha$-$\beta$ 전력이 심층에서 정점을 찍는 재현 가능한 패턴—는 확립되어 있으나, 그 기저에 있는 세포 및 분자적 기질은 여전히 미해결 상태이다. 기존 연구들은 전형적인 세포 마커(예: PV, NRGN)를 진동 전력과 연결 지었으나, 이러한 저차원적 접근 방식은 고차원적인 전사체 하위 부류, 지역적 세포 유형 구성, 그리고 유전자 발현 프로그램을 해상할 수 없다. 핵심적인 과제는 곡선 형태의 공간 전사체 데이터와 프로브 삽입 각도 및 피질 곡률에 따라 변하는 전기생리학적 깊이 프로파일을 정렬할 공통 좌표계의 부재이다. 본 연구는 일반적인 피질 구배(gradient)와 층 기하학을 넘어, 층별 전사체 프로그램이 피질 깊이에 따라 표현되는 스펙트럼 프로세스를 예측할 수 있는지 여부를 다룬다.

방법론
저자들은 12개의 일치된 마카크 피질 영역에서 국소 전위(LFP) 기록과 단일 세포 공간 전사체를 통합하는 "스펙트로-오믹(spectro-omic)" 프레임워크를 개발하였다.

1. 전기생리학적 정렬 및 분해 (LSE):

   • 정렬: LFP 프로브는 vFLIP-v2 프레임워크를 사용하여 정렬되었으며, 이는 $\alpha$-$\beta$/고주파 교차 지점(제4층 근처)에 깊이 축을 고정하여 공통의 $z=0$ 기준점을 생성한다.
   • 분해: 고정된 대역 전력 분석 대신, 저자들은 **국소 스펙트럼 확장(Local Spectral Expansion, LSE)**을 채택하였다. 이 가법 모델(additive model)은 주파수-깊이 전력 표면 $Y(f, z)$를 6개의 전형적인 스펙트럼 프로세스($\delta, \theta, \alpha, \beta, \text{low-}\gamma, \text{high-}\gamma$)로 분해하며, 국소 로렌츠 커널(Lorentzian kernels)을 사용한다. 이를 통해 프로세스별 주파수 중심, 깊이 중심 및 공간적 확산을 추정함으로써 스펙트로라미나 모티프의 기하학적 구조를 보존한다.

2. 전사체 정렬 (L4P):

   • 변환: 원래 곡선 형태의 2D 섹션에서 측정된 공간 전사체 세포들은 **제4층 투영(Layer-4 Projection, L4P)**을 통해 변환되었다. 이 절차는 경험적 제4층 중간 지점을 중심으로 피질 리본을 펼쳐, 세포들을 LSE 깊이 그리드와 일치하는 평탄한 부호가 있는 깊이 좌표계(표층은 음수, 심층은 양수)로 매핑한다.
   • 매칭: 전사체 데이터셋과 전기생리학 데이터셋은 개별 동물이 아닌 종, 피질 영역, 층별 깊이 수준에서 매칭되었으며, 프로브와 칩을 해부학적 영역 내의 독립적인 샘플로 취급하였다.

3. 예측 모델링:

   • 베이스라인 통제: 계층 구조 좌표, 층 좌표 및 이들의 상호작용을 사용하여 광범위한 피질 구배를 설명하기 위한 "계층-플러스-층(hierarchy-plus-layer)" 베이스라인 모델을 구축하였다.
   • 잔차 예측: 홀드아웃(held-out) 층별 스펙트럼 구성을 예측하기 위해 부분 최소 제곱(Partial Least Squares, PLS)을 사용한 **부분 로짓 모델(partwise logit model)**을 적용하였다. 전사체 예측 변수들은 일반적인 구배 이상의 분산을 설명하는지 테스트하기 위해 베이스라인에 대해 잔차화되었다.
   • 검증: 유의성 평가는 지역 제외 교차 검증(Leave-One-Region-Out, LORO), Freedman–Lane 지역 블록 잔차 순열 검정, 그리고 계층 보존 이동/반사 널(null) 모델을 사용하여 단순한 공유 기하학에 의한 결과가 아님을 확인하였다.

4. 기능적 디코딩:

   • 전체 깊이 PLS1: 안정적이고 부호가 있는 유전자 프로그램(PLS1+ 및 PLS1-)을 식별하기 위해 21개 빈(bin)의 전체 깊이 프로파일에 일성분 PLS 모델을 적용하였다.
   • 농축 분석: 안정적인 유전자 세트는 유전자 온톨로지(GO) 생물학적 과정 및 세포 유형 하위 부류 농축 분석에 사용되었다.

주요 결과

• 예측력: 전사체 구조는 계층-플러스-층 베이스라인을 넘어 홀드아웃 층별 스펙트럼 구성을 유의미하게 개선하였다. 전체 모델은 베이스라인의 $R^2 = 0.384, r = 0.692$와 비교하여 $R^2 = 0.621, r = 0.809$를 달성하였다($\Delta R^2 = +0.237$). 이 개선은 통계적으로 유의미하였으며($p_{FL} = 0.001$), 계층 보존 이동/반사 널 모델을 초과하였다($p_{shift} = 0.042$).
• 프로세스 특이성: 예측 이득은 더 빠른 리듬($\beta, \text{low-}\gamma, \text{high-}\gamma$)에서 가장 두드러졌으며, 느린 리듬($\delta, \theta$)은 더 지역 특이적인 이질성을 보였다.
• 분자적 시그니처:
  • $\alpha/\beta$ 리듬: 심층(L4/5/6, L6) 글루타메이트성(glutamatergic) 프로그램과 연관되었다. 이들은 부착, 축삭 유도, 시냅스 형성, 수초화 및 영양 신호 전달(예: NFASC, PTPRD, ERC2, EFNA5, RALYL, AFDN)과 관련된 유전자들이 농축되었다. 이는 느린 심층 리듬이 피드백 경로 및 구조적 통합과 결합되어 있다는 견해를 뒷받침한다.
  • $\gamma$ 리듬 (Low/High): 표층-중층 GABAergic 프로그램, 특히 PVALB가 풍부한 억제성 하위 부류와 연관되었다. 이들은 흥분성, 이온 항상성, 에너지 대사 및 시냅스 소포 방출(예: MEF2C, SCN1A, SLC12A5, ATP1A1, GABRA1, SNAP25)과 관련되었다.
  • 비신경세포 기여: 두 스펙트럼 영역 모두 비신경세포 프로그램(희소돌기아교세포, 별아교세포, 내피세포)의 농축을 보여, 조직 환경(수초, 대사 지원, 세포외 기질)이 LFP 스펙트럼 구조를 형성함을 시사한다.

의의 및 주장
본 논문은 스펙트로라미나 모티프를 순수하게 생리학적인 관찰에서 예측 가능하고 세포 유형이 해상된 스펙트로-오믹 프레임워크로 확장한다고 주장한다. 주요 의의는 층별 분자 구조가 마카크 피질 전반의 스펙트로라미나 리듬에 대한 타당한 기질을 제공한다는 점을 입증한 데 있다.

• 기하학적 기여: 본 연구는 이질적인 양상(LFP와 공간 전사체)을 공통의 제4층 참조 축에 정렬하는 것이 교차 모달리티 추론에 필수적임을 확립하였다.
• 기전적 통찰: 연구 결과는 피질 리듬에 대한 이해를 정교화한다: $\alpha/\beta$ 활동은 심층의 흥분성, 피드백 관련 및 구조적 아키텍처와 결합되어 있는 반면, $\gamma$ 활동은 억제성, 흥분성 및 대사 프로그램과 결합되어 있다.
• 한계점: 저자들은 본 연구가 관찰 연구이며 특정 유전자에 의한 리듬의 인과적 생성을 확립하는 것은 아니라고 명시하였다. 데이터셋은 영역과 깊이에 따라 매칭되었으나 동일한 조직에서 샘플링된 것은 아니므로 직접적인 인과 추론은 불가능하다. LSE 모델은 고유하게 식별 가능한 회로 모델이라기보다 구조화된 근사치이다.

결론적으로, 본 연구는 피질 기둥(cortical column)이 층별 스펙트럼 역동성이 지역 특이적인 세포 및 분자 아키텍처(신경 및 비신경 성분 포함) 내에 내재되어 있는 공동의 스펙트로-오믹 구조임을 제안한다.