YIKE安徽医科大学学报Acta Universitatis Medicinalis AnhuiActa Universitatis Medicinalis Anhui安徽医科大学学报1000-149234-1065/RAnhui Lianzhong Printing Limited CompanyEditorial Board of Acta Universitatis Medi-cinalis Anhui Meishan Road , Hefei 230032《安徽医科大学学报》编辑部安徽省合肥市安徽医科大学校内老图书馆三楼V271-李鹏-综合生物信息学1000–1492（2026）05–0861–1210.19405/j.cnki.issn1000–1492.2026.05 009R 77A基础医学研究综合生物信息学分析及实验验证糖尿病视网膜病变中血管新生相关基因的作用Integrated bioinformatics analysis and experimental validation of angiogenesis-related genes in diabetic retinopathy李鹏LiPeng1梁坤LiangKun2吴峰WuFeng1李佳LiJia3刘伦LiuLun4陶玉林TaoYulin2安徽医科大学附属阜阳医院，眼科，阜阳　236000Department of Ophthalmology，安徽医科大学附属阜阳医院，神经外科，阜阳　236000Department of Neurosurgery， Fuyang Hospital Affiliated to Anhui Medical University， Fuyang　236000安徽医科大学 第二附属医院眼科，合肥　230601Department of Ophthalmology， The Second Affiliated Hospital of Anhui Medical University， Hefei　230601安徽医科大学第一附属医院眼科，合肥　230022Department of Ophthalmology， The First Affiliated Hospital of Anhui Medical University， Hefei　230022

李鹏，男，硕士研究生，主治医师

陶玉林，男，博士，通信作者，E-mail：yelinto@163.com

刘伦，男，博士，硕士生导师，主任医师，通信作者，E-mail：13956956362@163.com

Tao Yulin， E-mail： yelinto@163.comLiu Lun， E-mail： 13956956362@163.com16032026 06022026 23052026615861871861-871目的

探讨糖尿病视网膜病变（DR）发生发展过程中血管新生相关分子机制。

方法

从Gencard网站获取血管新生相关基因，并与DR数据集（GSE60436和GSE94019）的差异表达基因取交集，随后通过功能富集、蛋白质-蛋白质相互作用（PPI）网络筛选候选基因并评估其诊断价值。基因集富集分析（GSEA）探讨候选基因潜在通路，免疫浸润分析揭示候选基因与免疫细胞的相关性。细胞实验验证纤维连接蛋白1（FN1）在高糖处理下的人视网膜微血管内皮细胞（HRMEC）中的作用。

结果

共筛选到237个血管新生相关差异基因，富集于磷脂酰肌醇3-激酶/蛋白激酶B信号通路（PI3K-Akt）、肿瘤抑制蛋白53（P53）、肿瘤坏死因子（TNF）和Janus激酶/信号转导及转录激活因子信号通路（JAK-STAT）等通路。其中胶原蛋白Ⅰ型α1链（COL1A1）、COL1A2、FN1、TNF和肿瘤抑制蛋白53（TP53）为关键基因，且具有较高诊断价值。GSEA提示这些基因参与P53等多条信号通路。CIBERSORTx分析发现其与多种免疫细胞浸润显著相关。高糖处理导致FN1表达升高。在高糖诱导的HRMEC中，与对照si-NC组相比，si-FN1组HRMEC的增殖、迁移及管腔形成能力均显著下降，同时P53蛋白表达升高。

结论

FN1在DR血管新生中具有重要作用，可能成为潜在的诊断和治疗靶点。

Objective

To investigate the molecular mechanisms related to angiogenesis during the development and progression of diabetic retinopathy （DR）.

Methods

Angiogenesis-related genes were obtained from the Gencard website and intersected with differentially expressed genes from DR datasets （GSE60436 and GSE94019）. Functional enrichment and protein-protein interaction （PPI） networks were then used to screen candidate genes and evaluate their diagnostic value. Gene set enrichment analysis （GSEA） was used to explore potential pathways underlying candidate genes， and immune infiltration analysis revealed associations between candidate genes and immune cells. Cellular experiments were conducted to validate the role of fibronectin 1 （FN1） in human retinal microvascular endothelial cells （HRMECs） under high glucose （HG） conditions.

Results

