Compartment-dependent selection regimes maintain the G-quadruplex paradox across the eukaryotic kingdoms

기술적 요약: 구획 의존적 선택 체계가 G-쿼드러플렉스 역설을 유지한다

문제 제기
G-쿼드러플렉스(G4) DNA 구조는 진핵생물 게놈 전반에 걸쳐 역설적인 분포를 보인다. 즉, 코딩 서열(CDS)에서는 체계적으로 결핍되어 있는 반면, 프로모터 및 비번역 부위(UTR)와 같은 조절 영역에서는 농축되어 있다. 이러한 "G4 역설"은 특정 계통(특히 동물계)에서 기록되었으나, 이것이 진핵생물 게놈 구조의 보편적 특징인지 아니면 역사적 샘플링 편향에 의한 인위적 결과인지는 여전히 불분명하다. 또한, 네 가지 진핵생물 왕국(동물계, 균계, 식물계, 원생생물계) 전반에서 이 구조를 유지하는 선택 체계는 완전히 이해되지 않았다. 특히, CDS에서의 결핍이 침묵 부위(silent sites)에 대한 구조적 제약에 의한 것인지, 아니면 코돈 적응(codon adaptation)과 같은 혼란 변수에 의한 것인지 명확하지 않으며, 비코딩 영역에서의 농축이 보존된 분자 기구에 의해 관리되는지 혹은 왕국 특이적인지에 대해서도 불확정적이다.

방법론
저자들은 기존 연구의 동물계 편향을 교정하기 위해 네 가지 진핵생물 왕국(동물계: 86, 균계: 49, 식물계: 33, 원생생물계: 30)으로 층화된 균형 잡힌 198개 참조 게놈 샘ku를 분석하였다. 본 연구는 다단계 계산 파이프라인을 사용하였다:

1. G4 탐지 및 안정성 분류: 게놈 영역(5′UTR, CDS, 인트론, 3′UTR, 프로모터, 인터제닉)을 regex 패턴과 G4Hunter 점수의 컨센서스를 사용하여 스캔하였다. 모티프는 Mergny–Lacroix 열역학적 안정성($\Delta G \le -5$ kcal/mol)에 의해 추가로 필터링되어 "컨센서스 안정" G4를 식출하였다.
2. 교차 왕국 농축 분석: 각 왕국에 대해 특정 영역에서의 안정적 G4의 log2 농축 비율(IRR)을 게놈 전체 배경 대비 계산하였다. 엄격한 기준(Rule B)과 타이브림이티브(tie-permissive) 기준(Rule A) 하에서 왕국 간 부호 일관성을 테스트하였다.
3. 코돈 수준 진화 분석 (dN/dS): 각 왕국 내 밀접하게 관련된 종 쌍에 대해, G4 중첩 코돈과 대조 코돈 사이의 동의 치환율($dS$) 및 비동의 치환율($dN$)에 대한 Nei–Gojobori 비교를 수행하였다. 결정적으로, 이 분석은 번역 선택에 의한 동의 치환율의 혼란 효과를 제어하기 위해 코돈 적응 지수(CAI)를 적용한 후 반복되었다.
4. 헬리케이스 공진화 (PGLS): 여섯 가지 G4 해소 헬리케이스 오솔로그 패밀리(BLM/WRN/RECQ, PIF1, FANCJ/BRIP1, RTEL1, DDX5, DDX17)의 복제수를 정량화하였다. 계통 발생적 비독립성을 교정하기 위해 계통 발생 일반화 최소 제곱(PGLS) 회귀를 사용하여 헬리케이스 레퍼토리와 다양한 게놈 구획 내 G4 풍부도 사이의 연관성을 테스트하였다. 견고성은 폴리토미가 없는 하위 집합과 독립적인 Open Tree of Life 백본을 통해 검증되었다.
5. 모티프-연령 층화: G4 모티프를 호스트 유전자 오솔로그 그룹에 매핑하여, 이것이 "고대"(최초 진핵생물 공통 조상(LECA) 이후 보존된) 유전자에 존재하는지 또는 "계통 특이적" 유전자에 존재하는지 결정하였다.

주요 결과

