張文靜，葛旭峰，歐茜文，蔡雨陽(yáng)，楊玲，劉芳*

1.200092 上海市，上海交通大學(xué)醫(yī)學(xué)院附屬新華醫(yī)院老年科

2.201822 上海市嘉定區(qū)精神衛(wèi)生中心精神科

3.200025 上海市，上海交通大學(xué)醫(yī)學(xué)院公共衛(wèi)生學(xué)院

4.200434 上海市，同濟(jì)大學(xué)附屬上海市第四人民醫(yī)院老年科

近年來(lái)抑郁癥患者明顯增多［1］。據(jù)估計(jì)全球約3.22 億抑郁癥患者，占總?cè)丝诒壤?.4%，其致殘生活時(shí)間遠(yuǎn)高于心血管疾病致殘生活時(shí)間，引起了廣泛關(guān)注［2］。2021 年中國(guó)精神衛(wèi)生調(diào)查顯示，我國(guó)抑郁癥患病率為3.4%，患病人數(shù)為9 500 萬(wàn)人，每年大約有28萬(wàn)人自殺，其中40%患有抑郁癥［3］。多項(xiàng)研究顯示神經(jīng)精神疾病與溫度、濕度、氣壓、日照時(shí)間等氣象元素有關(guān)［4-6］。中國(guó)的二十四節(jié)氣“春雨驚春清谷天，夏滿芒夏暑相連，秋初露秋寒霜降，冬雪雪冬小大寒”，是智慧的祖先對(duì)一年自然節(jié)律變化進(jìn)行的精煉總結(jié)［7］。目前關(guān)于二十四節(jié)氣的溫濕指數(shù)（Humidex）與抑郁癥關(guān)系的研究相對(duì)較少。既往研究多關(guān)注單一指標(biāo)如平均氣溫、濕度、日照時(shí)長(zhǎng)等［3-6］，忽略多個(gè)因素間相互作用對(duì)抑郁癥患者產(chǎn)生的影響。本研究通過(guò)分析上海市嘉定區(qū)精神衛(wèi)生中心2016-10-08—2019-10-07 抑郁癥患者的數(shù)據(jù)和嘉定區(qū)氣象數(shù)據(jù)，全面系統(tǒng)研究二十四節(jié)氣的Humidex 與抑郁癥的關(guān)系，得出抑郁癥患者人數(shù)隨二十四節(jié)氣變化的趨勢(shì)，為減少抑郁癥患病率、采取相應(yīng)的預(yù)防措施提供理論依據(jù)。

1 資料與方法

1.1 資料收集

收集上海市嘉定區(qū)精神衛(wèi)生中心精神科2016-10-08—2019-10-07 抑郁癥患者的就診資料，包括性別、年齡、就診時(shí)間、就診門診及專家門診數(shù)量。納入標(biāo)準(zhǔn)：抑郁癥患者符合《國(guó)際疾病分類（第10 版）》抑郁癥（ICD-10 編碼為F32）的診斷標(biāo)準(zhǔn)。排除標(biāo)準(zhǔn)：（1）雙相情感障礙抑郁發(fā)作；（2）既往診有除抑郁癥外的其他精神病者；（3）器質(zhì)性病變導(dǎo)致的抑郁癥；（4）嚴(yán)重頭創(chuàng)傷史。通過(guò)上海市氣象局公布數(shù)據(jù)收集嘉定區(qū)2016-10-08—2019-10-07 日均溫度、相對(duì)濕度、降雨量、氣壓和平均風(fēng)速數(shù)據(jù)。Humidex 計(jì)算公式：-10］，其中tem 為日均溫度，hum 為相對(duì)濕度。納入的樣本量估算采用描述性研究樣本量的計(jì)算公式N=Z2×P×（1-P）/ε2，其中N 為樣本量；Z 為統(tǒng)計(jì)量，置信區(qū)間（CI）為95%時(shí)，Z=1.96；ε為誤差值；P為概率值。中國(guó)精神障礙患病率橫斷面流行病學(xué)研究［1］顯示抑郁癥終生患病率為3.4%。N=11 364.7，考慮人數(shù)為整數(shù)得出N=11 365。最終實(shí)際納入樣本18 289 人次。二十四節(jié)氣從寒露節(jié)氣開始，到秋分節(jié)氣結(jié)束，繪制二十四節(jié)氣與Humidex、Humidex 和抑郁癥患者數(shù)的關(guān)系圖。

1.2 分組

共納入患者18 289 人次，根據(jù)參考文獻(xiàn)，將患者依據(jù)不同年齡分為≤18 歲組（n=157）、19～44 歲組（n=3 099）、45～59 歲組（n=4 848）、60～74 歲組（n=6 270）和≥75 歲組（n=3 915）［8-9］。

1.3 統(tǒng)計(jì)學(xué)方法