A total of 237 differentially expressed genes related to angiogenesis were identified， enriched in pathways such as phosphoinositide 3-kinase/protein kinase B signaling pathway （PI3K-Akt）， tumor suppressor protein 53 （P53）， tumor necrosis factor （TNF）， and Janus kinase （JAK）/signal transducer and activator of transcription signaling pathway （STAT）. Among them， collagen type I alpha 1 chain （COL1A1）， COL1A2， FN1， TNF， and tumor protein p53 （TP53） were key genes with high diagnostic value. GSEA indicated that these genes were involved in multiple signaling pathways， including P53. CIBERSORTx analysis revealed significant associations with the infiltration of multiple immune cells. HG treatment led to the upregulation of FN1. In HG-induced HRMECs， compared with the si-NC control group， si-FN1 significantly reduced cell proliferation， migration， and tube formation， while P53 protein expression was increased.

Conclusion

This study reveals the important role of FN1 in angiogenesis in DR and suggests that it may be a potential diagnostic and therapeutic target.

糖尿病视网膜病变血管生成FN1P53信号通路管腔形成diabetic retinopathyangiogenesisFN1P53 signaling pathwaytube formation安徽省自然科学基金项目2508085QH303安徽省自然科学基金项目（编号：2508085QH303）Natural Science Foundation of Anhui Province2508085QH303Natural Science Foundation of Anhui Province （No. 2508085QH303）version1.0.0.25090structure-time2026-06-30T13:59:53word-sourceFX

糖尿病视网膜病变 （diabetic retinopathy，DR） 是糖尿病最常见的慢性微血管并发症之一［1］。持续的高血糖、高血压和高血脂，以及吸烟和饮酒等生活方式可导致DR的发生和发展。目前，DR的治疗方法主要包括视网膜激光光凝术、玻璃体内注射血管内皮生长因子 （vascular endothelial growth factor，VEGF） 抑制剂以及玻璃体切除术等。研究［2］表明，作为VEGF下游的功能分子，纤维连接蛋白1（fibronectin 1，FN1）在病理性视网膜新生血管形成中发挥重要作用，能够促进内皮-间质转化并推动异常血管生成。FN1还被鉴定为原发性开角型青光眼（primary open-angle glaucoma，POAG）的关键基因，与免疫相关通路、氧化应激和内质网应激有关［3］。除了血管和神经炎症机制外，年龄相关性黄斑变性（age-related macular degeneration，AMD）的有关研究［4］表明，FN1在循环成纤维细胞和骨髓来源的巨噬细胞中表达显著上调，并在视网膜下纤维化诱导后进一步升高。该研究旨在利用整合生物信息学和实验验证方法，识别血管生成相关基因，重点分析FN1在DR血管生成中的调控作用。

材料与方法材料细胞株来源

人视网膜微血管内皮细胞（human retinal microvascular endothelial cells，HRMEC）来源于江苏海星生物科技有限公司（货号：PAHX-C115）。在HRMEC完全培养基中，于 37 ℃、5% CO₂ 培养条件下维持培养。

抗体和试剂

HRMEC完全培养基（货号：CM-H130）购自武汉普诺赛生命科技有限公司；葡萄糖（货号：HY-B0389）购自美国MedChemExpress有限责任公司；CCK-8 细胞活力检测试剂盒（货号：ml095229）和TRIzol 试剂（货号：ml095476）购自上海酶联生物科技有限公司；LipofectamineTM3000转染试剂（货号：L3000008）由美国Thermo Fisher Scientific公司提供；胎牛血清（FBS，货号：10099-141）购自美国Gibco公司；PrimeScriptTM RT reagent Kit及SYBR® Premix Ex TaqTM Ⅱ 购自大连Takara公司；Transwell 小室（8 μm 孔径，货号：3422）购自美国Corning 公司；Matrigel 基质胶（货号：C0372-5 mL）、结晶紫染色液（货号：C0121-100 mL）、RIPA 裂解液（货号：P0013B）、BCA蛋白定量试剂盒（货号：P0010）、ECL化学发光检测试剂盒（BeyoECL Plus，货号：P0018S）和PVDF膜（货号：FFP33）均购自上海碧云天生物技术有限公司；FN1抗体（兔多克隆，货号：15613-1-AP，1∶2 000）、P53抗体（兔多克隆，货号：10442-1-AP，1∶5 000）、GAPDH抗体（兔多克隆，货号：10494-1-AP，1∶5 000）以及HRP-标记山羊抗兔 IgG（H+L）（货号：SA00001-2，1∶5 000）购自武汉三鹰生物技术有限公司。

主要仪器