• G4 역설의 검증: 본 연구는 G4 역설이 교차 왕국적 특징임을 확인하였다. 엄격한 기준 하에서, CDS 결핍과 인터제닉 농축은 네 왕국 모두에서 부호 일관성을 보였다. 타이브림이티브 기준 하에서는 인트론 및 프로모터 농축 또한 부호 일치성을 보였으나, 이는 원생생물계와 균계의 일부 종이 해당 영역에서 안정적 G4를 전혀 보유하지 못한 것에 기인한다.
• CDS 결핍은 코돈 적응에 의해 주도된다 (침묵 부위 제약이 아님):
  • 초기 Nei–Gojobori 분석은 7개 중 6개 불포화 종 쌍에서 G4 중첩 코돈에서의 $dS$ 억제를 보여주었다.
  • 그러나, per-gene CAI를 적용하여 조정한 후에는 이러한 억제가 사라졌다. 7개 쌍 중 0개에서 조정된 G4 오즈비(odds ratio)가 1 미만으로 유의미하게 유지되었으며, 5개 쌍에서는 방향이 역전되었다 ($OR > 1$).
  • 결론: 관찰된 $dS$ 억제는 G4 중첩 코돈이 고발현 및 코돈 최적화된 유전자에 우선적으로 위치함에 따른 인위적 현상이다. 선택 압력은 G4 모티프 자체의 특정 구조적 제약이 아니라 코돈 적응에 있다.
• 비코딩 영역에서의 왕국 특이적 헬리케이스 공진화:
  • CDS와 달리, 특정 헬리케이스 연관성이 다중 테스트 교정을 통과하지 못한 것과 대조적으로, 인트론과 프로모터에서는 유의미한 연관성이 나타났다.
  • Benjamini–Hochberg 교정($q < 0.10$) 후에도 유의미하게 남은 다섯 가지 연관성은 다음과 같다:
    • 원생생물계: FANCJ/BRIP1은 프로모터 G4와 음의 상관관계를 가짐; RTEL1 및 BLM/WRN/RECQ는 인트론 G4와 연관됨.
    • 균계: PIF1은 프로모터 G4와 음의 상관관계를 가짐.
    • 식물계: FANCJ/BRIP1은 인트론 G4와 양의 상관관계를 가짐.
  • 주목할 점은 FANCJ/BRIP1이 맥락 의존적 행동을 보인다는 것이다. 즉, 식물에서는 인트론 G4 유지와 복제수가 상관관계가 있는 반면, 원생생물계에서는 프로모터 G4 해소와 상관관계가 있다.
• 인트론 G4의 진화적 깊이: 인트론 G4 모티프는 압도적으로 깊게 보존된 오솔로그 그룹(모든 왕국에서 $\ge 87.5\%$)에 집중되어 있으며, 이는 이러한 구조의 조절 맥락이 진핵생물의 초기 진화 단계부터 유지되어 왔음을 시사한다. 반면, 프로모터 G4는 주로 계통 특이적 유전자에서 발견된다.

의의 및 주장
본 논문은 G4 역설의 유지를 위한 **이중 체계 구조(two-regime architecture)**를 제안한다:

1. 코딩 서열(CDS)에서: 지배적인 선택 압력은 코돈 적응이다. CDS에서의 G4 결핍은 침묵 부위의 직접적인 구조적 제약 때문이 아니라, 고발현 유전자에서 특정 코돈을 선호하는 번역 선택의 부산물이다.
2. 비코딩 영역(인트론/프로모터)에서: 선택은 G4 구조에 직접 작적으로 작용한다. 이러한 구조의 유지는 왕국 특이적으로 배치된 헬리케이스 기구(예: FANCJ/BRIP1, PIF1, RTEL1)에 의해 관리되며, 이는 구획별로 공진화한다.

저자들은 교차 왕국적 G4 역설이 Animalia에 국한된 인위적 결과가 아니라 보존된 조절 구조라고 주장한다. 이 역설의 메커니즘적 핵심은 G4가 활발히 관리되는 비코딩 영역에 존재하며, 이는 아마도 RNA 가공 및 전사 결합(transcription coupling) 역할을 수행할 것이다. 본 연구는 CDS 결핍을 G4 허용성의 실패가 아닌 호스트 유전자의 발현 프로파일에 따른 결과로 재정의하는 한편, 인트론과 프로모터 구획을 기능적 G4 생물학 및 선택적 유지의 진정한 장소로 식별한다.

저자가 언급한 한계점
저자들은 프로모터 정의(고정된 -1000 bp 윈도우)가 키네토플라스티드(원생생물계)의 경우 생물학적으로 이질적일 수 있으며, 이것이 프로모터 특이적 연관성 해석에 영향을 미칠 수 있음을 인정하였다. 또한, 헬리케이스 분석은 특정 파라로그를 놓칠 수 있는 KEKG 오솔로그 주석에 의존하며, 비동물계 왕국에 대한 $dN/dS$ 분석은 단일 종 쌍에 의존하므로 왕국 전반의 일반성을 확인하기 위한 추가적인 재현이 필요하다. 본 연구는 비교 게놈학적 분석이므로, 헬리케이스 수와 G4 풍부도 사이의 상관관계를 보여주기는 하나, 분석된 비모델 종들에서 이들 헬리케이스의 구체적인 메커니즘적 역할을 실험적으로 검증하지는 않았다.