FTO基因SNP rs9939609, rs1421085多態性與兒童青少年肥胖及其代謝指標的相關性研究

曹凌峰 羅飛宏 支滌靜 程若倩 沈水仙 楊 毅

飲食和生活習慣等環境因素對肥胖的發生和發展起著十分關鍵的作用，但遺傳因素構成了肥胖產生的內在基礎。研究表明至少5%的肥胖兒童存在單基因的缺陷，而5%左右的肥胖存在復雜的多基因數量性狀效應[1]。在FTO(fat-mass and obesity associated)基因研究中，單核苷酸多態性(SNP)與肥胖及其代謝綜合征的關聯研究是近10多年的研究熱點和重點，既往研究報道了很多肥胖易感基因SNP，如漿細胞膜糖蛋白1(ENPP1)基因、脂聯素(APM1)基因和轉錄因子7類似物2(TCF7L2)等，與BMI、胰島素抵抗和肥胖高度相關，但研究結果不盡一致，提示基因多態性的效應存在人種和人群間的差異[2]。

在既往研究基因變異和肥胖及其相關代謝指標的關聯研究中，一直沿用的基因組連鎖分析和候選基因策略雖然成功地確定了某些與肥胖相關的基因或基因變異，但僅有少數能在不同人群和人種中重現。主要的原因是實驗設計欠缺，表型評估差異，基因標記的選擇不一和研究樣本量的不同[3]。

全基因組關聯研究(GWAS)是在全基因組范圍內選擇遺傳變異進行基因分型, 比較病例和對照間每種變異頻率的差異, 計算變異與疾病的關聯強度, 選出最相關的變異進行驗證并最終確認與疾病相關。通過對SNP與不同癥狀的相關性研究，最終確定這些SNP對于特定疾病患病風險的影響[4]。GWAS研究不需要預先對特定關聯基因進行假設，因此主觀理論對研究結果影響較小。

在上海市科技攻關基金項目(024119013)的系列研究中，對AMP1基因SNP276和漿細胞膜糖蛋白1(K121Q)與肥胖和(或)超重的相關性進行分析，采用的是國外候選基因策略所獲得的關聯位點[5]，未發現這2個基因多態性與中國漢族兒童代謝指標的相關性，可能與上述原因有關。而本文采用的是基于國外GWAS研究策略所獲得的關聯位點，以了解FTO基因SNP rs9939609、 rs1421085多態性與中國漢族兒童青少年單純性肥胖及代謝指標間的相關性，為今后類似研究提供方法學選擇的借鑒。

1 方法

1.1 研究對象、人體指標測量、代謝指標檢測、基因分型和統計學方法 FTO基因SNP rs9939609、 rs1421085多態性的研究人群(納入標準、排除標準、正常對照組的設計，倫理和知情同意)，空腹血糖受損的判斷方法，胰島素抵抗指數(HOMA-IR)、胰島素敏感指數(QUICKI)的計算方法， DNA抽提方法，人體指標測量數據、一般資料數據、代謝指標檢測數據、統計學方法均與參考文獻[5]相同，在此不贅述。

本研究采用的基因探針型號為：rs 9939609:C_30090620_10;rs1421085：C_8917103_10,2×ABI Taqman Universal PCR master mix,No AmpErase UNG (lot No. E09854),20 ng基因組DNA模板,儀器采用ABI7900HT Sequence Detection System。重復性實驗：隨機取70份樣品，約占總數的10%，基因分型結果與原先測定的結果完全符合。

rs9939609和rs1421085分型圖見圖1。

圖1 FTO基因SNP rs9939609和rs1421085多態性的檢測結果

Fig 1 FTO gene SNP rs9939609 and rs1421085 allele distribution plot

Notes A:rs9939609; B:rs1421085

2 結果

