溫州地區氣象因素與川崎病發病的相關性

周愛華，岳亞延，金增游，潘璐璐，張營營，王丹，錢燕，張園海，吳蓉洲，褚茂平

（1.溫州醫科大學附屬第一醫院 兒科，浙江 溫州 325000；2.江蘇大學附屬醫院 神經內科，江蘇鎮江 212000；3.溫州醫科大學附屬第二醫院育英兒童醫院 兒童心血管科，浙江 溫州 325027）

?論 著?

溫州地區氣象因素與川崎病發病的相關性

周愛華1，岳亞延2，金增游1，潘璐璐3，張營營1，王丹1，錢燕1，張園海3，吳蓉洲3，褚茂平3

（1.溫州醫科大學附屬第一醫院 兒科，浙江 溫州 325000；2.江蘇大學附屬醫院 神經內科，江蘇鎮江 212000；3.溫州醫科大學附屬第二醫院育英兒童醫院 兒童心血管科，浙江 溫州 325027）

目的：探討氣象因素與川崎?。↘D）發病的相關性。方法：收集2005年1月至2015年12月溫州地區住院的KD患兒資料及該地區氣象資料，對每月KD發病人數與同時期該地區氣象因素，包括降雨量、風速、日照時長、氣溫及氣壓進行相關性分析。結果：2005年1月至2015年12月溫州地區約1 956例KD住院患兒納入本研究，KD發病率以春夏季為高，冬季發病率較低。降雨量與KD發病率成正相關（r=0.217，P＜0.05），氣壓、風速與KD發病率成負相關（r=-0.209、-0.652，P＜0.05），氣溫、日照時數與KD發病率相關性不顯著（P＞0.05）。結論：KD的發病與氣象因素具有一定相關性。

川崎??；風速；降雨量

川崎?。↘awasaki disease，KD）是一種多發生于5歲以內兒童，以全身中小血管炎為主要病理特征的急性發熱出疹性疾病，具有一定的自限性[1-2]。自1967年川崎富作先生首次報道[3]以來，至今已有40余年，但其病因及發病機制仍不明確[4-9]。有研究表明，氣象因素對KD發病具有一定影響[10-12]，且KD的發病具有季節性、地域性、流行性和種族差異[13]。由于不同季節、不同地域、不同時間對應不同的氣象條件，據此推測KD發病可能與氣象因素存在一定關系，有必要對其展開研究。目前國內外這方面的研究比較少見。本研究通過對溫州市2005 年1月至2015年12月每月降雨量、風速、日照時數、氣溫和氣壓與KD月發病率的相關性分析，篩選出與KD發病相關的氣象因素，并探討氣象因素對KD發病產生影響的可能原因。

1 資料和方法

1.1 一般資料 收集2005年1月至2015年12月溫州市所有開展醫療衛生服務的二級甲等及以上醫院的全部KD患兒臨床資料，內容包括：姓名、性別、年齡、發病日期、就診日期、出院診斷、臨床表現等。病例納入采用2007年中國第八次KD專題討論會紀要的KD診斷標準[14]。典型KD（CKD）為發熱5 d以上，伴下列5項臨床表現中的4項，排除其他疾病后，即可診斷為KD：①眼結膜充血，非化膿性；②唇充血破裂，口腔黏膜彌漫充血，舌乳頭突起、充血，呈草莓舌；③多形性紅斑；④四肢變化：急性期掌跖紅斑，手足硬性水腫，恢復期指（趾）端膜狀脫皮；⑤頸部淋巴結腫大。如具備除發熱以外3項臨床表現，并超聲心動圖提示有冠狀動脈損害者，亦可確診。除發熱5 d或以上外，具備至少2項KD主要臨床表現，并具備炎癥指標明顯升高，除外其他疾病時，可考慮不完全KD（IKD），若超聲心動圖提示冠狀動脈損害，則可明確診斷。排除標準：①非溫州戶口居民；②溫州戶口居民但不在溫州生活者；③戶口或居住地信息不詳者；④資料不全且診斷不明確者；⑤KD冠狀動脈后遺癥期，因其他疾病入院者。

