2型糖尿病周圍神經(jīng)病變風(fēng)險的列線圖預(yù)測模型研究

李永生，張學(xué)良，李丞，封志煒，王凱*

2型糖尿?。═2DM）是世界范圍內(nèi)公認的重大的公共衛(wèi)生問題，2019年國際糖尿病聯(lián)盟（IDF）的糖尿病地圖（IDF Diabetes Atlas）（第9版）表明，全球約有4.63億糖尿病患者，估計2045年糖尿病患者將升至7.002億［1］。隨著糖尿病發(fā)生率的增加，并發(fā)癥的發(fā)生率也相應(yīng)增加［2］。據(jù)報道，作為糖尿病最常見的慢性并發(fā)癥，糖尿病周圍神經(jīng)病變（DPN）的患病率在中國為 35.34%［3］，緬甸為 33.7%［4］，印度為 29.2%［5］，也門為 56.2%［6］，約旦為 39.5%［7］，伊朗為 53%［8］。DPN是一種由長期高血糖引起的神經(jīng)病變［9-10］，可導(dǎo)致足部潰瘍，嚴重時需截肢，其致殘率、病死率極高，死亡風(fēng)險甚至超過部分癌癥，如乳腺癌和前列腺癌［11］。又因DPN發(fā)病機制復(fù)雜，目前尚無有效的治療方法。而DPN發(fā)病早期臨床癥狀缺乏特異性，多數(shù)患者一旦出現(xiàn)肢體麻木、疼痛等癥狀，周圍神經(jīng)已經(jīng)發(fā)生不可逆轉(zhuǎn)的病理變化［12］。如能早期發(fā)現(xiàn)和治療，可以延緩或預(yù)防DPN的不良后果，從而提高患者的生活質(zhì)量。因此，迫切需要及早發(fā)現(xiàn)并適當(dāng)控制可改變的DPN危險因素，以便早期預(yù)防和治療DPN。現(xiàn)有關(guān)于預(yù)測DPN的國內(nèi)外研究多基于社區(qū)型小樣本人群［13-16］，代表性方面有待商榷。本研究擬通過對新疆15 020例T2DM患者的臨床指標進行分析，探討DPN的危險因素，并以DPN作為結(jié)局變量，利用經(jīng)濟、有效且易于獲取的參數(shù)，在中國新疆地區(qū)建立一個精確的個性化DPN風(fēng)險預(yù)測模型，為臨床工作者在DPN的早期篩查和診斷治療方面提供參考。

1 資料與方法

1.1 研究對象 回顧性收集2010—2019年在新疆醫(yī)科大學(xué)第一附屬醫(yī)院就診的15 020例T2DM患者的住院資料，T2DM和DPN的診斷均依據(jù)《中國2型糖尿病防治指南（2010年版）》［17］中的診斷標準，本研究已獲新疆醫(yī)科大學(xué)第一附屬醫(yī)院倫理委員會批準（批號：K202108-24）。納入標準：（1）年齡≥18歲；（2）確診或新診斷為T2DM；（3）能夠自主交流，無精神障礙。排除標準：（1）資料不全；（2）未簽署知情同意書。

1.2 數(shù)據(jù)收集

1.2.1 一般資料 包括性別、年齡、糖尿病視網(wǎng)膜病變（DR）、糖尿病病程等。

1.2.2 生化資料 入院次日空腹狀態(tài)下取靜脈血3 ml，采用全自動生化分析儀（以3 000 r/min離心10 min，離心半徑22 cm）檢測白細胞計數(shù)（WBC）、嗜酸粒細胞分數(shù)、淋巴細胞計數(shù)、糖化血紅蛋白（HbA1c）、糖化血清蛋白（GSP）、總膽固醇（TC）、三酰甘油（TG）、高密度脂蛋白（HDL）、低密度脂蛋白（LDL）、載脂蛋白A1（ApoA1）、載脂蛋白B（ApoB）。

