Single-cell RNA-seq and bulk-seq identify RAB17 as a potential regulator of angiogenesis by human dermal microvascular endothelial cells in diabetic foot ulcers

Abstract

Background: Angiogenesis is crucial in diabetic wound healing and is often impaired in diabetic foot ulcers (DFUs). Human dermal microvascular endothelial cells (HDMECs) are vital components in dermal angiogenesis; however, their functional and transcriptomic characteristics in DFU patients are not well understood. This study aimed to comprehensively analyse HDMECs from DFU patients and healthy controls and find the potential regulator of angiogenesis in DFUs. Methods: HDMECs were isolated from skin specimens of DFU patients and healthy controls via magnetic-activated cell sorting. The proliferation, migration and tube-formation abilities of the cells were then compared between the experimental groups. Both bulk RNA sequencing (bulk-seq) and single-cell RNA-seq (scRNA-seq) were used to identify RAB17 as a potential marker of angiogenesis, which was further confirmed via weighted gene co-expression network analysis (WGCNA) and least absolute shrink and selection operator (LASSO) regression. The role of RAB17 in angiogenesis was examined through in vitro and in vivo experiments. Results: The isolated HDMECs displayed typical markers of endothelial cells. HDMECs isolated from DFU patients showed considerably impaired tube formation, rather than proliferation or migration, compared to those from healthy controls. Gene set enrichment analysis (GSEA), fGSEA, and gene set variation analysis (GSVA) of bulk-seq and scRNA-seq indicated that angiogenesis was downregulated in DFU-HDMECs. LASSO regression identified two genes, RAB17 and CD200, as characteristic of DFU-HDMECs; additionally, the expression of RAB17 was found to be significantly reduced in DFU-HDMECs compared to that in the HDMECs of healthy controls. Overexpression of RAB17 was found to enhance angiogenesis, the expression of hypoxia inducible factor-1α and vascular endothelial growth factor A, and diabetic wound healing, partially through the mitogen-activated protein kinase/extracellular signal-regulated kinase signalling pathway. Conclusions: Our findings suggest that the impaired angiogenic capacity in DFUs may be related to the dysregulated expression of RAB17 in HDMECs. The identification of RAB17 as a potential molecular target provides a potential avenue for the treatment of impaired angiogenesis in DFUs.

Keywords

Human dermal microvascular endothelial cells / Angiogenesis / Diabetic foot ulcers / RAB17 / Single-cell RNA-seq / Diabetic wound healing

References

[1] Chen H, Cheng Y, Tian J, Yang P, Zhang X, Chen Y, et al. Dissolved oxygen from microalgae-gel patch promotes chronic wound healing in diabetes. Sci Adv. 2020; 6: eaba4311. https://doi.org/10.1126/sciadv.aba4311.

[2] Fesseha BK, Abularrage CJ, Hines KF, Sherman R, Frost P, Langan S, et al. Association of Hemoglobin a(1c) and wound healing in diabetic foot ulcers. Diabetes Care. 2018; 41: 1478-85. https://doi.org/10.2337/dc17-1683.

[3] Okonkwo UA, DiPietro LA. Diabetes and wound angiogenesis. Int J Mol Sci. 2017; 18. https://doi.org/10.3390/ijms18071419.

[4] Wu X, He W, Mu X, Liu Y, Deng J, Liu Y, et al. Macrophage polarization in diabetic wound healing. Burns Dent Traumatol. 2022; 10: tkac051. https://doi.org/10.1093/burnst/tkac051.

[5] Ishihara J, Ishihara A, Starke RD, Peghaire CR, Smith KE, McKinnon TAJ, et al. The heparin binding domain of von Willebrand factor binds to growth factors and promotes angiogenesis in wound healing. Blood. 2019; 133: 2559-69. https://doi.org/10.1182/blood.2019000510.

[6] Cao XZ, Xie T, Lu SL. Receptor pathways of glycated basic fibroblast growth factor affecting the proliferation and vascularization of human dermal microvascular endothelial cells. Chin j Burns. 2021; 37: 17-24. https://doi.org/10.3760/cma.j.cn501120-20200916-00412.

[7] Bouïs D, Kusumanto Y, Meijer C, Mulder NH, Hospers GA. A review on pro-and anti-angiogenic factors as targets of clinical intervention. Pharmacol Res. 2006; 53: 89-103. https://doi.org/10.1016/j.phrs.2005.10.006.

