《Medical Sciences》：Multilocus Variation in Phylogenetic Support for HIV-1 Subtype B Assignment in an Institutional Cohort from Central Mexico
编辑推荐：
背景/目标：本研究在来自墨西哥中部的机构队列中，评估了pol、gag、vif、tat和vpr中临时HIV-1（1型人类免疫缺陷病毒）B亚型兼容性（provisional HIV-1 subtype B compatibility）的节点支持（nodal sup
背景/目标：本研究在来自墨西哥中部的机构队列中，评估了pol、gag、vif、tat和vpr中临时HIV-1（1型人类免疫缺陷病毒）B亚型兼容性（provisional HIV-1 subtype B compatibility）的节点支持（nodal support）变异。方法：研究人员对60名未接受治疗的成人进行了二次描述性生物信息学再分析（secondary descriptive bioinformatics reanalysis）。58份样本有配对FASTQ（测序读长文件）文件，其中42份在全部五个区域被保留。稳健支持（robust support）要求SH-aLRT（Shimodaira–Hasegawa近似似然比检验）≥80且UFBoot（ultrafast bootstrap，超快自举）≥95。其他分析改变了参考组成阈值、重新选择分配节点、统一参考面板并子采样比对列。结果：在原始记录的节点上，18/210份样本–区域记录（8.6%）为稳健支持；25/42份样本无稳健区域，16份有一个，1份有两个。在60%组成阈值下进行节点重选保留了全部18条记录，并新增S42/pol和S29/tat，得到20/210（9.5%）。两种分析均未在三个或更多区域中识别出稳健支持。使用共同参考面板时，pol与vif之间情境性内部支支持（contextual internal-branch support）的排序发生变化。结论：支持取决于所评估的节点和分析数据集。固定节点阈值稳定性并未在节点重选下确立稳定性。区域分析之间的差异不能独立归因于基因组区域。单片段分配仍为临时性，不能确立全基因组亚型分类。V3和ViReMa分析为探索性分析。
研究背景：HIV-1（1型人类免疫缺陷病毒）遗传多样性使分子分类解释复杂化；复制、重组和选择压力造成基因组区域差异，一个片段所得分类不必然代表整个病毒基因组。分子监测需区分序列在树中的位置与支持该分组的证据。B亚型在拉丁美洲占优势，但与其他谱系和重组型共存；墨西哥城都市区已描述B亚型、BG和BF重组及异质系统发育结构。比较基因组区域有助于认识单片段分配局限，但比较亚型标签与比较支支持回答不同问题：两个区域可能均拓扑兼容（topological compatibility，指样本相对于参考序列在树中的位置兼容）B亚型，却提供不同统计支持。本研究选择pol（酶功能）、gag（结构）、vif和vpr（附属）、tat（调控）区域，因检测设计和五区域共同样本可得；功能多样性提供框架但不证明基因功能解释支持。V3仅部分恢复，未作为主要分析等效区域。需区分拓扑兼容与节点支持（nodal support，指特定系统发育节点获得统计支持的程度），阈值达标须指支持分配的节点，而非邻近分支或树整体稳健性；区域差异本身不确认亚型间重组或重组型。前40份样本曾用于耐药和亚型多样性研究；本再分析以独立流程重处理可用reads，生成共识序列（consensus sequence，由比对结果生成的共有序列）和位点特异性系统发育重建，不设计比较DeepChek软件2.0版性能。研究目的是描述同一42份样本在五区域中临时B亚型分配节点支持变异，并确定一个区域稳健支持是否也见于其他区域。研究人员开展二次描述性生物信息学再分析，得出支持取决于评估节点和分析数据集、固定节点阈值稳定不等于节点重选稳定、区域差异不能独立归因于基因组区域、单片段分配仍为临时且不确立全基因组分类等结论；其意义在于强调报告所选节点、参考组成和两个支持值，不能将单片段临时分配视为确定性全基因组分类。
主要关键技术方法：队列来自墨西哥州Ixtapaluca的Hospital Regional de Alta Especialidad de Ixtapaluca 2024—2025年横断面研究，60名未治疗HIV-1成人，非概率便利抽样；58份有配对FASTQ，42份五区域均满足保留标准。研究人员行二次描述性生物信息学再分析：质控后去人源，比对HXB2，ViralConsensus生成共识序列，MAFFT位点特异性比对，IQ-TREE 3.1.2最大似然重建，ModelFinder选模型，1000次SH-aLRT和1000次UFBoot；稳健支持为SH-aLRT≥80且UFBoot≥95；并行节点重选、共同参考面板及343列长度子采样敏感性分析；V3与ViReMa为探索性分析。
研究结果：
3.1 队列特征（Characteristics of the Cohort）：通过队列描述性统计，60例中48男（80.0%）、12女（20.0%）；平均年龄38.9±12.0岁，中位35.5岁（IQR：30.8–46.5）；中位病毒载量392,143 copies/mL（IQR：187,422–1,584,822）。58/60有配对FASTQ，42/58纳入共同系统发育数据集，16/58未纳入；V3为34条记录的二次数据集。纳入与排除样本的描述性比较不证明等效或排除选择偏倚。
3.2 系统发育矩阵和序列质量（Phylogenetic Matrices and Sequence Quality）：通过五区域矩阵与质量指标汇总，五棵树均含相同42份样本；序列总数131（vpr）至356（gag），比对长度343至2844列；简约信息位点（parsimony-informative site，系统发育分析中具信息量的变异位点）177（vpr）至1409（pol）；完整面板信息位点数与稳健内部支比例的描述性Spearman系数为0.600，与稳健样本分配比例为0.791；队列-only信息位点为pol 727、gag 448、vif 213、tat 75、vpr 118；平均完整度96.29%（gag）至98.26%（pol），平均深度5648.8×（vpr）至10,083.2×（gag）。这些描述不替代覆盖广度和平均深度的保留标准，也不确立个体排除原因。
3.3 内部支支持的情境分析（Contextual Analysis of Internal-Branch Support）：通过统计带两种支持值的内部支（internal-branch support），pol 337、gag 339、vif 270、tat 195、vpr 125；同时达阈值比例最高为vif 37.41%，最低为tat 19.49%。分母为标记内部支而非参与者，说明情境树支持需与分配节点支持分开。
3.4 样本分配的记录支持（Recorded Support for Sample Assignment）：通过记录节点支持矩阵，210份样本–区域记录中18条达双阈值，189条支持不足，3条未解决；稳健记录比例在pol和vif均为16.7%，gag为4.8%，tat和vpr各为2.4%。S19、S20、S34在vpr未解决；支持不足不确立替代亚型。18条稳健记录参考B比例均为100%，固定节点六阈值下均保持稳健；其他类别可变化。节点重选在50%、60%、70%、80%、90%、100%阈值下分别得22、20、20、20、20、19条稳健记录；60%时保留原18条并新增S42/pol（重选节点SH-aLRT/UFBoot为82.3/95，原节点25.4/82）和S29/tat（97.1/98，原63.1/80）；S34/vpr从未解决变为支持不足但未增加稳健记录。结论是支持依赖评估节点和数据集。
3.5 跨位点稳健性一致性（Agreement in Robustness Across Loci）：通过每样本跨五区域双阈值计数，42份样本中25份（59.5%）无任何位点达双阈值，16份（38.1%）一个位点，1份（2.4%）两个位点；无样本在三个或更多区域达双阈值。S39是唯一在两个位点（pol和gag）稳健的样本；pol的7份与vif的7份不重合。60%节点重选下23份无稳健区域、18份一个、1份两个，S39仍为唯一两个稳健区域者；仍无三区域或更多。该结果描述稳健性一致性，不必然等同分配一致或不一致；其他区域支持丧失本身不构成替代分配；探索性相关未见一致模式。
3.6 二次分析（Secondary Analyses）：通过V3基因型规则和ViReMa读支持分析，V3与ViReMa均为探索性。34条V3源记录中28条记为R5、6条记为X4，不解释为验证性嗜性流行率；ViReMa信号为读支持，不确认生物学断点或重组；二者未用于解释系统发育支持变异。
3.7 参考面板和比对长度敏感性（Reference Panel and Alignment Length Sensitivity）：通过共同参考面板重建和343列长度匹配子采样，五区域全长共同面板重建均达UFBoot收敛；稳健标记内部支为pol 33/79（41.77%）、gag 27/79（34.18%）、vif 16/79（20.25%）、tat和vpr各18/79（22.78%）；pol与vif差异从原面板?4.77百分点变为共同面板+21.52百分点，属不同树间内部支摘要的描述性比较。长度匹配共40次尝试（pol、gag、vif、tat各10个343列无放回子样本；vpr全长已343列作比较），39次完成，1次pol达900 s限时；31/39完成运行达UFBoot收敛，7个最终相关<0.99，1个缺最终系数；描述性均值含全部完成运行为pol 15.19%（10.13–18.99%）、gag 21.52%（16.46–26.58%）、vif 15.95%（11.39–20.25%）、tat 14.94%（10.13–22.78%）。结论是内部支支持随参考面板和子采样变化，探索性控制不分离可归因生物学区域效应。
讨论部分总结：作者强调拓扑兼容B亚型不必然伴随稳健节点支持；记录节点分析为18/210，节点重选为20/210，均无三区域或以上稳健支持。固定节点阈值变化稳定不能推广至节点重选；支持与所评估二分体（bipartition，树中由一条分支定义的类群划分）相关，重根（rerooting）须保留该关联。历史B亚型优势及多样谱系和重组提供背景；前40份样本研究强调耐药和亚型多样性，本再分析问题不同，且不评估工作流一致性或性能。部分基因组分类局限已被认识；本研究以共同操作框架比较同42份样本五区域记录节点支持，区分拓扑兼容与稳健节点支持，但不确立确定性全基因组分类。vif内部支稳健比例最高37.41%，而分配仅7/42；pol为7/42、gag为2/42、tat和vpr各1/42，说明总体支支持不能替代分配节点评价。完整面板与队列-only信息位点描述不同序列集；支支持可依赖信息字符、分类群采样和模型；共同面板与长度子采样量化描述性敏感性，不识别独特生物区域效应，不将支持变异百分比因果归因于区域、长度或参考数；基因功能不独立解释支持。pol耐药监测意义在于部分片段结果可能不同于其他区域或更长序列，分配及节点支持应在所检区域解释，全基因组意义需额外证据。18条稳健记录是样本–区域组合而非18名独立参与者；25/42无区域达阈值不证明非B；报告评估节点、参考组成和两支持值使临时分配基础与范围明确。质量记录不能分离读深、覆盖广度、模糊、参考组成和信息内容贡献；无稳健支持不是重组证据。V3共识调用非表型确认；ViReMa读支持信号不确立生物断点或亚型间重组，二者未用于解释支持差异。限制包括单机构非概率便利抽样、五区域完整样本可能选择偏倚、描述性比较不排除选择偏倚、60%阈值为操作性、历史节点选择规则未独立重建、共同面板和长度控制不验证确定性全基因组分类或分离所有变异源、未收敛运行标记并纳入描述均值、相关和控制摘要为描述性不建立统计或因果关联、执行记录不完整限制独立重建、V3不捕获样本内全部分布且模糊残基限制基因型调用、150-nt reads和聚合输出不能验证生物重组。优势为明确分母、同一样本五区域比较、保留支持不足和未解决记录、报告节点重选变化及修订材料改善追溯。
研究结论翻译：与B亚型的拓扑兼容不必然与稳健节点支持一致。记录节点分析识别18/210条稳健样本–区域记录，而60%组成阈值下节点重选识别20/210条。两种分析均未识别在三个或更多区域有稳健支持的样本。情境内部支支持也随参考面板和比对子采样而变化。这些发现描述临时、区域特异性分配及对分析选择的敏感性；它们不建立独立生物学区域效应或确定性全基因组亚型分类。报告所选二分体、参考组成和两个支持值对解释每个分配至关重要。
