MTHFRC677T基因多态性与孤独症谱系障碍易感性关系的Meta分析

张宇昕，盖聪，杨璐平，马浩洁，张锦坤，孙红梅，胡蝶

【摘要】：目的探讨亚甲基四氢叶酸还原酶（MTHFR）基因多态性与孤独症谱系障碍（ASD）易感性的关系。方法检索中国知网、万方数据、PubMed、Embase、Proquest、Webofscience等数据库，收集MTHFRC677T基因多态性与ASD关联的病例对照研究，检索时间至2020年6月30日。应用Stata14.0软件进行Meta分析。结果共纳入16篇文献，其中病例组2149例，对照组2251例。Meta分析结果显示在等位基因模型（Tvs.C，OR=1.51，95%CI：1.15~1.99，P=0.003）、显性模型（TT+CTvs.CC，OR=1.67，95%CI：1.21~2.29，P=0.002）中，MTHFRC677T多态性与ASD发生具有相关性。而在隐性模型（TTvs.CT+CC，OR=1.20，95%CI：0.82~1.76，P=0.340）和纯合子模型（TTvs.CC，OR=1.43，95%CI：0.92~2.22，P=0.113）中MTHFRC677T多态性与ASD易感性无关。种族亚组分析显示，高加索人MTHFRC677T多态性仅在显性模型（OR=1.30，95%CI：1.04~1.63，P=0.024）中与ASD发生相关，而亚洲人在等位基因模型（OR=1.85，95%CI：1.13~3.01，P=0.014）和显性模型（OR=1.99，95%CI：1.18~3.36，P=0.010）中和ASD发病有关。结论MTHFRC677T基因多态性与ASD易感性相关。

【关键词】：孤独症谱系障碍；亚甲基四氢叶酸还原酶（NADPH）；多态性，单核苷酸；Meta分析

中图分类号：R748文献标志码：ADOI：10.11958/20202256基金项目：国家自然科学青年基金资助项目（81803857）作者单位：北京中医药大学中医学院（邮编102488）

作者简介：张宇昕（1996），女，硕士在读，主要从事脑病研究。

前言

孤独症谱系障碍（ASD）是一种常见的儿童神经发育障碍，其典型症状包括社交行为和沟通缺陷、兴趣受限和重复行为等［1-2］。ASD的病因尚不完全明确，其发病受遗传、环境等多种因素的影响［3］。DNA甲基化作为最重要的表观遗传学修饰之一，对大脑的生长、发育和认知功能至关重要［4］。亚甲基四氢叶酸还原酶（MTHFR）为叶酸代谢途径的关键酶，在DNA的合成及甲基化中具有重要作用［5］。MTHFR基因多态性可通过改变DNA甲基化和叶酸代谢降低MTHFR酶活性，可能与ASD等神经障碍有关［6］。目前已有较多关于MTHFR基因多态性与ASD风险之间的研究，但结论存在一定争议。因此，本研究采用Meta分析的方法对已有研究进行统计，以明确MTHFRC677T基因多态性与ASD易感性之间的关系。

1资料与方法

1.1文献检索对中国知网（CNKI）、万方数据、PubMed、Embase、Proquest、Webofscience等数据库进行全面的文献检索。中文检索词为“亚甲基四氢叶酸还原酶”、“孤独症谱系障碍”、“孤独症”、“自闭症”；英文检索词为“methylenetetrahydrofolatereductase”、“MTHFR”、“autismspectrumdisorders”、“autism”、“ASD”，检索2020年6月30日之前发表的所有MTHFR基因多态性与ASD易感风险关系的中英文相关出版物。

1.2文献纳入与排除标准纳入标准：（1）病例组患者符合《精神障碍诊断和统计手册》（diagnosticandstatisticalmanualofmentaldisorders，DSM）或儿童孤独症评定量表（childhoodautismratingscale，CARS）诊断标准。（2）有关MTHFRC677T基因多态性与ASD发生风险关系的公开发表文献。（3）病例组和健康对照组的基因型频率符合Hardy-Weinberg平衡

（HWE）。排除标准：（1）重复发表文献、动物研究及综述文献。（2）非病例对照研究文献。

