近克隆同一宿主粪便–血液沙门氏菌配对中经验证的小型Col样附属元件

《BMC Microbiology》：Validated small Col-like accessory elements in a near-clonal same-host stool–blood Salmonella pair

　　

**背景**：非伤寒沙门氏菌可在易感宿主中引起血流感染，但当配对粪便与血液分离株已经近乎克隆时，两者之间残留的细菌差异尚不明确。核心基因组SNP分析有助于评估亲缘关系，但可能忽略标准组装流程难以解析的小型环状附属元件。本研究调查了一对临床关联的同一宿主粪便–血液沙门氏菌配对是否含有经验证的小型Col样附属元件，而这些元件仅凭核心SNP分析无法显现。

**方法**：研究人员分析了一对时间相邻的同一宿主粪便–血液沙门氏菌指示配对FB1–XY1，以及两对短间隔粪便–粪便对照配对。该分析整合了严格reads支持的核心SNP过滤、附属元件分析、Illumina和ONT读长水平支持、PCR/Sanger/qPCR验证以及公共基因组背景筛查。

**结果**：严格reads支持过滤后，FB1–XY1配对仅保留4个核心SNP，而两对粪便–粪便对照配对均无保留的严格核心SNP。尽管核心基因组分歧极小，血液分离株XY1携带两个经验证的小型Col样附属元件。第一个是经组装解析的13.8-kb ColRNAI样质粒，含有ColRNAI样复制子命中序列和重复的cka/大肠杆菌素K样注释，并具有蛋白酶K敏感的抑制表型。第二个是5,241-bp的Col156样环状元件contig_40，该元件未在最终XY1闭合组装输出中表现为近全长复制子，但得到全长XY1 Illumina和ONT读长覆盖、116条ONT跨越连接点读长、PCR、Sanger测序、多菌落筛查和qPCR的支持。FB1 Illumina读长未显示对contig_40的相应支持。在筛查的16,456个公共沙门氏菌基因组背景中，两个元件的严格全长命中均罕见：13.8-kb质粒有2个命中，contig_40有15个命中。在402个筛查的Muenster基因组中，两个元件均无严格全长命中。

**结论**：这些发现将小型Col样环状附属元件确定为近克隆同一宿主粪便–血液沙门氏菌配对中一个经验证的、对组装敏感的基因组变异层。通过证明临床关联的粪便–血液分离株尽管核心基因组分歧极小，仍可在小型环状附属元件水平上存在差异，本研究将配对分离株解读扩展至核心SNP亲缘关系之外。整合读长水平证据、分子验证和公共基因组背景可能有助于改进临床微生物基因组学中小型环状元件的解读。

## 论文解读：近克隆同一宿主粪便–血液沙门氏菌配对中的小型Col样附属元件验证

### 一、研究背景与问题

非伤寒沙门氏菌（nontyphoidal Salmonella, NTS）是全球细菌性胃肠炎的主要病原体之一，并可在高危宿主中引起包括血流感染在内的侵袭性疾病。血流感染通常反映肠道或携带状态下细菌向肠道外播散，因此临床上不仅关注粪便与血液分离株是否相关，更关注二者在高度相近时究竟存在哪些细菌学差异。全基因组测序已成为沙门氏菌监测、暴发调查和临床分离株比较的核心工具，核心基因组SNP（core-genome SNP）和核心基因组多位点序列分型（cgMLST）广泛用于评估分离株亲缘关系。然而，当粪便与血液分离株在核心基因组水平已近乎克隆时，亲缘关系本身无法解释肠道与血液分离株之间残留的细菌差异。已有研究表明，同一患者的沙门氏菌分离株可能仅相差少数SNP，却仍存在基因含量、质粒含量或调控差异。小型质粒和环状附属元件作为组装敏感的基因组特征，在标准组装流程中可能被低估或遗漏，因此需要专门检验。

### 二、研究目的与设计

本研究旨在调查一对近克隆的同一宿主粪便–血液沙门氏菌配对是否携带经验证的小型Col样环状附属元件，从而超越单纯的核心SNP亲缘关系评估，探索组装敏感的精细附属基因组变异层。研究对象包括一对时间相邻的粪便–血液指示配对FB1（粪便）与XY1（血液），来自同一名急性淋巴细胞白血病儿科患者，两株菌均属于沙门氏菌Muenster血清型，间隔1天采集；另有两对短间隔粪便–粪便临床对照配对，分别为间隔5天的FB301–FB3（Infantis血清型）和间隔7天的FB4–FB401（Agona血清型），作为宿主内稳定性参照。研究整合了严格reads支持的SNP过滤、附属元件分析、Illumina和ONT读长支持、PCR/Sanger/qPCR分子验证以及公共基因组背景筛查。

### 三、关键研究方法

研究人员采用Illumina短读长测序完成全部6株菌的测序，并对XY1额外进行Oxford Nanopore长读长测序以生成闭合组装。核心SNP分析先经MUMmer/NUCmer进行组装级候选筛选，再通过自定义脚本进行reads支持过滤，排除短contig、contig边缘、高拷贝contig、重复区、低深度位点、非SNP及等位基因频率支持不足的位点，仅保留严格固定样核心SNP。小型附属元件通过高拷贝contig筛查、闭合组装比较、PlasmidFinder复制子注释和读长水平支持进行鉴定。contig_40的环状结构通过ONT读长映射至双拷贝参考序列并统计跨越首尾连接点的读长验证。分子验证包括内部PCR、连接点PCR、多菌落PCR、质粒富集DNA PCR、Sanger测序和基于SYBR的qPCR拷贝数估计。抑菌活性采用点种法测定，并用蛋白酶K处理评估蛋白质性质。公共背景筛查使用包含16,456个沙门氏菌基因组的SalmClinotype主目录，另设元数据丰富子集（11,036个基因组）、严格粪便–侵袭性主要子集（10,028个基因组）和Muenster背景（402个基因组），定义严格全长命中为≥90%查询覆盖率和≥95%加权核苷酸同一性。此外还选取106个公共记录进行定向raw-read敏感性分析。

