超深度HiFi测序揭示左心室辅助装置感染相关的多重耐药铜绿假单胞菌中罕见结构变异亚群

《Current Research in Microbial Sciences》：Ultra-deep HiFi sequencing reveals rare structural-variant subpopulations in multidrug-resistant Pseudomonas aeruginosa from an LVAD-associated infection

　　

超深度HiFi测序解析了与左心室辅助装置（LVAD）感染相关的多重耐药（MDR）铜绿假单胞菌中罕见的结构变异（SV）异质性。罕见的结构变异涉及噬菌体、菌毛、转运蛋白和应激反应相关基因座。耐药相关的小型变异与MDR表型的组成部分一致。标准深度的全基因组测序（WGS）难以有效采样大多数低频结构变异信号。长读长测序揭示了器械相关感染中隐藏的多样性。

研究背景、现存问题与研究目的

铜绿假单胞菌（Pseudomonas aeruginosa）是一种机会性革兰氏阴性病原体，在医院环境中构成重大威胁，尤其易定植于医疗器械如心室辅助装置（VAD）。VAD相关感染难以根除，常需长期抗生素抑制治疗，且多重耐药（MDR）菌株的增多严重限制了治疗选择。结构变异（SV），包括缺失、重复、倒位、插入和染色体重排，是驱动该菌适应性、持久性和致病性的重要遗传机制。然而，由于生物学、解释和技术层面的限制，SV的研究十分困难。标准深度的短读长测序无法可靠检测大型、复杂或低频率的SV，尤其是在重复序列或移动元件富集区域，导致与耐药性、噬菌体相互作用、生物膜形成等临床相关的变异可能被遗漏。高保真度（HiFi）长读长测序结合了单分子精度和足以解析复杂基因组区域的读长，有望克服这些局限。

本研究假设，超深度长读长测序能够揭示临床铜绿假单胞菌分离株中传统测序方法无法检测到的低频SV，这些SV可能是慢性感染背景下表型变异的一个未被充分认识的来源。为此，研究人员对从一名82岁男性左心室辅助装置（LVAD）感染中分离的多重耐药菌株NRD619进行了超深度HiFi测序，旨在解析隐蔽的结构异质性、注释与已知耐药表型相关的共有水平突变，并评估HiFi组装与混合组装的差异。研究的核心目标是展示深度、高精度测序在揭示临床相关病原体中显性和隐性基因组变异方面的能力。

主要关键技术方法

1. 1.
   细菌培养与DNA提取：将冷冻保存的NRD619菌株在LB琼脂上复苏，挑取单菌落进行液体培养，采用酚-氯仿法提取高分子量基因组DNA。
2. 2.
   文库制备与超深度测序：使用PacBio SMRTbell prep kit 3.0制备文库，并在PacBio Revio平台上利用环化共有测序（CCS）技术进行单细胞测序，获得超深度HiFi读段。
3. 3.
   基因组组装与比较：使用Flye软件对HiFi读段进行de novo组装，并通过Circlator进行环化，得到完整的单重叠群基因组。将该HiFi组装与先前发表的NRD619杂交组装（Illumina/PacBio CLR）进行比对，评估一致性。
4. 4.
   变异检测与注释：
   • •
     小型变异：将HiFi读段比对至PAO1参考基因组，使用Clair3软件调用小型变异，并用Ensembl VEP进行功能注释。
   • •
     同质性与异质性SV：将读段比对至自身组装基因组，使用Sniffles2软件以默认参数调用同源性SV，并以逐步降低最小等位基因频率（--mosaic-af-min）的镶嵌模式调用异质性SV。对低频SV信号进行IGV人工检查和点阵图分析。
5. 5.
   覆盖度需求计算：基于观察到的变异等位基因频率（VAF），使用公式计算在不同概率下至少采样到一个携带SV分子的所需测序深度。

研究结果

超深度测序与de novo基因组组装将NRD619解析为一个完整、单重叠群的基因组**

通过对预处理临床分离株NRD619进行超深度HiFi测序，研究人员成功将其基因组de novo组装成一个长度为6,438,253 bp的单重叠群，经抛光环化后获得了一个完成级、参考质量的基因组。将该HiFi组装与先前基于Illumina/PacBio CLR技术的杂交组装进行比较，发现两者长度完全相同，仅在四个单核苷酸多态性（SNP）位点存在差异，表明两种组装高度一致（质量值QV > 60）。

