松茸（Tricholoma matsutake, TM）是一种在欧亚大陆备受推崇的野生食用菌，其经济价值极高，全球零售市场估计达5亿美元。它通过与宿主树木根系形成外生菌根（ectomycorrhizal, EcM）共生而生存。然而，由于气候变化、寄主病害、森林砍伐和过度采伐等因素，其野生种群面临严重的可持续性挑战。因此，研究TM与其共生微生物的相互作用，对于理解其生态学并开发保育与栽培策略至关重要。
TM的菌根圈，即“shiro”，是一个独特的生态位，其中定殖着由TM结构化的特殊微生物群落。自20世纪60年代开创性研究以来，大量研究利用元条形码（Next Generation Sequencing）和条形码（Sanger测序）技术，从土壤（S-soil）、TM定殖的根系（S-root）和子实体（S-fruitbody）等不同生态位鉴定了大量shiro相关微生物。然而，由于测序平台、靶标区域、参考数据库和分类学分配工具等技术异质性，以及研究结果在不同分类水平上报告的不一致，导致我们对这些微生物的认知是零散且割裂的。为了克服这一碎片化问题，并进行全球性的比较，本研究进行了首次跨研究的整合分析。
研究方法
研究团队系统收集了已发表及新生成的细菌（16S rRNA基因）和真菌（ITS区域）扩增子序列数据，最终整合了来自芬兰、中国（西南和东北地区）、韩国和日本等地理区域的28项研究，共计44个数据集。为了整合技术异质性的序列，研究采用了创新的基于系统发育树的聚类方法。具体流程包括：序列预处理和质量过滤，构建高质量的细菌和真菌参考数据库及系统发育骨干树，将查询序列通过进化位置放置整合到骨干树中，随后利用uDance等工具对树进行分解和细化。最后，通过TreeCluster等算法将序列聚类为谱系单元，并利用自定义的“树行走”算法进行一致的分类学分配，从而定义了微生物“谱系型”。这种基于系统发育的聚类方法，有效解决了因技术差异导致的分类学不一致问题，使得跨研究比较成为可能。
此外，研究还进行了微生物多样性分析、生态位特异性鉴定、群落功能预测以及关键细菌谱系的功能验证实验。通过比较shiro土壤与非shiro土壤（NS-soil），以及shiro内部的土壤、根部和子实体等不同生态位，来识别核心和特异性的微生物谱系。细菌群落的功能通过PICRUSt2工具基于MetaCyc通路和酶委员会编号进行预测。为了在谱系水平上验证预测的功能，研究挑选了代表S-soil、S-root和S-fruitbody特异性谱系的8株细菌菌株，利用表型微阵列技术测试了它们对166种不同底物的利用能力。
研究结果
1. 跨研究序列数据集与微生物谱系
研究构建了一个涵盖广泛地理区域和生态位的全球shiro微生物序列数据集。经过处理，共获得59,770条细菌和15,241条真菌查询序列，最终聚类得到39,530个细菌谱系型和8,350个真菌谱系型。分析显示，shiro土壤样本的微生物α多样性显著低于非shiro土壤，并且群落组成存在明显差异。
2. 测序方法比较与地理分布模式
比较元条形码和培养条形码两种方法发现，两者在检测某些类群上存在重叠，但也存在明显的偏好性。例如，两种方法都能频繁检测到S-soil中的Caballeronia和Paraburkholderia，以及子实体中的Ewingella和Pseudomonas。然而，培养方法倾向于高估能形成抗性结构（如芽孢）的革兰氏阳性菌（如Bacillus, Paenibacillus, Streptomyces）以及真菌中的Penicillium和Trichoderma。相反，一些类群如Alphaproteobacteria, Acidobacteriia等在培养方法中被低估。这些方法学偏差突显了有针对性地分离未培养核心微生物的重要性。
尽管shiro微生物群落的系统发育结构在不同地理区域内部和之间存在变异，但细菌群落的系统发育差异与地理距离并无显著相关性。例如，相距超过3500公里的中国东北和西南地区的细菌群落相似性较高，而韩国境内短距离采集的样本却表现出高度异质性。
3. shiro核心与生态位特异性微生物谱系
尽管存在地理和技术差异，研究成功鉴定出在shiro中全局保守的“核心”微生物谱系型，包括39个细菌谱系型和12个真菌谱系型。它们虽然只占谱系型总数的很小比例（细菌0.4%，真菌1.1%），却贡献了shiro群落总丰度的很大一部分（细菌31.0%，真菌72.5%（不含TM时为10.4%））。细菌核心谱系主要属于Proteobacteria（22个）和Acidobacteria（11个），其中Rhizobiales谱系型（与Bradyrhizobium AXAI_s密切相关）和Paraburkholderia谱系型最为丰富和普遍。真菌核心谱系中，Umbelopsis dimorpha和Oidiodendron chlamydosporicum的普遍性最高。
通过多重分析方法，研究进一步鉴定出7个细菌和4个真菌谱系型为S-soil特异性（相对于NS-soil）。例如，Paraburkholderia谱系型 108.3和真菌Mucor zonatus在S-soil中 consistently表现出更高的丰度。尽管U. dimorpha在某一研究中因TM的绝对优势而丰度较低，但其与TM consistently的正相关性仍使其被鉴定为S-soil特异性。
通过比较shiro内部的生态位，研究揭示了高度特异性的谱系。在S-root细菌群落中，Acinetobacter谱系型 100.30和Ralstonia insidiosa占主导地位。在S-fruitbody细菌群落中，Pseudomonas谱系型 97.28、Cupriavidus pauculus和Ewingella americana最为丰富。在S-root真菌群落中，除了TM，Phialocephala谱系型 33.50占主导（12.0%），而在S-soil中非常罕见（<1%），显示出明确的S-root特异性。Thelephora ganbajun和Epicoccum谱系型 24.5被鉴定为S-fruitbody特异性真菌。
对S-fruitbody微生物组的跨界共现网络分析揭示了一个模块化结构，形成了两个特征鲜明的群组：一个与TM正相关（TM-positive），包含E. americana、Pseudomonas谱系型 97.28等；另一个与TM负相关（TM-negative），包含S-root特异性的Ra. insidiosa和Acinetobacter谱系型 100.30等作为主要枢纽。这两个群组彼此间存在强烈的负相关。
4. 跨生态位的功能分化
功能预测分析揭示了微生物群落在不同生态位中的功能分化。通过比较S-soil与NS-soil、S-root及S-fruitbody，鉴定出49条富集的代谢通路和417个富集的酶促反应。
S-soil富集功能：与氧化应激耐受相关的通路（如硫酸盐还原、乙醛酸循环、血红素生物合成）在S-soil中特异性富集。与膜合成相关的通路在S-soil和S-fruitbody中均富集。
S-fruitbody富集功能：S-fruitbody中特异性富集的通路包括2,3-丁二醇生物合成、L-组氨酸降解、烟酸酯降解I、TCA循环VII（乙酸生产者）、磷酸吡哆醛生物合成I、去甲亚精胺生物合成、葡萄糖氧化降解、鸟氨酸降解超通路以及多粘菌素抗性等。
跨生态位保守的芳香化合物降解通路：三条芳香化合物降解通路（儿茶酚降解为β-酮己二酸、儿茶酚降解III、通过β-酮己二酸的芳香化合物降解）在所有比较中（S-soil vs. NS-soil, S-root, S-fruitbody）均一致性地富集。这些通路属于“通过3-氧代己二酸的芳香化合物降解超通路”。其中，原儿茶酸和儿茶酚开环反应的关键酶（EC 1.13.11.3和EC 1.13.11.1）及其下游降解为3-氧代己二酸的反应在S-soil中 consistently富集，并且在S-root和S-fruitbody中进一步富集。此外，3-脱氢莽草酸降解反应（EC 4.2.1.118）在S-root中特异性富集，而水杨酸降解反应（EC 1.14.13.172）在S-root和S-fruitbody中均富集。
对于真菌群落，与NS-soil相比，S-soil中 consistently富集的功能包括线粒体生物发生、DNA复制蛋白、氧化磷酸化、心肌收缩和细胞粘附分子等通路。与S-soil相比，S-root中则富集了抗菌药耐药基因和细菌运动蛋白相关通路。
5. 细菌功能的实验验证
表型微阵列实验验证了不同生态位特异性细菌谱系的底物利用能力，结果与群落水平的功能预测相呼应。
S-soil特异性谱系：Paraburkholderia谱系型 108.3和Caballeronia菌株优先利用簇5中的有机酸，包括三种芳香化合物（4-羟基苯甲酸、间羟基苯乙酸、2-羟基苯甲酸）和三种非芳香化合物（L-酒石酸、癸二酸、乙醛酸）。其中，Paraburkholderia谱系型 108.3对4-羟基苯甲酸的利用能力最强。这与S-soil中芳香化合物降解和乙醛酸循环通路的富集相吻合。此外，这些菌株还偏好利用簇6中的底物，如戊糖醇和景天庚酮糖，这对应了磷酸戊糖途径的激活。
S-root和S-fruitbody特异性谱系：S-root特异性的Paenibacillus谱系型主要利用植物相关底物，如水杨苷、水苏糖、果胶、蔗糖和D-棉子糖，这些底物很少被其他测试细菌利用。有趣的是，虽然Paenibacillus能高效利用水杨苷，但不利用其衍生物水杨酸。S-fruitbody特异性的Pseudomonas谱系型 97.28则高效利用L-精氨酸、L-鸟氨酸和L-组氨酸等氨基酸，这与子实体中氨基酸降解通路的富集一致。该菌株还对衣康酸和D-海藻糖表现出最高的呼吸和生长活性。
讨论与意义
本研究首次在全球范围内，通过系统发育整合框架，系统比较了shiro微生物组。尽管存在技术和地理异质性，我们鉴定出了一组对维持群落稳定功能至关重要的shiro核心微生物谱系型，并揭示了S-soil、S-root和S-fruitbody中具有明确生态位分化的特异性谱系和功能。
研究揭示了方法学偏差，强调了未来需要针对性分离未培养的核心和生态位特异性谱系，以更好地合成人工shiro微生物组，这是实现TM人工栽培的关键一步。
核心谱系的时空一致性强化了这些微生物在TM生存和生长所需关键生态过程中发挥作用的观点。值得注意的是，历史上一些被误分类的微生物（如曾被归为Pseudomonas的Caballeronia/Paraburkholderia，曾被归为Mortierella的Umbelopsis dimorpha）在本研究中得到了更准确的鉴定。
生态位特异性谱系的功能暗示了它们与TM相互作用的多样性。例如，尽管Caballeronia和Paraburkholderia在含葡萄糖的培养基上显示出对TM菌丝的抑制作用，但它们在以有机酸为唯一碳源时可能起到菌根辅助细菌的作用。S-soil核心真菌U. dimorpha和Mu. zonatus的菌丝提取物能促进TM生长，暗示其作为菌根辅助真菌的潜力。此外，Acidobacteria门成员（如Acidipila KB906735_s和Edaphobacter aggregans）作为S-soil核心组分，可能通过分解顽固性植物材料和耐受shiro中的酸性及抗菌化合物环境，在养分循环中扮演关键角色。
S-root被Acinetobacter和Ralstonia insidiosa主导，它们与多种菌根相关，提示其在菌根发育中的潜在作用。S-root特异性真菌（如Phialocephala fortinii）很可能在降解衰老的菌根根系中发挥重要作用。在S-fruitbody中，与TM正相关的细菌（如E. americana, Pseudomonas谱系型 97.28）可能作为菌根辅助细菌促进菌丝生长并抑制有害霉菌。
总之，这项跨研究整合为理解松茸菌根圈微生物组的组成、功能和生态位分化提供了全面的视角。鉴定出的核心和特异性微生物谱系及其功能特征，为未来针对性的生态学研究、微生物分离培养以及开发基于微生物组的松茸保育和栽培策略奠定了坚实基础。