CO₂培养箱（型号：Heracell Vios 160i CR，美国Thermo Scientific公司）；倒置显微镜（型号：IXploreTM IX85，日本Olympus公司）；酶标仪（型号：Bio-Rad iMark，美国Microplate Absorbance Reader公司）；实时荧光定量 PCR 仪（ABI 7500 Real-Time PCR System，美国 Applied Biosystems公司）；电泳及转膜系统（型号：Bio-Rad Mini-PROTEAN，美国Tetra Cell公司）；凝胶成像系统（型号：Gel Doc XR+ Gel Documentation System，美国Bio-Rad公司）。

方法DR数据库的分析和差异表达基因（differentially expressed genes，DEGs）的鉴定

从基因表达综合数据库（gene expression omnibus，GEO；https：//www.ncbi.nlm.nih.gov/geo/）下载2个与DR相关数据集，分别是GSE60436（6个DR样本，3个对照样本）和GSE94019（9个DR样本，4个对照样本）。利用R语言（版本号 4.3.2）的“limma”包对2个数据集进行DEGs分析。其中，GSE60436 数据集的筛选标准为 |FC|≥1.3，GSE94019 的筛选标准为|FC|≥1.5，且都满足P<0.05， 结果均通过“ggplot”包进行可视化。

DR相关基因的功能富集分析以及蛋白质-蛋白质相互作用（protein-protein interaction，PPI）网络构建

在Genecard网站（https：//www.genecards.org）检索了5 928个“血管新生相关基因”，基于相关性评分>1的标准筛选得到1 861个候选基因。使用在线工具对1 861个血管新生相关基因以及GSE60436和GSE94019中上调表达基因和下调表达基因进行交集分析。随后DAVID数据库（https：//david.ncifcrf.gov/）对得到的交集基因进行京都基因与基因组百科全书（Kyoto encyclopedia of genes and genomes，KEGG）和基因本体论（gene ontology，GO）富集分析。 随后，基于 STRING 数据库（https：//string-db.org/）（highest confidence=0.900）构建 PPI 网络，并使用 Cytoscape 软件中的 CytoHubba 插件（Degree 算法）筛选前5个关键基因。最后对5个候选基因之间进行皮尔森相关性分析，通过R语言的lggcorrplot 包可视化。

DR中候选基因的验证和功能鉴定

为了进一步验证枢纽基因的诊断价值并探索其潜在功能，使用Sangerbox在线平台（版本3.0，http：//sangerbox.com/）进行了受试者工作特征（receiver operating characteristic，ROC）曲线分析，计算曲线下面积（area under the curve，AUC）和95%置信区间（confidence interval，CI），以评估GSE60436和GSE94019数据集中候选基因（COL1A1、COL1A2、FN1、TNF、TP53）的诊断性能。在GSE94019数据集中，按照候选基因中位表达水平将样本分成2组（高、低表达组）。基于MSigDB数据库（https：//www.gsea-msigdb.org/gsea/msigdb/）中的KEGG基因集进行基因集富集分析（gene set enrichment analysis，GSEA），以识别候选基因显著富集的信号通路。最后，比较2个数据集中DR组和对照组中候选基因的表达水平，并使用箱线图对结果进行可视化。

DR中候选基因的免疫浸润分析

将GSE94019数据集样本上传至CIBERSORTx平台（https：//cibersortx.stanford.edu/），以LM22基因集作为参考，分析每个样本中22种免疫细胞的浸润比例。随后评估5个血管新生相关候选基因［胶原蛋白Ⅰ型α1链（collagen type I alpha 1 chain， COL1A1）、COL1A2、FN1、肿瘤坏死因子（tumor necrosis factor， TNF）、肿瘤抑制蛋白53（tumor suppressor protein 53， TP53）］与免疫细胞浸润水平之间的Pearson相关性。

细胞处理和转染

在培养基中加入30 mmol/L葡萄糖处理HRMEC细胞24 h，以建立DR细胞模型。按照说明书说明，采用LipofectamineTM 3000转染试剂将siRNA对照（si-NC）或靶向FN1的siRNA（si-FN1，购自上海GenePharma公司）转染至HRMEC细胞，48 h后收集细胞用于后续实验。

qRT-PCR

采用TRIzol试剂提取细胞总RNA，利用PrimeScriptTM RT reagent Kit逆转录为cDNA。qRT-PCR在ABI 7500 Real-Time PCR System上进行，使用SYBR® Premix Ex TaqTM Ⅱ进行扩增。GAPDH作为内参基因，相对表达量采用2-ΔΔCT方法计算。引物序列见表1。