2.1 研究對象的基本情況 與正常對照組比較，肥胖+超重組的BMI、血清空腹葡萄糖(FPG)、空腹胰島素(FIns)和三酰甘油(TG)水平均顯著高于正常對照組(P<0.01)，具體數據見參考文獻[5]。

2.2 基因分布頻率的平衡性檢驗 rs9939609和rs1421085多態性的基因分布頻率差異均無統計學意義(rs9939609：χ2=0.144,P=0.704；rs1421085:χ2=0.878,P=0.349)，基因型期望值與觀察值吻合較好，符合Hardy-Weinberg平衡。

2.3 各組等位基因頻率及基因型頻率比較 肥胖、超重和正常對照組rs9939609分型成功率分別為94.9%(224/236例)、97.9(234/239例)和95.9%(231/241名)。rs9939609 AA基因型頻率：肥胖組為2.7%，超重組為0.4%，正常對照組為1.7%；肥胖+超重組A等位基因頻率顯著高于正常對照組，P=0.048,OR=1.437，95%CI: 1.006～2.078(表1)。

肥胖、超重和正常對照組rs1421085分型成功率分別為92.8%(219/236例)、97.1%(232/239例)和95.4%(230/241名)。rs1421085 CC基因型頻率：肥胖組為2.7%，超重組為0.9%，正常對照組為1.7%，肥胖+超重組C等位基因頻率高于正常對照組，但差異無統計學意義，P=0.076，OR=1.388,95%CI:0.965～1.200(表2)。

2.4 FTO基因SNP rs9939609、rs1421085多態性與肥胖代謝指標間的相關性 rs9939609 TT基因型BMI顯著低于TA+AA基因型(P=0.000 3)；rs1421085 TT基因型BMI顯著低于TC+CC基因型(P=0.000 5)。rs9939609 TT基因型FPG、FIns、TG、總膽固醇(TC)、HOMA-IR、QUICKI與TA+AA基因型差異均無統計學意義；rs1421085 TT基因型FPG、FIns、TG、TC、HOMA-IR、QUICKI與TC+CC基因型差異均無統計學意義(表3)。

表1 肥胖組、超重組和正常對照組rs9939609等位基因頻率及基因型頻率比較[n(%)]

Notes OB:obesity group;OW:overweight group;NC:normal control group;1)obesity+overweight groupvsnormal control group,χ2=3.901,P=0.048,OR=1.437(95%CI:1.006-2.078);2)obesity+overweight groupvsnormal control group,χ2=5.57,P=0.061 7

表2 肥胖組、超重組和正常對照組rs1421085等位基因頻率及基因型頻率比較[n(%)]

Notes OB:obesity group;OW:overweight group;NC:normal control group;1)obesity+overweight groupvsnormal control group,χ2=3.148,P=0.076,OR=1.388(95%CI:0.965-1.200);2)obesity+overweight groupvsnormal control group,χ2=4.026,P=0.136 6

表3 FTO基因SNP rs9939609、rs1421085多態性與肥胖及其相關代謝指標間的相關性

Notes M:male;F:female

3 討論

人類與肥胖相關的候選基因的數量有幾百種，其SNP數量也相當龐大。1996年，Risch和Merikangas率先提出使用GWAS方法來確定基因多態性與疾病的相關性。近年來，隨著人類基因組計劃的順利完成和國際人類基因組單倍型圖計劃(HapMap Project)的不斷發展，以及高通量基因分型技術的成熟，使得同時研究幾萬，幾十萬，甚至幾百萬個SNP成為可能[6]，基因連鎖分析或候選基因策略，僅能在單基因缺陷的疾病或特性中發揮重要作用，但在特定疾病的多基因遺傳模式中，對確定多基因與復雜性疾病或特性的相關分析時作用相對有限。相對而言，GWAS能更精確地確定導致復雜性疾病的易感性基因,從而闡明這些復雜性疾病的遺傳學發病機制,是目前科學界公認的最為有效的復雜疾病的研究方法。2006年后，相繼有40多種復雜性疾病及其數量性狀(如BMI和身高等)的GWAS研究報道，并發現了某些以往未檢測到但與疾病或特性密切相關的基因和多態性位點，如FTO基因。

