趙納 孟令惠 張英東 杜國帥 劉紅冉 高明龍 任貝貝 劉富榮

精神創(chuàng)傷性事件，又稱精神應激事件，包括親身經(jīng)歷或者目睹他人的死亡、遭受死亡威脅、嚴重傷害或性暴力等，影響著越來越多的人［1］。精神創(chuàng)傷性事件不僅是診斷創(chuàng)傷后應激障礙（posttraumatic stress disorder，PTSD）的先決條件，還與許多慢性疾病的發(fā)病有關，如糖尿病、高血壓、心臟疾病等［2］。所以，許多研究關注創(chuàng)傷后腦結構及功能的改變，尋找創(chuàng)傷后干預的生物學表征。但是，目前大部分研究關注的是大腦局部腦區(qū)或局部腦網(wǎng)絡連接的改變，缺乏從全腦水平大尺度分析腦網(wǎng)絡連接的變化。人腦網(wǎng)絡具有小世界屬性，腦區(qū)之間以低耗能的形式進行信息傳輸，具有高效的全腦和局部信息傳遞能力［3］。許多神經(jīng)精神疾病，包括PTSD均發(fā)現(xiàn)腦網(wǎng)絡的小世界屬性遭到破壞［4］。已有研究發(fā)現(xiàn)創(chuàng)傷后腦功能網(wǎng)絡有向隨機化發(fā)展的趨勢［5］，腦結構網(wǎng)絡的小世界屬性是否也發(fā)生了改變，目前尚不清楚。本研究采用基于彌散張量成像（diffusion tensor imaging，DTI）的圖論方法對交通事故受害者進行分析，在全腦水平探討急性創(chuàng)傷后早期腦結構網(wǎng)絡小世界屬性的變化，為尋找有效的創(chuàng)傷后早期干預方法提供理論基礎。

方法

1. 研究對象

招募2020年11月至2022年3月來自河北醫(yī)科大學第三醫(yī)院門診及住院的交通事故受害者31例（創(chuàng)傷組）和年齡、性別、受教育時間與創(chuàng)傷組相匹配的健康對照者31例（對照組）。創(chuàng)傷組由2名精神科醫(yī)生使用漢密爾頓焦慮量表（Hamilton anxiety scale，HAMA）進行焦慮測評，使用簡明國際神經(jīng)精神障礙訪談檢查表（Mini International Neuropsychiatric Interview，MINI）排除精神疾病。創(chuàng)傷組納入標準：①年齡大于18歲且小于60歲；②初中或以上文化程度；③右利手；④近1周內(nèi)親身經(jīng)歷交通事故，感受到自身生命受到威脅；⑤無頭外傷及意識喪失；⑥由精神科醫(yī)師使用MINI排除精神疾?。虎邿o神經(jīng)系統(tǒng)及其他器質(zhì)性疾??；⑧一級親屬中無精神疾病病史。對照組納入標準：①年齡大于18歲且小于60歲；②初中或以上文化程度；③右利手；④未經(jīng)歷嚴重精神創(chuàng)傷性事件；⑤無腦外傷史；⑥由精神科醫(yī)師使用MINI排除精神疾??；⑦無神經(jīng)系統(tǒng)及其他器質(zhì)性疾病；⑧一級親屬中無精神疾病病史。

排除標準：①嚴重的軀體疾病，既往神經(jīng)系統(tǒng)疾病史；②藥物及物質(zhì)濫用史；③女性處于妊娠、哺乳期；④有磁共振檢查禁忌證；⑤磁共振檢查有器質(zhì)性病變；⑥磁共振圖像質(zhì)量不符合要求。

本研究得到河北醫(yī)科大學第三醫(yī)院倫理委員會的批準，所有受試者均簽署知情同意書，并且在同一天內(nèi)完成臨床數(shù)據(jù)和影像學數(shù)據(jù)采集。創(chuàng)傷組數(shù)據(jù)采集時間在經(jīng)歷交通事故后1周內(nèi)。

2. 研究方法

2.1 影像數(shù)據(jù)采集

采用SIEMENS Verio 3.0 T 磁共振進行影像數(shù)據(jù)采集。使用12單元頭矩陣線圈。掃描過程中使用降噪耳機或棉球降低噪聲，同時使用泡沫軟墊固定頭部以防止頭動。磁共振掃描序列及參數(shù)如下：

（1）3D-T1序列：TR=1900 ms，TE=2.26 ms，視野=256 mm×256 mm，層厚1 mm，無間隔，層數(shù)176。

（2）DTI序列：TR=12800 ms，TE=85 ms， 視野=220 mm×220 mm，層厚2 mm，無間隔，層數(shù)65，b=1000 s/mm2，采集64個方向，1個b0。

為排除顱內(nèi)器質(zhì)性病灶，所有受試者均先行顱腦MRI平掃，未發(fā)現(xiàn)顱內(nèi)器質(zhì)性病變。

2.2 影像數(shù)據(jù)分析

