自噬相關基因頻率分布特征與肺結核相關性研究*

李 穎,黃 秦,陳維賢,李升錦△

(1.重慶醫科大學附屬第二醫院呼吸內科 400010；2.重慶醫科大學附屬第一醫院呼吸內科 400016；3.重慶醫科大學附屬第二醫院檢驗科 400010)

論著·臨床研究

自噬相關基因頻率分布特征與肺結核相關性研究*

李 穎1,黃 秦2,陳維賢3,李升錦1△

(1.重慶醫科大學附屬第二醫院呼吸內科 400010；2.重慶醫科大學附屬第一醫院呼吸內科 400016；3.重慶醫科大學附屬第二醫院檢驗科 400010)

目的 根據肺結核患者自噬相關基因(ATG)頻率分布特征探討ATG位點多態性與肺結核易感相關性。方法 應用SequenomMassArray質譜陣列技術對202例肺結核患者(病例組)及222例健康體檢者(對照組)ATG的18個單核苷酸多態性(SNP)位點進行基因分型，統計分析各位點基因型與肺結核易感相關性。結果 在肺結核患者和對照組中,校正了性別、年齡等因素后,通過二元logistic回歸分析發現ATG4A基因rs807185和 rs5973822位點基因型和等位基因頻率分布差異有統計學意義(P<0.05)。通過體質量指數(BMI)分層分析,發現這一差異在高BMI人群中更為顯著，且均有病例組人群分布頻率低于對照組的表現，其余16個SNP位點差異均無統計學意義。結論 ATG4A基因rs5973822和rs807185位點多態性可能與肺結核易感性相關聯，且在高BMI群體中差異更顯著。

結核,肺；自噬；基因；多態性,單核苷酸

肺結核是由結核分枝桿菌持續感染導致的慢性傳染病，機體發病與否主要取決于宿主、自然環境和菌株三因素。近來眾多研究表明，基因多態性可能影響肺結核的患病風險，暗示個體遺傳因素與結核病的發生與發展密切相關[1]。本課題組前期研究證實自噬相關基因(autophagy associated gene，ATG)4A基因單核苷酸多態性(single nucleotide polymorphism,SNP)位點rs807185的等位基因多態性與中國西南地區人群肺結核易感性呈負相關[2]；而自噬通路其他基因位點的遺傳變異可能同樣在結核病免疫中扮演重要角色。故本文應用SequenomMassArray質譜陣列技術對202例肺結核患者及222例健康者ATG的18個SNP位點進行基因分型，統計分析各位點基因型與肺結核易感相關性。

1 資料與方法

1.1 一般資料 本研究血液樣本由重慶醫科大學附屬第二醫院提供，共202例(男118例，女84例)肺結核患者(病例組)納入分析，平均(42±13)歲，同期在參加體檢的健康人群中挑選222例(男119例，女103例)進行對照研究,平均年齡(41±14)歲。個體之間均無血緣關系，為重慶市漢族人群；征得知情同意后簽署知情同意書并通過了重慶醫科大學倫理委員會的審批，符合醫學倫理學實驗標準。病例組均通過痰結核菌涂片或培養、影像學檢查并結合相應的臨床癥狀確診為肺結核，對照組均除外糖尿病、嚴重免疫系統疾病、惡性腫瘤及HIV等合并癥。

1.2 SNP位點來源 基于課題組前期對ATG4A及ATG16L1基因多態性與肺結核相關性研究的基礎上，在Ensembl數據庫中進行檢索，收集和整理相關研究和已有報道中指示與肺結核密切相關的SNP功能位點，采用pairwise-tagging方法，根據重慶地區漢族人群基因型數據挑選單體標簽的SNP，高分辨率溶解(HRM)方法再次進行挑選，并對少數樣本測序，最終采集肺結核典型患者生物樣本對挑選的基因及具體基因型進行驗證。

1.3 實驗標本、試劑及儀器

1.3.1 研究對象標本 為肺結核確診患者和健康人群的加入EDTA的外周靜脈血。

表1 自噬相關基因SNP位點引物序列設計