采用SPSS 16.0 與R 4.2.0 軟件進(jìn)行數(shù)據(jù)分析，符合正態(tài)分布的計(jì)量資料以（±s）表示。通過(guò)結(jié)合類Poisson 分布和分布滯后非線性模型估計(jì)二十四節(jié)氣的Humidex 與抑郁癥之間的關(guān)系。運(yùn)用Pearson 相關(guān)性分析計(jì)算Humidex 和氣象變量之間的相關(guān)性。當(dāng)Humidex 與單一氣象變量的相關(guān)系數(shù)（r）＞0.7，被認(rèn)為是強(qiáng)相關(guān)，模型需排除此變量。模型如下：Log［E（Yt）］=α+βXt，l+ns（Wind speedt，df=3）+ns（Time，df=5/year）+DOWt+Holidayt+Expertst，其中E（Yt）表示第t 天抑郁癥患者的預(yù)期計(jì)數(shù)，α為回歸模型截距，X 表示Humidex，l 表示最大滯后天數(shù)；β表示Humidex 的系數(shù)；Wind speedt表示第t 天的風(fēng)速；Time 表示時(shí)間的變量；ns 是自然三次樣條函數(shù)。ns（Time，df）指的是利用自然樣條函數(shù)控制季節(jié)和長(zhǎng)期趨勢(shì)的效果，包括星期幾（Dow），節(jié)假日（Holiday）和專家門診（Experts）在內(nèi)的分類變量也被控制。參考相關(guān)文獻(xiàn)，對(duì)于風(fēng)速指標(biāo)，使用具有3 df 的自然三次樣條曲線，使用自然三次樣條曲線以每年5 df 的時(shí)間控制長(zhǎng)期趨勢(shì)和季節(jié)性分布［10-12］。每個(gè)節(jié)氣平均時(shí)間為15 d，并且根據(jù)赤池信息準(zhǔn)則（AIC），本研究選擇最大滯后天數(shù)為15 d。以中位Humidex 作為對(duì)照，計(jì)算二十四節(jié)氣的Humidex 和抑郁癥的相對(duì)風(fēng)險(xiǎn)度（RR）。同時(shí)，探討性別與年齡分組對(duì)效應(yīng)的修飾作用，計(jì)算不同人群效應(yīng)值差的95%CI。以P＜0.05 為差異有統(tǒng)計(jì)學(xué)意義。95%CI 檢驗(yàn)差異的顯著性。以P＜0.05 為差異有統(tǒng)計(jì)學(xué)意義。

2 結(jié)果

2.1 資料統(tǒng)計(jì)描述

2016-10-08—2019-10-07共納入抑郁癥患者18 289人次，男6 900 人次，女11 389 人次；研究期間嘉定區(qū)氣象因素分布見(jiàn)表1。2016-10-08—2017-10-07 小寒節(jié)氣對(duì)應(yīng)的平均Humidex 最低，抑郁癥患者數(shù)最多；大暑節(jié)氣對(duì)應(yīng)的平均Humidex 最高，抑郁癥患者數(shù)最少。2017-10-08—2018-10-07 大寒節(jié)氣對(duì)應(yīng)的平均Humidex 最低，抑郁癥患者數(shù)最多；大暑節(jié)氣對(duì)應(yīng)的平均Humidex 最高，抑郁癥患者數(shù)最少。2018-10-08—2019-10-07 大寒節(jié)氣對(duì)應(yīng)的平均Humidex 最低，抑郁癥患者數(shù)最多；大暑節(jié)氣對(duì)應(yīng)的平均Humidex 最高，抑郁癥患者數(shù)最少。2016-10-08—2019-10-07 二十四節(jié)氣的平均Humidex 和抑郁癥患者數(shù)分布情況見(jiàn)圖1。

表1 2016-10-08—2019-10-07 嘉定區(qū)氣象因素分布數(shù)據(jù)Table 1 Distribution of meteorological factors in Jiading District from 2016-10-08 to 2019-10-07

圖1 2016-10-08—2019-10-07 二十四節(jié)氣的平均Humidex 和抑郁癥患者數(shù)分布圖Figure 1 Distribution of average Humidex and number of depression patients in 24 solar terms from 2016-10-08 to 2019-10-07

2.2 分布滯后非線性模型變量納入情況

相關(guān)性分析結(jié)果顯示，Humidex 與日均溫度（r=0.899）、氣壓（r=-0.832）強(qiáng)相關(guān)（P＜0.05），故不納入模型；Humidex 與相對(duì)濕度（r=0.028）、風(fēng)速（r=-0.023）、降雨量（r=0.044）無(wú)相關(guān)性（P＞0.05）。此外，因Humidex 為日均溫度與相對(duì)濕度計(jì)算得到的指標(biāo)，故相對(duì)濕度不納入模型。因降雨量數(shù)據(jù)離散性大，故不納入模型。最終模型納入Humidex 和風(fēng)速2 種變量。