使用PANDA軟件（http：//www.nitrc.org/projects/panda/）進行數(shù)據(jù)的預處理及網(wǎng)絡連接矩陣的建立。①去除頭皮；②頭動和渦流校正；③采用DTIFIT建立彌散張量模型，獲得各向異性分數(shù)（fractional anisotropy，F(xiàn)A）圖；④采用確定性纖維追蹤，使用自動解剖標記（Automated Anatomical Labeling，AAL）90圖集，建立全腦白質(zhì)結構網(wǎng)絡。

使用Gretna軟件（http：//www.nitrc.org/projects/gretna/）計算腦結構網(wǎng)絡屬性。

2.3 腦結構網(wǎng)絡屬性分析

（1）聚類系數(shù)（clustering coefficient，Cp）：衡量網(wǎng)絡的模塊化程度。值越大，說明網(wǎng)絡模塊化程度越高。

（2）標準聚類系數(shù)γ：聚類系數(shù)Cp與隨機網(wǎng)絡對應參數(shù)的比值，γ=Cnet/Crand（rand表示隨機網(wǎng)絡，net表示真實網(wǎng)絡）。

（3）最短路徑長度（shortest path length，Lp）：衡量網(wǎng)絡內(nèi)部信息的傳輸。網(wǎng)絡通過該路徑可以快速的傳遞信息，節(jié)省各種資源。

（4）標準最短路徑長度λ：最短路徑長度Lp與隨機網(wǎng)絡對應參數(shù)的比值，λ=Lnet/Lrand（rand表示隨機網(wǎng)絡，net表示真實網(wǎng)絡）。

（5）小世界性σ：同一個網(wǎng)絡中標準聚類系數(shù)與標準最短路徑長度的比值。σ=γ/λ。

（6）全局效率（global efficiency，Eg）：衡量網(wǎng)絡整體的信息傳輸能力。

（7）局部效率（local efficiency，Eloc）：衡量網(wǎng)絡局部的信息傳輸效率。

同時滿足γ＞1且λ≈1的網(wǎng)絡屬于小世界網(wǎng)絡［3］。

3. 統(tǒng)計學分析

使用SPSS21.0統(tǒng)計軟件對2組受試者的性別、年齡、受教育時間等資料進行統(tǒng)計學檢驗。使用獨立樣本t檢驗比較2組的年齡、受教育時間；性別比的比較采用χ2檢驗。腦結構網(wǎng)絡指標使用Gretna軟件，采用獨立樣本t檢驗進行分析?？刂颇挲g和受教育程度，采用偏相關分析研究HAMA評分與創(chuàng)傷組腦結構網(wǎng)絡指標的相關性。上述假設檢驗顯著性閾值設為0.05。

結果

1. 創(chuàng)傷組和對照組的人口學資料

創(chuàng)傷組31例和健康對照組31例被納入研究，2組之間年齡、性別比、受教育時間相匹配，組間比較差異無統(tǒng)計學意義（P＞0.05）。詳見表1。

表1 創(chuàng)傷組和對照組的人口學資料

2. 創(chuàng)傷組腦結構網(wǎng)絡屬性的改變

創(chuàng)傷組與對照組均滿足γ＞1且λ≈1，表明2組受試者的大腦結構網(wǎng)絡均具有小世界屬性。創(chuàng)傷組的σ值高于對照組，2組間差異具有統(tǒng)計學意義（t=3.037，P=0.004）；創(chuàng)傷組的γ值高于對照組，2組間差異具有統(tǒng)計學意義（t=2.980，P=0.004）。創(chuàng)傷組與對照組的Cp、Lp、λ、Eg、Eloc值間差異無統(tǒng)計學意義（P＞0.05），見表2。此外，創(chuàng)傷組腦結構網(wǎng)絡屬性未見與焦慮評分相關。

表2 創(chuàng)傷組與對照組腦結構網(wǎng)絡屬性比較

討論

本研究應用DTI構建創(chuàng)傷后早期的腦白質(zhì)結構網(wǎng)絡，從全腦尺度出發(fā)，采用圖論方法，探討創(chuàng)傷后早期腦白質(zhì)結構網(wǎng)絡的變化。研究結果顯示創(chuàng)傷組和對照組的腦白質(zhì)結構網(wǎng)絡均具有小世界屬性，創(chuàng)傷組的σ和γ值較對照組升高。

基于彌散磁共振成像和纖維追蹤技術構建的腦白質(zhì)結構網(wǎng)絡分析已經(jīng)廣泛應用于各種疾病［6-8］，揭示其腦網(wǎng)絡拓撲屬性，為進一步理解疾病的發(fā)病機制及早期診斷提供了新的線索。與基于靜息態(tài)fMRI的腦功能網(wǎng)絡反映不同腦區(qū)之間神經(jīng)活動在時間上的關聯(lián)性不同，腦結構網(wǎng)絡可直觀地反映腦區(qū)之間真實的結構連接［9］。腦結構網(wǎng)絡和腦功能網(wǎng)絡的拓撲屬性具有相關性，又有差異性。有研究發(fā)現(xiàn)，結構網(wǎng)絡比功能網(wǎng)絡具有更高的全局效率和局部效率［10］；功能連接較高的腦區(qū)，它們之間并不一定存在直接的結構連接［11］。所以，結合多模態(tài)成像，全面定量評價腦網(wǎng)絡的拓撲屬性仍然是當前研究的一個重要方向。