### 四、研究结果

**三个真实世界同一宿主临床配对定义粪便–血液指示比较和粪便–粪便对照**
三对分离株均来自去标识的临床记录。指示配对FB1–XY1为同一天住院期间分别从粪便和血液培养获得的Muenster血清型菌株；两对粪便–粪便对照分别为Infantis和Agona血清型，间隔5天和7天采集，构成短间隔的宿主内稳定性参照。

**严格reads支持的SNP分析显示有限的核心基因组变异**
FB1–XY1初始组装级候选变异有7个，经严格reads支持过滤后仅保留4个核心SNP，包括potH G15G同义变异、ssb_1邻近基因间变异、putA R1171T和sgrR Q437E错义变异；两对粪便–粪便对照均无保留的严格核心SNP。这表明指示配对核心基因组接近克隆，但仍有少数可验证的位点差异。

**在XY1中检测到一个13.8-kb ColRNAI样质粒**
XY1闭合组装中解析出一个13,795-bp的环状复制子，为ColRNAI样质粒，包含17个ORF特征，携带ColRNAI_1复制子命中序列和重复的cka/大肠杆菌素K样注释，未检测到获得性耐药基因。spot-on-lawn实验显示XY1诱导上清对指示菌产生可测量的抑制圈，蛋白酶K处理后抑制圈消失，FB1上清和LB对照无抑制圈，表明该质粒赋予蛋白酶K敏感的蛋白质类抑制表型。

**一个5.2-kb Col156样contig_40环状元件得到读长水平和分子证据支持**
XY1短读长组装中存在5,241-bp的高拷贝contig_40，包含Col156_1复制子命中序列和5个预测ORF（均为假设蛋白），但未在最终闭合组装中恢复为近全长复制子。FB1 Illumina读长无任何覆盖支持；XY1 Illumina读长覆盖全元件（平均深度4,682.87×），XY1 ONT读长覆盖全元件（平均深度123.96×），且双拷贝参考映射发现116条跨过原contig_40首尾边界（5,241/5,242）的连接点读长。内部PCR和环形连接点PCR在XY1中均产生预期产物，而FB1无产物；多菌落PCR证实XY1多个独立菌落阳性、FB1阴性；Sanger测序确认内部产物映射到contig_40位置76–592（99.61%同一性，517 bp），连接点产物映射到双拷贝参考位置4,878–5,242（97.26%同一性，365 bp）并跨越原始边界。qPCR显示XY1中contig_40靶标Ct为13.82±0.09，染色体参考Ct为16.45±0.06，估算约6.2拷贝/染色体当量；FB1中靶标仅出现晚期非可重复扩增。

**XY1相关Col样元件的公共基因组筛查**
在16,456个沙门氏菌基因组中，13.8-kb ColRNAI样质粒的严格全长命中仅2个（0.0122%），均属于Derby血清型；contig_40的严格全长命中为15个（0.0911%），分布于Typhimurium、I 1,4,[5],12:i:-、Javiana等血清型。在11,036个元数据丰富子集中分别为2个（0.0181%）和12个（0.1087%）；在10,028个严格粪便–侵袭性主要子集中，质粒无命中，contig_40有7个（0.0698%）。402个Muenster基因组中两个元件均无严格全长命中。额外的106个公共记录定向raw-read分析中，8个记录达到近全长读长支持阈值，且全部为已检出的严格全长组装命中；未发现组装阴性但raw-read阳性的新候选。

### 五、讨论与结论总结

讨论部分指出，FB1–XY1配对在严格核心SNP过滤后仅保留4个SNP，而两个粪便–粪便对照无保留SNP，这与既往研究一致：同一患者分离株在核心基因组水平高度相近，但附属含量可有差异。本研究的独特贡献在于，临床关联的粪便–血液配对在近克隆状态下，最明显的差异不是核心基因组分歧，而是小Col样环状附属元件：一个组装解析的13.8-kb ColRNAI样质粒和一个组装缺失但经读长与分子验证的5.2-kb Col156样环状元件。ColRNAI样质粒带有大肠杆菌素K相关模块和蛋白酶K敏感的抑制表型，提供了可解释的附属基因组差异，但不应将其归因于特定cka样位点或视为侵袭性决定因素。contig_40的组装缺失提示短环状元件可能在DNA提取、大小选择和文库制备中丢失，需要读长水平和分子验证才能确认。公共背景筛查表明这两个元件在全库及Muenster背景中均罕见，支持其作为样本特异或稀有元件，但应视为组装检测频率而非真实流行率。研究局限性包括单一指示对、表型未定位到特定基因、contig_40功能未知、qPCR拷贝数为染色体归一化值等。总体而言，研究支持将读长证据、分子验证和公共背景相结合来解读临床微生物基因组中的小环状元件。

结论：这些发现将小型Col样环状附属元件确定为近克隆同一宿主粪便–血液沙门氏菌配对中一个经验证的、对组装敏感的基因组变异层。通过证明临床关联的粪便–血液分离株尽管核心基因组分歧极小，仍可在小型环状附属元件水平上存在差异，本研究将配对分离株解读扩展至核心SNP亲缘关系之外。更广泛地说，这项工作支持一个整合读长水平证据、分子验证和公共基因组背景的框架，用于解读临床微生物基因组学中的小型环状元件。