FTO基因位于染色體16q12.2，在藻類和脊椎動物中高度保守，但在無脊椎動物、真菌和綠色植物中不表達[7，8]。FTO基因屬于非血紅色加雙氧酶家族，編碼2-酮戊二酸-依賴型核酸脫甲基酶[9，10]。FTO基因mRNA在大腦內，尤其下丘腦mRNA的表達最豐富，但在其他組織也有表達[7，9]。動物實驗研究發現，小鼠進食時下丘腦弓狀核FTO基因mRNA表達水平上調，禁食時則下調[9，11]，但在大鼠中調控模式與此相反[7]，FTO基因能夠使DNA去甲基化，并參與DNA修復、脂肪酸代謝和翻譯后修飾。

2007年，Frayling等[12]率先報道了位于FTO基因第1內含子的rs9939609與成人和兒童肥胖具高度相關性，隨后在歐洲幾項獨立研究中，也分別報道了該基因多態性與高加索人群、西班牙人群BMI和肥胖的相關性。FTO基因是至今為止在歐洲人群中唯一的、重復成功的肥胖易感基因[13]。在東亞地區，在有限的報道中結果并不一致。目前在中國僅見臺灣、北京和上海地區的成年人群FTO基因多態性與肥胖的關聯報道。Chang等[14]報道，FTO基因rs9939609多態性與中國臺灣成人肥胖和BMI存在相關性，但與2型糖尿病無相關性；但在上海和北京的成人研究中未發現其與肥胖、BMI和糖尿病之間的關聯性[15]。在日本人群的研究結果也不盡相同，Horikoshi等[16]未發現rs9939609多態性與BMI和肥胖的相關性；Omori等[17]研究得出陽性結果。提示研究結果可能與人群樣本的選擇有關。本研究發現，肥胖+超重組rs9939609 A等位基因頻率顯著高于正常對照組，rs1421085 C等位基因頻率高于正常對照組，但差異無統計學意義；rs1421085 TC +CC基因型和rs9939609 TA +AA 基因型的BMI顯著高于TT基因型。因此提示在中國漢族兒童青少年中rs1421085、rs9939609多態性均與肥胖的體脂形成有密切的關聯性。

本研究正常對照組rs9939609 A等位基因和rs1421085 C等位基因的頻率分別為9.7%和9.8%,其中A等位基因頻率略低于中國臺灣人群(12.6 %)[14]，而AA和CC的基因型頻率均為1.7%, 與文獻報道的結果一致[14]，但明顯低于歐洲人群(45%)[18]和非洲裔美國人(21%)[13]，提示FTO基因該2個多態性位點基因頻率存在顯著的種族差異。

本研究人群的rs9939609和rs1421085多態性與BMI顯著相關，但與肥胖其他相關代謝指標包括FPG、FIns、TC、TG和HOMA-IR等均無顯著相關性；與日本[15]和中國[16]的研究結果較為一致。1項納入17 037名高加索人群的Meta分析結果顯示，FTO基因SNP rs9939609和rs1421085多態性與FPG、TG、FIns存在相關性[19]。提示rs9939609和rs1421085多態性對肥胖患者代謝指標的影響有種族和地區差異。

本研究首次證實了FTO基因多態性與中國漢族兒童青少年BMI的相關性，雖然其rs9939609 A等位基因頻率遠低于歐洲人群，但該等位基因對BMI的作用與歐洲人群較為相似，但對代謝指標影響卻有所不同。由于中國是個多民族國家，在其他民族兒童青少年中FTO基因該位點多態性是否還有類似的作用，尚有待進一步研究證實。

本研究的不足之處與局限性：本研究對象僅限于上海某地區的漢族兒童青少年。作為全基因組相關研究，研究對象的數量尚顯不足。

