Seasonal Dynamics of Microbial Communities, Co-occurrence Networks, and Grazing-Mediated Carbon Flow Characteristics

기술 요약: 미생물 군집의 계절적 역학, 공생 네트워크 및 Grazing 매개 탄소 흐름 특성

문제 제기
연안만(Coastal bays)은 강한 수문학적 및 영양염 변화를 특징으로 하며, 이는 미생물 생물 다양성을 재편하고 미생물 먹이망을 통한 탄소 전달을 변화시킬 수 있다. 기존 연구들은 온도, 혼합, 영양염 가용성에 의해 구동되는 뚜렷한 계절적 순환을 미생물 군집이 보인다는 점을 확립해 왔으나, 여전히 중요한 지식 격차가 존재한다. 특히, 양식업의 영향을 받는 연안만에서 세 가지 핵심 차원—(1) 미생물 군집 구성 및 조립 메커니즘(결정론적 vs 확률적), (2) 통합된 박테리아-미세유기체(microeukaryotic) 공생 네트워크의 구조적 조직, (3) 직접적인 속도 기반의 grazing 매개 탄소 전달 측정—을 하나의 공유된 시간적 틀 안에서 통합하여 다룬 연구는 드물다. 또한, 계절적 환경 강제력이 식물 플랑크톤과 박테리아 플랑크톤 경로 모두로부터의 탄소 흐름 효율을 어떻게 동시에 조절하는지에 대한 실증적 증거도 제한적이다. 본 연구는 이러한 격차를 해소하기 위해 양식업의 영향을 크게 받는 동해만(Dinghai Bay, 중국 남동부)을 대상으로 조사하였다.

연구 방법론
본 연구는 4번의 계절별 크루즈(2020년 여름, 2020년 가을, 2021년 겨울, 2021년 봄)에 걸쳐 환경 DNA(eDNA) 메타바코딩, 군집 조립 모델링, 네트워크 분석 및 생리적 속도 측정을 결합한 다각적인 접근 방식을 채택하였다.

• 샘플링 및 시퀀싱: 간조대 양식 구역, 전이 지대(ecotones), 외해 지역의 24개 정점에서 표층 해수를 채집하였다. 0.2 µm 필터로부터 환경 DNA를 추출하였다. 박테리아 군집은 16S rRNA 유전자 앰플리콘 시퀀싱(V4 영역)을 통해, 미세유기체 군집은 18S rRNA 유전자 시퀀싱을 통해 특성화하였다.
• 군집 분석: 표준 지수와 Bray-Curtis 비유사도를 사용하여 알파 및 베타 다양성을 평가하였다. 군집 조립은 중립 군집 모델(Neutral Community Model, NCM)을 사용하여 확률적(중립적) 과정과 결정론적 과정을 구분하여 평가하였다. 부분 환형 다변량 분석(pRDA)과 맨텔 테스트(Mantel tests)를 사용하여 군집 구성을 수문학적 및 영양염 변수(예: 온도, pH, 용존 무기 질소)와 연결하였다.
• 네트워크 분석: 각 계절별로 박테리아 ASV, 미세유기체 ASV 및 환경 변수를 포함하는 통합 공생 네트워크를 구축하였다. 네트워크는 위상적 특성(모듈성, 클러스터링 계수, 차수 분포)을 분석하였으며, 이를 Erdős–Rényi 무작위 그래프와 비교하였다.
• Grazing 실험: 해수 희석 실험을 수행하여 식물 플랑크톤과 박테리아 플랑크톤 모두에 대한 고유 성장률($\mu$)과 미세 동물 플랑크톤의 grazing rate($g$)를 추정하였다. 탄소 흐름 효율은 일일 생산자의 생산량 중 미세 동물 플랑크톤에 의해 소비된 비율로 계산되었다.
• 통계적 종합: 환경 구배, 생물 다양성 및 탄소 흐름 과정 사이의 광범위한 연관성을 탐색하기 위해 부분 최소 제곱 구조 방정식 모델링(PLS-SEM)을 사용하였다.

주요 결과

