ultivar-resolved assembly and characterization of the mitochondrial genome of Medicago sativa cv. AH reveal repeat-mediated alternative junctions and extensive RNA editing

기술 요약: Medicago sativa cv. AH 미토콘드리아 게놈의 품종별 해상도 조립 및 특성 분석

문제 제기
식물 미토콘드리아 게놈은 반복 서열에 의한 재조합, 광범위한 RNA 편집, 그리고 엽록체-미토콘드리아 DNA 전이(MTPT)로 인해 높은 구조적 가소성을 특징으로 한다. Medicago sativa(알팔파)는 중요한 사료 작물이지만, 기존의 미토콘드리아 게놈 자원은 구조적 검증, RNA 편집 지형, 그리고 계통학적 맥락을 통합하는 품종별 해상도의 조립본이 부족하다. 특히, 특정 품종의 특정 긴 직립 반복 서열(long direct repeats)이 실험적으로 검증 가능한 대안적 게놈 접합부(alternative genome junctions)를 생성하는 정도와, 이러한 구조적 변이가 RNA 편집과 같은 전사체 수준의 특징과 어떻게 연관되는지는 여전히 충분히 해결되지 않은 상태이다.

방법론
본 연구는 롱리드 시퀀싱(long-read sequencing), 그래프 기반 조립, 그리고 직교 실험 검증을 결합한 종합적인 워크플로우를 채택하였다:

• 시퀀싱 및 조립: M. sativa cv. AH의 게놈 DNA를 PacBio HiFi(Circular Consensus Sequencing)를 사용하여 시퀀싱하여 고정밀 리드를 생성하였다. 이를 Oatk(복잡한 식물 세포 소기관 게놈을 위해 설계된 도구)를 사용하여 조립하였으며, 그 결과 그래피컬 단편 조립(Graphical Fragment Assembly, GFA)을 얻었다.
• 그래프 해상도 확보: 조립 그래프를 Bandage를 사용하여 분석하였다. 반복 서열로 인해 발생하는 분기를 해결하기 위해 리드 지원 연결(read-supported connections)을 우선순위에 두어, 그래프를 하나의 지배적인 원형 경로로 단순화하였다.
• 주석 달기 및 분석: PMGA, tRNAscan-SE, BLASTn을 사용하여 게놈 주석을 달았다. 코돈 사용 빈도는 MEGA12를 통해 분석하였다. 반복 서열 지형(SSR, 탄뎀 및 분산 반복 서열)은 MISA-web, Tandem Repeats Finder, ROUSFinder를 사용하여 매핑하였다.
• 구조의 실험적 검증: 반복 서벨에 의한 재조합을 검증하기 위해, 후보 긴 직립 반복 서열의 고유한 플랭킹 영역을 타겟팅하는 4가지 프라이머 조합을 사용하여 접합부 PCR(junction PCR)을 수행하였다. 생성된 산물은 생더(Sanger) 시퀀싱을 통해 대안적 접합부를 확인하는 데 사용되었다.
• MTPT 및 계통학: BLASTn을 사용하여 엽록체 게놈에 대한 엽록체-미토콘드리아 전이(MTPT)를 식별하였다. 계통학적 배치는 23개의 보존된 단백질 코딩 유전자를 바탕으로 최대 우도 추론(IQ-TREE 2)을 통해 결정되었다. 신테니(synteny)는 MCScanX로 분석하고 TBtools로 시각화하였다.
• RNA 편집: Deepred-Mt(확률 컷오프 0.9)를 사용하여 C-to-U 편집 사이트를 예측하였다. 대표적인 사이트들은 gDNA와 cDNA 크로마토그램을 비교하는 RT-PCR 및 생더 시퀀싱을 통해 검증되었다.

주요 결과