1.3 統(tǒng)計學(xué)方法 采用R 4.0.4軟件進行統(tǒng)計學(xué)分析，連續(xù)變量采用Kolmogorov-Smimovz正態(tài)檢驗，符合正態(tài)分布的計量資料采用（±s）描述，兩組間比較采用t檢驗；計數(shù)資料的分析采用χ2檢驗；通過嚴格的數(shù)據(jù)過濾與預(yù)處理，將符合條件的患者（n=15 020）按3∶1的比例隨機分為訓(xùn)練組（n=11 265）和驗證組（n=3 755）。以訓(xùn)練組數(shù)據(jù)中DPN作為結(jié)局變量，采用最小絕對收縮與選擇算子方法（Lasso回歸）［18-19］篩選出獨立預(yù)測因子，在此基礎(chǔ)上利用多因素Logistic回歸分析進一步探討并建立列線圖（Nomogram）預(yù)測模型，使用驗證組數(shù)據(jù)驗證構(gòu)建的預(yù)測模型，為了評價預(yù)測模型的效能，本研究采用受試者工作特征（ROC）曲線下面積驗證預(yù)測模型預(yù)測效果［20］，并通過Hosmer-Lemeshow擬合優(yōu)度檢驗判斷模型的擬合優(yōu)度，繪制校正曲線［21］，同時采用決策曲線分析（DCA）［22-23］預(yù)測DPN發(fā)生的風(fēng)險。以P＜0.05為差異有統(tǒng)計學(xué)意義。

2 結(jié)果

2.1 一般資料 15 020例T2DM患者中包括8 887例無DPN患者和6 133例DPN患者。無DPN患者和DPN患者嗜酸粒細胞分數(shù)、TC、LDL比較，差異均無統(tǒng)計學(xué)意義（P＞0.05）。無DPN患者和DPN患者性別、年齡、DR發(fā)生率、糖尿病病程、WBC、淋巴細胞計數(shù)、HbA1c、GSP、TG、HDL、ApoA1、ApoB比較，差異均有統(tǒng)計學(xué)意義（P＜0.05），見表1。

2.2 DPN的危險因素 以訓(xùn)練組是否發(fā)生DPN（賦值：是=1，否=0）為因變量，以性別（賦值：男=1，女=2）、年齡、DR（賦值：是=1，否=0）、糖尿病病程、WBC、嗜酸粒細胞分數(shù)、淋巴細胞計數(shù)、HbA1c、GSP、TC、TG、HDL、LDL、ApoA1、ApoB作為自變量（其余自變量賦值為連續(xù)變量）進行Lasso回歸分析（圖1），隨著懲罰系數(shù)λ的變化，模型納入的變量也在逐漸減少，最終選擇10倍交叉驗證誤差為最小λ+1（λ最小值的1 SE=0.023）作為模型的最優(yōu)值，此時篩選的預(yù)測變量包括年齡、DR、糖尿病病程、HbA1c和HDL。將篩選后的預(yù)測變量納入多因素Logistic回歸分析，結(jié)果顯示，年齡、DR、糖尿病病程、HbA1c和HDL是2型糖尿病患者發(fā)生DPN的影響因素（表2）。

表2 T2DM患者發(fā)生DPN影響因素的多因素Logistic回歸分析Table 2 Multivariate Logistic regression analysis of influencing factors of diabetic peripheral neuropathy in type 2 diabetics

圖1 采用Lasso回歸進行臨床特征的篩選Figure 1 Selection of potential predictors of diabetic peripheral neuropathy in type 2 diabetics by the least absolute shrinkage and selection operator（Lasso）regression

2.3 構(gòu)建列線圖預(yù)測模型 將上述多因素Logistic回歸分析結(jié)果篩選出的變量納入列線圖預(yù)測模型，結(jié)局指標選取DPN的發(fā)病風(fēng)險，繪制列線圖（圖2），根據(jù)每個風(fēng)險因素所對應(yīng)列線圖上方的標尺，從而得到該因素的單項評分，所有風(fēng)險因素評分相加得到總分，便可得到對應(yīng)患者DPN發(fā)生率，總分越高，意味發(fā)生DPN風(fēng)險的可能性越大。

圖2 T2DM患者發(fā)生DPN風(fēng)險預(yù)測列線圖Figure 2 Nomogram predicting diabetic peripheral neuropathy in type 2 diabetics

2.4 訓(xùn)練組和驗證組發(fā)生DPN的預(yù)測價值 繪制列線圖預(yù)測精度的ROC曲線，訓(xùn)練組列線圖預(yù)測模型的ROC曲線下面積為0.858〔95%CI（0.851，0.865）〕，診斷臨界值為0.464，靈敏度為76.5%，特異度為83.7%；通過驗證組數(shù)據(jù)集進行驗證，預(yù)測列線圖的ROC曲線下面積為0.852〔95%CI（0.840，0.865）〕（圖3、4），診斷臨界值為0.522，靈敏度為74.0%，特異度為87.4%。此外，驗證組的ROC曲線下面積較建模人群僅降低0.006，說明預(yù)測模型在訓(xùn)練組和驗證組人群均具有良好的預(yù)測區(qū)分度。

