# 中文标题 山毛榉树皮小蠹真菌群落组成的生态学特征：宿主、海拔与生境梯度下的真菌组成一致性

《Microbial Ecology》：The Ecology of Beech Bark Beetle Mycobiomes: Similar Fungal Composition Despite Host, Elevation and Site Gradients

　　

# 摘要翻译

树皮甲虫（bark beetles）可与真菌形成多种多样的关联，然而许多取食树皮的落叶树种甲虫的真菌群落（mycobiome）研究仍较为匮乏，其生态学意义亦有待深入探索。研究人员利用内转录间隔区2（internal transcribed spacer 2, ITS2）宏条形码技术（metabarcoding）对采自德国六个区域、61个森林斑块的山毛榉树皮小蠹（Ernoporicus fagi与Taphrorychus bicolor）的真菌群落进行了表征。尽管研究沿海拔梯度、实验性森林管理处理以及六个区域对两种甲虫进行了采样，但甲虫物种身份、处理措施与海拔对甲虫真菌群落组成变异的解释程度均较小，区域尺度仅呈现微弱效应，且在更细致的检验中消失。然而，研究发现真菌群落组成在更为精细的单个森林样地（forest site）尺度上变异更为一致，且斑块之间的地理距离具有弱但显著的影响，表明这些树皮甲虫所携带的真菌群落是由局地的、斑块水平的因素所塑造，而非更广泛的区域或宿主相关因素。两种甲虫的真菌群落总体高度相似，并以广泛分布的腐生真菌（saprotrophs）和酵母菌（yeasts）为优势类群。

# 论文解读

## 研究背景与意义

随着气候变化导致树皮甲虫（bark beetles）爆发频率增加、对森林的影响愈发严重，树皮甲虫与真菌的共生关系日益受到关注。以大小蠹属（Dendroctonus）和齿小蠹属（Ips）为代表的针叶树甲虫-真菌系统研究较为深入，但许多取食树皮的落叶树种甲虫的真菌群落（mycobiome）仍研究不足，其沿环境梯度的变异规律亦不明确。山毛榉树皮小蠹Ernoporicus fagi与Taphrorychus bicolor均取食欧洲山毛榉（Fagus sylvatica）的韧皮部，二者虽常在中欧山毛榉林中同域分布，但在交配系统、坑道结构与物候期上存在差异：E. fagi为单配制、雌虫钻蛀不规则纵向坑道，五月与七月飞行；T. bicolor为多雌制、多雌虫钻蛀星形坑道，三至六月飞行。既往对这两种甲虫真菌关联的认识完全依赖培养法，且仅涉及个别属，其真菌群落整体仍未得到表征。因此，阐明其真菌群落的组成与变异驱动因素，对于理解甲虫生物学、预测环境变化下昆虫-真菌共生关系的响应具有重要的基准意义。本研究据此提出四项假设：甲虫物种间真菌组成存在差异（H1）；真菌组成沿海拔梯度发生温度驱动的转变（H2）；异质林分中的真菌多样性更高（H3）；区域尺度显著影响真菌组成（H4）。

## 主要技术方法

研究人员依托德国BETA-FOR（BETA-FOR实验框架）的61个森林斑块（六区域：维尔茨堡大学森林、巴伐利亚森林、帕绍森林、洪斯吕克森林、萨尔兰森林与吕贝克森林）开展多区域野外研究，在每个斑块悬挂一根新鲜山毛榉原木，经自然定殖后以羽化陷阱（emergence traps）收集羽化的甲虫，基于形态特征鉴定并在必要时以COI条形码（COI barcoding）Sanger测序验证。对整虫进行DNA提取（ZymoBIOMICS? 96 MagBead DNA Kit），采用真菌特异性引物5.8 S-Fun_1与ITS4-mod扩增ITS2区，在MiSeq平台进行250 bp双端测序。生物信息学处理使用DADA2（v1.26.0）构建扩增子序列变体（amplicon sequence variants, ASVs），并采用仅正向读取（forward only）策略以保留长度相关分类群，基于UNITE数据库进行物种分类。统计分析采用Hill数（q=1即Hill-Shannon指数）量化alpha多样性，以置换多元方差分析（permutational multivariate analysis of variance, PERMANOVA）检验beta多样性组间差异，采用基于距离矩阵的多元回归（multiple regression on distance matrices, MRM）评估空间自相关，并运用多水平模式分析（multilevel pattern analysis）鉴定指示物种。

## 研究结果

**数据质量与稳健性**：研究共收集157头甲虫（E. fagi 66头、T. bicolor 91头），因单头虫体生物量不足而按原木合并为83个样本；多个甲虫合并未影响真菌群落组成（R2=0.012, p=0.792）。iNEXT稀疏曲线表明优势成分（q=1、q=2）已被充分采样。阴性对照与空白检出低读长，模拟群落可检出预期菌属，经decontam（阈值0.5）鉴定并剔除6个污染ASV（占读长0.11%）后，主要PERMANOVA结果未受影响；羽化陷阱效应不显著（Wilcoxon检验p=0.933），说明羽化陷阱未系统塑造真菌群落。

