Structure-informed evolutionary analysis of the meiotic recombination machinery

기술 요약: 감수 분열 재조합 기구의 구조 정보 기반 진화 분석

문제 제기
감수 분열에서의 상동 재조합은 성적 생식을 하는 진핵생물의 생식력과 게놈 안정성에 필수적이다. 역설적이게도, 이 과정은 보존되어 있음에도 불구하고 이를 매개하는 단백질 기구는 계통 특이적인 단백질 소실, 새로운 단백질의 출현, 가속화된 아미노산 서열 변화를 포함하여 급격한 진화적 다양성을 보인다. 이러한 가소성을 유도하는 진화적 힘은 여전히 잘 알려져 있지 않다. 양성 선택(positive selection)을 탐지하기 위한 전통적인 방법들은 단백으로를 선형 서열로 취급하여, 서열 변화와 구조적 제약(예: 매몰된 잔기 대 표면 잔기, 질서 있는 영역 대 무질서한 영역) 및 기능적 역할 사이의 상호작용을 무시하기 때문에 민감도가 떨어지는 경우가 많다. 또한, 표준 코돈 기반 분석은 특히 특정 도메인이나 잔기에 따라 선택 압력이 변할 때, 완화된 정화 선택(relaxed purifying selection)과 적응적 양성 선택을 구분하지 못하는 경우가 많다.

연구 방법론
저자들은 고해 resolution 구조 정보를 결합하여 선택의 징후를 더 민감하게 탐지할 수 있는 통합적 프레임워크를 개발하였다.

• 데이터 수집: 연구진은 네 가지 진핵생물 계통(영장류, 설치류, 조류, 출아 효모(Saccharomyces))에 걸쳐 25개의 핵심 감수 분열 재조합 단백질을 분석하였으며, 58종에서 1,156개의 서열을 수집하였다.
• 진화 분석: PAML의 codeml 패키지를 사용하여 동의 치환율($D_S$)과 비동의 치환율($D_N$)을 추정하였다. 저자들은 양성 선택을 받는 코돈을 식별하기 위해 사이트 모델(M1 vs. M2; M7 vs. M8)과 우도비 검정(LRTs)을 사용하였다.
• 구조적 통합: 경험적 모델 및 AlphaFold로 생성된 구조 모델을 진화율 추정치 위에 매핑하였다. 저자들은 잔기를 다음과 같이 분류하였다:
  • 구조적 질서(Structural Order): 구조화된 도메인과 본질적 무질서 영역(IDR)을 구분하기 위해 예측된 국소 거리 차이 테스트(pLDDT) 점수를 사용하였다.
  • 용매 접근성(Solvent Accessibility): 매몰된 잔기와 표면 잔기를 구별하기 위해 상대적 용매 노출 표면적(rel. SESA)을 사용하였다.
  • 단백질-단백질 인터페이스: AlphaFold3 예측을 사용하여 다단백질 복합체(예: SPO11-TOP6BL, MRX/N, MutS$\gamma$, MutL$\gamma$) 내의 접촉 잔기를 식별하였다.
• 파라로그 비교: 배경 선택 압력을 제어하기 위해, 저자들은 감수 분열 특이적 파라로그(예: DMC1, MSH4)를 "일반론적" 대응체(RAD51, MSH2)와 비교하여, 구조적으로 일치하는 잔기들 사이의 서로 다른 진화율을 식별하였다.
• 실험적 검증: 서열 변이의 기능적 관련성은 Saccharomyces cerevisiae를 이용한 교차 보완 분석(cross-complementation assays)을 통해 테스트되었다. 서로 다른 Saccharomyces 종(S. paradoxus, S. bayanus, S. arboricola)에서 유래한 FLAG 태그가 붙은 MSH4 유전자를 msh4$\Delta$ 배경에서 발현시켜 포자 생존력, 재조합 효율(유전적 거리), 그리고 감수 분열 비분리 현상을 평가하였다.

주요 결과