圖3 列線圖預(yù)測模型預(yù)測訓(xùn)練組T2DM患者發(fā)生DPN風(fēng)險的ROC曲線Figure 3 ROC curve of the predictive nomogram for the risk of diabetic peripheral neuropathy in type 2 diabetics of training group

圖4 列線圖預(yù)測模型預(yù)測驗證組T2DM患者發(fā)生DPN風(fēng)險的ROC曲線Figure 4 ROC curve of the predictive nomogram for the risk of diabetic peripheral neuropathy in type 2 diabetics of verification group

2.5 訓(xùn)練組和驗證組發(fā)生DPN的校正曲線 在訓(xùn)練數(shù)據(jù)集上，列線圖的校正曲線顯示預(yù)測結(jié)果與觀察吻合較好（圖5）。Hosmer-Lemeshow擬合優(yōu)度檢驗表明模型不顯著（P＞0.05），表明模型與觀察數(shù)據(jù)吻合較好。通過驗證組數(shù)據(jù)集進行驗證，同樣顯示預(yù)測與觀察吻合較好，Hosmer-Lemeshow擬合優(yōu)度檢驗表明模型不顯著（P＞0.05），表明模型與觀察數(shù)據(jù)吻合較好。

圖5 列線圖預(yù)測模型預(yù)測T2DM患者發(fā)生DPN風(fēng)險的校準曲線Figure 5 Calibration curve for predicting the risk of diabetic peripheral neuropathy in type 2 diabetics by the predictive nomogram

2.6 訓(xùn)練組和驗證組發(fā)生DPN的決策性曲線（DCA）在列線圖預(yù)測模型的基礎(chǔ)上，將篩選出的變量進行DPN的DCA（圖6），結(jié)果顯示，當(dāng)患者的閾值概率為0～0.9，使用列線圖預(yù)測T2DM發(fā)生DPN風(fēng)險的凈收益更高，同樣在驗證組得到證實（圖6），因此，廣泛的備選閾值概率表明該模型是一個很好的評估工具。

圖6 列線圖預(yù)測模型預(yù)測T2DM患者發(fā)生DPN風(fēng)險的決策性曲線分析Figure 6 Decision curve analysis for predicting the risk of diabetic peripheral neuropathy in type 2 diabetics by the predictive nomogram

3 討論

本研究結(jié)果發(fā)現(xiàn)年齡、DR、糖尿病病程、HbA1c和HDL是T2DM發(fā)生DPN的影響因素，其中，DPN發(fā)生率與年齡、是否發(fā)生DR、糖尿病病程和HbA1c呈正相關(guān)，與HDL呈負相關(guān)。

本研究結(jié)果顯示，DR是T2DM患者發(fā)生DPN的危險因素，其本身也是糖尿病微血管并發(fā)癥之一，雖然DPN和DR的病理生理機制有些不同，但其病理基礎(chǔ)均與糖代謝紊亂、微血管疾病及微循環(huán)障礙有關(guān)［24］，因此，在T2DM患者中，隨著DR的患病率增加，DPN的患病率也隨之增加，二者成平行關(guān)系［25］。在國內(nèi)外相關(guān)DPN風(fēng)險預(yù)測模型研究中［14-16］，少有文獻提出將DR納入DPN風(fēng)險預(yù)測模型的相關(guān)報道，而本研究結(jié)果顯示二者具有非常強的相關(guān)性，因而，在臨床工作中，醫(yī)務(wù)工作者對T2DM患者進行DPN篩查的過程中應(yīng)仔細詢問其是否患有DR。目前，DPN的發(fā)病機制尚不明晰，但已提出炎癥、氧化應(yīng)激和線粒體功能障礙主要參與DPN病理改變［26］。所有這些變化與衰老過程有關(guān)［27］。DPN已被證實與一些可改變和不可改變的危險因素有關(guān)［13，28-30］。在多數(shù)DPN的流行病學(xué)研究中，年齡和糖尿病病程是報道最多的不可改變的危險因素［2，31-32］。LI等［29］研究顯示，20～34 歲、35～49 歲、50～64歲和≥65歲的患者DPN發(fā)生率分別為8.4%、22.7%、33.0%和42.4%，年齡與DPN發(fā)生率顯著相關(guān)。POPESCU等［33］研究結(jié)果顯示，根據(jù)密歇根神經(jīng)病變篩查量表（MNSI）評分，DPN患病率為28.8%，與較高的年齡呈顯著正相關(guān)，證明了年齡影響DPN的發(fā)生；而SHUN等［34］和ZHANG等［35］發(fā)現(xiàn)，糖尿病病程越長，DPN患病率越高，可能是DPN患者神經(jīng)纖維病變導(dǎo)致纖維神經(jīng)密度的減少，加深了皮膚去神經(jīng)化的程度，進而增加了發(fā)生DPN患病風(fēng)險，因此在臨床工作中，需要提高對病程較長的T2DM患者的篩查力度，以減緩或減少DPN的發(fā)生，從而提升患者生活質(zhì)量。

