汪洋，嚴茂林，陳和平

(1電子科技大學附屬醫(yī)院 四川省人民醫(yī)院，成都610072；2電子科技大學醫(yī)學院)

慢性胰腺炎(CP)是胰腺長期進行性不可逆性炎性疾病，持續(xù)性胰腺炎癥會導致罹患胰腺癌的風險增加。CP的病因復雜，基因與環(huán)境相互作用在CP發(fā)生發(fā)展中發(fā)揮重要作用。近年來，與CP相關(guān)的易感基因和保護基因變異位點不斷被發(fā)現(xiàn)，糜蛋白酶C(CTRC)基因已被證實為CP易感基因，該基因部分突變可以使酶分泌減少并降低其酶活性[1，2]，從而增加CP的發(fā)生風險。CTRC基因編碼區(qū)和非編碼區(qū)有大量新的單核苷酸多態(tài)性(SNP)位點，這些SNP位點與CP發(fā)病有關(guān)。本研究旨在通過基因測序篩選西南地區(qū)漢族人群CTRC基因有意義的基因變異，從基因?qū)用嫔咸接慍P的發(fā)病機制。

1 資料與方法

1.1 臨床資料 選取2014年10月～2018年2月在四川省人民醫(yī)院和崇州市人民醫(yī)院住院治療的漢族CP患者108例(CP組)，男62例、女46例，年齡(57.0±15.4)歲，吸煙35例、飲酒24例。均符合中華醫(yī)學會《慢性胰腺炎診治指南》[3]中的診斷標準。排除合并其他胰腺疾病(如胰腺癌、胰腺囊腫等)；高血壓、糖尿病等代謝性疾??；遺傳性疾病；自身免疫性疾病患者。另選取同期健康體檢的漢族人群150例作為對照組，男88例、女62例，年齡(60.0±16.9)歲，吸煙32例、飲酒28例。兩組性別、年齡、吸煙、飲酒情況具有可比性。本研究經(jīng)四川省人民醫(yī)院倫理委員會批準，受試者均簽署知情同意書。

1.2 CTRC基因檢測方法

1.2.1 DNA提取 采集受試者外周靜脈血4 mL，EDTA抗凝。采用DNA抽提試劑盒提取白細胞核DNA。利用NanoDrop 2000分光光度計檢測DNA濃度與純度，OD值檢測在1.6～1.8保證樣本達標，校正樣本濃度為50 ng/μL，-20 ℃保存。

1.2.2 引物設(shè)計及合成 在UCSC數(shù)據(jù)庫檢索CTRC基因序列，截取外顯子及附近內(nèi)含子序列，使用Primer3.0在線設(shè)計引物，由上海生工公司進行序列合成。稀釋保存于-20 ℃。

1.2.3 PCR擴增 PCR擴增體積10 μL。反應采用MyCycler Thermal Cycler擴增儀。反應體系：2×Taq Master Mix 5 μL，正向引物0.5 μL，反向引物0.5 μL，DNA 1 μL，雙蒸水3 μL。反應條件：95 ℃ 5 min；95 ℃ 30 s，55～65 ℃ 30 s，72 ℃ 45 s，共35個循環(huán)；72 ℃ 7 min，12 ℃保溫。取擴增產(chǎn)物，采用2%瓊脂糖電泳，電壓120 V、時間20 min，獲取清晰特異的條帶，取對應的PCR產(chǎn)物進行測序。

1.2.4 PCR產(chǎn)物純化及測序 使用純化酶將PCR產(chǎn)物純化后，用上述擴增儀進行純化程序和測序程序。純化體積10 μL。純化反應體系：純化酶2 μL，PCR產(chǎn)物6 μL。純化條件：80 ℃15 min，37 ℃ 30 min，12 ℃保溫。測序體積10 μL。測序反應體系：BigDye 0.3 μL，5×Buffer 2 μL，引物(正向或反向)0.33 μL，純化后的DNA 2 μL，雙蒸水5.37 μL。測序條件：95 ℃ 5 min；95 ℃ 30 s，55～65 ℃ 30 s，72 ℃ 45 s，共29個循環(huán)；72 ℃ 7 min，12 ℃保溫。反應后產(chǎn)物進行乙醇純化離心，避光晾干30 min，加入10 μL雙蒸水，使用ABI3130基因分析儀進行測序，Chromas Lite Application軟件分析實驗結(jié)果，判定外顯子序列堿基是否有變化。

