引用本文: 石明巧, 陳捷, 王蘭蘭, 廖云, 張琦, 嚴琳, 劉新樂. TM6SF2基因rs58542926多態性與肝細胞損傷和肝纖維化嚴重程度相關性的Meta分析. 中國循證醫學雜志, 2016, 16(6): 655-662. doi: 10.7507/1672-2531.20160102 復制
非酒精性脂肪性肝病(non-alcoholic fatty liver disease,NAFLD)是一組進行性肝病的統稱,包括單純脂肪變性、非酒精性脂肪性肝炎、肝纖維化等,其以彌漫性肝細胞大泡性脂肪變性為主要特征,同時可以除外酒精和其他明確的肝損因素,患者肝中甘油三脂的含量增多,肝細胞脂肪變性的病理改變逐漸進展至纖維化甚至硬化。在肝硬化前期,往往沒有明顯癥狀出現,只能通過侵入性方法輔助診斷,而肝纖維化是肝細胞癌(hepatocellular carcinoma,HCC)衍變的關鍵環節,故有效控制此環節對預防肝硬化及HCC發生具有極為重要的意義 [1]。
近年來隨著精準醫療的發展,醫學實驗室越來越關注分子標志在疾病輔助診斷和預后判斷中的應用,在肝纖維化和硬化等肝細胞損傷中,最新的研究發現跨膜蛋白6超家族成員2(TM6SF2)的變異與NAFDL的發生存在相關性,TM6SF2基因的變異會使甘油三酯沉積在肝臟中而使肝臟中脂質含量增多,導致肝臟發生脂肪變性引起NAFLD,并逐漸發展為嚴重肝纖維化或者肝硬化 [2, 3]。但是,不同的研究對于TM6SF2基因的變異是否與脂質的代謝、肝纖維化嚴重程度以及肝細胞的損傷存在相關性還有爭議,本研究采用Meta分析的方法評價TM6SF2基因多態性對脂代謝指標以及肝纖維化嚴重程度的影響,為進一步探索新的分子標志物提供理論基礎。
1 資料與方法
1.1 納入與排除標準
1.1.1 研究類型
TM6SF2基因多態性與肝損傷及肝纖維化嚴重程度關系的橫斷面研究。
1.1.2 暴露因素 rs58542926位點多態性。
1.1.3 結局指標
嚴重肝纖維化的風險、AST、ALT水平。
1.1.4 排除標準
① 重復發表的文獻;② 沒有相關基因型的具體數據的文獻;③ 肝纖維化不是通過活檢確診;④ 非中、英文文獻。
1.2 檢索策略
計算機檢索PubMed、CNKI、WanFang Data和CBM數據庫,搜集有關TM6SF2基因多態性與肝纖維化以及肝細胞損傷的相關研究,檢索時限為建庫至2016年1月27日。同時手檢納入研究的參考文獻,以補充獲取相關文獻。中文檢索詞為:TM6SF2、rs58542926;英文檢索詞為:TM6SF2、rs58542926。以PubMed為例,其具體檢索策略見框 1。
框 1 PubMed 檢索策略
TM6SF2 rs58542926 #1 OR #2
1.3 文獻篩選、資料提取和方法學質量評價
由2位研究者獨立篩選文獻、提取資料,并交叉核對。如遇分歧,通過咨詢第三方協助解決。文獻篩選時首先閱讀文題和摘要,在排除明顯不相關的文獻后,進一步閱讀全文來確定是否納入。采用自制的資料提取表提取資料,提取內容主要包括:① 納入研究的基本特征,包括第一作者、發表時間、國家;② 納入患者的基本特征,包括診斷標準、種族、年齡等;③ 各種基因型對應的樣本量、基因型檢測方法、是否符合HWE平衡及其他指標。由于目前對于橫斷面研究的偏倚風險的評價沒有一個明確的標準,本研究根據AHRQ量表及相關文獻評價納入研究的方法學質量,共6個條目,分別作出“是”、“否”和“不清楚”的評價,評價為“是”得 1分,“否”或者“不清楚”得0分(表 1)。