1.2 人口資料 2005年1月至2015年12月溫州市常住居民兒童人口數據由溫州市統計局提供，根據其資料進行整理。KD患兒主要為＜5歲嬰幼兒[14-15]，采用＜5歲兒童發病率進行比較分析。根據溫州市每月＜5歲KD患兒發病數與該年份同年齡段的常住居民總數，計算出＜5歲兒童KD的每月發病率（/10萬）。

1.3 地理和氣象資料 溫州市位于中國東南部，介于東經119°37′-121°18′、北緯27°03′-28°36′之間，屬于典型的中亞熱帶季風氣候。本研究中季節劃分采用中國氣象季節劃分方式，3-5月為春季，6-8月為夏季，9-11月為秋季，12-2月為冬季。本研究所采用的2005年1月至2015年12月溫州市氣象數據由浙江省氣象局提供，是溫州氣象站（站臺號58659）以月為單位的常規地面氣象觀測數據。所有數據均符合世界氣象組織（World Meteorological Organization，WMO）全球觀測系統規范和中國氣象局系統觀測規范，可良好地代表溫州市的氣象環境。具體內容包括：每月降雨量、風速、日照時數、氣溫及氣壓。其中風速、氣溫、氣壓采用月平均值進行分析，降雨量和日照時數則采用每月總量進行分析。

1.4 數據質量控制 所有病例的收集、篩選和統計均由兒童心血管專業人員進行，篩選并剔除資料不全或診斷不詳者。采用雙輸入方式錄入，并經過雙份錄入的一致性比較，以確保數據的準確性及可靠性。數據輸入及管理采用Epidata 3.1。

1.5 統計學處理方法 采用SPSS19.0統計學軟件進行數據分析。采用spearman秩相關性檢驗將溫州市＜5歲兒童KD月發病率與同時段氣象因素（氣溫、氣壓、風速、降雨量及日照時數）進行相關性分析。采用多元線性回歸方法分析各氣象因素與KD發病率的依存關系。計量資料用±s表示。P＜0.05為差異有統計學意義。

2 結果

2.1 KD及氣象因素的一般情況 2005年1月至2015 年12月溫州地區KD患兒住院例數共1 956例，符合納入標準者1 951例，其余5例資料不全予以剔除。＜5歲KD患兒1 764例（占90.42%），其中CKD 1 473例（占83.5%），IKD 291例（占16.5%）。＜5歲患兒中男1 164例（占65.99%），女600例（占34.01%）。＜5歲KD患兒發病率呈逐年升高趨勢，2005至2015年平均發病率為40.80/10萬，男性患兒平均發病率為47.26/10萬，女性患兒平均發病率為32.44/10萬，男女發病率比例約1.46∶1，此期間未出現KD暴發流行現象，見圖1，與溫州既往流行病學調查結果一致[16]。根據我國季節劃分方式，溫州市KD發病率以春夏季為高，冬季發病率較低，見圖2。

圖1 溫州市2005至2015年＜5歲兒童KD發病率

圖2 2005至2015年溫州市＜5歲兒童KD月平均發病率

KD的發病具有季節性，而不同的季節具有不同的氣象特點。為明確不同的季節是否存在氣象差異，我們對溫州市2005年1月至2015年12月氣象因素（降雨量、風速、日照時數、氣溫、氣壓）進行了比較分析。結果提示除風速外，不同季節間的氣象因素總體差異有統計學意義（P＜0.01）。進一步采用Bonferroni檢驗對各季節的氣象因素進行多重兩兩比較，其中降雨量、日照時數及氣溫以夏季最高，冬季最低（P＜0.01），氣壓則是冬季最高，夏季最低（P＜0.01），而風速在四季中差異無統計學意義（P ＞0.05）。見表1。

表1 2005至2015年溫州市各季節氣象因素情況（±s）

2.2 氣象因素與KD的相關性分析 由于氣象因素及＜5歲KD月發病率均呈非正態分布，故采用spearman秩相關檢驗對兩者進行分析，結果顯示5項氣象因素中，風速、降雨量、氣壓與KD的發病率具有顯著相關性，其中風速（r=-0.652，P＜0.01），氣壓（r=-0.209，P＜0.05）與KD的發病率呈負相關，降雨量與KD的發病率呈正相關（r=0.217，P＜0.05），氣溫、日照時數與KD發病率的相關性無統計學意義（P＞0.05）。通過各氣象因素與不同性別KD發病率相關性分析得出，與男女性別KD發病率相關性有統計學意義的氣象因素均為風速、降雨量和氣壓，與合計發病率保持一致。見表2。

