IF50.5｜幽门螺杆菌存在古老生态种Hardy

幽门螺杆菌分化出两个古老生态种，Hardy与Ubiquitous共存于人类胃中。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

幽门螺杆菌长期定植人类胃部，与溃疡和胃癌相关，其毒力因子存在广泛地理变异。此前研究已根据多位点序列分型将全球菌株划分为多个种群，如hpAfrica2、hpSiberia等，并揭示了与人类迁移的共进化关系。然而，在特定地区如西伯利亚和美洲原住民人群中，部分菌株在系统发育树上形成独立分支，与同地区其他菌株明显不同，提示可能存在更深层的分化。

科学问题

卡在哪里： 这些独立分支菌株与同地区其他菌株的基因组差异程度、起源时间及维持机制尚不清楚。此前研究未系统比较其全基因组变异模式，也未阐明这种分化是否代表一个独立的生态种，以及其与人类饮食、宿主适应之间的关系。

本文要回答： 本文旨在鉴定并命名这一新生态种，解析其与Ubiquitous菌株的基因组分化特征、进化历史及生态适应机制。

技术路线

作者主要用了：全基因组测序、GWAS、FST、系统发育分析、泛基因组分析、dN/dS。

作者收集了全球9188个幽门螺杆菌基因组，经质控后保留6864个，利用染色体涂染法将菌株划分为13个种群。通过比较不同方法下的聚类结果，发现一组来自西伯利亚和美洲原住民的菌株在系统发育树和PCA中独立成群，命名为Hardy。随后，作者对Hardy与同地区Ubiquitous菌株进行全基因组关联分析（GWAS）和FST扫描，鉴定出100个高度分化基因，并分别基于分化基因和未分化基因构建系统发育树，以解析其进化历史。

图 1　Fig. 1。图1a的GWAS曼哈顿图显示，在1.50-1.55 Mb区域存在多个显著关联的SNP，且这些SNP的FST值很高，表明该区域在Hardy和Ubiquitous之间高度分化。图1b和1c分别基于未分化基因和分化基因构建的系统发育树，显示Hardy菌株在未分化基因树上与同地区Ubiquitous菌株聚类，但在分化基因树上形成独立深分支。

核心亮点

亮点 1｜Hardy与Ubiquitous在100个基因上几乎固定分化

GWAS和FST分析显示，Hardy菌株与同地区Ubiquitous菌株在100个基因上存在近乎固定的SNP差异，这些基因多编码外膜蛋白和宿主互作因子。在未分化基因的系统发育树上，Hardy菌株与同地区Ubiquitous菌株聚类，但在分化基因树上，Hardy菌株形成比hpAfrica2更深的独立分支，表明其分化早于hpAfrica2的分化。

图 2　Fig. 2。图2a展示了假设的生态种分化与全球传播情景：幽门螺杆菌在人类进化早期分化为Hardy和Ubiquitous两个生态种，随后随人类迁移传播至各大洲。图2b显示基于分化基因的种群水平树，Hardy和Ubiquitous的种群间关系相似，支持平行传播。

亮点 2｜Hardy与Ubiquitous的克隆祖先分裂古老

点图分析显示，Hardy与Ubiquitous菌株之间的基因组重排数量远高于同一生态种内部（约140个 vs 约10个），表明两者在细胞水平上的分裂非常古老。尽管存在频繁的同源重组，分化基因仍保持各自的身份，提示选择作用维持了这种分化。

图 3　Fig. 3。图3a显示分化基因的单倍型分布，Hardy和Ubiquitous菌株中均存在少量对方生态种的单倍型，表明近期基因流但被选择清除。图3b展示Hardy等位基因数量与遗传距离的关系，大部分菌株符合趋势线，但有11个来自Hardy分布区的菌株明显偏离，提示局部基因流。

亮点 3｜Hardy菌株携带铁辅因子脲酶和额外vacA拷贝

大多数Hardy菌株拥有第二个脲酶（UreA2B2），其辅因子为铁而非镍，这一特征通常见于食肉动物胃内的螺杆菌。此外，Hardy菌株还携带两个额外的vacA基因拷贝。这些基因内容差异可能反映了对不同宿主饮食的适应。

图 4　Fig. 4。图4a的泛基因组聚类显示Hardy菌株形成独立分支。图4b显示Hardy菌株中cagA、vacA和ureA/B的分布：大多数Hardy菌株缺乏完整cag致病岛，但拥有额外的vacA拷贝和铁辅因子脲酶基因ureA2B2。图4c展示了ureA/B和vacA的基因组排列。

亮点 4｜Hardy与Ubiquitous在分化基因上呈现平行传播历史

基于分化基因构建的种群水平系统发育树显示，Hardy和Ubiquitous菌株在非洲、萨胡尔、西伯利亚和美洲之间具有相似的拓扑结构和分支长度，表明两者可能由相同的人类迁移事件传播。Hardy菌株在美洲原住民中的分布可能对应更早的迁移。

亮点 5｜分化基因受到正选择且功能富集于外膜蛋白

dN/dS分析显示，分化基因在Hardy与Ubiquitous之间的比较中dN/dS升高，提示正选择驱动了分化。功能富集分析发现，100个分化基因中有21个编码外膜蛋白，显著富集，此外还包括铁镍获取和酸适应相关基因，支持生态位特化假说。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： 全基因组测序、GWAS、FST、系统发育分析、泛基因组分析、dN/dS
• 公开数据： 原文未给出公开组学登录号
• 代码： 公开仓库 https://github.com/EliseTourrette/HpEcospecies_Tourrette2023、https://github.com/univieCUBE/gepard；本地已克隆：gepard、HpEcospecies_Tourrette2023
• 谁该点开原文： 研究细菌进化、群体基因组学或幽门螺杆菌致病机制的研究者。

带走一句

细菌种群内可长期维持高度分化的生态种，即使存在广泛重组。分析时应区分基因组不同区域，避免用单一系统发育树掩盖生态分化信号。