A genus betrayed by its chloroplast: comparative plastomics of 54 Corymbia species exposes plastid paraphyly and microevolutionary signatures

*기술 요약: 54종 Corymbia의 비교 플래스토믹스(Comparative Plastomics)*

문제 제 statement
호주가 원산인 Corymbia 속(Myrtaceae)은 상당한 경제적 및 임업적 가치를 지닌다. 그러나 이들의 분자 분류학은 여전히 논쟁의 여지가 있다. 부분적인 엽록체 DNA(cpDNA) 마커를 이용한 기존 연구들은 Corymbia를 그 자매 속인 Angophora에 대해 측계통(paraphyletic)으로 일관되게 규정해 왔으며, Angophora 종들이 Corymbia 클레이드 내부에 위치하는 결과를 보여주었다. 형태학적 분류와 플래스티드 신호 사이의 이러한 불일치는 역사적인 엽록체 침입(introgression)의 결과로 간주되곤 한다. 그럼에도 불구하고, 선행 연구들은 샘플링의 부족과 부분적 마커에 대한 의존성으로 인해, 오랜 분류학적 불일치를 해결하고 종 간 구분에 진정으로 유익한 게놈 영역을 식별하는 데 한계가 있었다. 또한, 이 속 전반에 걸친 미세진화 역학, 구조적 다형성 및 선택 압위에 대한 포괄적인 이해도 부족한 실정이다.

방법론
본 연구는 GenBank에서 사용 가능한 공공 플래스토옴(plastomes)을 사용하여 54종의 Corymbia에 대한 전례 없는 비교 게놈 분석을 수행하였다. 방법론은 다음과 같은 다각적 접근 방식을 포함한다:

• 주석 달기 및 구조 분석: GeSeq와 CPGVAS2를 사용하여 게놈을 주석 처리하였으며, ChloroPlot을 통해 구조적 매핑을 수행하였다.
• 선택 압력 분석: MAFFT와 DnaSP를 사용하여 81개의 공유 단백질 코딩 유전자(PCGs)에 대한 비동의(Ka) 및 동의(Ks) 치환율을 계산하여 Ka/Ks 비율을 결정하였다.
• 다양성 및 핫스팟 식별: 전체 정렬에 대해 슬라이딩 윈도우 분석(600 bp 윈도우, 100 bp 스텝)을 통해 뉴클레오타이드 다양성($\pi$)을 평가하여 고변이 영역을 식별하였다.
• 구조적 역학: Progressive Mauve를 사용하여 공선성(collinearity)과 신테니(synteny)를 평가하였다. 역반복 서열(Inverted Repeat, IR) 경계의 확장 및 수축은 IRplus를 통해 분석되었다.
• 반복 요소: Vmatch를 사용하여 분산된 반복 서열을 식별하였고, MISA를 통해 단순 서열 반복(SSR)을 검출하였다.
• 코돈 사용: 상대적 동의 코돈 사용 빈도(RSCU)와 코돈 편향성을 RSCU Calculator를 통해 분석하였다.
• 계통 재구축: 두 가지 데이터셋을 사용하여 계통수를 추론하였다: (1) 77개의 공유 PCG, (2) 전체 엽록체 게놈. 분석은 Allosyncarpia ternata, Angophora costata, Eucalyptus grandis, Syzygium aqueum를 외집단으로 하여 IQ-TREE 3에서 1,000회의 부트스트랩 반복을 통해 수행되었다.

주요 결과