1.4 統計分析
采用Stata 12.0軟件進行Meta分析。計數資料采用比值比(OR)為效應指標,計量資料采用均數差(MD)為效應指標,各效應量均給出其點估計值和95%CI。納入研究結果間的異質性采用χ2檢驗進行分析(檢驗水準為α=0.1),同時結合I2定量判斷異質性的大小。若各研究結果間無統計學異質性,采用固定效應模型進行Meta分析;若各研究結果間存在統計學異質性,則進一步分析異質性來源,在排除明顯臨床異質性的影響后,采用隨機效應模型進行 Meta 分析。明顯的臨床異質性采用亞組分析或敏感性分析等方法進行處理。采用Begg’s檢驗評估納入研究的發表偏倚,檢驗水準設為α=0.05。
2 結果
2.1 文獻篩選流程及結果
初檢出相關文獻49篇,其中中文1篇,英文48篇,經逐層篩選后,最終納入12篇文獻 [3-14],共23個研究。文獻篩選流程及結果見圖 1。

2.2 納入研究的基本特征
納入研究的基本特征和方法學質量評價結果見表 2。

2.3 Meta分析結果
2.3.1 TM6SF2-rs58542926與脂代謝指標的相關性
2.3.1.1 TM6SF2-rs58542926與TC水平的相關性
共納入16個研究 [4-12, 14],Meta分析結果顯示:TM6SF2-rs58542926基因位點多態性與TC水平相關,且差異有統計學意義[MD= -9.628,95%CI(-12.019,-7.237),P < 0.000?1](圖 2)。

2.3.1.2 TM6SF2-rs58542926與TG水平的相關性
共納入18個研究[4-13],Meta分析結果顯示:TM6SF2-rs58542926基因位點多態性與TG水平相關[MD= -11.666,95%CI(-16.088,-7.244),P < 0.000?1](圖 3)。基于人種的亞組分析結果顯示: 在高加索人群中,TM6SF2-rs58542926基因位點多態性與TG水平相關[MD= -12.820,95%CI(-17.847,-7.793),P < 0.000?1];而在亞洲及其他人群中,該基因位點多態性與TG水平不具有相關性[MD= -5.698,95%CI(-13.744,2.348),P=0.165]。基于年齡的亞組分析結果提示:在成年人與兒童及青少年中,TM6SF2-rs58542926基因位點多態性與TG水平均相關[MD= -13.139,95%CI(-18.445,-7.833),P < 0.000?1;MD= -7.844,95%CI(-15.652,-0.035),P=0.049](表 3)。


2.3.1.3 TM6SF2-rs58542926與LDL-C水平的相關性
共納入14個研究 [5, 8-11, 13, 14] Meta分析結果顯示:TM6SF2-rs58542926基因位點多態性與LDL-C水平相關[MD= -3.821,95%CI(-5.912,-1.730),P < 0.000?1],(圖 4)。亞組分析結果顯示:在不同人群和不同年齡中,TM6SF2-rs58542926基因位點多態性與LDL-C水平均相關(P均 < 0.05)(表 3)。

2.3.2 TM6SF2-rs58542926與肝細胞損傷的相關性
2.3.2.1 TM6SF2-rs58542926與ALT水平的相關性
共納入9個研究 [4, 6, 10, 11, 13],Meta分析結果顯示:TM6SF2-rs58542926基因位點多態性與ALT水平相關[MD=1.185,95%CI(0.042,2.329),P=0.042](圖 5)。基于人種的亞組分析結果顯示:在高加索人群中,TM6SF2-rs58542926基因位點多態性與ALT水平相關[MD=1.298,95%CI(0.208,2.387),P=0.020];而在亞洲及其他人群中,該基因位點多態性與ALT水平不具有相關性[MD= -0.054,95%CI(-4.927,4.819),P=0.983]。基于年齡的亞組分析結果顯示:在成年人中,TM6SF2-rs58542926基因位點多態性與ALT水平相關[MD=1.083,95%CI(0.790,1.377),P < 0.000?1];而在兒童及青少年組中,該基因位點多態性與ALT水平不具有相關性[MD=4.397,95%CI(-3.361,12.156),P=0.267](表 3)。