10.19405/j.cnki.issn1000–1492.2026.05.001.T001

引物序列

Primer sequences
GeneForward primer （5′-3′）Reverse primer （5′-3′）
FN1AAAGACCCCTTTCGTCACCCTCTTGTCCTACATTCGGCGG
GAPDHCTAGCTGGCCCGATTTCTCCATGGAATTTGCCATGGGTGG
Western blot实验

采用RIPA裂解液提取总蛋白，并使用BCA试剂盒测定蛋白浓度。取等量蛋白进行SDS-PAGE电泳分离后转移至PVDF膜。随后在室温下用5%脱脂奶粉封闭1 h，并加入FN1、P53和GAPDH一抗，于4 ℃条件下孵育过夜。次日加入HRP标记的山羊抗兔IgG二抗，使用ECL发光试剂进行显色，条带信号通过Bio-Rad凝胶成像系统检测。

细胞增殖实验

转染处理后的HRMEC细胞接种于96孔板（2×10³个/孔），在0、24、48、72、96 h检测细胞增殖情况。每个时间点加入10 μL CCK-8试剂，继续孵育2 h后，于450 nm波长下使用酶标仪测定吸光度（absorbance， A）值。

细胞迁移实验

将5×10⁴个细胞悬液（200 μL，无血清）加入Transwell上室，下室加入600 μL含10% FBS培养基作为趋化因子。37 ℃孵育24 h后，棉签擦去未迁移细胞，甲醇固定15 min并用0.1%结晶紫染色30 min。随机选取5个视野，显微镜拍照并计数迁移细胞数。

管腔形成实验

在96孔板中预铺Matrigel基质胶，于37 ℃孵育30 min凝胶化。将2×10⁴个处理后的HRMEC细胞接种于孔内，37 ℃孵育6 h后拍照。采用倒置显微镜观察，并用ImageJ软件分析管腔分支数与总长度。

统计学处理

所有实验均至少重复3次，数据以x¯±s表示。组间差异采用Student’s t检验（两组比较）或单因素方差分析（one-way ANOVA）（多组比较），并进行Tukey事后检验。数据分析在GraphPad Prism 9.0软件上完成。相关性分析采用 Pearson 相关系数。P<0.05为差异有统计学意义。

结果DR相关基因的富集分析

从GSE60436数据集中共筛选到2 550个上调基因和2 007个下调基因（图1A）；在GSE94019数据集中则筛选到4 029个上调基因和2 502个下调基因（图1B）。将这些DEGs 与血管新生相关基因取交集后，共得到237个重叠基因（其中227个上调，10个下调）（图2）。GO 富集分析结果显示，这些基因主要富集在生物过程（biological process，BP），分子功能（molecular function，MF），细胞组分（cellular component，CC）中的条目，包括血管生成（angiogenesis）、炎症反应（inflammatory response）以及蛋白结合（protein binding）等（图3A-3C）。在KEGG通路分析中，这些基因富集于PI3K-Akt、p53及JAK-STAT等经典信号通路（图3D）。

10.19405/j.cnki.issn1000–1492.2026.05.001.F001DR相关数据集中差异表达基因的鉴定Identification of differentially expressed genesin DR-related datasets

A： Volcano plot of DEGs in the GSE60436 dataset； B： Volcano plot of DEGs in the GSE94019 dataset.

10.19405/j.cnki.issn1000–1492.2026.05.001.F002GSE60436、GSE94019与血管新生相关基因的鉴定Identification of GSE60436， GSE94019， and angiogenesis-related genes

A： Venn diagram of upregulated genes； B： Venn diagram of downregulated genes.

10.19405/j.cnki.issn1000–1492.2026.05.001.F003糖尿病视网膜病变中血管新生相关基因的功能富集分析Functional enrichment analysis of angiogenesis-related genes in diabetic retinopathy

A： GO BP analysis indicated enrichment in angiogenesis， inflammatory response， cell adhesion， and hypoxia response； B： GO CC analysis revealed enrichment in extracellular matrix， basement membrane， and collagen-containing extracellular matrix； C： GO MF analysis showed enrichment in protein binding， collagen binding， and integrin binding； D： KEGG pathway enrichment indicated significant pathways including PI3K-Akt， p53， TNF， TGF-β， JAK-STAT， and AGE-RAGE signaling pathways.

DR中血管新生相关的关键基因的筛选和临床诊断价值分析