表2 自噬相關基因SNP位點詳細信息

1.3.2 實驗試劑簡介 人淋巴細胞分離液(挪威Axis-shield公司)、血基因組DNA提取試劑盒(北京天恩生物技術公司)、多種DNA酶(日本TaKaRa公司)、PCR引物(北京奧科生物有限責任公司)等。

1.3.3 實驗儀器簡介 GeneAmp9700型PCR儀(美國Applied Biosystems公司)、臺式離心機(雷博爾離心機有限公司)、全自動凝膠成膠系統(ViberLounmat公司)、高速臺式冷凍離心機(Coulter Beckman公司)等。

1.4 實驗方法

1.4.1 SNP篩選及引物設計 所有研究對象均采集靜脈血，2 mL容量一式三份，分別用于實驗室檢測、復檢和備查，采用EDTA-K2抗凝，置于-80 ℃低溫保存。用血液基因組提取試劑盒(北京天根公司)提取外周血基因組DNA，測定DNA純度和濃度，見表1。

1.4.2 基因分型 利用SequenomMassArray質譜陣列技術對ATG4A基因的18個SNP位點進行基因分型。聚合酶鏈式反應(Polymerase chain reaction，PCR)用于擴增位于兩段已知序列之間的DNA區域，堿性磷酸酶(shrimp alkaline phosphatase，SAP)反應,去除未反應完的dNTP，最后進行單堿基延伸PCR反應和點樣及Mass ARRAY分析。

1.5 統計學處理 采用SPSS19.0進行分析。人口統計學中的年齡差異采用獨立樣本t檢驗進行分析，性別差異采用χ2檢驗。病例組與對照組的基因型、等位基因分布比較采用χ2檢驗，應用二元Logistic回歸模型計算不同基因型相對風險度的比值比(OR)及95%可信區間(CI)，檢驗水準α取值0.05，雙側概率。

2 結 果

2.1 調查對象基本情況 兩組人群在性別、年齡、吸煙史、飲酒史等一般資料上實行了適當的匹配，保證研究對象一般資料基本一致。

2.2 ATG 18個SNP多態性位點具體位置 所有調查對象ATG的SNP多態性位點基因型頻率分布符合Hardy-Weinberg平衡定律，基因型頻率統計概率P>0.05。本次研究所檢測的ATG 18個SNP多態性位點具體位置見表2。

2.3 ATG的SNP多態性位點基因型和肺結核易感性的關系 分段比較兩組對象18個SNP位點等位基因及基因型的差異，結果顯示病例組rs5973822位點GG基因型(G等位基因)的頻率(3.5%)顯著低于對照組(8.6%)，差異有統計學意義(P=0.016)；病例組rs807185位點AT基因型(A等位基因)的頻率(24.8%)顯著低于對照組(35.6%),差異有統計學意義(P=0.009)。其他位點基因型及等位基因頻率分布在兩組對象中差異無統計學意義(P>0.05)。見表3。

2.4 兩組對象中rs5973822和rs807185位點等位基因及基因型頻率在不同體質量指數(BMI)對象中的分布 BMI<18.5 kg/m2的對象比較，rs807185和rs5973822基因位點頻率分布差異均無統計學意義(P>0.05)。BMI>18.5 kg/m2的人群rs807185基因位點的A等位基因(P=0.021，OR=0.580，95%CI0.365～0.922)和AT基因型(P=0.009，OR=0.434，95%CI0.230～0.816)差異有統計學意義，rs5973822基因位點的G等位基因(P=0.034，OR=0.606，95%CI0.381～0.965)和GG基因型(P=0.039，OR=0.302，95%CI0.092～0.989)差異有統計學意義。見表4。

表3 ATG的SNP位點等位基因及基因型在兩組對象中的頻率分布

NS:無意義

表4 兩組對象中rs5973822和rs807185位點等位基因及基因型頻率在不同BMI人群中的分布

續表4 兩組對象中rs5973822和rs807185位點等位基因及基因型頻率在不同BMI人群中的分布

-:表示無計算意義

3 討 論

