湖北省水痘发病情况与气象因素的关系研究

罗玲，黎明

水痘是由水痘-带状疱疹病毒（VZR）感染引起的急性传染病，是儿科临床常见病［1］。其特征为皮肤黏膜分批出现斑疹、丘疹、水疱和痂，传染性较强［2］，针对免疫力低下的婴幼儿，严重时可引发呼吸道、神经系统损伤等并发症［3］。大量流行病学研究表明，水痘的患病与人群特征、气象因素等密切相关，且呈现高度的空间异质性［4-5］。因此，研究华中区域气候因素对水痘发病流行的影响，可以有效评价该区域水痘流行特征，从而为预防控制水痘流行提供科学依据。

1 资料与方法

1.1 资料来源 资料来源于湖北省孝感市疾病预防控制中心2008年以来水痘疫情报告卡，收集2009年1月—2012年12月患水痘的6 324例患者的临床资料。同期气象数据包括气温、相对湿度、降水量及日照时数，由湖北省孝感市气象局提供。

1.2 研究地点气象资料 湖北省孝感市地处中纬度地带，属亚热带季风气候。境内四季分明，冬季盛行偏北风，夏季盛行偏南风，严寒酷暑时间短，春、初夏、秋温和气候时间长；雨量充沛，年平均降水量为1 133.800 mm；光热充足，年平均日照时数为1 996.700 h，年平均气温为16.000 ℃，≥10.000 ℃的积温为5 091.900 ℃。研究期间孝感市气象资料见表1。

1.3 方 法 将 患 者 按 照 年 龄 分 为 0～＜1、1～＜4、4～＜7、7～＜10、10～＜15、15～＜20、20～ 岁 7 个年龄段，观察水痘的高发年龄，并分析各年龄段的性别分布情况。统计2009—2012年不同月份水痘发病例数的动态变化，比较不同季节水痘发病例数。分析2009—2012年气象因素与水痘发病情况的相关性，以及不同季节气象因素与水痘发病情况相关性。

1.4 统计学方法 采用Levene's test 进行方差齐性检验，季节性差异采用Welch's t-test检验，多重比较采用Dunnet's T3进行季节间比较。水痘患病率与单一气象因素的线性回归模型采用SPSS 20.0（IBM）进行。选择SPSS回归分析模式中Curve Estimation模块，勾选线性回归、二次曲线、三次曲线、指数函数、幂函数、生长曲线和复合函数这7种模型进行气象因子与水痘患病例数之间关系的曲线估计，选择最大相关系数（R2）模型评估气象因素与水痘发病率的相关性。

2 结果

2.1 2009—2012年不同年龄段水痘患者性别分布情况 本研究所划分的7个年龄段中，各年龄段均有水痘病例出现，以4～＜7岁人群为主，＜15岁的人群占总人群的比例高达85.29%。在6 324例水痘患者中，男4 017例（63.52%），女2 307例（36.48%），男女比例为1.74:1，其中0～＜1岁男女比例高达2.43:1（详见表2）。

2.2 2009—2012年不同月份水痘发病例数的动态变化2009—2012年水痘年平均发病率为32.838/10万，无死亡病例。2009—2012年的48个月中，共发水痘6 324例，每月最低发病24例，最高发病297例。2009年水痘发病峰值出现在6月和12月，且维持到了2010年的1月，自2010年开始，其余年份水痘发病峰值均出现在6月和12月，呈现典型的双锋曲线（见图1）。发病高峰出现在夏季和冬季，分别占报告病例总数的34.69%（2 194/6 324）和32.84%（2 077/6 324）。

2.3 2009—2012年不同季节水痘发病例数比较 不同季节水痘发病例数比较，差异有统计学意义（P＜0.001）；夏季水痘发病例数高于春季和秋季，差异有统计学意义（P＜0.01）；冬季水痘发病例数高于春季和秋季，差异有统计学意义（P＜0.05，见表3）。

2.4 2009—2012年气象因素与水痘发病情况的相关性分析水痘发病例数与气温、相对湿度、降水量及日照时数均不存在相关性（P＜0.05，见图2）。

表1 2009—2012年湖北省孝感市气象资料汇总Table 1 Summary of meteorological data from 2009 to 2012 from Xiaogan City, Hubei Province

表2 2009—2012 年孝感市水痘病例年龄和性别分布Table 2 Number of cases of chickenpox and gender distribution in Xiaogan from 2009 to 2012

表3 2009—2012年孝感市不同季节水痘发病例数比较（x ±s，例/月）Table 3 Comparison of the incidence among different seasons in Xiaogan

2.5 2009—2012年不同季节气象因素与水痘发病情况的相关性分析 数据分析表明，气象因素存在显著共线性，多元线性回归模型未能通过F检验（P=0.589）。因此，为了深入探讨季节性差异的成因，拟通过建立单一气象因子的回归模型来评估气象因素对患病情况的影响。复合模型结果显示，春季与冬季气温、日照时数均与水痘发病率呈显著相关（P＜0.05），夏季各气象因子均与水痘发病率呈显著相关（P＜0.05），而秋季仅有日照时数与水痘发病率呈显著相关（P＜0.05，见表4）。

3 讨论