STRING数据库构建了DR中237个血管新生相关基因的PPI网络，结果显示162个基因之间存在紧密的相互作用关系（图4A）。进一步通过CytoHubba插件分析，筛选出前5个血管新生相关的关键基因（COL1A1、COL1A2、FN1、TNF、TP53）（图4B）。相关性分析显示，这5个关键基因之间呈显著的正相关关系（图4C）。为评估关键基因的诊断价值，利用GSE60436和GSE94019两个独立数据集绘制ROC 曲线。结果显示，在GSE60436数据集中，COL1A1、COL1A2、FN1、TNF和TP53的AUC值均为1.000（图5A-5E）。在 GSE94019 数据集中，COL1A1、COL1A2、FN1、TNF和TP53的AUC分别为1.000、0.889、0.917、0.806和0.889，进一步验证了这些基因在DR中具有较高的诊断价值（图5F-5J）。

10.19405/j.cnki.issn1000–1492.2026.05.001.F004PPI网络及候选基因相关性分析PPI network and correlation analysis of candidate genes

A： PPI network of angiogenesis-related genes in DR， consisting of 162 nodes and 298 edges； B： The Degree algorithm identified five candidate genes （COL1A1， COL1A2， FN1， TNF， and TP53）； C： Heat map showed the Pearson correlation between the five candidate genes.

10.19405/j.cnki.issn1000–1492.2026.05.001.F005GSE60436和GSE94019数据集中5个基因的ROC曲线分析ROC curve analysis of five genes in the GSE60436 and GSE94019 datasets

A-E： ROC curves for COL1A1， COL1A2， FN1， TNF， and TP53 in the GSE60436 dataset； F-J： ROC curves for COL1A1， COL1A2， FN1， TNF， and TP53 in the GSE94019 dataset.

GSEA结果揭示关键通路

基于GSE94019数据集，对5个候选基因的高低表达组进行了 GSEA分析。结果表明，COL1A1显著富集于Toll样受体（Toll like receptor，TLR）信号通路，细胞黏附分子（cell adhesion molecules，CAMs）等通路（图6A）；COL1A2与甘油脂代谢（glycerolipid metabolism，GLM）显著关联（图6B）；FN1富集于细胞周期、凋亡及ECM受体相互作用等通路（图6C）；TNF与JAK-STAT、NOD-like receptor 及TLR信号通路密切相关（图6D）；而TP53则富集于神经营养因子及P53信号通路等中（图 6E）。这些结果提示这5个基因可能通过多条信号通路参与DR的发生发展。

10.19405/j.cnki.issn1000–1492.2026.05.001.F006基于GSE94019数据集的候选基因的GSEAGSEA of candidate genes basedon the GSE94019 dataset

A-E： GSEA plots showed the enriched KEGG pathways associated with COL1A1， COL1A2， FN1， TNF， and TP53 expression in the GSE94019 dataset； The X-axis represents the rank of genes in the ordered dataset， while the Y-axis indicates the enrichment score （ES）； The normalized enrichment score （NES） and nominal P values are displayed in each panel. ES represents the enrichment score， and NP represents the nominal P value.

5个血管新生相关基因在DR中高表达

在独立数据集GSE60436和GSE94019中，对5个候选基因的表达水平进行验证。结果显示，这些基因在GSE60436数据集的DR组中较对照组显著上调（图 7A-7E，P<0.05）。在GSE94019数据集中，COL1A1、COL1A2、FN1和TP53也在DR样本中表达显著上调（图 7F-7H、7J，P<0.05），而TNF的差异无统计学意义（图7I）。

10.19405/j.cnki.issn1000–1492.2026.05.001.F007糖尿病视网膜病变中5个候选基因的表达水平验证Validation of the expression levels of five candidate genesin diabetic retinopathy

A-J： Boxplots showed the expression levels of COL1A1 （A， F）， COL1A2 （B， G）， FN1 （C， H）， TNF （D， I）， and TP53 （E， J） in DR samples and controls from the GSE60436 and GSE94019 datasets； *P<0.05， **P<0.01 vs Con group.

血管新生相关基因的免疫浸润分析

使用CIBERSORTx算法评估了GSE94019数据集样本中22种免疫细胞的浸润比例。结果显示，免疫细胞在不同样本中的分布差异明显（图8A）。随后对 5 个血管新生相关关键基因（COL1A1、COL1A2、FN1、TNF、TP53）与免疫细胞浸润水平进行Pearson相关性分析。结果表明，这些基因与多种免疫细胞显著相关，其中COL1A1、COL1A2与浆细胞，初始型CD4 T细胞，调节性T细胞（T cell regulatory，Tregs）等呈正相关；FN1和TP53与浆细胞，初始型CD4 T细胞呈显著正相关；而TNF与活化型自然杀伤（natural killer，NK）细胞、活化型肥大细胞呈正相关。结果提示这些基因可能通过免疫细胞浸润参与DR的发生发展（图8B）。

