廣州市2013—2015年氣象因子與手足口病發(fā)病關(guān)聯(lián)的時(shí)間序列研究

劉靜雯 謝華萍

[摘要] 目的 探究氣象因子與廣州市手足口病發(fā)病的關(guān)聯(lián)性。方法 收集2013—2015年廣州市手足口病發(fā)病例數(shù)及同期的氣溫、氣壓、相對(duì)濕度、風(fēng)速、降雨量等氣象數(shù)據(jù)，使用分布滯后非線性模型分析氣象因素與手足口病發(fā)病的關(guān)系及其滯后效應(yīng)。結(jié)果 廣州市2013—2015年共有手足口病病例5 919例，發(fā)病高峰期在夏秋季節(jié)。各氣象因子與手足口病發(fā)病均呈非線性關(guān)系。氣溫、降雨量、濕度均與手足口病的發(fā)病呈正相關(guān)，而平均氣壓與其呈負(fù)相關(guān)。最高氣溫與不同性別的手足口病發(fā)病均呈非線性關(guān)系。降雨量對(duì)手足口病發(fā)病有影響。結(jié)論 氣象因素可能是影響廣州市手足口病發(fā)病的一個(gè)重要因素，在以后的手足口病防控中應(yīng)重視這一發(fā)病指征。

[關(guān)鍵詞] 手足口病;氣象因子;分布滯后非線性模型

[中圖分類號(hào)] R725.1 [文獻(xiàn)標(biāo)識(shí)碼] A [文章編號(hào)] 1672-5654（2019）01（a）-0178-05

手足口?。╤and，foot and mouth disease，HFMD） 是一種常見(jiàn)的急性傳染病，多發(fā)于5歲（特別是3歲）及以下兒童[1]。該病可由多種腸道病毒引起，而柯薩奇病毒A16型（Cox A16）和腸道病毒71型（EV71）是其常見(jiàn)的兩種病原，它的潛伏期為2～10 d，病程一般為7～10 d[1]。手足口病的主要癥狀是以發(fā)熱和手、足、口腔等部位出現(xiàn)皮疹或皰疹，少數(shù)患兒可發(fā)生多種嚴(yán)重并發(fā)癥，個(gè)別重癥病情進(jìn)展很快甚至死亡[2]。2013—2015年，廣州市手足口病高發(fā)期從3月份開(kāi)始，至7月份有所緩解。據(jù)文獻(xiàn)報(bào)道，手足口病發(fā)病與氣象因素之間存在一定的聯(lián)系[3-4]。因此，該研究根據(jù)廣州市2013—2015年氣象數(shù)據(jù)和手足口病病例資料，利用分布滯后非線性模型（DLNM），研究氣象因素對(duì)手足口病發(fā)病的關(guān)系，為今后對(duì)該病防控提供科學(xué)的參考依據(jù)。

1? 資料與方法

1.1? 一般資料

從中國(guó)疾控傳染病報(bào)告信息管理系統(tǒng)中取得2013年1月—2015年12月廣州市手足口病病例信息（包括臨床診斷和實(shí)驗(yàn)室確診病例）。并從廣州市氣象局獲取包括氣溫、氣壓、相對(duì)濕度、風(fēng)速等的同期數(shù)據(jù)。

1.2? 標(biāo)本采集

2013年1月—2015年12月廣州市哨點(diǎn)監(jiān)測(cè)醫(yī)院的疑似手足口病病例及疑似暴發(fā)疫情的標(biāo)本，共有5 919例，采集疑似及確診病例的糞便、肛拭子或皰疹液。各哨點(diǎn)醫(yī)院或區(qū)疾控中心采集標(biāo)本后凍存于-20℃，并將標(biāo)本及臨床資料送往廣州市疾病預(yù)防控制中心實(shí)驗(yàn)室檢測(cè)[5]。

1.3? 統(tǒng)計(jì)方法

用Spearman等級(jí)相關(guān)分析廣州市各氣象因子與手足口病日發(fā)病人數(shù)的相關(guān)性以及各氣象因子間的相關(guān)程度。利用分布滯后非線性模型對(duì)氣象因子與手足口病日發(fā)病人數(shù)的關(guān)聯(lián)進(jìn)行擬合，基本模型如下[6]：

上式中Yt表示第t天的手足口病發(fā)病人數(shù);Xi是除M以外的其他與手足口病發(fā)病人數(shù)相關(guān)的解釋變量;DOWt是第t天為星期幾的啞變量;Time表示季節(jié)性和長(zhǎng)期趨勢(shì)。根據(jù)赤池信息量準(zhǔn)則（AIC），最終把日最高氣溫納入模型中計(jì)算。參考值為廣州市各氣象因子的中位數(shù)，并計(jì)算在整個(gè)滯后時(shí)間內(nèi)各氣象因子對(duì)手足口病發(fā)病的相對(duì)危險(xiǎn)度（RR，總體效應(yīng)），本研究采用R 3.3.0中的DLNM軟件進(jìn)行統(tǒng)計(jì)分析[6]。