耐药相关的小型变异与NRD619的表型耐药谱部分一致

研究人员分析了与抗菌药物耐药性相关的基因中的突变，并与表型药敏谱进行了关联。研究发现，mexZ基因的早期移码突变（E21X）可能导致MexXY-OprM外排泵去抑制，这与对氟喹诺酮类和头孢吡肟的耐药性一致。oprD基因的移码突变可能减少碳青霉烯类药物的摄取。而ampC的T105A替换和parE的D533E替换则被认为不具耐药性或作用不确定。这些变异为MDR表型提供了文献支持的遗传学背景，但并非完整的因果耐药组。

超深度HiFi测序捕获了低频率的异质性结构变异

通过逐步降低Sniffles2的最小等位基因频率阈值，研究人员回收了五个异质性SV信号。其中，SV2（重复）由95条读段支持，观察到的VAF为1.2119%；SV1（缺失）由3条读段支持，VAF为0.0383%；SV3、SV4和SV5（均为插入）各由一条读段支持，VAF约为0.0122%-0.0125%。由于事件间未进行分型，这五个SV不一定代表五个不同的细菌亚群，但其极低的频率说明了为何常规测序深度难以捕获此类稀有结构状态。

异质性结构变异的潜在表型效应

研究人员评估了每个异质性SV的潜在生物学和临床意义。

• •
  SV1/SV2：影响同一约32 kb的区域，包含42个基因，主要为噬菌体相关基因（如尾部蛋白、裂解酶、重组酶、末端酶等）。其互斥的缺失/重复架构提示可能存在噬菌体介导的重组过程。
• •
  SV3：一个3,475 bp的插入，破坏了pilQ基因，同时提供了pilO和pilP的额外拷贝。这可能影响IV型菌毛的组装、运动性和生物膜形成。
• •
  SV4：一个2,987 bp的插入，破坏了yhiL开放阅读框，该基因编码一个ABC转运体。其潜在效应尚不清楚，但可能改变转运相关表型。
• •
  SV5：一个4,175 bp的插入，破坏了PA5339（RidA同源物）并引入了oxyR等五个基因的额外拷贝。PA5339的破坏可能影响运动和生物膜形成，而oxyR剂量的增加可能增强细胞对氧化应激的反应。

覆盖低频率SV分子的测序深度需求

基于观察到的VAF，研究人员计算了以不同概率采样到至少一个SV支持分子所需的测序深度。结果显示，除了SV2，其余所有SV都需要超过1,800倍覆盖度才能有50%的概率采样到一个支持分子。例如，对于VAF为0.0122%的SV3，需要5,667倍覆盖度才有50%的采样概率。

讨论与结论

讨论部分总结：本研究强调了超深度HiFi测序在检测临床分离株内基因组变异方面的强大能力。HiFi组装与杂交组装的高度一致性验证了其准确性，而HiFi技术在单分子精度上的优势使其在分析低频变异时尤为突出。研究中发现的耐药相关小型变异为部分耐药表型提供了合理解释，但并非全部。超深度测序揭示的低频SV（最低VAF约0.0122%）凸显了常规WGS的局限性，因为临床上重要的耐药谱系可能源自稀有细胞并在选择压力下扩增。尽管SV3-SV5仅由单条分子支持而被视为推定事件，但它们影响的基因座（如菌毛、转运体和应激反应）具有明确的生物学意义，为后续功能研究提供了可验证的假设。研究也指出了局限性，包括仅分析单个治疗前分离株、无法确定SV的体内起源（可能在体内存在或在实验室传代中出现）、以及缺乏表型和转录组验证。

研究结论翻译：总之，本研究揭示了临床铜绿假单胞菌分离株中被常规测序采样极限所掩盖的基因组重排和罕见结构异质性。超深度长读长测序揭示了一个影响噬菌体相互作用、宿主适应、表面行为和医疗器械持久性相关基因座的SV携带基因型的嵌合体。单克隆的培养历史限制了对个体稀有事件体内起源和频率的推断，但这并不否定其相关性：相同的结构状态可能在体内存在，而培养相关的多样化本身也可能改变用于基因组和表型表征的群体。这些发现挑战了共识基因组能完全代表生物学相关多样性的假设，并支持将高分辨率基因组监测作为研究慢性和器械相关感染中微生物复杂性的工具。