鹿科动物肠道微生物组中定植因子谱系的基因组解析揭示生态分层

《Frontiers in Microbiology》：Genome-resolved analysis of colonization factor repertoires reveals ecological stratification in cervid gut microbiomes

　　

摘要

引言：定植因子（Colonization Factors, CFs）是与肠道微生物持久性和宿主适应相关的重要微生物性状，然而它们在鹿科动物肠道微生物组中的大规模组织方式仍不清楚。方法：研究人员分析了来自15个公共可用项目和1个内部数据集的688份鹿科动物肠道宏基因组样本中衍生的共3,311个非冗余高质量宏基因组组装基因组（Metagenome-Assembled Genomes, MAGs）。通过与GHA CF数据库比对鉴定CF相关基因，并在基因组、宿主物种和胃肠道区段水平上表征CF谱系。结果：共鉴定出跨越71个CF家族的138,729个CF相关基因。来自鹿亚科（Cervinae）的MAGs比羊亚科（Caprinae）含有更丰富的CF谱系，其中CF47（Peptidase_C69）、CF24_29（QueH）和CF18（Glycos_transf_2）是最常见的家族之一。CF谱系受分类学强烈结构化，与细菌系统发育距离呈中度关联，并形成两种具有不同KEGG功能谱的重复性基因组水平构型。整合样本元数据进一步揭示了CF谱系在宿主物种和胃肠道区段间的分化，这代表了本研究考察的主要生态维度。与区段相关的CF变异伴随更广泛功能谱的重新分布，包括空肠中碳水化合物和脂质代谢的富集、回肠中膜运输的富集、盲肠中外源生物降解的富集以及瘤胃中环境适应的富集。讨论：这些发现提供了鹿科动物肠道微生物组中CF谱系组织的基因组解析视图，并证明定植相关功能在微生物谱系和生态背景中结构化分布。本研究强调了在解释哺乳动物肠道生态系统中CF谱系分布时，考虑微生物分类学和宿主相关环境的重要性。

论文解读：鹿科动物肠道微生物组中定植因子谱系的生态分层——基于基因组的比较分析

**研究背景与问题**

微生物在肠道中的稳定定植是宿主-微生物长期互作的基础。定植因子（Colonization Factors, CFs）包括营养获取、环境适应、转运、黏附及与宿主表面相互作用等相关基因，是影响微生物持久性和宿主适应的重要功能模块。然而，以往对CF的研究多集中于单一菌株或有限实验系统，对复杂肠道微生物组中CF的大规模组织方式缺乏系统认识。鹿科动物（Cervidae）因其宿主多样性与反刍动物高度区室化的消化系统相叠加，为研究CF在宿主物种和胃肠道区段间的分布提供了理想模型。已有研究表明，肠道微生物群落在宿主物种间、近缘宿主谱系间及更广泛反刍动物类群间存在显著差异，同时胃肠道不同区段在理化条件、底物可利用性和宿主分泌物等方面各异，可能对微生物持久性特征施加不同的选择压力。因此，CF谱系不太可能作为均一的功能背景存在，而可能随宿主物种和肠道区段系统变化。然而，仅记录生态差异不足以解释CF组织的机制，需要明确CF谱系在微生物基因组中的分布是否受分类学和系统发育约束、是否形成重复的基因组水平构型，以及这些构型如何随宿主相关生态背景变化。为此，研究人员提出假设：CF谱系在鹿科动物肠道微生物组中呈层级结构化，而非随机分布，具体表现为微生物分类学和系统发育约束基因组水平的CF组织，而宿主物种和胃肠道区段驱动CF组成的生态分化。

**研究内容与结论**

为验证上述假设，研究人员整合了来自15个公共宏基因组项目和1个内部数据集的688份鹿科动物肠道样本，经质量过滤和去冗余后获得3,311个非冗余高质量宏基因组组装基因组（MAGs，完整性≥90%，污染≤5%）。通过比对GHA CF数据库，共鉴定出138,729个CF相关基因，归属71个CF家族。研究结果显示：CF谱系在鹿科动物肠道微生物组中呈现显著的生态分层模式。在基因组水平，CF组织与微生物分类学强相关，与细菌系统发育距离呈中度正相关，并形成两种重复性基因组构型，这两种构型具有不同的分类组成和KEGG功能谱。在生态水平，CF组成和多样性在宿主物种间及胃肠道区段间显著分化，且这种分化与更广泛的功能谱差异相伴随。该研究首次提供了鹿科动物肠道微生物组中CF谱系组织的基因组解析视图，证明了定植相关功能在微生物谱系和生态背景中的结构化分布，强调了在解释哺乳动物肠道生态系统中CF谱系分布时考虑微生物分类学和宿主相关环境的重要性。研究论文发表在《Frontiers in Microbiology》。

**主要关键技术方法**

研究人员采用了以下关键技术方法：（1）宏基因组组装与基因组分箱：使用fastp进行质量过滤，Bowtie2去除宿主序列，MEGAHIT组装，MetaBAT2分箱生成MAGs；（2）基因组质量评估与去冗余：采用CheckM2评估质量，dRep以99%平均核苷酸一致性去冗余；（3）分类学注释与系统发育重建：使用GTDB-Tk进行分类学分配，PhyloPhlAn3基于蛋白质序列重建系统发育树；（4）CF基因鉴定：利用DIAMOND BLASTP将预测蛋白与GHA CF参考数据库比对，以E-value≤1×10^–10且覆盖率≥70%为阈值；（5）功能注释：使用Pfam-A HMM数据库和KEGG数据库对CF相关蛋白进行结构域和直系同源注释；（6）统计与可视化：采用Jaccard和Bray–Curtis距离、PCoA、PAM聚类、PERMANOVA、Kruskal–Wallis检验及Spearman相关分析等。样本队列来源包括15个公共项目和1个内部数据集，共688份鹿科动物肠道样本。

