Does climate shape plant immune gene repertoires? A phylogenetically controlled test across angiosperms

기술 요약: 기후가 식물 면역 유전자 레퍼토리를 형성하는가? 피사계통학적으로 통제된 속씨식물 전반에 걸친 테스트

문제 제기
속씨식물(angiosperms)의 면역 수용체 유전자 수는 종마다 수십 개에서 천 개 이상의 NLR(Nucleotide-Binding Leucine-Rich Repeat) 및 PRR(Pattern Recognition Receptor)에 이르기까지 한 자릿수 이상의 차이를 보이며 매우 다양하다. 따뜻하고 습하며 저위도인 환경에 처한 식물들이 더 높은 병원균 압력을 경험하고, 결과적으로 더 큰 면역 레퍼토리를 유지할 것이라는 광범위한 거시생태학적 기대가 존재한다. 그러나 이러한 경향성에 대한 실증적 증거는 엇갈린다. 결정적인 혼란 요인은 면역 유전자 수와 기후적 니치(niche) 모두 계통 발생적으로 상속된다는 점이다. 즉, 밀접하게 연관된 종들은 유사한 기후와 유사한 유전자 수를 공유하는 경향이 있다. 명시적인 계통 발생적 교정 없이 관찰된 상관관계는 환경에 대한 적응적 반응이라기보다 공유된 조상의 특성을 반영하는 것일 수 있다. 또한, 식물의 기능적 형질, 특히 목본(나무/관목)과 초본 성장 형태 사이의 구분이 면역 투자 형성에 미치는 역할은 거시생교학적 규모에서 아직 충분히 탐구되지 않았다.

연구 방법론
본 연구는 최근 발표된 808종의 속씨식물을 대상으로 균일하게 주석(annotation)이 달린 면역 수용체 유전자 카탈로그를 활용하였다. 분석은 6개의 NLR 패밀리, 8개의 상위 수준 RLK(Receptor-Like Kinase) 패밀리, 21개의 LRR-RLK 하위 클래스, 4개의 RLP(Receptor-Like Protein) 패밀리를 포함한 39개의 주석 처리된 면역 수용체 패밀리에 집중하였다.

• 데이터 필터링: 특수 생활 양식/서식지(SLH) 종(예: 수생 식물, 기생 식물, 식충 식물)은 면역 레퍼토리의 수렴적 감소가 거시생태학적 패턴을 왜곡할 수 있으므로 주요 분석에서 제외하였다. 그 결과 기후 분석에는 584종, 형질 분석에는 176종의 데이터셋이 사용되었다.
• 예측 변수:
  • 기후/지리: WorldClim v2.1에서 유도된 18개 변수(19개 생물 기후 변수의 중앙값 및 표준 편차)와 지리적 좌표(위도, 범위).
  • 기능적 형질: GIFT의 범주형 형질(목질화, 성장 형태)과 TRY 데이터베이스의 연속형 형질(잎/줄기 특성).
  • 병원균 노출: GBIF 출현 기록과 GlobalFungi 토양 메타바코딩 데이터를 통해 유도된 주변 환경 곰팡이 병원균 압력의 독립적 추정치.
• 통계적 프레임워크:
  • 이변량 스크리닝: 모든 수용체 패밀리와 예측 변수에 대해 Spearman 순위 상관관계(연속형) 및 Kruskal-Wallis 검정(범주형)을 수행한 후, 공동 Benjamini-Hochberg(FDR) 교정을 적용하였다.
  • 계통 발생적 교정: 유의미한 연관성은 Pagel's $\lambda$를 최대 우도법으로 추정한 PGLS(Phylogenetic Generalized Least Squares) 모델을 사용하여 재평가하였다. 게놈 조립 품질을 제어하기 위해 BUSCO 완결성을 공변량으로 포함하였다.
  • 다중 검정 전략: 두 단계의 교정 접근법을 채택하였다. 사전 지정된 핵심 패밀리인 LysM-RLK를 포함한 12개 주요 수용체 패 family 내에서 확인적 테스트를 수행하였고, 전체 39개 패밀리를 대상으로 탐색적 스윕(sweep)을 실시하였다.
  • 민감도 분석: 모델을 브랜치 길이 변환, 분류군 서브샘플링, 가축화 상태, 게놈 크기 및 조상 상태 재구성에 대해 테스트하였다.

주요 결과