**H1不成立——甲虫物种间真菌组成相似**：尽管两种甲虫在属水平上比例略有差异（E. fagi中枝孢属（Cladosporium）与青霉属（Penicillium）平均相对丰度较高，T. bicolor中耶氏酵母属（Yamadazyma）与曲霉属（Aspergillus）较高），但群落组成无显著差异（PERMANOVA：完整模型p=0.1185；仅物种p=0.0521），alpha多样性亦无差异（Kruskal–Wallis, p=0.085）。指示种分析在错误发现率（false discovery rate, FDR）校正后未发现与任一宿主显著关联的ASV，Hyphopichia、Phallus、Trechispora等属仅在T. bicolor样本中出现，但整体上两种甲虫真菌群落均被腐生真菌与酵母菌主导。

**H2不成立——海拔无显著影响**：任一甲虫的alpha多样性（Shannon多样性与ASV丰富度）与海拔均无显著相关关系（p≥0.378），拒绝温度驱动的海拔过滤假说。

**H3不成立——森林经营处理无显著影响**：处理斑块与对照斑块间alpha多样性（p=0.902）与群落组成（PERMANOVA, R2=0.0110, p=0.8435）均无差异，表明死木与林冠操纵所营造的结构异质性并未改变甲虫真菌群落。

**H4不成立——区域效应微弱且不稳健**：区域尺度仅呈现微弱效应（R2=0.105, p=0.026），但在两项敏感性检验中消失：排除样本量仅1–2个的区域后不再显著（R2=0.071, p=0.094）；两个样本量最大的区域（维尔茨堡大学森林与巴伐利亚森林）直接比较亦无差异（R2=0.027, p=0.435）。

**空间结构存在于更精细尺度**：以森林样地为空间单位时，真菌群落组成呈现显著的空间结构化（PERMANOVA, R2=0.208, p=0.005），且在两个采样最完整的区域内独立于区域-样地混杂效应依然显著（R2=0.174, p=0.028）；MRM分析显示地理距离是真菌群落不相似性的显著预测因子（coefficient=0.0100, p=0.011），解释2.2%的方差（R2=0.022, p=0.033），表明真菌群落的空间结构存在于比区域更精细的尺度。

**优势真菌类群**：枝孢属、青霉属、耶氏酵母属与曲霉属为两物种共有优势属，其中枝孢属与青霉属在几乎所有样本中检出，曲霉属出现于约75%样本；Hypholoma fasciculare高流行；Talaromyces rugulosus多见于T. bicolor；Aureobasidium pullulans在两物种中均常见但在T. bicolor中略高。在已知甲虫共生菌中，Geosmithia属在23个样本中检出、两物种水平相当，但Geosmithia cnesini几乎仅见于萨尔兰区域的3个E. fagi样本，Geosmithia microcorthyli仅为微量检出；Ophiostoma属仅在单个T. bicolor样本中检出一次；昆虫病原菌Beauveria仅见于两个T. bicolor样本。

## 讨论与结论

研究指出，两种甲虫重叠的真菌群落与森林管理、海拔的微弱效应一致于如下观点：取食落叶树的韧皮部取食甲虫大多不具有专性的、物种特异的核心共生菌，而是从基质中获取真菌并沿水平与种间交换大量分类群，兼性的、生态灵活的关联策略可能是此类系统的普遍特征。需要承认的是，所有原木运至同一设施统一羽化一年，后期发育在部分标准化条件下进行，可能均质化了区域间差异，但定殖与早期幼虫发育仍发生在源地林分，故所测群落仍反映局地的真实定殖组装；甲虫未表面消毒、羽化陷阱向环境空气开放，无法判定单个分类群的来源（原木、陷阱环境或甲虫自身），但因所有样本经历相同的野外与陷阱条件，该局限主要影响单类群生物学解释而不影响组成比较；ITS2扩增对Ophiostoma等类群存在已知的检测偏差，其近乎缺失可能部分源于方法学局限而非自然缺失。未来研究应增加单样地样本量、整合局域微气候数据、关注单坑道尺度及跨代垂直传播分析，以揭示潜在的特异性专性共生菌。

**研究结论**：研究表明，山毛榉树皮小蠹E. fagi与T. bicolor与多样的、几乎一致的森林腐生真菌及酵母菌真菌群落相关联，唯一例外是几乎仅出现于E. fagi的Geosmithia cnesini。区域、处理与海拔等广泛的景观尺度梯度仅微弱影响这些关联，而真菌群落组成在单个森林样地及地理距离尺度上呈现更一致的结构，提示区域尺度的采样设计可能遗漏甲虫真菌群落中真实的空间结构。对于取食落叶树的树皮甲虫而言，以灵活、腐生真菌为主的关联而非专性共生菌，可能是这类山毛榉树皮小蠹的典型特征。