10.19405/j.cnki.issn1000–1492.2026.05.001.F008GSE94019 数据集中的免疫浸润分析Immune infiltration analysis in the GSE94019 dataset

A： Stacked barplot showed the fraction of 22 immune cell subsets in each sample， estimated by CIBERSORTx using the LM22 signature matrix； The X-axis lists the GSM sample IDs， and the Y-axis indicates the proportion of each immune cell type； B： Heatmap of Pearson correlations between the expression of five genes （COL1A1， COL1A2， FN1， TNF， TP53） and the infiltration levels of the 22 immune cell types.

<italic>FN1</italic>沉默抑制高糖（high glucose， HG）诱导的HRMEC的增殖

通过综合生物信息学分析筛选出5个潜在关键基因。综合差异表达分析、PPI网络结果及诊断效能曲线后，FN1表现出最高的连接度与诊断价值。qRT-PCR和Western blot结果表明，FN1在HG处理的HRMEC细胞中显著上调（图9A–9C）。进一步通过siRNA敲低FN1发现，其mRNA和蛋白水平均明显下降，验证了转染效率（图9D–9F）。此外，CCK-8检测显示，FN1沉默显著抑制了高糖处理后HRMEC的增殖能力（图9G）。

10.19405/j.cnki.issn1000–1492.2026.05.001.F00924 h高糖刺激下<italic>FN1</italic>敲低对HRMEC增殖的影响The effects of <italic>FN1</italic> knockdown on HRMEC proliferationunder high glucose stimulation

A： Quantitative real-time polymerase chain reaction （qRT-PCR） analysis of FN1 mRNA expression levels； B，C： Western blot analysis of FN1 protein expression levels； D： qRT-PCR analysis of FN1 knockdown efficiency in the DR cell model； E，F： Western blot analysis of FN1 knockdown efficiency in the DR cell model； G： After FN1 knockdown in the DR cell model， CCK-8 assay was used to detect the proliferation of HRMECs at 0， 24， 48， 72， and 96 h in different treatment groups； *P<0.05， **P<0.01， ***P<0.001 vs NG group； #P<0.05， ##P<0.01 vs HG+si-NC group.

<italic>FN1</italic>沉默抑制HG诱导的HRMEC的迁移和管形成并激活P53信号通路

Transwell实验结果显示，HG处理促进了HRMEC迁移，而FN1沉默显著减弱了这种效应（图10A、10B）。管腔形成实验进一步证实，HG诱导的血管生成能力显著增强，而FN1敲低显著抑制了管腔数量和分支长度（图10C–10E）。Western blot结果显示，HG处理下P53蛋白水平显著下调，而沉默FN1后P53表达部分恢复（图10F、10G）。提示FN1可能通过调控P53信号通路参与DR疾病发展。

10.19405/j.cnki.issn1000–1492.2026.05.001.F010<italic>FN1</italic>敲低对HRMEC迁移、血管生成和P53表达的影响The effects of <italic>FN1</italic> Knockdown on HRMEC migration， angiogenesis， and P53 expression

A，B： Transwell assays showed the migration of HRMECs under NG， HG， HG + si-NC， or HG + si-FN1 treatment； migrated cells were stained with crystal violet and quantified × 200； C-E： Tube formation assays on Matrigel， including representative images （C）， the number of tubes formed （D）， and quantification of branch length （E）× 200； F，G： Representative Western blot images of tumor suppressor P53 expression （F） and quantification normalized to GAPDH （G）； *P<0.05， **P<0.01， ***P<0.001 vs NG group； #P<0.05， ##P<0.01 vs HG + si-NC group.

讨论

随着DR病程进展，视网膜缺血区域会分泌大量促血管生成因子，从而诱导异常血管的生成［5］。血管新生程度与DR严重程度和视力丧失高度相关。本次研究得到与DR相关的237个血管新生相关基因。这些基因与血管生成、PI3K-Akt、P53、TNF、TGF-beta及JAK-STAT等通路相关。

