Phylogeny-associated codon evolution and conditional boundary dinucleotide structure in Platycerium plastomes

기술 요약: Platycerium 플라스트로옴(plastomes)의 계통 연관 코돈 진화 및 조건부 경계 디뉴클레오타이드 구조

문제 제기
플라스트 게놈은 세포 소기관의 조상 정보와 계통 다양성을 보존하고 있으나, 동의 치환(synonymous substitutions, 아미노산 정체성을 변화시키지 않는 치환)은 흔히 진화적으로 특징이 없거나 단순히 돌연변이 편향에 의해 결정되는 것으로 간-주된다. 플라스트 게놈 연구의 핵심 과제는 코돈 사용 패턴이 공유된 계통 역사(lineage history)를 반영하는 것인지, 아니면 단순히 국소적인 서열 제약 및 조성적 보존(compositional conservation)에 의한 것인지를 구별하는 것이다. 또한, 디뉴클레오타이드 빈도(CpG 및 TpA 등)가 식물 게놈에서 비무작위적이라는 사실은 알려져 있으나, 이러한 패턴이 코돈 인벤토리 제약(아미노산 요구량 및 뉴클레오타이드 조성)에서 기인하는지, 아니면 코돈 경계에서의 특정 동의 코돈 배열에서 기인하는지는 여전히 불분명하다. 본 연구는 양치식물 속인 Platycerium을 대상으로 두 가지 규모의 동의 진화를 다룬다: (1) 종 수준의 코돈 프로파일이 좁은 범위의 게놈 조성에도 불구하고 계통적 거리를 추적하는지 여부, 그리고 (2) 엄격한 조건부 널 모델(단백질 서열과 정확한 코돈 수를 고정한 모델)과 비교했을 때 코돈 경계에서 특정 디뉴클레오타이드가 결핍되어 있는지 여부이다.

연구 방법론
본 연구는 세 가지 지리적 계통(Afro-American, Javan-Australian, Malayan-Asian)을 대표하는 18종의 Platycerium을 분석하였다. 연구진은 발표된 전장 플라스트로옴 정렬(whole-plastome alignment)을 활용하였으며, 84개의 고유한 단백질 코딩 특징을 추출하여 무결성을 검토하였다. 주요 분석 매트릭스는 18종 전체에 걸쳐 편집이 필요 없는 온전한 44개의 유전자로 구성되었다(총 792개 종-유전자 서열).

1. 계통적 코돈 분석: 59개 센스 코돈에 대한 상대적 동의 코돈 사용도(RSCU)를 계산하였다. RSCU 벡터 간의 유클리드 거리를 계산하여 참조 플라스트로옴 트리와 맨텔 검정(Mantel test)을 통해 비교하였다. RSCU 기반의 최근접 이웃(neighbor-joining) 표현형도(phenogram)를 구축하여 참조 트리와의 위상학적 중첩을 평가하였다.
2. 유전자 수준 진단: 돌연변이 대 선택적 영향력을 평가하기 위해 유효 코돈 수(ENC), 3번째 위치의 GC 함량(GC3), 그리고 패리티 규칙-2(parity-rule-2) 좌표를 포함한 지표들을 분석하였다. 샘플링 효과와 진화적 분화를 구분하기 위해 유전자 길이와 RSCU 발산 사이의 관계를 테스트하였다.
3. 정확한 코돈 조건부 무작위화: 코돈 배열의 효과만을 격리하기 위해, 각 아미노산 그룹 내에서 동의 코돈을 9,999회 셔플링하는 널 모델을 적용하였다. 이 절차는 아미노산 서열, 모든 코돈의 정확한 개수, 총 뉴클레오타이드 조성, 위치별 GC 함량, 그리고 코돈 내부의 디뉴클레오타이드를 그대로 보존하였다. 오직 동의 코돈의 배치와 그에 따른 코돈 경계 디뉴클레오타이드만이 변화하도록 하였다.
4. 경계 디뉴클레오타이드 분석: 코돈 경계에 존재하는 가능한 모든 16가지 DNA 디뉴클레오타이드를 계수하였다. 관측된 횟수를 무작위화된 기대값과 비교하여 결핍 또는 풍부도를 계산하였다. 특히 CpG와 TpA에 주목하였으며, 계층적 부트스트래핑(hierarchical bootstrapping)과 부호 검정(sign tests)을 통해 종 및 유전자 전반에 걸친 통계적 유의성을 평가하였다.

주요 결과