2.3.2.2 TM6SF2-rs58542926與AST水平的相關性
共納入9個研究 [4, 6, 10, 11, 13],Meta分析結果顯示:TM6SF2-rs58542926基因位點多態性與AST水平存在相關性[MD=0.860,95%CI(0.560,1.161),P < 0.000?1](圖 6)。

2.3.3 TM6SF2-rs58542926與肝纖維化嚴重程度的相關性
共納入4個研究 [3, 5, 8],Meta分析結果顯示:TM6SF2-rs58542926基因位點多態性與肝纖維化嚴重程度相關[OR=1.491,95%CI(1.131,1.851),P < 0.000?1](圖 7)。

2.4 敏感性分析
逐一剔除納入研究進行敏感性分析,Meta分析結果的方向并未發生改變,僅在TM6SF2-rs58542926基因位點多態性與ALT水平相關性的分析中,在剔除了Goffredo等 [10]的研究后,研究間的異質性減少,I2為43.3%。
2.5 發表偏倚
各組研究在漏斗圖兩側的分布對稱性較好,提示無明顯發表偏倚(圖 8~10),Begg's檢驗也提示本研究無發表偏倚(P > 0.05)。



3 討論
全基因組關聯研究(GWAS)發現TM6SF2在分泌富含TG的脂蛋白以及富含TG的脂滴中起重要作用,參與介導脂質在肝臟的沉積,而大量的脂質沉積在肝中會引起肝臟損傷,導致NAFLD,并逐步進展為肝纖維化,同時也有研究報道該基因介導的脂質沉積與心血管疾病的風險密切相關。隨著精準醫療的發展,醫學實驗室越來越關注分子標志在疾病輔助診斷和預后判斷中的應用。對于肝纖維化,目前只能通過有創的肝活檢來確定,尋找新的輔助檢測方法來預防或明確肝纖維化顯得尤為重要。
在本研究中,我們采用Meta分析方法,對TM6SF2基因的錯意突變與脂質代謝、肝細胞損傷及肝纖維化的相關性進行了評價。結果顯示,TM6SF2-rs58542926的變異與脂質的形成及代謝密切相關,并且與肝纖維化進展以及肝細胞損傷也存在相關性。這與2015年Akuta等 [15]對日本人群的研究和Arslanow等 [16]對174例患者的研究結論不一致,這可能是由于監測肝纖維化的方法不同或結局指標不一樣。本研究還發現T的攜帶者(EK+KK的個體)發生嚴重肝纖維化的風險大于C的純合子(EE),然而他們的脂代謝指標如TC、TG、LDL-C的水平卻較低,分別為-9.628、-11.67和-3.821 mg/dL。采用隨機效應模型的分析發現TG、LDL-C與TM6SF2變異具有相關性。該結果表明TM6SF2突變型的個體血清中脂質的含量低于野生型個體,這可能是由于突變型個體大量的脂質沉積于肝臟中,而不能從肝臟中分泌出來,大量脂質沉積于肝臟中必然會對肝臟造成一些損傷,進而發展為肝纖維化。此外,本研究還根據種族以及年齡進行了亞組分析,結果提示,在高加索人群中,TM6SF2-rs58542926基因位點多態性與TG和LDL-C水平相關,亞洲人群僅與LDL-C水平相關。在年齡結構方面,兒童及青少年除在ALT水平方面與該多態性不具有相關性,其他指標在年齡上均未體現差別。
本研究的局限性:① 部分納入研究樣本量較小,因而合并分析的檢驗效能仍可能不足;② 部分研究因缺乏活檢結果的支持,或缺乏具體的數據,而未被納入;③ 納入研究人群高加索人群占絕大多數,亞洲人群較少,從而導致本研究結果的外推性會受到一定程度的影響;④ 通過年齡和種族分層分析,研究的異質性依然存在,提示研究結果可能受其他因素影響,如個體生活環境的不同。
綜上所述,TM6SF2的基因變異與不同人群脂質形成及代謝密切相關,突變型個體血脂含量水平較低,發生嚴重肝纖維化及肝損傷的風險更大。受納入研究數量和質量的限制,上述結論尚需開展更多高質量的研究予以驗證。
非酒精性脂肪性肝病(non-alcoholic fatty liver disease,NAFLD)是一組進行性肝病的統稱,包括單純脂肪變性、非酒精性脂肪性肝炎、肝纖維化等,其以彌漫性肝細胞大泡性脂肪變性為主要特征,同時可以除外酒精和其他明確的肝損因素,患者肝中甘油三脂的含量增多,肝細胞脂肪變性的病理改變逐漸進展至纖維化甚至硬化。在肝硬化前期,往往沒有明顯癥狀出現,只能通過侵入性方法輔助診斷,而肝纖維化是肝細胞癌(hepatocellular carcinoma,HCC)衍變的關鍵環節,故有效控制此環節對預防肝硬化及HCC發生具有極為重要的意義 [1]。
近年來隨著精準醫療的發展,醫學實驗室越來越關注分子標志在疾病輔助診斷和預后判斷中的應用,在肝纖維化和硬化等肝細胞損傷中,最新的研究發現跨膜蛋白6超家族成員2(TM6SF2)的變異與NAFDL的發生存在相關性,TM6SF2基因的變異會使甘油三酯沉積在肝臟中而使肝臟中脂質含量增多,導致肝臟發生脂肪變性引起NAFLD,并逐漸發展為嚴重肝纖維化或者肝硬化 [2, 3]。但是,不同的研究對于TM6SF2基因的變異是否與脂質的代謝、肝纖維化嚴重程度以及肝細胞的損傷存在相關性還有爭議,本研究采用Meta分析的方法評價TM6SF2基因多態性對脂代謝指標以及肝纖維化嚴重程度的影響,為進一步探索新的分子標志物提供理論基礎。
1 資料與方法
1.1 納入與排除標準
1.1.1 研究類型
TM6SF2基因多態性與肝損傷及肝纖維化嚴重程度關系的橫斷面研究。
1.1.2 暴露因素 rs58542926位點多態性。
1.1.3 結局指標
嚴重肝纖維化的風險、AST、ALT水平。
1.1.4 排除標準
① 重復發表的文獻;② 沒有相關基因型的具體數據的文獻;③ 肝纖維化不是通過活檢確診;④ 非中、英文文獻。
1.2 檢索策略
計算機檢索PubMed、CNKI、WanFang Data和CBM數據庫,搜集有關TM6SF2基因多態性與肝纖維化以及肝細胞損傷的相關研究,檢索時限為建庫至2016年1月27日。同時手檢納入研究的參考文獻,以補充獲取相關文獻。中文檢索詞為:TM6SF2、rs58542926;英文檢索詞為:TM6SF2、rs58542926。以PubMed為例,其具體檢索策略見框 1。
框 1 PubMed 檢索策略
TM6SF2 rs58542926 #1 OR #2
1.3 文獻篩選、資料提取和方法學質量評價
由2位研究者獨立篩選文獻、提取資料,并交叉核對。如遇分歧,通過咨詢第三方協助解決。文獻篩選時首先閱讀文題和摘要,在排除明顯不相關的文獻后,進一步閱讀全文來確定是否納入。采用自制的資料提取表提取資料,提取內容主要包括:① 納入研究的基本特征,包括第一作者、發表時間、國家;② 納入患者的基本特征,包括診斷標準、種族、年齡等;③ 各種基因型對應的樣本量、基因型檢測方法、是否符合HWE平衡及其他指標。由于目前對于橫斷面研究的偏倚風險的評價沒有一個明確的標準,本研究根據AHRQ量表及相關文獻評價納入研究的方法學質量,共6個條目,分別作出“是”、“否”和“不清楚”的評價,評價為“是”得 1分,“否”或者“不清楚”得0分(表 1)。