**研究结果**

*基因组解析的鹿科相关微生物组系统发育结构*：从15,494个MAGs中经质量过滤保留4,995个高质量MAGs，去冗余后获得3,311个非冗余MAGs。系统发育树显示这些MAGs主要归属于Bacillota_A和Bacteroidota，两者合计占比超过70%，另有Pseudomonadota、Verrucomicrobiota等门类。属水平上Alistipes、HGM04593、Cryptobacteroides等最为丰富。

*CF谱系的基因组解析编目与分类学结构化*：共鉴定138,729个CF相关基因，归属71个CF家族，其中CF47（Peptidase_C69）、CF24_29（QueH）和CF18（Glycos_transf_2）最为普遍。稀释曲线表明CF发现尚未饱和。CF流行率在属和种水平上均表现出显著的分类学结构化，Bacillota_A和Bacillota的属种通常具有更高且变异更大的CF流行率，而Campylobacterota和Pseudomonadota则较低。

*鹿亚科与羊亚科MAGs的CF谱系跨亚科比较*：Jaccard排序显示两数据集间清晰分离（PERMANOVA，R²=0.114，P=0.001），鹿亚科MAGs比羊亚科含有显著更多的CF家族（Wilcoxon检验，P<2.2×10^–16）。两亚科共享高流行率核心CF家族如CF43、CF35、CF15等，但额外家族流行模式存在差异。

*CF谱系的分类学结构化与聚类*：CF谱系受微生物分类学强烈结构化，与细菌系统发育距离呈中度正相关（Spearman ρ=0.341）。PAM聚类将基因组水平CF谱系分为两种主要构型。Cluster 1主要关联营养摄取、氨基酸和辅因子代谢、PTS介导的碳水化合物利用、转运系统和生物膜形成等通路；Cluster 2则富集碳和能量代谢、细胞包膜生物合成和应激相关代谢适应等通路。聚类成员与门水平（Cramér's V=0.801）和科水平（Cramér's V=0.868）组成强相关，Bacillota_A优势谱系主要归入Cluster 1，Bacteroidota相关谱系集中于Cluster 2。

*宿主物种相关的CF谱系生态分化*：四种鹿属宿主（Cervus nippon、Cervus albirostris、Cervus canadensis、Cervus elaphus）的CF基因组成显著不同（PERMANOVA，R²=0.206，P=0.001），CF基因丰富度和Shannon多样性均存在显著差异（Kruskal–Wallis，P<0.001）。聚类水平上，Cluster 1在所有宿主中占主导，但相对丰度在宿主间仍有差异，Cervus elaphus的Cluster 1比例最高。多种优势CF家族在特定宿主中优先富集，不同宿主关联部分不同的CF家族组合。

*胃肠道区段相关的CF谱系分化*：瘤胃、空肠、回肠、盲肠和结肠五个区段间CF基因组成存在显著差异（Bray–Curtis PCoA显示分离），CF基因丰富度和Shannon多样性均显著不同（Kruskal–Wallis，P<0.001）。空肠和瘤胃样本CF基因丰富度最高，盲肠最低，回肠和结肠居中。聚类水平上，Cluster 1在所有区段占绝对优势，但Cluster 2在瘤胃和空肠的相对丰度高于回肠、盲肠和结肠。

*区段相关KEGG功能分层*：KEGG功能谱在胃肠道区段间重新分布，整体组成存在显著差异（PERMANOVA，R²=0.223，P=0.001）。空肠群落富集碳水化合物代谢、聚糖生物合成与代谢及脂质代谢；回肠以膜运输为特征；盲肠显示更强的外源生物降解与代谢；瘤胃则具有较高的环境适应和翻译功能。结肠通常处于中间位置。多个优势CF家族与这些区段变化的KEGG类别呈显著正或负相关，表明CF组成与更广泛功能谱在区段间共变。

**讨论与结论**

讨论部分指出，CF谱系的强分类学结构化表明其可能遵循基因组组织的普遍规则，即生态相关功能受进化历史和宿主相关选择共同塑造。鹿亚科-羊亚科比较为鹿科CF组织提供了更广泛的反刍动物背景，但差异也可能反映宿主生态和微生物群落分类结构的变化，不应仅归因于进化分歧。两种重复CF构型表明CF内容变异是结构化的而非随机分布，但尚不能确证实验验证的定植策略。宿主物种差异需在微生物谱系依赖性背景下解读，宿主相关生态过滤可能重新加权分类学结构化的微生物谱系库。Cervus albirostris因其局限于青藏高原而提供了可解释的生态案例，其CF水平分化与能量代谢富集共现，提示持久性功能与更广泛代谢适应可能响应相同环境压力。胃肠道区段差异中，瘤胃和空肠的高CF丰富度及Cluster 2的较高贡献与各自的消化功能和代谢条件相关，盲肠的窄CF谱系与其外源生物降解功能背景相一致，但均为相关性而非因果证据。

研究结论为：鹿科动物肠道微生物组中CF谱系呈现层级的生态分层模式。在基因组水平，CF组织与微生物分类学强相关，与细菌系统发育距离中度相关，并通过两种重复CF构型表达。在生态水平，CF组成和多样性在宿主物种和胃肠道区段间分化，伴随区段相关功能变异。这些数据库定义的谱系模式是否转化为微生物持久性、宿主适应或微生物组稳定性的差异，需进一步的实验和纵向验证。