马凡综合征关联药物及相关基因

王天毅　牛兆倬　黄强　生伟　吴建涛　李丹丹

[摘要]目的 探讨某些药物以及所涉及的基因引起马凡综合征的可能机制。

方法通过比较毒理基因组学数据库（CTD）分析引起马凡综合征关联的药物及其相互作用的基因。选取药物与马凡综合征高度相关的基因进行基因本体分析和KEGG信号通路等基因富集分析。

结果有24种关键基因/蛋白参与凋亡过程的正向调节、血压的负向调节、细胞对有机氮的反应、转化生长因子β信号、多细胞的生物过程、肾脏发育、脂多糖反应等生物学进程，以及与蛋白多糖反应、肾素分泌、激活丝裂原活化蛋白激酶信号通路，其中与基因有关的药物有阿替洛尔、多四环素、氯沙坦、硫代巴比妥酸。

结论药物可能通过影响基因表达从而引起某些信号通路的改变等内在分子生物学机制与马凡综合征产生关联，通过分析引起马凡综合征疾病相关基因/蛋白通路，可为临床应用新药、生物疗法和手术疗法提供科学依据。

[关键词]马凡综合征;分子生物学;药物基因组学检测;比较基因学

[中图分类号]R682

[文献标志码]A

[文章编号]2096-5532（2020）01-0035-05

马凡综合征（MFS）是最常见的常染色体显性遗传性结缔组织疾病之一[1]，其发病机制不明[2]。目前，国内相关研究报道MFS发病率為0.04%～0.10%[3]，其临床表现主要累及心血管系统、骨骼系统、眼和神经系统等[4]。该疾病一旦累及心血管系统[5]，则病人极有可能死于心脏衰竭或主动脉瘤破裂[6-7]。为更深入了解MFS遗传学特征，本研究通过利用专业数据库查阅相关文献，尝试分析某些药物以及基因/蛋白引起MFS的可能性及其可能机制。现将结果报告如下。

1资料和方法

1.1基因筛选

以MFS（MESH：D008382）作为目标疾病，应用比较毒理基因组学数据库（CTD）[8]与STRING数据库筛选与MFS关联的药物及基因的数据集。

1.2通路和生物进程功能富集分析

分别利用http：//string-db.org/中的string[9]和http：//metascape.org中的metascape[10]对关键基因进行注释和进一步分析。参考的数据资源包括KEGG通路、GO生物学进程、Reactome基因组、经典通路和哺乳动物蛋白复合物数据库（CORUM）。

根据相关文献方法[11-13]进行通路和生物进程功能富集分析。

1.3蛋白-蛋白互作富集分析

应用MetascapeBioGrid[14]、InWeb_IM以及OmniPath[15]数据库对蛋白-蛋白互作进行富集分析。应用MCODE算法[16]对互作蛋白进行模块化分析，每个模块的蛋白进一步进行通路和生物进程富集分析。根据P值确定最佳前3位的通路或生物进程作为所在模块的功能描述。

2结果

2.1与MFS相关的基因/蛋白

以MFS作为目标疾病，通过CTD数据库检索显示，有258种基因与MFS有关，其中基因评价分值≥5.5有26种基因，见表1。其中GJA1基因与MFS联系最为密切，其次为ABCB11基因。

2.2MFS相关基因/蛋白对应药物

CTD数据库检索显示，1546种药物与MFS有关。其中基因评价分值≥5.5的26种基因对应的药物见表2，有阿替洛尔、多四环素、氯沙坦、硫代巴比妥酸等4种药物与26种基因密切相关。其中氯沙坦与26种基因相关程度达84.62%。

2.326种基因/蛋白相互作用关系

STRING数据库分析结果显示，有24种基因在STRING数据库识别（除SNORA49、PINK-1）（图1）。24种基因之间存在着密不可分的互作关系，基因与基因、蛋白与蛋白的相互作用可以通过基因数据库和实验来验证，还可能通过基因间相邻作用、基因间相互结合与基因间相互同现而发生相互联系，一些基因转录文本、基因间共表达及基因同源性都可以导致基因与基因、蛋白与蛋白相互作用。

2.424种关键基因/蛋白富集分析

应用metascape进一步对24种关键基因富集分析，将富集后的结果按功能分簇，结果见表3。24种关键基因主要集中在循环系统、超氧化物代谢过程、凋亡过程的正向调节、血压的负向调节、细胞对有机氮的反应、转化生长因子β信号、多细胞的生物过程、肾脏发育、脂多糖反应、钠离子运输过程以及有机羟基化合物转运等GO生物学进程，与蛋白多糖反应、乙型肝炎、肾素分泌以及激活丝裂原活化蛋白激酶信号通路等信号通路相关。

