白細胞介素相關基因多態性與塵肺易感性系統評價與meta分析

鄭 玲，蘇晶瑩，伍定輝，姚向陽

(1.福建醫科大學第三臨床醫學院, 福建 福州 350000；2.廈門大學附屬第一醫院 肺科，福建 廈門 361003；3.廈門大學附屬第一醫院杏林分院 兒科，福建 廈門 361003)

塵肺病是一種在職業活動中長期吸入生產性礦物性粉塵并在肺內潴留而引起的以肺組織彌漫性纖維化為主的疾病[1]，是目前我國職業病中最常見和危害最嚴重的疾病[1-2]。根據全國職業報告，我國職業性塵肺病新易感例2014年報告26 873例，2018年報告19 468例，2020年報告14 367例，分別占當年新發職業病病例總數89.66%、82.85%、84.19%，2020年因塵肺病死亡6 668例[3-5]。可見我國職業性塵肺病新易感例仍在全國各類職業病病例中居高不下。塵肺病早期可無明顯癥狀，但晚期常常導致殘疾甚至死亡。塵肺病的易感原因有很多，包括工作場所特征、宿主免疫反應和基因環境的相互作用等[6-7]。

近幾年來，基因多態性與塵肺病易感性的研究越來受關注[8-10]。白細胞介素(interleukin，IL)是指在白細胞或免疫細胞間相互作用的一種細胞因子。細胞因子在廣泛的生物過程中發揮作用，特別是在炎癥和免疫反應中，是在接觸二氧化硅和煤塵等礦物粉塵的人類中觀察到的有毒和致病作用的重要介質[11]。許多研究表明，IL相關基因多態性與人類許多疾病包括肺結核[12-14]、惡性腫瘤[15-17]、哮喘[18-19]、類風濕性關節炎[20-21]、膿毒癥[22-23]及塵肺[24-25]等的發生有關。近年來，有關研究對IL-1[26-27]、IL-17[24-25]等相關基因與塵肺易感的遺傳易感性進行了相關探討，然而研究結果尚存在爭議。王煒等[28]曾對IL-1基因多態位點與塵肺病易感性進行系統評價，但該研究排除了少于3篇同時研究一個IL-1多態位點的文獻，且僅對顯性模型進行meta分析。郝艷星等[29]曾對IL-6基因-174C/G和-634C/G多態性與塵肺病易感性進行系統評價，但研究基因有限。目前關于IL相關基因與塵肺病易感風險的大樣本量相關研究較少，本文將通過全面收集相關文獻進行綜合分析IL相關基因多態性與塵肺病的易感關聯風險。

1 資料與方法

1.1文獻檢索策略

1.1.1計算機檢索數據庫 檢索PubMed、the Cochrane Library、Embase、中國知網(CNKI)、WanFang Data、VIP、CBM等數據庫，同時手工檢索，并回溯納入文獻的參考文獻，檢索時間限定為建庫至2022年1月8日。檢索方式為主題詞和自由詞相結合。

1.1.2中文檢索詞 塵肺，矽肺，煤工塵肺，煤矽肺，石棉肺，肺塵埃沉著病，白介素，白細胞介素，IL，基因。以萬方數據庫為例檢索式：(塵肺OR矽肺OR 煤工塵肺OR石棉肺OR煤矽肺OR肺塵埃沉著病)AND (白介素OR白細胞介素OR IL)AND 基因。