[8] Dixon D, Edmonds M.Managing diabetic foot ulcers: pharmacotherapy for wound healing. Drugs. 2021; 81: 29-56. https://doi.org/10.1007/s40265-020-01415-8.

[9] Tsou PS, Campbell P, Amin MA, Coit P, Miller S, Fox DA, et al. Inhibition of EZH2 prevents fibrosis and restores normal angiogenesis in scleroderma. Proc Natl Acad Sci U S A. 2019; 116: 3695-702. https://doi.org/10.1073/pnas.1813006116.

[10] Fandaros M, Joseph K, Kaplan AP, Rubenstein DA, Ghebrehiwet B, Yin W. gC1qR antibody can modulate endothelial cell permeability in angioedema. Inflammation. 2022; 45: 116-28. https://doi.org/10.1007/s10753-021-01532-w.

[11] Peng L, Li Q, Wang H, Wu J, Li C, Liu Y, et al. Fn14 deficiency ameliorates psoriasis-like skin disease in a murine model. Cell Death Dis. 2018; 9: 801. https://doi.org/10.1038/s41419-018-0820-6.

[12] Lawley TJ, Kubota Y. Induction of morphologic differentiation of endothelial cells in culture. J Invest Dermatol. 1989; 93: 59s-61. https://doi.org/10.1111/1523-1747.ep12581070.

[13] Clark PR, Manes TD, Pober JS, Kluger MS. Increased ICAM-1 expression causes endothelial cell leakiness, cytoskeletal reorganization and junctional alterations. J Invest Dermatol. 2007; 127: 762-74. https://doi.org/10.1038/sj.jid.5700670.

[14] Beckman JA, Doherty SP, Feldman ZB, Banks ES, Moslehi J, Jaffe IZ, et al. Comparative Transcriptomics of ex vivo, patient-derived endothelial cells reveals novel pathways associated with type 2 diabetes mellitus. JACC Basic Transl Sci. 2019; 4: 567-74. https://doi.org/10.1016/j.jacbts.2019.05.012.

[15] Singh K, Sinha M, Pal D, Tabasum S, Gnyawali SC, Khona D, et al. Cutaneous epithelial to mesenchymal transition activator ZEB1 regulates wound angiogenesis and closure in a Glycemic status-dependent manner. Diabetes. 2019; 68: 2175-90. https://doi.org/10.2337/db19-0202.

[16] Botusan IR, Sunkari VG, Savu O, Catrina AI, Grünler J, Lindberg S, et al. Stabilization of HIF-1alpha is critical to improve wound healing in diabetic mice. Proc Natl Acad Sci U S A. 2008; 105: 19426-31. https://doi.org/10.1073/pnas.0805230105.

[17] Rivard A, Silver M, Chen D, Kearney M, Magner M, Annex B, et al. Rescue of diabetes-related impairment of angiogenesis by intramuscular gene therapy with adeno-VEGF. Am J Pathol. 1999; 154: 355-63. https://doi.org/10.1016/s0002-9440(10)65282-0.

[18] Qi J, Zhao P, Li F, Guo Y, Cui H, Liu A, et al. Down-regulation of Rab17 promotes tumourigenic properties of hepatocellular carcinoma cells via Erk pathway. Int J Clin Exp Pathol. 2015; 8: 4963-71.

[19] Wang M, Wang W, Ding J, Wang J, Zhang J. Down-regulation of Rab17 promotes cell proliferation and invasion in non-small cell lung cancer through STAT3/HIF-1α/VEGF signaling. Thorac Cancer. 2020; 11: 379-88. https://doi.org/10.1111/1759-7714.13278.

[20] Striz AC, Tuma PL. The GTP-bound and Sumoylated form of the rab17 small molecular weight GTPase selectively binds Syntaxin 2 in polarized hepatic WIF-B cells. J Biol Chem. 2016; 291: 9721-32. https://doi.org/10.1074/jbc.M116.723353.

[21] El-Remessy AB, Abou-Mohamed G, Caldwell RW, Caldwell RB. High glucose-induced tyrosine nitration in endothelial cells: role of eNOS uncoupling and aldose reductase activation. Invest Ophthalmol Vis Sci. 2003; 44: 3135-43. https://doi.org/10.1167/iovs.02-1022.