表2 spearman相關分析氣象因素與KD發病率相關性

表3 多元線性回歸分析氣象因素與KD發病率相關性（r）

2.3 擬合預測模型 采用多元線性回歸方法分析各氣象因素與KD發病率的數量依存關系。經分析可知各氣象因素間存在一定相關性（見表3），為避免氣象因素對KD發病的共線性作用對擬合模型的影響，首先對各氣象因素進行共線性診斷，以容忍度＜0.1，方差膨脹因子＞10為多重共線性標準。結果提示各氣象因素間不存在多重共線性（見表4），故均被納入回歸方程。采用逐步選擇法中的前進法進行回歸方程擬合，預測方程為：?=8.726-0.011X1-5.919X2+0.093X3+0.078X4-0.089X5（R2=0.386，P＜0.01）（?：發病率，/10萬；X1：降雨量，mm；X2：風速，m/s；X3：日照時數，h；X4：氣溫，℃；X5：氣壓，hpa）。

表4 氣象因素的共線性分析

2.4 風速與發病率 為了進一步了解各氣象因素如何影響KD的發病，更好地服務于臨床實踐，我們將各氣象因素分別與KD發病率進行曲線擬合。發現風速與發病率的擬合度較好（見圖3），結合表3的spearman相關系數及上述多元線性回歸方程的標準化回歸系數均以風速為大，據此認為氣象因素中風速對KD發病的影響最大。溫州市2005年1月至2015 年12月的月平均風速范圍為[0.4，1.4]m/s。根據圖3A曲線，將風速分為[0.4，0.7]m/s，[0.7，1.2]m/s及[1.2，1.4]m/s，分別與相應月份的發病率進行相關性分析表明，在風速[0.4，0.7]m/s或[1.2，1.4]m/s時，發病率隨風速的增大而增加，但趨勢沒有統計學意義，在風速介于[0.7，1.2]m/s時，發病率與風速呈顯著負相關（r=-0.652，P＜0.01）。不同性別KD發病率與風速關系的與總體基本一致，見圖3。

3 討論

本研究將降雨量、風速、日照時數、氣溫和氣壓5項氣象因素分別與KD的發病率進行相關性分析，顯示風速、降雨量、氣壓與KD的發病率具有相關性（P＜0.01），其中風速、氣壓與KD的發病率呈負相關，降雨量與KD的發病率呈正相關，氣溫及日照時數與KD發病率相關性不顯著。為進一步分析各氣象因素與發病率的數量依存關系，采用多元線性回歸方法進行統計分析。盡管spearman相關性分析中氣溫及日照時數與KD發病率的相關性不顯著，但考慮人們所處的氣象環境是由多種因素組成的統一整體，且在氣象因素的共線性診斷中，氣溫、日照時數與其他氣象因素間不存在多重共線性，故仍被納入回歸方程。現就各氣象因素與KD發病率的相關性研究結果，分別分析如下：

圖3 風速與KD發病率的相關性

3.1 風速 相關系數、多元線性回歸方程的標準化回歸系數及各氣象因素與KD發病率曲線擬合圖形均提示在氣象因素中，風速對＜5歲兒童KD發病的影響最大。本研究結果顯示，KD的發病與風速呈顯著負相關（見表2），尤其在風速為[0.7，1.2]m/s時，發病率隨風速的增大而降低（r=-0.652，P＜0.01）。本研究認為風速較小時KD的發病率高，可能因為風速較小，病原微生物不易隨風擴散，導致單位空間內的病原微生物密度較大，增加KD的致病風險。雖然圖3顯示在風速[0.4，0.7]m/s或[1.2，1.4]m/s時，發病率隨風速的增大而增加，但趨勢沒有統計學意義，原因可能是在這2種風速范圍內的KD樣本量偏少，誤差較大所致。由于我們沒有更低的或更高的月平均風速，故缺少對風速小于0.4 m/s或大于1.4 m/s時與KD發病率的關系分析。

