식물 종 간 대사산물 분포의 전역 구조와 형성 메커니즘

초록

**
플라보노이드와 같은 2차 대사산물이 식물 종 사이에 어떻게 분포하는지를 이중 네트워크로 분석하였다. 실험 데이터는 6개 대형 과에서 4 725종·6 846개의 플라보노이드를 포함한다. 네트워크는 이질적인 연결도(멱법칙), 중첩성(nestedness), 모듈성(modularity)을 동시에 보이며, 기존의 특성 기반 BC 모델보다 단순한 진화 모델(복제‑발산)로 더 정확히 재현된다. 이는 대사산물 분포가 복잡한 특성 상호작용이 아니라 기본적인 진화 과정에 의해 형성될 수 있음을 시사한다.

**

상세 분석

**
본 연구는 식물 종‑플라보노이드 관계를 이중(bipartite) 그래프로 모델링하고, 네트워크 구조의 세 가지 핵심 특성인 이질적 연결도, 중첩성, 모듈성을 정량적으로 평가하였다. 데이터는 Metabolomics.JP에서 추출한 14 378개의 종‑플라보노이드 쌍을 기반으로 하며, 특히 Fabaceae, Asteraceae, Lamiaceae, Rutaceae, Moraceae, Rosaceae 등 여섯 개 대형 과에 초점을 맞추었다. 연결도 분포는 플라보노이드와 종 모두에 대해 멱법칙 형태에 지수적 절단을 보였으며, 이는 기존 연구에서 보고된 ‘power‑law‑like’ 특성과 일치한다. 중첩성(N)은 0.5~0.9 사이의 높은 값을 나타냈고, 무작위 널 모델(null model 2) 대비 통계적으로 유의미(p < 0.0001)하였다. 모듈성(Q) 역시 높은 값을 기록했으며, N과 Q 사이에 상관관계가 거의 없다는 점에서 두 구조가 독립적인 현상임을 확인했다.

두 가지 모델을 비교하였다. 첫 번째는 저자들이 제안한 진화 모델로, (i) 새로운 종이 기존 종을 복제하고, 각 플라보노이드를 상속 확률 q로 유지하거나 소실하고, (ii) 새로운 플라보노이드는 기존 종‑플라보노이드 쌍을 무작위 선택해 추가한다는 두 단계 프로세스로 구성된다. 모델 파라미터 p와 q는 실제 데이터의 종 수(S), 플라보노이드 수(F), 상호작용 수(L)로부터 식(1),(2)를 이용해 추정하였다. 두 번째는 기존의 BC(Bipartite Cooperation) 모델로, 종과 플라보노이드의 특성(trait)과 외부 요인(environmental factor)을 사전에 지정하고, 고정된 종 수와 상호작용 수를 기반으로 네트워크를 생성한다.

구조 예측 정확도는 Pearson 상관계수(CC)와 RMSE로 평가했으며, 진화 모델은 N에 대해 CC = 0.949, RMSE = 0.0039, Q에 대해 CC = 0.770, RMSE = 0.0282로 BC 모델보다 우수했다. 또한, 종 및 플라보노이드의 차수 분포에 대한 Kolmogorov‑Smirnov(wKS) 거리에서도 진화 모델이 더 작은 값을 기록, 특히 플라보노이드 차수(P(k_F))에 대한 적합도가 현저히 높았다. 이러한 결과는 복제‑발산 메커니즘만으로도 복잡한 네트워크 구조를 재현할 수 있음을 강력히 뒷받침한다.

결론적으로, 플라보노이드와 같은 2차 대사산물의 종 간 분포는 특성 기반 상호작용보다 기본적인 진화 과정—즉, 종 복제와 대사산물의 점진적 변이—에 의해 형성되는 것이 보다 타당하다. 이는 생물학적 네트워크 연구에서 ‘특성‑기반’ 모델이 반드시 필요하지 않으며, 간단한 진화 규칙만으로도 실제 관찰되는 복합 구조를 설명할 수 있음을 시사한다.

**