Genome-wide identification and comparative analysis of endogenous retroviruses across five breed-representative pig assemblies

기술 요약: 5종의 품종 대표 돼지 어셈블리를 통한 전장 게놈 규모의 내인성 레트로바이러스 식별 및 비교 분석

문제 제기
돼지 이종 이식은 인간의 장기 부족 문제에 대한 잠재적 해결책이지만, 내인성 레트로바이러스(ERV)의 종간 전파 위험은 여전히 중요한 안전 장벽으로 남아 있다. 현재의 위험 평가는 주로 단일 참조 게놈에 의존하고 있어, 다양한 유전적 배경에 따른 품종별 ERV 다양성, 구조적 변이 및 삽입 연령에 관한 "사각지대"를 생성한다. 기존 연구들은 주로 생물 의학적 안전성(예: PERV 불활성화)이나 특정 레트로트랜스포존 조사에 집중해 왔으며, 돼지의 ERV 풍부도, 구조적 완전성 및 게놈 맥락에 대한 통합적이고 다중 어셈블리 기반의 기술이 부족했다. 본 연구는 분석 프레임워크를 단일 참조 모델에서 비교 가능한 다중 품종 경관(multi-breed landscape)으로 전환함으로써 품종별 변이 문제를 해결하고자 한다.

연구 방법론
연구진은 서로 다른 유전적 배경을 나타내는 5종의 염색체 수준 돼지 게놈 어셈블리(Bama 미니 돼지, Duroc, Meishan, Ningxiang, USMARC)를 대상으로 비교 게놈 분석을 수행하였다.

• 식별 및 주석 달기: LTRharvest/LTRdigest 파이프라인을 사용하여 ERV 로커스(loci)를 식별하였다. 후보군들은 쌍을 이룬 긴 말단 반복(LTR) 경계(길이 100–2,000 bp, 유사도 80% 이상)와 표적 삽입 부위 중복(TSD, 4–6 bp)을 기준으로 필터링되었다.
• 구조적 분류: 쌍을 이룬 LTR, 프라이머 결합 부위(PBS), 폴리퓨린 트랙트(PPT) 및 내부 레트로바이러스 도메인(예: gag, pol, env, RVT_1)의 보존 여부에 따라 로커스를 완전(complete), 준완전(near-complete) 또는 **부분(partial)**으로 분류하였다.
• 계통 발생 및 연령 분석: 완전한 형태의 ERV들을 정렬(MAFFT)하고, 24종의 레트로바이러스 참조 서열 패널을 대상으로 최대 우도 계통 발생 추론(IQ-TREE 2)을 실시하였다. 상대적인 삽입 연령은 쌍을 이룬 LTR 간의 Kimura 2-parameter 발산($T = K / 2r$)을 사용하여 추정하였다.
• 기능 및 게놈 맥락: ERV 관련 유전자(10 kb 이내 인접 영역)를 분석하여 유전자 온톨로지(GO) 및 교토 유전자 및 경로(KEGG) 농축도를 확인하였다. 조직별 발현 프로파일은 GSE93855 데이터셋(6개 조직)에서 도출하였으며, 양적 형질 위치(QTL)와의 공국소화(colocalization)를 평가하였다.

주요 결과

