《Microbiology Spectrum》：Longitudinal characterization of mixed-genotype SARS-CoV-2 infections in a military cohort reveals compartmentalized viral populations
编辑推荐：
混合基因型严重急性呼吸综合征冠状病毒2（SARS-CoV-2）感染因其可能产生导致新变异株出现的新型重组病毒而备受关注。为更好地理解宿主内病毒动态，研究人员分析了来自美国军事卫生系统“大流行潜在新兴传染病流行病学、免疫学与临床特征”（EPICC）COVID-1
混合基因型严重急性呼吸综合征冠状病毒2（SARS-CoV-2）感染因其可能产生导致新变异株出现的新型重组病毒而备受关注。为更好地理解宿主内病毒动态，研究人员分析了来自美国军事卫生系统“大流行潜在新兴传染病流行病学、免疫学与临床特征”（EPICC）COVID-19队列中24名疑似混合基因型SARS-CoV-2感染参与者的标本。在最初24例疑似病例中，研究人员确认17例为真实共感染，并按证据强度进行分级：7例为“强”证据，4例为“中等”证据，6例为“弱”证据，另有7例被判定为不太可能属于真正的混合基因型感染。研究人员获取了感染过程中多个身体部位的拭子标本，从而能够观察到区室化现象和变异株优势地位的转换，这些在单时间点分析中无法被发现，同时还发现了一个重组Omicron BA.1/BA.2基因组。通过采用基于证据的生物信息学框架对来自特征明确的临床病例的测序数据进行评估，研究人员区分了真实共感染与生物信息学伪影。研究结果强调了广泛标本采集和审慎生物信息学方法在确定双重基因型感染中的重要性。
论文解读
一、研究背景与研究意义
SARS-CoV-2自出现以来持续通过快速突变进化，尤其是在编码刺突蛋白的区域。共感染（coinfection）被认为是产生新型重组病毒、进而催生新变异株的重要机制之一。早在2020年，共感染即被提出作为新变异株出现的潜在机制而受到关注。在疫情早期，Alpha、Beta、Gamma和Delta等多种关切变异株（variants of concern, VOCs）相继出现；自2021年末以来，Omicron及其后代谱系成为全球主要流行类型。
尽管文献中已对SARS-CoV-2双重基因型共感染有所探讨，但大多数研究受限于依赖单一时间点采集的单一标本的基因组数据。单次拭子采集无法解决宿主内病毒动态的关键问题，例如多种病毒谱系在感染过程中的不同时间点是否存在竞争或共存，或者呼吸道不同部位（如鼻腔和咽喉）是否含有不同的病毒群体，即区室化（compartmentalization）现象。此外，常规分析过程中可能过滤掉能够回答这些问题的数据，因为现有方法难以区分由低频变异代表的真实共感染与序列伪影。这些空白并非由于缺乏分析工具，而是由于实际限制以及缺乏捕捉单一人类宿主内SARS-CoV-2随时间变化动态的数据。
为应对上述问题，研究人员利用了EPICC研究这一独特的纵向队列数据集。该队列涵盖美国军事卫生系统（Military Health System, MHS）受益人群，其关键特征在于感染过程中采集了多种类型的标本，包括鼻咽拭子（nasopharyngeal swab, NP）、口咽拭子（oropharyngeal swab, OP）、前鼻孔拭子（anterior nares swab, NS）和直肠拭子（rectal swab, RS）。研究的主要目标是在严格区分真实共感染与序列伪影的前提下，定义混合基因型SARS-CoV-2感染的时间和空间动态。
二、主要技术方法
研究人员对EPICC队列中24名疑似混合基因型感染参与者的93份标本进行了全基因组测序。标本主要采集于2020年4月至8月，另有2名参与者的标本采集于2021年4月。测序采用扩增子平铺策略，使用ARTIC v3或Midnight引物，在Illumina NextSeq 550平台上完成；部分标本使用ARTIC v4引物在Illumina MiSeq平台进行重测序。生物信息学处理采用病毒扩增子Illumina工作流程（Viral Amplicon Illumina Workflow, VAIW）生成共有基因组，使用iVar进行宿主内变异分析，Pango分类工具和Nextclade分别进行谱系和进化枝（clade）分类。研究人员建立了一套基于证据的分级框架，综合变异分析、谱系和进化枝分类三方面证据，将混合基因型感染证据分为强、中等、弱和不可能四级。此外，使用Bolotie工具检测重组事件，并将研究结果与Freyja和Kallisto两种单样本反卷积工具进行比较。
三、研究结果
队列特征与发现总结
研究人员对24名参与者的93份标本进行了测序，这些参与者占EPICC整个队列（n=7,911）的0.3%。最终确定7例具有强证据的混合基因型感染，5例中等证据，8例弱证据，4例被判定为不太可能属于真正的混合基因型感染。共有基因组分析产生了26份编码区完整的标本和67份存在至少一个零覆盖度核苷酸的标本。Pango谱系分配结果涵盖19个独特谱系，Nextclade将全部93份标本分配到7个进化枝之一。研究面临的主要挑战之一是疫情早期许多共循环谱系共享大量定义突变，这需要基于多重证据的方法来可靠判定混合基因型感染。
纵向和多部位采样揭示混合基因型感染可能性
研究人员观察到谱系持续存在、区室化或时间性转换等复杂模式。在6例“弱”证据病例中，各时间点和各标本的谱系和进化枝均保持一致而非混合状态。4例中等证据病例主要表现出变异和谱系均不同但Nextclade进化枝相同。7例强证据病例中有5例在标本间存在10个或更多差异变异。以参与者BDRD-13为例，第0天鼻咽拭子为B.1.369谱系/20C进化枝，而后续5份标本均为B.1.2谱系/20G进化枝，表明感染过程中变异株从混合状态向单一优势谱系收敛。参与者BDRD-18的三份标本虽属于同一谱系和进化枝，但存在5个独特突变，且变异按感染天数而非拭子类型分隔，提示可能存在低水平或序贯性继发感染。
重组病毒的调查
研究人员在参与者BDRD-21中鉴定出一个单一重组SARS-CoV-2病毒，该病毒为Omicron BA.1和BA.2的重组体。该参与者采集的三份拭子标本的共有序列在单核苷酸变异方面基本一致，重组断点位于非刺突基因区域，其中ORF1a区域（8,000-13,000 nt）来源于BA.1，基因组其余部分由BA.2组成。该重组病毒的出现必然要求两种亲本谱系曾发生共感染。
NMRC框架与单样本工具的比较
研究人员将分级框架的结果与Freyja和Kallisto两种常用单来源反卷积工具进行比较。三种方法仅在24名参与者中的6名（25%）产生一致判定。NMRC框架与Freyja的一致性较高（16/24，66.7%），与Kallisto的一致性较低（7/24，29.2%）。Freyja检测到了除BDRD-24之外的所有被分级为“中等”或“强”证据的混合基因型感染。Kallisto由于依赖预定义数据库，在谱系未知的发现场景中表现不够可靠。
四、讨论与结论
本研究利用独特的纵向、多部位队列，详细描述了17名患者SARS-CoV-2共感染的宿主内动态。研究结果表明，虽然真正的混合基因型感染较为罕见，但存在区室化和谱系优势时间性转换等可观察要素，这些在单样本分析中无法显现。研究人员观察到从混合基因型感染向单一优势谱系的短期转换，以及低水平或未成功建立的继发感染。重组Omicron BA.1/BA.2病毒的发现与同期意大利报道的BA.1和BA.2共感染一致。
研究人员建立的保守分级框架强调多重证据来源，与现有单样本工具相比具有独特优势。NMRC框架适用于复杂数据集的调查研究，旨在数据质量存疑时减少假阳性；Freyja则适合作为单标本高通量筛查的敏感替代方案。方法选择涉及权衡：仅采用强证据判定最为准确但可能损失灵敏度，而将中等和强证据均纳入判定则可兼顾灵敏度和准确性。
研究存在若干局限性：各参与者的可用拭子标本数量不等，限制了参与者间的直接比较；部分病例存在低质量或低覆盖度测序数据问题；不同研究间的方法学多样性使共感染发现的直接比较复杂化。尽管如此，EPICC队列中混合基因型感染的发生频率支持此前关于其相对罕见的结论，同时证明共感染并非仅见于免疫功能低下患者或关切变异株，可发生于多种基因型组合中。区分同质性与多毒株SARS-CoV-2感染对于指导单克隆抗体治疗等医疗对策的选择、密切接触环境（如舰船或兵营）中感染聚集的调查以及军事行动准备度的维护具有重要意义。研究所产生的数据和方法不仅补充了围绕该病原体的知识基础，其基于多重遗传证据和证据优势评分判定混合基因型感染可能性的方法学思路，也可能推广应用于SARS-CoV-2之外的其他病毒感染研究。