1.3文献筛选及资料提取通过阅读文献的标题和摘要进行初筛，而后阅读全文进行二次筛查，最终根据纳入标准决定文献是否纳入，意见不同时通过讨论或第3位研究者协助判断。从每项确定的研究中提取出以下信息：第一作者、发表年份、国家、种族、样本量、基因型信息和所有研究等位基因的频率。如果文献中没有给出上述信息，则通过联系作者获得。

1.4文献质量评价由2名研究者独立进行评价，结论出现不一致时讨论协商或第3研究者协助判断解决。计算纳入文献中对照组的基因型频率是否符合HWE，若P≤0.05，代表文献质量较低；若P＞0.05，则该文献质量较高。同时使用纽卡斯尔-渥太华量表（NOS）进行文献质量评价，从研究对象的选择、组间可比性和暴露因素的测量3个方面进行评价，共包含8个条目，共计9分，得分≥5分为高质量研究。

1.5统计学方法所有的统计分析使用Stata14.0软件，MTHFR多态性与ASD之间的关系采用OR值和95%CI进行评估。分别计算等位基因模型（Tvs.C），显性模型（TT+CTvs.CC），纯合子模型（TTvs.CC）和隐性模型（TTvs.CT+CC）的OR值及其95%CI，并进行种族亚组分析。采用χ2检验判定基因型在对照组的分布是否符合HWE。采用I2检验文献的异质性，检验水准α=0.05，若P＞0.05且I2≤40%，说明异质性小，选择固定效应模型；若P≤0.05且I2>40%，说明各研究之间存在较高异质性，选择随机效应模型进行分析。通过敏感性分析来检验结果的稳定性，为了控制潜在的发表偏倚，采用漏斗图并通过Egger线性回归检验进行分析；P＜0.05为差异有统计学意义。

2结果

2.1纳入文献基本特征 经文献筛选后共纳入16篇文献［7-22］，其中病例组2149例，对照组2251例，文献筛选流程及结果见图1，纳入文献的基本资料见表1，NOS质量评分均≥7分。图1文献筛选流程图

2.2MTHFR分析结果

2.2.1MTHFRT等位基因与ASD发病之间的Meta分析 各研究结果存在异质性（I2=85.7%，P＜0.001），采用随机效应模型进行Meta分析。ASD组和对照组携带T等位基因分别为1528例、2770例，C等位基因分别为2770例、3084例。以T等位基因为暴露因素，C等位基因为非暴露因素进行分析，结果显示，携带T等位基因者发生ASD的风险是携带C等位基因的1.51倍（OR=1.51，95%CI：1.15~1.99，P=0.003）。种族亚组分析显示，亚洲人携带T等位基因者发生ASD的风险是携带C等位基因人群的

1.85倍（P=0.014）；而在高加索人及非洲人组中，携带T等位基因者发生ASD的风险与携带C等位基因无关（高加索人OR=1.20，P=0.102；非洲人OR=2.17，P=0.354），见图2。表1MTHFRC677T纳入研究的基本特征。图2MTHFR基因C677T多态性（Tvs.C）与ASD发病的Meta分析森林图

2.2.2显性模型TT+CTvs.CC与ASD发病的Meta分析 各研究结果存在异质性（I2=80.8%，P＜0.001），采用随机效应模型进行Meta分析。ASD组和对照组TT+CT基因型分别为1248例、1149例，CC基因型分别为901例、1102例。以TT+CT基因型为暴露因素，CC基因型为非暴露因素进行分析，结果显示，携带TT+CT基因型人群发生ASD的风险是携带CC基因型的1.67倍（P=0.002）。种族亚组分析显示，高加索人携带TT+CT基因型人群发生ASD的风险是携带CC基因型的1.30倍（P=0.024）；亚洲人携带TT+CT基因型人群发生ASD的风险是携带CC基因型人群的1.99倍（P=0.010）；而在非洲人组中，携带TT+CT基因型人群发生ASD的风险与携带CC基因型无关（OR=2.51，P=0.324），见图3。