• 계절적 다양성 및 구성: 박테리아와 미세유기체 군집 모두 뚜렷한 계절적 전환(turnover)을 보였다. 알파 다양성(풍부도 및 Shannon 지수)은 겨울에 정점에 달하고 여름에 최저치를 기록했다. 박테리아 군집은 Pseudomonadota가 우점하였고, 미세유기체는 Metazoa, Ciliophora, Dinoflagellata가 우점하였다.
• 환경적 동인: 군집 비유사성은 온도, pH, 염도 및 용존 무기 질소(DIN)와 강하게 연관되어 있었다. pRDA는 이러한 변수들이 정점 효과를 고려한 후에도 유의미한 주요 동인임을 확인하였다.
• 군집 조립: 중립 군집 모델 분석 결과, 박테리아 군집은 미세유기체 군집(95% 신뢰 구간 내 ASV 약 27%)에 비해 중립적 기대치와 더 높은 일치성(약 54%의 ASV가 신뢰 구간 내 존재)을 보였다. 미세유기체는 중립적 예측으로부터 더 강한 이탈을 보였으며, 이는 결정론적 선택 또는 기타 비중립적 과정의 영향이 더 크다는 것을 시사한다.
• 공생 네트워크: 통합 네트워크는 높은 연결성, 강한 모듈성을 보였으며, 양의 상관관계(80–90%의 edge)가 지배적이었다. 미세유기체-미세유기체 간의 연결이 가장 큰 비율의 edge를 차지했다. 네트워크 연결성(edge 밀도 및 평균 차수)은 여름에 정점에 달했다. 모든 네트워크는 두꺼운 꼬리를 가진 차수 분포를 가진 small-world 위상을 나타냈으며, 이는 핵심-주변부(core-periphery) 구조를 시사한다.
• 탄소 흐름 역학:
  • 생산: 식물 플랑크톤 1차 생산과 박테리아 플랑크톤 생산는 여름에 가장 높았다.
  • 효율: 생산율과는 대조적으로, grazing 매개 탄소 흐름 효율은 식물 플랑크톤과 박테리아 플랑크톤 모두에서 여름에는 일반적으로 낮았고 겨울/가을에는 높았다.
  • 경로: 모든 계절에 걸쳐 식물 플랑크톤 생산은 박테리아 플랑크톤 생산보다 grazing 매개 탄소 전달에 실질적으로 더 많이 기여하였다.
• 구조와 기능의 디커플링(Decoupling): PLS-SEM 분석 결과, 계절이 생물 다양성에는 강력한 영향을 미쳤으나 성장 및 grazing 과정에 미치는 직접적인 영향은 약했다. 생물 다양성에서 탄소 전달율로 가는 직접적인 경로는 거의 0에 가까웠으며, 이는 속도 기반의 탄소 전달이 광범위한 군집 패턴보다 밀접하게 결합되어 있지 않음을 시사한다.

주요 기여 및 의의
본 연구는 양식업의 영향을 받는 만(bay)에서 계절적 환경 변화가 미생물 시스템을 어떻게 재편하는지에 대한 프로세스 중심의 평가를 제공한다. 주요 기여는 다음과 같다:

1. 통합적 프레임워크: eDNA 기반의 군집 구조, 중립적 조립 신호, 계절적 공생 위상을 독립적으로 측정된 grazing 기반 탄소 흐름과 성공적으로 연결하였으며, 이는 연안 연구에서 드물게 결합된 형태이다.
2. 교차 도메인 조립 통찰: 박테리아와 미세유기체 간의 뚜렷한 조립 메커니즘을 강조하며, 이 환경에서 미세유기체 군집이 박테리아 군집보다 비중립적 과정에 의해 더 강하게 형성됨을 입증하였다.
3. 탄소 흐름의 디커플링: 연구 결과, 미생물 군집 구조와 구성이 계절에 따라 평행하게 변함에도 불구하고, 탄소 전달 효율(기능적 지표)은 이러한 구조적 변화를 단순히 반영하지 않는다는 것을 밝혔다. 특히, 높은 여름 생산량이 높은 전달 효율로 이어지지 않는다는 점은 탄소 흐름이 단순히 생산자 생물량이나 광범위한 군집 다양성이 아닌, 복잡한 속도 제한 메커니즘(예: grazer 조성, 먹이 품질)에 의해 제어됨을 시사한다.
4. 네트워크 위상: 미생물 네트워크의 계절적 재편을 규명하였으며, 모듈형 골격은 안정적으로 유지되지만 모듈 내 결합 강도는 환경 조건에 따라 변동하며 여름에 정점에 달한다는 점을 주목하였다.

결론
저자들은 계절적 환경 강제력이 동해만의 미생물 생물 다양성, 조립 신호, 공생 조직 및 grazing 매개 탄소 전달를 재구조화한다고 결론짓는다. 그러나 저자들은 속도 기반의 탄소 전달이 군집 구성 자체보다는 광범위한 군집 패턴과 덜 밀접하게 결합되어 있다는 점을 완곡하게 언급하였다. 본 연구는 역동적인 연안 환경에서 생태계 기능을 완전히 이해하기 위해 eDNA 기반의 네트워크 구조와 독립적인 생태학적 측정을 병행해야 할 필요성을 강조한다.