Diabetic kidney disease induces transcriptome alterations associated with angiogenesis activity in human mesenchymal stromal cells
Stem Cell Research & Therapy volume 14, Article number: 49 (2023)
Therapeutic interventions that optimize angiogenic activities may reduce rates of end-stage kidney disease, critical limb ischemia, and lower extremity amputations in individuals with diabetic kidney disease (DKD). Infusion of autologous mesenchymal stromal cells (MSC) is a promising novel therapy to rejuvenate vascular integrity. However, DKD-related factors, including hyperglycemia and uremia, might alter MSC angiogenic repair capacity in an autologous treatment approach.
To explore the angiogenic activity of MSC in DKD, the transcriptome of adipose tissue-derived MSC obtained from DKD subjects was compared to age-matched controls without diabetes or kidney impairment. Next-generation RNA sequencing (RNA-seq) was performed on MSC (DKD n = 29; Controls n = 9) to identify differentially expressed (DE; adjusted p < 0.05, |log2fold change|> 1) messenger RNA (mRNA) and microRNA (miRNA) involved in angiogenesis (GeneCards). Paracrine-mediated angiogenic repair capacity of MSC conditioned medium (MSCcm) was assessed in vitro using human umbilical vein endothelial cells incubated in high glucose and indoxyl sulfate for a hyperglycemic, uremic state.
RNA-seq analyses revealed 133 DE mRNAs (77 upregulated and 56 down-regulated) and 208 DE miRNAs (119 up- and 89 down-regulated) in DKD-MSC versus Control-MSC. Interestingly, miRNA let-7a-5p, which regulates angiogenesis and participates in DKD pathogenesis, interacted with 5 angiogenesis-associated mRNAs (transgelin/TAGLN, thrombospondin 1/THBS1, lysyl oxidase-like 4/LOXL4, collagen 4A1/COL4A1 and collagen 8A1/COL8A1). DKD-MSCcm incubation with injured endothelial cells improved tube formation capacity, enhanced migration, reduced adhesion molecules E-selectin, vascular cell adhesion molecule 1 and intercellular adhesion molecule 1 mRNA expression in endothelial cells. Moreover, angiogenic repair effects did not differ between treatment groups (DKD-MSCcm vs. Control-MSCcm).
MSC from individuals with DKD show angiogenic transcriptome alterations compared to age-matched controls. However, angiogenic repair potential may be preserved, supporting autologous MSC interventions to treat conditions requiring enhanced angiogenic activities such as DKD, diabetic foot ulcers, and critical limb ischemia.
Tissue hypoxemia from loss of microvascular capillary density is a major contributor to morbidity and mortality in individuals with diabetic kidney disease (DKD), one of the most common causes of kidney failure worldwide . Sustained hyperglycemia and uremia alter vascular function and integrity in the kidney and limbs wherein disease pathogenesis is further fueled by microvascular rarefaction, tissue hypoxemia, sterile inflammation, cellular senescence, endothelial dysfunction, and extracellular matrix deposition [2,3,4,5,6]. These effects culminate in complications common in DKD including end-stage kidney disease (ESKD), cardiovascular disease, critical limb ischemia, and amputations of lower extremities . Despite a 2-decade decline in lower extremity amputations between 2009 and 2015, rates rebounded by 50% in adults with diabetes . Furthermore, the odds of major amputations are nearly twofold higher in individuals with chronic kidney disease (CKD) or ESKD compared to persons without CKD/ESRD . Therefore, therapeutics aimed at optimizing vascular repair and relieving endothelial dysfunction in the DKD patient population are urgently needed.
Regenerative cell-based therapies, such as mesenchymal stem/stromal cells (MSC), have considerable potential to boost endogenous repair of injured tissues. MSC release soluble factors, such as extracellular vesicles (EVs), growth factors, and cytokines, which enhance tissue regeneration and stimulate pro-angiogenic signaling [10,11,12]. In preclinical models of DKD, renovascular disease, and ischemic limb disease, MSC-based approaches improve kidney function, reduce hypoxemia, and restore microvascular density [2, 11, 13,14,15]. Collectively, these encouraging studies provided the basis for clinical application of MSC in phase I and II studies in DKD (NCT05362786, NCT02585622, NCT03840343, NCT04869761, NCT04125329, NCT04216849, NCT02008851, NCT03270956, NCT02836574)  and diabetic foot ulcers (NCT04464213, NCT04104451, NCT00955669, NCT04466007, NCT02375802, NCT01686139) .
Despite the promise of MSC in DKD, numerous challenges remain. Chronic hyperglycemia and uremia may damage endogenous repair systems like MSC, disturbing the coordinated network of pro- and anti-angiogenic factors activated by these cells. Autologous (versus allogeneic) MSC source may be preferable, given reduced risk of allosensitization, particularly for patients that might ultimately require kidney transplantation . However, multiple DKD-related factors including age, hyperglycemia, uremia, cellular senescence abundance, or obesity may alter MSC functionality and number [13, 19,20,21]. We and others demonstrated that immunomodulatory and paracrine activities of MSC are maintained despite kidney impairment with and without diabetes [21,22,23]. However, it is unclear whether the microenvironment of DKD negatively influences the angiogenic potential of MSC. Elucidation of pro-angiogenic and tissue-repairing potential of DKD-MSC in the development of MSC-based regenerative medicine strategies for this population is crucial.
Our previous phase 1 clinical trial showed that intra-arterial infusion of autologous adipose-derived MSC (AD-MSC) improved kidney blood flow and reduced kidney tissue hypoxia as assessed through blood oxygen level dependent magnetic resonance imaging (BOLD-MRI) and improved glomerular filtration rate in atherosclerotic renovascular disease subjects . These findings reflected AD-MSC angiogenic capabilities to form new vessels and restore the microcirculation, consistent with our experimental observations in swine atherosclerotic renal artery stenosis [2, 10, 25]. AD-MSC remain a viable cell option for clinical applications given their high abundance in adipose tissue and ease of harvest when compared to bone marrow-derived MSC (BM-MSC). Interestingly, AD-MSC proved more effective than BM-MSC in promoting neovascularization in animal models of ischemic diseases [26, 27]. Given the need for regenerative investigations in CKD, we are currently conducting two phase 1 clinical trials testing AD-MSC (NCT04869761) and BM-MSC (NCT05362786) in the treatment of DKD/CKD. Therefore, based on our prior experiences and the literature, we compared the transcriptome of AD-MSC obtained from DKD subjects and age-matched controls focusing on angiogenesis-related mRNA and miRNA targeting these genes. We further compared the angiogenic repair capacity in vitro of AD-MSC from DKD subjects and controls.
The expression of angiogenesis-related mRNAs and miRNAs and angiogenic potential was characterized in AD-MSC isolated from diabetic kidney disease (DKD) subjects. Eligible individuals age ≥ 18 years were recruited from Mayo Clinic in Rochester, Minnesota, in the Nephrology clinic between November 2015–2019 (IRB: 15-000933), as previously reported . Informed consent was obtained by trained clinical research coordinators. DKD was defined as estimated glomerular filtration rate (eGFR) < 60 mL/min/1.73 m2 and/or abnormal albuminuria with preserved kidney function (eGFR ≥ 60) in the setting of pharmacologically treated Type 1 or Type 2 diabetes mellitus. Controls were age-matched individuals, without diabetes or CKD, undergoing laparoscopic nephrectomy for kidney donation. Adipose tissue was collected from DKD (n = 29) participants in an outpatient surgical suite, while adipose tissue from controls (n = 9) was sampled during their surgical procedure. To minimize bias, those with dialysis dependency, kidney transplant, immunosuppressive therapy, hemoglobin A1c (HbA1c) > 11% (97 mmol/mol), or malignancy were excluded. The Mayo Clinic Institutional Review Board approved all experimental study procedures and all participants provided written informed consent prior to participating.
MSC harvest and characterization
Abdominal subcutaneous adipose tissue (0.5–2.0 g) was collected from control and DKD subjects to isolate MSC as previously described [20, 21, 25]. MSC were cultured in standard conditions (37 °C with 5% CO2) in advanced minimum essential medium supplemented with 5% platelet lysate  (PLTmax, Mill Creek Life Sciences, Rochester, MN). Primary MSC at passages 3–5 were phenotyped in both groups to confirm cell surface marker expression by imaging flow cytometry (Amnis® FlowSight®, Austin, TX) [19,20,21]. The following fluorescently labeled antibodies were used: cluster of differentiation (CD)73 (ab106677; Abcam, Waltham, MA), CD90 (ab124527; Abcam), and CD105 (ab53321; Abcam). Conversely, MSC lack expression of CD34 (340441; BD Biosciences, San Jose, CA) or CD14 (ab82012; Abcam). For surface marker expression, 0.25–1.0 × 106 cells were incubated with the antibodies for 30 min at 4 °C according to the manufacturer’s instructions. Following incubation, cell suspensions were washed, pelleted, and resuspended in staining buffer (00-4222-57; Thermo Fisher Scientific, Waltham, MA). Thereafter, samples were run on the FlowSight® and analyzed using IDEAS software version 6.2. Initial gates were selected for single cells (aspect ratio intensity versus brightfield), followed by gates of positively stained cells. Thresholds for positivity were established from scatterplots generated by single-stained compensation beads (A10497, Molecular Probes, Eugene, OR).
Collection of MSC conditioned medium (MSCcm)
MSC (passages 3–5) were cultured to 70% confluence in T-75 flasks, then washed extensively with phosphate-buffered saline, and replenished with serum-free advanced minimum essential medium for 36 h. After serum starvation, MSCcm was centrifuged at 2000×g for 20 min at 4 °C to remove cellular debris, aliquoted and stored at − 80 °C for in vitro studies. MSCcm was normalized by MSC density per sample (1 mL conditioned medium for every 0.13 × 106 MSC) . Thrombospondin 1 (TSP1) protein was measured in MSCcm using a commercially available kit (ELISA; R&D DTSP10) following manufacturer’s protocol.
MSC mRNA sequencing and data analysis
To explore MSC angiogenic regulation in DKD, high-throughput RNA sequencing was performed as previously reported . Briefly, RNA libraries were generated from 1 µg of total RNAs using an Illumina TruSeq RNA Sample Prep Kit v2 and loaded onto flow cells (8–10 pM) to generate cluster densities of 700,000/mm2. MSC were sequenced on an Illumina HiSeq 2000. Comprehensive analysis including alignment, gene and exon quantification, and QC of raw RNA sequencing reads was performed using the Mayo Analysis Pipeline for RNA Seq (MAP-RSeq) version 3.1.3 , an in-house bioinformatics workflow. R-bioinformatics package (edgeR) was used for differential gene expression analysis between MSC from DKD and control participants. Genes with p adjusted value < 0.05 and |log2FC|> 1 were considered significantly differentially expressed (DE) genes. GeneCards® database (http://www.genecards.org/) was utilized to screen genes associated with angiogenesis. DE angiogenesis mRNAs-of-interest were visualized on volcano plots and heatmaps created using R package ggplot2.
MSC miRNA sequencing and data analysis
Secondary analyses of high-throughput miRNA analyses on MSC were performed utilizing CAP-miRSeq version 1.1 workflow . CAP-miRSeq generates raw and normalized expression counts for both known mature miRNAs and predicted novel miRNAs from unaligned reads (FASTQ). Differential miRNA expression analysis was performed using bioinformatics R package edgeR to identify miRNA upregulated in DKD-MSC compared to Control-MSC (|log2FC|> 1 | and Benjamini–Hochberg adjusted p value < 0.05).
Integrated mRNA/miRNA analysis
The miRNA/mRNA interactions were performed by crossing DE MSC mRNAs and miRNAs to predict targets analyzed using microRNA-Target analysis module in Ingenuity Pathway Analysis software (Ingenuity® Systems, www.ingenuity.com).
Validation of mRNA and miRNA by quantitative polymerase chain reaction (qPCR)
To validate the expression levels of representative upregulated and down-regulated mRNAs and miRNAs, qPCR was performed on MSC. RNA was isolated from MSC (0.5–1.0 × 106 cells) and ran using an Applied Biosystems ViiA7 Real-Time PCR system as previously mentioned . Fold change of gene expression was calculated using 2−ΔΔCT method. All probes were from Thermo Fisher Scientific (bone morphogenetic protein 2/BMP2: Hs00154192, proenkephalin/PENK: Hs00175049, vascular cell adhesion molecule 1/VCAM1: Hs01003372, insulin-like growth factor-binding protein 2/IGFBP2: Hs01040719_m1, thrombospondin 1/THBS1: Hs00962908, integrin subunit beta 8/ITGB8: Hs00174456, glyceraldehyde 3-phosphate dehydrogenase/GAPDH: Hs02786624, miR-let-7a-5p: hsa-let-7a-5p, miR-30c-5p: hsa-miR-30c and U6: U6 small nuclear RNA/snRNA, 715680. mRNAs and miRNAs were normalized to GAPDH and U6, respectively.
MSC angiogenic activity
Angiogenic potential of MSC was assessed in vitro. Briefly, commercially available human umbilical vein endothelial cells (HUVEC, 200K-05f, Cell Applications Inc., San Diego, CA) were grown in endothelial cell growth medium [20, 31]. HUVEC were divided into 4 groups for experimentation: group 1 cells were cultured under normal conditions (non-injured), while group 2–4 cells were co-incubated with high glucose (HG, 25 mmol) and indoxyl sulfate (IS, 1 mmol) for 6 h to simulate a uremic, diabetic milieu. Groups 3 and 4 HUVEC were simultaneously incubated with Control-MSCcm or DKD-MSCcm (medium from 1.3 × 105 cells), respectively, at the time of HG + IS injury. After 6 h, all HUVEC underwent functional assessments in vitro (migration, proliferation, or tube formation assay) or were lysed and prepared for qPCR analyses.
HUVEC (1.6 × 105/well) were seeded onto a 6-well plate and treated as described above. To assess proliferation of HUVEC in groups 1–4, the cells were fixed in 4% paraformaldehyde for 15 min, then stained for Ki67 (Abcam, ab15580) according to the manufacturers protocol. Treated HUVEC were also assayed for migration capacity via the QCM 24-well Colorimetric Cell Migration (Millipore, ECM 508) following manufacturer’s protocol. Samples were observed and imaged using an EVOS M5000 imaging system and analyzed via Image J. HUVEC (1.4 × 104/well) were also seeded onto a 96-well plate coated with GelTrex™ matrix (A1413201; Thermo Fisher Scientific) to assess tube formation. Plates were placed into a 37 °C incubator with 5% CO2. After 4 h, HUVEC were stained with calcein-AM (10 µg/mL, for 30 min at 37 °C) to assess capillary-like structures and observed under an EVOS M5000 imaging system (Objective 40×; lens EVOS_AMEP4683; Camera 3.2 megapixels, monochrome, CMOS; Detectors and filters model DAPI—AMEP4650, F1420-1404-0373, TX Red—AMEP4655, F1620-1475-0197). A total number of tubular structures were manually counted. Gene expression was evaluated in HUVEC using the following probes from Thermo Fisher Scientific (E-selectin/SELE: Hs00174057, VCAM1, intercellular adhesion molecule 1/ICAM-1: Hs00164932, VEGF: HS00900055_m1, and THBS1:HS00962908_m1). Expression was normalized to TATA-binding protein (TBP: HS00427620).
Data analysis was performed using GraphPad Prism 9 statistical software. Shapiro–Wilk normality test was performed on data from each group before statistical evaluation. Normally distributed data were represented as mean ± standard deviation, and non-normal data as median and interquartile range. Comparisons between Control-MSC and DKD-MSC were performed using either a two-sample t test with a 5% type-I error rate or Wilcoxon Rank Sum, as appropriate. One-way analysis of variance (ANOVA) and post hoc pairwise testing were employed to detect differences in co-culture experiments. Statistical significance was accepted at p ≤ 0.05.
In total, 38 participants were included (n = 9 controls and n = 29 DKD subjects). Clinical characteristics of each cohort are summarized in Table 1. No significant differences were evident between groups for age, sex, or race. Body mass index (BMI) was higher (35.4 ± 5.6 vs. 28.8 ± 3.6 kg/m2; p < 0.0001), and by design, eGFR was lower in DKD participants (38.9 ± 15.4 vs. 80.5 ± 13.3 mL/min/1.73 m2; p < 0.0001) compared to non-diabetic, non-CKD control subjects.
AD-MSC express putative surface markers
Undifferentiated AD-MSC obtained from 38 study participants produce a relatively homogenous cell population exhibiting plastic adherence, fibroblast-like morphology, and multi-lineage potential as previously shown . Culture-expanded MSC demonstrated robust MSC-specific cell surface positivity to CD90, CD105, and CD73, while low expression of hematopoietic markers, CD34 and CD14, was observed in MSC from control and DKD subjects (Fig. 1). Flow cytometric analyses revealed that MSC from DKD subjects maintain surface marker characteristics meeting the criteria required for MSC .
mRNA alterations in DKD-MSC
A total of 4611 mRNAs associated with angiogenesis were identified from the GeneCards® database (https://www.genecards.org/). Of these mRNAs, a total of 133 DE genes were related to angiogenesis. The distribution of DE upregulated and down-regulated genes is displayed in a volcano plot (Fig. 2A). Among them, 77 mRNAs were upregulated (Fig. 2B) and 56 down-regulated in DKD-MSC (Fig. 2C). Among the top 15 DE upregulated mRNAs were pro-angiogenic genes, including IGFBP2, lysyl oxidase-like 4 (LOXL4) and VCAM1, and anti-angiogenic genes, including THBS1 and ITGB8. Among the top 15 DE down-regulated genes several were pro-angiogenic genes, including interleukin (IL)13 receptor subunit alpha-2, C-X-C motif chemokine ligand (CXCL)3, IL33, PENK, matrix metalloproteinase (MMP)10, CXCL8, BMP2, and MMP3 as well as an anti-angiogenic gene, MMP12. Expression patterns of several mRNAs were subsequently confirmed by qPCR (Additional file 1: Fig. S1A and C). In addition, DKD-MSCcm tended to release higher levels of the anti-angiogenic protein TSP1 compared to Control-MSCcm (Additional file 1: Fig. S1B).
miRNA alterations in DKD-MSC
A total of 575 annotated miRNAs were identified, of which 208 (119 upregulated and 89 down-regulated) were DE in DKD-MSC. Expression patterns of miRNAs were subsequently confirmed by qPCR (Additional file 1: Fig. S1D). DE miRNAs were further inspected to identify mRNA targets involved in angiogenesis. Target prediction analysis resulted in 14 unique DE miRNAs regulating 18 unique DE angiogenesis-related genes (Fig. 3A) and 25 mRNA-miRNA interactions (Fig. 3B). Surprisingly, only one miRNA was down-regulated (miR-148a-3p), which interacted with ribosomal protein S6 kinase A5 (RPS6KA5), involved in phosphorylating cAMP-response element-binding protein (CREB) that promotes angiogenesis by regulating vascular endothelial growth factor (VEGF) receptor-1(VEGFR1). Notably, let-7a-5p targeted five mRNAs (transgelin/TAGLN, THBS1, LOXL4, collagen 4A1/COL4A1, and collagen 8A1/COL8A1; Fig. 3C).
DKD-MSCcm angiogenic functionality
Apart from the influence of hyperglycemia, uremic conditions present in DKD may further impact gene expression, intrinsic properties of MSC, and MSC functionality. Therefore, we assessed the reparative capacity of DKD-MSCcm on injured endothelial cells by establishing a uremic and diabetic microenvironment using HG + IS. Pre-incubation of HUVEC with HG + IS diminished tube formation capacity. Co-incubation of HUVEC with DKD-MSCcm prior to HG + IS exposure restored formation of tube-like networks (Fig. 4A, B), reflecting sustained pro-angiogenic properties. In other functional studies, incubation of DKD-MSCcm restored the migration potential of injured HUVEC, while proliferation activity was not different between groups (Additional file 2: Fig. S2A, B). Co-culture of MSCcm with injured HUVEC increased THBS1 mRNA expression, though DKD-MSCcm exhibited lower THBS1 levels when compared with Control-MSCcm (p = 0.07) (Additional file 2: Fig. S2C). Furthermore, exposure to HG + IS led to upregulation of adhesion molecules SELE, VCAM1, ICAM-1 in HUVEC, while incubation of DKD-MSCcm attenuated mRNA expression, suggesting that DKD-MSCcm prevents endothelial dysfunction (Fig. 4C). No differences were observed between DKD-MSCcm and Control-MSCcm groups. Furthermore, VEGF mRNA levels increased in HG + IS-exposed HUVEC, but MSCcm from DKD and Controls prevented the rise in VEGF expression (Fig. 4D).
Our study used comprehensive sequencing analysis to compare the enrichment and suppression of mRNA and miRNA expression in adipose tissue-derived MSC in DKD and control subjects. We further explored the pro-angiogenic potential of DKD-MSCcm in vitro. Compared to MSC from controls, DKD-MSC exhibited heterogeneity in the expression profiles of mRNA associated with angiogenic activity. Furthermore, 14 miRNAs targeted these angiogenesis-related genes modifying mRNA activity. In particular, miRNA let-7a-5p, which regulates angiogenesis, has anti-thrombospondin activity, and is associated with DKD pathogenesis, was upregulated and had the highest number of mRNA targets. In our in vitro studies mimicking the DKD microenvironment, DKD-MSCcm incubation restored endothelial cell function and halted overexpression of adhesion molecules associated with vascular insult. These observations suggest that despite alterations of the transcriptome, the pro-angiogenic reparative capacity of MSC in vitro appears relatively preserved in DKD subjects.
To limit potential immunological reactions and minimize the risk of allosensitization in subjects who may eventually desire kidney transplantation, an autologous MSC source is preferable. However, DKD patient factors such as advanced age, hyperglycemia, and uremia may diminish MSC reparative ability. Our previous ex vivo studies in MSC from patients with DKD showed that DKD-MSC anti-inflammatory, immunomodulatory, anti-apoptosis, and anti-fibrosis functionality remained intact compared with Control-MSC. However, the migratory capacity of DKD-MSC was lower than that of Control-MSC . Therefore, we further explored functional differences in angiogenic capacity between DKD-MSC and Control-MSC in this study.
Angiogenic activity is dependent on the intricate balance between pro-angiogenic factors (such as VEGF, fibroblast growth factor-2/FGF2, transforming growth factor-β/TGF-β, and angiopoietins) and factors that inhibit angiogenesis (angiostatin, endostatin, thrombospondins) . Diabetes is associated with angiogenesis dysregulation leading to structurally immature blood vessel formation after injury, impaired wound healing, and excessive angiogenesis—particularly in DKD and diabetic retinopathy [33, 34]. We previously found no difference in secretion of VEGF-A between DKD-MSC and Control-MSC, while hepatocyte growth factor (HGF), maintaining anti-apoptosis, anti-fibrotic, and pro-angiogenic activity, was substantially higher in DKD-MSC . VEGF is the key mediator of angiogenesis and vessel repair. Xin et al.  found that the combination of HGF and VEGF increased neovascularization greater than either growth factor alone. In our studies, DKD-MSCcm restored HUVEC formation of tube-like networks and migration ability in a uremic, diabetic milieu, thus demonstrating that DKD-MSC maintain angiogenic repair capability in vitro. Previous studies identified MSC subsets possessing greater angiogenic paracrine activity. Du et al.  found that VCAM1+ expressing MSC had increased gene expression and release of pro-angiogenic factors and cytokines (VEGF, HGF, MMP2, and CXCL1) compared to VCAM1−-MSC. In in vivo studies, exogenous delivery of VCAM1+-MSC increased perfusion in the preclinical ischemic hindlimb model versus VCAM1−-MSC. Likewise, VCAM1 mRNA was upregulated in our DKD-MSC, which could potentially influence pro-angiogenic repair activity observed in vitro.
Contrarily, anti-angiogenic genes were upregulated (e.g. THBS1) and pro-angiogenic genes were down-regulated (e.g. BMP2) in DKD-MSC in our study, which could disturb angiogenic activity. Thrombospondin 1 (TSP1), encoded by the THBS1 gene, is a secreted anti-angiogenic protein that modulates cell migration and adhesion by regulating vascular nitric oxide signaling [37, 38]. In a study by Dzhoyashvili et al. , MSC harvested from subjects with type 2 diabetes and coronary artery disease showed impaired ability to stimulate angiogenesis in vitro with significant increase in mRNA level of THBS1, as well as closely negative correlation between its mRNA level and MSC angiogenic activity measured by tube formation assay. Similar to these findings, our MSC harvested from participants with DKD also had trends of higher THBS1 mRNA expression and protein release in condition media, but no statistical differences compared with Control-MSC for angiogenic capability to restore formation of tube-like networks in injured endothelial cells. Similarly, we found that DKD-MSCcm led to increased migration of HG + IS injured HUVEC versus injured controls with no difference observed between DKD- and Control-MSCcm groups. BMP2 mRNA expression was lower in DKD-MSC in our study. BMP2 has been proven to promote angiogenesis by facilitating endothelial migration , invasion, and proliferation  and induce HUVEC tube formation . Overall, these observations suggest that the ultimate functional outcome may depend on a weighted effect of these genes and their inhibitors, in concert with the impact of the microenvironment.
Given that miRNA regulate gene activity, we further assessed the relationship between miRNA and angiogenesis-related mRNA in DKD-MSC. In our study, 208 differentially expressed miRNAs were identified among which 14 miRNAs targeted differentially expressed angiogenesis-related mRNA. Interestingly, let-7a-5p was upregulated in RNA-seq analysis (and had a trend toward higher gene expression, supplemental Fig. 1D) and targets the highest number of mRNAs, including TAGLN, THBS1, LOXL4, COL4A1, and COL8A1. Let-7a-5p has pro-angiogenic properties with increasing sprout formation of hypoxic HUVEC [43, 44] by directly targeting the 3′UTR of the THBS1 mRNA . Hence, we speculate that DKD-MSC restore angiogenetic abilities partly because of increased let-7a-5p. Zhu et al.  conducted miRNA sequencing revealing that let-7 family (let-7a-5p) was enriched in EVs derived from human AD-MSC, and blocking of let-7 impaired tube formation of HUVEC in vitro. In addition, platelet-derived EVs can deliver let-7a into HUVEC in vitro with inhibition of HUVEC tube formation by down-regulating TSP1 protein expression . Consistent with these findings, we used DKD-MSCcm, containing MSC-EVs and found trends of higher tube formation and lower THBS1 mRNA expression in HUVEC compared to Control-MSCcm. Further investigation is required to explore the mechanism and outcomes of altered let-7a-5p in vivo and in vitro. In addition, our studies found that TSP1 targets included: TAGLN, a regulator of angiogenesis , LOXL4, a contributor to the vascular permeability in diabetes , types of collagen (COL4A1 and COL8A1) that influence angiogenesis plus THSB1 an inhibitor of angiogenesis . Further observations are needed to explore how let-7a-5p adjusts these mRNA. Taken together, these observations suggest that differentially expressed miRNAs, particularly let-7a-5p, in DKD-MSC might contribute to angiogenesis dysregulation.
We further examined MSC angiogenic repair capacity in vitro. After an in vitro injury mimicking a DKD microenvironment, endothelial cells lost tube formation capacity and increased adhesion molecule (SELE, VCAM, ICAM-1) expression, key molecules associated in early vascular injury [49,50,51]. These effects were reversed following co-incubation with DKD-MSCcm. Additionally, HUVEC migration was reduced after HG + IS treatment but improved with DKD-MSCcm co-culture but was not different versus Control-MSCcm. Angiogenic function impairment in vivo was shown by Kim et al.  wherein exogenous delivery of MSC from diabetic rats failed to improve blood flow compared to MSC from non-hyperglycemic rats in a model of hindlimb ischemia. Interestingly, higher expression levels of pro-angiogenic factors (HGF, FGF2, PGF, insulin-like growth factor-1, and angiopoietin-2 were evident in both control- and diabetic MSC infused into ischemic tissue 2 weeks-post administration, indicating that diabetic MSC still possess paracrine angiogenic repair capacities albeit limited. Hence, a therapeutic boost of autologous MSC infusion could still help enhance angiogenic repair capacity in patients with DKD.
This study has limitations. First, the predominantly Caucasian and relatively older study cohort may limit generalizability of our findings to other races and younger individuals with DKD, but remains reflective of the DKD population in the USA . Secondly, AD-MSC may differ from other autologous cell sources, such as bone marrow or umbilical cord , though AD-MSC may possess more potent pro-angiogenic activity . Furthermore, in vivo studies are needed to confirm these in vitro findings in DKD-MSC. Finally, future mechanistic studies employing loss of function (e.g., siRNA or miRNA inhibitors) are necessary to establish which DE genes are involved in modification of DKD-MSC angiogenic function.
Compared to Control-MSC, AD-MSC harvested from DKD participants possess modified genetic messages related to angiogenesis. Yet, DKD-MSC maintained the intrinsic capacity to release pro-angiogenic paracrine factors to avert further endothelial dysfunction. Identification of disease-related dysfunction and interindividual variation in MSC may allow for development of regenerative strategies to improve MSC function and thwart barriers that prevent successful MSC therapeutic outcomes in individuals with DKD.
Availability of data and materials
The datasets supporting the conclusions of this article are available from the corresponding author on reasonable request, without undue reservation. Sequencing data were deposited in NCBI’s Gene Expression Omnibus and are publicly available [accession number: GSE217712].
Adipose tissue-derived mesenchymal stem cells
Analysis of variance
Bone morphogenetic protein 2
Body mass index
Cluster of differentiation
Chronic kidney disease
Control mesenchymal stromal cells
CAMP-response element-binding protein
C–X–C motif chemokine ligand
Diabetic kidney disease
Diabetic kidney disease mesenchymal stromal cells
Diabetic kidney disease mesenchymal stromal cells condition medium
Estimated glomerular filtration rate
Fibroblast growth factor-2
Glyceraldehyde 3-phosphate dehydrogenase
Hepatocyte growth factor
Human umbilical vein endothelial cells
Intercellular adhesion molecule 1
Insulin-like growth factor-binding protein 2
Interleukin 13 receptor subunit alpha-2
Integrin subunit beta 8
Lysyl oxidase-like 4
Mayo analysis pipeline for RNA-seq
Mesenchymal stromal cells
Mesenchymal stromal cells condition medium
Adjusted p value
Plasminogen activator inhibitor-1
Placental growth factor
Ribosomal protein S6 kinase A5
Small nuclear RNA
Transforming growth factor-β
Thrombospondin 1 (gene)
Thrombospondin 1 (protein)
Vascular cell adhesion molecule 1
Vascular endothelial growth factor
Vascular endothelial growth factor receptor-1
Saran R, Robinson B, Abbott KC, Bragg-Gresham J, Chen X, Gipson D, et al. US renal data system 2019 annual data report: epidemiology of kidney disease in the United States. Am J Kidney Dis. 2020;75(1 Suppl 1):A6–7.
Kim SR, Zou X, Tang H, Puranik AS, Abumoawad AM, Zhu XY, et al. Increased cellular senescence in the murine and human stenotic kidney: effect of mesenchymal stem cells. J Cell Physiol. 2021;236(2):1332–44.
Yu S, Kim SR, Jiang K, Ogrodnik M, Zhu XY, Ferguson CM, et al. Quercetin reverses cardiac systolic dysfunction in mice fed with a high-fat diet: role of angiogenesis. Oxid Med Cell Longev. 2021;2021:8875729.
Bian X, Griffin TP, Zhu X, Islam MN, Conley SM, Eirin A, et al. Senescence marker activin A is increased in human diabetic kidney disease: association with kidney function and potential implications for therapy. BMJ Open Diabetes Res Care. 2019;7(1):e000720.
Hickson LJ, Langhi Prata LGP, Bobart SA, Evans TK, Giorgadze N, Hashmi SK, et al. Senolytics decrease senescent cells in humans: preliminary report from a clinical trial of Dasatinib plus Quercetin in individuals with diabetic kidney disease. EBioMedicine. 2019;47:446–56.
Palmer AK, Xu M, Zhu Y, Pirtskhalava T, Weivoda MM, Hachfeld CM, et al. Targeting senescent cells alleviates obesity-induced metabolic dysfunction. Aging Cell. 2019;18(3):e12950.
Hickson LJ, Rule AD, Thorsteinsdottir B, Shields RC, Porter IE, Fleming MD, et al. Predictors of early mortality and readmissions among dialysis patients undergoing lower extremity amputation. J Vasc Surg. 2018;68(5):1505–16.
Geiss LS, Li Y, Hora I, Albright A, Rolka D, Gregg EW. Resurgence of diabetes-related nontraumatic lower-extremity amputation in the young and middle-aged adult U.S. population. Diabetes Care. 2019;42(1):50–4.
Anantha-Narayanan M, Sheikh AB, Nagpal S, Jelani QU, Smolderen KG, Regan C, et al. Systematic review and meta-analysis of outcomes of lower extremity peripheral arterial interventions in patients with and without chronic kidney disease or end-stage renal disease. J Vasc Surg. 2021;73(1):331–40.
Eirin A, Zhu XY, Krier JD, Tang H, Jordan KL, Grande JP, et al. Adipose tissue-derived mesenchymal stem cells improve revascularization outcomes to restore renal function in swine atherosclerotic renal artery stenosis. Stem Cells. 2012;30(5):1030–41.
Liew A, O’Brien T. Therapeutic potential for mesenchymal stem cell transplantation in critical limb ischemia. Stem Cell Res Ther. 2012;3(4):28.
Ferguson CM, Farahani RA, Zhu XY, Tang H, Jordan KL, Saadiq IM, et al. Mesenchymal stem/stromal cell-derived extracellular vesicles elicit better preservation of the intra-renal microvasculature than renal revascularization in pigs with renovascular disease. Cells. 2021;10(4):763.
Zhu XY, Klomjit N, Conley SM, Ostlie MM, Jordan KL, Lerman A, et al. Impaired immunomodulatory capacity in adipose tissue-derived mesenchymal stem/stromal cells isolated from obese patients. J Cell Mol Med. 2021;25:9051–9.
Hickson LJ, Abedalqader T, Ben-Bernard G, Mondy JM, Bian X, Conley SM, et al. A systematic review and meta-analysis of cell-based interventions in experimental diabetic kidney disease. Stem Cells Transl Med. 2021;10(9):1304–19.
Zou X, Jiang K, Puranik AS, Jordan KL, Tang H, Zhu X, et al. Targeting murine mesenchymal stem cells to kidney injury molecule-1 improves their therapeutic efficacy in chronic ischemic kidney injury. Stem Cells Transl Med. 2018;7(5):394–403.
Packham DK, Fraser IR, Kerr PG, Segal KR. Allogeneic mesenchymal precursor cells (MPC) in diabetic nephropathy: a randomized, placebo-controlled, dose escalation study. EBioMedicine. 2016;12:263–9.
Bura A, Planat-Benard V, Bourin P, Silvestre JS, Gross F, Grolleau JL, et al. Phase I trial: the use of autologous cultured adipose-derived stroma/stem cells to treat patients with non-revascularizable critical limb ischemia. Cytotherapy. 2014;16(2):245–57.
Hickson LJ, Herrmann SM, McNicholas BA, Griffin MD. Progress toward the clinical application of mesenchymal stromal cells and other disease-modulating regenerative therapies: examples from the field of nephrology. Kidney360. 2021;2(3):542–57.
Isik B, Thaler R, Goksu BB, Conley SM, Al-Khafaji H, Mohan A, et al. Hypoxic preconditioning induces epigenetic changes and modifies swine mesenchymal stem cell angiogenesis and senescence in experimental atherosclerotic renal artery stenosis. Stem Cell Res Ther. 2021;12(1):240.
Conley SM, Hickson LJ, Kellogg TA, McKenzie T, Heimbach JK, Taner T, et al. Human obesity induces dysfunction and early senescence in adipose tissue-derived mesenchymal stromal/stem cells. Front Cell Dev Biol. 2020;8:197.
Hickson LJ, Eirin A, Conley SM, Taner T, Bian X, Saad A, et al. Diabetic kidney disease alters the transcriptome and function of human adipose-derived mesenchymal stromal cells but maintains immunomodulatory and paracrine activities important for renal repair. Diabetes. 2021;70(7):1561–74.
Reinders ME, Roemeling-van Rhijn M, Khairoun M, Lievers E, de Vries DK, Schaapherder AF, et al. Bone marrow-derived mesenchymal stromal cells from patients with end-stage renal disease are suitable for autologous therapy. Cytotherapy. 2013;15(6):663–72.
van Rhijn-Brouwer FCC, van Balkom BWM, Papazova DA, Hazenbrink DHM, Meijer AJ, Brete I, et al. Paracrine proangiogenic function of human bone marrow-derived mesenchymal stem cells is not affected by chronic kidney disease. Stem Cells Int. 2019;2019:1232810.
Abumoawad A, Saad A, Ferguson CM, Eirin A, Herrmann SM, Hickson LJ, et al. In a phase 1a escalating clinical trial, autologous mesenchymal stem cell infusion for renovascular disease increases blood flow and the glomerular filtration rate while reducing inflammatory biomarkers and blood pressure. Kidney Int. 2020;97(4):793–804.
Saad A, Dietz AB, Herrmann SMS, Hickson LJ, Glockner JF, McKusick MA, et al. Autologous mesenchymal stem cells increase cortical perfusion in renovascular disease. J Am Soc Nephrol JASN. 2017;28(9):2777–85.
El-Badawy A, Amer M, Abdelbaset R, Sherif SN, Abo-Elela M, Ghallab YH, et al. Adipose stem cells display higher regenerative capacities and more adaptable electro-kinetic properties compared to bone marrow-derived mesenchymal stromal cells. Sci Rep. 2016;6:37801.
Adolfsson E, Helenius G, Friberg O, Samano N, Frobert O, Johansson K. Bone marrow- and adipose tissue-derived mesenchymal stem cells from donors with coronary artery disease; growth, yield, gene expression and the effect of oxygen concentration. Scand J Clin Lab Investig. 2020;80(4):318–26.
Crespo-Diaz R, Behfar A, Butler GW, Padley DJ, Sarr MG, Bartunek J, et al. Platelet lysate consisting of a natural repair proteome supports human mesenchymal stem cell proliferation and chromosomal stability. Cell Transplant. 2011;20(6):797–811.
Kalari KR, Nair AA, Bhavsar JD, O’Brien DR, Davila JI, Bockol MA, et al. MAP-RSeq: mayo analysis pipeline for RNA sequencing. BMC Bioinform. 2014;15:224.
Sun Z, Evans J, Bhagwate A, Middha S, Bockol M, Yan H, et al. CAP-miRSeq: a comprehensive analysis pipeline for microRNA sequencing data. BMC Genom. 2014;15:423.
Yigitbilek F, Conley SM, Tang H, Saadiq IM, Jordan KL, Lerman LO, et al. Comparable in vitro function of human liver-derived and adipose tissue-derived mesenchymal stromal cells: implications for cell-based therapy. Front Cell Dev Biol. 2021;9:641792.
Dominici M, Le Blanc K, Mueller I, Slaper-Cortenbach I, Marini F, Krause D, et al. Minimal criteria for defining multipotent mesenchymal stromal cells. The International Society for Cellular Therapy position statement. Cytotherapy. 2006;8(4):315–7.
Tahergorabi Z, Khazaei M. Imbalance of angiogenesis in diabetic complications: the mechanisms. Int J Prev Med. 2012;3(12):827–38.
Sorrells WS, Mao SA, Taner T, Jadlowiec CC, Farres H, Davila V, et al. Endarterectomy for iliac occlusive disease during kidney transplantation: a multicenter experience. Int J Angiol. 2021;30(2):91–7.
Xin X, Yang S, Ingle G, Zlot C, Rangell L, Kowalski J, et al. Hepatocyte growth factor enhances vascular endothelial growth factor-induced angiogenesis in vitro and in vivo. Am J Pathol. 2001;158(3):1111–20.
Du W, Li X, Chi Y, Ma F, Li Z, Yang S, et al. VCAM-1+ placenta chorionic villi-derived mesenchymal stem cells display potent pro-angiogenic activity. Stem Cell Res Ther. 2016;7:49.
Isenberg JS, Ridnour LA, Perruccio EM, Espey MG, Wink DA, Roberts DD. Thrombospondin-1 inhibits endothelial cell responses to nitric oxide in a cGMP-dependent manner. Proc Natl Acad Sci USA. 2005;102(37):13141–6.
Isenberg JS, Martin-Manso G, Maxhimer JB, Roberts DD. Regulation of nitric oxide signalling by thrombospondin 1: implications for anti-angiogenic therapies. Nat Rev Cancer. 2009;9(3):182–94.
Dzhoyashvili NA, Efimenko AY, Kochegura TN, Kalinina NI, Koptelova NV, Sukhareva OY, et al. Disturbed angiogenic activity of adipose-derived stromal cells obtained from patients with coronary artery disease and diabetes mellitus type 2. J Transl Med. 2014;12:337.
Garcia de Vinuesa A, Abdelilah-Seyfried S, Knaus P, Zwijsen A, Bailly S. BMP signaling in vascular biology and dysfunction. Cytokine Growth Factor Rev. 2016;27:65–79.
Dyer LA, Pi X, Patterson C. The role of BMPs in endothelial cell function and dysfunction. Trends Endocrinol Metab. 2014;25(9):472–80.
Finkenzeller G, Hager S, Stark GB. Effects of bone morphogenetic protein 2 on human umbilical vein endothelial cells. Microvasc Res. 2012;84(1):81–5.
Aday S, Hazan-Halevy I, Chamorro-Jorganes A, Anwar M, Goldsmith M, Beazley-Long N, et al. Bioinspired artificial exosomes based on lipid nanoparticles carrying let-7b-5p promote angiogenesis in vitro and in vivo. Mol Ther. 2021;29(7):2239–52.
Zhu Y, Zhang J, Hu X, Wang Z, Wu S, Yi Y. Extracellular vesicles derived from human adipose-derived stem cells promote the exogenous angiogenesis of fat grafts via the let-7/AGO1/VEGF signalling pathway. Sci Rep. 2020;10(1):5313.
Anene C, Graham AM, Boyne J, Roberts W. Platelet microparticle delivered microRNA-Let-7a promotes the angiogenic switch. Biochim Biophys Acta Mol Basis Dis. 2018;1864(8):2633–43.
Tsuji-Tamura K, Morino-Koga S, Suzuki S, Ogawa M. The canonical smooth muscle cell marker TAGLN is present in endothelial cells and is involved in angiogenesis. J Cell Sci. 2021. https://doi.org/10.1242/jcs.25492.
Song B, Kim D, Nguyen NH, Roy S. Inhibition of diabetes-induced lysyl oxidase overexpression prevents retinal vascular lesions associated with diabetic retinopathy. Investig Ophthalmol Vis Sci. 2018;59(15):5965–72.
Zhou L, Isenberg JS, Cao Z, Roberts DD. Type I collagen is a molecular target for inhibition of angiogenesis by endogenous thrombospondin-1. Oncogene. 2006;25(4):536–45.
Qiu S, Cai X, Liu J, Yang B, Zugel M, Steinacker JM, et al. Association between circulating cell adhesion molecules and risk of type 2 diabetes: a meta-analysis. Atherosclerosis. 2019;287:147–54.
Siddiqui K, George TP, Nawaz SS, Joy SS. VCAM-1, ICAM-1 and selectins in gestational diabetes mellitus and the risk for vascular disorders. Future Cardiol. 2019;15(5):339–46.
Zhito AV, Iusupova AO, Kozhevnikova MV, Shchendrygina AA, Privalova EV, Belenkov YN. E-selectin as a marker of endothelial dysfunction in patients with coronary artery disease including those with type 2 diabetes mellitus. Kardiologiia. 2020;60(4):24–30.
Kim H, Han JW, Lee JY, Choi YJ, Sohn YD, Song M, et al. Diabetic mesenchymal stem cells are ineffective for improving limb ischemia due to their impaired angiogenic capability. Cell Transplant. 2015;24(8):1571–84.
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. Next generation sequencing data analysis. EMBnetjournal. 2011;17(1):10. https://doi.org/10.14806/ej.17.1.200.
Russell AL, Lefavor R, Durand N, Glover L, Zubair AC. Modifiers of mesenchymal stem cell quantity and quality. Transfusion. 2018;58(6):1434–40.
We would like to thank the following individuals for their contributions: Shahrzad Tehranian MD (data collection); Andre van Wijnen, PhD (overall RNA-seq advisement) Sunil Krishnan, MD and Kartik Lakhotiya (cell studies); Zhuo Li, MS (statistical support); Caroline Skaff (revision analyses), and Mrunal Dehankar, MS, Pritha Chana, MS, Asha A. Nair, PhD, and Jaime I. Davila, PhD (Bioinformatics Core; RNA Sequencing).
This research was funded by Burroughs Wellcome Fund Postdoctoral Enrichment Program (SMC); Mayo Clinic Robert D. and Patricia E. Kern Center for the Science of Health Care Delivery (LJH); NIDDK Diabetes Complications Consortium, Grant/Award Numbers: DK115255 and DK076169 (LJH); National Institutes of Health (NIH), Grant/Award Numbers: DK118120 (SMH), DK106427 (AE), DK122137 (AE), DK123492 (LJH), DK109134 (LJH); DK120292 (LOL); Mayo Clinic Florida Center for Regenerative Medicine; and Regenerative Medicine Minnesota, Grant/Award Number: RMM 091718CT (LJH). The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
The studies involving human participants were reviewed and approved by Mayo Clinic Institutional Review Board (IRB: 15-000933; “Stem Cells in Chronic Kidney Disease”) in May 2015. This study adheres to the Declaration of Helsinki. Study participants provided their written informed consent to participate in this study.
Consent for publication
Lilach Lerman, MD, PhD, received grant funding from Novo Nordisk and is an advisor to Weijian Technologies and AstraZeneca. The remaining authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Validation of upregulated (A) and down-regulated (C) differentially expressed (DE) messenger RNAs (mRNAs) in diabetic kidney disease (DKD)-mesenchymal stromal cells (MSC), as well as microRNAs (miRNAs) (D). ELISA of TSP1 in MSC conditioned medium (MSCcm) (B). BMP2, bone morphogenetic protein 2; PENK, proenkephalin; VCAM1, vascular cell adhesion molecule 1; IGFBP2, insulin-like growth factor-binding protein 2; THBS1/TSP1, thrombospondin 1; ITGB8, integrin subunit beta 8.
Representative images of capillary-like tubes formed by non-injured human umbilical vein endothelial cells (HUVEC) or in the presence of high glucose (HG) plus indoxyl sulfate (IS) injury and either Control-MSC conditioned medium (cm) or DKD-MSCcm. Images were acquired at 40X resolution and were not enhanced. (A). 4′,6-diamidino-2-phenylindole (DAPI) nuclear DNA (A). Nuclear staining (blue), proliferation marker Ki67 protein staining (red). Quantification of proliferation ability (A) and migratory function (B). Thrombospondin (THSB1) gene expression in HUVEC groups (C).
About this article
Cite this article
Bian, X., Conley, S.M., Eirin, A. et al. Diabetic kidney disease induces transcriptome alterations associated with angiogenesis activity in human mesenchymal stromal cells. Stem Cell Res Ther 14, 49 (2023). https://doi.org/10.1186/s13287-023-03269-9