2.2.3隐性模型TTvs.CT+CC与ASD发病的Meta分析 见图4。各研究结果存在异质性（I2=61.7%，PFig.3TheMetaanalysisforestmapbetweenMTHFRgeneC677Tpolymorphism(TT+CTvs.CC)andASD图3MTHFR基因C677T多态性（TT+CTvs.CC）与ASD的Meta分析森林图图4MTHFR基因C677T多态性（TTvs.CT+CC）与ASDMeta分析森林图0.001），因此采用随机效应模型进行Meta分析。ASD组、对照组TT基因型分别为280例、269例，ASD组、对照组CT+CC基因型分别为1869例、1982例。以TT基因型为暴露因素，CT+CC基因型为非暴露因素进行分析，结果显示，携带TT基因型人群发生ASD的风险与携带CT+CC基因型无关（OR=1.20，P=0.340）。在种族亚组分析中，3组Meta分析结果均显示携带TT基因型人群发生ASD的风险与携带CT+CC基因型无关（高加索人：OR=1.19，P=0.376；亚洲人：OR=1.46，P=0.309；非洲人：OR=1.69，P=0.618）。

2.2.4

纯合子模型TTvs.CC与ASD发病的Meta分析 各研究结果存在异质性（I2=67.5%，P＜0.001），采用随机效应模型进行Meta分析。ASD组、对照组TT基因型分别为280例、269例，ASD组、对照组CC基因型分别为901例、1102例。以TT基因型为暴露因素，CC基因型为非暴露因素进行分析，结果显示，携带TT基因型人群发生ASD的风险与携带CC基因型人群差异无统计学意义（OR=1.43，P=0.113）。种族亚组分析均显示携带TT基因型人群发生ASD的风险与携带CC基因型人群差异无统计学意义（高加索人：OR=1.36，P=0.205；亚洲人：OR=1.72，P=0.165；

非洲人：OR=2.47，P=0.490），见图5。

敏感性分析通过逐一排除研究的方法进行敏感性分析，Meta分析结果显示各效应量无明显改变，证明结果具有稳定性，见图6。

2.4发表偏倚 通过漏斗图和Egger线性回归法检验发表偏倚，各基因型漏斗图可观察到各独立研究OR在坐标轴上围绕漏斗图的中心线散开，各基因型Egger检验结果均P＞0.05（Tvs.C，t=2.040，P=0.061；TT+CTvs.CC：t=1.923，P=0.076；TTvs.CT+CC，t=1.711，P=0.109；TTvs.CC，t=1.730，P=0.106），提示本研究无明显发表偏倚。见图7。

a：等位基因模型（Tvs.C）；b：显性模型（TT+CTvs.CC）；c：隐性模型（TTvs.CT+CC）；d：纯合子模型（TTvs.CC）Fig.7Funnelplotforpublicationbiastestinthefourgeneticmodels图7各遗传模型的发表偏倚漏斗图

3讨论

化为丙氨酸［27］，显著，降／低MTHFR酶的活性和血浆同型半胱氨酸的浓度，导致神经和血管功能障碍，继而影响大脑功能［28］，而这种改变可能会增加ASD的易感性。

近年来关于MTHFR基因多态性与ASD易感性的研究较多，但结论尚不完全统一，因此笔者收集了相关文献进行Meta分析。本研究共纳入了16篇文献，其中6篇表明MTHFRC677T多态性与ASD发生

献，并进行了敏感性分析及发表偏倚检测，其结果也叶酸主要参与甲基化循环以及DNA和RNA的生物合成，叶酸代谢障碍可导致机体氧化应激易感性，增加ASD的发病风险［23］。MTHFR是单碳代谢过程中叶酸和同型半胱氨酸代谢的关键酶。它将5，10-亚甲基四氢叶酸‘转化为5-甲基四氢叶酸，并参与DNA甲基化相关的叶酸和，同／，型半胱氨酸的转化［24-25］。遗传因素与ASD的发生密切相关，人类MTHFR基因定位于染色体1p36.3［26］，经鉴定MTHFR基因有14个与酶缺陷相关的单核苷酸多态性位点，其中C677T基因位点是研究最多、在临床上最重要的变异之一。C677T的变化是由于胞嘧啶被胸腺嘧啶取代，导致在密码子222位点的缬氨酸转提示分．、析、的、结论是较可靠的。