• 게놈 조립: M. sativa cv. AH의 완전한 미토콘드리아 게놈은 300,797 bp 크기의 원형 분자로 조립되었으며, GC 함량은 45.36%였다. 이는 33개의 고유한 단백질 코딩 유전자(PCG), 16개의 tRNA, 3개의 rRNA를 인코딩한다.
• 반복 서열 지형 및 구조적 역동성: 게놈은 59개의 단순 반복 서열(SSR), 9개의 탄뎀 반복 서열, 112쌍의 분산 반복 서열을 포함한다. 특정 4,354 bp 직립 반복 서열(R1)이 재조합 핫스팟으로 식별되었다.
  • 검증: 접합부 PCR과 생더 시퀀싱은 이 반복 서열이 상동 재조합을 매개하여 "마스터 서클(master-circle)" 구성과 대안적 서브게놈 분자를 모두 생성함을 확인하였다. 마스터 서클 접합부가 더 강하게 증폭되었는데, 이는 지배적인 이소형(isoform)이 존재함을 시사하며, 재조합된 접합부는 낮은 빈도로 존재하였다.
• MTPT: 게놈의 0.24%(725 bp)를 차지하는 6개의 엽록체 유래 단편이 검출되었다. 여기에는 5개의 온전한 엽록체 유래 tRNA 유전자(trnD-GUC, trnH-GUG, trnM-CAU, trnN-GUU, trnW-CCA)가 포함된다.
• 계통 및 신테니: 계통 분석 결과, cv. AH는 다른 M. sativa 조립본들과 함께 높은 부트스트랩 값(bootstrap = 96)을 보이며 M. sativa 클레이드 내에 배치되었다. 그러나 신테니 분석 결과, 다른 M. sativa 품종 및 아종과 비교했을 때 유전자 블록의 순서와 방향에서 상당한 재배열이 관찰되었으며, 이는 밀접한 계통학적 관계에도 불구하고 현저한 종내 구조적 변이가 존재함을 보여준다.
• RNA 편집: 33개의 PCG 전체에서 총 468개의 고신뢰도 C-to-U RNA 편집 사이트가 예측되었다.
  • 특징: 편집의 95.73%가 비동의(nonsynonymous) 편집이었다. 가장 빈번한 아미노산 전환은 Ser→Leu (107회) 및 Pro→Leu (104회)였다. nad4 유전자가 가장 많은 사이트(45개)를 보유하였다.
  • 검증: 실험적 검증을 통해 nad4L의 정형 AUG 개시 코돈 복원(ACG-to-AUG 편집을 통해) 및 다수의 비동의 전환을 포함한 편집 사건들이 확인되었다.

의의 및 주장
저자들은 본 연구가 구조적 이소형 검증, MTPT 특성 분석, 계통 배치, 그리고 실험적으로 뒷받침된 RNA 편집 지형을 단일 프레임워크 내에 통합한 M. sativa cv. AH의 품종별 해상도 미토콘드리아 게놈 자원을 제공한다고 주장한다.

본 논문은 다음과 같은 점을 통해 알팔파의 미토콘드리아 게놈 역동성에 대한 이해를 진전시킨다고 단언한다:

1. (접합부 PCR을 통해) 알팔파 품종에서 긴 반복 서열에 의한 대안적 접합부에 대한 직접적인 실험적 증거를 제공하였다.
2. M. sativa 품종들이 보존된 계통학적 신호를 공유하면서도, 게놈 구조 측면에서는 뚜렷한 종내 재배열을 보인다는 것을 입증하였다.
3. 표적 검증(정형 번역 신호의 기능적 복원 포함)을 통해 뒷받침되는 품종 특이적 RNA 편집 지형을 전달하였다.

저자들은 이 자원이 향후 세포질 다양성 및 육종 응용 연구를 위한 토대가 될 것으로 보고 있으며, 특히 미토콘드리아 변이가 세포질 수컷 불임(CMS)에 갖는 관련성을 언급하였다. 다만, 특정 구성을 농업적 표현형과 연결하기 위해서는 추가적인 전용 표현형 분석 및 유전적 테스트가 필요함을 명시함으로써 연구의 범위를 겸허히 제한하였다. 본 연구는 단일 품종 분석 및 이소형의 양적 검증 부족과 같은 한계를 인정하며 다소 제한적인 범위 내에서 수행되었다.