Distinct Structuring Mechanisms of Extrinsic and Intrinsic Factors on the Wild Marine Fish Gut Microbiome

기술 요약: 야생 해양 어류 장내 미생물 군집에 대한 외적 및 내적 요인의 구별된 구조화 메커니즘

문제 제기
어류의 장내 미생물 군집(GM)은 영양분 흡수, 면역 항상성 및 숙주 생리에 매우 중요하다. 외적 요인(예: 계절, 먹이, 수온)과 내적 요인(예: 숙주 종, 유전) 모두가 GM을 형성한다는 점은 확립되어 있으나, 이러한 요인들이 GM 구조를 유도하는 구체적인 메커니즘은 여전히 불분명하다. 기존 연구들은 주로 단일 지표(예: 분류학적 조성 또는 알파/베타 다양성)에 의존하여 특정 요인이 GM에 영향을 미친다는 것을 입증해 왔다. 그러나 이러한 접근 방식은 요인들이 유사하게 작용하는지, 아니면 서로 다르게 작용하는지—즉, 이들이 주로 계통적 유사성 없이 미생물 풍부도(abundance)를 재편하는 것인지, 아니면 뚜렷한 조성 변화 없이 계통적 유사성을 형성하는 것인지—를 구분하지 못한다. 본 연구는 외적 및 내적 요인이 서로 다른 과정을 통해 어류 GM을 어떻게 차별적으로 구조화하는지에 대한 이해의 간극을 메우고자 한다.

연구 방법론
연구진은 한국 남해안에서 채집한 21종의 야생 해양 어류의 GM 특성을 분석하였다 (2021년 여름과 2022년 겨울 두 계절).