1.1.3英文檢索詞 “pneumoconiosis”、“Pneumoconiose”、“Bagassosis”、“coalworker pneumoconiosis”、“CWP”、“silicosis”、“asbestosis” 、“lung fibrosis”、“interleukins”、“cytokine”、“interleukin”、“IL”、“Polymorphism, Genetic”、“Genetic Polymorphism”、“Gene Polymorphism”、“Polymorphism, Gene”、“Polymorphism (Genetics)”。以PubMed為例檢索式：(“pneumoconiosis”[MeSH Terms]) OR “Pneumoconiose”[All Fields]OR“Bagassosis”[All Fields]OR“coalworker pneumoconiosis”[All Fields]OR“CWP”[All Fields]OR“silicosis ”[All Fields]OR“asbestosis”[All Fields]OR“lung fibrosis”[All Fields]AND(“interleukins”[MeSH Terms]) OR“cytokine”[All Fields]OR“interleukin”[All Fields]OR“IL”[All Fields]) AND(“Polymorphism, Genetic”[MeSH Terms]OR“Genetic Polymorphism”[All Fields]OR“Gene Polymorphism”[All Fields]OR“Polymorphism, Gene”[All Fields]OR“Polymorphism (Genetics)”[All Fields])。

1.2文獻篩選標準

1.2.1文獻納入標準 (1)研究類型為病例對照研究；(2)研究的病例組均為塵肺人群，均符合塵肺診斷標準，對照組為接觸職業性粉塵的健康人群；(3)研究給出明確的病例組樣本量和對照組樣本量；(4)研究的指標為IL相關基因，從文獻中可直接或間接獲得不同基因型、等位基因的分布數據；(5)研究語種為中英文。

1.2.2文獻排除標準 (1)合并其他疾?。喝缃诟腥?、免疫系統疾病、遺傳代謝疾病等的患者；(2)研究數據無法完整獲取的文獻；(3)重復發表的文獻或數據；(4)綜述等二手文獻、個案報告。

1.3文獻篩選、資料提取與質量評價

1.3.1文獻篩選 由鄭玲和蘇晶瑩分別獨立完成文獻篩選、資料提取和文獻質量評價，并交由伍定輝主任醫師核對，出現意見分歧時共同協商討論解決或由主任醫師分析判定。通過閱讀文獻標題和摘要排除明顯不符的文獻，再對初篩后的文獻通過閱讀全文進一步排查。

1.3.2文獻資料提取 納入文獻的提取資料包括：(1)一般資料：第一作者、發表年份、塵肺種類、國家、種族、樣本量(病例組/對照組)、檢測的單核苷酸多態性(SNPs)、基因檢測方法、病例組的年齡和接塵時間、對照組的年齡和接塵時間等。(2)結局測量指標：病例組和對照組的各基因型數和等位基因數。

1.3.3文獻質量評價 由鄭玲和蘇晶瑩參照遺傳關聯性研究報告規范行質量評價[30]以及依據文獻[31]選取的14條標準對遺傳關聯性研究進行質量評價, 以“是”(滿足標準) 、“否”(不滿足標準) 和“不清楚”(文獻未描述) 進行文獻質量評價。由2名評價員獨立完成，意見不統一時，與伍定輝主任醫師或姚向陽主任醫師共同討論，直到達成一致意見。

1.4統計學方法 應用Review Manager 5.2以及Stata 14.0軟件進行統計學數據分析。關聯強度用風險比(odds ratio，OR)表示，效應量用點估計及其95%可信區間(confidence intervals，CI)表示；異質性分析采用Q檢驗與I2檢驗判斷，若I2≤50%，表示無統計學上異質性，采用固定效應模型；若I2>50%，表示存在統計學上異質性，則進一步分析異質性來源，在排除明顯臨床異質性的影響后，采用隨機效應模型。明顯的臨床異質性采用亞組分析、敏感性分析或meta回歸分析等方法進行處理。采用漏斗圖、Egger’s法檢驗分析評估文獻的發表偏倚。以P<0.05為差異有統計學意義。

2 結 果

