Domain- and Lineage-Specific Heterogeneity in Selective Constraint on Gliomedin (GLDN), a Node of Ranvier Gene Restricted to Jawed Vertebrates

기술적 요약: Gliomedin(GLDN)의 도메인 및 계통 특이적 이질성

문제 및 동기
정화 선택(purifying selection)은 거의 모든 단백질 코딩 유전자에 작용하지만, 단일 유전자 전반의 비동의 치환 대 동의 치환 비율($\omega = dN/dS$)은 중요한 진화적 세부 사항을 가릴 수 있다. 진화 분석의 효용은 선택적 제약이 어디에 존재하는지를 규명하는 데 있다: 즉, 특정 단백로 도메인 전반, 특정 계통 내, 그리고 유전자의 깊은 진화적 기원과 관련하여 말이다. Gliomedin(GLDN)은 이러한 접근 방식에 대한 이상적인 테스트 케이스가 된다. GLDN은 슈반 세포 유래 당단백질로서 랑비에 결절(Nodes of Ranvier)의 분자적 조립에 필수적이며, 축삭 부착 분자인 neurofascin-186(NF186) 및 NrCAM과 결합한다. 이 단백질은 N-말단 콜라겐 유사 도메인과 C-말단 olfactomedin(OLF) 도메인이라는 뚜렷한 2개 도메인 구조를 가지고 있다. 저자의 이전 임상 연구에서는 GLDN 유전자 좌위와 방광암 위험 사이의 암시적인 서브 게놈 수준의 연관성과 공격적인 종양에서 핵 내 GLDN 발현이 감소하는 현상을 확인한 바 있다. 그러나 GLDN이 척추동물 진화 과정에서 어떻게 선택 압력을 받아 형성되었는지, 그리고 이러한 압력이 기능적 도메인에 따라 어떻게 분포되는지에 대한 근본적인 질문은 체계적이고 정량적인 분자 진화 분석을 통해 해결되지 않은 상태로 남아 있었다.

방법론
본 연구는 539개의 턱이 있는 척추동물 GLDN 코딩 서열(CDS)로 구성된 큐레이션된 데이터셋을 사용하여 비교 유전체학 프레임워크를 채택하였다. 방법론은 다음과 같은 엄격한 단계들을 거쳤다:

1. 데이터셋 큐레이션 및 상동성 검증: NCBI에서 검색된 937개의 GLDN 주석 서열을 길이, 완전성 및 분류학적 대표성을 기준으로 필터링하였다(과잉 대표된 지느러미가 있는 어류의 다운샘플링 포함). 결정적으로, HMMER과 계통 발생적 거리 분석을 사용하여 9개의 병렬 olfactomedin 패밀리 유전자(예: MYOC, OLFM1–4, OLFML1–3) 패널에 대해 상동성을 검증하였다. 이를 통해 539개의 후보가 모두 진정한 GLDN 상동 유전자임을 확인하였다.
2. 계통수 및 정렬 구축: FastTree 2를 사용하여 최대 가능도(maximum-likelihood) 계통수를 생성하였다. OLF 도메인의 단백질 서열을 MAFFT로 정렬하였고, 그에 따른 코돈 정렬을 구축하였다.
3. 선택 분석:
   • PAML (codeml): $\omega$를 추정하기 위해 최대 가능도 코돈 모델을 사용하였다. 모델에는 양성 선택을 탐지하기 위한 one-ratio 모델(M0), site 모델(M1a, M2a, M7, M8), 그리고 포유류와 비포유류를 비교하기 위한 two-ratio branch 모델이 포함되었다. 계산적 제약을 관리하기 위해 브랜치 길이는 FastTree 추정치를 기반으로 고정되었다.
   • HyPhy 교차 검증: 견고성을 보장하기 위해 독립적인 방법인 Fixed Effects Likelihood (FEL) 및 Mixed Effects Model of Evolution (MEME)을 사용하여 사이트 수준의 선택 신호를 교차 검증하였다.
4. 구조 매핑: PAML에 의해 식별된 후보 사이트들을 인간 GLDN OLF 도메인의 결정 구조(PDB 5YBY)에 매핑하여 구조적 맥락(B-factor, 2차 구조, 리간드 근접성)을 평가하였다.
5. 무악류 척추동물 탐색: 칠성치(Petromyzontidae) 및 혜골(Myxini)에 대한 이름 기반 및 BLASTP 검색을 유전체 자원을 통해 수행하여, 무악류 척추동물 외부에 GLDN 상동 유전자가 존재하는지 확인하였다.
6. 도메인 비교: OLF 도메인과 인접한 콜라겐 유사 도메인에 대해 각각 별도로 선택적 제약을 추정하였다.

