Mammalian TRIM21 Exhibits Gene-Wide Episodic Diversifying Selection with Statistically Supported and Structurally Positioned Evolutionary Candidate Residues

기술 요약: 포유류 TRIM21의 진화 유전체학

문제 제기
TRIM21(Ro52)은 세포질 내 IgG Fc 수용체이자 E3 유비퀴틴 연결효소(E3 ubiquitin ligase)로서, 항체 의존적 세포 내 중화에 결정적인 역할을 한다. 면역 효과 유전자들은 종종 숙주-병원체 간의 군비 경쟁으로 인해 양성 선택(positive selection)의 징후를 보이지만, 포유류 전반에 걸친 TRIM21의 진화 역학은 아직 해결되지 않은 상태이다. 구체적으로, 다양성 선택(diversifying selection)이 유전자 전체 수준, 특정 계통 수준, 또는 개별 잔기(residue) 수준에서 작용하는지, 그리고 그러한 선택이 구조적 맥락이나 반복적인 진화적 사건에 의해 뒷받침되는지는 알려져 있지 않다. 본 연구는 구조 분석 및 조상 상태 분석을 통합하여, 유전자 수준, 계통 수준, 잔기 수준의 증거를 명확히 구분함으로써 이러한 질문들을 해결하고자 한다.

방법론
본 연구는 MAFFT를 통해 정렬되고 PAL2NAL을 통해 코돈으로 역번역된 96개의 포유류 TRIM21 코딩 서열(631개 코돈) 데이터셋을 활용하였다. 최대 가능도 계통수(maximum-likelihood phylogeny)는 JTT+G4 모델 하에서 IQ-TREE2를 사용하여 추론되었다.

분석 프레임워크는 다층적 접근 방식을 채택하였다:

1. 유전자 전체 선택: BUSTED (HyPhy)를 사용하여 유전자 전체에 걸친 에피소드적 다양성 선택(episodic diversifying selection)을 테스트하였다. 전역적 dN/dS 비율은 SLAC을 통해 추정되었다.
2. 계통 특이적 선택: 16개 분류군 서브셋에 aBSREL (HyPhy)을 적용하여 에피소드적 선택이 일어난 가지(branch)를 식별하였으며, Holm–Bonferroni 교정을 사용하였다. PAML branch-site 모델도 시도되었으나 검증 문제로 인해 최종 결과에서는 제외되었다.
3. 잔기 수준 선택: 96개 분류군 정렬에 대해 FEL, FUBAR, SLAC (HyPhy)을 실행하였다. FEL p-값은 Benjamini–Hochberg (FDR)를 사용하여 조정하였다. 보완적 검증을 위해 포유류 서브셋에 PAML site 모델(M7 vs. M8)을 적합시켰다.
4. 구조 및 조상 분석: 잔기 좌표를 인간 TRIM21–IgG Fc 결정 구조(PDB 2IWG)에 매핑하였다. 특히 잔기 453에 대해 조상 상태 재구성(IQ-TREE2, JTT 모델)과 네 가지 별도의 귀무 시뮬레이션(N=100,000)을 사용하여 관찰된 G→R 전이의 통계적 특이성을 조사하였다.
5. 검증: 분류군 샘플링에 따른 잔기 수준 호출(call)의 성능과 민감도를 평가하기 위해 일치하는 분류군 교차 프레임워크 테스트(8개 분류군)를 수행하였다.

주요 결과