1.4 統計分析
采用Stata 12.0軟件進行Meta分析。計數資料采用比值比(OR)為效應指標,計量資料采用均數差(MD)為效應指標,各效應量均給出其點估計值和95%CI。納入研究結果間的異質性采用χ2檢驗進行分析(檢驗水準為α=0.1),同時結合I2定量判斷異質性的大小。若各研究結果間無統計學異質性,采用固定效應模型進行Meta分析;若各研究結果間存在統計學異質性,則進一步分析異質性來源,在排除明顯臨床異質性的影響后,采用隨機效應模型進行 Meta 分析。明顯的臨床異質性采用亞組分析或敏感性分析等方法進行處理。采用Begg’s檢驗評估納入研究的發表偏倚,檢驗水準設為α=0.05。
2 結果
2.1 文獻篩選流程及結果
初檢出相關文獻49篇,其中中文1篇,英文48篇,經逐層篩選后,最終納入12篇文獻 [3-14],共23個研究。文獻篩選流程及結果見圖 1。

2.2 納入研究的基本特征
納入研究的基本特征和方法學質量評價結果見表 2。

2.3 Meta分析結果
2.3.1 TM6SF2-rs58542926與脂代謝指標的相關性
2.3.1.1 TM6SF2-rs58542926與TC水平的相關性
共納入16個研究 [4-12, 14],Meta分析結果顯示:TM6SF2-rs58542926基因位點多態性與TC水平相關,且差異有統計學意義[MD= -9.628,95%CI(-12.019,-7.237),P < 0.000?1](圖 2)。