[22] Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15: 550. https://doi.org/10.1186/s13059-014-0550-8.

[23] Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012; 16: 284-7. https://doi.org/10.1089/omi.2011.0118.

[24] Korotkevich G, Sukhov V, Budin N, Shpak B, Artyomov MN, Sergushichev A. Fast gene set enrichment analysis. bioRxiv. 2021; 060012. https://doi.org/10.1101/060012.

[25] Hänzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics. 2013; 14: 7. https://doi.org/10.1186/1471-2105-14-7.

[26] Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015; 43: e47. https://doi.org/10.1093/nar/gkv007.

[27] Hao Y, Hao S, Andersen-Nissen E, Mauck WM, 3rd, Zheng S, Butler A, et al. Integrated analysis of multimodal single-cell data. Cell. 2021; 184: 3573-87.e29. https://doi.org/10.1016/j.cell.2021.04.048.

[28] Squair JW, Gautier M, Kathe C, Anderson MA, James ND, Hutson TH, et al. Confronting false discoveries in single-cell differential expression. Nat Commun. 2021; 12: 5692. https://doi.org/10.1038/s41467-021-25960-2.

[29] Langfelder P, Horvath S. Fast R functions for robust correlations and hierarchical clustering. J Stat Softw. 2012; 46: i11.

[30] Simon N, Friedman J, Hastie T, Tibshirani R. Regularization paths for Cox’s proportional hazards model via coordinate descent. J Stat Softw. 2011; 39: 1-13. https://doi.org/10.18637/jss.v039.i05.

[31] Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez JC, et al. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics. 2011; 12: 77. https://doi.org/10.1186/1471-2105-12-77.

[32] Zhao J, Yang S, Shu B, Chen L, Yang R, Xu Y, et al. Transient high glucose causes persistent vascular dysfunction and delayed wound healing by the DNMT1-mediated Ang-1/NF-κB pathway. J Invest Dermatol. 2021; 141: 1573-84. https://doi.org/10.1016/j.jid.2020.10.023.

[33] Theocharidis G, Thomas BE, Sarkar D, Mumme HL, Pilcher WJR, Dwivedi B, et al. Single cell transcriptomic landscape of diabetic foot ulcers. Nat Commun. 2022; 13: 181. https://doi.org/10.1038/s41467-021-27801-8.

[34] Shim J, Oh SJ, Yeo E, Park JH, Bae JH, Kim SH, et al. Integrated analysis of single-cell and spatial Transcriptomics in keloids: highlights on Fibrovascular interactions in keloid pathogenesis. J Invest Dermatol. 2022. https://doi.org/10.1016/j.jid.2022.01.017.

[35] LoGerfo FW, Coffman JD. Current concepts. Vascular and microvascular disease of the foot in diabetes. Implications for foot care. N Engl J Med. 1984; 311: 1615-9. https://doi.org/10.1056/nejm198412203112506.

[36] Falanga V. Wound healing and its impairment in the diabetic foot. Lancet. 2005; 366: 1736-43. https://doi.org/10.1016/s0140-6736(05)67700-8.

[37] Vorwald CE, Murphy KC, Leach JK. Restoring vasculogenic potential of endothelial cells from diabetic patients through spheroid formation. Cell Mol Bioeng. 2018; 11: 267-78. https://doi.org/10.1007/s12195-018-0531-1.

[38] Huang X, Sun J, Chen G, Niu C, Wang Y, Zhao C, et al. Resveratrol promotes diabetic wound healing via SIRT1-FOXO1-c-Myc Signaling pathway-mediated angiogenesis. Front Pharmacol. 2019; 10: 421. https://doi.org/10.3389/fphar.2019.00421.

[39] Chen C, Zhao C, Gu C, Cui X, Wu J. MiRNA-144-3p inhibits high glucose induced cell proliferation through suppressing FGF16. Biosci Rep. 2019; 39. https://doi.org/10.1042/bsr20181788.

[40] Henrot P, Laurent P, Levionnois E, Leleu D, Pain C, Truchetet ME, et al. A method for isolating and culturing skin cells: application to endothelial cells, fibroblasts, keratinocytes, and melanocytes from punch biopsies in systemic sclerosis skin. Front Immunol. 2020; 11: 566607. https://doi.org/10.3389/fimmu.2020.566607.