2.1文獻檢索結果 初步檢索獲得1 486篇相關文獻，102篇文獻進入全文閱讀篩選，最終29篇文獻[24-27，32-56]共包含33組數據符合納入和排除標準進入本研究。共計病例組5 315例，對照組5 332例。中國人群的研究20篇[24, 26-27, 32-41, 46-47, 52-56], 國外人群的研究為9篇[25, 42-45, 48-51]。其中矽肺研究16項[25-27, 32, 34, 38-41, 47-51, 54, 56]，煤工塵肺研究14項[27，32, 34-37, 39, 42-44, 46, 52-53, 55]，煤矽肺研究1項[45]，塵肺未分類研究2項[24, 33]。文獻篩選流程情況見圖1, 納入文獻基本情況見表1。同一個基因位點有2項及以上研究數據時納入meta分析，僅有單項數據報告時不能進行meta分析。納入meta分析文獻中各基因位點的基因型分布見表2。其他有29種IL基因僅有單篇文獻研究，各基因型分布見表1。

圖1 文獻篩選流程圖

表1 納入研究基本情況表

2.2文獻質量評價結果 29篇納入文獻[24-27, 32-56]未進行計算把握度、連鎖不平衡檢測、計算錯誤率、多重檢驗校正、重復性以及功能驗證性試驗；29篇文獻均未描述是否采用盲法，所有病例組和對照組均來自同一人群，同時描述了病例組人口學及臨床信息。9篇文獻[25, 27, 32-34, 38-39, 43, 56]未提及HWE檢驗情況，其余文獻均符合HW平衡檢測。18篇文獻[25-26, 32, 34-38, 40-41, 43, 45, 47, 52-56]進行了多態性的鑒定，18篇文獻[25, 40-56]描述了人群分層分析。見圖2。

圖2 29篇納入文獻的偏倚風險評價

2.3Meta分析結果 本研究對17種IL相關基因分別采用等位基因模型、顯性模型、隱性模型、共顯性模型、超顯性模型等進行meta分析探討其與塵肺易感性的相關性。

2.3.1IL-1β(-511)C/T多態位點與塵肺易感性的關聯 29篇納入文獻[24-27，32-56]中有6篇文獻[27, 32, 36, 41-42, 45]包括8組數據提供了IL-1β(-511)C/T多態位點的各基因型和等位基因，包括病例組949例，對照組1 051例。對各基因模型數據進行異質性檢驗，I2>50%時說明存在明顯異質性，采用隨機效應模型，I2<50%，無明顯異質性，采用固定效應模型。Meta分析結果見表3。對各基因模型進一步行敏感性分析，隱性模型(CC+CT vs TT)和共顯性模型(TT vs CC)分別剔除Volobaev等[45]等的研究后，異質性檢驗I2分別降低為50.3%、62.8%, meta分析結果發生明顯改變(OR=1.36, 95%CI=0.97～1.92,P=0.079)、(OR=1.43, 95%CI=0.89～2.31,P=0.137)，剔除其他研究meta分析結果明顯改變。而其他模型，剔除任意一項研究對meta分析結果無明顯改變。對各基因模型進一步回歸分析顯示，塵肺種類、種族P值均<0.05，可能是異質性來源。

表2 納入meta分析文獻的IL基因多態位點的基因型分布匯總

2.3.2IL-1RA(+2018)T/C多態位點與塵肺易感性的關聯 5篇文獻[26-27，40, 43-44]包括6組數據提供了IL-1RA(+2018)T/C多態位點的各基因型和等位基因，包括病例組660例，對照組711例。對各基因模型數據進行異質性檢驗，I2>50%時說明存在明顯異質性，采用隨機效應模型，I2<50%，無明顯異質性，采用固定效應模型。Meta分析結果見表4。對各基因模型進一步行敏感性分析，剔除任意一項研究對meta分析結果無明顯改變。進一步回歸分析顯示，除共顯性模型、顯性模型外，其余模型顯示塵肺種類、種族P值均<0.05，提示塵肺種類、種族可能為異質性來源。

表3 IL-1β(-511)C/T多態性與塵肺易感性的meta分析和亞組分析結果

表4 IL-1RA(+2018)T/C多態性與塵肺易感性的meta分析和亞組分析結果