本研究結(jié)果顯示，除年齡、DR和糖尿病病程外，血糖、血脂與DPN的發(fā)生同樣密切相關(guān)。T2DM患者常伴有脂代謝紊亂，表現(xiàn)為HDL水平降低，LDL水平升高。楊輝等［36］研究顯示HDL低水平是T2DM患者發(fā)生DPN的危險因素，其原因可能是HDL對逆向轉(zhuǎn)運膽固醇起到了促進作用，加速了TG和膽固醇代謝，從而穩(wěn)定血糖濃度。另有研究證明，HDL與神經(jīng)傳導(dǎo)速度呈正相關(guān)，DPN患者HDL水平明顯降低，神經(jīng)傳導(dǎo)速度受限提示HDL可能是DPN的危險因素［37］。本研究結(jié)果顯示，HDL與DPN呈負相關(guān)，起到了保護作用。HbA1c是葡萄糖與血紅蛋白反應(yīng)產(chǎn)物，能反映患者過去1～2個月血糖水平，李美［38］在研究中指出，DPN患者的HbA1c水平明顯高于非DPN患者，HbA1c與DPN發(fā)病率有關(guān)。HbA1c水平越高，DPN患病率也越高［39］，本研究結(jié)果與之一致，HbA1c水平與DPN呈顯著正相關(guān)，HbA1c升高是DPN發(fā)生的獨立危險因素之一，有研究發(fā)現(xiàn)長期高血糖可促使神經(jīng)髓鞘上的蛋白非酶促糖基化，導(dǎo)致神經(jīng)纖維信號傳導(dǎo)功能障礙［40］，而HbA1c水平的升高便表明T2DM患者處于高血糖狀態(tài)［41］，因此針對T2DM患者，需要嚴格控制其血糖，降低HbA1c水平，從而可以有效預(yù)防、延緩DPN的發(fā)生。

列線圖是臨床事件個體化預(yù)測分析的統(tǒng)計模型。與其他預(yù)測性統(tǒng)計學(xué)方法相比，列線圖分析可以通過直觀、可視化的方式提供更好的個體化預(yù)測風(fēng)險評估［14，42］。本研究通過Lasso回歸結(jié)合多因素Logistic回歸分析確定篩選出的5個危險因素是臨床醫(yī)生容易獲得的常規(guī)臨床變量，從而使其易于在實踐中采用。此外，使用隨機抽樣將符合條件的患者分為訓(xùn)練組和驗證組，以評價預(yù)測模型的效能。在訓(xùn)練組中，預(yù)測列線圖的ROC曲線下面積為0.858，相比SANTOS等［14］建立的模型（ROC曲線下面積=0.73）及寧改君等［16］建立的模型（ROC曲線下面積=0.789）具有更高的效能。將建立的預(yù)測模型在外部驗證時，不論是預(yù)測精度，還是校正曲線和DCA均具有良好的效能，由此說明該預(yù)測模型是一個很好的評估工具。本研究的局限性在于本研究為單中心研究，并非多中心大樣本流行病學(xué)調(diào)查結(jié)果，數(shù)據(jù)來自同一醫(yī)療中心，亦存在病例選擇偏倚，在后續(xù)的研究工作中，本課題組計劃與其他中心進行合作，利用其臨床數(shù)據(jù)資料進一步完善并提高模型的預(yù)測價值。

作者貢獻：李永生、張學(xué)良、李丞、封志煒和王凱進行文章的構(gòu)思與設(shè)計；張學(xué)良、李丞進行研究的實施與可行性分析；李丞、封志煒進行數(shù)據(jù)收集；李永生、封志煒進行數(shù)據(jù)整理；李永生、王凱進行統(tǒng)計學(xué)處理；李永生、張學(xué)良進行結(jié)果的分析與解釋；李永生撰寫論文；張學(xué)良進行論文的修訂；王凱負責(zé)文章的質(zhì)量控制及審校；張學(xué)良、王凱對文章整體負責(zé)，監(jiān)督管理。

本文無利益沖突。