[41] Tsou PS, Amin MA, Campbell PL, Zakhem G, Balogh B, Edhayan G, et al. Activation of the thromboxane A2 receptor by 8-Isoprostane inhibits the pro-Angiogenic effect of vascular endothelial growth factor in scleroderma. J Invest Dermatol. 2015; 135: 3153-62. https://doi.org/10.1038/jid.2015.323.

[42] Hu SC, Lan CE. High-glucose environment disturbs the physiologic functions of keratinocytes: focusing on diabetic wound healing. J Dermatol Sci. 2016; 84: 121-7. https://doi.org/10.1016/j.jdermsci.2016.07.008.

[43] Zhang S, Ke Z, Yang C, Zhou P, Jiang H, Chen L, et al. High glucose causes distinct expression patterns of primary human skin cells by RNA sequencing. Front Endocrinol (Lausanne). 2021; 12: 603645. https://doi.org/10.3389/fendo.2021.603645.

[44] Wimmer RA, Leopoldi A, Aichinger M, Wick N, Hantusch B, Novatchkova M, et al. Human blood vessel organoids as a model of diabetic vasculopathy. Nature. 2019; 565: 505-10. https://doi.org/10.1038/s41586-018-0858-8.

[45] Dong L, Zhang Z, Liu X, Wang Q, Hong Y, Li X, et al. RNA sequencing reveals BMP4 as a basis for the dual-target treatment of diabetic retinopathy. J Mol Med (Berl). 2021; 99: 225-40. https://doi.org/10.1007/s00109-020-01995-8.

[46] Li Q, Zhu Z, Wang L, Lin Y, Fang H, Lei J, et al. Single-cell transcriptome profiling reveals vascular endothelial cell heterogeneity in human skin. Theranostics. 2021; 11: 6461-76. https://doi.org/10.7150/thno.54917.

[47] Paik DT, Tian L, Williams IM, Rhee S, Zhang H, Liu C, et al. Single-cell RNA sequencing unveils unique transcriptomic signatures of organ-specific endothelial cells. Circulation. 2020; 142: 1848-62. https://doi.org/10.1161/circulationaha.119.041433.

[48] Rustagi Y, Abouhashem AS, Verma P, Verma SS, Hernandez E, Liu S, et al. Endothelial phospholipase Cγ2 improves outcomes of diabetic ischemic limb rescue following VEGF therapy. Diabetes. 2022; 71: 1149-65. https://doi.org/10.2337/db21-0830.

[49] Ding J, Adiconis X, Simmons SK, Kowalczyk MS, Hession CC, Marjanovic ND, et al. Systematic comparison of single-cell and single-nucleus RNA-sequencing methods. Nat Biotechnol. 2020; 38: 737-46. https://doi.org/10.1038/s41587-020-0465-8.

[50] Chung KB, Oh J, Roh WS, Kim TG, Kim DY. Core gene signatures of atopic dermatitis using public RNA-sequencing resources: comparison of bulk approach with single-cell approach. J Invest Dermatol. 2022; 142: 717-21.e5. https://doi.org/10.1016/j.jid.2021.07.169.

[51] Luo Q, Chen Q, Wang W, Desrivières S, Quinlan EB, Jia T, et al. Association of a Schizophrenia-Risk Nonsynonymous Variant with putamen volume in adolescents: a Voxelwise and genome-wide association study. JAMA Psychiatry. 2019; 76: 435-45. https://doi.org/10.1001/jamapsychiatry.2018.4126.

[52] Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008; 9: 559. https://doi.org/10.1186/1471-2105-9-559.

[53] Liu W, Huang X, Liang X, Zhou Y, Li H, Yu Q, et al. Identification of key modules and hub genes of keloids with weighted gene Coexpression network analysis. Plast Reconstr Surg. 2017; 139: 376-90. https://doi.org/10.1097/prs.0000000000003014.

[54] Wang H, Zhou Z, Liu Y, Wang P, Chen L, Qi S, et al. Identification and validation of HOXD3 and UNC5C as molecular signatures in keloid based on weighted gene co-expression network analysis. Genomics. 2022; 114: 110403. https://doi.org/10.1016/j.ygeno.2022.110403.