• 게놈 보존성 및 구조: 54종의 플래스토옴은 높은 구조적 보존성을 보였으며, 크기는 158,745~160,839 bp 범위였고 GC 함량은 약 36%였다. 모든 종은 전형적인 4부 구조(quadripartite structure)를 유지하였다. 그러나 게놈 경계에서 유의미한 미세진화 역학이 관찰되었다: 오직 27.78%의 종만이 IR 영역(IRa 및 IRb)에서 동일한 크기를 보유하고 있었다.
• 선택 압력: 플래스토옴은 정화 선택(purifying selection)에 의해 지배된다. 대부분의 유전자는 돌연변이가 없었다 (Ka = Ks = 0). ccsA 유전자는 측정 가능한 Ka/Ks 비율(~0.28)을 보인 유일한 로커스였으며, 이는 단백질 기능을 보존하기 위한 강력한 음성 선택을 나타낸다.
• 뉴클레오타이드 다양성: 전체 뉴클레오타이드 다양성은 낮았다 ($\pi$ = 0.00119). 그러나 5개의 유전자 간 영역과 ycf1 유전자 복사본이 고변이 핫스팟으로 나타났다. C. gummifera에서 독특한 구조적 거대 다형성(macropolymorphism)이 확인되었는데, 이는 첫 번째 ycf1 복사본과 ycf1/ndhF 중첩 영역에 걸친 독점적인 ~1.2 kb의 인델(indel)이다.
• 반복 요소: 분산된 반복 서열은 회문(palindromic) 및 직접형(direct) 타입으로 엄격히 제한되었으며, 역방향 및 상보적 반복 서열은 존재하지 않았다. SSR 분석 결과, 모노뉴클레오타이드 반복이 우세했으며, 그 뒤를 이어 트리(tri-) 및 디(di-)뉴클레오타이드가 나타났다.
• 코돈 사용: 번역 기구의 매우 높은 보존성이 관찰되었다. 20개의 인코딩된 아미노산에 대해 54종 중 51종이 동일한 선호 코돈을 공유하였으며, 이는 세 번째 위치에서의 강한 AT 편향을 특징으로 한다. 유일한 예외는 발린(Valine)과 관련된 것으로, 3종(C. umbonata, C. semiclara, C. dichlomophloia)은 대다수가 선호하는 GTA 코돈 대신 GTT를 선호하였다.
• 계통 추론: 두 가지 계통적 접근 방식(77개 유전자 및 전체 게놈) 모두 높은 지지도로 주요 Corymbia 클레이드를 복원하였다. 결정적으로, 두 분석 모두 Angophora costata가 높은 통계적 지지도로 Corymbia 클레이드 내부에 위치함을 보여주었으며(전체 게놈 100, 공유 유전자 83), 이는 Corymbia의 플래스티드 측계통성을 확인시켜 준다.

의의 및 주장
저자들은 본 연구가 현재까지 Corymbia 그룹에 대한 가장 강력한 플래스티드 참조 자료를 제공한다고 주장한다. 주요 의의는 다음과 같다:

1. 측계통성 확인: 본 연구는 플래스티드 게놈의 관점에서 Corymbia가 Angophora에 대해 측계통임을 보여주는 고해상도 증거를 제공하며, 이는 최근의 핵 게놈 분석 결과와 일치함과 동시에 Corymbia의 단계통성에 도전하는 결과이다.
2. 미세진화적 특징 식별: 연구는 게놈이 전반적으로는 안정적이지만, 특정 영역(IR 경계, ycf1, 유전자 간 스페이서)과 구조적 사건(예: C. gummifera의 인델 및 특정 계통의 IR 수축)이 중요한 미세진화적 마커 역할을 한다는 점을 강조한다.
3. 계통 및 육종을 위한 자원: 고변이 핫스팟과 독특한 구조적 특징의 식별은 향후 분류학적 진단 및 분자 진단을 위한 토대를 제공한다. 저자들은 현재 공공 서열의 품질(모호한 염기)로 인해 프라이머 설계가 제한되어 있으나, 식별된 마커들(예: *C. g-indels)은 고정밀 재서열 분석이 적용될 경우 종 식별 및 임업 및 육종 프로그램에서의 추적성을 위한 잠재력을 지니고 있다고 언급하였다.

결론적으로, Corymbia의 엽록체 게놈은 안정적이면서도 역동적인 시스템이며, 특정 구조적 및 서열 변이가 계통 관계를 정교화하고 역사적인 분류학적 불일치를 해결하는 데 필요한 해상도를 제공한다.