研究目的进食障碍是指通常起病于青少年期或成年早期的进食行为的严重不协调,患者以女性为主,常表现为因过度关注体重或身材而导致严重的进食减少或过度饮食,并因此感到心里极度痛苦,并造成患者的社会功能受损。依据美国精神病学会最新的诊断标准,进食障碍主要分为神经性厌食、神经性贪食、暴食障碍和其他未分型的进食障碍四种主要亚型。进食障碍的患者常共病其他精神疾病,如强迫障碍、抑郁障碍及焦虑障碍等。进食障碍患者的预后不良,长期的过度营养不良引起的心律失常等并发症,或是由于共病的抑郁障碍造成自杀等不良后果均可能是进食障碍的死因。病因学研究一直是进食障碍的研究热点之一。依据目前的研究结果,进食障碍的病因除了与社会、家庭环境和人格特质等因素有关外,遗传因素在其发病中也存在较大影响。按照遗传学的观点,目前认为进食障碍是多条神经病理通路的基因异常共同作用导致的,其中五羟色胺(serotonin,5-HT)信号通路、多巴胺(dopamine,DA)信号通路、脑源性神经营养因子(brain-derived neurotrophic factor,BDNF)、Ghrelin和雌激素(estrogen)等受到较多关注。本研究的目的是通过meta分析来评价在进食障碍发病过程中某些单核苷酸多态性(Single Nucleotide Polymorphism,SNP)的功能,选取的单核苷酸多态性分别来自五羟色胺、多巴胺、脑源性神经营养因子和Ghrelin通路,以期为进食障碍的病因、诊断、治疗和预后提供有效的参考依据。研究方法检索截止到2017年12月31日之前公开发表的关于SNP与进食障碍相关性的收录在Web of Science、PubMed或中国知网(China National Knowledge Internet,CNKI)数据库内的所有文献,文献的写作语言限制在中文和英文。纳入的文献需是病例对照研究,且以健康人群为对照组被试。采用纽卡斯尔-渥太华量表对文献质量进行评估,将得分在6星及6星以上的文献纳入研究进行数据提取和分析。将纳入文章按照SNP位点进行汇总。采用RevMan 5.3软件分别进行OR值的合并和95%置信区间及相关效应量的估计。对于未明确遗传模式的SNP分别采用:等位基因模型、显性模型、隐性模型、共显性模型和超显性模型进行meta分析;已明确遗传模式地采用等位基因模型和相应模型进行分析。继而按照被试人种以及疾病亚型不同分别进行亚组分析。采用逐一剔除法进行敏感性分析,漏斗图用于估计纳入文献的发表偏倚。结果最终纳入141篇针对不同基因位点与进食障碍相关性的研究,将同一多态性位点的文献数量大于3篇的研究进行meta分析。本文共选取5个多态性位点进行分析,分别是:五羟色胺转运体基因连锁多态性区域(serotonin transporter gene linked polymorphic region,5-HTTLPR)、5-HT2A受体(serotonin 2A receptors,5-HT2ARs)启动子基因(-1438A/G)、儿茶酚氧位甲基转移酶(catechol-O-methyltransferase,COMT)基因(rs4680)、BDNF基因(rs6265)和Ghrelin基因(408C/A)。5-HTTLPR基因相关文献共有13篇,共包含2515例进食障碍患者和健康对照组2847例。按照等位基因模型,5-HTLPPR与高加索人种AN(OR=1.09,P=0.23)或BN(OR=1.12,P=0.57)患病无相关性,但增加蒙古利亚人种AN的患病风险(OR=1.55,P<0.001);按照显性模型,5-HTLPPR与高加索人AN(OR=1.06,P=0.55)和蒙古利亚人AN患病无关(OR=1.54,P=0.05),也与高加索人BN(OR=1.23,P=0.46)的患病无关。5-HT2A(-1438A/G)基因相关文献共检索到18篇,共有进食障碍患者2463例和健康对照组2713例。在高加索人种中-1438位点A等位基因突变是AN的患病的风险因子[等位基因模型(OR=1.36,P=0.003)、显性模型(OR=1.39,P=0.02)、隐性模型(OR=1.76,P=0.0006)、共显性模型(OR=2.02,P=0.003)];但是与蒙古利亚人种AN的患病无关[等位基因模型(OR=0.91,P=0.43)、显性模型(OR=0.87,P=0.46)、隐性模型(OR=0.88,P=0.41)、共显性模型(OR=0.81,P=0.42)、超显性模型(OR=0.99,P=0.96)]。5-HT2A(-1438A/G)是高加索人BN的患病风险因子[等位基因模型(OR=1.42,P=0.02)、显性模型(OR=1.81,P<0.000)、隐性模型(OR=1.66,P<0.000)],但关于蒙古利亚人的研究过少,无法分析。DA通路COMT基因rs4680(Val158Met)相关文献共有10篇,共有进食障碍患者2637例和健康对照组3776例。亚组分析表明,高加索人种AN发病与COMT基因rs4680 A等位基因突变无关[等位基因模型(OR=0.84,P=0.09)、显性模型(OR=0.90,P=0.47)、隐性模型(OR=0.99,P=0.90)、共显性模型(OR=1.04,P=0.62)、超显性模型(OR=0.96,P=0.50)];该基因突变将降低高加索人种BN的发病风险[等位基因模型(OR=0.76,P=0.02)、显性模型(OR=0.65,P=0.03)、共显性模型(OR=0.58,P=0.03)]。蒙古利亚人中相关文献数量太少无法进行meta分析。BDNF基因rs6265(Val66Met)相关文献纳入14篇,共有进食障碍患者4754例,健康对照组5215例。根据亚组分析结果,BDNF基因rs6265位点携带A(Met)等位基因将增加高加索人种AN[在等位基因模型(OR=1.13,P=0.005)、显性模型(OR=1.17,P=0.004)]和BN[等位基因模型(OR=1.36,P<0.000)、隐性模型(OR=1.74,P=0.03)、共显性模型(OR=1.91,P=0.01)]的患病风险,但不增加蒙古利亚人AN的患病风险。蒙古利亚人种有关BN的相关研究过少,无法进行meta分析。Ghrelin基因408C/A(Leu72Met)相关文献纳入4篇,共有进食障碍患者713例,健康对照组702例。根据疾病亚型进行亚组分析,Ghrelin基因408C/A(Leu72Met)位点A(Met)等位基因突变并不增加AN的患病风险[等位基因模型(OR=1.44,P=0.17)、显性模型(OR=1.30,P=0.12)、隐性模型(OR=1.09,P=0.85)、共显性模型(OR=1.18,P=0.71)、超显性模型(OR=0.76,P=0.11)];但与BN[等位基因模型(OR=1.46,P=0.02)、显性模型(OR=1.78,P=0.003)]和BED[等位基因模型(OR=2.31,P=0.005)、显性模型(OR=2.55,P=0.003)]的患病风险增加有关。结论在5-HT通路中,5-HTTLPR与进食障碍的发病无相关性;但5-HT2A(-1438A/G)基因多态性与进食障碍的发病风险增高有关,尤其是在高加索人种中。在DA通路rs4680位点单核苷酸多态性与进食障碍发病无关。BDNF基因rs6265基因多态性与进食障碍发病有关,尤其是高加索人种。Ghrelin基因408C/A与AN发病无关,但可能增加BN、BED的患病风险。未来需要更多针对亚洲人种的有关研究以为了解进食障碍的遗传学病因提供更多的循证依据。