《Microbiology Spectrum》：Analysis of temporal virus evolution and intra-host diversity in long-term non-progressors by bulk next-generation sequencing
编辑推荐：
本研究聚焦长期不进展者（LTNPs）这一特殊HIV-1感染群体，旨在探究其体内病毒演化是否持续以及宿主内病毒准种的多样性特征。研究团队通过对两名LTNPs进行约两年的纵向追踪，利用近全长二代测序技术，系统分析了血浆RNA和前病毒DNA的演化、多样性及CTL表位突变。结果表明，LTNPs体内HIV-1演化有限，前病毒DNA的宿主内多样性高于血浆RNA，且CTL逃逸突变呈高频、稳定状态。这些发现为理解HIV-1自然控制机制和功能性治愈策略提供了关键的病毒学新视角。
艾滋病，一种由人类免疫缺陷病毒（HIV）感染引起的疾病，至今仍是全球公共卫生的重大挑战。对于绝大多数感染者而言，若不进行抗逆转录病毒治疗（ART），病毒会持续复制，导致免疫系统功能逐渐崩溃。然而，在庞大的感染者群体中，存在一个极为特殊的亚群——长期不进展者（Long-term non-progressors, LTNPs）。他们仿佛拥有“与生俱来”的抵抗力，在无需药物治疗的情况下，能够多年（通常7-10年以上）维持正常的CD4+T细胞计数和较低（甚至检测不到）的病毒载量，病情不向艾滋病阶段发展。这些“精英控制者”为我们理解人体如何天然控制HIV、以及如何实现“功能性治愈”（即在不持续用药的情况下长期控制病毒）提供了极其珍贵的研究模型。
然而，关于LTNPs体内病毒的“真实状态”，科学界仍存在争论。一个核心问题是：在这些看似“静止”的感染者体内，病毒是真的停止了演化，还是在免疫系统的强大压力下，以一种我们尚未完全理解的、极其缓慢或隐秘的方式继续变异和进化？此外，在感染者体内，循环于血液中的病毒（血浆RNA）和整合在细胞基因组中的“病毒档案库”（前病毒DNA）之间，是否存在遗传差异？谁更具有多样性？这些差异又反映了怎样的病毒与宿主免疫系统之间复杂的“军备竞赛”？过去的研究多局限于对HIV-1单个基因（如env, pol）的片段进行分析，难以从全基因组层面窥见病毒准种（由许多遗传学上略有差异的病毒变体组成的群体）的动态全景。为了回答这些悬而未决的问题，一项发表于《Microbiology Spectrum》的研究，采用了更强大的“侦察工具”——近全长二代测序技术，对两名LTNPs进行了为期约两年的“密集观测”，旨在从核苷酸水平深入揭示病毒准种的时空演化轨迹、宿主内多样性以及关键的免疫逃逸突变特征。
为开展研究，研究人员主要运用了以下关键技术方法：首先，建立了符合严格标准的LTNPs队列，并从两名研究对象（LZ-L-02和QZ-L-06）在2017至2020年间的四个时间点采集了配对的血浆和外周血单个核细胞样本。其次，通过近全长基因组扩增和基于Illumina NovaSeq平台的二代测序，同时获得了血浆HIV-1 RNA和PBMC来源的前病毒DNA的高深度序列数据。最后，利用一系列生物信息学工具（如SmaltAlign、CliqueSNV、MEGA7等）进行序列比对、共识序列生成、系统发育树构建、遗传距离计算、宿主内单核苷酸变异分析和CTL表位突变预测，系统评估了病毒的演化、多样性及免疫逃逸特征。
研究结果
参与者的特征
研究纳入的两名LTNPs均符合长期（>10年）无ART治疗且CD4+T细胞计数稳定>500 cells/μL的标准。两人均携带保护性HLA I类等位基因（LZ-L-02: HLA-B58:01；QZ-L-06: HLA-B51:01），但不携带CCR5Δ32缺失突变。对四个时间点（T1-T4）样本的分析中，成功获得了LZ-L-02四个时间点的血浆RNA序列，以及QZ-L-06四个时间点配对的RNA和DNA序列。
系统发育和进化分析
基于测序获得的共识序列构建的最大似然系统发育树显示，两名受试者不同时间点的序列都各自形成高度支持的单系群，表明在采样期间感染的是单一亚型变异株。遗传距离计算显示，与T1相比，后续时间点的RNA和DNA序列均表现出轻微增加的遗传差异。线性回归分析表明，三个病毒序列集（LZ-L-02 RNA, QZ-L-06 RNA, QZ-L-06 DNA）的遗传差异随时间推移均无统计学显著性，提示随访期间病毒演化有限。值得注意的是，来自QZ-L-06的RNA和DNA序列在系统发育树上显示出遗传上的区室化（即分为不同的分支）。对病毒准种单倍型的分析进一步显示，血浆RNA的序列按采样时间聚类，而前病毒DNA的序列则在不同时间点相互混杂。
宿主内多样性
研究人员量化了宿主内单核苷酸变异（iSNV，频率>0.01）。结果发现，RNA序列中的iSNV数量相对较低，而DNA序列中的数量是RNA的5到10倍。RNA序列中iSNV的频率均未超过0.5（即未成为主导变异），而DNA序列在T4时间点有部分突变频率超过0.5，成为主导并可能持续存在。通过香农熵（Shannon Entropy）评估的宿主内多样性也呈现相同模式：DNA序列的多样性高于其匹配的RNA序列。在结构基因中，env基因的iSNV数量和多样性最高，gag基因最低。氨基酸突变分析显示，env基因的突变数量也最多，且RNA与DNA序列之间的突变数量存在显著差异。
env基因可变区特征
对包膜蛋白（Env）可变区（V1-V5）的分析显示，两名LTNPs血浆病毒V1区长度适中，V5区长度较短且潜在N-连接糖基化位点（PNGS）较少。与HIV-1参考株HXB2相比，两人在CD4结合位点（CD4bs）的接触区域（如内域、Loop D、β23/loop V5/β24）均未出现与广泛中和抗体（bNAb）广谱耐药相关的特征性替代突变（如K97E, A281G等）。
CTL/CD8+T细胞表位逃逸突变
根据受试者的HLA I类基因型，研究人员鉴定了最佳定义的CTL表位，并检测了其中的突变。在LZ-L-02的RNA序列中发现了8个突变表位，在QZ-L-06的RNA和DNA序列中分别发现了11个和10个突变表位。其中，IW9（LZ-L-02）和AK11（QZ-L-06）是已知的CTL逃逸表位。绝大多数CTL表位突变在所有时间点都呈现高频率（>95%）且频率变化极小。令人惊讶的是，在QZ-L-06中，RNA序列中检测到的11个突变有10个也同时存在于DNA序列中，且类型和频率高度一致，尽管两者的系统发育关系是分开的。DNA序列中还观察到更多低频突变以及少数表位内变异体频率随时间波动甚至发生回复突变（即变回野生型）的现象。
研究结论与意义
这项研究通过高分辨率的近全长二代测序，深入描绘了两名LTNPs体内HIV-1的基因组动态图谱。主要结论包括：首先，在约两年的观察期内，LTNPs体内的HIV-1病毒（无论是血浆RNA还是前病毒DNA）演化极为有限，支持了“病毒相对静止”的观点，但前病毒DNA显示出比血浆RNA更高的遗传距离和演化迹象。其次，前病毒DNA库的宿主内遗传多样性显著高于循环的血浆病毒，这主要由大量低频突变贡献，且env基因是多样性最高的区域。第三，血浆RNA与前病毒DNA序列之间存在明显的遗传差异（区室化），暗示外周血细胞中的病毒库可能并非循环血浆病毒的直接来源。第四，尽管存在显著的遗传差异，但RNA和DNA序列中却共同维持着大量高频、稳定的CTL逃逸突变，这表明这些逃逸很可能发生在感染早期、在序列发生区室化之前，是急性感染期强大CTL免疫压力选择的结果。最后，血浆病毒env基因可变区（特别是V5区）较短且缺乏已知的bNAb广谱耐药突变，提示其中和抗体敏感性得以维持。
这些发现具有多重重要意义。它们挑战了LTNPs体内“病毒完全停滞”的简单认知，揭示了即使在卓越的免疫控制下，病毒库（前病毒DNA）仍保持着更高的复杂性和潜在的演化能力。研究清晰展示了病毒在“档案库”（DNA）和“前线部队”（RNA）之间的差异演化格局，为追溯病毒起源和库的动态提供了新线索。尤为关键的是，研究揭示了“高频、稳定的CTL逃逸突变”与“持续病毒控制”可以共存，这表明LTNPs可能并非通过完全清除逃逸变体，而是通过与这些变体达成一种“平衡”来实现长期不进展。这种平衡状态，以及相关保护性HLA等位基因、有限的病毒演化、特定的env基因特征等，共同构成了LTNPs表型的多维基础。尽管该研究样本量小、观察期有限，但其高分辨率的分析方法提供了前所未有的细节。它表明，实现功能性治愈或许并非依赖于单一机制，而是多种因素共同作用的结果。这为未来旨在模拟自然控制者免疫状态的治疗策略和疫苗设计，例如靶向关键保守表位或利用特定的免疫应答模式，提供了重要的理论依据和潜在的生物标志物。