2.3.1.2 TM6SF2-rs58542926與TG水平的相關性
共納入18個研究[4-13],Meta分析結果顯示:TM6SF2-rs58542926基因位點多態性與TG水平相關[MD= -11.666,95%CI(-16.088,-7.244),P < 0.000?1](圖 3)。基于人種的亞組分析結果顯示: 在高加索人群中,TM6SF2-rs58542926基因位點多態性與TG水平相關[MD= -12.820,95%CI(-17.847,-7.793),P < 0.000?1];而在亞洲及其他人群中,該基因位點多態性與TG水平不具有相關性[MD= -5.698,95%CI(-13.744,2.348),P=0.165]。基于年齡的亞組分析結果提示:在成年人與兒童及青少年中,TM6SF2-rs58542926基因位點多態性與TG水平均相關[MD= -13.139,95%CI(-18.445,-7.833),P < 0.000?1;MD= -7.844,95%CI(-15.652,-0.035),P=0.049](表 3)。


2.3.1.3 TM6SF2-rs58542926與LDL-C水平的相關性
共納入14個研究 [5, 8-11, 13, 14] Meta分析結果顯示:TM6SF2-rs58542926基因位點多態性與LDL-C水平相關[MD= -3.821,95%CI(-5.912,-1.730),P < 0.000?1],(圖 4)。亞組分析結果顯示:在不同人群和不同年齡中,TM6SF2-rs58542926基因位點多態性與LDL-C水平均相關(P均 < 0.05)(表 3)。

2.3.2 TM6SF2-rs58542926與肝細胞損傷的相關性
2.3.2.1 TM6SF2-rs58542926與ALT水平的相關性
共納入9個研究 [4, 6, 10, 11, 13],Meta分析結果顯示:TM6SF2-rs58542926基因位點多態性與ALT水平相關[MD=1.185,95%CI(0.042,2.329),P=0.042](圖 5)。基于人種的亞組分析結果顯示:在高加索人群中,TM6SF2-rs58542926基因位點多態性與ALT水平相關[MD=1.298,95%CI(0.208,2.387),P=0.020];而在亞洲及其他人群中,該基因位點多態性與ALT水平不具有相關性[MD= -0.054,95%CI(-4.927,4.819),P=0.983]。基于年齡的亞組分析結果顯示:在成年人中,TM6SF2-rs58542926基因位點多態性與ALT水平相關[MD=1.083,95%CI(0.790,1.377),P < 0.000?1];而在兒童及青少年組中,該基因位點多態性與ALT水平不具有相關性[MD=4.397,95%CI(-3.361,12.156),P=0.267](表 3)。