3.2 降雨量 Bronstein等[11]將1994至1998年美國圣迭戈市KD月平均發病率與月降雨量進行相關性分析，顯示KD的發病率與降雨量呈正相關，與本研究結果一致。溫州市冬季降雨量相較其他季節少，KD的發病率也較低（P＜0.05），支持上述相關性結果。有學者[17]報道KD的發病與室內加濕器的應用有關，經常清洗地毯的家庭的兒童和居住在河邊的兒童更加容易患上KD，提示空氣濕度也許為KD發病的影響因素之一[13，18]。據此我們推測，降雨量可能通過增加空氣濕度來影響KD的發病。降雨量多時空氣濕度較大，有利于大部分微生物的滋生和繁殖，增加了感染的風險。

3.3 氣壓 在大自然中，氣壓的明顯變化一般伴隨著降水、大風、低溫或者高溫、高濕等天氣現象[19]，提示氣壓改變受其他氣象因素的影響較大。本研究結果顯示KD的發病率與氣壓成負相關，考慮氣壓較低時，空氣中的粉塵顆粒、化學污染物、過敏原、病原菌等流動性差，不易擴散，滯留于此環境中，易增加過敏原及病原菌感染的風險。因此，低氣壓可能通過增加病原微生物的感染風險來增加KD的發病。

3.4 氣溫 Bronstein等[11]認為KD的發病率與月平均氣溫呈負相關，而本實驗研究結果則提示兩者的相關性不顯著，研究結果存在差異。原因可能有以下幾點：①冬季天氣寒冷時，居民可能存在開空調現象，氣象臺監測氣溫與兒童實際所處空調氣溫差別較大，致使研究結果出現誤差。②在不同氣溫范圍，KD發病率與氣溫的相關性結果可能存在不同。溫州與圣迭戈氣溫范圍不同，導致兩者相關性結果的差異。③可能因圣迭戈和溫州KD病例樣本量偏少，影響統計結果，有待繼續隨訪增大樣本量進行分析。所以，盡管spearman相關性結果提示氣溫與KD發病率相關性不顯著，亦不能完全排除其對KD發病的影響。

3.5 日照時數 日照可以通過太陽輻照影響人體健康，與輻射的強度、被機體吸收的程度和它的生物學效應相關[20]。本研究中，日照時數與KD發病率相關性不顯著，可能由于溫州市嬰幼兒大多數時間居于室內，暴露于日照的時間與溫州真正的日照時數不成正比，致使研究結果出現誤差，所以亦不能完全排除其對KD發病的影響。

3.6 性別差異 KD的發病具有性別差異。為了解氣象因素對不同性別兒童KD的發病的影響，本研究將＜5歲不同性別KD發病率分別與氣象因素進行相關性分析，結果提示氣象因素對KD發病的影響沒有明顯的性別差異。

由于溫州市地域跨度相對較小，缺乏極端氣候，本研究的結論是否適用于極端氣候條件尚有待進一步研究。人類所處的氣象環境是多種氣象因素綜合起來的復雜整體，各氣象因素間相互影響，本研究也只涉及了其中常見的一小部分，對于氣象因素與KD發病的關系，仍需要大家共同努力探索。

[1]NEWBURGER J W, FULTON D R.Kawasaki disease[J].Curr Opin Pediatr, 2004, 16(5): 508-514.

[2]DANIELS L B, TJAJADI M S, WALFORD H H, et al.Prevalence of Kawasaki disease in young adults with suspected myocardial ischemia[J].Circulation, 2012, 125(20): 2447-2453.

[3]KAWASAKI T.Acute febrile mucocutaneous syndrome with lymphoid involvement with specific desquamation of the fingers and toes in children[J].Arerugi, 1967, 16(3): 178-222.

[4]WOOD L E, TULLOH R M.Kawasaki disease in children [J].Heart, 2009, 95(10): 787-792.

[5]NAKAMURA Y, YANAGAWA H, et al.The worldwide epidemiology of Kawasaki disease[J].Prog Pediatr Cardiol, 2004, 19(2): 99-108.

[6]BURNS J C.Kawasaki disease update[J].Indian J Pediatr, 2009, 76(1): 71-76.

[7]CURTIS N.Kawasaki disease[J].BMJ, 1997, 315(7104): 322-323.

[8]KUO H C, YANG K D, CHANG W C, et al.Kawasaki disease: an update on diagnosis and treatment[J].Pediatr Neonatol, 2012, 53(1): 4-11.

