栖息地干扰下生态特化灵长类的微生物与代谢灵活性：黑白领狐猴的整合多组学证据

《Molecular Ecology》：Microbial and Metabolic Flexibility in Response to Habitat Disturbance in an Ecologically Specialist Primate

　　

： 在人类世，理解何种特征使物种易于灭绝对野生动物管理至关重要。生态特化种（ecological specialists）因其形态、生理和/或行为上的限制被认为难以应对快速栖息地退化，因而被假设面临特别高的灭绝风险。然而，已有研究对此假设的支持结果不一，由此引出一个关键问题："特化种对环境变化的适应弹性（resilience）究竟有多大？"本研究检验了如下假设：受行为和生理限制的生态特化种在退化栖息地中可能处于能量不利地位。具体而言，研究人员在12个月期间检验了生活在次生林中的极度濒危食性特化灵长类——黑白领狐猴（Varecia variegata）——相较于生活在原始林栖息地的同种个体，是否遭受营养和能量亏缺。为此，研究人员采用混合模型方法，检验了行为、营养化学、16S测序、宏基因组功能预测（metagenome functional predictions）、代谢物谱（metabolite profiles）与能量结局之间的关系。与原始林中的同种个体相比，研究人员发现退化森林中的动物从多样性较低的饮食中摄入了略少的卡路里。这些动物表现出较低的肠道微生物组（gut microbiota）多样性、降低的微生物功能潜力以及改变的代谢组学谱（metabolomic profiles）。然而，尽管存在明显的营养限制，不同栖息地之间的能量结局大致相似。这些发现表明，生物体的肠道微生物组（gut microbiome）可能能够调节微生物代谢潜力，以促进在次优条件下的适应弹性。这些发现凸显了宿主-微生物组互作作为生态特化种适应弹性的重要组成部分，对预测持续环境变化中物种存续具有广泛意义。

论文解读

一、研究背景与科学问题

随着全球人为干扰日益加剧，动物正面临不断累积的生态胁迫，进化错配（evolutionary mismatch）的可能性随之上升。生态特化种（ecological specialists）因生态位宽度狭窄，可利用的资源、栖息地及环境条件范围有限，被认为对环境变化尤为脆弱，其特化特征长期与灭绝风险相关联。尽管特化种常因行为或食性灵活性有限而被假定为脆弱，但适应性可能通过其他生物学尺度的补偿机制浮现，例如生理或微生物可塑性可能通过提高消化效率、重新分配代谢途径或稳定宿主能量平衡来部分抵消能量代价。然而，这类补偿机制能否缓冲特化种应对退化栖息地的影响，仍知之甚少。

栖息地退化预期会改变可利用食物的丰度和营养组成，这可能通过限制特化种获取偏好资源而对其构成约束。营养摄入的改变进而可影响肠道微生物群落组成和功能潜力，因为肠道微生物对底物可利用性变化高度敏感，并在消化、发酵和宿主能量获取中发挥核心作用。在此框架下，宿主能量状况并非直接响应栖息地退化，而是源于资源可利用性、饮食行为、微生物功能与宿主代谢之间的相互作用。本研究以马达加斯加雨林特化种黑白领狐猴（Varecia variegata）为对象，比较了生活在原始林与再生次生林中的两个种群，检验肠道微生物组可塑性是否有助于缓冲生态特化种面临的环境与营养挑战。

二、研究设计与主要技术方法

研究人员在Ranomafana国家公园（RNP）对两个黑白领狐猴种群开展了为期12个月（2022年8月至2023年7月）的纵向研究：Mangevo（MGV）为原始雨林样地，含20只佩戴无线电信标的成年个体；Talatakely（TLK）为曾受伐木干扰的再生次生林样地，含5只成年个体。研究整合了多维度数据：通过物候监测计算每月食物可利用性指数（FAI_m）；通过全天焦点动物跟踪收集全发生摄食行为数据（共1727小时）；采集摄食植物进行营养化学分析，测定粗蛋白（CP）、可利用蛋白（AP）、粗脂肪、中性洗涤纤维（NDF）和总非结构性碳水化合物（TNC）等宏量营养素；采集粪便样本（N=185）进行16S rRNA基因V4-V5区域扩增子测序，经QIIME2流程处理，使用Deblur去噪和SILVA数据库分类，并通过PICRUSt2预测宏基因组功能；对40个粪便样本子集进行非靶向代谢组学（Q-TOF MS）和短链脂肪酸（SCFAs）定量分析（LC-MS/MS）；采集尿液样本（N=124）测定C-肽（uCP）浓度作为能量平衡指标。统计分析采用混合模型、PERMANOVA、平衡重采样、贝叶斯广义线性混合模型（GLMMs）等方法。

三、研究结果

3.1 跨栖息地的营养相似性

次生林中的狐猴摄入热量较少（MGV均值410 kcal，TLK均值337 kcal；t=1.898, df=212, p=0.059），但未达显著水平，且两样地的季节性能量摄入模式相似：12月至3月最高，5月至9月最低。尽管次生林食物可利用性和饮食多样性较低，狐猴在两个栖息地表现出高度相似的营养平衡模式，均摄入高碳水化合物、中低脂肪、低蛋白的饮食，仅在非蛋白能量（NPE）组分（即碳水化合物与脂肪）比例上呈现季节性变化。这表明这些生态特化种在营养上高度受限，在次生或退化栖息地中可能难以摄入足够热量。