2.3 二十四節(jié)氣的Humidex 對(duì)抑郁癥的影響

以中位Humidex（21.5）作為比較基準(zhǔn)，二十四節(jié)氣的Humidex 與總?cè)巳阂钟舭Y風(fēng)險(xiǎn)呈非線性關(guān)系。隨著Humidex 逐漸升高，總?cè)巳阂钟舭Y風(fēng)險(xiǎn)分布呈現(xiàn)先下降后上升再下降的趨勢(shì)。2016-10-08—2019-10-07 總?cè)巳阂钟舭Y風(fēng)險(xiǎn)分布前2 位Humidex 對(duì)應(yīng)的節(jié)氣分別為大寒、立春；總?cè)巳阂钟舭Y風(fēng)險(xiǎn)分布最低Humidex 對(duì)應(yīng)的節(jié)氣均為大暑。繪制Humidex 與不同性別和年齡抑郁癥人群的劑量反應(yīng)圖，可見(jiàn)隨著Humidex 逐漸升高，不同人群抑郁癥風(fēng)險(xiǎn)分布均大致呈現(xiàn)先下降后上升再下降的趨勢(shì)。男性、女性、≤18 歲、19～44 歲、45～59歲、60～74 歲和≥75 歲抑郁癥風(fēng)險(xiǎn)分布前2 位Humidex對(duì)應(yīng)的節(jié)氣分別為大寒、立春，抑郁癥風(fēng)險(xiǎn)分布最低Humidex 對(duì)應(yīng)的節(jié)氣為大暑，見(jiàn)圖2。

圖2 2016-10-08—2019-10-07 二十四節(jié)氣的Humidex 與抑郁癥人群的劑量反應(yīng)圖Figure 2 Dose-response graph of Humidex and depression patients in twenty-four solar terms from 2016-10-08 to 2019-10-07

2.4 Humidex 對(duì)抑郁癥的滯后影響

抑郁癥總?cè)巳涸诖蠛?、立春?jié)氣后第10 天出現(xiàn)最大滯后效應(yīng)（RR=1.020，95%CI=1.001～1.040），在大暑節(jié)氣后第8 天出現(xiàn)滯后效應(yīng)，持續(xù)至第9 天。分析不同人群發(fā)現(xiàn)，僅60～74 歲年齡組人群在大寒、立春節(jié)氣存在滯后效應(yīng)，60～74 歲年齡組人群在大寒、立春節(jié)氣后第7 天出現(xiàn)滯后效應(yīng)，第9 天出現(xiàn)最大滯后效應(yīng)（RR=1.054，95%CI=1.007～1.104），持續(xù)至第11 天。僅男性和60～74 歲年齡組人群在大暑節(jié)氣存在滯后效應(yīng)。男性抑郁癥人群在大暑節(jié)氣后第1 天和第8 天出現(xiàn)滯后效應(yīng)，分別持續(xù)至第2 天和第10 天，于第10 天出現(xiàn)最大滯后效應(yīng)（RR=0.952，95%CI=0.911～0.994）；60～74 歲年齡組人群在大暑節(jié)氣后第8 天出現(xiàn)滯后效應(yīng)，持續(xù)至第9 天，見(jiàn)圖3。

圖3 二十四節(jié)氣Humidex 對(duì)抑郁癥人群影響的3D 圖Figure 3 3D figure of the effect of Humidex in twenty-four solar terms on depression in the total population

3 討論

本研究結(jié)果顯示，二十四節(jié)氣的Humidex 與抑郁癥呈非線性相關(guān)關(guān)系，總?cè)巳阂钟舭Y風(fēng)險(xiǎn)分布前2 位Humidex 對(duì)應(yīng)的節(jié)氣分別為大寒、立春，抑郁癥風(fēng)險(xiǎn)分布最低Humidex 對(duì)應(yīng)的節(jié)氣為大暑。二十四節(jié)氣具有重要的科學(xué)價(jià)值與文化價(jià)值，也體現(xiàn)了中國(guó)獨(dú)具特色的生產(chǎn)方式和生活態(tài)度［13］。1961—2014 年華北地區(qū)二十四節(jié)氣氣溫?cái)?shù)據(jù)顯示，大暑節(jié)氣平均氣溫最高（30.4 ℃），在大寒節(jié)氣最低（-7.5 ℃）［14］。本研究發(fā)現(xiàn)在大寒、立春節(jié)氣的Humidex 低，大暑節(jié)氣Humidex 高，Humidex 受溫度影響更大。本研究結(jié)果與光照治療抑郁癥的研究結(jié)果具有相同點(diǎn)，抑郁癥在冬季加重，夏季減輕［15］。中醫(yī)認(rèn)為肝陽(yáng)乃春生少陽(yáng)之氣，始萌未盛，易受戕伐。如寒邪正盛，肝陽(yáng)受損，肝失升發(fā)、暢達(dá)之性而處于瘀滯狀態(tài)［16］。在冬季抑郁癥患者的肝疏泄功能降低，更容易導(dǎo)致抑郁癥。