1. 계통 특이적 양성 선택: 효모는 (동의 치환 포화로 인해) 유의미한 양성 선택을 보이지 않았으나, 설치류, 조류, 영장류 계통에서는 각각 15, 7, 6개의 단백질에서 유의미한 양성 선택의 징후가 나타났다. 주요 단백질로는 SPO11, TOP6BL, IHO1, REC114, MRE11, RAD50, MSH4, HFM1, RNF212, BLM, MUS81, 및 MLH3가 포함된다.
2. 무질서 영역에서의 선택: MLH3에서 양성 선택을 받는 잔기들은 구조화된 ATPase 도메인과 이량체화 도메인을 연결하는 본질적 무질서 링커 도메인에 유의미하게 농축되어 있었다. 이 영역은 모든 계통에서 구조화된 도메인보다 높은 $D_N/D_S$ 비율을 보였으며, 이는 완화된 구조적 제약이 빠른 진화를 허용함을 시사한다.
3. 구조적 제약에도 불구하고 나타나는 선택: 연구는 양성 선택을 받는 잔기들이 매몰되어 있고 구조적으로 질서 정연하다는 사실을 확인하였는데, 이는 구조적 영역은 엄격히 보존되어야 한다는 가설에 도전하는 것이다. 주목할 만한 사례는 다음과 같다:
   • MEI1: 구조적 핵심부 전체에 분포된 46개의 $D_N/D_S > 1$ 잔기.
   • TOP6BL: 설치류에서 특히 두드러지는, GHKL 및 트랜스듀서 도메인 내의 매몰된 잔기들의 양성 선택.
   • RAD50: 조류 계통에서 매몰되어 있음에도 불구하고 유의미한 양성 선택을 보인 4개의 잔기(Walker A 모티프 및 코일드 코일 도메인 포함).
4. 단백질-단백질 인터페이스 역학: 단백질-단백질 인터페이스의 접촉 잔기들은 일반적으로 강한 정화 선택을 경험한다. 그러나 특정 인터페이스 잔기들(예: SPO11-TOP6BL 복합체 내)은 완화된 선택 또는 양성 선택의 징후를 보였다. 설치류의 경우, SPO11-TOP6BL 인터페이스에서의 가변적인 아미노산 조합이 상호작용을 변화시키거나 방해할 것으로 예측되었으며, 이는 계통 특이적 공진화를 시사한다.
5. 파라로그 특이적 발산: 감수 분열 특이적 파라로그(DMC1, MSH4)는 일반론적 대응체(RAD51, MSH2)보다 일관되게 높은 $D_N/D_S$ 비율을 보였다. 구조적 매핑 결과, 감수 분열 특이적 단백질들은 일반론적 단백질에서 보존되는 특정 표면 루프(예: 설치류 DMC1의 경우)에 변이를 축적하는 것으로 나타났다.
6. 서열 변이의 기능적 영향: 교차 보완 실험은 MSH4의 서열 발산이 기능에 영향을 미친다는 것을 입증하였다. S. arboricola 및 S. bayanus 유래의 이종(heterologous) MSH4는 S. cerevisiae 또는 S. paradoxus 대조군에 비해 감소된 포자 생존력, 감소된 유전적 재조합 거리, 그리고 증가된 감수 분열 비분리 현상을 초래하였다. 결함의 심각도는 진화적 거리에 상관관계가 있었다. 흥미롭게도, S. bayanus의 MSH4와 이종 MSH5를 공동 발현했을 때 결함이 악화되었는데, 이는 이종 하위 단위 간의 부적합성을 시사한다.

의의 및 주장
본 논문은 진화적 가소성 자체가 감수 분열 재조합 기구의 핵심적인 보존 특성이라고 주장한다. 저자들은 이 기구가 엄격한 서열 보존을 통해서가 아니라, 구조적 제약과 적응적 혁신 사이의 역동적인 균형을 통해 기능을 유지한다고 제안한다.

이 연구의 주요 기여는 다음과 같은 기계론적 프레임워크를 통합한 데 있다:

• 구조적 제약 (질서/무질서, 매몰/표면).
• 단백질-단별 상호작용 인터페이스.
• 서로 다른 기능적 맥락을 가진 파라로그의 비교 분석.

이러한 접근 방식은 단백질을 선형 서열로 분석하거나 유전자 전체의 비율을 평균화할 때 가려질 수 있는 "잠재적(cryptic)" 양성 선택 및 완화된 제약을 탐지할 수 있게 한다. 저자들은 이러한 구조 정보 기반 전략이 감수 분열 이외의 다른 단백질 시스템에서도 기능적으로 유의미한 진화적 변화를 탐지하는 데 일반화되어 적용될 수 있다고 상정한다. 결론적으로, 재조합의 핵심 메커니즘은 고대부터 존재해 왔으나, 그 구체적인 분자적 구현은 계통 특이적 적응과 특정 구조적 영역에서의 제약 완화에 의해 구동되는 매우 가소성이 높은 체계이다.