然而本研究还存在不足之处：（1）文献合并分析后异质性较大，进行种族亚组分析后不能明显降低其异质性。Park等［16］发现男性患病率大于女性，提示异质性来源可能与性别有关，但大部分研究未报道不同性别各基因型的详细数据，无法进行性别亚组分析。（2）纳入的部分研究样本量偏小。

综上所述，MTHFRC677T多态性与ASD易感性具有一定相关性，高加索人仅在显性遗传基因型中MTHFRC677T多态性与ASD发生密切相关，而亚洲人在等位基因模型和显性模型中均具有相关性。

参考文献

［1］YamasueHDomesG.OxytocinandautismspectrumdisordersJ］.CurrTopBehavNeurosci，2018，35：449-465.doi：10.1007/7854201724.

［2］梁颖，李楠，郭园园，等.不同干预模式下孤独症谱系障碍儿童干预效果的比较［J］.天津医药，2016，44（7）：880-

3］LyallK，CroenL，DanielsJ，etal.Thechangingepidemiologyofautismspectrumdisorders［J］.AnnuRevPublicHealth，2017，38：81-102.doi：10.1146/annurev-publhealth-031816-04431＼8.

［4］HiraokaM，Kagaw●aY.Geneticpolymorphismsandfolatestatus［J］.CongenitAnom（Kyoto），2017，57（5）：142-149.doi：10.1111/CC）］进行综合分析，结果显示在等位基因模型和显性模型中MTHFRC677T多态性与ASD发生具有相

［5］WanL，LiYZhangZ，etal.Methylenetetrahydrofolatereductaseandpsychiatricdiseases［J］.TranslPsychiatry，2018，8（1）：242.doi：10.1038/s41398-018-0276-6.

［6MuftinNQJubairSHadiSM.IdentificationofMTHFRgeneticpolymorphisminIraqiautisticchildren［J］.GeneRep，2020，18：100585.doi：10.1016/j.genrep.2019.100585.

［7］JamesSJ，MelnykS，JerniganS，etal.Metabolicendophenotypeandrelatedgenotypesareassociatedwithoxidativestressinchildrenwithautism［J］.AmJMedGenetBNeuropsychiatrGenet，2006，141B（8）：947-956.doi：10.1002/ajmg.b.30366.

［9PaşcaSPDroncaEKaucsárTetal.OnecarbonmetabolismdisturbancesandtheC677TMTHFRgenepolymorphisminchildrenwithautismspectrumdisorders［J］.JCellMolMed，2009，13（10）：4229-4238.doi：10.1111/j.1582-4934.2008.00463.x.

［10dosSantosPALongoDBrandalizeAPetal.MTHFRC677TisnotariskfactorforautismspectrumdisordersinSouthBrazil［J］.PsychiatrGenet，2010，20（4）：187-189.doi：10.1097/YPG.0b013e32833a2220.

［11LiuXSolehdinFCohenILetal.PopulationandfamilybasedstudiesassociatetheMTHFRgenewithidiopathicautisminsimplexfamilies［J］.JAutismDevDisord，2011，41（7）：938-944.doi：10.1007/s10803-010-1120-x.

［12SchmidtRJHansenRLHartialaJetal.Prenatalvitaminsonecarbonmetabolismgenevariants，andriskforautism［J］.Epidemiology，2011，22（4）：476-485.doi：10.1097/EDE.0b013e31821d0e30.

［13GuoTChenHLiuBetal.MethylenetetrahydrofolatereductasepolymorphismsC677TandriskofautismintheChineseHanpopulation［J］.GenetTestMolBiomarkers，2012，16（8）：968-973.doi：10.1089/gtmb.2012.0091.

［14DivyakoluSTejaswiniYThomasWetal.EvaluationofC677Tpolymorphismofthemethylenetetrahydrofolatereductase（MTHFR）geneinvariousneurologicaldisorders［J］.JNeurolDisord，2013，2（1）：142-146.doi：10.4172/2329-6895.1000142.