[55] Li Z, Huang B, Yi W, Wang F, Wei S, Yan H, et al. Identification of potential early diagnostic biomarkers of sepsis. J Inflamm Res. 2021; 14: 621-31. https://doi.org/10.2147/jir.S298604.

[56] Rütsche D, Michalak-Micka K, Zielinska D, Moll H, Moehrlen U, Biedermann T, et al. The role of CD200-CD200 receptor in human blood and lymphatic endothelial cells in the regulation of skin tissue inflammation. Cells. 2022; 11. https://doi.org/10.3390/cells11061055.

[57] Lütcke A, Jansson S, Parton RG, Chavrier P, Valencia A, Huber LA, et al. Rab17, a novel small GTPase, is specific for epithelial cells and is induced during cell polarization. J Cell Biol. 1993; 121: 553-64. https://doi.org/10.1083/jcb.121.3.553.

[58] Beaumont KA, Hamilton NA, Moores MT, Brown DL, Ohbayashi N, Cairncross O, et al. The recycling endosome protein Rab17 regulates melanocytic filopodia formation and melanosome trafficking. Traffic. 2011; 12: 627-43. https://doi.org/10.1111/j.1600-0854.2011.01172.x.

[59] Mori Y, Matsui T, Furutani Y, Yoshihara Y, Fukuda M. Small GTPase Rab17 regulates dendritic morphogenesis and postsynaptic development of hippocampal neurons. J Biol Chem. 2012; 287: 8963-73. https://doi.org/10.1074/jbc.M111.314385.

[60] Yin C, Kim Y, Argintaru D, Heit B. Rab17 mediates differential antigen sorting following efferocytosis and phagocytosis. Cell Death Dis. 2016; 7: e2529. https://doi.org/10.1038/cddis.2016.431.

[61] Haobam B, Nozawa T, Minowa-Nozawa A, Tanaka M, Oda S, Watanabe T, et al. Rab17-mediated recycling endosomes contribute to autophagosome formation in response to group a streptococcus invasion. Cell Microbiol. 2014; 16: 1806-21. https://doi.org/10.1111/cmi.12329.

[62] Diaz-Vera J, Palmer S, Hernandez-Fernaud JR, Dornier E, Mitchell LE, Macpherson I, et al. A proteomic approach to identify endosomal cargoes controlling cancer invasiveness. J Cell Sci. 2017; 130: 697-711. https://doi.org/10.1242/jcs.190835.

[63] Zhou X, Xia G, Liu Y, Xin X, Shi R, An L, et al. Analysis of carcinogenic signaling networks in endometrial cancer identifies RAB17 as a potential target. J Cell Physiol. 2021; 236: 328-39. https://doi.org/10.1002/jcp.29845.

[64] Arnaoutova I, George J, Kleinman HK, Benton G. The endothelial cell tube formation assay on basement membrane turns 20: state of the science and the art. Angiogenesis. 2009; 12: 267-74. https://doi.org/10.1007/s10456-009-9146-4.

[65] Mace KA, Yu DH, Paydar KZ, Boudreau N, Young DM. Sustained expression of Hif-1alpha in the diabetic environment promotes angiogenesis and cutaneous wound repair. Wound Repair Regen. 2007; 15: 636-45. https://doi.org/10.1111/j.1524-475X.2007.00278.x.

[66] Gao Y, Xie Z, Ho C, Wang J, Li Q, Zhang Y, et al. LRG1 promotes keratinocyte migration and wound repair through regulation of HIF-1α stability. J Invest Dermatol. 2020; 140: 455-64.e8. https://doi.org/10.1016/j.jid.2019.06.143.

[67] Carmeliet P, Jain RK. Molecular mechanisms and clinical applications of angiogenesis. Nature. 2011; 473: 298-307. https://doi.org/10.1038/nature10144.

[68] Wu S, Luo M, To KKW, Zhang J, Su C, Zhang H, et al. Intercellular transfer of exosomal wild type EGFR triggers osimertinib resistance in non-small cell lung cancer. Mol Cancer. 2020; 19: 145. https://doi.org/10.1186/s12943-020-01264-2.

[69] Liu H, Tang L, Jiang X, He X, Liang J, Liao X, et al. Autophagy deficiency in endothelial cells accelerates wound healing in mice possibly in a paracrine-dependent fashion. J Invest Dermatol. 2021. https://doi.org/10.1016/j.jid.2021.05.012.