2.5蛋白质互作网络模块

应用Cytoscape中MCODE模块化分析结果显示，24种基因/蛋白共有2个基因模块（MCODE1和MCODE2）（图2）。NOS2、REN、AGT、MMP9、MYC、CREB1、EDN1与ACTA2、CASP3、MMP2分别通过结合、同现、共表达以及同源性相互作用。MCODE1中7种基因/蛋白（NOS2、REN、AGT、CREB1、EDN1、MMP9、MYC）参与相关生物过程，包括转化生长因子β信号生物过程（CREB1、MYC以及EDN1）、激活丝裂原活化蛋白激酶信号通路（AGT、CREB1、MYC、EDN1、REN）、NO正向调节（AGT、EDN1、NOS2）及纤维蛋白原正向调节通路（AGT、MYC）;而模块MCODE2中3种基因/蛋白（ACTA2、CASP3、MMP2）不参与相关生物过程。

3讨论

本文通过对MFS与人类相关的24种主要基因/蛋白富集显示，MFS相关基因/蛋白主要集中在循环系统、超氧化物代谢过程、凋亡过程的正向调节、血压的负向调节、细胞对有机氮的反应、转化生长因子β信号、多细胞的生物过程、肾脏发育、脂多糖反应、钠离子运输过程以及有机羟基化合物转运、蛋白多糖反应、乙型肝炎、肾素分泌以及激活丝裂原活化蛋白激酶信号通路等信号通路。其中转化生长因子超家族信号通路对维持细胞外基质和组织稳态起着重要的作用[17-18]。

国内外相关学者研究认为，原纤维蛋白-1基因（FBN1）是MFS致病的主要基因[19-21]。FBN1是一种抗胶原蛋白酶的非硫酸化的糖蛋白，参与蛋白多糖反应生物过程，为细胞外基质微纤维的重要组成成分[22];同时还可以调节转化生长因子β（TGF-β）的释放[23-25]。有研究表明，FBN1突变会引起结构改变，不仅可以导致弹性纤维组装障碍，而且还会引起TGF-β的过度激活，从而引起TGF-β介导的信号传导增强[26]。KIELTY等[27]研究发现，由于基因突变引起的微纤维蛋白异常可能导致血管弹性纤维改变，从而使病人更容易发生主动脉夹层或形成主动脉瘤[28-29]。

TGF-β的过度激活与细胞外基质中MMPs的异常表达有关，MMPs作为一种锌内肽酶，其中的MMP2和MMP9是经典TGF-β通路的下游调控产物[30]。CHUNG等[31]研究显示，TGF-β针对性诱导MMP2和MMP9高表达，在实验诱导以及自然形成的动脉瘤组织中都可以观察到MMPs高表达[24]，而MMPs高表达使得主动脉相关蛋白质分解从而诱导动脉瘤形成[32-33]。DIETZ[34]对MFS小鼠模型研究发现，主动脉壁组织中确实存在TGF-β信号转导通路过度激活的现象，使用TGF-β信号转导通路阻断剂氯沙坦，可阻断MFS小鼠主动脉病变的进展;同样，FORTEZA等[7]研究显示，血管紧张素Ⅱ受体阻滞剂氯沙坦能抑制TGF-β介导的细胞外信号调节激酶激活，进而延缓甚至逆转主动脉扩张等病变的进展，从而为MFS血管病变的治疗提供有力临床证据[37]。本文通过CTD数据库初步分析显示，氯沙坦与26种基因密切相关，因此，氯沙坦可应用于治疗MFS病人。

本文通过基因/蛋白富集分析显示，TGF-β还可激活其他信号传导通路，该通路主要为丝裂原活化的蛋白激酶（MAPK）信号转导通路，涉及AGT、CREB1、MYC、EDN1、REN等相关基因，这与国内外报道的结果相符[36-39]。这些基因通过结合、同现、共表达以及同源性相互作用共同作用于MAPK信号转导通路[40]，该通路上还有细胞外信号调节激酶（ERK）、Jun氨基端激酶（JNK）和P38。在MFS血管组织病变过程中主要是ERK起决定作用[41]，这可能为治疗MFS血管病变提供新的靶点。

同时我们的研究结果也显示，环境因素一氧化氮（NO）也可能对MFS血管病变有一定的影响。这与CHUNG等[42]应用MFS小鼠模型的研究结果相似，MFS与内皮NO合成减少有关，并推测血管内皮功能紊乱可能是MFS病人主动脉并发症发生机制之一。