2.3.2.2 TM6SF2-rs58542926與AST水平的相關性
共納入9個研究 [4, 6, 10, 11, 13],Meta分析結果顯示:TM6SF2-rs58542926基因位點多態性與AST水平存在相關性[MD=0.860,95%CI(0.560,1.161),P < 0.000?1](圖 6)。

2.3.3 TM6SF2-rs58542926與肝纖維化嚴重程度的相關性
共納入4個研究 [3, 5, 8],Meta分析結果顯示:TM6SF2-rs58542926基因位點多態性與肝纖維化嚴重程度相關[OR=1.491,95%CI(1.131,1.851),P < 0.000?1](圖 7)。

2.4 敏感性分析
逐一剔除納入研究進行敏感性分析,Meta分析結果的方向并未發生改變,僅在TM6SF2-rs58542926基因位點多態性與ALT水平相關性的分析中,在剔除了Goffredo等 [10]的研究后,研究間的異質性減少,I2為43.3%。
2.5 發表偏倚
各組研究在漏斗圖兩側的分布對稱性較好,提示無明顯發表偏倚(圖 8~10),Begg's檢驗也提示本研究無發表偏倚(P > 0.05)。



3 討論
全基因組關聯研究(GWAS)發現TM6SF2在分泌富含TG的脂蛋白以及富含TG的脂滴中起重要作用,參與介導脂質在肝臟的沉積,而大量的脂質沉積在肝中會引起肝臟損傷,導致NAFLD,并逐步進展為肝纖維化,同時也有研究報道該基因介導的脂質沉積與心血管疾病的風險密切相關。隨著精準醫療的發展,醫學實驗室越來越關注分子標志在疾病輔助診斷和預后判斷中的應用。對于肝纖維化,目前只能通過有創的肝活檢來確定,尋找新的輔助檢測方法來預防或明確肝纖維化顯得尤為重要。
在本研究中,我們采用Meta分析方法,對TM6SF2基因的錯意突變與脂質代謝、肝細胞損傷及肝纖維化的相關性進行了評價。結果顯示,TM6SF2-rs58542926的變異與脂質的形成及代謝密切相關,并且與肝纖維化進展以及肝細胞損傷也存在相關性。這與2015年Akuta等 [15]對日本人群的研究和Arslanow等 [16]對174例患者的研究結論不一致,這可能是由于監測肝纖維化的方法不同或結局指標不一樣。本研究還發現T的攜帶者(EK+KK的個體)發生嚴重肝纖維化的風險大于C的純合子(EE),然而他們的脂代謝指標如TC、TG、LDL-C的水平卻較低,分別為-9.628、-11.67和-3.821 mg/dL。采用隨機效應模型的分析發現TG、LDL-C與TM6SF2變異具有相關性。該結果表明TM6SF2突變型的個體血清中脂質的含量低于野生型個體,這可能是由于突變型個體大量的脂質沉積于肝臟中,而不能從肝臟中分泌出來,大量脂質沉積于肝臟中必然會對肝臟造成一些損傷,進而發展為肝纖維化。此外,本研究還根據種族以及年齡進行了亞組分析,結果提示,在高加索人群中,TM6SF2-rs58542926基因位點多態性與TG和LDL-C水平相關,亞洲人群僅與LDL-C水平相關。在年齡結構方面,兒童及青少年除在ALT水平方面與該多態性不具有相關性,其他指標在年齡上均未體現差別。
本研究的局限性:① 部分納入研究樣本量較小,因而合并分析的檢驗效能仍可能不足;② 部分研究因缺乏活檢結果的支持,或缺乏具體的數據,而未被納入;③ 納入研究人群高加索人群占絕大多數,亞洲人群較少,從而導致本研究結果的外推性會受到一定程度的影響;④ 通過年齡和種族分層分析,研究的異質性依然存在,提示研究結果可能受其他因素影響,如個體生活環境的不同。
綜上所述,TM6SF2的基因變異與不同人群脂質形成及代謝密切相關,突變型個體血脂含量水平較低,發生嚴重肝纖維化及肝損傷的風險更大。受納入研究數量和質量的限制,上述結論尚需開展更多高質量的研究予以驗證。