2.3.3IL-1α(-889)C/T多態位點與塵肺易感性的關聯 3篇文獻[27, 41, 51]包括4組數據提供了IL-1α(-889)C/T多態位點的各基因型和等位基因，包括病例組400例，對照組371例。對各基因模型數據進行異質性檢驗，I2>50%時說明存在明顯異質性，采用隨機效應模型，I2<50%無明顯異質性，采用固定效應模型。Meta分析結果見表5。對異質性較高的基因模型進行回歸分析，塵肺種類、種族P值均>0.05，未找到異質性來源。進一步行敏感性分析，剔除任意一項研究對meta分析結果無明顯改變。

2.3.4IL-1α(+4845)G/T多態位點與塵肺易感性的關聯 3篇文獻[42-44]包括3組數據提供了IL-1α(+4845)G/T多態位點的各基因型和等位基因，包括病例組696例，對照組600例。對各基因模型數據進行異質性檢驗，I2>50%時說明存在明顯異質性，采用隨機效應模型，I2<50%，無明顯異質性，采用固定效應模型。Meta分析結果見表2。

2.3.5IL-1β(+3953)C/T多態位點與塵肺易感性的關聯 4篇文獻[41-44]包括4組數據提供了IL-1β(+3953)C/T多態位點的各基因型和等位基因，包括病例組879例，對照組 711例。對各基因模型數據進行異質性檢驗，I2>50%時說明存在明顯異質性，采用隨機效應模型，I2<50%，無明顯異質性，采用固定效應模型。Meta分析結果見表6。對于異質性較高的基因模型進一步回歸分析，塵肺種類、種族P值均>0.05，未找到異質性來源。進一步敏感性分析，剔除任意一項研究對meta分析結果無明顯改變。

2.3.6IL-6(-174)G/C多態位點與塵肺易感性的關聯 5篇文獻[33, 42, 44-46]包括5組數據提供了IL-6(-174)G/C多態位點的各基因型和等位基因，包括病例組719例，對照組880例。對各基因模型數據進行異質性檢驗，I2>50%時說明存在明顯異質性，采用隨機效應模型，I2<50%，無明顯異質性，采用固定效應模型。Meta分析結果見表7。對異質性較高的基因模型進一步回歸分析顯示，塵肺種類、種族P值均<0.05，未找到異質性來源。進一步敏感性分析，共顯性模型(GC vs GG) 分析中剔除Volobaev等[45]的研究后，meta分析結果發生明顯改變(OR=0.72, 95%CI=0.52～0.98,P=0.079)，而逐一剔除任意一項研究，其他基因模型的meta分析結果無明顯改變。

表5 IL-1α(-889)C/T多態性與塵肺易感性的meta分析和亞組分析結果

表6 IL-1β(+3953)C/T多態性與塵肺易感性的meta分析和亞組分析結果

表7 IL-6(-174)G/C多態性與塵肺易感性的meta分析和亞組分析結果

2.3.7IL-6(-634)C/G多態位點與塵肺易感性的關聯 有3篇文獻[33-35]包括4組數據提供了IL-6(-634)C/G多態位點的各基因型和等位基因，包括病例組344例，對照組362例。對各基因模型數據進行異質性檢驗，I2>50%時說明存在明顯異質性，采用隨機效應模型，I2<50%，無明顯異質性，采用固定效應模型。Meta分析結果見表3。