综上所述，MFS的发病机制目前尚不明确，而国内外的研究多集中在实验动物模型，但仅凭动物研究结果并不能说明其在人类中的机制。本文借助CTD数据库和其他在线分析工具研究显示，药物可能通过影响基因表达从而引起某些信号通路的改变等内在分子生物学机制与MFS产生关联，而该种关联尚需要在人群中进一步验证。

[参考文献]

[1]武玉敏，崔建，阮银镜，等.一家四代马凡综合征6例报告[J].实用医学杂志，2005，21（9）：999.

[2]冯彤.马凡氏综合征心血管病变的临床与影像学表现[J].影像研究与医学应用，2018，2（4）：64-65.

[3]那开宪.马凡氏综合征[J].首都食品与医药，2017（15）：58-59.

[4]贺晶，高凌根.马凡综合征的分子遺传学研究进展[J].中华保健医学杂志，2015，17（5）：421-424.

[5]苏育琳，陈子怡.马凡氏综合征继发脑血管疾病的临床分析[J].中国神经精神疾病杂志，2019，45（1）：47-50.

[6]覃勇民，伍广伟.马凡综合征的临床研究进展[J].疑难病杂志，2010，9（10）：797-801.

[7]FORTEZAA，EVANGELISTAA，SNCHEZV，etal.EfficacyoflosartanvsatenololforthepreventionofaorticdilationinMarfanSyndrome：arandomizedclinicaltrial[J].EuropeanHeartJournal，2015，37（12）：978-985.

[8]DAVISAP，GRONDINCJ，JOHNSONRJ，etal.Thecomparativetoxicogenomicsdatabase：update2019[J].NucleicAcidsResearch，2019，47（D1）：D948-D954.

[9]SZKLARCZYKD，MORRISJH，COOKH，etal.TheSTRINGdatabasein2017：quality-controlledprotein-proteinassociationnetworks，madebroadlyaccessible[J].NucleicAcidsResearch，2017，45（D1）：D362-D368.

[10]ZHOUYY，ZHOUB，PACHEL，etal.Metascapeprovidesabiologist-orientedresourcefortheanalysisofsystems-leveldatasets[J].NatureCommunications，2019，10（1）：1523-1531.

[11]JERROLDZH.Biostatisticalanalysis[M].India：PearsonEducationIndia，1999.

[12]HOCHBERGY，BENJAMINIY.Morepowerfulproceduresformultiplesignificancetesting[J].StatisticsinMedicine，1990，9（7）：811-818.

[13]COHENJ.Acoefficientofagreementfornominalscales[J].EducationalandPsychologicalMeasurement，1960，20（1）：37-46.

[14]STARKC，BREITKREUTZBJ，REGULYT，etal.BioGRID：ageneralrepositoryforinteractiondatasets[J].NucleicAcidsResearch，2006，34（Databaseissue）：D535-D539.

[15]LITB，WERNERSSONR，HANSENRB，etal.Ascoredhumanprotein-proteininteractionnetworktocatalyzegenomicinterpretation[J].NatureMethods，2017，14（1）：61-64.

[16]BADERGD，HOGUECW.Anautomatedmethodforfin-

dingmolecularcomplexesinlargeproteininteractionnetworks[J].BMCBioinformatics，2003，4（1）：2-10.

[17]BERTOLI-AVELLAAM，GILLISEM，MORISAKIH，etal.MutationsinaTGFβligand，TGFB3，causesyndromicaorticaneurysmsanddissections[J].JournaloftheAmericanCollegeofCardiology，2015，65（3）：1324-1336.

[18]PEPEG，GIUSTIB，STICCHIE，etal.Marfansyndrome：currentperspectives[J].ApplClinGenet，2016，9：55-65.

[19]VERSTRAETENA，ALAERTSM，VANLAERL，etal.Marfansyndromeandrelateddisorders：25yearsofgenediscovery[J].HumanMutation，2016，37（6）：524-531.

[20]张琳，张莹，范鹏，等.FBN1基因钙结合转化生长因子结构域终止突变基因型表型研究[J].中国分子心脏病学杂志，2018，18（3）：2513-2516.

[21]FRANKENR，GROENINKM，DEWAARDV，etal.GenotypeimpactssurvivalinMarfansyndrome[J].EuropeanHeartJournal，2015，36（1）：523-527.

[22]GUGLIELMINAP，GIUSTIB，STICCHIE，etal.Marfansyndrome：currentperspectives[J].TheApplicationofClinicalGenetics，2016，9：55-65.

[23]張琳，张莹，范鹏，等.2型马凡综合征TGFBR2基因新致病突变研究[J].中国心血管杂志，2018，23（4）：288-293.