[9]GOMARD-MENNESSON E, LANDRON C, DAUPHIN C, et al.Kawasaki disease in adults report of 10 cases[J].Medicine, 2010, 89(3): 149-158.

[10]TSE S M, SILVERMAN E D, MCCRINDLE B W, et al.Early treatment with intravenous immunoglobulin in patients with Kawasaki disease[J].J Pediatr, 2002, 140(4): 450-455.

[11]BRONSTEIN D E, DILLE A N, AUSTIN J P, et al.Relationship of climate, ethnicity and socioeconomic status to Kawasaki disease in San Diego County, 1994 through 1998 [J].Pediatr Infect Dis J, 2000, 19(11): 1087-1091.

[12]RODO X, BALLESTER J, CAYAN D, et al.Association of Kawasaki disease with tropospheric wind patterns[J].Sci Rep, 2011, 1(10): 152.

[13]BURNS J C, HERZOG L, FABRIL O, et al.Seasonality of Kawasaki disease: A global perspective[J].PLoS One, 2013, 8(9): e74529.

[14]《中華兒科雜志》編輯委員會.川崎病專題討論會紀要[J].中華兒科雜志, 2007, 45(11): 826-829.

[15]HARNDN A, TULLOH R, BURGNER D.Kawasaki disease[J].BMJ, 2014, 349: g5336.

[16]張營營, 張園海, 周愛華, 等.2001-2010年溫州市川崎病住院患兒流行病學調查[J].中華風濕病學雜志, 2012, 16(11): 763-766.

[17]TREADWELL T A, MADDOX R A, HOLMAN R C, et al.Investigation of Kawasaki syndrome risk factors in Colorado [J].Pediatr Infect Dis J, 2002, 21(10): 976-978.

[18]RAUCH A M, GLODE M P, WIGGINS J W, et al.Outbreak of Kawasaki syndrome in Denver, Colorado: association with rug and carpet cleaning[J].Pediatrics, 1991, 87(5): 663-669.

[19]畢家順.高原氣候環境對人體生理機能影響的研究[G].中國氣象學會2005年年會論文集, 2005: 4922-4924.

[20]譚斐, 楊映霞.廣西地面太陽輻射分布特征以及對人體健康的影響[J].氣象研究與應用, 2012, 33(2): 45-48.

（本文編輯：吳昔昔）

Correlative research on meteorological factors and incidence of Kawasaki disease in Wenzhou

ZHOU Ai-hua1, YUE Yayan2, JIN Zengyou1, PAN Lulu3, ZHANG Yingying1, WANG Dan1, QIAN Yan1, ZHANG Yuanhai3, WU Rongzhou3, CHU Maoping3.
1.Department of Paediatrics, the First Affi liated Hospital of Wenzhou Medical University, Wenzhou, 325000; 2.Department of Neurology, the Affi liated Hospital of Jiangsu University, Zhenjiang, 212000; 3.Department of Pediatric Cardiovascular Disease, the Second Affi liated Hospital & Yuying Children’s Hospital of Wenzhou Medical University, Wenzhou, 325027

Objective: To explore the correlation between meteorological factors and incidence of Kawasaki disease (KD) in Wenzhou.Methods: The medical records of 1 956 KD patients and meteorological factors in Wenzhou from January 2005 to December 2015 were retrieved.The correlation between the number of KD patients and the meteorological factors was studied.Results: The incidence rate of KD was high in spring and summer but low in winter.Amount of rainfall was a positive correlation factor on the incidence of KD (r=0.217, P＜0.05), while atmospheric pressure and wind speed were negative correlation factors (r=-0.209, -0.652, P＜0.05).Temperature and sunshine duration had no signifi cant correlation with KD incidence (P＞0.05).Conclusion: The results suggest that meteorological factors are correlated with the pathogenesis of KD.

Kawasaki disease; wind speed; rainfall

R725

A

10.3969/j.issn.2095-9400.2016.12.008

2016-07-09

溫州市科技局科研基金資助項目（Y20140056，Y2015 0015）；浙江省醫藥衛生科技項目（2016KYB197）。

周愛華（1980-），女，湖北荊門人，住院醫師，碩士。

褚茂平，主任醫師，碩士生導師，Email：chmping@126.com。