通过对2009—2012年湖北省孝感市水痘患病情况的分析，发现湖北省水痘发患病数呈现典型的双峰曲线动态变化。结合季节差异分析，夏季与冬季为水痘的患病高发期，分别占报告病例总数的34.69%和32.84%。已有的研究表明，东南亚等热带区域的相对较冷时期水痘发病率较高［6］，而温带区域特别是发达国家则夏季水痘发病率较高［7］，而本研究表明，亚热带区域可能呈现与温带及热带区域不同的流行模式。

图1 2009—2012年孝感市不同月份水痘发病例数的动态变化Figure 1 Dynamic change inincidence of chickenpox from 2009 to 2012 in Xiaogan

图2 2009—2012年孝感市气象因素与水痘发病情况的拟合曲线Figure 2 Linear regression between the case number and the meteorological factors from 2009 to 2012 in Xiaogan

表4 2009—2012年孝感市不同季节气象因素与水痘发病情况的相关性分析Table 4 Correlation analysis between chickenpox incidence and meteorological factors in different seasons in Xiaogan

众多研究指出，气候因素可能是影响VZR扩散及水痘感染皮肤并发症（VSC）的关键因素［8］。本研究通过对不同季节水痘发病数的多重比较分析，发现湖北省水痘发病率存在显著的季节性差异，表明气象因素可能对水痘发例数造成影响，但采用总体数据分析表明，水痘发病例数与各气象因子不存在显著相关性。但通过分析季节对水痘发病例数与单一气象因子进行线性拟合，可以发现春季、夏季及冬季气温是影响水痘发病例数的关键因子之一，这与已有的研究结果相一致［8-9］。但相对湿度与水痘发病例数仅在夏季表现出了显著的相关性，这也说明，华中区域表现出与沿海区域不同的发病特点［10］。值得提出的是，日照时数在四季中均与水痘发病例数呈现显著相关，通过复合方程及指数方程等形式较好地阐释了水痘患病时间的方差变动。

综上所述，水痘相对其他烈性传染病来讲，其发病比较温和且有很高的自愈率［11］。但由于其潜在空气传播途径［12］，根据不同的人群特征、环境因素及其他因素，存在突然大规模流行的可能［1］。同时，由于个人免疫体系的差异及其他并发症等情况，水痘仍然有一定的致死可能［7］。因此，深入研究其染病机制及环境要素之间的相互关系，仍然是医疗工作者及卫生部门努力的方向。本研究分析了气象因素中的4个因素对水痘发病的影响，但就人群特征及个体差异尚未展开研究，尚不能排除其他因素对水痘发病率的影响，因此本文研究结论存在一定的局限性。同时，需要注意各气象因素之间高度的多重共线性特征，因此，本研究建立结构性的非线性模型可能可为相关部门制定水痘疫情流行的防治措施提供更为科学的依据。

作者贡献：罗玲进行文章的构思与设计，数据收集；黎明进行研究的实施与可行性分析，数据整理，统计学处理，论文修订；罗玲、黎明进行结果分析与解释，撰写论文；黎明负责文章的质量控制及审校，对文章整体负责，监督管理。

本文无利益冲突。

［1］LEE B R，FEAVER S L，MILLER C A，et al. An elementary school outbreak of varicella attributed to vaccine failure: policy implications［J］.J Infect Dis，2004，190（3）：477-483.

［2］MARCHETTO S，DE BENEDICTIS F M，DE MARTINO M，et al.Epidemiology of hospital admissions for chickenpox in children：an Italian multicentre study in the pre-vaccine era［J］.Acta Paediatrica，2007，96（10）：1490-1493.

［3］LOSURDO G，BERTOLUZZO L，CANALE F，et al.Varicella and its complications as cause of hospitalization［J］.Infez Med，2005，13（4）：229-234.

［4］FAIRLEY C K，MILLER E.Varicella-zoster virus epidemiology：a changing scene?［J］.J Infect Dis，1996，174：S314-319.

［5］YAWN B，YAWN R，LYDICK E.Community impact of childhood varicella infections［J］.J Pediatr，1997，130（5）：759-765.

［6］LEE B W.Review of varicella zoster seroepidemiology in India and Southeast Asia［J］.Trop Med Int Health，1998，3（11）：886-890.

［7］BOZZOLA E，BOZZOLA M，KRZYSZTOFIAK A，et al.Varicella skin complications in childhood：a case series and a systematic review of the literature［J］.Int J Mol Sci，2016，17（5）：E688.DOI：10.3390/ijms17050688.

［8］KOKAZE A，YOSHIDA M，SEKINE Y，et al.The magnitude of variation in temperature within a year has an effect on the seasonal variations of chickenpox incidence in Japan［J］.Epidemiol Infect，2001，126（2）：269-277.

［9］ALMUNEEF M，MEMISH Z A，BALKHY H H，et al.Chickenpox complications in Saudi Arabia：is it time for routine varicella vaccination? ［J］.Int J Infect Dis，2006，10（2）：156-161.

［10］CHAN J Y，TIAN L，KWAN Y，et al.Hospitalizations for varicella in children and adolescents in a referral hospital in Hong Kong，2004 to 2008：a time series study［J］.BMC Public Health，2011，11：366.DOI：10.1186/1471-2458-11-366.

［11］BOELLE P Y，HANSLIK T.Varicella in non-immune persons：incidence，hospitalization and mortality rates［J］.Epidemiol Infect，2002，129（3）：599-606.

［12］SUZUKI K，YOSHIKAWA T，IHIRA M，et al.Spread of varicellazoster virus DNA to the environment from varicella patients who were treated with oral acyclovir［J］.Pediatr Int，2003，45（4）：458-460.