2.3.8其他基因多態位點與塵肺易感性的關聯 1篇文獻[27]包括2組數據提供了IL-1α VNTR多態位點的各基因型和等位基因，包括病例組125例，對照組125例。1篇文獻[39]包括2組數據提供了IL-8(781)C/T、IL-8 Met31Arg T/G、IL-8RA(+860)G/C、IL-8(-251)A/T多態位點的各基因型和等位基因，包括病例組213例，對照組251例。2篇文獻[32, 47]包括3組數據提供了IL-10(-592)A/C多態位點的各基因型和等位基因，包括病例組271例，對照組199例。2篇文獻[37, 45]包括2組數據提供了IL-12β(+1188)A/C多態位點的各基因型和等位基因，包括病例組209例，對照組232例。2篇文獻[25, 54]包括2組數據提供了IL-17A(-832)A/G和IL-17F(+7488)A/G多態位點的各基因型和等位基因，包括病例組161例，對照組178例。2篇文獻[24, 55]包括2組數據提供了IL-17A rs4711998 A/G多態位點的各基因型和等位基因，包括病例組916例，對照組936例。對各基因模型數據進行異質性檢驗，I2>50%時說明存在明顯異質性，采用隨機效應模型，I2<50%，無明顯異質性，采用固定效應模型。Meta分析結果見表4。

2.3.9發表偏倚 對納入研究≥4項的各基因位點不同基因模型進行發表偏倚檢測，IL-1β(-511)C/T隱性模型漏斗圖不對稱，顯示可能存在發表偏倚, Egger法檢驗進一步分析顯示，t=3.33，P=0.016，均顯示發表偏倚可能性大，對其進行剪補法結果顯示，剪補前后合并效應值和P值無明顯變化，說明此meta分析結論相對穩定；其余研究漏斗圖大致對稱，Egger檢驗P值大于0.05，顯示發表偏倚可能性小。見圖3。

圖3 IL-1β(-511)C/T隱性模型與塵肺易感關聯的發表偏倚 a.漏斗圖；b. Egger圖；c. 剪補法圖

3 討 論

塵肺病是一種肺纖維化疾病，按照《職業病分類和目錄》，主要包括矽肺、煤工塵肺、石墨塵肺、炭黑塵肺、石棉肺、滑石塵肺，水泥塵肺、云母塵肺、陶工塵肺、鋁塵肺、電焊工塵肺、鑄工塵肺等十二種[5]。其中矽肺和煤工塵肺是我國最常見的塵肺種類。塵肺病發生大多需要10～20年甚至更長的接塵工齡，且脫離粉塵接觸后仍然可以易感，一旦易感，將難以治愈，給個人和社會造成嚴重的經濟負擔[5, 57]。塵肺病是可防可控的疾病[5]。塵肺的發生發展及病變程度主要受工人接觸粉塵的種類、濃度、分散度、工齡、個體因素、以及防護措施的影響[58]。雖然職業粉塵暴露是塵肺易感的首要因素，但塵肺易感并不完全取決于粉塵暴露，據報道，遺傳在塵肺的病因和易感機制中起著重要作用[59-61]。近年來，基因多態性與塵肺病易感風險的相關研究越來越多，包括IL[24-27，32-56]、腫瘤壞死因子-α(tumor necrosis factor-α，TNF-α)[62]、轉化生長因子-β(transforming growth factor box1，TGF-β)[63]、高遷移率族蛋白 B1 (high mobility group box-1 protein，HMGB-1)[64]、基質金屬蛋白酶(matrix metalloproteinases，MMP-2)[65]、熱休克蛋白70(heat shock protein，HSP70)[66]、免疫相關鳥苷三磷酸酶 M(immunity related GTPase M，IRGM)[67]等。IL屬于細胞因子，是一種具有重要的免疫調節的炎癥介質，在參與肺纖維化中扮演重要角色。IL-1在塵肺發生過程中主要刺激成纖維細胞，促使其分泌膠原及膠原酶，并且可以調節免疫介導的肺損傷[68]。矽塵導致矽肺的易感進程中肺泡巨噬細胞、T細胞等通過IL-4Ra通道發揮免疫功能[69]，IL-4在塵肺易感中起重要作用。IL-8 是一種白細胞趨化因子，通過趨化作用使炎癥細胞向肺組織遷移并釋放活性物質參與肺纖維化，可作為單獨的細胞因子對肺纖維化發揮作用[70]。目前研究與塵肺易感相關的IL基因主要有：IL-1、IL-4、IL-6、IL-8、IL-10、IL-12、IL-13、IL-17、IL-18等[24-27，32-56，68]。王煒等[28]曾對IL-1基因多態位點與塵肺病易感性進行系統評價和meta分析，但該研究僅對一種基因模型進行分析，不夠全面。本文通過全面搜索目前已發表的IL相關基因與塵肺病易感風險的相關研究，以多種遺傳基因模型進行進行綜合分析，結果可能更加全面可靠。