3.2 跨栖息地的微生物与代谢差异

次生栖息地的狐猴种群表现出中等程度的α多样性降低（Shannon熵和Faith系统发育多样性）、弱分化的β多样性差异及中等程度的核心微生物组差异。平衡重采样分析支持Shannon多样性在MGV普遍高于TLK（中位差=0.24，95% CI: ?0.20–0.61），但仅20.3%的比较显著。微生物相对丰度在栖息地间差异显著（Poisson GLMM后验均值=?0.463，pMCMC=0.006），最大差异出现在Firmicutes和Bacteroidota两个最丰富门类中，若干在原始栖息地存在的稀有门类在次生栖息地完全缺失。宏基因组功能预测显示，最丰富的代谢途径（碳水化合物代谢23.1%、氨基酸代谢19.3%、能量代谢13.7%等）在两样地相当，但18条KEGG和62条MetaCyc代谢途径存在显著差异，且几乎所有差异途径均在原始栖息地富集。四种主要SCFAs（乙酸盐、丁酸盐、丙酸盐、异戊酸盐）产量在两样地大致相似，但MGV均略高。代谢组学分析显示显著差异：在61,608种代谢物中，1430种化合物代表的88条代谢途径在栖息地间表现出显著强度差异；次生林中脂肪酸氧化代谢物强度较低，而淀粉和蔗糖代谢物强度较高。

3.3 营养摄入、微生物代谢与能量结局的关联

狐猴肠道微生物组整体群落组成在营养摄入类别间无显著变化，但特定营养素的摄入与特定微生物ASV亚群的丰度变化相关：NDF、TNC和脂肪摄入与若干Bacteroidota和Cyanobacteria ASV的丰度变化相关，AP摄入则与任何ASV无关联。营养摄入虽不能预测微生物群落组成，但显著预测了推定功能途径丰度，AP和脂肪摄入是整体途径相对丰度的最佳预测因子（后验均值_AP=?27,985，pMCMC=0.049；后验均值_Fat=?30,409，pMCMC=0.028）。宏量营养素摄入减少分别与不同功能途径的富集相关：AP摄入减少与多种氨基酸代谢途径丰度升高相关，TNC摄入减少与多种碳水化合物代谢途径丰度升高相关，脂肪摄入减少与多种脂质代谢途径丰度升高相关。能量平衡方面，样地对能量平衡（C-肽ng/mL）有边际效应（后验均值=?4.871，pMCMC=0.054），TLK仅略低于MGV；MGV在9月至12月出现明显的C-肽升高，而TLK缺乏此季节性升高。宏基因组代谢途径丰度显著正向预测能量平衡变异（后验均值=2.148，pMCMC=0.004），解释22.2%的总方差；代谢物强度也显著正向预测能量平衡变异（后验均值=0.006，pMCMC<0.001），解释73%的总方差（CI 7%–96%）。

四、讨论与结论

与预期相反，生活在再生次生栖息地中的黑白领狐猴——极度濒危的生态特化种——相对于原始林中的同种个体并未处于显著能量劣势。相反，数据表明通过微生物功能的快速转变实现的代谢灵活性可能在代谢缓冲中发挥关键作用，有助于在次优环境条件下维持能量状况。尽管次生栖息地中热量摄入略低且微生物组成改变，但两样地的能量结局大致相似。微生物功能途径丰度和代谢物强度均正向预测尿C-肽变异，表明微生物代谢响应可能与宿主能量稳定性相关。值得注意的是，尽管次生栖息地中宏基因组功能预测和代谢组学谱显示功能潜力降低——这常被视为肠道微生物组健康下降的指标——但这些功能差异未伴随广泛的分类更替。此外，虽然次生栖息地中整体途径丰富度和代谢物多样性降低，但两样地在较低营养摄入下均富集了与营养利用相关的特定代谢途径，提示存在靶向功能响应而非代谢能力的均匀丧失。这些模式表明环境变化可能更容易改变微生物群落的功能而非其成员组成，凸显功能可塑性作为适应弹性的潜在机制。

然而，这些发现不应被解释为跨栖息地完全能量等效的证据。次生种群的尿C-肽边际较低（pMCMC=0.054），且缺乏原始种群中观察到的季节性C-肽显著升高，这可能表明能量缓冲不完全或资源丰盛期积累能量储备的能力降低，可能对繁殖或长期适合度产生重要下游后果。两个种群在食物可利用性降低时均转向节能行为，提示宿主摄食决策与肠道微生物响应可能协同作用，使特化种在资源景观改变时维持能量稳定。研究设计上，单一原始林样地与单一再生次生林样地的比较使栖息地类型与样地身份部分混淆，无法完全排除未测量的样地特异性因素的影响，因此研究结果应被视为两个具有已知干扰历史、资源可利用性和食性生态差异的狐猴种群的局部比较。

综上，研究结果展示了肠道微生物组在缓解宿主面临的能量和营养挑战中的关键作用。这种代谢灵活性可能是包括生态特化种在内的部分动物应对环境条件改变的重要组成。整合行为、营养、微生物、代谢组学和生理数据的多因素机制框架，为评估生态特化种在栖息地退化下的适应弹性提供了有力工具，对预测持续环境变化中的物种存续具有广泛意义。从应用角度看，维持或增强饮食多样性和关键功能性食物资源的栖息地恢复策略，可能有助于支持生态特化种的肠道微生物和代谢灵活性。