Comparative plastome genomics reveals phylogeny, biogeography, and character evolution of Primula (Primulaceae)

기술 요약: Primula 속의 비교 엽록체 게놈학

문제 제기
500종 이상의 종을 포함하는 Primula(앵초과) 속은 급격한 방사 진화, 형태적 수렴 진화(homoplasy), 빈번한 잡종 형성, 그리고 불완도 계통 분류(incomplete lineage sorting, ILS)로 인해 계통 분류학적 측면에서 상당한 어려움을 안겨준다. 170년에 걸친 분류학적 논쟁에도 불구하고, 널리 받아들여지는 7개 아속, 38개 절(section) 체계(Richards, 2003)는 게놈 규모의 데이터에 의해 거의 검증되지 않은 상태로 남아 있다. 이전의 분자 연구들은 희소한 분류군 샘플링과 소수의 마커(예: rbcL, matK, ITS)에 의존했기에, 밀접하게 연관된 계통들을 해결하는 데 필요한 계통 발생적 신호가 부족했다. 또한, 이형꽃술(heterostyly, 타가수분 촉진)에서 동형꽃술(homostyly, 자가수분)로의 전이 과정에 대한 진화적 동인과 칭하이-티베트 고원(QTP) 상승과 관련된 본 속의 생물지리학적 역사는 견고한 계통 발생적 맥락 내에서 여전히 명확히 밝혀지지 않았다.

연구 방법론
본 연구는 현재까지 수행된 Primula에 대한 가장 포괄적인 엽록체 계통 유전체 분석을 제시한다.

• 분류군 샘플링: 분석에는 106개 종을 대표하는 180개의 개체가 포함되었으며, 이는 7개 아속 중 6개와 38개 절 중 22개를 포괄한다. 여기에는 중국 운남과 티베트에서 새로 시퀀싱된 18개의 샘플과 GenBank에서 검색된 162개의 서열이 포함되었다. 네 종의 Androsace가 외집단(outgroup)으로 사용되었다.
• 게놈 시퀀싱 및 조립: 게 genomic DNA를 추출하고, short-insert 라이브러리를 Illumina HiSeq X Ten 플랫폼을 통해 시퀀싱하였다. 엽록체 게놈(plastomes)은 GetOrganelle을 사용하여 조립되었으며 Geneious Prime에서 주석을 달았다.
• 비교 게놈학: 29개의 대표 종을 대상으로 rps19, rpl2, ndhF, ycf1을 포함한 역반복(inverted repeat, IR) 경계의 구조적 특징을 분석하였다. 뉴클레오타이드 다양성($\pi$)은 슬라이딩 윈도우 방식(1,000 bp 윈도우, 500 bp 스텝)을 사용하여 초변이 영역을 식별하는 데 계산되었다.
• 계통 발생 재구성: 184개의 엽록체 게놈 서열(Primula 180개 + 외집단 4개)을 사용하여 최적 적합 치환 모델(ModelFinder를 통해 선택된 TVM+F+R2) 기반의 최대 우도법(Maximum Likelihood, ML) 트리를 구축하였다. 가지 지지도는 Ultrafast Bootstrap(UFBoot)과 SH-유사 근사 우도비 검정(SH-aLRT)을 통해 평가되었다.
• 분기 시간 및 생물지리학: 분기 시간은 두 개의 화석 보정 지점(Primula–Androsace 분기점(마이오세 종자 화석) 및 Primula의 크라운 그룹(P. riosacea 화석))을 사용하여 BEAST v2.5.2로 추정하였다. 조상 영역 재구성은 RASP의 S-DEC 모델을 사용하여 수행되었으며, 꽃술 유형(이형꽃술 vs. 동형꽃술)에 대한 조상 상태 재구성은 Mk1 모델을 활용하였다.

주요 결과

• 엽록체 게놈 구조: Primula의 엽록체 게놈은 115개의 뚜렷한 유전자를 가진 매우 보존된 사분 구조(quadripartite structure, 150,164–154,897 bp)를 보인다. 전체적인 구조는 안정적이지만, IR 경계 분석 결과 rps19, rpl2, ndhF, ycf1을 포함하는 특정 확장/수축이 나타났다. 슬라이딩 윈도우 분석을 통해 ndhF–rpl32와 ycf1이 초변이 영역($\pi > 0.068$)임을 확인하였으며, 이는 표준 바코드보다 우수한 종 식별 능력을 제공한다.
• 계통 발생 체계: 본 속은 단계통으로 확인되었으며, 세 개의 강력하게 지지되는 주요 클레이드(P1, P2, P3)와 하나의 독특한 계통(P. palinuri)으로 구분된다.
  • 클레이드 P1: Subg. Pinnatae, Subg. Auriculastrum (특히 Sect. Amethystina 포함), 그리고 Subg. Aleuritia의 4개 절을 포함한다.
  • 클레이드 P2: Subg. Primula, Subg. Auriculastrum (Sect. Auricula 포함), 그리고 Subg. Aleuritia의 8개 절을 포함한다.
  • 클레이드 P3: Subg. Auganthus와 Subg. Carolinella를 포함한다.
• 분류학적 함의: 본 연구는 전통적 분류 체계에서 광범위한 비단계통성(non-monophyly)을 발견하였다. 오직 Subg. Pinnatae와 Subg. Primula만이 단계통이었다. 저자들은 다음과 같은 구체적인 개정을 제안한다:
  • Sect. Amethystina는 Subg. Aleuritia 내에 위치하므로 그곳에 머물러야 한다.
  • Subg. Pinnatae는 Subg. Aleuritia 내에 위치하므로 절(sectional) 단계로 격하되어야 한다.
  • Sect. Carolinella는 Subg. Auganthus 내에 위치하며, 그 정의적 형질(circumscissile capsules)이 수렴 진화된 것이므로 절 단계로 격하되어야 한다.
• 분기 및 생물지리학: Primula의 크라운 그룹은 초기 마이오세(~19.9 Mya)에 기원하였다. 주요 다양화는 칭하이-티베트 고원(QTP)의 상승 및 지구 냉각기와 일치하는 플라이오세(3–5 Mya)에 발생하였다. 조상 영역 재구성은 횡단산맥-동부 히말라야, 중국 중부/동부, 그리고 유럽의 중첩 지역에서 기원했음을 시사하며, 이후 산맥 형성에 의해 주도된 분산 사건들이 일어났다.
• 형질 진화: 이형꽃술이 조상 상태임이 확인되었다. 동형꽃술은 고산 지대(2,400–5,000 m)에서만 독점적으로, 그리고 QTP 상승기인 3–5 Mya 기간 동안 최소 8회 독립적으로 진화하였다.

의의 및 주장
저자들은 본 연구가 Primula 속의 계통 분류학적 논란을 해결하는 견고한 게놈 규모의 계통 발생 체계를 제공한다고 주장한다. 본 속의 다양화 맥동을 마이오세-플라이오세의 지질학적 사건(특히 QTP 상승)과 기후 변화에 연결함으로써, 이 연구는 고산 식물의 방사 진화에 대한 새로운 관점을 제시한다. 또한, 동형꽃술의 반복적이고 독립적인 진화는 매개체가 부족한 고산 환경에 대한 명확한 적응 반응으로 제시되며, 이는 "생식 보장(reproductive assurance)" 가설을 뒷받침한다. 본 연구는 향후 분류학적 개정, 보전 계획 및 진화 연구를 위한 토대를 마련하며, 이 빠르게 방사 진화하는 속의 복잡한 진화 역사를 완전히 해결하기 위해서는 핵 게놈 데이터의 통합이 필요함을 인정하고 있다.