주요 결과
분석 결과, 단일 유전자 전반의 지표로는 보이지 않는 세 가지 뚜렷한 선택적 제약의 이질성이 드러났다:

1. 도메인 특이적 제약: 리간드 결합을 담당하는 OLF 도메인은 인접한 반복적인 콜라겐 유사 도메인($\omega = 0.365$)에 비해 유의미하게 강한 정화 선택($\omega = 0.165$)을 받는다. 이 차이는 2배 이상이며, 이는 변이성이 높은 콜라겐 반복 서열에 비해 OLF 도메인이 갖는 구조적 경직성과 기능적 필수성을 반영한다.
2. 계통 특이적 제약: 포유류 계통은 비포유류 턱 있는 척추동물($\omega = 0.170$)에 비해 OLF 도메인에서 유의미하게 강한 정화 선택($\omega = 0.130$)을 보인다. 이는 단순한 시간 경과에 따른 제약 완화 모델과 상반되며, 포유류에서 기능적 중요성이나 복잡성이 증가했음을 시사한다.
3. 진화적 기원 및 존재 여부: 무악류 척추동물(칠성치 또는 혜골)에서는 GLDN 상동 유전자가 검출되지 않았다. 다른 olfactomedin 패밀리 병렬 유전자들은 이들의 게놈에 존재하지만, GLDN은 압축된 수초(compact myelin)를 독특하게 보유한 그룹인 턱 있는 척추동물(Gnathostomata)에 국한되어 나타난다.
4. 양성 선택 신호: PAML M7/M8 모델은 단일 후보 사이트(참조 서열의 Leu545, 정렬 사이트 288)가 $\omega > 1$일 확률이 99% 이상의 후험 확률을 가짐을 지적하였다. 그러나 독립적인 HyPhy 방법(FEL 및 MEME)은 이 신호를 재현하지 못했다; 오히려 FEL은 이 사이트에서 정화 선택을 나타냈다. 구조 매핑 결과, 이 잔기는 OLF 도메인의 무질서한 C-말단 경계에 위치하며, 높은 B-factor를 보이고 리간드 결합 인터페이스와 근접하지 않았다. 저자는 이 신호가 진정한 양성 선택이라기보다는 유연한 도메인 경계에서의 구조적 제약 완화를 반영하는 것으로 결론지었다.

의의 및 주장
본 논문은 GLDN의 기원 이후 작용해 온 선택 압력에 대한 최초의 정량적 설명을 제공한다고 주장한다. 본 연구의 주요 의의는 "유전자 전반"의 보존성 진술만으로는 필수 유전자의 진화 역사를 이해하기에 불충분하다는 점을 입증했다는 데 있다. 대신, 본 연구는 다음을 확립하였다:

• 선택적 제약은 매우 이질적이며, 기능적 OLF 도메인과 구조적 콜라겐 도메인 사이에서 2배 이상의 차이를 보인다.
• 제약은 정적인 것이 아니라 특정 계통(포유류)에서 강화되며, 이는 신경계의 복잡성과 상관관계가 있을 수 있다.
• GLDN의 출현은 압축된 수초를 보유하지 않은 무악류 척추동물과 달리, 턱 있는 척추동물의 기원과 진화적으로 일치한다.
• 양성 선택의 후보 사이트는 적응적 핫스팟이라기보다 구조적 유연성을 나타내는 영역으로 해석하는 것이 더 적절하며, 이는 독립적인 방법 및 구조 데이터와 함께 선택 신호를 교차 검증할 필요성을 강조한다.

저자는 GLDN의 암 관련 역할이 조직 맥락에 따라 달라질 수 있음을 언급하며, 핵심 OLF 도메인의 강력한 보존성을 바탕으로 질병 관련 메커니즘이 빠르게 진화하는 잔기보다는 정형적인 구조적 아키텍처에 의존할 것임을 명시하며 신중한 입장을 유지한다. 본 연구는 새로운 치료적 응용을 제안하는 것이 아니라, 인간 질병에서 GLDN의 역할을 해석하기 위한 정량적 진화 프레임워크를 제공한다.