盐城市2012—2023年手足口病流行特征及病原学变化分析

[摘 " 要] " 目的：分析盐城市2012—2023年手足口病流行特征和病原学变化，为制定有效防治措施提供依据。方法：选取中国疾病预防控制信息系统报告数据，描述性分析2012—2023年盐城市手足口病发病及病原学检测结果，采用Cox-Stuart进行时间序列趋势检验，中断时间序列法评价EV-A71疫苗应用前后手足口病发病变化趋势。结果：2012—2023年盐城市累计报告手足口病例69 172例，年均发病率（46.87～147.63）/10万。发病高峰多位于春夏季，发病居于前3位的县（市、区）分别为亭湖区、盐都区、东台市，发病年龄主要集中在5岁及以下，男女性别比为1.56 ∶ 1。2012—2023年优势病原逐步由EV-A71型转为CVA6、16型。EV-A71疫苗接种措施实施前后EV-A71感染手足口估算发病率、手足口病重症发病率均呈下降趋势。结论：江苏省盐城市应进一步加强手足口病春夏季防控工作，重点关注5岁及以下儿童，重点防控地区为大市区，当前优势病原以CVA6、16型为主。

[关键词] " 手足口病；流行特征；病原学；EV-A71疫苗

[中图分类号] " R181.3 [文献标志码] " B [DOI] " 10.19767/j.cnki.32-1412.2024.06.024

手足口病（hand， foot and mouth disease，HFMD）由肠道病毒感染引起，常见于5岁以下儿童，以手、足、口腔等部位皮肤黏膜斑丘疹、疱疹、溃疡为典型表现，重症病例病情进展快，甚至可引起死亡[1]。2008年5月手足口病被纳入我国丙类法定传染病管理，并先后发布诊疗指南、重症病例临床救治专家共识、聚集性和暴发疫情处置工作规范，以稳步推进防控工作。手足口病在较长时间内处于盐城市法定传染病报告发病率首位，学校和托幼机构聚集性疫情频发。本研究对盐城市2012—2023年手足口病流行特征及变化趋势进行分析，以期为制定精准防制措施提供科学依据。

1 " 资料与方法

1.1 " 资料来源 " 本研究数据来源于“中国疾病预防控制信息系统”传染病监测信息系统，统计时间为2012年1月1日至2023年12月31日，内容包括地区、年龄、性别、发病时间等传染病报卡信息以及发病数、发病率。依据《手足口病诊疗指南（2018年版）》[2]，结合临床特征、流行病学调查以及实验室检测结果进行病例判定。

1.2 " 病原学监测 " 依据《手足口病预防控制指南（2009版）》[3]，采集手足口病新发病例咽拭子、粪便或肛拭子标本，采用反转录聚合酶链反应（reverse transcription-polymerase chain reaction，RT-PCR）进行EV-A71、CVA16、6、10等肠道病毒核酸检测。

1.3 " 统计学处理 " 运用Excel 2016软件录入、整理数据，SAS 9.4统计学软件进行数据统计分析。计数资料以n（%）表示，组间比较采用χ2检验；疾病时间分布、地区分布，人群分布及病原学监测结果采用描述性流行病学方法进行分析；时间序列趋势分析采用Cox-Stuart检验；有序分类变量的线性趋势分析使用Cochran-Armitage检验。P＜0.05为差异有统计学意义。采用R 4.3.2软件进行中断时间序列（interrupted time-series，ITS）分析并进行绘图。

2 " 结 " " "果

2.1 " 发病概况 " 2012—2023年盐城市累计报告手足口病例69 172例，其中重症病例275例，死亡1例。手足口病年均发病率（46.87～147.63）/10万，2012—2017年发病率基本呈“大小年”隔年高发，2018—2021年发病率逐渐降低，2021—2023年发病率逐步增高。见图1。

应用Cox-Stuart趋势检验，分析2012—2023年手足口病发病率的时间序列是否存在单调趋势。观察年份N=12，C=N/2=6，前半段为xi，后半段为xi+C，i=1，2，…，C，统计xi-xi+C的“+”“-”个数，得出6个“+”，0个“-”，故K=0，n=6+0=6。Binomial检验结果为P=0.03125＜0.05，即盐城市手足口病发病序列在12年间有显著下降趋势。

2.2 " 流行病学特征

2.2.1 " 季节分布：2012—2023年手足口病发病存在季节性高发现象（图2），季节指数法[4]分析显示，4～7月和10～12月是手足口病的高发季节，其中5、6、7月发病强度较大（图3）。

分段式圆形分布法[5]分析显示，春夏季Z=267.678，P＜0.05，差异有统计学意义，即平均角存在，x=204.37°，S=91.60°，对应高峰日为6月13日，高峰期为4月27日至7月30日。秋冬季Z=154.352，P＜0.05，差异有统计学意义，平均角存在，x=141.36°，S=94.29°，对应高峰日为11月11日，高峰期为9月24日至12月28日。见表1。

2.2.2 " 地区分布：2012—2023年盐城市各县区均有手足口病病例报告，年均报告发病率最高为亭湖区（200.58/10万），最低为滨海县（20.66/10万）。大市区（亭湖区、盐都区）病例数（36 144例）多于县（市、区）病例数（33 028例），大市区年均报告发病率（179.51/10万）高于县区报告发病率（50.19/10万），差异具有统计学意义（χ2=32 057.464，P＜0.001）。各县（市、区）发病分布情况见图4。

2.2.3 "人口学分布：2012—2023年男性手足口病报告病例42 184例（60.98%），女性报告病例26 988例（39.02%），性别比为1.56 ∶ 1。男性平均报告发病率（98.17/10万）高于女性（63.47/10万），差异有统计学意义（χ2=3 184.144，P＜0.001）。发病年龄主要集中在5岁及以下（91.79%），1岁组报告病例数最多，为17 637例，发病率为2 048.28/10万。随年龄增加报告发病率呈显著下降趋势（Cochran-Armitage检验：Z=-631.401，P＜0.001）。见表2。各人口学分类均有病例报告，以散居儿童为主，2012—2023年共报告45 833例（66.26%），其次为托幼儿童20 429例（29.53%）和学生2 584例（3.74%）。

2.3 " 病原学监测 " 2012—2023年共对5 622份手足口病病例的咽拭子/肛拭子标本进行病原学检测，检出阳性2 957份，阳性检出率52.60%，其中CVA16型占32.23%，CVA6型占30.88%，EV-A71型占17.89%，CVA10型占1.01%，其它肠道病毒占17.99%。EV-A71、CVA16、其它肠道病毒在2012—2017年均有检出，优势病原并不显著。随着EV-A71疫苗接种的不断开展，2018—2023年病原构成逐步转为CVA6、16型交替流行。见图5。

对275份重症病例标本进行病原学检测，阳性检出率56.00%，其中EV-A71型占64.29%，CVA6、16型各占9.74%，其它肠道病毒占16.23%，自2018年起重症病例优势病原由EV-A71型转变为柯萨奇病毒，且重症病例数明显减少。盐城市2012—2023年有1例手足口病死亡病例，为EV-A71型感染。见图6。

2.4 " EV-A71疫苗接种对手足口病发病的影响 " EV-A71感染手足口病估算发病率=（该年EV-A71型阳性检出率×该年手足口病报告发病数）/该年同期人口数×K，K=100 000/10万[6]。2012—2023年盐城市各年度EV-A71感染手足口估算发病率、重症发病率及ITS模型所需变量（X1、X2、X3）赋值见表3[7-8]。

使用Durbin Waston（DW）法对相关数据进行自相关检验，其中2012—2023年EV-A71感染手足口病估算发病率回归分析与重症发病率回归分析中DW统计量分别为2.46、2.80，提示不存在1阶自相关。

对EV-A71感染手足口估算发病率分析结果显示，中断线性模型方程为：Y=28.893-7.350X1+16.937X2+5.304X3+ε，在EV-A71疫苗接种措施实施前（2012—2016年）估算发病率每年下降7.350/10万（P=0.018）；措施实施当年（2017年）估算发病率瞬间变化量为16.937/10万（P=0.122）；措施实施后（2018—2022年）估算发病率斜率改变量为5.304/10万（P=0.103），斜率为-2.046（β1+β3=-2.046）。

对手足口病重症发生率分析结果显示，中断线性模型方程为：Y=1.10609-0.28952X1+0.76089X2+0.20440X3+ε，在EV-A71疫苗接种措施实施前重症发生率斜率为-0.28952，呈下降趋势，差异有统计学意义（P＜0.001）；措施实施当年重症发生率较接种前上升0.761%，差异有统计学意义（P＜0.01）；措施实施后重症发生率斜率为-0.08512（β1+β3=-0.08512），呈下降趋势，差异有统计学意义（P=0.014）。见表4、图7、图8。

3 " 讨 " " "论

盐城市2012—2023年全市累计报告手足口病例69 172例，12年间年均发病率为80.90/10万，低于全国平均水平，并呈下降趋势。2017年前全市手足口病发病基本符合“大小年”的规律，与周边浙江[9]等省及全国的调查结果类似，这一现象可能与肠道病毒感染后自身产生抗体有关。2017年后发病“大小年”现象不再明显，可能与我国在2016年全面开放二孩政策，新生儿数量增多带来易感人群扩大有关。2020年全市手足口病发病高峰大幅后移，2021年手足口病发病率为12年来最低，主要原因为2020年开始受新型冠状病毒疫情的影响，托幼机构与学校聚集性活动大大减少，居民保持社交活动安全距离，多方强化手卫生与健康宣教，在一定程度上阻断手足口病的传播。

全市手足口病发病呈较明显季节分布特征，季节指数法与分段式圆形分布法结果均显示，手足口病发病呈双峰趋势，主峰位于5～7月，次峰位于10～12月。赵霞等[10]研究显示，温湿度越高，手足口病发病率越高。盐城市属于亚热带季风气候，具有雨热同期的特点，高温炎热、多雨潮湿的季节有利于肠道病毒的繁殖与传播[11]。

从地区分布来看，大市区发病率显著高于各县（市、区）发病率。可能与大市区人口较密集、人员流动频繁、聚集性场所较多、经济水平较高、诊疗水平与传染病监测系统敏感性较强等特点有关。报告病例主要集中在5岁及以下婴幼儿群体，与散居儿童为主的人群分布特征相一致。发病人群以男性为主，可能与男童户外活动更多，个人卫生习惯养成欠佳，更易接触致病菌有关。

病原学监测结果显示，盐城市手足口病病原型别构成不断波动变化。2012—2023年优势病原由EV-A71型逐渐过渡为柯萨奇病毒，以CVA16、CVA6交替成为主要病原。2016年起盐城市在连续监测到其他肠道病毒阳性占比居高不下的情况下，对病原学检测工作进行适应性调整，逐步增加对CVA6、CVA10的分离。盐城市2017年开始EV-A71疫苗接种工作，自2018年起EV-A71型检出占比显著减少，重症病例数量显著减少，证实EV-A71疫苗接种可有效预防EV-A71型手足口病，降低手足口病的危害[12-15]。

通过中断时间序列方法分析显示，EV-A71疫苗接种措施实施前后EV-A71手足口病估算发病率、手足口病重症发生率均呈下降趋势，但措施实施前下降趋势较措施实施后更明显，这是由于在整体下降的情况下，基数越小（如接种措施实施后估算发病率与重症发生率在部分年份为0），下降幅度会随之减缓。此外，2017年夏天盐城市开始EV-A71疫苗接种，当年干预时间较短，引起估算发病率与重症发生率瞬时变化量增加现象。

本研究存在局限性，首先本研究报告发病数据来源于被动监测网络，可能存在低估实际发病情况的问题；其次在探究季节性因素对于手足口病发病的影响时，缺乏实际气象数据的支持与验证，且未对冬季第二峰的发生驱动因素进行探索[16]。

综上所述，盐城市应进一步加强手足口病春夏季防控工作，重点防控地区为大市区，重点关注5岁及以下儿童，当前优势病原以CVA6、16型为主。通过多渠道健康教育措施以强化幼童父母与托幼机构从业者掌握手足口病预防知识，应持续强化手足口病原学监测能力，以应对肠道病毒病原谱的变化。

[参考文献]

[1] FU X，WAN Z，LI Y，et al. National epidemiology and evolutionary history of four hand，foot and mouth disease-related enteroviruses in China from 2008 to 2016［J］. Virol Sin，2020，

35（1）：21-33.

[2] 国家卫生健康委员会. 手足口病诊疗指南（2018年版）［J］.中国病毒病杂志，2018，8（5）：347-352.

[3] 中华人民共和国卫生部. 手足口病预防控制指南（2009版）［J］. 中国实用乡村医生杂志，2012，19（19）：3-6.

[4] 王若琳，李亚飞，王文华，等. 2008—2023年河南省手足口病季节分布特征［J］. 现代疾病预防控制，2024，35（7）：504-508.

[5] 段陈林，周银柱，陈水连，等. 应用集中度与圆形分布法分析2012—2018年长沙市手足口病季节性特征［J］. 实用预防医学，2023，30（1）：30-33.

[6] 姚红岑，周杰，董兆鹏，等. 上海市金山区肠道病毒71型灭活疫苗接种情况及接种前后手足口病流行病学和病原学变化分析［J］. 中国生物制品学杂志，2023，36（3）：315-320，

329.

[7] 徐文杰，尹志英，郑灿杰，等. 基于中断时间序列评价EV71疫苗上市应用对衢州市手足口病报告发病情况影响［J］. 现代预防医学，2024，51（9）：1702-1707.

[8] 陈纯，李意兰，蔡文锋. 肠道病毒A组71型疫苗接种对广州市肠道病毒A组71型手足口病发病影响分析［J］. 中国病毒病杂志，2024，14（2）：130-134.

[9] 马婉婉，龚磊，肖永康，等. 2015—2022年安徽省手足口病流行病学特征及病原谱变化［J］. 安徽医科大学学报，2023，

58（10）：1763-1768.

[10] 赵霞，崔劲松，张奕，等. 2014—2018年上海市闵行区手足口病流行特征及其与气象因子相关分析［J］. 上海预防医学，2022，34（3）：219-222，234.

[11] 郑艳妮，李盛，赵玉锐，等. 武威市手足口病流行病学特征分析及ARIMA预测［J］. 疾病预防控制通报，2023，38（4）：10-14，68.

[12] 袁小亮，吴尤佳. 南通市区重症手足口病患儿早期预警指标及其评估价值研究［J］. 交通医学，2020，34（1）：80-82.

[13] 李明，朱梦洁，李海朋，等. 连云港市2015—2020年手足口病流行病学特征分析［J］. 安徽预防医学杂志，2023，29（3）： 222-224，245.

[14] DU Z，HUANG Y，BLOOM M S，et al. Assessing the vaccine effectiveness for hand，foot，and mouth disease in Guangzhou，China： a time-series analysis［J］. Hum Vaccin Immunother，2021，17（1）：217-223.

[15] 王瑶，李峰，王海燕，等. 儿童手足口病流行特征及肠道病毒核酸检测分析［J］. 交通医学，2024，38（5）：540-542.

[16] 韩春英，苟发香，张晓曙，等. 2010—2018年甘肃省手足口病季节性流行第二峰气象驱动因素分析［J］. 中国公共卫生，2024，40（2）：218-222.

[收稿日期] 2024-10-07

（本文编辑 " 缪宏建）