• 유전자 전체 선택: BUSTED는 유전자 전체의 에피소드적 다양성 선택에 대한 강력한 증거(LRT = 266.78, p < 0.0001)를 제공하였으며, 정화 선택(purifying selection, ω ≈ 0.06)을 배경으로 하는 강한 양성 선택(ω ≈ 134)을 보이는 소수의 부위를 식별하였다. PAML M7 vs. M8 테스트가 서브셋에서 이를 뒷받침하였다 (2ΔlnL = 80.28, p = 3.69 × 10⁻¹⁸).
• 계통 이질성: aBSREL은 다양한 계통(Perissodactyla, Afrotheria, Lagomorpha, Strepsirrhini, Cetacea, Eulipotyphla)에 걸친 7개의 가지에서 유의미한 에피소드적 선택을 식별하였다. 그러나 신호는 분류군 샘PLING에 민적인 영향을 받았는데, 코끼리(elephant) 신호는 접근성(accession)에 의존적이었으며, 하나의 내부 가지 신호는 위상적 아티팩트(topological artifact)로 표시되었다.
• 잔기 수준 후보: 96개 분류군 분석에서 FEL(FDR 교정 q ≤ 0.05)을 통과한 잔기는 R57, G453, H336, E252, K455, Q308 총 6개였다.
  • R57: 가장 강력한 통계적 후보로 식별되었으며, FEL (q = 0.0028), FUBAR (posterior = 0.983), SLAC (p = 0.0022) 모두에서 일치하는 결과를 보였다. 이 잔기는 결정화된 PRY/SPRY 도메인 외부인 RING–B-box 링커 영역에 위치한다.
  • G453: 중간 정도의 통계적 지지(FEL q = 0.0145; SLAC p = 0.042; FUBAR는 임계값 미달)를 보였다. 그러나 이는 PDB 2IWG에서 IgG Fc 영역과의 백본 매개 접촉(3.50 Å)으로서 구조적으로 유의미하다. 또한 가장 높은 엔트로피(2.433 bits)를 가진 인터페이스 잔기이다.
• G453에서의 반복 진화: 조상 상태 재구성을 통해 두 개의 고신뢰도 독립적 G→R 전이(Bovidae 및 Ursidae 계통)를 확인하였다. 사이트 동질적 JTT 과정 하에서의 귀무 시뮬레이션(N=100,000)은 관찰된 ≥2회의 그러한 사건이 통계적으로 드문 현상임을 나타냈다 (귀무 모델에 따라 P = 0.00031 ~ 0.0329).
• 검증 능력: 일치하는 8개 분류군 검증은 96개 분류군 분석에서 발견된 유의미한 잔기 수준 호출을 재현하는 데 실패하였으며, 이는 잔기 수준의 탐지가 분류군 샘플링 및 통계적 검정력에 크게 의존함을 확인시켜 주었다.

의의 및 주장
본 논문은 유전자, 계통, 잔기 수준에서 포유류 TRIM21의 진화 역학을 최초로 규명했다고 주장한다. 주요 결론은 포유류 TRIM21이 잘 뒷받침되는 유전자 전체의 에피소드적 다양성 선택을 받는다는 것이다.

잔기 수준에서 저자는 통계적 강도와 구조적 근거를 구분한다:

• R57은 가장 강력한 통계적 후보로 제시되나, 알려진 Fc 결합 도메인 외부에 위치하므로 기능적 메커니즘은 정의되지 않았다.
• G453은 반복 진화에 대한 가장 강력한 구조적 근거를 가진 가설로 제시된다. 이의 식별은 통계적 증거, Fc 인터페이스의 백본 접촉으로서의 위치, 높은 변이성, 그리고 관찰된 G→R 치환 패턴의 통계적 희귀성의 결합에 기초한다.

저자는 본 결과가 기능적 적응, 병원체 유도 선택, 또는 Fc 친화력의 변화를 입증하는 것이 아님을 명시적으로 밝힌다. 대신, 본 연구는 검증 가능한 예측을 생성한다: R57 변이는 단백질 안정성 및 RING 매개 신호 전달에 미치는 영향을 테스트해야 하며, G453 변이는 IgG Fc 결합 및 항체 의존적 세포 내 중화의 변화를 테스트해야 한다. 본 논문은 진화적 신호는 견고하지만, 기능적 해석은 실험적 검증이 필요한 가설임을 강조한다.