[1]Bouchard C. Childhood obesity: are genetic differences involved? Am J Clin Nutr,2009,89(5):1494-1501

[2]Bell CG, Walley AJ, Froguel P. The genetics of human obesity. Nat Rev Genet,2005,6(3):221-234

[3]Rankinen T, Zuberi A, Chagnon YC, et al. The human obesity gene map: the 2005 update. Obesity (Silver Spring)，2006，14(4): 529-644

[4]Kuniholm MH. Replicating genome-wide association studies. Science，2007，318(5849):390-391

[5]Cao LF(曹凌峰),Luo FH,Zhi DJ,et al.The relationship between SNP276 polymorphism of adiponectin gene and obesity in children and adolescents.Chin J Evid Based Pediatr(中國循證兒科雜志),2009，4(6):525-529

[6]Goodman E, Daniels SR, Morrison JA, et al. Contrasting prevalence of and demographic disparities in the World Health Organization and National Cholesterol Education Program Adult Treatment Panel Ⅲ definitions of metabolic syndrome among adolescents. J Pediatr,2004,145(4):445-451

[7]Fredriksson R, Hagglund M, Olszewski PK, et al. The obesity gene, FTO, is of ancient origin, upregulated during food deprivation and expressed in neurons of feeding-related nuclei of the brain. Endocrinology， 2008，149(5):2062-2071

[8]Robbens S, Rouzé P, Cock JM, et al. The FTO gene, implicated in human obesity, is found only in vertebrates and marine algae. J Mol Evol，2008，66(1): 80-84

[9]Gerken T, Girard CA, Tung YC, et al. The obesity-associated FTO gene encodes a 2-oxoglutarate-dependent nucleic acid demethylase.Science，2007，318(5855): 1469-1472

[10]Fernández AM, Kim JK, Yakar S, et al.Functional inactivation of the IGF-I and insulin receptors in skeletal muscle causes type 2 diabetes. Genes Dev，2001，15(15): 1926-1934

[11]Stratigopoulos G, Padilla SL, LeDuc CA, et al.Regulation of Fto/Ftm gene expression in mice and humans. Am J Physiol Regul Integr Comp Physiol ，2008，294(4):1185-1196

[12]Frayling TM, Timpson NJ, Weedon MN，et al. A common variant in the FTO gene is associated with body mass index and predisposes to childhood and adult obesity. Science，2007，316(5826):889-894

[13]Loos RJ, Bouchard C. FTO: the first gene contributing to common forms of human obesity. Obes Rev，2008，9(3):246-250

[14]Chang YC, Liu PH, Lee WJ,et al. Common variation in the fat mass and obesity-associated (FTO) gene confers risk of obesity and modulates BMI in the Chinese population.Diabetes，2008，57(8):2245-2252

[15] Li H, Wu Y, Loos RJ,et al. Variants in the fat mass- and obesity-associated (FTO) gene are not associated with obesity in a Chinese Han population. Diabetes，2008，57(1):264-268

[16]Horikoshi M, Hara K, Ito C,et al. Variations in the HHEX gene are associated with increased risk of type 2 diabetes in the Japanese population. Diabetologia，2007，50(12):2461-2466

[17]Omori S, Tanaka Y, Takahashi A,et al. Association of CDKAL1, IGF2BP2, CDKN2A/B, HHEX, SLC30A8, and KCNJ11 with susceptibility to type 2 diabetes in a Japanese population. Diabetes，2008，57(3):791-795

[18]Peeters A, Beckers S, Verrijken A,et al. Variants in the FTO gene are associated with common obesity in the Belgian population. Mol Genet Metab，2008，93(4)：481-484

[19]Freathy RM, Timpson NJ, Lawlor DA,et al. Common variation in the FTO gene alters diabetes-related metabolic traits to the extent expected given its effect on BMI. Diabetes，2008，57(5):1419-1426