《Microorganisms》：Integrating Phenotypic Testing, Targeted PCR and Long-Read Genomics in Staphylococcus aureus: Defining the Scope of Genotype–Phenotype Concordance
编辑推荐：
表型抗菌药物敏感性试验（antimicrobial susceptibility testing, AST）、靶向PCR及全基因组测序（whole-genome sequencing, WGS）可提供关于金黄色葡萄球菌（Staphylococcus aureu
表型抗菌药物敏感性试验（antimicrobial susceptibility testing, AST）、靶向PCR及全基因组测序（whole-genome sequencing, WGS）可提供关于金黄色葡萄球菌（Staphylococcus aureus）的互补信息，但其结果并不能直接互换。研究人员在一组来自门诊患者的450株鼻咽及咽喉部金黄色葡萄球菌分离株（2016–2020年）诊断性菌株库中整合了上述三种方法。采用EUCAST 2019判定标准对17种抗菌药物进行AST。对10株经刻意挑选的分离株（4株耐甲氧西林金黄色葡萄球菌MRSA，6株甲氧西林敏感金黄色葡萄球菌MSSA）进行长读长WGS，并针对MEGARes及VFDB数据库进行筛查；将靶向PCR结果与其基因组进行比较，并应用于40株随机选取的分离株。苄青霉素及头孢西丁耐药率分别为86.0%及38.9%；三个及以上类别内的耐药（30.4%）仅限于MRSA。MEGARes返回31条耐药相关条目，其中13条存在于全部10个基因组中。100个可直接评估的分离株-药物配对中，一致性为92.0%（bootstrap 95% CI 85.0–98.0%），不一致结果仅限于4株分离株；氟喹诺酮耐药位于可评估范围之外，且未发现典型的靶位点替代。VFDB鉴定出75个毒力相关基因；lukF-PV在4株分离株中被检出，但其中3株实际为lukD，仅1株携带PVL位点。在180个配对观察中，PCR与基因组位点匹配179个，尽管18个靶标中有10个为不变位点。基因型-表型一致性应在明确界定的分析范围内进行评估。
论文解读
一、研究背景与研究意义
金黄色葡萄球菌（Staphylococcus aureus）是一种高度适应性的人类共生菌和机会性致病菌，可定植于黏膜及皮肤表面，同时也能引起从浅表感染到肺炎、菌血症、心内膜炎等多种侵袭性疾病。耐甲氧西林金黄色葡萄球菌（methicillin-resistant S. aureus, MRSA）是全球最重要的耐药细菌病原体之一，持续加重抗菌药物耐药的全球负担。前鼻孔是金黄色葡萄球菌的主要生态储库，咽部及其他解剖部位也可发生定植。定植群体具有流行病学意义，因其参与传播并可能成为后续内源性感染的来源。
传统上，金黄色葡萄球菌的耐药性和毒力通过互补的湿实验室方法进行表征。表型药敏试验（antimicrobial susceptibility testing, AST）提供临床可解释的表达耐药信息，而常规及多重PCR可快速检测预设的耐药和毒力决定因子。然而，这两种方法均有固有局限：AST不能直接鉴定耐药性的完整遗传基础，PCR则局限于预设靶标，对等位基因多样性、基因组背景或额外非预期决定因子提供的信息有限。全基因组测序（whole-genome sequencing, WGS）通过同时表征耐药组、毒力组、克隆背景及移动遗传元件提供了更广泛的框架。长读长测序进一步促进从头组装，并改善对质粒、重复区域及结构复杂位点的解析。
然而，基因型与表型不可互换。耐药相关基因的存在并不总是导致表型耐药，而计算机检测可能受序列质量、组装、等位基因变异、数据库组成及分析阈值的影响。近期对mecA阳性但表型头孢西丁敏感的MRSA分离株的长读长分析进一步表明，SCCmec内部的结构变异可能构成基因型-表型不一致的基础，而常规PCR无法解决这一问题。这些发现支持一种整合性解读：表型、靶向分子检测和全基因组分析提供互补而非等同的信息。
二、研究设计与主要结论
本研究对2016至2020年间从门诊上呼吸道感染患者采集的鼻咽和咽喉金黄色葡萄球菌分离株，在三个分析层面进行了表征：对450株分离株进行表型药敏检测，对50株进行靶向PCR检测，对10株进行长读长全基因组测序。研究目的包括：(i) 表征由MRSA和MSSA分离株组成的测序子集的耐药组、毒力组及克隆背景；(ii) 预先为每种受试抗菌药物界定其主要耐药机制是否属于获得性基因筛查的分析范围，并在该声明范围内评估基因型-表型一致性；(iii) 评估靶向PCR panel与同一批分离株组装基因组的一致性，随后应用于独立菌株集以确定标记频率；(iv) 为保加利亚诊断菌株库中的鼻咽和咽喉金黄色葡萄球菌分离株提供试点基因组参考点。
主要结论包括：表型检测报告在限定条件下表达的耐药性，PCR报告预设序列的存在，基因组筛查报告参考数据库内容在组装中的回收情况。在声明范围之外，这些方法并非相互矛盾而是回答了不同问题；在范围之内，100个观察中有8个不一致，基因组数据能够描述但无法解释这些不一致。
三、关键技术方法
研究人员采用了三种互补的分析方法：表型AST应用于450株分离株，涉及17种抗菌药物，按EUCAST v9.0判定标准重新解读；长读长牛津纳米孔（Oxford Nanopore Technology, ONT）测序应用于10株经刻意挑选的分离株（4株MRSA、6株MSSA），采用Flye v2.9.2进行从头组装，并以MEGARes、CARD、ResFinder及NCBI AMRFinderPlus数据库筛查耐药决定因子，以VFDB数据库筛查毒力相关基因；靶向PCR panel包含20个毒力相关标记及mecA，应用于10株测序分离株进行PCR-WGS一致性评估，并应用于另外40株未测序分离株以确定标记频率。菌株来源于2016至2020年间保加利亚普罗夫迪夫SMDL "Synevo Bulgaria" Ltd.门诊上呼吸道感染患者的鼻咽和咽喉拭子。
四、研究结果
3.1 测序分离株的基因组特征与克隆背景
10株临床分离株经长读长测序表征，染色体基因组大小为2.77至2.92 Mb，G+C含量为32.82%至32.99%。多位点序列分型（multilocus sequence typing, MLST）将分离株归入6个序列型（sequence type, ST）：ST30三株，ST5和ST398各两株，ST15、ST22和ST88各一株。spa分型解析出9种类型。SCCmec元件在全部4株MRSA中检出，6株MSSA均未检出，头孢西丁表型分类与SCCmec携带完全一致。
3.2 耐药组
3.2.1 表型药敏结果
450株分离株中，苄青霉素耐药率最高（86.0%），头孢西丁耐药率38.9%，未观察到头孢他林耐药。红霉素耐药率26.2%，克林霉素16.9%。三个及以上类别内耐药占30.4%，全部为MRSA。万古霉素和利奈唑胺耐药率各为0.7%，均视为临时结果。
3.2.2 全基因组测序检出的耐药决定因子
MEGARes筛查返回31条耐药相关条目，13条存在于全部10个基因组中，18条在分离株间变异。β-内酰胺酶操纵子（blaZ、blaR、blaI）在9株中检出；mecA在4株中检出，与头孢西丁耐药及SCCmec携带完全对应。质粒定位分析显示β-内酰胺酶操纵子在316M和726D中仅存在于染色体外contig，在670M中为四拷贝（三拷贝染色体、一拷贝质粒）。
3.2.3 基因组决定因子与表型药敏的一致性
在100个可直接评估的分离株-药物配对中，基因型与表型一致率为92.0%（bootstrap 95% CI 85.0–98.0%），Cohen's κ = 0.81。8个不一致观察中，2个为表型耐药但无可预测决定因子（R/?），6个为携带预测决定因子但无表型耐药（S/+）。不一致仅限于4株分离株。β-内酰胺类完全一致。6个S/+不一致涉及ermA、ermT、ant(4′)及aac(6′)-Ie–aph(2″)-Ia，均为完整开放阅读框。2个R/?不一致为316M红霉素耐药但无erm基因，205M庆大霉素耐药但无氨基糖苷修饰酶基因。氟喹诺酮耐药位于可评估范围之外，4株耐药分离株的GyrA和GrlA QRDR均为野生型。
3.3 毒力组
3.3.1 全基因组测序检出的毒力相关基因
VFDB筛查鉴定出75个毒力相关基因，48个存在于全部10个基因组中，27个在分离株间变异。lukF-PV在4株中检出，但仅670M在99.80%同一性且紧邻lukS-PV处携带真正的PVL位点，其余3株的低同一性hit（82.06–82.37%）经系统发育分析确认为lukD。荚膜位点中，cap8H–cap8K在5株中缺失，对应荚膜5型。β-溶血素基因在7株携带scn的分离株中覆盖率为81.9–83.0%，在3株无免疫逃避基因的分离株中为99.9–100%，对应φSa3前噬菌体整合位点。
3.3.2 靶向PCR与测序基因组的一致性
18个可评估靶标产生180个配对观察，PCR与基因组结果一致179个。唯一不一致为316M的sea阳性，该基因组携带相关肠毒素基因sep而非sea。8个变异靶标贡献45个阳性及35个阴性观察，PCR结果全部与基因组发现匹配。
3.3.3 未测序分离株中靶向标记的频率
40株未测序分离株中，mecA在全部15株表型MRSA中检出，25株MSSA均未检出。7个标记在全部40株中检出：clfB、icaA、isdA、tarS、hla、hlgA和psmβ1。超抗原和肠毒素标记频率最低：tsst-1为30.0%，sec为27.5%，sea为22.5%，seb为15.0%。
五、讨论与结论
讨论部分指出，测序分离株的克隆多样性描述了panel本身而非整个菌株库。苄青霉素耐药率86.0%和头孢西丁耐药率38.9%反映了诊断菌株库的组成而非健康携带者中的MRSA定植流行率。耐药组分析显示，13个不变条目代表物种核心基因组组分而非获得性决定因子；aac(3)因保守染色体位置被排除为获得性预测因子；dfrS1仅在CARD和NCBI AMRFinderPlus中检出，说明数据库组成影响hit list。长读长组装使Tn554、pUB110等已知元件的物理排列可恢复。基因型-表型一致性92.0%在预先界定的分析范围内评估，该范围声明了哪些药物可通过获得性基因筛查评估。毒力组分析表明，lukF-PV的低同一性hit对应lukD，荚膜5型通过cap8H–cap8K缺失模式识别，φSa3整合解释了hlb的部分覆盖。PCR与WGS在179/180个配对观察中一致，但该一致性仅证明检测方法扩增了目标序列，不涉及表达。
研究结论为：表型检测报告在限定条件下表达的耐药性，PCR报告预设序列的存在，基因组筛查报告参考数据库内容在组装中的回收情况。在声明范围之外，这些方法并非相互矛盾而是回答了不同问题；在范围之内，100个观察中有8个不一致，基因组数据能够描述但无法解释这些不一致。基因型-表型一致性应在明确界定的分析范围内进行评估。
