Epidemiological characteristics and temporal-spatial clustering analysis of hand, foot and mouth disease in Chengdu City from 2014 to 2023

doi:10.20199/j.issn.1672-2302.2025.02.006

Abstract

Abstract:

Objective To analyze the epidemiological and spatiotemporal distribution characteristics of hand, foot and mouth disease (HFMD) in Chengdu area so as to provide evidence for targeted prevention and control of this infection. Methods The data on HFMD cases and pathogen surveillance in Chengdu from 2014 to 2023 were retrieved from the Chinese Disease Prevention and Control Information System. Descriptive epidemiological methods were used to analyze the characteristics of temporal, population and regional distribution. Spatial autocorrelation analysis and spatiotemporal scan statistics were used to analyze the clustering of HFMD in time and space. Results A total of 337 633 HFMD cases were reported in Chengdu area from 2014 through 2023, during which the incidence reported ranged from 83.53/100 000 to 383.86/100 000, with an annual incidence being 195.82/100 000 on average. The prevalence peaked in April-July and November-December in general. Cases were predominant in population ≤6 years old (96.44%), and most involved in children living at home (59.35%). The annual incidence rate was 224.52/100 000 for males and 166.81/100 000 for females. Global spatial autocorrelation analysis revealed a positive spatial correlation in 2015 (Moran's I = 0.19, Z=2.03, P < 0.05). Local spatial autocorrelation identified “high-high” clusters in peripheral urban districts (e.g., Xindu, Pidu) and “low-low” clusters in core urban districts (e.g., Gaoxin, Qingyang). Spatiotemporal scan analysis indicated primary clusters in peripheral districts (e.g., Chongzhou, Dayi, Longquanyi) during May-August and October-November. Other enteroviruses (untyped) were dominated from 2017 to 2020, whereas CVA6 became predominant from 2021 to 2023. Conclusion HFMD in Chengdu area exhibited spatiotemporal clustering, primarily affecting children ≤6 years old, with CVA6 as dominant pathogen in past years. Our findings suggest that prevention efforts should focus on high-risk seasons, key areas, and vulnerable populations.

Key words: Hand, foot and mouth disease, Epidemiological characteristics, Pathogen, Spatiotemporal clustering, Chengdu City

CLC Number:

R512.5，R181.3

XU Kai, WANG Yao, SU Liyuan, FENG Songqi, CHENG Yue, LIU Hui. Epidemiological characteristics and temporal-spatial clustering analysis of hand, foot and mouth disease in Chengdu City from 2014 to 2023[J]. Journal of Tropical Diseases and Parasitology, 2025, 23(2): 91-97.

Figures/Tables 7

Table 1

Figure 1

Table 2

Table 3

Incidence rate of HFMD in different regions in Chengdu area from 2014 to 2023（/100 000）

地区	2014年	2015年	2016年	2017年	2018年	2019年	2020年	2021年	2022年	2023年
成华区	294.23	246.47	316.38	167.08	456.24	257.02	156.84	133.40	71.60	191.83
崇州市	310.63	167.08	239.15	215.63	270.10	228.72	151.49	110.61	91.46	164.58
大邑县	195.06	147.37	206.25	161.03	344.63	290.60	159.84	162.69	96.96	276.11
东部新区							49.87	96.76	62.52	134.98
都江堰市	143.78	92.65	167.31	130.13	284.20	133.86	75.27	168.92	63.06	158.65
高新区	233.88	201.80	297.56	173.20	492.91	243.18	126.86	134.65	68.34	204.20
简阳市			91.94	79.44	226.66	135.30	120.55	80.74	84.40	89.60
金牛区	268.03	202.51	284.88	147.35	305.27	164.17	109.32	111.98	64.99	160.13
金堂县	153.79	77.65	133.48	150.98	228.38	267.23	100.20	147.64	87.63	161.05
锦江区	275.17	190.80	293.82	154.98	500.59	213.70	96.61	125.50	58.14	183.41
龙泉驿区	306.20	269.92	348.06	186.57	624.38	265.32	140.55	232.77	82.99	265.14
彭州市	139.35	88.64	115.00	110.81	226.89	139.59	100.50	204.31	85.32	142.02
郫都区	229.16	170.78	312.95	198.93	410.16	257.41	160.10	142.45	90.69	179.97
蒲江县	210.61	198.69	275.07	166.47	412.05	246.36	159.99	287.95	205.82	193.56
青白江区	375.26	301.93	392.53	209.10	483.06	269.93	148.11	164.11	140.25	246.36
青羊区	253.50	207.80	326.33	163.35	380.73	193.03	107.13	114.70	66.55	166.37
邛崃市	305.34	101.54	185.46	125.02	256.19	155.09	70.93	123.37	63.94	181.90
双流区	306.13	205.14	348.55	215.77	499.15	343.34	234.23	139.32	88.37	221.25
天府新区	289.40	175.51	292.96	213.66	560.97	210.35	145.09	184.62	139.42	280.70
温江区	343.18	237.41	337.05	188.28	539.79	382.66	159.56	174.59	93.22	256.90
武侯区	230.33	153.01	199.33	113.90	272.32	136.10	100.17	103.29	57.00	171.27
新都区	194.60	204.77	358.69	206.51	599.19	249.84	267.65	133.49	86.08	194.32
新津区	241.44	163.87	246.83	152.94	228.89	100.91	98.01	138.53	109.91	166.95

Table 3

Table 4

Table 5

Spatiotemporal scanning results of HFMD in Chengdu area from 2014 to 2023

聚集年份	聚集类型	聚集时间	聚集区域	实际病例数（例）	期望病例数（例）	RR值	LLR值	P值
2014	主集群	3/1—4/30	崇州市，大邑县，温江区，新津区，郫都区，双流区，邛崃市	3 819	1 461.90	2.79	1 387.79	<0.001
	次集群1	10/1—11/30	龙泉驿区，锦江区，成华区，天府新区，高新区	2 945	1 750.28	1.74	357.56	<0.001
	次集群2	4/1—4/30	青羊区，金牛区，武侯区	1 224	751.59	1.65	127.53	<0.001
2015	主集群	11/1—12/31	新都区，青白江区，金牛区，青羊区，成华区	2 289	1 129.53	2.13	484.89	<0.001
	次集群	11/1—12/31	崇州市，大邑县，温江区，新津区，郫都区，双流区	2 061	1 158.62	1.85	301.36	<0.001
2016	主集群	6/1—7/31	新都区，青白江区，金牛区，青羊区，成华区	3 845	1 643.61	2.48	1 130.29	<0.001
	次集群	5/1—6/30	崇州市，大邑县，温江区，新津区，郫都区	2 108	1 132.43	1.91	346.53	<0.001
2017	主集群	4/1—5/31	崇州市，大邑县，温江区，新津区，郫都区	1 667	717.82	2.41	473.80	<0.001
	次集群	10/1—11/30	新都区，青白江区，金牛区，青羊区，成华区	1 696	1 043.33	1.67	180.14	<0.001
2018	主集群	7/1—8/31	龙泉驿区，锦江区，成华区，天府新区，高新区	9 916	2 411.51	4.71	7 013.03	<0.001
	次集群	7/1—8/31	郫都区，温江区，金牛区，青羊区，武侯区	7 765	2 922.47	2.90	2 950.99	<0.001
2019	主集群	5/1—6/30	金堂县，青白江区，新都区，龙泉驿区，成华区	3 178	1 365.42	2.46	921.80	<0.001
	次集群	5/1—6/30	崇州市，大邑县，温江区，新津区，郫都区，双流区	2 824	1 285.65	2.30	719.39	<0.001
2020	主集群	10/1—11/30	金堂县，青白江区，新都区，龙泉驿区，成华区	3 957	889.64	5.20	3 070.76	<0.001
	次集群	10/1—11/30	崇州市，大邑县，温江区，新津区，郫都区，双流区	3 779	834.16	5.26	2 977.78	<0.001
2021	主集群	6/1—7/31	都江堰市，彭州市，郫都区，温江区，崇州市，大邑县	3 163	1 232.59	2.75	1 116.12	<0.001
	次集群	6/1—7/31	龙泉驿区，锦江区，成华区，天府新区，东部新区，高新区	3 504	1 489.92	2.53	1 054.74	<0.001
2022	主集群	6/1—6/30	龙泉驿区，锦江区，成华区，天府新区，东部新区，高新区	1 402	426.42	3.49	721.31	<0.001
	次集群1	5/1—6/30	彭州市，郫都区，新都区，都江堰市，青白江区，金牛区	1 884	872.20	2.30	470.31	<0.001
	次集群2	5/1—6/30	大邑县，崇州市，邛崃市，双流区，新津区，温江区，蒲江县	1 576	688.98	2.41	440.64	<0.001
2023	主集群	6/1—7/31	双流区，邛崃市，蒲江县，大邑县，新津区，崇州市，温江区，天府新区	5 113	1 896.78	2.94	1 989.11	<0.001
	次集群	7/1—8/31	金堂县，青白江区，新都区，龙泉驿区，成华区	4 616	1 871.07	2.65	1 521.37	<0.001

Table 5

Table 6

References 28

[1]	国家卫生健康委员会办公厅. 关于印发手足口病诊疗指南(2018年版)的通知:国卫办医函〔2018〕327号[A/OL]. (2018-05-15)[2024-08-21]. https://www.gov.cn/zhengce/zhengceku/2018-12/31/content_5435156.htm.
[2]	张静, 李秀惠, 李丽, 等. 手足口病病原学和流行病学研究进展[J]. 中华流行病学杂志, 2022, 43(5):771-783.
[3]	齐孝旗, 欧剑鸣, 陈武, 等. 2016—2020年福建省手足口病流行特征及病原学分析[J]. 中国预防医学杂志, 2023, 24(9):960-963.
[4]	李佳霖, 刘莹莹, 于秋丽, 等. 2016—2019年河北省手足口病流行特征及时空聚集性分析[J]. 疾病监测, 2022, 37(7):941-948.
[5]	刘青, 刘露, 张态, 等. 2016—2021年云南省大理白族自治州手足口病流行特征及时空聚集性分析[J]. 疾病监测, 2023, 38(1):75-81.
[6]	蔡靖靖, 栾荣生. 2017—2022年四川省法定传染病流行特征分析[J]. 现代预防医学, 2023, 50(9):1566-1571.
[7]	刘雅琼, 童文彬, 马小珍, 等. 2017—2021年四川省手足口病流行病学和病原学特征分析[J]. 现代预防医学, 2024, 51(10):1742-1747.
[8]	徐凯, 戴映雪, 王瑶, 等. 2011—2021年成都市其他感染性腹泻流行特征及时空聚集性分析[J]. 热带病与寄生虫学, 2023, 21(1):13-19.
[9]	张艺馨, 代炳芹, 王舒姿, 等. 2019—2023年山东省其他感染性腹泻病流行特征及时空聚集性分析[J]. 疾病监测, 2025, 40(2):196-201.
[10]	梁小洁, 张智芳, 王晓欢. 中国大陆手足口病流行特征及空间聚集性分析[J]. 中国人兽共患病学报, 2023, 39(10):1009-1016. doi: 10.3969/j.issn.1002-2694.2023.00.119
[11]	顾云, 康伟俐, 娄冬华, 等. 扬州市2012—2021年手足口病流行特征分析及发病趋势预测[J]. 医学动物防制, 2024, 40(9):833-838.
[12]	马婉婉, 龚磊, 肖永康, 等. 2015—2022年安徽省手足口病流行病学特征及病原谱变化[J]. 安徽医科大学学报, 2023, 58(10):1763-1768.
[13]	陈炜. 手足口病非EV71及非CVA16的肠道病毒分子流行病学及CVA10与EV71差异性致重症机制研究[D]. 福州: 福建医科大学, 2019.
[14]	关恒云, 王春荣, 刘岚铮, 等. 济南市2012—2017年手足口病肠道病毒谱变化及流行特征[J]. 山东大学学报(医学版), 2019, 57(11):90-97.
[15]	陈纯, 李意兰, 蔡文锋. 肠道病毒A组71型疫苗接种对广州市肠道病毒A组71型手足口病发病影响分析[J]. 中国病毒病杂志, 2024, 14(2):130-134.
[16]	朱也凡, 张书志, 段春梅, 等. 新型冠状病毒肺炎疫情防控措施对上海市嘉定区法定传染病发病率的影响[J]. 疾病监测, 2023, 38(2):181-184.
[17]	黄建英, 李静静, 刘妙玲, 等. 新型冠状病毒肺炎疫情对儿童传染病报告影响[J]. 中国公共卫生, 2021, 37(10):1527-1530.
[18]	罗艳, 杜中华, 邱丽, 等. 2017—2021年海南省手足口病重复感染病例流行特征分析[J]. 现代预防医学, 2023, 50(12):2132-2136.
[19]	谢玲, 王宏卫, 刘素红, 等. 气象因素影响下中国手足口病时空演化特征研究[J]. 地球信息科学学报, 2021, 23(3):431-442. doi: 10.12082/dqxxkx.2021.200024
[20]	杨雅斯, 卢亚陵, 方莅媛, 等. 气象因素对四川省手足口病发病率的影响及预测模型构建[J]. 四川大学学报(医学版), 2020, 51(5):685-690.
[21]	刘莹莹, 苏通, 薛卫聪, 等. 2017—2019年某市气象因素与手足口病发病相关性及滞后效应[J]. 中华疾病控制杂志, 2024, 28(11):1355-1360.
[22]	张元元, 王建醒, 袁璐, 等. 山东省2017—2021年手足口病流行病学特征及病原学监测分析[J]. 病毒学报, 2023, 39(2):453-459.
[23]	刘青, 刘露, 张态, 等. 大理州手足口病的病原学特征及其时间序列预测[J]. 中华疾病控制杂志, 2023, 27(5):574-580.
[24]	董阿莉, 刘元, 安建华, 等. 2017—2020年山西省手足口病流行特征及空间分布特点[J]. 预防医学论坛, 2022, 28(6):401-404,410.
[25]	柴向琪, 王骥涛, 高瑞红, 等. 2013—2018年太原市手足口病流行特征及病原学分析[J]. 中国生物制品学杂志, 2022, 35(1):74-78.
[26]	姜晓曼, 岳勇, 陈恒, 等. 成都市托幼机构手足口病高发危险因素研究[J]. 中国预防医学杂志, 2016, 17(2):86-89.
[27]	刘竹, 安寒, 邓小华, 等. 成都市20家托幼机构手卫生管理调查[J]. 中国消毒学杂志, 2016, 33(6):562-563,566.
[28]	安寒, 刘竹, 邓小华, 等. 成都市托幼机构卫生消毒管理现状调查[J]. 中国消毒学杂志, 2014, 31(11):1187-1189.

年份	病例数（例）	发病率（/10万）	重症数（例）	重症率（%）	死亡数（例）
2014	38 018	266.12	439	1.15	1
2015	25 867	180.36	251	0.97	2
2016	40 310	256.73	503	1.25	4
2017	25 483	160.09	173	0.68	1
2018	61 591	383.86	89	0.14	0
2019	35 075	214.79	23	0.07	0
2020	22 164	133.67	6	0.03	0
2021	30 238	145.11	21	0.07	0
2022	17 632	83.53	2	0.01	0
2023	41 255	193.98	5	0.01	0
合计	337 633	195.82	1 512	0.45	8

特征	报告病例数（例）											年均发病率（/10万）
特征	2014年	2015年	2016年	2017年	2018年	2019年	2020年	2021年	2022年	2023年	合计	年均发病率（/10万）
性别
男	21 997	15 059	23 448	14 694	35 632	20 102	12 329	17 349	10 039	24 359	195 008	224.52
女	16 021	10 808	16 862	10 789	25 959	14 973	9 835	12 889	7 593	16 896	142 625	166.81
年龄（岁）
<1	2 143	2 003	3 085	1 974	6 723	2 293	1 872	2 370	1 022	3 692	27 177	2 025.26
1	11 083	7 870	12 139	7 744	23 688	9 366	7 968	9 243	3 872	10 482	103 455	6 169.84
2	10 570	6 588	9 166	5 625	13 091	7 909	4 475	6 171	3 464	7 063	74 122	4 685.49
3	8 383	5 800	8 875	5 479	9 337	7 117	4 109	5 657	3 811	7 525	66 093	5 671.56
4	3 277	1 981	3 947	2 502	4 101	3 772	1 940	3 338	2 288	5 162	32 308	2 809.90
5	1 262	811	1 517	1 104	1 944	2 016	725	1 525	1 469	2 666	15 039	1 299.43
6	492	312	658	467	907	1 096	311	763	685	1 713	7 404	602.29
7	227	130	252	193	382	547	173	341	347	979	3 571	233.39
8	121	89	178	91	282	269	129	217	218	512	2 106	132.44
9	103	60	100	65	186	157	84	172	146	353	1 426	141.93
≥10	357	223	393	239	950	533	378	441	310	1 108	4 932	3.12
职业
散居儿童	21 533	15 405	24 440	15 177	44 030	18 985	14 126	17 325	8 179	21 188	200 388
幼托儿童	15 560	9 894	14 827	9 629	15 481	14 287	7 061	11 509	8 277	16 356	122 881
学生	774	471	866	571	1 573	1 579	790	1 191	1 067	3 284	12 166
其他	151	97	177	106	507	224	187	213	109	427	2 198

年份	Moran′s I	Z值	P值	是否聚集
2014	0.08	1.15	0.13	否
2015	0.19	2.03	0.03	是
2016	0.16	1.81	0.06	否
2017	-0.00	0.43	0.30	否
2018	-0.08	-0.39	0.37	否
2019	-0.06	-0.14	0.47	否
2020	-0.05	-0.02	0.48	否
2021	-0.01	0.38	0.35	否
2022	0.06	1.10	0.14	否
2023	-0.08	-0.32	0.40	否

年份	CVA16	EV-A71	其他肠道病毒	CVA4	CVA5	CVA6	CVA10
2017	131（40.31）	45（13.85）	149（45.85）
2018	83（31.80）	1（0.38）	177（67.82）
2019	74（52.86）	1（0.71）	65（46.43）
2020	0（0）	0（0）	11（100.00）
2021	54（12.71）	0（0）	53（12.47）	38（8.94）	3（0.71）	180（42.35）	97（22.82）
2022	122（28.84）	0（0）	37（8.75）	9（2.13）	64（15.13）	173（40.90）	18（4.26）
2023	15（2.02）	0（0）	38（5.11）	135（18.17）	3（0.40）	537（72.27）	15（2.02）