2? 結(jié)果

2.1? 手足口病基本情況

2013—2015年廣州市共報(bào)告手足口病病例5 919例，3年報(bào)告病例數(shù)分別為1 632、2 289和1 998例;其中柯薩奇A16型（Cox A16）1 483例（25.05%，1 483/5 919），腸道病毒71型（EV71）1 073例（18.13%，1 073/5 919）。發(fā)病高峰集中夏秋季節(jié)，5～7月份病例最多。日平均、最高和日最低氣溫中位數(shù)分為23.50℃、28.20℃、20.10℃，日降雨量中位數(shù)為0 mm，日均風(fēng)速中位數(shù)為3.32 m/s。見(jiàn)圖1和表1。

2.2? 氣象因子與手足口病發(fā)病間的關(guān)系

2.2.1 各氣象因子與每日手足口病發(fā)病的相關(guān)性? 由Spearman等級(jí)相關(guān)分析得出，（最高、平均、最低）氣溫、降雨量、平均相對(duì)濕度與手足口病發(fā)病呈正相關(guān)（r>0），平均氣壓與其發(fā)病呈負(fù)相關(guān)（r<0）。Cox A16、EV71、男性病例均與各氣象因子相關(guān) （P<0.05），其他型別的腸道病毒（除風(fēng)速和相對(duì)濕度外）、女性病例（除風(fēng)速外）均與氣象因子相關(guān)（P<0.05）。見(jiàn)表2。

2.2.2 各氣象因子在不同滯后時(shí)間的手足口病發(fā)病風(fēng)險(xiǎn)及總體效應(yīng)? 通過(guò)分布滯后非線性模型DLNM擬合各氣象因子在不同滯后時(shí)間與手足口病發(fā)病的情況，氣壓與最高氣溫有強(qiáng)相關(guān)關(guān)系（r=0.77），根據(jù)AIC排除氣壓。結(jié)果顯示，各個(gè)氣象因子都與手足口病的發(fā)病呈非線性關(guān)系。當(dāng)最高氣溫為37℃，滯后為25 d時(shí)，RR值最大為1.326（95%CI：1.013～1.735）;當(dāng)平均相對(duì)濕度為35%，滯后為30 d時(shí)，RR值最大為1.216（95%CI：1.022～1.447），均對(duì)手足口發(fā)病的影響有統(tǒng)計(jì)學(xué)意義。當(dāng)降雨量135 mm，滯后21 d的累積效應(yīng)值為5.996（95%CI：1.045～34.416），對(duì)手足口病發(fā)病的影響有統(tǒng)計(jì)學(xué)意義。最高氣溫37℃，平均相對(duì)濕度35%的累積效應(yīng)，均對(duì)手足口病發(fā)病的影響無(wú)統(tǒng)計(jì)學(xué)意義。見(jiàn)圖2。

2.2.3 氣溫在不同滯后時(shí)間的各型別手足口病發(fā)病風(fēng)險(xiǎn)及總體效應(yīng)? 通過(guò)DLNM擬合氣溫（平均、最高、最低）在不同滯后時(shí)間與各型別手足口病發(fā)病的情況，根據(jù)AIC，將日最高氣溫納入模型中計(jì)算。結(jié)果顯示，最高氣溫與各型別手足口病發(fā)病均呈非線性關(guān)系。當(dāng)最高氣溫為37℃，滯后為25 d時(shí)，Cox A16型的RR值最大為1.472（95%CI：1.121～1.933）;當(dāng)最高氣溫為37℃，滯后為當(dāng)天時(shí)，EV71型的RR值最大為1.691（95%CI：1.198～2.387），對(duì)兩種型別的手足口病發(fā)病的影響均有統(tǒng)計(jì)學(xué)意義。見(jiàn)圖3。

2.2.4 氣溫在不同滯后時(shí)間與不同性別的手足口病發(fā)病風(fēng)險(xiǎn)及總體效應(yīng)? 氣溫（平均、最高、最低）通過(guò)分布滯后非線性模型將不同滯后時(shí)間與不同性別的手足口病發(fā)病情況進(jìn)行擬合，根據(jù)AIC，將日最高氣溫納入模型中計(jì)算。結(jié)果顯示，最高氣溫與不同性別的手足口病發(fā)病均呈非線性關(guān)系。當(dāng)最高氣溫為37℃，滯后為25 d時(shí)，女性病例的RR值最大為1.474（95%CI：1.051～2.067），只對(duì)手足口病發(fā)病為女性的影響有統(tǒng)計(jì)學(xué)意義。見(jiàn)圖4。