自噬是真核細胞中普遍存在的防御機制[3]，能依賴溶酶體直接殺傷胞內病原體，并協助抗原提呈，在天然免疫和特異性免疫中發揮重要作用[4]。自噬功能障礙繼發的免疫反應和結核病的發生與發展密切相關[5]。自噬體的形成是自噬過程的關鍵環節，依賴眾多自噬基因協調配合，其中包括ATG12-ATG5和微管相關蛋白輕鏈-磷脂酰乙醇胺(ATG8-PE)兩套泛素化連接系統[6]。ATG4A是C-54家族中的一種半胱氨酸蛋白酶，通過控制ATG8的裂解進而調控其與ATG8-PE的共價結合與解離而影響自噬體的形成[7]。本課題組前期研究已證實，ATG4A基因rs807185位點多態性與肺結核易感性呈負相關，但自噬通路其他基因多態性與結核病的相關性仍少見報道。因此，筆者擴大研究，深入探討ATG頻率分布特征與肺結核易感相關性。

本次研究通過對比分析健康者與肺結核患者18個SNP位點等位基因及基因型頻率，發現ATG4A基因rs807185和rs5973822位點多態性與肺結核的發病風險相關聯，其中rs5973822 SNP 位于3′-UTRs(非翻譯區)， 其多態性可能影響基因表達水平， rs807185 SNP位于內含子區域。目前有多項研究證實，內含子可能涉及mRNA轉錄、翻譯及組織特異性表達[8-10]。因此，筆者推測rs807185和rs5973822位點多態性可能通過影響自噬功能進而影響細胞免疫反應；其余16個SNP位點多態性與肺結核發病風險無明顯相關性，考慮與樣本數量不足，檢測方法局限及種族地域性差異有關。

分層研究顯示，ATG4A基因rs807185和rs5973822位點多態性與高BMI人群結核病患病風險關聯性更大。BMI是營養狀況的反映指標，高BMI人群日常可能攝取高能量和高蛋白，免疫功能增強，降低機體對結核病的易感性[11-12]。而BMI<18.5 kg/m2被認定為營養不良[13]，常伴隨常量營養素及微量營養素缺乏；通過加速胸腺、淋巴組織萎縮，擾亂單核巨噬細胞、T淋巴細胞和細胞因子的相互作用[14]，從而降低免疫應答，增加宿主對結核桿菌的感染風險，加快潛伏性感染進展為活動性肺結核進程[15]。

綜上所述，筆者考慮ATG4A rs807185和rs5973822位點多態性可能降低結核病發病風險，且在高BMI群體中更為顯著。但本次研究仍有許多不足之處，后期筆者將會擴大樣本量，檢測更多的自噬相關基因，設計有效的功能試驗，為結核病的風險預測、診斷及治療提供理論基礎。

[1]Ghamari E,Farnia P,Saif S， et al.Susceptibility to pulmonary tuberculosis:host genetic deficiency in tumor necrosis factor alpha (TNF-α) gene and tumor necrosis factor receptor 2(TNFR2)[J].Int J Mycobacteriol，2016，5(Suppl 1):S136-137.

[2]黃秦,趙鳳容,李霞，等．ATG4A 及 ATG16L1 基因多態性與肺結核易感性的關系[J]．重慶醫科大學學報,2015,40(4):599-603.

[3]Xu Y,Eissa NT.Autophagy in innate and adaptive immunity[J].Proc Am Thorac Soc,2010,7(1):22-28.

[4]Ní Cheallaigh C,Keane J,Lavelle C,et al.Autophagy in the immune response to tuberculosis:clinical perspectives[J].Clin Exp Immunol,2011,164(3):291-300.

[5]Ponpuak M,Davis AS,Roberts EA,et al.Delivery of cytosolic components by autophagic adaptor protein p62 endows autophagosomes with unique antimicrobial properties[J].Immunity,2010,32(3):329-341.

[6]Romanov J,Walczak M,Ibiricu I,et al.Mechanism and functions of membrane binding by the Atg5-Atg12/Atg16 complex during autophagosome formation[J].EMBO J,2012,31(22):4304-4317.

