[论文解读] Linking statistical and ecological theory: Hubbell's unified neutral theory of biodiversity as a hierarchical Dirichlet process

Keith Harris, Todd L. Parsons+4 more|arXiv (Cornell University)|Oct 15, 2014

Gut microbiota and health参考文献 53被引用 4

一句话总结

本文提出了一种贝叶斯拟合方法来处理Hubbell的统一中性生物多样性理论（UNTB），通过证明多地点UNTB在大种群极限下收敛于分层狄利克雷过程（HDP），从而实现在多个地点上高效且精确的参数估计。该方法在低迁徙情景下优于现有近似方法，并揭示人类肠道微生物组在较高分类等级上表现出非中性动态，但在特定属内接近中性，且身体质量指数（BMI）与Ruminococcaceae中迁徙率呈负相关。

ABSTRACT

Neutral models which assume ecological equivalence between species provide null models for community assembly. In Hubbell's Unified Neutral Theory of Biodiversity (UNTB), many local communities are connected to a single metacommunity through differing immigration rates. Our ability to fit the full multi-site UNTB has hitherto been limited by the lack of a computationally tractable and accurate algorithm. We show that a large class of neutral models with this mainland-island structure but differing local community dynamics converge in the large population limit to the hierarchical Dirichlet process. Using this approximation we developed an efficient Bayesian fitting strategy for the multi-site UNTB. We can also use this approach to distinguish between neutral local community assembly given a non-neutral metacommunity distribution and the full UNTB where the metacommunity too assembles neutrally. We applied this fitting strategy to both tropical trees and a data set comprising 570,851 sequences from 278 human gut microbiomes. The tropical tree data set was consistent with the UNTB but for the human gut neutrality was rejected at the whole community level. However, when we applied the algorithm to gut microbial species within the same taxon at different levels of taxonomic resolution, we found that species abundances within some genera were almost consistent with local community assembly. This was not true at higher taxonomic ranks. This suggests that the gut microbiota is more strongly niche constrained than macroscopic organisms, with different groups adopting different functional roles, but within those groups diversity may at least partially be maintained by neutrality.We also observed a negative correlation between body mass index and immigration rates within the family Ruminococcaceae.

研究动机与目标

• 开发一种计算高效且精确的贝叶斯拟合策略，用于多地点统一中性生物多样性理论（UNTB），以解决一般情况下似然计算不可行的问题。
• 通过证明在大种群极限下，生态中性理论与机器学习中的统计模型——分层狄利克雷过程（HDP）——存在正式数学关联，建立两者之间的联系。
• 通过检验HDP近似对实证数据的拟合程度，区分非中性群落分布与完全中性的UNTB动态。
• 将该方法应用于真实世界数据集，包括热带树木群落和人类肠道微生物组，以评估中性在宏观与微观生态体系中的作用。

提出的方法

• 作者证明，具有大陆-岛屿结构和固定局部群落大小的一类广泛多地点中性模型，无论局部动态如何，只要物种在生态上等价，在大种群极限下均会收敛至分层狄利克雷过程（HDP）。
• 他们将现有的HDP贝叶斯推断技术（特别是基于‘中式餐厅特许经营’比喻的吉布斯采样）进行适配，以估计UNTB参数，包括基本生物多样性数（θ）和各站点特异的迁徙率（m_i）。
• 该方法采用非信息先验，并使用吉布斯采样算法，通过迭代重采样物种丰度和迁徙率，生成参数的完整后验分布。
• 通过在不同迁徙率和θ值下生成的模拟数据验证该方法，与Etienne的近似和精确似然方法进行性能比较。
• 将HDP近似应用于真实数据：巴拿马的热带树木数据集和包含570,851条序列、278个样本的人类肠道微生物组数据集，并在多个分类水平上测试分类分辨率。
• 该方法可通过比较HDP模型在不同假设下的拟合优度，实现中性检验：完全中性（元群落和局部群落均为中性）与非中性元群落但局部群落中性组装。

实验结果

研究问题

• RQ1 在大种群极限下，多地点UNTB是否无论局部群落动态如何，都会收敛至分层狄利克雷过程（HDP）？
• RQ2 HDP近似是否能为完整多地点UNTB提供比现有基于似然的近似方法更准确且计算更高效的贝叶斯拟合策略？
• RQ3 人类肠道微生物组在群落水平上是否与完全中性一致，还是功能生态位过程占主导地位？
• RQ4 在肠道微生物组中是否存在某些分类群，其局部群落组装近似中性，即使整个群落并非如此？
• RQ5 身体质量指数（BMI）是否与特定肠道微生物类群的迁徙率相关，提示肥胖状态下生态动态发生变化？

主要发现

• HDP近似在参数估计上比Etienne的近似方法更准确，尤其在低迁徙情景下（例如当I << θ时），HDP优于近似方法，并与Etienne的‘精确’方法性能相当。
• 对于巴拿马的热带树木数据集，HDP近似得到的参数估计（θ ≈ 231 ± 22，BCI站点迁徙率 ≈ 65.5 ± 5.9）与Etienne的‘精确’方法高度一致，验证了该方法在高θ条件下的有效性。
• 在人类肠道微生物组数据集中，完全群落中性被拒绝，表明在较高分类阶元上功能生态位过程占主导。
• 在某些属内——特别是Ruminococcaceae家族中——物种丰度分布几乎与局部中性组装一致，表明中性可能在维持功能群内多样性方面发挥作用。
• 在Ruminococcaceae家族中，发现身体质量指数（BMI）与迁徙率之间存在显著负相关，暗示肥胖可能降低外部迁徙的影响，转而有利于该关键碳水化合物代谢类群的本地生长。
• 基于HDP的方法具有可扩展性和高效性，能够对大规模微生物组数据集（如570,851条序列）实现完整后验推断，在复杂、低迁徙情景下的计算可行性与准确性均优于现有方法。