• ERV 경관 및 구조: 5개 품종 전체에서 총 23,767개의 ERV 로커스가 식별되었다. 이 경관은 구조적으로 퇴화된 요소들이 지배하고 있다: 98.33%가 부분적(partial)이었고, 1.35%가 준완전(near-complete)이었으며, 오직 0.32%(75개 로커스)만이 완전한(complete) 상태로 분류되었다.
• 계통 발생적 다양성: 완전한 ERV들은 네 가지 레트로바이러스 클레이드(Betaretrovirus, Gammaretrovirus, Spumaretrovirus, Epsilonretrovirus)로 군집을 이루었다. 품종별 차이가 관찰되었는데, 예를 들어 Betaretrovirus와 Spumaretrovirus 후보는 Ningxiang에서 더 빈번하게 나타난 반면, Gammaretrovirus 후보는 Duroc와 USMARC에서 더 흔하게 발견되었다.
• 삽입 연령의 이례적 패턴: 상대적 삽입 연령 추정 결과 역설적인 패턴이 드러났다. 작은 유효 집단 크기($N_e$)를 가진 중국 토착 품종인 Ningxiang이 가장 낮은 중앙값 삽입 연령(0.81 Myr)을 보인 반면, 복합 외래 혈통인 USMARC는 가장 높은 연령(7.85 Myr)을 보였다. 이는 작은 유효 집단 크기에서 정화 선택(purifying selection)의 효율성이 감소함에 따라, 오래된 삽입 기록이 반드시 다양한 상업적 혈통에만 국한되는 것이 아니라 오히려 더 젊은 완전한 요소들이 지속될 수 있음을 시사한다.
• 게놈 맥락 및 기능:
  • 공유 유전자: 5개 어셈블리 모두에서 278개의 ERV 관련 유전자가 공유되었다.
  • 농축도: 이러한 공유 유전자들은 시냅스 신호 전달, 소낭 수송 및 세포 접합 경로(GO)와 글루타메이트성 시냅스 및 옥시토신 신호 전달(KEKG)과 같은 경로에 농축되어 있었다.
  • 발현: 공유 유전자 세트는 조직 특이적 발현을 보였으며, 비장, 신장 또는 폐에서 정점 발현이 자주 관찰되었다.
  • QTL 공국소화: ERV 관련 유전자와 육질, 성장, 비육 및 번식에 관한 주석 처리된 QTL 사이의 상당한 중첩이 확인되었다. 특히, 278개의 공유 유전자 중 56.5%가 QTL 구간과 중첩되었는데, 이는 품종 특이적 세트보다 높은 비율이다.

의의 및 주장
저자들은 본 연구를 돼지 ERV에 대한 포괄적이고 품종별로 해상도가 높은 범-게놈(pan-genome) 프레임워크의 구축으로 규정한다. 본 연구의 의의는 다음과 같다:

1. PERV 중심적 관점의 탈피: 본 연구는 고전적인 PERV(Porcine Endogenous Retrovirus) 프레임워크를 넘어 이해를 확장하며, 다양한 품종에 걸친 더 넓은 범위의 ERV 후보군과 그 구조적 상태에 대한 목록을 제공한다.
2. 인구 통계학적 통찰: 삽입 연령 프로파일에 관한 연구 결과는 품종의 인구 통계학적 역사(창시자 효과, 선택 강도, 혼혈)가 검출 가능한 ERV 구성(특히 완전한 요소의 유지에 미치는 유효 집단 크기의 영향)을 어떻게 형성하는지에 대한 새로운 관점을 제공한다.
3. 조절 잠재력: ERV 관련 유전자가 경제적으로 중요한 형질(육질, 번식)에 관한 QTL과 공국소화되고 시냅스/수송 경로에 농축되어 있다는 점은, 통합된 ERV 로커스가 단순한 게놈 잔해가 아니라 기능적 조절 요소로서 작용할 수 있음을 시사한다.
4. 실용적 적용: 생성된 카탈로그는 이종 이식 맥락에서 기증자 계통 선택 및 정밀한 CRISPR 매개 바이러스 불활성화를 위한 실질적인 벤치마크와 타겟 리스트를 제공한다.

한계 및 겸허한 태도
본 논문은 인과관계에 대해 명확히 겸허한 입장을 유지한다. 저자들은 다음과 같이 밝힌다:

• 본 분석은 참조 어셈블리에 기반하며, 집단 수준의 삽입 다형성(다수의 개체 샘플링 필요)을 포착하지 않는다.
• 탐지 방법(쌍을 이룬 LTR)은 젊고 덜 분화된 후보를 선호하는 결정적 편향(ascertainment bias)을 생성하므로, 연령 추정치는 절대적인 것이 아니라 상대적이다.
• QTL과의 공국소화 및 기능적 농축은 후보 게놈 맥락을 나타내는 것이지, 특정 ERV 로커스가 인접한 유전자나 표현형을 조절한다는 직접적인 증거는 아니다. 기능적 활성을 확립하기 위해서는 후성유전체 프로파일링이나 표적 섭동과 같은 실험적 검증이 필수적이다.