相關(guān)研究發(fā)現(xiàn)抑郁癥高發(fā)的原因可能與冬季日照時(shí)間改變和褪黑素水平變化有關(guān)［17］。有研究表明冬季抑郁癥高發(fā)［18］，本研究結(jié)果與之一致。有關(guān)臨床試驗(yàn)發(fā)現(xiàn)，季節(jié)性抑郁癥在使用光療時(shí)療效顯著［19］。神經(jīng)影像學(xué)研究表明光療具有抗抑郁作用，可通過(guò)視網(wǎng)膜的雙突觸回路影響抑郁樣行為，神經(jīng)回路受光照劑量多少、波長(zhǎng)長(zhǎng)短影響。研究發(fā)現(xiàn)，情緒受環(huán)境光照時(shí)間長(zhǎng)短、光照溫度高低和波長(zhǎng)長(zhǎng)短的影響，光照效果也與不同個(gè)體及主觀感受有關(guān)［20］。褪黑素通過(guò)多種途徑改善抑郁癥，對(duì)于治療抑郁癥起著非常重要的作用。褪黑素來(lái)源于腦松果體，可以調(diào)節(jié)晝夜節(jié)律，同時(shí)具有抗抑郁功能，其作用機(jī)制可能與單胺神經(jīng)遞質(zhì)，下丘腦-垂體-腎上腺軸等有關(guān)［21］。

本研究通過(guò)對(duì)不同性別和各年齡組進(jìn)行分層分析發(fā)現(xiàn)，男性和女性抑郁癥患者在二十四節(jié)氣的風(fēng)險(xiǎn)分布趨勢(shì)相同，大寒、立春為其風(fēng)險(xiǎn)分布最高的2 個(gè)節(jié)氣，大暑為其風(fēng)險(xiǎn)分布最低的節(jié)氣。不同年齡組抑郁癥患者在二十四節(jié)氣的風(fēng)險(xiǎn)分布趨勢(shì)相同，表現(xiàn)為大寒、立春節(jié)氣的風(fēng)險(xiǎn)分布最高，大暑節(jié)氣的風(fēng)險(xiǎn)分布最低。不同性別、年齡組抑郁癥患者最大滯后風(fēng)險(xiǎn)不同，滯后時(shí)間也不同。

本研究存在以下不足：（1）本研究針對(duì)上海市嘉定區(qū)精神衛(wèi)生中心抑郁癥患者門診就診數(shù)據(jù)進(jìn)行分析，而上海市屬于亞熱帶季風(fēng)性氣候，與中國(guó)其他地區(qū)結(jié)果可能存在地區(qū)差異，后續(xù)需要進(jìn)行全國(guó)范圍抑郁癥人群數(shù)據(jù)分析。（2）本研究探討二十四節(jié)氣的Humidex 的短期影響而非長(zhǎng)期影響；有必要在更長(zhǎng)時(shí)間范圍探討二十四節(jié)氣的Humidex 與抑郁癥人群的關(guān)系。（3）研究納入模型的抑郁癥患者可能存在定期前往嘉定區(qū)精神衛(wèi)生中心門診就診情況，就診原因可能與二十四節(jié)氣的Humidex 因素?zé)o關(guān)，但由于數(shù)據(jù)限制，無(wú)法對(duì)此種抑郁癥患者數(shù)據(jù)進(jìn)行排除，后續(xù)研究將會(huì)收集更多與個(gè)體特征有關(guān)數(shù)據(jù)進(jìn)行更全面分析，也需要收集更長(zhǎng)時(shí)間樣本數(shù)據(jù)進(jìn)行檢驗(yàn)。

綜上所述，大寒、立春節(jié)氣的Humidex 是抑郁癥的危險(xiǎn)因素，大暑節(jié)氣的高Humidex 是抑郁癥的保護(hù)因素，提示臨床上應(yīng)根據(jù)不同節(jié)氣的Humidex，對(duì)抑郁癥患者人群做到精準(zhǔn)診療，制定個(gè)體化的干預(yù)方案。

作者貢獻(xiàn)：張文靜、葛旭峰負(fù)責(zé)論文撰寫、數(shù)據(jù)收集、調(diào)查研究；歐茜文、蔡雨陽(yáng)負(fù)責(zé)統(tǒng)計(jì)分析；楊玲、劉芳對(duì)研究進(jìn)行監(jiān)督管理，對(duì)本文負(fù)責(zé)。

本文無(wú)利益沖突。