GPC6基因的生物信息学分析

绪论Glypican家族研究背景硫酸乙酰肝素蛋白多糖家族（Glypicans，GPCs）是一类通过GPI锚定在细胞膜表面的硫酸乙酰肝素蛋白聚糖（HSPGs），其核心蛋白包含保守的Glypican结构域，可结合多种生长因子、细胞因子及细胞外基质成分，参与调控细胞增殖、分化、迁移及信号传导REF_Ref14701\r\h[2]REF_Ref14707\r\h[3]。硫酸乙酰肝素蛋白聚糖家族的结构高度保守，其核心特征包括：核心蛋白、硫酸乙酰肝素（HS）链和通过糖基磷脂酰肌醇（GPI）锚定将其连接至细胞膜外表面REF_Ref7969\r\h[4]。哺乳动物中已鉴定出6个Glypican家族成员（GPC1-GPC6），它们通过结构差异与配体特异性结合，在胚胎发育、组织稳态及疾病发生中发挥多样功能REF_Ref14773\r\h[5]。例如，GPC1在胰腺癌中高表达并促进肿瘤细胞代谢REF_Ref14812\r\h[6]，GPC3是肝细胞癌的重要诊断标志物REF_Ref14887\r\h[7]，而GPC6作为家族中较晚被关注的成员，其生物学功能仍需深入探究。GPC6基因简介GPC6基因是糖蛋白家族（Glypicans）的第六个成员，其编码的Glypican-6是一种糖基磷脂酰肌醇（GPI）锚定的肝素硫酸蛋白聚糖（HSPG），通过细胞表面受体功能调控生长因子信号转导（如Hh、Wnt等）和细胞-基质相互作用。该基因定位于13号染色体13q31.3-q32.1区域，全长约1.19Mb，包含12个外显子，可产生4种转录本及蛋白变体。GPC6蛋白由555个氨基酸构成，其核心蛋白分子量约为60-65kD，其结构包括N端信号肽、14个保守的半胱氨酸残基形成的二硫键结构域，以及C端GPI锚定位点，与GPC4的同源性最高，高达63%REF_Ref8181\r\h[8]。GPC6在胎儿及成人组织中广泛表达，尤其高表达于卵巢、肝脏、肾脏、小肠和结肠，但在胸腺和外周血白细胞中缺失，在胚胎发育中，主要定位于间充质组织，如血管平滑肌和肺间质细胞REF_Ref8181\r\h[8]，其在胎儿胚胎发育过程中的高表达暗示了其在发育调控过程中具有重要的作用。GPC6的作用机制GPC6作为细胞表面核心受体，GPC6通过硫酸肝素链(HS)介导Wnt、Hedgehog（Hh）、FGF、BMP等信号通路，并且调控这些信号通络的传导，从而控制细胞外环境与细胞或者细胞之间的相互作用。调控Hedgehog（Hh）信号通路（1）促进Hh信号传导Hh是生长因子/形态发生素，在胚胎发育过程中的组织模式化和保持成人体内平衡中起关键作用。Shh和Ihh是哺乳动物两种主要的Hh因子REF_Ref8295\r\h[9]，GPC6通过核心蛋白与Hh配体（如Shh、Ihh）结合，并通过HS链与Hh受体Ptc1相互作用，稳定Hh-Ptc1复合物，增强信号传导REF_Ref14949\r\h[10]。初级纤：在无Hh时，GPC6定位于纤毛外；Hh刺激后，GPC6与Ptc1共定位至纤毛膜，促进信号级联REF_Ref14949\r\h[10]。功能缺失：GPC6突变或HS链缺失（ΔGPC6）导致Hh信号活性丧失，引发长骨生长迟缓（如OMOD1）REF_Ref14949\r\h[10]。调控Wnt信号通路经典Wnt通路

