The whole profiling and competing endogenous RNA network analyses of noncoding RNAs in adipose-derived stem cells from diabetic, old, and young patients
Stem Cell Research & Therapy volume 12, Article number: 313 (2021)
Mesenchymal stem cells including adipose-derived stem cells (ASCs) have a considerable potential in the field of translational medicine. Unfortunately, multiple factors (e.g., older age, co-existing diabetes, and obesity) may impair cellular function, which hinders the overall effectiveness of autologous stem cell therapy. Noncoding RNAs—including microRNAs (miRNAs), long ncRNAs (lncRNAs), and circular RNAs (circRNAs)—have been shown to play important roles in stem cell biology. However, the overall diabetes-related and aging-related expression patterns and interactions of these RNAs in ASCs remain unknown.
The phenotypes and functions of ASCs isolated from diabetic (D-ASCs), old (O-ASCs), and young (Y-ASCs) donors were evaluated by in vitro assays. We conducted high-throughput RNA sequencing (RNA-seq) in these ASCs to identify the differentially expressed (DE) RNAs. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, and protein-protein interaction (PPI) analyses were performed to investigate mRNAs with significant differences among groups. The lncRNA- or circRNA-associated competing endogenous RNA (ceRNA) networks were constructed based on bioinformatics analyses and real-time polymerase chain reaction (RT-PCR) results. The miR-145-5p mimics were transfected into O-ASCs and verified by PCR.
ASCs from diabetic and old donors showed inferior migration ability and increased cellular senescence. Furthermore, O-ASCs have decreased capacities for promoting endothelial cell angiogenesis and fibroblast migration, compared with Y-ASCs. The DE miRNAs, mRNAs, lncRNAs, and circRNAs were successfully identified by RNA-seq in O-ASCs vs. Y-ASCs and D-ASCs vs. O-ASCs. GO and KEGG analyses demonstrated that DE mRNAs were significantly enriched in aging and cell senescence terms separately. PPI networks revealed critical DE mRNAs in the above groups. RNAs with high fold changes and low p values were validated by PCR. ceRNA networks were constructed based on bioinformatics analyses and validated RNAs. Additionally, the lncRNA RAET1E-AS1–miR-145-5p–WNT11/BMPER axis was validated by PCR and correlation analyses. Finally, the overexpression of miR-145-5p was found to rejuvenate O-ASCs phenotype and augment the functionality of these cells.
Our research may provide insights regarding the underlying mechanisms of ASC dysfunction; it may also offer novel targets for restoring therapeutic properties in ASCs.
Stem cell-based therapies provide an alternative option for tissue regeneration and have been widely studied for the treatment of various human diseases in the 21st century. The most recent clinical successes involving limbal stem cells in corneal restoration and transgenic stem cells for epidermis regeneration have been compelling [1, 2], which presumably reflect the considerable potential of this field. However, most of these studies remain in preclinical or clinical trial stages. Although substantial investments have been made, clinical applications remain limited because of the poor efficacy and safety characteristics . Notably, there are certain obstacles we must concern when applying stem cells to treat diseases. First, there is a risk of allosensitization when transplanting allogenic stem cells . Second, intravascular stem cell therapy may trigger instant blood-mediated inflammatory reaction, as well as thrombosis and embolization . Third, the hostile host microenvironment hinders stem cell vitality and functionality . Finally, the therapeutic effect of stem cells can be easily impaired by physical conditions within donors (e.g., high age, diabetes, obesity, and chronic kidney disease) [7,8,9]. Given the high morbidity of diabetes and the rapidly aging global population, the donor condition merits extensive consideration.
Stem cells dysfunction causes localized failed tissue maintenance and regeneration, while limiting overall effectiveness in repairing distant tissues and organs [8, 9]. Adipose-derived stem cells (ASCs) have been recognized as a promising tool for the treatment of many disorders (e.g., nerve injury, skin wound healing, cardiac diseases, and autoimmune disorders) [10,11,12]. ASC functions are primarily mediated by their paracrine factors such as growth factors, cytokines, and extracellular vesicles [13, 14]. However, previous studies demonstrated that diabetic or aged ASCs had reduced capacities to promote neovascularization and skin wound healing [15, 16]. Other studies showed that diabetes and aging led to the upregulation of inflammatory markers and changes in ASC immunomodulatory properties [17, 18]. Furthermore, ASCs isolated from donors with diabetes showed impaired effects on the treatment of critical limb ischemia and diabetes-related complications [19, 20]. These cells showed inferior proliferative and proangiogenic capacities, as well as increased intracellular reactive oxygen species (ROS) accumulation; they also exhibited a pro-inflammatory secretome and senescence-associated phenotype, compared with ASCs from normal donors [21,22,23]. Thus, it is important to consider the negative influences of diabetes and aging on ASCs before implantation. Thus far, the precise molecular mechanisms of these alterations and impairments in the niche remain largely unknown. Identification of the underlying mechanisms may provide opportunities to reverse ASC malfunctions and achieve improved stem cell-based therapies.
Regulatory non-coding RNAs—including microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs)—have been demonstrated to play essential roles in many cellular processes and influence various human diseases . Recent studies have shown that these RNAs extensively participate in the proliferation, migration, differentiation, quiescence, and senescence of stem cells [25, 26]. miRNAs are small non-coding RNAs (22-26 nucleotides) that can repress gene expression post-transcriptionally . Emerging evidences elucidated that diabetes and aging manipulated the expression levels of some miRNAs in stem cells, which might impair their functionalities [22, 28]. Moreover, modifying the levels of miR-122 , miR-5591-5p , and miR-34a-5p  could enhance the therapeutic efficacies of ASCs against liver fibrosis, diabetic wound healing, and ischemic myocardial infarction, respectively. lncRNAs are a cluster of transcripts with lengths greater than 200 nucleotides, which can function as signals, decoys, guides, and scaffolds . Li et al.  found that the overexpression of lncRNA Bmncr could reverse aging-related damage in bone marrow mesenchymal stem cells (BMMSCs) and promote their ability to engage in bone formation. circRNAs are covalently closed, endogenous transcripts that function by regulating transcription and protein production, by interfering with splicing or by translating themselves . Recent studies have demonstrated the important roles of circRNAs in aging, as well as various cancers and diabetes. Moreover, circRNAs and lncRNAs can function as competing endogenous RNAs (ceRNAs), which compete with miRNAs for binding to mRNAs, thus promoting the translation of target mRNAs . ceRNA networks have been implicated in stem cell function . Thus, the mechanism by which non-coding RNAs affect the therapeutic capacities of ASCs in relation to diabetes and aging warrants comprehensive investigation.
In this study, we conducted high-throughput RNA sequencing (RNA-seq) in ASCs isolated from type 2 diabetic, old and young patients to systematically explore dysregulated miRNAs, mRNAs, lncRNAs and circRNAs. To the best of our knowledge, this is the first establishment of lncRNA- and circRNA-associated ceRNA networks for simultaneous interpretation of poor ASC function in the context of diabetes and aging. Furthermore, we demonstrated that the lncRNA RAET1E-AS1–miR-145-5p–WNT11/BMPER axis may play a vital role in regulating ASC senescence and function. Therefore, we presume that this study will provide a theoretical basis and novel methods to restore the function of ASCs, thus promoting future ASC-based therapies.
Methods and materials
ASCs were isolated and incubated using our previously described protocols . ASCs were cultured in Dulbecco’s modified Eagle’s medium (DMEM, Gibco, USA) supplemented with 10% fetal bovine serum (FBS, Serapro, USA) and 1% penicillin/streptomycin. ASCs were identified and characterized using a flow cytometric analysis and multi-lineage differentiation method, in accordance with previous published articles [11, 13, 37]. The patients under skin flap operations were selected from hand surgery department of Wuhan union hospital at the Tongji Medical College of Huazhong University of Science and Technology. These donors were consented at first, and then assigned into three groups: young group (Y-ASCs), old group (O-ASCs), and diabetic group (D-ASCs). The detailed information of these groups was described in Table 1. Human umbilical vein endothelial cells (HUVECs) (#GDC166, CCTCC) were obtained from the China Center for Type Culture Collection (CCTCC, Wuhan, China). Human foreskin fibroblasts were isolated and cultured following our previous protocols .
ASC proliferation and migration
ASCs at passage 4 from the above groups were grown in 96-well plates (5000 cells per well) for 24 h. After 2 h incubation with the EdU, the proliferation rates of ASCs from each group were evaluated with Cell-Light EdU Apollo In Vitro Kit (Ribobio, Guangzhou, China).
The 24-well Transwell Chamber (8.0 μm pore size, Corning, USA) was used for assessing the migration ability of ASCs from above groups. Briefly, ASCs (20,000 cells/well) suspended in DMEM without serum were added to the upper compartment, and then incubated in complete culture medium containing 10% FBS for 24 h. Then, ASCs migrated to the bottom surface were stained with crystal violet staining (Solarbio, Beijing, China) and counted under microscopy.
Senescence-associated β-galactosidase staining
Senescence-associated β-galactosidase staining was performed in accordance with a previously published method . ASCs (10,000 cells/well) at passage 4 from above groups were seeded in 48-well plates. After incubation for 3 days, the activity of senescence-associated β-Galactosidase was detected using Senescence β-gal Staining kit (#9860, CST) according to the manufacturer’s instruction.
Wound scratch and tube formation
ASCs from each group were grown in six-well plates (2 × 105 cells per well) for 2 days. The complete culture medium was then replaced with serum-free DMEM. After another day incubation, the medium was collected and considered as the conditioned medium.
Human fibroblasts (5 × 104 cells per well) were seeded in 12-well plates and cultured until 100% confluence. Then, the monolayer was scored with a 200-μl sterile pipette tip. After that, the culture medium was replaced with conditioned media from above groups and maintained for 24 h. The cells were then stained with crystal violet and counted.
HUVECs (2 × 104 cells per well) suspended in conditioned media were seeded in 96-well plates previously coated with 50 μl Matrigel Basement Membrane Matrix (BD Biosciences, USA). After 4 h incubation, tube formation was detected with the microscopy.
ASCs at passage 3 were used for RNA extraction and sequencing. Standard cDNA libraries were conducted and sequenced using the BGISEQ-500 platform and Illumina HiSeq 4000 platform (BGI-Shenzhen, China). The sequence data were filtered to get clean data, and then mapped to the reference genome with Bowtie2  (version 2.2.5). Further, differential expression analysis was performed using the DEGseq . The expected fragments per kilobase of transcript per million fragments sequenced (FPKM)  was used to determine the expression level of mRNAs, lncRNAs and circRNAs. The expression level of miRNAs was determined by the transcript per million (TPM) . The significantly dysregulated RNAs must meet the following criteria: (1) fold changes ≥ 2 or ≤ − 2, (2) p and q value < 0.001, and (3) TPM value ≥ 10 for miRNA or FPKM value ≥ 0.1 for other RNAs. Detailed RNA-seq data have been deposited at Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo; GSE174502).
RT-PCR was performed in accordance with a previously published method . Total RNA was isolated with miRNeasy Mini Kit (Qiagen, Germany) according to the manufacturer’s instructions. MiRNAs and other RNAs were reversely transcribed into cDNA using the Mir-X™ miRNA First-Strand Synthesis kit (#RR638315, TaKaRa) and Prime-Script® RT reagent Kit (#RR047A, TaKaRa) separately. Real-time PCR was performed on the StepOnePlus™ platform (Applied Biosystems, USA) using TB Green® Premix Ex Taq™ II kit (#RR820A, TaKaRa). The primer sequences were presented in Tables S2, S3, S4, S5. The relative expression levels of targeted genes were calculated using the 2−ΔΔCt method and normalized to β-actin, gapdh, or u6.
Western blotting was performed in accordance with a previously published method . Total proteins were extracted by ripa lysis buffer with proteinase inhibitor (Roche, Switzerland). In total, 30 μg aliquots of protein were separated in 10% SDS-PAGE gels, then transferred onto PVDF membranes (Millipore, USA). After 1.5 h of blocking with 5% w/v nonfat dry milk buffer, the membrane was incubated overnight with primary antibodies for fibronectin (1:1000, #A12932, ABclonal, China), cyclin D1 (1:1000, #A11022, ABclonal, China), cyclin A1 (1:1000, #ab206746, Abcam), NANOG (1:1000, #14295-1-AP, Proteintech, China), OCT4 (1:1000, #11263-1-AP, Proteintech, China), p21 (1:1000, #10355-1-AP, Proteintech, China), and GAPDH (1:1000, #10494-1-AP, Proteintech, China). Then, the membrane was incubated with secondary antibodies (1:5000, Aspen, China) and exposed to X-ray film (UVP, USA).
ceRNA and PPI network analysis
The significant differentially expressed mRNAs and ncRNAs between each group were used for ceRNA network constructions. The RNAs which could be predicted by at least two of these databases (RNAhybrid [https://bibiserv.cebitec.uni-bielefeld.de/rnahybrid], miRanda [http://www.microrna.org/microrna/home.do] and TargetScan [http://www.targetscan.org]) were considered as the miRNA targets. The sequences of mRNAs, lncRNAs, and circRNAs were screened to get the potential MREs. Protein-protein interaction (PPI) analysis of differentially expressed mRNAs was based on the STRING database (https://string-db.org). These networks were illustrated using Cytoscape 3.7.1. The degree centrality of involved genes was calculated by Cytoscan.
GO and KEGG enrichment analyses
The differentially expressed mRNAs were analyzed by Gene Ontology (GO) database and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database. The GO categories (http://geneontology.org) were used for defining the molecular function and biological process involved of candidate genes. The biological function of these genes was further annotated by KEGG database (http://www.genome.jp/kegg). We used the hypergeometric distribution test to find out significantly enriched gene sets. A p value < 0.05 was considered statistically significant.
The mimic and inhibitor of hsa-miR-145-5p were designed and synthesized by Ribobio (Guangzhou, China). ASCs were transfected using riboFECT™CP Reagent according to the manufacturer’s instruction. After 48 h of incubation, ASCs were harvested for subsequent experiments. The sequences of these oligonucleotides are described in Table S5.
All statistical analysis was conducted using Graph Pad prism v 7.0 software. All data were expressed as the mean ± SEM. Unpaired Student’s t test was used for evaluating the statistical discrepancy between two groups. For group ≥ 3, one-way ANOVA with Bonferroni post hoc test was employed. The Pearson correlation coefficient was used to calculate the expressional correlation of two RNAs. Statistical significance was set at p or q < 0.05.
ASCs from diabetic or old donors showed inferior proliferation and migration ability, exhibited a senescent phenotype
To identify whether diabetes and aging affected ASC phenotype, a series of fundamental experiments was conducted. First, the EdU assay showed that the proliferation rate of ASCs from old donors was 24.48%, which was significantly lower than the proliferation rate of ASCs from young donors (38.88%) (Fig. 1a, b). Notably, diabetic conditions did not reduce the proliferation rate of ASCs (Fig. 1a, b). Next, Transwell assays revealed that the migration ability of ASCs was best in the young group, followed by the old group; the worst ability was observed in the diabetic group (Fig. 1c, d). Last, β-galactosidase staining showed that the rate of senescent stem cells in diabetic group was significantly higher in the diabetic group than in the old group, while ASCs in the young group exhibited the lowest number of senescent cells (Fig. 1e, f). Overall, diabetes and aging both have detrimental effects on ASC phenotype.
Diabetic conditions and aging impaired the capacities of ASCs for modulating fibroblast and HUVEC functions
Given the above results, we investigated whether diabetic conditions and aging impaired the therapeutic abilities of ASCs. Thus, we evaluated the differences in ASCs from these three groups in terms of their abilities to optimize HUVEC and fibroblast functions. In vitro tube formation assays showed that HUVECs treated with conditioned media from Y-ASCs exhibited twofold increases in the number of closed tubular structures, compared with HUVECs treated with conditioned media from O-ASCs (Fig. 1g, h). Moreover, the number of tubular structures was twofold greater in the O-ASC group than in the D-ASC group, but this difference was not statistically significant (Fig. 1g, h). Notably, wound scratch assays revealed that the migration of fibroblasts treated with conditioned media from Y-ASCs was significantly greater than the migration of fibroblasts treated with conditioned media from O-ASCs (Fig. 1i, j). Additionally, the number of migratory fibroblasts was slightly lower in the D-ASC group than in the O-ASC group (Fig. 1i, j). Therefore, old donor age has substantially deleterious effects on the abilities of ASCs to modulate cellular functions, while the presence of diabetes in donors may have a slight harmful effect on ASC function.
Differential expression analyses: Y-ASCs vs. O-ASCs and D-ASCs vs. O-ASCs
Following identification of the negative influences of diabetes and aging on ASC phenotype and function, we performed RNA-seq analysis on ASCs from these three groups to uncover the underlying mechanisms. The differentially expressed miRNAs, mRNAs, lncRNAs, and circRNAs from all three groups (n = 3) were illustrated using hierarchical clustering heat maps (Fig. 2). In total, 34 (27 up and seven down) differentially expressed miRNAs, 553 (291 up and 262 down) differentially expressed mRNAs, 132 (39 up and 92 down) differentially expressed lncRNAs, and 196 (123 up and 73 down) differentially expressed circRNAs were detected between Y-ASCs and O-ASCs (Fig. 2a–d, Data S1, S2, S3, S4). In this context, up indicates RNAs upregulated in Y-ASCs, compared with O-ASCs. Comparing D-ASCs with O-ASCs, we found 43 (24 up and 19 down) differentially expressed miRNAs, 926 (483 up and 443 down) differentially expressed mRNAs, 535 (496 up and 39 down) differentially expressed lncRNAs, and 1966 (1806 up and 160 down) differentially expressed circRNAs (Fig. 2e–h, Data S5, S6, S7, S8). In this context, up indicates RNAs upregulated in D-ASCs, compared with O-ASCs. Overall, we found that four miRNAs, 22 mRNAs, eight circRNAs were upregulated in O-ASCs, compared with D-ASCs; these RNAs were all upregulated in Y-ASCs, compared with O-ASCs (Data S9, S10, S11). Moreover, 15 mRNAs, 14 lncRNAs, and 17 circRNAs were downregulated in O-ASCs, compared with D-ASCs; these RNAs were all downregulated in Y-ASCs, compared with O-ASCs (Data S10, S11, S12).
Functional enrichment analysis of differentially expressed mRNAs
GO and KEGG analyses were conducted to investigate the biological functions of dysregulated mRNAs. Comparing Y-ASCs with O-ASCs, GO biological process analysis showed that 750 GO terms were significantly enriched (Data S13). Terms such as aging, positive regulation of cell proliferation, positive regulation of cell migration and angiogenesis were potentially associated with donor age-related ASC dysfunction (Fig. 3a). KEGG pathway analysis showed that 47 pathways were statistically enriched (Data S14). Among them, the cellular senescence, TGF-beta signaling, p53 signaling, and PPAR signaling pathways may be involved in impaired ASC functions (Fig. 3b). Comparing O-ASCs with D-ASCs, 553 GO terms were significantly enriched in GO biological process analysis (Data S15). Terms such as oxidation-reduction process, aging, cell migration, wound healing and angiogenesis may be associated with decreased efficacy in ASCs derived from patients with diabetes (Fig. 3c). KEGG pathway analysis revealed that 37 pathways were statistically enriched (Data S16). Among them, the AGE-RAGE signaling, TGF-beta signaling, p53 signaling, mTOR signaling, and cellular senescence pathways may be involved in the impaired functions of ASCs isolated from diabetic donors (Fig. 3d).
PPI network construction
PPI networks were constructed to identify critical genes among the differentially expressed mRNAs. Comparing Y-ASCs with O-ASCs, the established network comprised 107 nodes and 373 edges (Fig. 3e). In this network, the top 15 genes with the highest core degree were APOE, BMP4, EGR1, GPC3, CXCL8, SPP1, FOS, HGF, TNC, IGFBP7, JUN, APOL1, GBP2, HLA-F, and OAS1 (Fig. 3e). Comparing O-ASCs with D-ASCs, the established network comprised 131 nodes and 783 edges (Fig. 3f). In this network, the top 15 genes with the highest core degree were C3, CXCL8, IL6, CXCL1, CCNB1, BUB1, CXCL12, MCHR2, CXCL2, CXCL6, CXCR4, CXCL5, QSOX1, S1PR1, and TAS2R31 (Fig. 3f).
RT-PCR confirmation of differentially expressed miRNAs
We randomly selected 14 miRNAs to validate the reliability of RNA-seq data using RT-PCR. Comparing Y-ASCs with O-ASCs, PCR analysis showed that miR-145-3P, miR-145-5p, miR-126-3p, miR-126-5p, miR-214-3p, miR-181a-3p, and miR-210-3p were upregulated in Y-ASCs; these findings were consistent with the RNA-seq results (Fig. 4a, b). Furthermore, PCR analysis showed miR-766-3p was downregulated in Y-ASCs, which conflicted with the RNA-seq findings (Fig. 4a, b). Comparing D-ASCs with O-ASCs, PCR analysis showed that miR-214-3p, miR-193a-3p, and miR-145-3p were downregulated in D-ASCs, miR-3529-3p, miR-302a-3p, and miR-302b-3p were upregulated in D-ASCs; these findings were consistent with the RNA-seq results (Fig. 4c, d). However, PCR analysis of two miRNAs (miR-615-3p and miR-12136) showed findings that contrasted with the RNA-seq results (Fig. 4c, d). The remaining six miRNAs showed no significant difference between groups according to PCR analysis (Fig. 4e, f).
Construction of circRNA– or lncRNA–miRNA–mRNA network
ceRNA networks were constructed based on the differentially expressed RNAs. miRNAs verified by PCR analysis, which exhibited expression patterns consistent with RNA-seq results, were selected as the cores of ceRNA networks. Comparing Y-ASCs with O-ASCs, the lncRNA-associated ceRNA network comprised four miRNAs, 171 mRNAs, and 46 lncRNAs (Fig. 5a); the circRNA-associated ceRNA network comprised three miRNAs, 160 mRNAs, and 36 circRNAs (Fig. 5b). Comparing D-ASCs with O-ASCs, the lncRNA-associated ceRNA network comprised four miRNAs, 118 mRNAs, and 97 lncRNAs (Fig. 6a); the circRNA-associated ceRNA network comprised four miRNAs, 119 mRNAs, and 103 circRNAs (Fig. 6b).
RT-PCR confirmation of differentially expressed mRNAs, lncRNAs, and circRNAs
mRNAs, lncRNAs, and circRNAs contained in the above PPI and ceRNA networks were selected for further RT-PCR analyses. Comparing Y-ASCs with O-ASCs, RT-PCR results showed that 14 mRNAs were significantly dysregulated (Fig. 7a). Among these dysregulated mRNAs, MMP3, EGR1, JUNB, and BMP4 were involved in the PPI network (Figs. 3e and 7a); ANK2, PODN, NOVA1, CLEC3B, RARRES3, IGFBP6, and MFAP5 were involved in the ceRNA network (Figs. 5a, b and 7a); and ITGB8, WNT11, and BMPER were involved in both PPI and ceRNA networks (Figs. 3e; 5a, b; and 7a). Comparing D-ASCs with O-ASCs, 11 mRNAs were confirmed to be differentially expressed by PCR analyses (Fig. 7b); of these, CXCL8, CXCL1, CXCL6, and CXCL5 were all upregulated in O-ASCs and involved in the PPI network (Figs. 3f and 7b). Furthermore, SOX4, ANGPT1, SLC5A3, NDUFV3, and LUC7L3 were involved in the ceRNA network (Figs. 6a, b and 7b); COL18A1 and FMOD were involved in both PPI and ceRNA networks (Figs. 3f; 6a, b; and 7b). Moreover, we found that five lncRNAs (RAET1E-AS1, LOC102723591, LOC105373230, LOC112267883, and LINC02595) and five circRNAs (hsa_circ_0017534, hsa_circ_0092630, hsa_circ_0080906, hsa_circ_0075045, and hsa_circ_0080909) were significantly upregulated in O-ASCs, compared with Y-ASCs (Fig. 7c). Additionally, PCR analysis revealed that two lncRNAs (LOC102724087 and LOC101928000) and three circRNAs (hsa_circ_0088199, hsa_circ_0088195, and hsa_circ_0058158) were significantly upregulated in O-ASCs, compared with D-ASCs (Fig. 7d). Conversely, two lncRNAs (LOC105377989 and NEAT1) were significantly downregulated in O-ASCs, compared with D-ASCs (Fig. 7d). Other RNAs without statistically significant differences between these groups are displayed in Figure S1A and B.
Establishment of ceRNA sub-networks based on RNAs verified by PCR analyses
Based on the above PCR results, we constructed circRNA– and lncRNA–miRNA–mRNA sub-networks for further investigation. Comparing Y-ASCs with O-ASCs, the network included five miRNAs, 10 mRNAs, five lncRNAs, and five circRNAs (Fig. 7e). Comparing D-ASCs with O-ASCs, the network included six miRNAs, seven mRNAs, four lncRNAs, and three circRNAs (Fig. 7f). We randomly selected hsa-miR-145-5p for further analysis; it had ceRNA relationships with four mRNAs (ITGB8, WNT11, BMPER, and RARRES3), two lncRNAs (RAET1E-AS1 and LOC102723591), and one circRNA (hsa_circ_0092630) (Fig. 7e). Correlation analyses revealed that the expression levels of RAET1E-AS1, WNT11, and BMPER had significantly negative correlations with the expression level of hsa-miR-145-5p (Fig. 7g–i), while the other RNAs did not exhibit such relationships (Figure S1C, D). Furthermore, the expression level of RAET1E-AS1 had a significantly positive correlation with the expression levels of WNT11 and BMPER (Fig. 7j, k). Finally, we found that overexpression of has-miR-145-5p in ASCs could significantly downregulated the expression levels of WNT11, BMPER, RAET1E-AS1, LOC102723591, and hsa_circ_0092630 (Fig. 7l). Therefore, we speculated that the RAET1E-AS1–hsa-miR-145-5p–WNT11/BMPER network might play a vital role in the phenotype and function of ASCs.
Overexpression of hsa-miR-145-5p improved the phenotype and function of ASCs, while its inhibition adversely affected the cells.
To evaluate the characteristics of the RAET1E-AS1–hsa-miR-145-5–WNT11/BMPER axis, we modulated the expression level of hsa-miR-145-5p in ASCs from old donors to investigate changes in their phenotypes and functions. miR-145-5p mimic, mimic-NC, miR145-5p inhibitor, or inhibitor NC was transfected into ASCs separately, and transfection efficiency was evaluated by PCR analysis (Figure S2). EdU and Transwell assays showed that the proliferation and migration abilities of ASCs were promoted by more than twofold after hsa-miR-145-5p overexpression, while its inhibition significantly reduced these abilities (Fig. 8a–d). Furthermore, the number of β-galactosidase-stained cells was twofold less in the miR-145-5p group than in the mimic-NC group, while inhibition of miR-145-5p significantly increased the number of stained cells (Fig. 8e, f). Additionally, in vitro tube formation and wound scratch assays revealed that cellular supernatants of cultured ASCs with miR-145-5p overexpression could remarkably enhance the functions of HUVECs and fibroblasts, as shown by increases in closed tubular structures and migratory cells; the inhibition of miR-145-5p led to opposite results (Fig. 8g–j). Additionally, western blot analysis showed that overexpression of miR-145-5p in ASCs could significantly improve the expression of migration–associated protein FN1, proliferation–associated proteins CCNA1 and CCND1, and pluripotent markers NANOG and OCT4; it could decrease the expression of senescence–associated protein p21 (Fig. 8k). PCR analysis showed that overexpression of miR-145-5p in ASCs could significantly improve the expression of MMP9, NANOG, and growth factors PDGFA and FGF; it could decrease the expression of IGF1, HIF-1, CXCL8, and senescent markers TP53 and CDKN2A (Fig. 8l). Overall, the overexpression of miR-145-5p in ASCs could ameliorate the unsatisfactory phenotype and function of ASCs isolated from old donors.
Stem cells have an important role in the field of translational medicine, largely through the secretion of trophic factors rather than direct involvement in tissue reconstruction . Clinical trials often use autologous stem cells to treat human diseases without consideration of immune rejection . Unfortunately, patients who may benefit from stem cell therapies often exhibited accompanying systemic diseases, which may compromise the functionalities of their own stem cells and cause uncertainty of autologous transplantation . ASCs have been recognized as a promising stem cell source in recent years because of a series of specific advantages: convenient acquisition of adipose tissue, lack of ethical concerns, and effortless in vitro isolation and expansion . However, poor donor physical conditions (e.g., older age and diabetes) can impair the overall benefits of ASCs and limit their applications [22, 23]. In the present study, we hypothesized that diabetic conditions and aging might change the epigenome, transcriptome, proteome, and metabolome of ASCs. We mainly focused on the transcriptome and conducted RNA-seq to uncover the corresponding underlying mechanisms.
First, we assessed the functional and phenotypic differences of ASCs isolated from diabetic, old, and young donors. The results were generally similar to the findings of previous studies. We found that ASCs isolated from old donors exhibited impaired proliferation and migration abilities; moreover, they easily became senescent after in vitro expansion, as described previously . Previous studies showed that O-ASCs had weakened capacities for accelerating skin wound healing compared to Y-ASCs; similarly, we found that O-ASCs had inferior capacities for modulating HUVEC and fibroblast functions . However, we found that diabetes did not influence the proliferation of ASCs, which considerably differed from the findings in most previous studies [20, 22]. This discrepancy might be related to differences in original sources, and expansion methods, as well as other factors. Importantly, the migratory and senescent phenotypes of diabetic ASCs in our study were consistent with the findings in previous research . Diabetic conditions reportedly impair the therapeutic effects of ASCs in wound healing and critical limb ischemia [19, 44]. Our research also revealed that diabetic ASCs were slightly impaired in their capacities to promote HUVEC angiogenesis and fibroblast migration, which might explain the results of the prior studies. Thus, donor physical conditions must be considered before stem cell therapies.
To the best of our knowledge, the pathogenesis and progression of many diseases (e.g., cancers and metabolic diseases) are often accompanied by the transcriptome alterations. Non-coding RNAs, which are vital components of the transcriptome, have been shown to participate in the pathophysiologic mechanisms underlying those diseases [24, 45]. Previous microarray analyses have identified differences in the expression levels of miRNAs and lncRNAs in BMMSCs isolated from young and old donors [46, 47]. However, the overall ASC transcriptome alterations related to diabetes and aging have not been fully characterized. In the present study, we used RNA-seq to simultaneously explore expression differences in miRNAs, mRNAs, lncRNAs, and circRNAs among ASCs isolated from diabetic, old, and young donors. Our study found that large amounts of RNAs were differentially expressed in Y-ASCs vs. O-ASCs and in D-ASCs vs. O-ASCs. GO analyses showed biological processes such as aging, angiogenesis, cell adhesion, positive regulation of cell proliferation, and cell-cell signaling were both enriched in the two comparisons. Furthermore, the telomere organization term was significantly enriched in GO analyses that compared D-ASCs with O-ASCs. Telomere shortening is recognized as a hallmark of stem cell senescence , and a previous study showed that diabetes could significantly affect telomere length . KEGG analyses revealed that the AGE-RAGE signaling, cellular senescence, and p53 signaling pathways were both enriched in the two comparisons. Notably, the p53 signaling pathway reportedly involved in the proliferation, senescence, apoptosis, and differentiation of stem cells .
Moreover, we conducted PPI network analyses to further characterize the differentially expressed mRNAs. Comparing Y-ASCs with O-ASCs, we found that MMP3, EGR1, JUNB, BMP4, ITGB8, WNT11, and BMPER were involved in the network; these findings were verified by PCR analyses. A recent study revealed that the induction of EGR1 could trigger glioblastoma cell dedifferentiation into a stem-like state, which involved the expression of pluripotent markers NANOG and OCT4 . JUNB, a member of the AP-1 transcription factor family, has been shown to regulate epidermal stem cells by balancing progenitor proliferation and differentiation . BMP4, a member of the BMP–SMAD signaling axis, could suppress p16/INK4A-mediated cell senescence . In our study, EGR1, JUNB, and BMP4 were all downregulated in O-ASCs, compared with Y-ASCs; this might explain the dysfunction observed in O-ASCs. Additionally, PCR analyses showed that the network members CXCL8, CXCL1, CXCL6, CXCL5, COL18A1, and FMOD were all downregulated in D-ASCs, compared with O-ASCs. These results contrasted markedly from the findings of a previous study, in which ASCs from patients with diabetes or atherosclerosis strongly expressed pro-inflammatory markers such as IL-1β, IL6, and IL-8/CXCL8 [17, 18]. However, a recent study demonstrated that the overexpression of dual chemokines CXCL6 and SDF-1α in ASCs could promote angio-vasculogenesis in those cells . Because these chemokines exhibit multiple functions , their specific roles in ASCs require further investigation.
ceRNA networks have been described as an intricate interplay among diverse RNA species, in which mRNAs, lncRNAs, pseudogenes, and circRNAs compete for binding to shared miRNAs . This crosstalk has been shown to participate in many biological processes and human diseases, especially involving oncogenesis and disease progression . However, it has received minimal attention in recent research. A previous study showed that LINC00707 could promote osteogenesis in human BMMSCs by acting as a ceRNA to upregulate WNT2B via miR-370-3p inhibition . circFOXP1 could sustain mesenchymal stem cell identity and differentiation by competing with WNT5A for miR-17-3p/miR-127-5p binding . Because miRNAs are the cores of the ceRNA network, we first evaluated the miRNAs that were differentially expressed in RNA-seq analysis by a secondary PCR analysis approach. Our results revealed that eight miRNAs were differentially expressed between Y-ASCs and O-ASCs. Of these eight miRNAs, miR-210-3p  and miR-126  could inhibit apoptosis in ASCs and BMMSCs, while miR-214-3p could anticipate in the therapeutic application of mesenchymal stem cell-derived extracellular vesicles for attenuating radiation-induced lung injury . However, all three miRNAs were downregulated in O-ASCs, which might have contributed to their poor function. Additionally, we found that eight miRNAs were differentially expressed between D-ASCs and O-ASCs. Among these eight miRNAs, miR-302 clusters were strongly expressed in human embryonic stem cells and associated with stem cell pluripotency . Surprisingly, most of the above miRNAs have not been extensively investigated in terms of mesenchymal stem cell phenotypes and functions. Thus, further studies are needed to investigate their characteristics.
Next, we constructed ceRNA networks using bioinformatics methods based on the differentially expressed mRNAs, lncRNAs, circRNAs, and verified miRNAs. The four network maps illustrated extensive potential ceRNA relationships that might play vital roles in ASC biology. Subsets of RNAs in the networks were analyzed by PCR to establish sub-networks. Most involved genes comprised extracellular proteins, and some genes (e.g., COL18A1, MFAP5, ITGB8, FMOD, and PODN) were extracellular matrix-related proteins. SOX4 and ANGPT1, which participate in stemness and angiogenesis, were downregulated in D-ASCs in our study. Liu et al.  reported that BMMSCs transfected with ANGPT1 had therapeutic effects on hyperoxia induced optic nerve injury. With respect to non-coding RNAs in the networks, NEAT1 has been proven to act as an oncogene in multiple cancers by binding to various miRNAs [63, 64]. LINC02595 could promote colorectal cancer progression by inhibiting miR-203b-3p activity . Notably, the biological functions of other lncRNAs and circRNAs have not been reported thus far.
The RAET1E-AS1–miR-145-5p–WNT11/BMPER axis was conformed on the basis of ceRNA relationships and expression correlation analyses. miR-145-5p has been recognized as a tumor repressor in many types of cancers [66, 67], a suppressor in the osteogenic differentiation of ASCs , and a biomarker for the diagnosis of type 2 diabetes because of its low level in plasma from patients with diabetes . Our study found that the overexpression of miR-145-5p in O-ASCs could rejuvenate cellular senescence and enhance their function. Thus, the modulation of miR-145-5p in ASCs might be a promising method for promoting stem cell function; further investigations are needed to explore the therapeutic effects of modified ASCs in animal models. Chen et al.  found that WNT11 overexpression could inhibit ASC proliferation and induce their differentiation into nucleus pulposus cells and osteoblasts. Pérez et al.  found that BMPER was upregulated in diet- and obese-derived ASCs; moreover, it could decrease the migration abilities of those cells. However, the role of the lncRNA RAET1E-AS1 in cell biology has not been explored. Further fundamental studies are needed to investigate the specific function of the RAET1E-AS1–WNT11–BMPER network in ASCs and its direct interactions with miR-145-5p.
In the present study, we simultaneously examined changes in the expression patterns of miRNAs, mRNAs, lncRNAs, and circRNAs in ASCs from older donors, younger donors, and donors with diabetes. GO and KEGG pathway analyses were conducted to identify the possible functions of differentially expressed mRNAs. PPI networks were established to find protein genes with critical roles in our disease models. ceRNA networks that included lncRNA–miRNA–mRNA and cirRNA–miRNA–mRNA interactions were successfully constructed based on the bioinformatics analyses and PCR results. Overexpression of miR-145-5p could rejuvenate the O-ASC phenotype and augment their abilities to modulate endothelial cell and fibroblast functions (Fig. 9). Thus, this study may contribute to the broader understanding of the underlying mechanisms of ASC instability and provide novel targets to reverse dysfunction in ASCs isolated from diabetic and old patients.
Availability of data and materials
Advanced glycation end products
Adipose-derived stem cells
Bone marrow mesenchymal stem cells
Bone morphogenetic protein 4
BMP binding endothelial regulator
Cyclin-dependent kinase inhibitor 2A
Competing endogenous RNA
C-type lectin domain family 3 member B
Collagen type XVIII alpha 1 chain
C-X-C motif chemokine ligand 1
C-X-C motif ligand 12
C-X-C motif chemokine ligand 2
C-X-C motif chemokine ligand 5
C-X-C motif chemokine ligand 6
C-X-C motif ligand 8
C-X-C motif chemokine receptor 4
Early growth response protein 1
Fibroblast growth factor
Fragments per kilobase of transcript per million fragments sequenced
Guanylate binding protein 2
Hepatocyte growth factor
Hypoxia inducible factor 1 subunit alpha
Major histocompatibility complex, class I, F
Human umbilical vein endothelial cells
Instant blood-mediated inflammatory reaction
Insulin-like growth factor 1
Insulin-like growth factor binding protein 6
Insulin-like growth factor binding protein 7
Integrin subunit beta 8
Kyoto Encyclopedia of Genes and Genomes
Long non-coding RNA
LUC7 like 3 pre-mRNA splicing factor
Melanin concentrating hormone receptor 2
Microfibril-associated protein 5
Matrix metallopeptidase 9
NADH: ubiquinone oxidoreductase subunit V3
Nuclear paraspeckle assembly transcript 1
NOVA alternative splicing regulator 1
2′-5′-oligoadenylate synthetase 1
Organic cation/carnitine transporter4
Platelet-derived growth factor subunit A
Peroxisome proliferators-activated receptor
Quiescin sulfhydryl oxidase 1
RAET1E antisense RNA 1
Receptor for advanced glycation end products
Retinoic acid receptor responder 3
Reactive oxygen species
Real-time polymerase chain reaction
Sphingosine-1-phosphate receptor 1
Solute carrier family 5 member 3
SRY-box transcription factor 4
Secreted phosphoprotein 1
Taste 2 receptor member 31
Blau HM, Daley GQ. Stem cells in the treatment of disease. N Engl J Med. 2019;380(18):1748–60. https://doi.org/10.1056/NEJMra1716145.
Cossu G, Birchall M, Brown T, De Coppi P, Culme-Seymour E, Gibbon S, et al. Lancet Commission: Stem cells and regenerative medicine. Lancet. 2018;391(10123):883–910. https://doi.org/10.1016/S0140-6736(17)31366-1.
Trounson A, McDonald C. Stem cell therapies in clinical trials: progress and challenges. Cell Stem Cell. 2015;17(1):11–22. https://doi.org/10.1016/j.stem.2015.06.007.
Welsch CA, Rust WL, Csete M. Concise review: lessons learned from islet transplant clinical trials in developing stem cell therapies for type 1 diabetes. Stem Cells Transl Med. 2019;8(3):209–14. https://doi.org/10.1002/sctm.18-0156.
Moll G, Ankrum JA, Kamhieh-Milz J, Bieback K, Ringdén O, Volk H, et al. Intravascular mesenchymal stromal/stem cell therapy product diversification: time for new clinical guidelines. Trends Mol Med. 2019;25(2):149–63. https://doi.org/10.1016/j.molmed.2018.12.006.
Mazini L, Rochette L, Admou B, Amal S, Malka G. Hopes and limits of adipose-derived stem cells (ADSCs) and mesenchymal stem cells (MSCs) in wound healing. Int J Mol Sci. 2020;21(4):1306. https://doi.org/10.3390/ijms21041306.
Hickson LJ, Eirin A, Lerman LO. Challenges and opportunities for stem cell therapy in patients with chronic kidney disease. Kidney Int. 2016;89(4):767–78. https://doi.org/10.1016/j.kint.2015.11.023.
Duscher D, Barrera J, Wong VW, Maan ZN, Whittam AJ, Januszyk M, et al. Stem Cells in wound healing: the future of regenerative medicine? A Mini-Review. Gerontology. 2016;62(2):216–25. https://doi.org/10.1159/000381877.
Louwen F, Ritter A, Kreis NN, Yuan J. Insight into the development of obesity: functional alterations of adipose-derived mesenchymal stem cells. Obes Rev. 2018;19(7):888–904. https://doi.org/10.1111/obr.12679.
Kang Y, Liu Y, Liu Z, Ren S, Xiong H, Chen J, et al. Differentiated human adipose-derived stromal cells exhibit the phenotypic and functional characteristics of mature Schwann cells through a modified approach. Cytotherapy. 2019;21(9):987–1003. https://doi.org/10.1016/j.jcyt.2019.04.061.
Ren S, Chen J, Duscher D, Liu Y, Guo G, Kang Y, et al. Microvesicles from human adipose stem cells promote wound healing by optimizing cellular functions via AKT and ERK signaling pathways. Stem Cell Res Ther. 2019;10(1). https://doi.org/10.1186/s13287-019-1152-x.
Ceccarelli S, Pontecorvi P, Anastasiadou E, Napoli C, Marchese C. Immunomodulatory effect of adipose-derived stem cells: the cutting edge of clinical application. Front Cell Dev Biol. 2020;8. https://doi.org/10.3389/fcell.2020.00236.
Chen J, Ren S, Duscher D, Kang Y, Liu Y, Wang C, et al. Exosomes from human adipose-derived stem cells promote sciatic nerve regeneration via optimizing Schwann cell function. J Cell Physiol. 2019;234(12):23097–110. https://doi.org/10.1002/jcp.28873.
Shukla L, Yuan Y, Shayan R, Greening DW, Karnezis T. Fat therapeutics: the clinical capacity of adipose-derived stem cells and exosomes for human disease and tissue regeneration. Front Pharmacol. 2020;11. https://doi.org/10.3389/fphar.2020.00158.
Rennert RC, Sorkin M, Januszyk M, Duscher D, Kosaraju R, Chung MT, et al. Diabetes impairs the angiogenic potential of adipose-derived stem cells by selectively depleting cellular subpopulations. Stem Cell Res Ther. 2014;5(3):79. https://doi.org/10.1186/scrt468.
Duscher D, Rennert RC, Januszyk M, Anghel E, Maan ZN, Whittam AJ, et al. Aging disrupts cell subpopulation dynamics and diminishes the function of mesenchymal stem cells. Sci Rep. 2014;4:7144.
Serena C, Keiran N, Ceperuelo-Mallafre V, Ejarque M, Fradera R, Roche K, et al. Obesity and type 2 diabetes alters the immune properties of human adipose derived stem cells. Stem Cells. 2016;34(10):2559–73. https://doi.org/10.1002/stem.2429.
Kizilay MO, Lora M, Shum-Tim D, Nadeau S, Rodier F, Colmegna I. A proinflammatory secretome mediates the impaired immunopotency of human mesenchymal stromal cells in elderly patients with atherosclerosis. Stem Cells Transl Med. 2017;6(4):1132–40. https://doi.org/10.1002/sctm.16-0221.
Inoue O, Usui S, Takashima S, Nomura A, Yamaguchi K, Takeda Y, et al. Diabetes impairs the angiogenic capacity of human adipose-derived stem cells by reducing the CD271+ subpopulation in adipose tissue. Biochem Bioph Res Co. 2019;517(2):369–75. https://doi.org/10.1016/j.bbrc.2019.07.081.
Wang M, Song L, Strange C, Dong X, Wang H. Therapeutic effects of adipose stem cells from diabetic mice for the treatment of type 2 diabetes. Mol Ther. 2018;26(8):1921–30. https://doi.org/10.1016/j.ymthe.2018.06.013.
Marędziak M, Marycz K, Tomaszewski KA, Kornicka K, Henry BM. The influence of aging on the regenerative potential of human adipose derived mesenchymal stem cells. Stem Cells Int. 2016;2016:1–15. https://doi.org/10.1155/2016/2152435.
Alicka M, Major P, Wysocki M, Marycz K. Adipose-derived mesenchymal stem cells isolated from patients with type 2 diabetes show reduced “stemness” through an altered secretome profile, impaired anti-oxidative protection, and mitochondrial dynamics deterioration. J Clin Med. 2019;8(6):765. https://doi.org/10.3390/jcm8060765.
Alicka M, Kornicka-Garbowska K, Kucharczyk K, Kępska M, Rӧcken M, Marycz K. Age-dependent impairment of adipose-derived stem cells isolated from horses. Stem Cell Res Ther. 2020;11(1):4. https://doi.org/10.1186/s13287-019-1512-6.
Slack FJ, Chinnaiyan AM. The role of non-coding RNAs in oncology. Cell. 2019;179(5):1033–55. https://doi.org/10.1016/j.cell.2019.10.017.
Shenoy A, Blelloch RH. Regulation of microRNA function in somatic stem cell proliferation and differentiation. Nat Rev Mol Cell Bio. 2014;15(9):565–76. https://doi.org/10.1038/nrm3854.
Hu S, Shan G. LncRNAs in stem cells. Stem Cells Int. 2016;2016:1–8. https://doi.org/10.1155/2016/2681925.
Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136(2):215–33. https://doi.org/10.1016/j.cell.2009.01.002.
Choi SW, Lee JY, Kang K. miRNAs in stem cell aging and age-related disease. Mech Ageing Dev. 2017;168:20–9. https://doi.org/10.1016/j.mad.2017.08.013.
Lou G, Yang Y, Liu F, Ye B, Chen Z, Zheng M, et al. MiR-122 modification enhances the therapeutic efficacy of adipose tissue-derived mesenchymal stem cells against liver fibrosis. J Cell Mol Med. 2017;21(11):2963–73. https://doi.org/10.1111/jcmm.13208.
Li Q, Xia S, Yin Y, Guo Y, Chen F, Jin P. miR-5591-5p regulates the effect of ADSCs in repairing diabetic wound via targeting AGEs/AGER/JNK signaling axis. Cell Death Dis. 2018;9:566.
Weng C, Wu C, Kao S, Chen J, Lin H. Down-regulation of miR-34a-5p potentiates protective effect of adipose-derived mesenchymal Stem cells against ischemic myocardial infarction by stimulating the expression of C1q/tumor necrosis factor-related protein-9. Front Physiol. 2019;10. https://doi.org/10.3389/fphys.2019.01445.
Ponting CP, Oliver PL, Reik W. Evolution and functions of long noncoding RNAs. Cell. 2009;136(4):629–41. https://doi.org/10.1016/j.cell.2009.02.006.
Li C, Xiao Y, Yang M, Su T, Sun X, Guo Q, et al. Long noncoding RNA Bmncr regulates mesenchymal stem cell fate during skeletal aging. J Clin Investig. 2018;128(12):5251–66. https://doi.org/10.1172/JCI99044.
Kristensen LS, Andersen MS, Stagsted LVW, Ebbesen KK, Hansen TB, Kjems J. The biogenesis, biology and characterization of circular RNAs. Nature reviews. Genetics. 2019;20(11):675–91. https://doi.org/10.1038/s41576-019-0158-7.
Tay Y, Rinn J, Pandolfi PP. The multilayered complexity of ceRNA crosstalk and competition. Nature. 2014;505(7483):344–52. https://doi.org/10.1038/nature12986.
Cherubini A, Barilani M, Rossi RL, Jalal MMK, Rusconi F, Buono G, et al. FOXP1 circular RNA sustains mesenchymal stem cell identity via microRNA inhibition. Nucleic Acids Res. 2019;47(10):5325–40. https://doi.org/10.1093/nar/gkz199.
Farahzadi R, Fathi E, Vietor I. Mesenchymal stem cells could be considered as a candidate for further studies in cell-based therapy of Alzheimer’s disease via targeting the signaling pathways. Acs Chem Neurosci. 2020;11(10):1424–35. https://doi.org/10.1021/acschemneuro.0c00052.
Fathi E, Vietor I. Mesenchymal stem cells promote caspase expression in Molt-4 leukemia cells via GSK-3alpha/beta and ERK1/2 signaling pathways as a therapeutic strategy. Curr Gene Ther. 2021;21(1):81–8. https://doi.org/10.2174/1566523220666201005111126.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9. https://doi.org/10.1038/nmeth.1923.
Wang L, Feng Z, Wang X, Wang X, Zhang X. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26(1):136–8. https://doi.org/10.1093/bioinformatics/btp612.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. Bmc Bioinformatics. 2011;12(1):323. https://doi.org/10.1186/1471-2105-12-323.
Fathi E, Valipour B, Sanaat Z, Nozad CH, Farahzadi R. Interleukin-6, -8, and TGF-beta secreted from mesenchymal stem cells Show functional role in reduction of telomerase activity of leukemia cell via Wnt5a/beta-catenin and P53 pathways. Adv Pharm Bull. 2020;10(2):307–14. https://doi.org/10.34172/apb.2020.037.
Fathi E, Farahzadi R, Javanmardi S, Vietor I. L-carnitine extends the telomere length of the cardiac differentiated CD117(+)- expressing stem cells. Tissue Cell. 2020;67:101429. https://doi.org/10.1016/j.tice.2020.101429.
Cianfarani F, Toietta G, Di Rocco G, Cesareo E, Zambruno G, Odorisio T. Diabetes impairs adipose tissue-derived stem cell function and efficiency in promoting wound healing. Wound Repair Regen. 2013;21(4):545–53. https://doi.org/10.1111/wrr.12051.
Sallam T, Sandhu J, Tontonoz P. Long noncoding RNA discovery in cardiovascular disease. Circ Res. 2018;122(1):155–66. https://doi.org/10.1161/CIRCRESAHA.117.311802.
Li CJ, Xiao Y, Yang M, Su T, Sun X, Guo Q, et al. Long noncoding RNA Bmncr regulates mesenchymal stem cell fate during skeletal aging. J Clin Invest. 2018;128(12):5251–66. https://doi.org/10.1172/JCI99044.
Li J, Dong J, Zhang Z, Zhang D, You X, Zhong Y, et al. miR-10a restores human mesenchymal stem cell differentiation by repressing KLF4. J Cell Physiol. 2013;228(12):2324–36. https://doi.org/10.1002/jcp.24402.
Fathi E, Charoudeh HN, Sanaat Z, Farahzadi R. Telomere shortening as a hallmark of stem cell senescence. Stem Cell Investig. 2019;6:7. https://doi.org/10.21037/sci.2019.02.04.
Wang J, Dong X, Cao L, Sun Y, Qiu Y, Zhang Y, et al. Association between telomere length and diabetes mellitus: a meta-analysis. J Int Med Res. 2016;44(6):1156–73. https://doi.org/10.1177/0300060516667132.
Fu X, Wu S, Li B, Xu Y, Liu J. Functions of p53 in pluripotent stem cells. Protein Cell. 2020;11(1):71–8. https://doi.org/10.1007/s13238-019-00665-x.
Almairac F, Turchi L, Sakakini N, Debruyne DN, Elkeurti S, Gjernes E, et al. ERK-mediated loss of miR-199a-3p and induction of EGR1 act as a “toggle switch” of GBM cell dedifferentiation into NANOG- and OCT4-positive cells. Cancer Res. 2020;80(16):3236–50. https://doi.org/10.1158/0008-5472.CAN-19-0855.
Singh K, Camera E, Krug L, Basu A, Pandey RK, Munir S, et al. JunB defines functional and structural integrity of the epidermo-pilosebaceous unit in the skin. Nat Commun. 2018;9(1):3425. https://doi.org/10.1038/s41467-018-05726-z.
Hayashi Y, Hsiao EC, Sami S, Lancero M, Schlieve CR, Nguyen T, et al. BMP-SMAD-ID promotes reprogramming to pluripotency by inhibiting p16/INK4A-dependent senescence. Proc Natl Acad Sci. 2016;113(46):13057–62. https://doi.org/10.1073/pnas.1603668113.
Min Y, Han S, Aae Ryu H, Kim S. Human adipose mesenchymal stem cells overexpressing dual chemotactic gene showed enhanced angiogenic capacity in ischaemic hindlimb model. Cardiovasc Res. 2018;114(10):1400–9. https://doi.org/10.1093/cvr/cvy086.
Morein D, Erlichman N, Ben-Baruch A. Beyond cell motility: the expanding roles of chemokines and their receptors in malignancy. Front Immunol. 2020;11. https://doi.org/10.3389/fimmu.2020.00952.
Zhu K, Zhang C, Ma X, Hu J, Cai T, Zhang L. Analyzing the interactions of mRNAs and ncRNAs to predict competing endogenous RNA networks in osteosarcoma chemo-resistance. Mol Ther. 2019;27(3):518–30. https://doi.org/10.1016/j.ymthe.2019.01.001.
Jia B, Wang Z, Sun X, Chen J, Zhao J, Qiu X. Long noncoding RNA LINC00707 sponges miR-370-3p to promote osteogenesis of human bone marrow-derived mesenchymal stem cells through upregulating WNT2B. Stem Cell Res Ther. 2019;10(1):67. https://doi.org/10.1186/s13287-019-1161-9.
Lai CH, Barik P, Hsieh DJ, Day CH, Ho TJ, Chen RJ, et al. Inhibition of cell death-inducing p53 target 1 through miR-210-3p overexpression attenuates reactive oxygen species and apoptosis in rat adipose-derived stem cells challenged with angiotensin II. Biochem Biophys Res Commun. 2020;532(3):347–54. https://doi.org/10.1016/j.bbrc.2020.07.052.
Kong R, Gao J, Ji L, Zhao D. MicroRNA-126 promotes proliferation, migration, invasion and endothelial differentiation while inhibits apoptosis and osteogenic differentiation of bone marrow-derived mesenchymal stem cells. Cell Cycle. 2020;19(17):2119–38. https://doi.org/10.1080/15384101.2020.1788258.
Lei X, He N, Zhu L, Zhou M, Zhang K, Wang C, et al. Mesenchymal stem cell-derived extracellular vesicles attenuate radiation-induced lung injury via miRNA-214-3p. Antioxid Redox Signal. 2020. https://doi.org/10.1089/ars.2019.7965.
Gu KL, Zhang Q, Yan Y, Li TT, Duan FF, Hao J, et al. Pluripotency-associated miR-290/302 family of microRNAs promote the dismantling of naive pluripotency. Cell Res. 2016;26(3):350–66. https://doi.org/10.1038/cr.2016.2.
Liu FY, Li GW, Sun CH, Chen S, Cao JF, Ma QQ, et al. Effects of bone marrow mesenchymal stem cells transfected with Ang-1 gene on hyperoxia-induced optic nerve injury in neonatal mice. J Cell Physiol. 2018;233(11):8567–77. https://doi.org/10.1002/jcp.26501.
Liu X, Yao W, Xiong H, Li Q, Li Y. LncRNA NEAT1 accelerates breast cancer progression through regulating miR-410-3p/ CCND1 axis. Cancer Biomark. 2020;29(2):277–90. https://doi.org/10.3233/CBM-190721.
Bu FT, Wang A, Zhu Y, You HM, Zhang YF, Meng XM, et al. LncRNA NEAT1: shedding light on mechanisms and opportunities in liver diseases. Liver Int. 2020;40(11):2612–26. https://doi.org/10.1111/liv.14629.
Yang Z, An Y, Wang N, Dong X, Kang H. LINC02595 promotes tumor progression in colorectal cancer by inhibiting miR-203b-3p activity and facilitating BCL2L1 expression. J Cell Physiol. 2020;235(10):7449–64. https://doi.org/10.1002/jcp.29650.
Zhang Z, Wang C, Zhang Y, Yu S, Zhao G, Xu J. CircDUSP16 promotes the tumorigenesis and invasion of gastric cancer by sponging miR-145-5p. Gastric Cancer. 2020;23(3):437–48. https://doi.org/10.1007/s10120-019-01018-7.
Chen J, Chen T, Zhu Y, Li Y, Zhang Y, Wang Y, et al. circPTN sponges miR-145-5p/miR-330-5p to promote proliferation and stemness in glioma. J Exp Clin Canc Res. 2019;38:398.
Yu L, Qu H, Yu Y, Li W, Zhao Y, Qiu G. LncRNA-PCAT1 targeting miR-145-5p promotes TLR4-associated osteogenic differentiation of adipose-derived stem cells. J Cell Mol Med. 2018;22(12):6134–47. https://doi.org/10.1111/jcmm.13892.
Shahrokhi SZ, Saeidi L, Sadatamini M, Jafarzadeh M, Rahimipour A, Kazerouni F. Can miR-145-5p be used as a marker in diabetic patients? Arch Physiol Biochem. 2020:1–6. https://doi.org/10.1080/13813455.2020.1762657.
Chen H, Huang A, He Y, Bian J, Li H. Wnt11 overexpression promote adipose-derived stem cells differentiating to the nucleus pulposus-like phenotype. Eur Rev Med Pharmaco. 2017;21:1462.
Pérez LM, de Lucas B, Gálvez BG. BMPER is upregulated in obesity and seems to have a role in pericardial adipose stem cells. J Cell Physiol. J Cell Physiol. 2021;236:132–45.
This work was supported by the National Natural Science Foundation of China (Grant No. 81772094, 81974289) and the Fundamental Research Funds for the Central Universities (HUST. 2019JYCXJJ051).
Ethics approval and consent to participate
The use of human adipose tissue and foreskin was approved by the Ethics Committee of Tongji Medical College, Huazhong University of Science and Technology (No. 2018-S288). All the young patients’ parents and adult patients have been informed in advance and consented to participate in this study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
–S4. Primers used for real-time polymerase chain reaction. Table S5. Sequences of mimics and inhibitors.
This figure illustrated the negative results of Fig. 7. The selected mRNAs, lncRNAs and circRNAs without significant difference by PCR analyses were shown in (A, B). (C, D) The correlation analysis found the expression level of miR-145-5p was not significantly correlated with the expression levels of LOC102723591 and hsa_circ_0092630.
Differentially expressed miRNAs between O-ASCs and Y-ASCs.
Differentially expressed mRNAs between O-ASCs and Y-ASCs.
Differentially expressed lncRNAs between O-ASCs and Y-ASCs.
Differentially expressed circRNAs between O-ASCs and Y-ASCs.
Differentially expressed miRNAs between D-ASCs and O-ASCs.
Differentially expressed mRNAs between D-ASCs and O-ASCs.
Differentially expressed lncRNAs between D-ASCs and O-ASCs.
Differentially expressed circRNAs between D-ASCs and O-ASCs.
Differentially expressed miRNAs between D-ASCs, O-ASCs and Y-ASCs.
Differentially expressed mRNAs between D-ASCs, O-ASCs and Y-ASCs.
Differentially expressed circRNAs between D-ASCs, O-ASCs and Y-ASCs.
Differentially expressed lncRNAs between D-ASCs, O-ASCs and Y-ASCs.
GO Biological Process Analysis between O-ASCs and Y-ASCs.
KEGG Pathway Analysis between O-ASCs and Y-ASCs.
GO Biological Process Analysis between D-ASCs and O-ASCs.
KEGG Pathway Analysis between D-ASCs and O-ASCs.
The efficiency of mimic and inhibitor transfection were examined. (A, B) Representative images of ASCs treated with 50 nM mimic-NC labeled with cy3 (red) and 200 nM inhibitor-NC labeled with 5-FAM (green). Scale bar = 50 μm. (C, D) The miR-145-5p expression level of ASCs transfected with miR-145-5p mimic or inhibitor was detected by PCR. n = 3. ***p < 0.001, ****p< 0.0001.
About this article
Cite this article
Ren, S., Xiong, H., Chen, J. et al. The whole profiling and competing endogenous RNA network analyses of noncoding RNAs in adipose-derived stem cells from diabetic, old, and young patients. Stem Cell Res Ther 12, 313 (2021). https://doi.org/10.1186/s13287-021-02388-5