3? 討論

手足口病全年均可發(fā)生，多發(fā)生于夏秋季，本病以糞口、呼吸道和密切接觸等途徑傳播，該病發(fā)生不僅與年齡、機(jī)體免疫力等因素有關(guān)，而且與氣象因素存在復(fù)雜的關(guān)系[7]。廣州市全年平均氣溫為20～22℃，以溫暖多雨、溫差小、光熱充足、夏季長(zhǎng)為特征。廣州市手足口病的高發(fā)可能與其獨(dú)特的氣候特點(diǎn)有關(guān)。

氣象因素對(duì)手足口病的發(fā)生和流行起著重要作用。該研究對(duì)廣州市2013—2015年手足口病的監(jiān)測(cè)結(jié)果顯示，病例主要集中在每年的4～7月，而9～10月有一個(gè)小高峰，提示氣溫升高可能是手足口病高發(fā)的原因之一[8]。據(jù)該研究顯示，氣溫、氣濕、氣壓和降雨量等的變化對(duì)廣州手足口病的發(fā)病有影響，而氣溫、降雨量、平均相對(duì)濕度與手足口病發(fā)病之間呈正相關(guān)，對(duì)手足口病總體的影響體現(xiàn)在氣溫漸高，手足口病的發(fā)病增多;降雨量越大，發(fā)病越多;平均氣壓與手足口病發(fā)病率呈負(fù)相關(guān)，這與深圳市[9]、日本[10]等的研究結(jié)果類似。近年來(lái)該市手足口病呈現(xiàn)上升趨勢(shì)，也有夏秋季高發(fā)的特點(diǎn)，雖然目前尚不清楚氣象因素作用于手足口病的機(jī)理，但是相關(guān)等研究[9-10]也發(fā)現(xiàn)，不同的溫度和濕度等氣象條件對(duì)手足口病發(fā)病影響較大，高溫高濕的天氣更有利該病病原體在環(huán)境中的生存和繁殖，低氣壓環(huán)境容易引起人體機(jī)體缺氧而使抵抗力下降。同時(shí)不同氣候會(huì)影響兒童的日?；顒?dòng)、接觸的環(huán)境，進(jìn)而影響手足口病傳播的幾率[9]。

氣象因素的健康效應(yīng)有一定的持續(xù)性和滯后性，人群健康指標(biāo)（如發(fā)病率） 可能與當(dāng)天、昨天或者十多天前的暴露水平有關(guān)[11]。并且手足口病病毒常在體內(nèi)潛伏3～7 d才發(fā)病，采用DLNM可以檢測(cè)氣象因素與其發(fā)病的相對(duì)的延遲效應(yīng)。該研究通過(guò)用DLNM分析本市的日均氣溫、日均相對(duì)濕度、日均風(fēng)速、日均氣壓和降雨量分別與手足口病的關(guān)聯(lián)性，結(jié)果顯示氣象因素與手足口病發(fā)病呈非線性關(guān)系，與潘潔[12]的研究結(jié)果一致。其中，氣溫在同樣的滯后時(shí)間對(duì)不同型別的影響體現(xiàn)在Cox A16和EV71型都受到氣溫的影響，而其他型別的腸道病毒受氣溫影響不明顯。Cox A16與EV71相比，EV71滯后時(shí)間較短，受氣溫的影響更大。并且氣溫對(duì)不同性別手足口病發(fā)病均有影響，對(duì)女性發(fā)病的影響更大。提示夏季氣溫較高、相對(duì)濕度較低、風(fēng)速較大時(shí)，手足口病的發(fā)病人數(shù)會(huì)增加，這與徐曉楠[13]的研究結(jié)果類似。由于廣州獨(dú)特的天氣特點(diǎn)，無(wú)明顯的冬季，夏季高溫天氣長(zhǎng)，在高溫的條件下，兒童受到家長(zhǎng)的管制可能會(huì)減少在室外活動(dòng)，而增加室內(nèi)活動(dòng)（如室內(nèi)游樂(lè)場(chǎng)，培訓(xùn)室等場(chǎng)所），進(jìn)而增加與同伴接觸，接觸傳染源的機(jī)會(huì)更大，而且手足口病常通過(guò)被病毒污染的手、生活用品等或與患者密切接觸或通過(guò)空氣（飛沫）傳播而感染[14]，于是更容易因增加接觸患者而感染手足口病的機(jī)會(huì)。