• 시료 채집: 21종으로부터 118개의 장 시료를 얻었다. 물고기를 마취시킨 후 후장(hindgut) 내용물을 추출하여 -80°C에 보관하였다.
• 시퀀싱: Illumina MiSeq 플랫폼을 사용하여 16S rRNA V3-V4 앰플리콘 시퀀싱을 수행하였다.
• 데이터 처리: QIQI2 (v2024.05)를 사용하여 시퀀스를 처리하였으며, DADA2를 통해 데노이징(denoising) 및 키메라 제거를 수행하였다. 분류학적 분류는 SILVA 데이터베이스 (v138.1)를 기준으로 수행되었다.
• 통계 분석:
  • 다양성 지표: 알파 다양성(Shannon, Simpson, Faith's PD)과 베타 다양성을 계산하였다.
  • 비유사도 지표: 군집 구조의 다양한 측면을 분석하기 위해 세 가지 지표를 사용하였다: 분류군 풍부도를 위한 Bray-Curtis (BC), 계통적 존재/부재를 위한 Unweighted UniFrac, 그리고 계통적 풍부도를 위한 Weighted UniFrac.
  • 동인(Drivers): 연구는 계절(외적), 영양 단계(먹이의 대리 지표로서의 외적 요인), 그리고 어류 종(내적)의 영향을 분석하였다.
  • 기능 예측: PICRUSt2를 사용하여 기능적 경로(KEGG orthologs)를 예측하였다. Random Forest 분류기(5-fold 교차 검증, 10회 반복)를 사용하여 클래스 불균형을 처리하기 위한 AUPRC(Area Under the Precision-Recall Curve)를 통해 이러한 요인들이 기능적 프로파일에 미치는 예측력을 평가하였다.
  • 핵심 ASV: 숙주 특이적 필터링을 평가하기 위해 어류 종별 핵심 앰플리콘 서열 변이(core ASVs)를 식별하였다.

주요 결과

1. 외적 요인은 풍부도 변화를 유도함 (Bray-Curtis):

   • 계절: 계절은 분류군 풍부도(BC 비유사도)를 기준으로 GM 조성에 유의미한 영향을 미쳤다. Catenococcus와 Enterovibrio는 여름의 지표 분류군이었으며, Aliivibrio와 Shewanella는 겨울의 특징을 나타냈다. 주목할 점은 Vibrio와 Photobacterium 같은 속(genera)은 두 계절 모두 존재했지만, 해당 속 내의 특정 ASV들은 뚜렷한 계절적 농축을 보였다는 것이다.
   • 영양 단계: 영양 단계 또한 GM 풍부도에 유의미한 영향을 미쳤다. Photobacterium은 높은 영양 단계(>3.4)의 지표였으며, Vibrio는 낮은 영양 단계(<3.4)와 연관되었다. 맨텔 테스트(Mantel test)는 영양 단계(연속 변수)와 BC 비유사도 사이의 유의미한 상관관계를 확인하였다.
   • 다양성: 계절이나 영양 단계 모두 알파 다양성(Shannon/Simpson)을 유의미하게 변화시키지 않았는데, 이는 이러한 요인들이 전체적인 다양성 풍부도를 변화시키기보다는 기존의 군집을 재편함을 의미한다.

2. 내적 요인은 계통적 유사성을 유도함 (Unweighted UniFrac):

   • 어류 종: 풍부도 기반의 차원 축소(PCoA)에서는 뚜렷한 종별 클러스터링이 나타나지 않았으나, Unweighted UniFrac(계통적 거리)을 통해 분석했을 때 어류 종 정체성은 유의미한 특이성을 보였다.
   • 핵심 ASV: "어류 종-핵심 ASV" 분석(숙주 종 내에서 지속적으로 검출되는 ASV) 결과, 숙주 종 내에서 계통적 클러스터링이 나타났다. 예를 들어, 특정 종의 핵심 ASV는 Vibrio/Catenococcus 클레이드 내에 클러스터링된 반면, 다른 종들은 Enterovibrio 내에 클러스터링되었다.
   • 메커니즘: 이는 내적 요인인 숙주 특이적 필터링이 상대적 풍부도와 관계없이 계통적으로 유사한 미생물을 선택한다는 것을 시사하며, 이러한 패턴은 계절적 변화가 어떤 특정 ASV가 "핵심"으로 검출되는지를 변화시키더라도 보존된 상태로 유지되었다.

3. 계절에 따른 기능적 층화:

   • 예측력: Random Forest 분류기는 계절이 기능적 경로 프로파일의 가장 강력한 예측 인자(AUPRC = 0.8571)임을 입증한 반면, 영양 범주는 더 약한 예측 인자(AUPRC = 0.5791)였다. 어류 종은 샘플 크기 제한으로 인해 평가할 수 없었다.
   • 기능적 변화:
     • 여름: 탄수화물 소화/흡수, 스티렌 분해, 그리고 다환 방향족 탄화수소(PAH) 분해와 관련된 경로가 풍로하였다. 항균 경로(포스포네이트 대사, 시데로포어 생합성) 또한 여름철 병원균 증식과 맞물려 풍로하였다.
     • 겨겨울: 단백질 소화/흡수, 아트라진 분해(잠재적으로 농업 유출수와 연관됨), 그리고 스테로이드 생합성이 풍로하였다.

주요 기여

• 메커니즘의 차별화: 본 연구는 외적 요인(계절, 영양 단계)과 내적 요인(어류 종)이 근본적으로 다른 메커니즘을 통해 어류 GM을 구조화한다는 증거를 제공한다. 외적 요인은 주로 미생물의 풍부도 변화(Bray-Curtis로 포착)를 유도하는 반면, 내적 요인은 계통적 유사성과 숙주 특이적 필터링(Unweighted UniFrac 및 핵심 ASV 분석으로 포착)을 유도한다.
• 분류학적 신호의 해상도: 본 연구는 특정 속(예: Vibrio, Photobacterium)에 대한 계절적 및 영양 단계적 신호가 속(genus) 수준이 아닌 ASV 수준에서만 해상도가 확보된다는 점을 강조하며, 고해상도 시퀀싱의 필요성을 역설한다.
• 기능적 역동성: 환경적 변화가 어류 GM의 예측된 기능적 능력(특히 탄수화물 대 단백질 소화 및 방어 기제와 관련된 기능)을 형성하는 주요 동력임을 보여준다.

의의 및 주장
저자들은 자신들의 연구 결과가 외적 및 내적 요인이 서로 다른 과정을 통해 어류 GM을 차별적으로 구조화한다는 것을 밝혀냈다고 주장한다. 구체적으로 다음과 같이 결론짓는다:

1. 숙주 특이적 필터링은 어류 종 내에서 GM의 계통적 구조를 주도하며, 계절적 변동이 특정 미생물 조성을 변화시키더라도 보존된 진화적 관계를 유지한다.
2. 계절적 변화는 어류 GM의 예측된 기능적 역량을 형성하는 지배적인 힘이며, 이는 먹이 가용성(해조류 대 동물성 먹이) 및 병원균 압력과 같은 환경적 변화에 의해 유도되는 것으로 보인다.
3. 이러한 구별된 메커니즘을 이해하는 것은 환경 변화에 대한 어류 GM의 반응을 예측하고, GM 변동의 주요 동인을 고려한 효과적인 프로바이오틱스를 설계하는 데 필수적이다.

저자들은 숙주 형질을 직접 측정하지 않았고 기능적 경로는 실험적으로 검증된 것이 아니라 추론된 것이라는 한계점을 언급하며 신중한 태도를 유지하였다. 따라서 이러한 메커니즘을 확증하기 위해서는 숙주의 생리학적 측정과 직접적인 기능 분석을 포함하는 추가 연구가 필요하다고 명시하였다.