GPC6通过HS链与Wnt配体（如Wnt3a）结合，增强其与受体复合物（Frizzled-LRP5/6）的亲和力，促进β-catenin稳定和核转位，激活经典Wnt靶基因,在经典Wnt通路中Wnt配体与Frizzled受体结合后，抑制GSK-3β活性，阻止β-catenin降解，使其进入核内激活靶基因（如细胞增殖相关基因）REF_Ref8527\r\h[11]。可以通过化学激活剂（如氯化锂、BIO、SB216763等）通过抑制GSK-3β活性，阻止β-catenin降解，增强经典Wnt信号，间接上调GPC6表达，影响细胞分化。非经典Wnt通路GPC6通过HS链结合Wnt5a，抑制其在上皮细胞的活性，但促进间质细胞释放Wnt5a，激活非经典通路（如Ca²⁺-CAMKII通路），从而激活JNK（c-JunN-末端激酶）和p38MAPK（丝裂原活化蛋白激酶）,调控细胞迁移和极性REF_Ref15099\r\h[12]。在鼻咽癌中，GPC6通过激活非经典Wnt5a信号，上调JNK和p38MAPK通路，促进癌细胞侵袭REF_Ref15099\r\h[12]。其他信号通路FGF通路：PC6通过HS链结合FGF2，增强其与FGFR的结合，促进MAPK/ERK信号传导，调控细胞增殖和分化。BMP通路：与BMPR2相互作用，参与毛囊干细胞分化REF_Ref15129\r\h[13]。GPC6基因与疾病关联研究进展肺癌肺癌是全球范围内发病率和死亡率极高的恶性肿瘤之一，其发生与发展涉及复杂的分子机制和病理过程。黄志文的研究发现，GPC6基因在肺腺癌细胞中的表达异常，上调GPC6基因的表达可以促进肺腺癌细胞的增殖，而通过敲减GPC6基因则能够抑制肺腺癌细胞的增殖并且诱导其凋亡，这表明GPC6在肺腺癌细胞的增殖过程中发挥着促进作用REF_Ref14342\r\h[1]。崔莹莹在研究GPC1表达对A549细胞的影响时发现，GPC1表达沉默时抑制A549增殖并且显著促进其增殖，研究结果同样提示了糖基磷脂酰肌醇蛋白聚糖家族成员在肺癌发生发展中可能存在的的潜在作用REF_Ref15207\r\h[14]。龚庆豪等人的研究显示，GPC6促进肺腺癌细胞增殖的机制可能与调节细胞周期相关蛋白的表达有关，GPC6的mRNA和蛋白在GIST组织中显著高于癌旁组织，也可能通过影响相关信号通路,如ERK1/2通路的活性，GPC6通过增强ERK1/2磷酸化，调控EMT进程，从而促进癌细胞的侵袭REF_Ref15233\r\h[15]。胃肠道间质瘤胃肠道间质瘤（GIST）是一类起源于胃肠道间叶组织的肿瘤，占胃肠道恶性肿瘤的1%~3%，是最常见的胃肠道间叶源性肿瘤。龚庆豪等人的研究表明，Glypican-6在胃肠道间质瘤组织中呈高表达，并且通过调节ERK1/2通路促进细胞的增殖、迁移及侵袭REF_Ref15233\r\h[15]。在正常生理状态下，ERK1/2通路参与细胞的增殖、生长和分化等多种生物学功能的调节，但在肿瘤细胞中，GPC6的异常高表达会持续激活该通路，导致细胞增殖失控、迁移和侵袭能力增强，从而促进胃肠道间质瘤的发生发展。因此，通过靶向抑制ERK1/2通路磷酸化可能是抗肿瘤的一个新的治疗策略REF_Ref15233\r\h[15]。此外，研究还发现，GPC6的表达水平与胃肠道间质瘤患者的临床病理特征及预后相关，高表达患者的预后较差REF_Ref15233\r\h[15]。鼻咽癌鼻咽癌是一种起源于鼻咽部黏膜上皮的恶性肿瘤，其发病部位位于鼻腔后方、颅底下方的鼻咽腔。利用CEO数据库分析GPC6表达水平，发现GPC6在鼻咽癌组织中的表达显著高于正常组织，GPC6可能通过调控多种信号通路来促进鼻咽癌细胞的增殖、迁移和侵袭REF_Ref15383\r\h[16]。GPC6与Wnt5A结合后可以激活非经典Wnt通路，诱导NFAT转录因子活化，从而促进癌细胞迁移和侵袭REF_Ref15099\r\h[12]，从而在鼻咽癌的发生发展中发挥重要作用。肾透明细胞癌肾透明细胞癌（ccRCC）是肾细胞癌（RCC）中最常见的亚型，约占所有肾细胞癌的70%～80%。研究发现，GPC6在肾透明细胞癌组织中的显著高表达，且与肿瘤的大小、分期及转移有关REF_Ref15449\r\h[17]。张琨通过研究体外实验表明，沉默GPC6基因可以抑制肾透明细胞癌细胞的增殖、迁移和侵袭，并诱导细胞凋亡，而体内实验显示，敲低GPC6基因能够抑制肿瘤的生长，肿瘤的转移能力降低REF_Ref15449\r\h[17]。胰腺癌胰腺癌是一种起源于胰腺导管上皮或腺泡细胞的恶性肿瘤，因其早期诊断困难、进展迅速、预后极差，被称为“癌中之王”。研究指出，GPC6在胰腺癌患者癌细胞中高表达，发挥着促癌作用REF_Ref14701\r\h[2]REF_Ref15553\r\h[18]。一方面，GPC6可以促进胰腺癌细胞的增殖、迁移和侵袭，抑制细胞凋亡；另一方面，GPC6可以通过与相关信号分子相互作用，调节肿瘤微环境，促进肿瘤血管生成和免疫逃逸REF_Ref14701\r\h[2]REF_Ref16703\r\h[19]。GPC6能够通过Kras/MEK/ERK/FOXO3a等信号通路发挥在胰腺癌中的作用REF_Ref15553\r\h[18]，增强胰腺癌细胞的干性，使其更具侵袭性和耐药性，从而促进胰腺癌的发生发展REF_Ref14701\r\h[2]REF_Ref14949\r\h[10]REF_Ref15553\r\h[18]。GPC6基因与其他疾病心脏功能障碍心肌细胞中GPC6过表达，会增加\o"LearnmoreaboutproteinsynthesisfromScienceDirect'sAI-generatedTopicPages"蛋白质合成并诱导肥大和心力衰竭特征基因REF_Ref9455\r\h[20],GPC6通过ERK1/2信号通路促进心肌细胞生长和纤维化，在心肌细胞中，GPC6水平升高诱导心力衰竭特征分子的表达，并通过ERK1/2激活诱导生长REF_Ref9510\r\h[21]。心脏成纤维细胞中的Glypican-6表达受ANGII和BMP4的调节，BMP4可以激活GPC6表达，加剧心脏重塑REF_Ref9510\r\h[21]。同时在压力超负荷诱导的小鼠模型中，GPC6缺失可减轻心肌肥厚和纤维化，提示其作为治疗靶点的潜力REF_Ref9510\r\h[21]。皮肤相关疾病在皮肤组织中也存在GPC6的表达。通过红外光谱成像研究了Glypican-4和Glypican-6在毛发生长周期中的作用，发现GPC6可能参与调节毛囊细胞的增殖和分化，在毛囊的