本研究鉴定出5个关键基因（COL1A1、COL1A2、FN1、TNF和TP53）在DR样本中持续高表达，并具有较高的临床诊断价值。HG刺激通过激活TP53并促其核转位，加速铁死亡与内皮损伤，导致血视网膜屏障破坏［6］。这5个基因在多个经典通路中富集，包括JAK-STAT信号通路、Toll样受体信号通路和P53通路。在糖尿病小鼠模型中，JAK/STAT通路持续激活，导致异常血管生成并促进DR进展［7］。5'-N-乙基酰胺基腺苷通过抑制树突状细胞中的Toll样受体信号通路减轻DR的严重程度［8］。异槲皮素通过抑制P53信号通路减轻了糖尿病小鼠的视网膜损伤［9］。综上所述，5个关键基因可能通过多种信号通路驱动疾病进展。然而，未能进一步通过蛋白组数据进行交叉验证。

在DR进展过程中循环免疫细胞可浸润进入视网膜，其中先天免疫细胞的异常激活和浸润会引发血管及神经退化［10］。中性粒细胞、T细胞、B淋巴细胞以及单核/巨噬细胞等多类免疫细胞也在DR中表现出异常表达。此外，在DR样本中，B细胞、浆细胞、记忆性CD4 T细胞、Treg、M0巨噬细胞、M1巨噬细胞和中性粒细胞这7种类型的免疫细胞均显著过表达［11］。本次研究通过免疫浸润分析发现样本间免疫细胞分布存在显著差异。并且，血管新生相关关键基因与多种免疫细胞显著相关。这提示这些基因可能通过调控免疫细胞浸润参与 DR 的发生发展。

FN1已证实能够促进内皮细胞活性和促血管生成因子的表达［12］。本研究中，HG处理导致FN1表达上调。沉默FN1抑制HRMEC的增殖、迁移和管腔形成。研究［13］发现P53可通过抑制VEGF信号转导发挥血管生成负调控因子的作用。在机制层面，P53加剧了DR患者的视网膜内皮细胞凋亡［14］。高血糖条件下LINC00673通过负调控P53发挥视网膜色素上皮细胞保护作用［15］。在本次研究中，P53在HG诱导的HRMEC细胞中高表达，而沉默FN1可部分逆转这一结果，这表明FN1可能通过抑制P53信号传导来促进DR中的病理性血管生成。值得注意的是，磷酸化P53在细胞内更为稳定，并同样能够促进细胞衰老及抑制增殖。未来研究应进一步探究FN1对磷酸化P53的调控及其在糖尿病视网膜病变血管新生中的具体机制，以完善FN1-P53信号轴的作用模型。

综上所述，本研究鉴定了5个血管生成相关基因（COL1A1、COL1A2、FN1、TNF和TP53），作为DR的潜在关键调控因子，它们表现出较高的诊断价值，在多条经典信号通路中富集，并与免疫细胞浸润密切相关。功能实验进一步证实，FN1 在HG条件下促进HRMEC增殖、迁移和管腔形成，部分机制是通过调控 P53 信号通路实现的。总而言之，这些发现凸显了FN1作为DR病理性血管生成的关键介质的作用，并提示其可能成为疾病干预的有前景的生物标志物和治疗靶点。