1. 기후 및 지리:

   • 이변량 스크리닝 결과 702개의 테스트 중 56개의 유의미한 상관관계가 식별되었다. 그러나 모든 효과 크기는 완만했다($|\rho| \le 0.24$).
   • 결정적으로, 그 어떤 기후-수용체 연관성도 계통 발생 교정(PGLS)을 통과하지 못했다. 예를 들어, 연평균 기온과 TNL 유전자 사이의 상관관계($\rho = -0.237$)는 계통 발생을 고려한 후 비유의적($p = 0.51$)이 되었다.
   • 가장 강한 원시 상관관계의 방향은 가설과 상충하는 경우가 많았다. 즉, 따뜻한 기후는 더 많은 헬퍼-NLR(TNL, NRG1, RNL)이 아니라 더 적은 수의 헬퍼-NLR과 연관되어 있었다.
   • CNL 유전자와 위도 사이의 음의 상관관계는 단자엽 식물(monocots, TNL 유전자가 결핍됨)의 저위도 과잉 대표성에 의한 분류학적 오류(Simpson's paradox)로 확인되었다.

2. 기능적 형질 (성장 형태):

   • 기후와 대조적으로, 성장 형태는 면역 유전자 함량과 더 일관된 연관성을 보였다.
   • LysM-RLK: 목본 및 나무 종은 비목본 및 초본 종보다 더 많은 LysM-RLK 유전자를 보유하고 있었다. 이 신호는 비모수 검정(FDR $p = 0.004$)과 12개 주요 패밀리에 대한 PGLS 모델($p = 0.007$) 모두에서 견고하게 나타났다.
   • 기타 수용체: Malectin-RLK와 여러 LRR-RLK 하위 클래스 역시 비모수 검정에서 목본/나무 종에서 높은 수치를 보였다.
   • 통계적 유의성: LysM-RLK 신호는 혼란 변수(브랜치 길이, 가축화, 게놈 크기)에 대해 견고했으나, 39개 전체 패밀리에 대한 엄격한 공동 FDR 교정은 통과하지 못했다(FDR $p = 0.13$). 저자들은 이를 실제 존재하지만 미미한 효과이며, 사전 지정된 가설로서는 통계적으로 유의미하지만 게놈 전체 규모에서는 그렇지 않다고 해석한다.

3. 병원균 노출:

   • GBF 및 GlobalFungi를 통한 주변 곰팡이 병원균 노출 분석 결과, 목본 종이 초본 종보다 더 큰 병원균 압력을 경험하지 않는 것으로 나타났다. 일부 지표에서는 목본 종의 노출도가 오히려 더 낮았다.
   • PGLS 모델에 병원균 노출 지표를 추가해도 성장 형태 계수가 변하지 않았는데, 이는 목본성과 LysM-RLK 사이의 연관성이 주변 환경 곰팡이 병원균에 대한 차별적 노출에 의해 매개되는 것이 아님을 나타낸다.

의의 및 주장
본 논문은 계통 발생적 구조를 고려할 때, 기후와 지리는 속씨식물의 면역 수용체 유전자 함량을 예측하는 견고한 변수가 아님을 결론짓는다. 관찰된 면역 유전자 수의 거시생태학적 경향은 직접적인 환경적 병원균 압력에 대한 적응적 반응이라기보다 주로 계통 발생적 구조에 의한 인공적 결과(artifact)이다.

식별된 가장 일관된 생태적 상관관계는 성장 형태이다. 목본 및 나무 종은 특히 LysM-RLK 및 Malectin-RLK와 같은 특정 세포 표면 수용체의 확장된 레퍼토리를 보유하는 경향이 있다. 그러나 본 연구는 집합적인 유전자 수가 병원균 면역에 대한 투자 증가를 반영하는지, 아니면 공생 신호 전달(균근/리조비아)에 대한 것인지를 명확히 구분할 수 없다고 명시한다. 이는 LysM 수용체가 두 경로 모두에서 기능하기 때문이다. 저자들은 하위 클래스 수준의 주석(LYK vs. LYR)이나 표적 샘플링이 필요하며, 이는 현재 데이터셋의 범위를 벗어난다고 언급하였다.

이 연구는 계통 발생적 통제가 엄격하지 않을 경우, 유전체 특성과 환경 변수 간의 상관관계가 얼마나 오해를 불러일으킬 수 있는지를 보여주는 거시생태학적 유전체학의 경고 사례로서 의미를 갖는다. 성장 형태가 기후보다 더 명확한 예측 변수임을 보여주지만, 효과 크기는 미미하며, 관찰된 패턴의 기계론적 기초는 여전히 미해결 과제로 남아 있다.