• RSCU의 계통 신호: RSCU 거리는 플라스트로옴의 경로적 거리(patristic distance)와 유의미하게 증가하였다(Mantel $\rho = 0.6685$, $P = 0.0001$). 게놈 GC 함량(0.4090–0.4128)과 ENC 값(51.68–52.34)의 범위가 매우 좁음에도 불구하고, RSCU 표현형도는 세 가지 지리적 계통을 단일 계통 클러스터로 성공적으로 복구하였다. 이러한 연관성은 계통 내 쌍(within-lineage pairs)으로 분석을 제한했을 때도 지속되었으며, 이는 동의적 발산이 깊은 분지와 최근의 분지 모두를 따라 축적됨을 나타낸다.
• 유전자 길이 민감도: 유전자 수준의 RSCU 발산은 코딩 서열 길이와 강한 음의 상관관계를 보였다($\rho = -0.6142$). 짧은 유전자는 더 큰 RSCU 거리를 나타냈는데, 이는 비율 기반 지표가 짧은 로커스(loci)에서의 적은 수의 치환에 매우 민감함을 시사한다.
• 조건부 경계 결핍: 조건부 무작위화 결과, 단백질 서열과 코돈 인벤토리가 고정된 상태에서도 코돈 경계는 교환 가능하지 않음이 밝혀졌다.
  • CpG: CpG는 코돈 경계에서 유의미하게 결핍되었다(관측/기대 비율 = 0.887). 중앙값 효과는 18종 전체와 44개 유전자 중 34개 유전자에서 음수였다. 이러한 결핍은 세 가지 계통, 다양한 재구성 방식, 그리고 유전자 길이 임계값(900 nt 이상 유전자 포함) 전반에서 견고하게 나타났다.
  • TpA: TpA는 통합된 결핍(비율 = 0.933)과 음수의 종 수준 중앙값을 보였으나, 유전자 수준에서는 변동성이 컸으며 44개 유전자 중 26개 유전자에서만 음의 효과를 보였다.
  • 전역적 경 landscape: 16가지 모든 디뉴클레오타이드 분석은 CpG나 TpA에 국한된 효과라기보다 코돈 맥락의 광범위한 재분포를 나타냈다. GpT가 가장 많이 결핍되었으며(0.865), 그다음이 CpG였다.

핵심 기여

1. 조성으로부터 계통의 분리: 본 연구는 조성적 지표(GC, GC3, ENC)가 계통 간에 매우 보존되어 있음에도 불구하고, 코돈 프로파일이 감지 가능한 계통 구조를 유지한다는 것을 입증하였다.
2. 코돈 순서에 대한 조건부 널 모델: 저자들은 정확한 코돈 조건부 무작위화를 구현함으로써 아미노산 조성, 총 코돈 수, 뉴클레오타이드 편향과 같은 혼란 변수로부터 코돈 순서의 신호를 격리하였다. 이를 통해 CpG의 경계 결핍이 단순히 코돈 인벤토리 때문이 아니라 코돈 배치의 특성임을 확인하였다.
3. 계통 신호와 유전자 신호의 구별: 본 연구는 계통 연관 동의 발산(RSCU 표현형도로 포착)과 보존된 국소 서열-맥락 경향(경계 디뉴클레오타이드 결핍으로 포착)을 구분하였다.
4. 짧은 유전자의 통계적 엄밀성: 본 연구는 짧은 유전자가 낮은 치환 횟수로 인해 RSCU 공간에서 매우 높은 발산을 보이는 통계적 인공물(artifact)을 강조하며, 비교 플라스트 게놈 연구에서 길이 기반의 민감도 필터를 사용할 것을 권고한다.

의의 및 주장

본 논문은 Platycerium의 플라스트 코돈 진화가 중첩된 모델에 의해 형성된다고 주장한다: 즉, 조상으로부터 물려받은 조성적 제약(A/T-rich)이 안정적인 틀을 제공하는 가운데, 계통의 역사가 동의 코돈 빈도의 점진적인 변화를 주도한다. 동시에, 특정 코돈 순서(특히 코딩 경계에서의 CpG 결핍)에 대한 보존된 경향성이 단백질 서열이나 총 코돈 사용량과 무관하게 지속된다.

저자들은 이러한 발견이 공유된 플라스트 조성, 조상, 유전자 수준의 제약, 그리고 국소 서열 맥락이 결합하여 코돈 구조를 형성한다는 진화 모델을 뒷받침한다고 주장한다. 모든 샘플링된 계통에서 나타나는 일관된 CpG 결핍은, 이것이 독립적인 진화적 기원이 아니라 보존된 조상의 서열-맥락적 경향이거나 유사한 플라스트 제약 하에서 반복적으로 유지된 결과임을 시사한다. 결론적으로, 동의적 변이는 계통의 역사와 국소 서열-맥락 제약을 모두 기록하며, 이 두 신호를 구별하기 위해서는 종 수준의 프로파일과 유전자 수준의 조건부 순서를 분리해야 한다고 밝히고 있다.