《Revista Argentina de Microbiología》：Tracking SARS-CoV-2: Molecular epidemiology and surveillance in Health Region VIII, Buenos Aires, Argentina (2020–2021)
编辑推荐：
研究人员分析了来自第八卫生区（Health Region VIII, HRVIII）的120 421份呼吸道样本，用于检测SARS-CoV-2（severe acute respiratory syndrome coronavirus 2，严重急性呼吸综合征冠
研究人员分析了来自第八卫生区（Health Region VIII, HRVIII）的120 421份呼吸道样本，用于检测SARS-CoV-2（severe acute respiratory syndrome coronavirus 2，严重急性呼吸综合征冠状病毒2）。2020年病毒阳性率更高；旅客与居民之间的趋势不同。第一波由谱系B.1.1.33、N.3和B.1.499主导。夏季复燃包括谱系N.5以及新出现的Alpha（B.1.1.7）和Lambda（C.37）变异株。第二波显示Gamma（B.1.1.28/P.1）和Lambda共同循环，随后被Delta（B.1.617.2）取代。去中心化地区的基因组监测可追踪病毒引入和演化。
论文解读文章
一、研究背景与问题
COVID-19（coronavirus disease 2019，2019冠状病毒病）由SARS-CoV-2（severe acute respiratory syndrome coronavirus 2，严重急性呼吸综合征冠状病毒2）引起，自2019年末出现以来构成全球公共卫生挑战。阿根廷受到显著影响，布宜诺斯艾利斯省第八卫生区（Health Region VIII, HRVIII）被识别为国家疫情中心之一。2020年3月至2021年12月，HRVIII经历多波感染，反映国家及全球趋势。大流行期间出现多种SARS-CoV-2变异株，部分具有更高传播力或毒力，影响流行动态并增加防控难度。阿根廷报告多种病毒谱系以及关注变异株（variants of interest, VOIs）和关切变异株（variants of concern, VOCs），包括Alpha、Beta、Gamma、Delta和Lambda。这些变异株塑造了大流行进程，因此需要持续基因组监测以支持及时公共卫生干预。2021年初阿根廷启动全国疫苗接种，至2021年12月，超过93%成人至少接种一剂，82.4%完成全程接种。尽管疫苗显著降低疾病负担，病毒演化和区域传播异质性仍是关键问题。HRVIII位于阿根廷东部，涵盖沿海、农村和工业区域，居民超过1.34百万，经济依赖农业、制造和旅游，Mar del Plata为主要城市和文化中心。因此，研究人员开展本研究，整合流行病学与病毒基因组数据，描述该地区谱系时空分布、突变模式和突变谱，为公共卫生规划和疾病控制提供证据。
二、主要关键技术与样本来源
本研究为回顾性研究，2020年3月至2021年12月采集HRVIII 16个市镇多家医院中有急性呼吸道感染（acute respiratory infection, ARI）且符合COVID-19定义患者的呼吸道样本，样本在症状出现5天内获得。SARS-CoV-2检测采用RT-qPCR（reverse transcription quantitative polymerase chain reaction，逆转录定量聚合酶链反应）。HRVIII居民且循环阈值（cycle threshold, Ct）小于30的阳性样本考虑测序。合格样本随机选送至Proyecto PAIS、Hospital Prof. Garrahan或Instituto Malbrán-ANLIS。全基因组测序采用ARTIC方案，部分2021年样本进行刺突基因部分测序。序列提交GISAID（Global Initiative on Sharing All Influenza Data）。谱系划分采用PANGOLIN（Phylogenetic Assignment of Named Global Outbreak LINeages）和CoVSurver。系统发育分析采用MAFFT比对、BioEdit校正、IQ-TREE最大似然树、ModelFinder模型选择及1000次超快自举。贝叶斯溯祖分析采用BEAST（Bayesian Evolutionary Analysis Sampling Trees），估计tMRCA（time to the most recent common ancestor，最近共同祖先时间）、演化速率和人口历史。突变分析采用COVDB Stanford和Outbreak.info，并进行刺突蛋白结构映射。统计学采用STATA 12中的卡方检验，显著性设为p < 0.05。
三、研究结果
（一）总体检测与阳性率
研究人员分析120 421份呼吸道样本，总体SARS-CoV-2阳性率为34.7%，2020年阳性率37.4%显著高于2021年的31.9%（p < 0.001）。样本中76.2%来自居民，23.8%来自旅客。2020年旅客阳性率41.7%高于居民36.4%；2021年居民阳性率33.4%高于旅客28.0%。居民样本主要来自General Pueyrredón区（54.0%），其中97.1%来自Mar del Plata。其他贡献市镇包括Tandil、Ayacucho、Necochea和La Costa。16个市镇中位阳性率为34.5%，最低为General Guido的20.0%，最高为Lobería的42.2%。
（二）疫情波次与谱系分布
第一波始于2020年第12流行病学周，第35周达峰。2020年第51周至2021年第12周出现显著复燃，随后第二波于2021年第13周开始，第16周达峰。2020年3月实施强制预防性社会隔离（Mandatory Preventive and Compulsory Social Isolation, ASPO）延迟了阳性病例指数增长；2020年11月实施预防性强制社会距离（Preventive and Compulsory Social Distancing, DISPO）后两周出现复燃，研究人员将其归因于夏季旅游和休闲活动导致的人口流动性增加。共选取450份阳性样本测序，2020年48份，2021年402份；2020年获得22份完整基因组，2021年获得235份完整基因组，并获114份2021年刺突基因部分序列。贝叶斯系统发育分析显示八个主要簇：B.1.499、B.1.1.33、N.3、N.5、Alpha（B.1.1.7）、Lambda（C.37）、Gamma（B.1.1.28及P.1衍生）和Delta（B.1.617.2及AY亚谱系）。部分刺突测序还确认Epsilon（B.1.427）、B.1和B.1.1.519。第一波以B.1.499为主（76%），其次为N.3（19%）和B.1.1.33（5%）。2020–2021夏季复燃以N.5为主（73%），并检测到B.1.499（14%）、B.1.1.33（7%）以及早期出现的Alpha和Lambda（各3%）。第二波谱系多样性增加，Gamma最常见（58%），其次为Delta（18%）、Lambda（17%），以及较低比例的Alpha（4%）、N.5（1%）、Epsilon（1%）、B.1（1%）、B.1.499（0.3%）和B.1.1.519（0.3%）。Gamma和Lambda在第15周至第42周共循环；Delta于第31周首次检测，并在2021年末逐渐占优。Gamma亚谱系包括P.1（52%）、P.1.15（38%）和P.1.2（10%）；Delta多样性最高，含13个不同亚谱系，其中AY.103最常见（30%），其次为AY.25（25%）、AY.43（12%）和AY.100（9%）。
（三）早期谱系的本地多样化
B.1.499和B.1.1.33的贝叶斯系统发育重建显示本地多样化证据。B.1.499形成两个主要分支A和B，分支A进一步分为A1和A2。B.1.1.33及其衍生谱系N.3和N.5呈现类似结构。B.1.499分支A中，本研究若干序列与后来在Santa Fe、Santa Cruz和Tierra del Fuego报告的毒株聚集。N.5分支A中也观察到类似模式，本研究序列与随后在Santiago del Estero、Santa Fe和Santa Cruz检测的毒株聚集。这些区域簇提示本地循环病毒的独特遗传特征，并提示可能从HRVIII向其他地区传播。这些分支的替换率为4.2 × 10
?4至5.6 × 10
?4 substitutions/site/year（s/s/y）。
（四）Lambda、Gamma和Delta的系统发育结构
Lambda、Gamma和Delta呈现不同系统发育结构，具有明确分支和亚分支，反映各自时空多样化模式。Lambda分为三个分支A–C，分支A进一步分为四个亚分支A1–A4。Gamma更复杂，包含五个分支A–E和八个亚分支，其中B1和D1序列密度高。Delta系统发育结构最碎片化，包含六个分支A–F和九个亚分支。与早期谱系B.1.499和B.1.1.33不同，研究人员未发现Lambda、Gamma或Delta从HRVIII区域传播的类似模式。HRVIII序列与其他省份序列持续交错，提示全国同期循环，而非早期本地出现。多数情况下，这些变异株先在其他地区报告，随后在HRVIII检测，提示该地区在此阶段主要作为广泛传播变异株的接受者。Delta替换率最高，在E2亚分支达6.83 × 10
?4 s/s/y，高于Lambda最高6.32 × 10
?4 s/s/y和Gamma最高4.36 × 10
?4 s/s/y。
（五）刺突蛋白突变与结构映射
各谱系具有定义性突变组合以及不同频率的附加替换。Lambda序列中，ORF3a:A110S和N:A199P在HRVIII序列中分别占48.4%（15/31）和38.7%（12/31），而全球GISAID数据中约10%的Lambda基因组报告这些突变。还观察到若干全球数据库中不常见的罕见突变，作为单次出现于不同谱系个体序列。跨谱系常见突变包括D614G，与增强ACE2（angiotensin-converting enzyme 2，血管紧张素转换酶2）受体结合和增加病毒传播性相关；P681H改变S1/S2弗林蛋白酶切割位点，可能影响病毒进入。N501Y可增加ACE2亲和力，E484K和L452Q与部分逃逸中和抗体相关。部分突变如B.1.1.33中的T716I以及Lambda中的G75V、T76I和T856N目前未与已知功能效应关联。Lambda中的Δ246–252缺失、Delta中的Δ157–158缺失以及Delta中的G142D和Gamma中的L18F等氨基末端结构域（N-terminal domain, NTD）替换，可能干扰靶向NTD的抗体结合，值得进一步研究。
四、讨论与结论翻译
讨论部分指出，本研究综合分析HRVIII在大流行前两年的时空动态、基因组多样性和分子演化。HRVIII首例SARS-CoV-2于2020年3月检测到，与阿根廷其他地区初始报告一致。2020年3月至2021年12月总体阳性率34.7%，从2020年37.4%降至2021年31.9%，与Sch?nfeld等报告的全国平均38.5%可比。疫情模式与国家趋势一致：2020年第35周初峰，随后夏季复燃，2021年4月第16周出现更大第二波。ASPO早期实施可能延迟广泛传播；Larrosa等也报告封锁在早期遏制中通过降低流动性有效。限制逐渐放松后，病例自2020年8月持续增加，9月中旬达峰，第50周下降。与其他南美国家相比，阿根廷初始应对可能延迟流行峰值，并争取时间加强医疗基础设施和诊断能力。研究优势包括纳入居民和旅客呼吸道样本。2020年旅客阳性率高于居民，2021年反转，反映HRVIII社会经济和旅游特征。季节性病例激增与夏季假期及ASPO向DISPO转变等流动性增加时期一致。多数样本来自General Pueyrredón，尤其Mar del Plata，其城市集中可能影响HRVIII传播动态；低人口密度和流动性市镇阳性率较低，而农工业活动较强、人口流动性较高的Lobería阳性率较高，与Marro等一致。第一波以B.1.499为主，两个主要分支共循环，提示多次引入，并在本地多样化，替换率5.6 × 10
?4 s/s/y。第二波由Lambda和Gamma主导，Delta后期出现。Lambda于2021年3月首次检测并迅速在多个区占优；Gamma显示更长时间持续和遗传多样性；Delta在2021年末急剧增加。系统发育显示，早期谱系B.1.499和B.1.1.33具有本地多样化，并随后在其他阿根廷省份检测到相关序列；而Lambda、Gamma和Delta显示多次引入和快速多样化，HRVIII序列与其他地区交错，提示全国同期循环，该地区在第二波主要作为接受者而非传播源。刺突突变结构映射显示D614G、P681H、N501Y、E484K和L452Q等反复出现，涉及ACE2结合、病毒进入和免疫逃逸。Lambda、Delta和Gamma中的缺失及NTD替换可能干扰抗体结合。研究局限包括：许多样本流行病学和临床数据不完整；多数样本来自Mar del Plata，成功测序样本多对应2021年，存在空间和时间偏倚；有限资金限制测序能力。尽管如此，数据集提供了大流行前两年区域SARS-CoV-2流行病学和演化的代表性概览。研究结论翻译：总之，本研究揭示了Omicron出现前阿根廷SARS-CoV-2异质且快速演化的图景。研究人员的发现强调了协调和去中心化基因组监测对于实时检测、监测和应对变异株出现的重要性。作为一个以高连通性和流动性为特征的地区，HRVIII说明了去中心化环境如何能够有意义地促进对国家传播动态的理解，并应被积极纳入大流行应对策略。