[7]Weidberg H,Shvets E,Shpilka T,et al.LC3 and GATE-16/GABARAP subfamilies are both essential yet act differently in autophagosome biogenesis [J].EMBO J，2010,29(11):1792-1802.

[8]ENCODE Project Consortium,Birney E,Stamatoyannopoulos JA,et al.Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project[J].Nature,2007,447(7146):799-816.

[9]Ortiz-Cuaran S,Cox D,Villar S,et al.Association between TP53 R249S mutation and polymorphisms in TP53 intron 1 in hepatocellular carcinoma[J].Genes Chromosomes Cancer,2013,52(10):912-919.

[10]Cheng J,Kapranov P,Drenkow J,et al.Transcriptional maps of 10 human chromosomes at 5-nucleotide resolution[J].Science,2005,308(5725):1149-1154.

[11]Hanrahan F,Golub E,Mohapi L,et al.Body mass index and risk of tuberculosis and death[J].AIDS,2010,24(10):1501-1508.

[12]Leung CC,Lam TH,Chan WM,et al.Lower risk of tuberculosis in obesity[J].Arch Intern Med,2007,167(12):1297-1304.

[13]Cegielski JP,Arab L,Cornoni-Huntley J.Nutritional risk factors for tuberculosis among adults in the United States,1971-1992[J].Am J Epidemiol,2012,176 (5):409-422.

[14]Ezeamama AE,Mupere E,Oloya J,et al.Age,sex,and nutritional status modify the CD4+ T-cell recovery rate in HIV-tuberculosis co-infected patients on combination antiretroviral therapy[J].Int J Infect Dis,2015,35(1):73-79.

[15]Elias D,Mengistu G,Akuffo H,et al.Are intestinal helminths risk factors for developing active tuberculosis? [J].Trop Med Int Health，2006,11(4):551-558.

Research on correlation between autophagy related gene frequency distribution characteristics and pulmonary tuberculosis*

LiYing1，HuangQing2，ChenWeixian3，LiShengjing1△

(1.DepartmentofRespiration,SecondAffiliatedHospital,ChongqingMedicalUniversity,Chongqing400010,China;2.DepartmentofRespiration,FirstAffiliatedHospital,ChongqingMedicalUniversity,Chongqing400016,China;3.DepartmentofLaboratoryMedicine,SecondAffiliatedHospital,ChongqingMedicalUniversity,Chongqing400010,China)

Objective To investigate the correlation between autophagy related gene(ATG) locus polymorphism and the pulmonary tuberculosis(PTB) susceptibility according to the frequency distribution characteristics of autophagy related gene in the patients with PTB.Methods Eighteen single nucleotide polymorphism (SNP) loci of autophagy related genes in 202 patients with PTB as the case group and 222 healthy controls were genotyped by SequenomMassArray mass spectrometry array technology.The correlation between the each locus genotype and the PTB susceptibility was statistically analyzed.Results In the PTB patients group and healthy control group,after correcting the factors of sex and age,the binary Logistic regression analysis found that the frequency distribution of genotype and allele had statistical difference between rs5973822 and rs807185 sites in ATG4A gene (P<0.05).The stratified analysis by body mass index (BMI) found that this difference was more significant in the high BMI population，moreover the distribution frequency in the patient group was lower than that in the control group,and the other 16 SNP loci had no statistical difference.Conclusion ATG4A gene rs5973822 and rs807185 loci polymorphism may be negatively correlated with PTB susceptibility,moreover which is more significant in the high BMI group.

tuberculosis,pulmonary;autophagy;genes;polymorphism,single nucleotide

10.3969/j.issn.1671-8348.2017.22.007

國家自然科學基金資助項目(81171628)；重慶市衛計委基金資助項目(2012-2-370)。作者簡介：李穎(1990-)，住院醫師，碩士，主要從事呼吸系統免疫防御機制研究。△

，E-mail：lsj1025@163.com。

A

1671-8348(2017)22-3046-04

2017-02-28

2017-04-16)