人腦網(wǎng)絡是高效的網(wǎng)絡系統(tǒng)，它能夠?qū)⒉煌墓δ苓M行分化與整合，具有較高的Cp和較小的Lp，即小世界屬性，這使得大腦以較低的能耗實現(xiàn)高效的局部和全局水平上的信息傳輸［9］。σ=γ/λ，衡量網(wǎng)絡的小世界特性，其中γ衡量網(wǎng)絡局部信息傳輸能力，λ衡量網(wǎng)絡的全局傳輸能力，同時滿足γ＞1且λ≈1的網(wǎng)絡屬于小世界網(wǎng)絡［3］。本研究顯示創(chuàng)傷組和對照組的γ＞1且λ≈1，說明其腦白質(zhì)結構網(wǎng)絡均具有小世界屬性。創(chuàng)傷組與對照組相比，其σ值和γ值升高，提示創(chuàng)傷組局部信息的處理能力可能有所增強，一定程度上反映了腦網(wǎng)絡的防御能力提高。研究表明，腦結構網(wǎng)絡可以根據(jù)行為和腦功能的需求出現(xiàn)動態(tài)變化，尤其是在需要快速應對外部變化的關鍵時期［9］。創(chuàng)傷暴露者在經(jīng)歷創(chuàng)傷事件后短時間內(nèi)就會產(chǎn)生痛苦的體驗，出現(xiàn)焦慮等癥狀；同時，創(chuàng)傷相關腦神經(jīng)環(huán)路的功能也會出現(xiàn)改變［12］，所以，此時腦結構網(wǎng)絡小世界性的增高，有可能是創(chuàng)傷后早期的一種應激反應，以增強腦網(wǎng)絡的防御能力。

但是，本研究的結果與前期腦功能網(wǎng)絡的研究結果并不完全一致，Du等［5］研究發(fā)現(xiàn)，地震后25天創(chuàng)傷者的腦功能網(wǎng)絡也具有小世界性，但較對照者有向隨機化發(fā)展的趨勢。分析其原因可能有如下幾點：第一，可能與創(chuàng)傷類型及創(chuàng)傷程度有關，既往研究表明不同的創(chuàng)傷類型及創(chuàng)傷程度引起腦結構及功能的改變不同［13-14］。本研究從醫(yī)學倫理方面考慮，沒有納入造成嚴重損傷的車禍受害者，納入人群創(chuàng)傷程度相對較輕，但Du等納入的是2008年汶川地震震源區(qū)的幸存者，他們親身經(jīng)歷了8.0級的地震，目睹了親人或者朋友的死亡，創(chuàng)傷程度很重。第二，可能與磁共振模態(tài)有關。采用不同模態(tài)的成像技術對大腦網(wǎng)絡拓撲屬性的分析可能會得到截然不同的結果［9］。例如，有研究報道精神分裂癥的腦結構網(wǎng)絡與腦功能網(wǎng)絡改變就截然相反［15］。不同模態(tài)的影像數(shù)據(jù)反映的腦網(wǎng)絡機制不同，功能腦網(wǎng)絡根據(jù)某一時刻腦活動的同步性構建腦網(wǎng)絡，而結構腦網(wǎng)絡則反映較穩(wěn)定的解剖結構。腦結構網(wǎng)絡和腦功能網(wǎng)絡拓撲屬性之間錯綜復雜的聯(lián)系還需要進一步研究。

本研究仍然存在著一些不足之處。第一，DTI的空間分辨率明顯大于腦白質(zhì)神經(jīng)纖維的實際大小，前者是毫米，而后者只有微米大小。第二，腦結構網(wǎng)絡構建采用確定性纖維追蹤技術，難以重建交叉纖維以及較長的纖維，導致腦區(qū)之間部分連接缺失，未來需要更先進的彌散成像技術來提高纖維追蹤的準確性。第三，本研究的樣本量較少，較大的樣本量有助于探討腦白質(zhì)結構網(wǎng)絡拓撲屬性與臨床癥狀之間的關系。第四，采用不同模態(tài)的成像技術可能會影響大腦網(wǎng)絡拓撲屬性的分析結果，未來需要融合多模態(tài)成像技術分析來理解神經(jīng)精神疾病的病理生理機制。

本研究從全腦大尺度出發(fā)，基于DTI采用圖論方法，發(fā)現(xiàn)創(chuàng)傷后早期腦結構網(wǎng)絡小世界性提高，加深了對創(chuàng)傷后大腦信息處理模式的理解，為揭示創(chuàng)傷后腦的病理生理機制提供了新的啟示。