本文納入29篇關于IL相關基因多態位點與塵肺易感的病例對照設計研究，包括33組數據，共計病例組5 315例，對照組5 332例。中國人群的研究20篇[24，26-27，32-41, 46-47，52-56], 國外人群的研究為9篇[25，42-45, 48-51]。除了Salum等[51]研究為混合人種以外，其他國外研究的人群均為高加索人種。王煒等[28]曾對IL-1相關基因進行了meta分析研究結果顯示，IL-1β(-511)C/T、IL-1α(+4845)G/T等多態位點與矽肺及煤工塵肺的易感無明顯關聯。但本研究meta分析發現，IL-1β(-511)C/T、IL-1α(+4845)G/T等多態位點可能與塵肺病的易感風險有關。本研究納入的6篇[27, 32, 36, 41-42, 45]有關IL-1β(-511)C/T多態位點與塵肺病易感的研究meta分析顯示，隱性模型(CC+CT vs TT)、共顯性模型(TT vs CC)均提示IL-1β(-511)C/T多態位點可能增加塵肺病的發生風險，這與范雪云等[27]、Volobaev等[45]的研究結果一致，與王煒等[28]及大多數研究[32, 36, 41-42]結果不一致。結合塵肺種類的亞組分析結果，剔除Volobaev等[45]的研究后，meta分析結果發生明顯改變，提示Volobaev等[45]研究的煤矽肺種類可能為異質性來源。因煤矽肺研究只有1篇，樣本量有限，對本研究產生了局限性。而等位基因模型(T vs C)、顯性模型(TT+CT vs CC)、共顯性模型(CT vs CC)、超顯性模型(CC+TT vs CT)等遺傳模式分析結果均提示IL-1β(-511)C/T多態位點可能與塵肺病的易感性無關聯，這與大多數研究結果一致[28，32, 36, 41-42]。本研究納入的3項有關IL-1α(+4845)G/T多態位點與塵肺病易感的等位基因模型(T vs G)分析結果提示等位基因T可能會增加塵肺病的易感風險，雖其異質性高，但3項研究均為高加索人種煤工塵肺，且研究較少，無法進行亞組分析[42-44]。由于研究數和樣本量少，對研究結果的可靠性有待未來更多相關研究進一步證實。因此，結合多種遺傳模型的分析結果，提示IL-1β(-511)C/T、IL-1α(+4845)G/T等多態位點與塵肺病的易感風險可能有關。本研究結果顯示IL-1RA(+2018)T/C多態位點中，帶TC或CC基因型的個體可增加煤工塵肺或矽肺的易感風險，同時等位基因模型(C vs T)、隱性模型(TT+TC vs CC)、共顯性模型(CC vs TT)等的陽性結果更加證實了IL-1RA(+2018)T/C多態位點與塵肺病的易感風險相關，這與大多數研究結果一致[26，28，40, 43，71]。此外，高加索人種IL-1RA(+2018)T/C多態位點與塵肺易感風險無關，但是該位點可增加矽肺或中國人群的塵肺病易感風險，這提示塵肺種類和種族可能是臨床異質性來源。IL-1α(-889)C/T多態位點的顯性模型(TT+TC vs CC)、共顯性模型(TC vs CC) 均提示攜帶TC基因型可能增加塵肺病的易感風險，這與范雪云等[27]研究結果一致，但與Wu等[41]、Salum等[51]研究結果不一致。本研究亞組分析表明攜帶TC基因型能增加煤工塵肺的易感風險，攜帶TC基因型不能增加矽肺的易感風險，提示塵肺種類可能是其臨床異質性來源。