手足口病已成為廣州重點(diǎn)防治的腸道傳染病，EV71和Cox A16型是手足口病的常見(jiàn)病原，而EV71更是引起手足口病重癥及死亡的主要病原，并且近年來(lái)其他腸道病毒中的Cox A6、Cox A10型也有一定比例存在，在手足口病高峰期要及時(shí)進(jìn)行各種病原檢測(cè)尤為重要，特別是發(fā)現(xiàn)由EV71引起的病例要及時(shí)救治，減少重癥的發(fā)生[5，15]。廣州市手足口病發(fā)病與氣溫、濕度、降雨量等氣象因子有關(guān)，有可能是本市兒童患手足口病的危險(xiǎn)因素之一，因此可以將其作為預(yù)測(cè)本市手足口病發(fā)病的指標(biāo)。目前，市面上只有預(yù)防EV71型手足口病疫苗，但早發(fā)現(xiàn)、早治療才是預(yù)防手足口病的重要措施。相關(guān)部門(mén)應(yīng)加強(qiáng)對(duì)公共場(chǎng)所設(shè)施的清潔和消毒的監(jiān)管;家庭應(yīng)做到勤洗手、勤開(kāi)窗、曬衣被、喝開(kāi)水、食熟食;如果有兒童感染，應(yīng)在家休息隔離，嚴(yán)防交叉感染。

[參考文獻(xiàn)]

[1]? 中華人民共和國(guó)衛(wèi)生部.手足口病預(yù)防控制指南[M].北京：中華人民共和國(guó)衛(wèi)生部，2010.

[2]? Xing W，Liao Q，Viboud C，et al.Hand，foot and mouth disease in China，2008-2012：an epidemiological study[J].Lancet Infect Dis，2014，14（4）：308-318.

[3]? 劉峰，劉鳳仁，李剛.深圳市龍崗區(qū)2011-2014年手足口病特征及與氣象因素相關(guān)性[J].中國(guó)熱帶醫(yī)學(xué)，2016，16（7）：688-691.

[4]? Huang Y，Deng T，Yu S，et al. Effect of meteorological variables on the incidence of hand，foot，and mouth disease in children：a time-series analysis in Guangzhou，China[J].BMC Infect Dis，2013，13：134.

[5]? 劉靜雯，耿進(jìn)妹，陳純，等.2011-2015年廣州市手足口病流行病學(xué)特征及病原學(xué)監(jiān)測(cè)分析[J].現(xiàn)代預(yù)防醫(yī)學(xué)，2017，44（5）：772-777.

[6]? 董航，李曉寧，劉華章，等.2012—2013年廣州市心血管病逐日死亡數(shù)與氣象因子的時(shí)間序列分析[J].熱帶病與寄生蟲(chóng)學(xué)，2016，14（3）：155-159.

[7]? 朱建明，李瀾，莫平華，等.基于氣象因素的上海市金山區(qū)手足口病多元線性回歸預(yù)測(cè)模型[J].實(shí)用預(yù)防醫(yī)學(xué)，2016， 23（1）：115-116.

[8]? 呂世峰，王惠君，劉慧玲，等.山東日照市2009—2012年手足口病與氣象因素的時(shí)間序列分析[J].華南預(yù)防醫(yī)學(xué)，2014， 40（4）：312-317.

[9]? Lin HL，Zou H，Nie J，et al.Short term effects of El Nino-Southern Oscillation on hand，foot，and mouth disease in Shenzhen，China[J]. PloS One，2013，8（7）：65585.

[10]? Onozuka D，Hashizume M.The influence of temperature and humidity on the incidence of hand，foot，and mouth disease in Japan[J].Sci Total Environ，2011，410-411;119-125.

[11]? 楊瓊英.廣州市流行性腮腺炎流行病學(xué)特征及與氣象因子關(guān)聯(lián)性研究[D].廣州：南方醫(yī)科大學(xué)，2015.

[12]? 潘潔，楊軍，李萌萌，等.分布滯后非線性模型在氣象因素對(duì)手足口病發(fā)病影響研究中的應(yīng)用[J].環(huán)境與健康雜志，2016，33（5）：377-380.

[13]? 徐曉楠，王金玉，李盛.嘉峪關(guān)市氣象因素與手足口病滯后效應(yīng)研究[J].中國(guó)公共衛(wèi)生管理，2016，32（3）：312-314.

[14]? 陳宏標(biāo)，楊偉康，王金明.一起托幼機(jī)構(gòu)手足口病暴發(fā)的病原學(xué)診斷及其傳播途徑分析[J].現(xiàn)代預(yù)防醫(yī)學(xué)，2015， 42（13）：2323-2325.

[15]? 王春榮，關(guān)恒云，楊國(guó)樑，等.濟(jì)南地區(qū)2009—2013年手足口病分子流行病學(xué)特征[J].中國(guó)公共衛(wèi)生，2016，32（1）：53-56.

（收稿日期：2018-09-26）