生发区（germinativezone）和内根鞘（IRS）中持续存在，还在退行期和休止期的基质分化区中检测到高相关性，推测其在毛发的正常生长及皮肤稳态的维持中起到一定作用REF_Ref14773\r\h[4]。虽然目前关于GPC6与常见皮肤疾病（如银屑病、特应性皮炎等）的直接关联研究较少，但基于其在皮肤细胞中的功能，推测其可能与某些皮肤疾病的发生发展相关，值得进一步深入研究REF_Ref14773\r\h[4]。生物信息学分析的必要性生物信息学是生物学、计算机科学、数学和统计学交叉形成的学科，旨在通过计算技术管理和分析生物分子数据，揭示生命活动的规律。其核心功能包括数据管理,存储与整合海量基因组、转录组、蛋白质组数据（如GenBank、UniProt等数据库）；序列分析，比对、组装和注释基因序列，识别功能元件（如启动子、外显子）；功能预测，通过算法预测基因功能、蛋白质结构及相互作用网络；进化研究，构建系统发育树，解析物种分化与适应性等。生物信息学通过整合计算技术与生物学数据，现在已经了成为基因分析中不可或缺的核心手段，随着未来AI的发展，生物信息学将进一步推动基因研究的发展。GPC6基因研究目的与意义GPC6基因在多种疾病的发生发展中扮演着重要角色，尤其是在肿瘤疾病中，其异常表达可以通过多种信号通路影响肿瘤细胞的增殖、迁移、侵袭和凋亡等过程,在非肿瘤疾病，如骨骼发育异常和皮肤相关疾病中，GPC6也参与了关键的生理病理过程。然而目前的研究仍存在一些不足之处，例如GPC6在不同疾病中的具体分子调控网络尚未完全明确，针对GPC6的靶向治疗策略在临床应用中还需要进一步探索和优化,通过对GPC6基因作用机制和功能的研究，以期待能够通过影响其在人体中的作用来治疗疾病。本研究旨在通过实验研究人源GPC6基因，并结合生物信息学方法，系统解析其序列特征、三维结构、进化关系及功能网络。具体目标包括，进行多序列对比，分析不同物种间的保守结构域；构建GPC6蛋白三维模型，并进行可视化分析；预测GPC6蛋白结构域、功能位点及配体结合区域；解析GPC6在肿瘤及正常组织中的表达模式，挖掘其参与的关键信号通路；构建蛋白互作网络，筛选潜在的功能协同蛋白；构建进化树，分析多物种中GPC6蛋白的进化关系。研究结果将为GPC6的功能机制研究提供理论支撑，尤其为肿瘤靶向治疗提供新的候选分子及作用靶点，具有重要的基础研究与临床转化意义。GPC6基因的生物信息学分析实验材料细胞HEK293T细胞（用于RNA提取）试剂TRIzol试剂、逆转录试剂盒生物信息学工具序列分析：BLAST、ClustalOmega、CDD结构预测：SWISS-MODEL、PyMOL表达分析：GEPIA2、UALCAN功能富集：DAVID、STRING本实验采用的glypican6（GPC6）基因和蛋白序列数据主要来源于NCBI数据库（/），由HGNC上传的智人GPC6基因（GeneID:10082），进行生物信息学分析。实验方法GPC6基因克隆（1）RNA提取与cDNA合成在含10%胎牛血清的DMEM（高糖）培养基中培养HEK293T细胞至细胞汇合率达80-90%