1.2.5 序列比對及預測分析 測序峰圖中不同顏色代表不同堿基，綠色表示脫氧腺嘌呤核苷酸，黑色表示脫氧鳥嘌呤核苷酸，紅色表示脫氧胸腺嘧啶核苷酸，藍色表示脫氧胞嘧啶核苷酸。單峰表示基因型為純合，雙峰表示基因型出現(xiàn)雜合。使用DNAMAN軟件將測序結(jié)果與數(shù)據(jù)庫獲得的目的基因相應序列進行比對，確定判讀結(jié)果的可靠性。

1.3 統(tǒng)計學方法 采用SPSS17.0統(tǒng)計軟件。使用Fisher確切概率檢驗頻率變化。使用優(yōu)勢比(OR)和相應95%置信區(qū)間(95%CI)表示遺傳模型和環(huán)境與CP的相關(guān)性。顯性模型為AB+AAvs.BB(A為最小等位基因，B為主要基因)。P<0.05為差異有統(tǒng)計學意義。

2 結(jié)果

2.1 兩組CTRC基因變異檢測結(jié)果 CP組CTRC基因第6外顯子檢測到1例c.611G>A雜合錯義突變，導致第204位甘氨酸改變?yōu)楣劝彼?，預測突變性質(zhì)為可能致病；第7外顯子檢測到1例c.703G>T雜合錯義突變，導致第703位纈氨酸改變?yōu)楸奖彼?，預測突變性質(zhì)為可能致病；第7外顯子檢測到1例c.736C>T雜合錯義突變，導致第736位精氨酸改變?yōu)榘腚装彼?，預測突變性質(zhì)為良性。對照組外顯子區(qū)域未檢出CTRC基因突變。另外發(fā)現(xiàn)5個SNP位點，均為于內(nèi)含子區(qū)域，分別為第1內(nèi)含子的c.40+92T>A、c.40+133G>A，第5內(nèi)含子的c.493+49G>C、c.493+51C>A，第7內(nèi)含子的c.793-137A>G；CP組c.40+133G>A、c.493+49G>C位點發(fā)生頻率高于對照組(P<0.05或<0.01)。見表1。

表1 兩組CTRC基因內(nèi)含子區(qū)域SNP位點發(fā)生情況比較[例(%)]

2.2 兩組CTRC基因SNP位點基因型頻率與等位基因頻率比較 經(jīng)上述統(tǒng)計分析篩選有意義的SNP位點，進一步行基因型頻率和等位基因頻率比較。CP組與對照組c.40+133G>A位點基因型頻率和等位基因頻率比較、c.493+49G>C位點基因型頻率和等位基因頻率比較，差異均有統(tǒng)計學意義(P<0.05或<0.01)。見表2、3。

表2 兩組CTRC基因SNP位點基因型頻率比較[例(%)]

表3 兩組CTRC基因SNP位點等位基因頻率比較[例(%)]

2.3 CTRC基因SNP遺傳模型與CP的關(guān)系 在顯性模型中,c.40+133G>A位點在兩組間的分布比較，CP組患病風險是對照組的1.088倍(P=0.030，95%CI為1.002～1.181)；c.493+49G>C位點在兩組間的分布比較，CP組患病風險是對照組的1.150倍(P=0.004，95%CI為1.038～1.275)。

2.4 CTRC基因SNP位點與吸煙、飲酒的關(guān)系 兩組飲酒者比較，c.40+133G>A位點和c.493+49G>C位點的分布差異均無統(tǒng)計學意義(P均>0.05)；兩組吸煙者比較，c.493+49G>C位點分布差異有統(tǒng)計學意義(OR=4.667，95%CI=1.341～16.239，P=0.014)。見表4、5。

表4 兩組吸煙者的CTRC基因SNP位點分布比較(例)

表5 兩組飲酒者的CTRC基因SNP位點分布比較(例)

3 討論