IL-6是機體重要的免疫-神經-內分泌調節因子，既可以通過誘導膠原合成又可以抑制膠原降解從而影響塵肺發生發展[68]。本研究中IL-6(-174)G/C多態位點各遺傳基因模型與塵肺病的易感風險無關，這與大多數研究結果一致[29，33, 42, 44-46]。相比于郝艷星等[29]的研究結果，本研究顯示除了等位基因模型、顯性模型和超顯性模型提示IL-6(-634)C/G多態位點可能與塵肺病易感有關外，共顯性模型(GG vs CC) 、共顯性模型(CG vs CC)也提示CC基因型可能增加塵肺的易感風險，等位基因C可能為塵肺等位基因的風險因素，而等位基因G可能為塵肺病的保護因素。Chen 等[54]研究發現，IL-17F 7488A/G與矽肺易感性有關，G等位基因可能具有保護作用。但本研究分析顯示IL-17F(+7488)A/G多態性與矽肺的易感風險可能無關，這與Hassani 等[25]研究一致。IL-8(781)C/T、IL-8 Met31Arg T/G多態性可能增加塵肺病的易感風險，IL-8(-251)A/T多態性可能降低塵肺病的易感風險。IL-1α(VNTR)、IL-8RA(+860)G/C、IL-10(-592)A/C、IL-12β(+1188)A/C、IL-17A(-832)A/G、IL-17A rs4711998 A/G、IL-17F(+7488)A/G等IL相關基因多態性的研究均未發現與塵肺病的易感風險有關。而其他如Wang 等[52]研究發現IL-4C -590 T的功能性多態性與中國黃種人群煤工塵肺的易感風險降低相關；Han 等[55]研究發現，IL-17 rs3748067 G>A和rs8193036 C>T多態性降低了煤工塵肺易感風險等。由于這些基因多數僅有單篇文獻報告，缺乏多中心數據，證據等級相對不足，期待未來有更多相關研究進一步探討分析。

本文部分納入研究的異質性較高，未能找到異質性來源。遺傳基因關聯性研究異質性來源可能是種族、民族、對照組來源、基因分型檢測方法，疾病表型、樣本量等等。本文僅對塵肺種類、種族進行亞組分析，由于納入研究較少，總體樣本量相對不多等種種局限性，未能進一步對民族、基因分型檢測方法、樣本量等進一步亞組分析。因此，本meta分析結果的異質性來源可能與臨床多種因素相關，期待未來有更多相關研究、更大樣本量，以及更好的研究設計，以進一步探討IL相關基因與塵肺易感關聯性。

IL-1RA(+2018)T/C多態位點與塵肺的易感風險顯著相關，等位基因C可能是塵肺發病的風險基因，攜帶CC或CT基因型可能會增加塵肺的發病風險，對矽肺的發病風險可能更為顯著。IL-6(-634)C/G多態位點與塵肺的易感風險有關，等位基因C可能是塵肺發病的風險基因，突變等位基因G可能為塵肺的保護因素，攜帶GG或GC基因型可能會減少塵肺的發病風險。IL-1β(+3953)C/T、IL-6(-174)G/C多態位點與塵肺病的易感風險無關聯。IL-1α(+4845)G/T多態位點與塵肺病的易感風險可能有關，等位基因T可能是塵肺發病的風險基因；IL-1β(-511)C/T多態位點與煤矽肺的易感風險可能有關；IL-1α(-889)C/T 多態位點與煤工塵肺的易感風險可能有關；IL-8(781)C/T、 IL-8 Met31Arg T/G、IL-8(-251)A/T多態性與塵肺的易感風險可能有關，需要更多相關研究進一步證實，期待后續多中心大樣本量、以及更好的研究設計進一步探討分析遺傳基因多態位點與塵肺易感性關聯。