[24]公兵，杨秀滨.马凡氏综合征血管病变的基础研究进展[J].心脏杂志，2012，24（6）：760-762.

[25]CANNAERTSE，BERTOLI-AVELLAA，GILLISE，etal.MutationsinaTGFβligand，TGFB3，causesyndromicaorticaneurysmsanddissections[J].JournaloftheAmericanCollegeofCardiology，2015，65（13）：1324-1336.

[26]SENGLEG，TSUTSUIK，KEENEDR，etal.Microenvironmentalregulationbyfibrillin-1[J].PLoSGenetics，2012，8（1）：e1002425.

[27]KIELTYCM，SHERRATTMJ，MARSONA，etal.Fibrillinmicrofibrils[J].AdvancesinProteinChemistry，2005，70：405-436.

[28]FAIVREL，MASUREL-PAULETA，COLLOD-BEROUDG，etal.Clinicalandmolecularstudyof320childrenwithmarfansyndromeandrelatedtypeⅠfibrillinopathiesinaseriesof1009probandswithpathogenicFBN1mutations[J].Pediatrics，2009，123（1）：391-398.

[29]LI-WAN-POA，LOEYSB，FARNDONP，etal.PreventingtheaorticcomplicationsofMarfansyndrome：acase-exampleoftranslationalgenomicmedicine[J].BritishJournalofClinicalPharmacology，2011，72（1）：6-17.

[30]刘鼎乾，李军，赖颢.TGF-β信号通路与马方综合征及Loeyz-Dietz综合征研究进展[J].中国分子心脏病学杂志，2018，18（1）：2386-2389.

[31]CHUNGAW，YANGHH，RADOMSKIMW，etal.Long-termdoxycyclineismoreeffectivethanatenololtopreventthoracicaorticaneurysminmarfansyndromethroughtheinhibitionofmatrixmetalloproteinase-2and-9[J].CirculationResearch，2008，102（8）：e73-e85.

[32]WUD，SHENYH，RUSSELLL，etal.Molecularmechanismsofthoracicaorticdissection[J].JournalofSurgicalResearch，2013，184（2）：907-924.

[33]FRANKENR，DENHARTOGAW，RADONICT，etal.BeneficialoutcomeoflosartantherapydependsontypeofFBN1mutationinmarfansyndrome[J].Circulation-CardiovascularGenetics，2015，8（2）：383-388.

[34]DIETZHC.TGF-betainthepathogenesisandpreventionofdisease：amatterofaneurysmicproportions[J].TheJournalofClinicalInvestigation，2010，120（2）：403-407.

[35]BROOKEBS，HABASHIJP，JUDGEDP，etal.AngiotensinⅡblockadeandaortic-rootdilationinMarfanssyndrome[J].NewEnglandJournalofMedicine，2008，358（26）：2787-2795.

[36]LOEYSB，DEPAEPEA.Newinsightsinthepathogenesisofaorticaneurysms[J].Verhandelingen-KoninklijkeAcademievoorGeneeskundeVanBelgie，2008，70（2）：69-84.

[37]NGCM，CHENGAL，MYERSLA，etal.TGF-beta-dependentpathogenesisofmitralvalveprolapseinamousemo-

delofMarfansyndrome[J].TheJournalofClinicalInvestigation，2004，114（11）：1586-1592.

[38]HABASHIJP，DOYLEJJ，HOLMTM，etal.AngiotensinⅡtype2receptorsignalingattenuatesaorticaneurysminmicethroughERKantagonism[J].Science（NewYork，N.Y.），2011，332（627）：361-365.

[39]LEMAIRER，FARINAG，BAYLEJ，etal.AntagonisticeffectofthematricellularsignalingproteinCCN3onTGF-beta-andWnt-mediatedfibrillinogenesisinsystemicsclerosisandMarfansyndrome[J].TheJournalofInvestigativeDermatology，2010，130（6）：1514-1523.

[40]LINCY，CHINCH，WUHH，etal.Hubba：hubobjectsanalyzer-aframeworkofinteractomehubsidentificationfornetworkbiology[J].NucleicAcidsResearch，2008，36（S）：W438-W443.

[41]FRANKENR，DENHARTOGAW，DEWAARDV，etal.Circulatingtransforminggrowthfactor-βasaprognosticbiomarkerinMarfansyndrome[J].InternationalJournalofCar-

diology，2013，168（3）：2441-2446.

[42]CHUNGAW，AUYEUNGK，CORTESSF，etal.EndothelialdysfunctionandcompromisedeNOS/AktsignalinginthethoracicaortaduringtheprogressionofMarfansyndrome[J].BritishJournalofPharmacology，2007，150（8）：1075-1083.

（本文編辑 黄建乡）