CP是一種由遺傳和環(huán)境因素共同作用所致的進展性疾病，病因復雜。研究發(fā)現(xiàn)，CTRC基因的功能喪失性突變會增加CP患病風險[4～7]。Grabarczyk等[8]發(fā)現(xiàn)，SNP位點與CP有很強的相關(guān)性。CTRC基因位于染色體1p36.21上，長度8.2 kb，包含8個外顯子。該基因編碼CTRC蛋白，非活化狀態(tài)的CTRC蛋白由16個氨基酸的分泌信號肽、13個氨基酸的前肽(活化肽)和239個氨基酸的糜蛋白酶樣酶組成，活化后可以促進胰蛋白酶原的自動激活，降解胰蛋白酶和胰蛋白酶原，并輔助活化CPA1和CPA2[9]。與CP相關(guān)的CTRC基因突變目前共發(fā)現(xiàn)54個，包括26個錯義突變、5個移碼突變、4個同義突變、一個框內(nèi)微突變及18個非編碼區(qū)的突變[10]，如主要發(fā)生在印度人的p.A73T和p.V235I突變以及主要發(fā)生在歐洲人種的p.R254W和p.K247_R254del突變[11，12]等。突變譜的不一致表明疾病的遺傳學病因與地域、人種等有關(guān)。此外CTRC基因的SNP改變也是CP發(fā)生的一個重要關(guān)聯(lián)因素。目前共發(fā)現(xiàn)3 209個CTRC基因SNP改變，包括rs121909294等致病性改變。Masson等[13]研究發(fā)現(xiàn)，兩個常見的同義CTRC基因多態(tài)性c.285C>T和c.180C>T，且c.180C>T多態(tài)性位點與CP發(fā)生有相關(guān)性。目前CP中研究最多的SNP位點是c.180C>T。

本研究利用Sanger直接測序法，對CTRC基因進行位點變異分析。結(jié)果顯示，CP組第6外顯子檢出1例c.611G>A雜合錯義突變，外顯子7上分別檢出1例c.703G>T雜合錯義突變和1例c.736C>T雜合錯義突變，對照組未檢出上述突變。檢索國內(nèi)外文獻，c.611G>A和c.703G>T雜合錯義突變未見報道。軟件預測突變c.611G>A和c.703G>T為致病性突變，c.736C>T突變?yōu)榉侵虏⌒酝蛔?。生物信息學領(lǐng)域的預測尚需基礎(chǔ)的功能驗證上述突變是否為致病性突變以及其具體的致病機制。此外還檢測到5個位于內(nèi)含子的SNP位點改變：c.40+133G>A、c.40+92T>A、c.493+49G>C、c.493+51C>A和c.793-137A>G。位點c.40+133G>A為新發(fā)現(xiàn)的SNP改變，未在NCBI數(shù)據(jù)庫上找到相應的rs號。CTRC基因c.40+133G>A和c.493+49G>C位點具有組間顯著差異。c.40+133G>A位點基因型分布和等位基因頻率在CP組分布高于對照組，在顯性模型下CP患病風險是對照組的1.088倍。c.493+49G>C位點基因型分布和等位基因頻率在CP組也顯著高于對照組，在顯性模型下CP患病的風險是對照組的1.150倍。表明CTRC基因c.40+133G>A和c.493+49G>C位點與中國西南漢族CP患者的發(fā)病可能存在一定關(guān)系，提示基因與疾病關(guān)聯(lián)的異質(zhì)性與種族和地區(qū)等因素有關(guān)。CP與吸煙和飲酒有很強的關(guān)聯(lián)性。LaRusch等[14]研究表明，伴有CTRC基因G60G位點單體型的吸煙者或酗酒者CP患病率更高。因此我們進一步探討了SNP位點與吸煙和飲酒的關(guān)系。在吸煙因素存在的情況下，具有c.493+49G>C位點改變的CP患病風險是對照組的4.667倍。其機制可能是吸煙引起的應激導致胰蛋白酶損傷而致病。多項研究表明c.180C>T位點與CP發(fā)病相關(guān)，但本研究未通過測序找到該SNP位點改變，其原因可能是樣本量小，未能驗證到該位點與CP的關(guān)系；或該位點基因型存在人群和地域差異。

綜上所述，本研究豐富了西南地區(qū)CTRC基因位點改變數(shù)據(jù)，表明不同地域、種族間位點變異分布有差異性。雜合錯義突變c.611G>A、c.703G>T和c.736C>T可能與西南地區(qū)漢族人群CP發(fā)病有關(guān)，SNP位點c.40+133G>A和c.493+49G>C可能是西南地區(qū)漢族人群CP易感位點，有位點c.493+49G>C改變的吸煙者可能有CP發(fā)病傾向。