关于《Porphyromonas gingivalis - Porphyromonas gulae Duo：探讨牙周炎起源的系统基因组分析》的学术研究报告
本研究由来自智利多所大学的研究人员共同完成。第一作者（并列）为Mauricio Morales-Olavarría和Josefa Núñez-Belmar，通讯作者为Cristian Cortez和Juan P. Cárdenas。参与机构包括：Universidad Mayor（智利圣地亚哥）的基因组学与生物信息学中心以及生物技术学院，以及Pontificia Universidad Católica de Valparaíso（智利瓦尔帕莱索）的医学技术学院。该研究成果以原创研究论文的形式发表于2023年7月19日的《Frontiers in Microbiology》期刊上。
一、 学术背景 本研究属于微生物基因组学、进化生物学及口腔医学的交叉领域，聚焦于牙周病的关键病原体——牙龈卟啉单胞菌（Porphyromonas gingivalis）。牙周炎是一种由菌斑生物膜引发的慢性免疫炎症性疾病，导致牙齿支持组织的破坏。P. gingivalis在其中扮演“基石病原体”（keystone pathogen）的角色，即其少量存在即可破坏微生物群落平衡并驱动疾病进程。
该研究的动因在于，尽管P. gingivalis的致病性已被广泛研究，但其在所属卟啉单胞菌属（Porphyromonas）中的进化历程，以及其如何获得作为人类牙周炎关键病原体的特性，仍不甚清晰。值得注意的是，与P. gingivalis亲缘关系极近的P. gulae存在于多种哺乳动物的口腔中，并同样与牙周病相关。这一“P. gingivalis-P. gulae”对为研究病原体宿主适应性进化提供了理想模型。
因此，本研究旨在利用公开的基因组数据，通过比较基因组学方法，达成四个主要目标：(a) 通过系统基因组分析确认P. gingivalis在卟啉单胞菌属中的系统发育位置；(b) 定义并比较P. gingivalis和P. gulae的泛基因组；© 评估P. gingivalis和P. gulae在分化过程中的基因家族获得/丢失事件；(d) 通过计算Tajima’s D值和dN/dS比率，比较这两个物种基因家族的进化压力（即自然选择强度）。
二、 详细工作流程 本研究完全基于生物信息学分析，其工作流程严谨而系统，主要包括以下几个步骤：
基因组数据集构建与物种界定：首先，研究者于2022年8月从NCBI GenBank下载了所有标注为卟啉单胞菌属的基因组。他们运用了两个关键过滤标准以确保数据质量：一是使用GTDB-Tk工具进行精确的物种分类学鉴定；二是使用CheckM工具评估基因组的完整度（≥90%）和污染度（%），最终从233个初始基因组中筛选出高质量数据集。为了界定“基因组物种”集群，他们进一步计算了所有基因组间的平均核苷酸一致性（Average Nucleotide Identity, ANI），并以95%为界进行聚类分析，最终确定了36个基因组物种集群。其中，P. gingivalis有84个高质量基因组，P. gulae有14个，这构成了后续分析的核心数据集。
基因注释、直系同源群鉴定与系统发育树构建：使用Prokka软件对所有选定基因组进行从头注释。随后，使用OrthoFinder工具在全数据集（233个基因组）中鉴定直系同源群（Orthogroups, OGs），即推测起源于单一祖先基因的一组基因。为了构建可靠的系统发育树，研究者从所有基因组中提取了38个单拷贝且保守的直系同源群，将其蛋白质序列进行多序列比对（使用MAFFT），然后基于比对结果利用IQ-TREE软件构建最大似然系统发育树，并以Tannerella forsythia作为外群。该树揭示了整个卟啉单胞菌属的进化关系。
泛基因组分析：为了深入比较P. gingivalis和P. gulae，研究者对这两个物种的基因组集合（84个和14个）分别独立运行了OrthoFinder，以构建各自的泛基因组。他们计算了泛基因组累积曲线，并用幂律方程拟合以判断其开放性（γ>0表示开放）。根据基因家族在菌株中的存在频率，将泛基因组划分为核心基因组（存在于100%菌株）、软核心基因组（90-99.9%）、壳基因组（15-89%）和云基因组（<15%）。同时，使用eggNOG-mapper对所有预测蛋白质进行功能分类（COG类别），并比较了不同泛基因组组分中功能类别的分布差异。
基因获得/丢失模型构建：为了追溯进化过程中基因内容的变化，研究者将P. gingivalis、P. gulae的全部基因组以及一个P. loveana基因组（作为外群）一起进行直系同源群分析。利用获得的直系同源群存在/缺失矩阵和系统发育树，使用Count软件基于Wagner简约性原则，推断了在P. gingivalis和P. gulae分化节点及其最近共同祖先（Last Common Ancestor, LCA）中发生的基因获得与丢失事件的数量。
进化压力分析：此项分析包含两个部分。一是计算dN/dS比率（ω）：研究者鉴定了P. gingivalis和P. gulae之间共有的单拷贝直系同源群，并以P. loveana为外群，使用PAML软件中的codeml程序分别计算每个基因家族在这两个物种配对比较中的非同义替换率（dN）与同义替换率（dS）的比值。ω值通常用于推断自然选择类型（ω<1：纯化选择；ω>1：正选择）。二是计算Tajima’s D值：该统计量基于群体内DNA序列的多态性，用于检测偏离中性进化模型的情况。研究者分别计算了P. gingivalis核心基因家族（1,335个）和P. gulae核心基因家族（1,367个）的Tajima’s D值分布。负值表示罕见等位基因过多，可能暗示种群扩张或正选择；正值则表示等位基因频率稳定或种群收缩。
毒力因子进化分析：为了探究关键毒力因子的进化起源，研究者使用隐马尔可夫模型（HMM）搜索，针对特定的Pfam数据库条目（代表不同毒力因子结构域），在整个卟啉单胞菌属数据集中进行扫描。这些因子包括：主要菌毛FimA（Pfam PF06321）、次要菌毛Mfa1（PF15495）、血凝素/牙龈蛋白酶共有结构域（PF07675, PF10365）、牙龈蛋白酶特有结构域（PF01364）以及肽基精氨酸脱亚胺酶（Peptidyl-arginine deiminase, PAD; PF04371）。获取相关蛋白质序列后，进行多序列比对并构建系统发育树，以分析这些毒力因子在属内的分布与进化关系。
三、 主要研究结果 1. 卟啉单胞菌属系统发育与物种关系：系统发育树证实了P. gingivalis与P. gulae是高度相关的姐妹谱系。一个关键且新颖的发现是，这两个物种与P. loveana（一种从麝鼠袋鼠口腔中分离的物种）构成了一个紧密的进化支。研究还发现了许多尚未被正式命名的基因组物种，揭示了该属存在大量未被充分认识的遗传多样性。这84个P. gingivalis和14个P. gulae基因组被用于后续深入分析。
P. gingivalis与P. gulae的开放泛基因组：对两个物种的泛基因组分析显示，它们的累积曲线符合开放泛基因组特征（γ值分别为0.194和0.204）。功能分析表明，核心基因组显著富集了与管家功能（如能量产生、翻译、细胞壁合成等）相关的COG类别。相反，附属基因组（壳、云、独特基因）则显著富集了功能未知的基因以及属于“防御机制”和“可移动元件”类别的基因。这表明，虽然维持基本生命活动的功能核心相对稳定，但两个物种都在通过获得或丢失与适应、防御和未知功能相关的基因来动态调整其基因组内容。
物种分化中的基因获得与丢失：基因获得/丢失模型推断，P. gingivalis与P. gulae分化节点的共同祖先拥有约1,662个基因家族。在P. gingivalis谱系的最远共同祖先中，预测获得了94个基因，丢失了95个基因，净变化很小。P. gulae谱系的最远共同祖先则获得了35个基因，丢失了65个基因。值得注意的是，两个物种祖先获得的基因中，很大一部分是功能未知的“假定蛋白”。在功能已知的获得基因中，包括与IX型分泌系统（T9SS，对毒力因子分泌至关重要）、硝基还原酶（可能与口腔硝化应激适应有关）、转录调控因子以及可移动遗传元件（如转座酶、松弛酶）相关的基因。这些发现提示，在物种分化过程中，除了功能未知基因的流动，与特定环境适应（如硝化应激）和遗传物质转移相关的基因可能扮演了重要角色。
进化压力与选择信号：dN/dS分析显示，P. gingivalis和P. gulae共享的核心基因绝大多数受到强烈的纯化选择（ω中位数远小于1）。通过比较两个物种各自与P. loveana的ω值，研究者识别出少数可能在一个谱系中受到更强正选择的基因家族。例如，P. gingivalis中的4-羟基丁酰辅酶A脱水酶（参与厌氧代谢）和UvrABC系统蛋白B（DNA修复），以及P. gulae中的D-丙氨酸-D-丙氨酸连接酶（细胞壁合成）等。这暗示了两个物种在分化后，各自在不同的功能通路上经历了特化的适应性进化。Tajima‘s D值分析揭示了更显著的差异：P. gingivalis核心基因的Tajima’s D值分布显著偏向负值（比P. gulae更低）。这种大量罕见等位基因的频率升高，通常解释为经历了种群瓶颈效应后的扩张事件。研究者据此提出一个假设：P. gingivalis在适应并定植于人类口腔这一特定生态位的过程中，可能经历了一个遗传瓶颈，随后种群扩张，导致其遗传多样性呈现出现今观察到的模式。而宿主范围更广的P. gulae则表现出相对更稳定的等位基因频率。
毒力因子的保守性与进化起源：对关键毒力因子的系统发育分析取得了极具启发性的结果：
四、 结论与意义 本研究的核心结论是：导致人类牙周炎的“基石病原体”P. gingivalis，与其动物口腔中的近亲P. gulae以及新近识别出的近亲P. loveana，构成了一个进化上紧密相关的“三重奏”。这个进化支的共同祖先很可能已经具备了包括菌毛、血凝素、牙龈蛋白酶（至少是Rgp）和PAD酶在内的多种关键毒力因子“工具箱”。因此，P. gingivalis作为人类牙周炎关键病原体的“能力基础”并非其独有或近期全新进化而来，而是从其与动物病原体共享的祖先那里继承的。
这一发现具有重要的科学价值：它改变了我们对牙周炎病原体起源的认知，将视角从单一的“人类病原体”扩展到了“人畜共患病原体进化支”的背景之下。研究提示，牙周炎相关病原体的进化历史可能远比想象的古老，可能追溯至哺乳动物早期分化的时代。此外，研究发现P. gingivalis在适应人类宿主过程中可能经历了遗传瓶颈，而P. gulae和P. loveana也拥有类似的毒力装备，这为未来研究不同宿主（人类、其他哺乳动物、有袋类）中卟啉单胞菌引发的牙周病差异及其分子机制提供了全新的进化框架和比较模型。
五、 研究亮点 1. 新颖的系统发育关系发现：首次通过大规模基因组数据分析，明确将P. loveana定位为与P. gingivalis和P. gulae紧密相关的第三个物种，形成了“P. gingivalis-P. gulae-P. loveana”进化支，为理解牙周炎相关病原体的进化提供了关键节点。 2. 毒力因子进化史的重新阐述：通过属范围的系统发育分析，有力证明了FimA、Mfa1、血凝素、Rgp牙龈蛋白酶、PAD等核心毒力因子在该进化支分化前就已存在，挑战了这些因子是P. gingivalis为适应人类宿主而特化获得的传统观点。 3. 多维度进化压力分析：综合运用泛基因组、基因获得/丢失模型、dN/dS和Tajima’s D值等多种进化生物学分析工具，从基因内容动态、自然选择类型和群体历史多个角度，全面刻画了P. gingivalis和P. gulae的进化轨迹，并提出了P. gingivalis经历宿主适应瓶颈的合理假设。 4. 高质量数据与严谨流程：研究基于大量高质量基因组（84个P. gingivalis, 14个P. gulae, 总计233个属内基因组），采用严格的生物信息学流程，确保了分析结果的可靠性，并揭示了卟啉单胞菌属中存在大量未描述物种的现状。
六、 其他有价值的发现 研究还指出，P. gulae的核心基因组中含有更多与CRISPR-Cas系统（一种原核生物防御机制）相关的“防御机制”基因，而P. gingivalis中此类基因更多存在于附属基因组中。这可能反映了二者在各自微生物群落中面临的噬菌体或质粒压力不同。此外，许多在物种分化中获得的基因功能未知，这些“暗物质”基因可能在宿主适应和致病过程中扮演尚未被认识的角色，是未来功能研究的重要靶点。