,在冰上弃培养基，用预冷PBS清洗细胞3次，加入TRIzol试剂1000μl，刮下细胞，移至离心管，在冰上裂解细胞20min，每管加入氯仿溶液200μl，剧烈震荡15秒，冰上静置10min后，4℃，12000×g离心20min。取离心后上层水相到新的EP管中，加入等体积异丙醇沉淀RNA，颠倒5次。在冰上静置10min，4℃，12000×g离心10分钟，弃去离心后的上清，用75%乙醇洗涤沉淀2次，后用4℃，8000×g离心3分钟，弃去离心后的上清，空气干燥后溶于DEPC水，使用琼脂糖凝胶电泳对提取的RNA质量进行检测。使用BeyoRT™IIcDNA第一链合成试剂盒进行反转录。（2）PCR扩增以cDNA为模板，进行扩增GPC6基因基因序列的扩增，扩增结束后，使用琼脂糖凝胶电泳（1%）进行验证。扩增体系如下：表2-SEQ表2-\*ARABIC1PCR扩增体系Table2-SEQTable2-\*ARABIC1PCRamplificationsystemTagPCRMasterMix25μLDNAtempiate1μLPrimerF(10μmol/L)2μLPrimerR(10ymol/L)2μLddH2O补至50μL表2-SEQ表2-\*ARABIC2PCR扩增程序Table2-SEQTable2-\*ARABIC2PCRamplificationprotocol预变性94℃4min1个循环变性9430s35个循环退火56−7030s延伸721kb/min终延伸72°C10min1个循环4°C4°C表2-SEQ表2-\*ARABIC3引物设计Table2-SEQTable2-\*ARABIC3Primerdesign引物名称序列上游引物5'-CCGCTCGAGGCatgccttcttggatcggggc-3'下游引物5'-CGGGGTACCtctgcacagtctctgcagtgc-3'序列特征分析（1）多序列比对从NCBIProtein数据库中以FASTA格式下载人、小鼠、大鼠的蛋白质序列，使用ClustalOmega(https://www.ebi.ac.uk/Tools/msa/clustalo/)对人、小鼠、大鼠等物种的GPC6基因在氨基酸水平上对该基因编码蛋白的同源序列进行比对，拥有相似的蛋白序列，可能在结构功能方面有着相似，基于多序列比对数据，GPC6基因在不同物种间的同源性及特异性得以明确展示，为后续的进化分析和功能研究方面提供线索。利用ESPript3.0(https://espript.ibcp.fr/ESPript/ESPript/)对多序列比对后的结果进行保守结构域分析，保守结构域是蛋白质序列中经受长期进化选择压力仍保持高度相似性的功能单元，其核心结构与关键活性位点在不同物种间高度稳定，这种进化稳定性往往与其承担的必需生物学功能（如催化中心或配体结合位点）存在深入关联，因此，我们能够通过对已知相同结构域的蛋白家族，对GPC6的蛋白结构和功能进行合理推测。（2）结构域预测使用NCBICDD（/Structure/cdd/wrpsb.cgi），以FASTA格式上传

GPC6的蛋白序列，对GPC6的结构域进行分析，GPC6的功能主要依赖于GPC6的结构域，结构域的突变或表达异常可导致疾病发生，确定结构域的位置有助于了解GPC6蛋白的功能调控机制以及深入研究GPC6作为治疗靶点的可能性。使用PredGPI（http://gpcr.biocomp.unibo.it/predgpi/）软件，输入GPC6蛋白质的序列，预测GPI锚定位点，GPC6蛋白的GPI锚定不仅介导GPC6的膜定位，还影响其与信号通路的动态调控，对GPC6蛋白至关重要，确定GPC6蛋白的GPI锚定位点，对GPC6蛋白功能作用的研究至关重要。三维结构建模（1）同源建模采用SWISS-MODEL在线平台（/

），选择高相似性模板构建GPC6蛋白结构模型，SWISS-MODEL通过同源建模的方法，根据目标蛋白与已知结构蛋白的序列相似性，选择合适的模板构建三维结构模型，建模完成后，通过评估模型的GMQE（GlobalModelQualityEstimation）值和QMEAN（QualitativeModelEnergyAnalysis）值，选择质量较好的模型进行后续分析。GMQE值范围为0-1，越接近1表示模型质量越好；QMEAN值范围为-4-0，越接近0表示模型质量越高。三级结构模型能直观展示蛋白质中活性位点（如酶催化中心）和结构域的三维排布，可用于预测蛋白质与其他分子（如配体、核酸）的结合界面，可以为药物靶点的开发提供线索。（2）结构分析用PyMOL可视化GPC6三级结构，标记GPC6的GPI锚定信号位点、HS链结合位点、二级结构和表面结构，PyMOL能够将蛋白质、核酸、配体等生物大分子的三维结构以多种形式呈现，帮助研究者直观理解原子空间排列、化学键网络及构象动态变化，通过表面渲染可观察蛋白质的静电势分布，而卡通模型能清晰展示二级结构（α螺旋、β折叠）REF_Ref15952\r\h[22]。功能富集与互作网络（1）互作蛋白预测蛋白质的功能协同性通常依赖于其参与的分子互作网络才能发挥最大作用。为系统解析GPC6蛋白的潜在功能模块，本研究基于STRING数据库（/）构建其蛋白-蛋白相互作用（PPI）网络。在STRING数据库中，输入GPC6蛋白的序列信息，选择相应的物种，设置互作置信度阈值，即可获取与GPC6蛋白与其他蛋白之间的相互作用网络，在该模型中，节点（Nodes）代表蛋白质，边（Edges）表示相互作用，通过分析网络中的节点和边，可以挖掘关键节点蛋白和信号通路，进一步揭示GPC6蛋白在细胞内的功能和调控网络。（2）通路分析将互作蛋白导入DAVID（），进行GO（如“细胞外基质组织”）和KEGG（如“Hippo信号通路”）富集分析，再将得到的结果绘制成气泡图进行分析，了解目的基因在特定生物学过程中的作用，分析代谢通路和信号网络，有助于深入GPC6引起疾病的分子机制。（3）PPI网络构建：通过Cytoscape的NetworkAnalyzer工具对GPC6互作网络进行参数化表征，筛选关键蛋白。通过可视化网络，可更加直观的发现与GPC6直接或间接相互作用的蛋白质，揭示其参与的生物学过程，确定GPC6可能参与的通路，为研究其调控机制提供线索。关键节点通常具有高连接度（度中心性）或中介作用（中介中心性），可能是信号传递的核心，关键节点功能异常可能导致网络失衡，引发疾病。进化树构建进化树（又称系统发育树或生命树）是一种以树状分支结构展示生物类群（如物种、基因、种群等）进化关系的图形工具。它通过节点、分支和距离标尺等要素，揭示生物多样性的起源、分化路径及亲缘关系。在NCBI数据库中获取多物种的GPC6基因序列，以FASTA格式保存，利用MEGA7.0软件，构建邻接法(Neighbor-Joining,NJ)进化树，分析不同物种间GPC6基因的进化关系。实验结果分析GPC6基因克隆分析使用琼脂糖凝胶电泳对提取的RNA质量进行检查，得到结果如下图()所示，在图中可以看见清晰的28SrRNA和18SrRNA两条主带，说明RNA未显著降解，适合用于后续实验。图2-SEQ图2-\*ARABIC1琼脂糖凝胶电泳（1%）验证RNA质量Fig2-SEQFig_2-\*ARABIC1Agarosegelelectrophoresis(1%)wasusedtoverifythequalityofRNA使用琼脂糖凝胶电泳（1%）对PCR扩增结果进行验证，得到结果如下图(REF_Ref17927\h图2-1)所示，我们在图中可以清晰的看到目标高亮条带，PCR扩增成功。图2-SEQ图2-\*ARABIC2琼脂糖凝胶电泳（1%）验证PCR扩增结果Fig2-SEQFig_2-\*ARABIC2Agarosegelelectrophoresis(1%)wasusedtoverifytheresultsofPCRamplification.序列特征分析（1）多序列对比根据图(REF_Ref18375\h图2-2)对人、小鼠、大鼠GPC6蛋白的多序列比对和保守结构域的分析，可以得到GPC6蛋白在整体上高度保守，核心结构域（22-548aa）的氨基酸序列一致性达90%以上，表明这个结构域在维持GPC6蛋白功能中起关键作用，可能在进化或功能上具有高度保守性。图2-SEQ图2-\*ARABIC3GPC6蛋白多序列比对和保守结构域分析Fig2-SEQFig_2-\*ARABIC3MultiplesequencealignmentandconserveddomainanalysisoftheGPC6protein（3）结构域预测利用NCBICDD对基因编码蛋白质的氨基酸fasta序列进行分析保守结构域分析，分析结果如下表（REF_Ref18688\h表2-4）所示，分析结果显示GPC6蛋白有一个名称为Glypican的核心结构域，其结构域区间为22-548aa。表2-SEQ表2-\*ARABIC4GPC6结构域位点Table2-SEQTable2-\*ARABIC4GPC6domainsites名称位置电子值Glypican22-5480e+00使用PredGPI预测GPI锚定位点，预测得到GPC6蛋白的GPI锚定位点位置为530aa。三维结构建模（1）同源建模分析采用SWISS-MODEL在线平台（/

），选择高相似性模板构建GPC6蛋白结构模型，得到结果如下（REF_Ref18972\h图2-3），基于全局模型质量评估（GMQE）评分≥0.8的标准筛选最优模型，选择GMQE值达0.82的高置信度模型，得到结构模型如下所示，经过分析，该蛋白三级结构整体上处于无规则卷曲状态，整体呈现出特定且稳定的三维形态，处于折叠后的天然状态。图2-SEQ图2-\*ARABIC4SWISS-MODEL构建GPC6蛋白的三级结构Fig2-SEQFig_2-\*ARABIC4SWISS-MODELwasusedtobuildthetertiarystructureoftheGPC6protein（2）结构分析用PyMOL可视化GPC6蛋白结构，标记GPC6蛋白二级结构，得到结果如下(REF_Ref23149\h图2-4)，经过分析，GPC6蛋白二级结构由α-螺旋和无规则卷曲结构组成；标记的GPC6蛋白的GPI锚定信号位点，得到结果如图(REF_Ref23172\h图2-6)；标记GPC6蛋白表面结构，展示分子间相互作用界面(REF_Ref19269\h图2-7)，对GPC6蛋白进行静电势分析，得到结果如下（REF_Ref19269\h图2-8），蓝色区域带正电，红色区域负电，GPC6蛋白的硫酸乙酰肝素结合位点作为带负电分子可能结合正电区域，推测GPC6蛋白的乙酰肝素结合位点可能处于图中的蓝色区域。图2-SEQ图2-\*ARABIC5PyMOL可视化GPC6蛋白的二级结构Fig2-SEQFig_2-\*ARABIC5PyMOLwasusedtovisualizethesecondarystructureoftheGPC6protein图2-SEQ图2-\*ARABIC6图PyMOL标记GPC6蛋白表面结构Fig2-SEQFig_2-\*ARABIC6PyMOLwasusedtolabelthesurfacestructureoftheGPC6protein图2-SEQ图2-\*ARABIC7PyMOL标记GPC6蛋白GPI锚定信号位点Fig2-SEQFig_2-\*ARABIC7PyMOLwasusedtolabeltheGPI-anchoredsignalsiteoftheGPC6protein图2-SEQ图2-\*ARABIC8PyMOL对GPC6蛋白进行静电势分析Fig2-SEQFig_2-\*ARABIC8PyMOLperformselectrostaticpotentialanalysisontheGPC6protein功能富集与互作网络（1）互作蛋白预测分析利用STRING工具(/)分析GPC6蛋白是否与其他种类蛋白存在相互作用，得到的结果如下（REF_Ref23221\h图2-8）所示。通过分析得到，GPC6蛋白和GPC5、SDC1、SDC3、HSPG2、EXT1、GLCE、HS6ST3、、SDC2、SDC4和GPC4蛋白间相互作用。EXT1、GLCE和LS6ST3蛋白参与硫酸乙酰肝素（HS）链合成与修饰，HS链的结构修饰决定GPC6与配体（如Wnt3a、FGF2）的结合能力及信号传导效率。GPC4和GPC5与GPC6同为Glypican家族成员，通过共享HS链结合生长因子（如Wnt3a、FGF2），形成信号复合物，增强或抑制配体与受体的结合效率，影响干细胞维持或肿瘤发生。SDC1、SDC2、SDC3和SDC4蛋白同属于Syndecan家族（跨膜HS蛋白聚糖），参与细胞黏附、迁移及生长因子信号传导，通过HS链结合细胞外基质分子。HSPG2蛋白是基底膜蛋白聚糖（Perlecan），通过基底膜HS链与GPC6结合，调控细胞迁移和组织稳态。图2-SEQ图2-\*ARABIC9STRING预测的蛋白的相互作用Fig2-SEQFig_2-\*ARABIC9ProteininteractionspredictedbySTRING（2）通路分析利用DAVID工具（/home.jsp），进行进行GO和KEGG通路富集分析，得到结果如下图(REF_Ref23267\h图2-9、REF_Ref23286\h图2-10)，从图中分析可以得到，GO:0098552~sideofmembrane和GO:0005886~plasmamembrane这两条通路和GPC6的功能显著相关，GPC6蛋白可能参与质膜动态调控、突触信号传递或细胞间通讯。分析得到，hsa04820:Cytoskeletoninmusclecells、hsa0524:Glycosaminoglycanbiosynthesis–heparansulfate/heparin、hsa05205:Proteoglycansincancer和hsa04514:Celladhesionmolecules(CAMs)通路与GPC6蛋白显著性最高，GPC6蛋白与肌肉细胞（心肌、骨骼肌）中细胞骨架的结构与功能调控、ECM重塑、细胞黏附和癌症相关通路显著相关。图2-SEQ图2-\*ARABIC10GO通路富集分析气泡图Fig2-SEQFig_2-\*ARABIC10BubblechartofGOpathwayenrichmentanalysis图2-SEQ图2-\*ARABIC11KEGG通路富集分析气泡图Fig2-SEQFig_2-\*ARABIC11BubblechartofKEGGpathwayenrichmentanalysis（3）PPI网络构建采用Cytoscape对GPC6蛋白-蛋白相互作用（PPI）网络进行拓扑分析，根据蛋白间的Degree，节点直径与其Degree值呈正相关，颜色梯度映射采用Jet色谱（蓝色:低度值→红色:高度值），筛选出关键节点蛋白，得到结果如下(REF_Ref23322\h图2-11)，从图中分析得到，GPC6与多个其他蛋白质有直接的相互作用，每个蛋白在相互作用的生物过程中都十分重要。图2-SEQ图2-\*ARABIC12蛋白互作网络可视化Fig2-SEQFig_2-\*ARABIC12Proteininteractionnetworkvisualization进化树构建分析选取一些物种的GPC6基因在NCBI数据库中以以FASTA格式下载，利用MEGA7.0软件，构建邻接法(NJ)进化树,构建结果如下图(REF_Ref23351\h图2-12)所示。在灵长类动物中，人类（Homosapiens）、黑猩猩（Panpaniscus）和大猩猩（Gorillagorillagorilla）的GPC6转录变体（分别为标准mRNA、X2和X1）聚集于同一分支，表明它们具有高度保守的GPC6基因结构；马科动物中，斑马（Equusquagga）与驴（Equusasinus）有着很近的亲缘关系；鲸类的进化地位尤为显著，虎鲸（Orcinusorca）、小须鲸（Balaenopteraacutorostrata）与北大西洋露脊鲸（Eubalaenaglacialis）的GPC6mRNA构成独立分支，与陆生偶蹄类骆驼（Camelusferus）、绵羊（Ovisaries）明显分离。图2-SEQ图2-\*ARABIC13多物种的GPC6基因进化树分析Fig2-SEQFig_2-\*ARABIC13MultispeciesphylogenetictreeanalysisoftheGPC6gene实验讨论本研究通过多序列比对和保守结构域分析，发现GPC6蛋白在人类、小鼠、大鼠等多个物种间整体高度保守，尤其是其Glypican核心结构域（22-548aa）在进化过程中保持了极高的稳定性。这一结果不仅表明GPC6基因在进化上具有重要的保守性，也暗示了该结构域在执行生物功能中可能具有关键的作用。利用生物信息学工具的进行分析，我们对GPC6蛋白的三维结构进行建模分析，结果显示该蛋白三级结构处于无规则卷曲状态，整体呈现出特定且稳定的三维形态，处于折叠后的天然状态。我们对GPC6蛋白的结构域和GPI锚定位点进行预测，发现GPC6蛋白具有典型的Glypican家族特征，拥有Glypican结构域，结构域位置为22-548aa，GPI锚定位点位置为530aa。对GPC6蛋白三维结构进行可视化分析，标记GPC6蛋白的二级结构、表面结构和GPI锚定位点，并且进行GPC6蛋白静电势分析并且将结果可视化，发现该蛋白二级结构主要由α-螺旋和无规则卷曲结构组成，并且更加直观的观察到GPI锚定位点以及硫酸乙酰肝素结合位点可能存在的位置。通过STRING工具分析GPC6蛋白与其他蛋白的相互作用网络，我们发现GPC6与GPC5、SDC1、SDC3等多个蛋白存在相互作用，这些蛋白共同参与硫酸乙酰肝素链的合成与修饰、细胞间通讯等关键生物过程，这些互作关系不仅揭示了GPC6蛋白在细胞内的复杂功能网络，也为其在疾病发生发展中的作用提供了新的线索。利用GPC6蛋白和其他蛋白的互作网络进行GO通路和KEGG通路进行反洗，发现GPC6蛋白可能通过结合细胞外基质成分和调控蛋白酶活性，参与细胞黏附、迁移及生长因子信号传导等生物过程。通过构建GPC6基因的进化树，我们进一步揭示了不同物种间GPC6基因的进化关系。例如，灵长类动物中人类、黑猩猩和大猩猩的GPC6基因高度保守，而马科动物中斑马与驴的亲缘关系较近，这些进化关系的研究有助于我们理解GPC6基因在物种进化中的适应性和功能多样性。实验总结与心得体会文章总结本文主要对GPC6基因进行了全面的生物信息学分析，内容涵盖了基因克隆、进化关系、表达模式及潜在功能等多个方面。研究利用了多种生物信息学工具和方法，包括多序列比对、结构域预测、进化树构建、三维结构预测、功能位点预测、互作网络分析和表达模式分析等。通过多序列比对和结构域预测，研究发现GPC6蛋白具有典型的Glypican家族特征，在多个物种间具有高度保守的Glypican核心结构域，表明其在进化过程中保持了功能的稳定性。三维结构预测结果显示，GPC6蛋白具有特定的三维形态，且其二级结构由α-螺旋和无规则卷曲组成。功能富集和互作网络分析揭示了GPC6蛋白与其他多个蛋白的相互作用关系，这些蛋白共同参与硫酸乙酰肝素链的合成与修饰、细胞间通讯等关键生物过程。进化树构建则揭示了不同物种间GPC6基因的进化关系，为理解其进化历史提供了线索。综上所述，本文通过对GPC6基因的生物信息学分析，揭示了其在结构、进化、功能和表达模式等方面的多个重要特征，为进一步深入研究其功能及在相关疾病中的作用机制提供了理论依据。心得体会在参与GPC6基因生物信息学分析的过程中，我深刻体会到了生物信息学在基因功能研究中的重要作用。通过多序列比对、三维结构建模、互作网络构建和功能富集以及进化树构建等生物信息学分析工具，我们能够快速获取基因的结构特征和进化关系，为后续的功能研究提供重要线索。通过这次研究，我不仅加深了对生物信息学工具使用的理解，同时更了解生物信息学工具功能便捷强大，对我的研究提供了很大的帮助，同时通过这次对GPC6基因进行系统的生物信息学分析，也对GPC6基因的功能和研究前景有了更全面的认识，了解到GPC6基因在疾病产生过程中的作用，特别是在癌症方面，期待GPC6基因的研究在未来能够有更大的进展，为癌症的提供帮助。同时我也认识到本次实验还有很大不足，对实验结果的认识较为浅薄，分析较为浅显，如果以后还有机会进行相关的研究，希望能够更加深入进行分析。参考文献[1]黄志文.Glypican6基因对肺腺癌细胞增殖的影响及其机制的初步研究[D].苏州大学.[2]姚艳丽。多糖HH1-1抗胰腺癌和糖基磷脂酰肌醇蛋白聚糖GPC6促胰腺癌的功能及机制研究[D].中国科学院大学，2019.[3]葛静，薛丹，王军.Glypican-3和Glypican-6在巨大儿胎盘组织中的表达及意义[J].中国优生与遗传杂志，2013,21(01):19-21.SukhneerajP.Kaur,BrianS.Cummings,Roleofglypicansinregulationofthetumormicroenvironmentandcancerprogression,BiochemicalPharmacology,Volume168,2019,Pages108-118,ISSN0006-2952,/10.1016/j.bcp.2019.06.020.[4]Colin-PierreC,UntereinerV,SockalingumGD,etal.InvestigationofGlypican-4and-6byInfraredSpectralImagingduringtheHairGrowthCycle.IntJMolSci.2023Mar3;24(5):4528.[5]HanX,LiY,ZhangY,etal.Glypican-1promotespancreaticcancerprogressionbyenhancingWarburgeffect.CellDeathDis.2018Jan25;9(2):189.[6]ZhouY,ChenX,LiX,etal.Glypican-3:apromisingbiomarkerforhepatocellularcarcinomadiagnosisandtherapy.ExpertRevMolDiagn.2019Mar;19(3):247-256.MarkVeugelers,BartDeCat,HelgaCeulemans,Anne-MarieBruystens,ChristienCoomans,JoachimDürr,JorisVermeesch,PeterMarynen,GuidoDavid,Glypican-6,aNewMemberoftheGlypicanFamilyofCellSurfaceHeparanSulfateProteoglycans*,JournalofBiologicalChemistry,Volume274,Issue38,1999,Pages26968-26977,ISSN0021-9258,/10.1074/jbc.274.38.26968.WenShi,TomoyukiKaneiwa,MarzenaCydzik,JeanGariepy,JorgeFilmus,Glypican-6stimulatesintestinalelongationbysimultaneouslyregulatingHedgehogandnon-canonicalWntsignaling,MatrixBiology,Volume88,2020,Pages19-32,ISSN0945-053X,/10.1016/j.matbio.2019.11.002.[7]MarianaCapurro,

TomomiIzumikawa,

PhilippeSuarez,

WenShi,

MarzenaCydzik,

TomoyukiKaneiwa,

JeanGariepy,

LuisaBonafe,

JorgeFilmus;Glypican-6promotesthegrowthofdevelopinglongbonesbystimulatingHedgehogsignaling.

JCellBiol

4September2017;216(9):2911–2926.石鹏双,房静远.Wnt信号通路与胃癌发生、发展关系的研究进展[J].胃肠病学,2021,26(06):353-357