劉靜雯 謝華萍

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

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

[中圖分類號] R725.1 [文獻標識碼] A [文章編號] 1672-5654（2019）01（a）-0178-05

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

1? 資料與方法

1.1? 一般資料

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

1.2? 標本采集

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

1.3? 統(tǒng)計方法

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

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

2? 結(jié)果

2.1? 手足口病基本情況

2013—2015年廣州市共報告手足口病病例5 919例，3年報告病例數(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。見圖1和表1。

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

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

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

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

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

3? 討論

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

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

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

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

[參考文獻]

[1]? 中華人民共和國衛(wèi)生部.手足口病預(yù)防控制指南[M].北京：中華人民共和國衛(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].中國熱帶醫(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]? 劉靜雯，耿進妹，陳純，等.2011-2015年廣州市手足口病流行病學(xué)特征及病原學(xué)監(jiān)測分析[J].現(xiàn)代預(yù)防醫(yī)學(xué)，2017，44（5）：772-777.

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

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

[8]? 呂世峰，王惠君，劉慧玲，等.山東日照市2009—2012年手足口病與氣象因素的時間序列分析[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]? 潘潔，楊軍，李萌萌，等.分布滯后非線性模型在氣象因素對手足口病發(fā)病影響研究中的應(yīng)用[J].環(huán)境與健康雜志，2016，33（5）：377-380.

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

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

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

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