参考文献AnsariP， TabasummaN， SnigdhaN N， et al. Diabetic retinopathy： an overview on mechanisms， pathophysiology and pharmacotherapy［J］. Diabetology， 2022， 3（1）： 159-75. doi：10.3390/diabetology3010011.LiuD， MengZ， JinC， et al. Fibronectin mediates endothelial-to-mesenchymal transition in retina angiogenesis［J］. Invest Ophthalmol Vis Sci， 2025， 66（3）： 10. doi：10.1167/iovs.66.3.10.WuJ， LinC， YangC， et al. Identification and validation of key biomarkers and potential therapeutic targets for primary open-angle glaucoma［J］. Sci China Life Sci， 2023， 66（12）： 2837-50. doi：10.1007/s11427-022-2344-5.YiC， LiuJ， DengW， et al. Old age promotes retinal fibrosis in choroidal neovascularization through circulating fibrocytes and profibrotic macrophages［J］. J Neuroinflammation， 2023， 20（1）： 45. doi：10.1186/s12974-023-02731-y.段　梅， 曹　凡， 桂衍超， 等. 塞来昔布下调JAML抑制糖尿病视网膜病变大鼠VEGF的表达及其机制［J］. 安徽医科大学学报， 2023， 58（8）： 1293-9. doi：10.19405/j.cnki.issn1000-1492.2023.08.008.DuanM， CaoF， GuiY C， et al. Mechanism of celecoxib inhibiting the expression of retinal VEGF in diabetic retinopathy rats via JAML［J］. Acta Univ Med Anhui， 2023， 58（8）： 1293-9. doi：10.19405/j.cnki.issn1000-1492.2023.08.008.LuC， LanQ， SongQ， et al. Identification and validation of ferroptosis-related genes for diabetic retinopathy［J］. Cell Signal， 2024， 113： 110955. doi：10.1016/j.cellsig.2023.110955.ChoC H， RohK H， LimN Y， et al. Role of the JAK/STAT pathway in a streptozotocin-induced diabetic retinopathy mouse model［J］. Graefes Arch Clin Exp Ophthalmol， 2022， 260（11）： 3553-63. doi：10.1007/s00417-022-05694-7.LiL， ChenJ， WangZ， et al. NECA alleviates inflammatory responses in diabetic retinopathy through dendritic cell toll-like receptor signaling pathway［J］. Front Immunol， 2024， 15： 1415004. doi：10.3389/fimmu.2024.1415004.CaiY， PengS， DuanB， et al. Isoquercetin alleviates diabetic retinopathy via inhibiting p53-mediated ferroptosis［J］. Cell Biol Int， 2025， 49（7）： 852-64. doi：10.1002/cbin.70027.ZhaoB， ZhaoY， SunX. Mechanism and therapeutic targets of circulating immune cells in diabetic retinopathy［J］. Pharmacol Res， 2024， 210： 107505. doi：10.1016/j.phrs.2024.107505.袁琳慧，张立军，刘 新，等.基于加权基因共表达网络识别糖尿病视网膜病变与免疫相关的关键基因［J］.国际眼科杂志， 2023， 23（8）： 1343-1351. doi：10.3980/j.issn.1672-5123. 2023.8.20.YuanL H， ZhangL J， LiuX， et al. Identification of key immune related genes in diabetes retinopathy based on weighted gene co-expression network ［J］. International Eye Science， 2023： 1343-51. doi：10.3980/j.issn.1672-5123.2023.8.20.ChenC， WangF， ChengC， et al. Cancer-associated fibroblasts-derived exosomes with HOXD11 overexpression promote ovarian cancer cell angiogenesis via FN1［J］. Reprod Sci， 2025， 32（5）： 1530-44. doi：10.1007/s43032-024-01716-3.ZhangC， LiuJ， WangJ， et al. The interplay between tumor suppressor p53 and hypoxia signaling pathways in cancer［J］. Front Cell Dev Biol， 2021， 9： 648808. doi：10.3389/fcell. 2021.648808.ChengY， ZhangM， XuR， et al. p53 accelerates endothelial cell senescence in diabetic retinopathy by enhancing FoxO3a ubiquitylation and degradation via UBE2L6［J］. Exp Gerontol， 2024， 188： 112391. doi：10.1016/j.exger.2024.112391.ChengY， ZhuY， MaL. LncRNA LINC00673 is downregulated in diabetic retinopathy and regulates the apoptosis of retinal pigment epithelial cells via negatively regulating p53［J］. Diabetes， Metabolic Syndrome and Obesity， 2021， 14： 4233-40. doi：10.2147/DMSO.S298185.

李鹏, 梁坤, 吴峰, 等. 综合生物信息学分析及实验验证糖尿病视网膜病变中血管新生相关基因的作用[J]. 安徽医科大学学报, 2026, 61(05): 861-871.

Li Peng, Liang Kun, Wu Feng, et al. Integrated bioinformatics analysis and experimental validation of angiogenesis-related genes in diabetic retinopathy[J]. Acta Universitatis Medicinalis Anhui, 2026, 61(05): 861-871.

1， 梁坤2，吴峰1，李佳3，刘伦4，陶玉林2"?>1眼科、3神经外科，阜阳　236000；2安徽医科大学第二附属医院眼科，合肥　230601；4安徽医科大学第一附属医院眼科，合肥　230022）"?>1， Liang Kun2， Wu Feng1， Li Jia3， Liu Lun4， Tao Yulin 2"?>1Department of Ophthalmology， 3Department of Neurosurgery， Fuyang Hospital Affiliated toAnhui Medical University， Fuyang　236000；2Department of Ophthalmology， The Second AffiliatedHospital of Anhui Medical University， Hefei　230601； 4Department of Ophthalmology，The First Affiliated Hospital of Anhui Medical University， Hefei　230022）"?>