- Open Access
Multiomics global landscape of stemness-related gene clusters in adipose-derived mesenchymal stem cells
Stem Cell Research & Therapy volume 11, Article number: 310 (2020)
Adipose-derived mesenchymal stem cells (AD-MSCs) are a type of stem cell that is abundant and widely used. The molecular characteristics of AD-MSCs from different passages from donors of different ages have not been well elucidated.
Six kinds of AD-MSCs ((E1, E2, E3, Y1, Y2, and Y3) with E denoting cells derived from an elderly patient, Y denoting cells derived from a young patient, and 1, 2, and 3 representing passages 3, 6, and 10) were obtained from human abdominal adipose tissue. We obtained the protein expression profile, the mRNA expression profile, the lncRNA expression profile, and the methylation profile of each kind of AD-MSC by sequencing. After calculating the stemness indices, genes related to stemness were extracted. The multiomics correlation analysis was performed in the stemness-related genes. In addition, short time-series expression miner (STEM) analysis was performed for all cell passages and donor ages. To further explore the biological functions of the stemness-related genes, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Finally, the lncRNA-KEGG network and transcription factor (TF)-KEGG network were constructed based on the RNAInter database and TRRUST v2 database.
The stemness of the Y1, E1, and Y2 cells was higher than that of the E2, Y3, and E3 cells. The stemness was the highest for Y1 cells and the lowest for E3 cells. STEM analysis showed that five stemness-related gene clusters were associated with the cell passages, and only one gene cluster was associated with age. The enrichment analysis results showed that the biological processes (BPs) and KEGG pathways were mainly involved in the proliferation, differentiation, and migration of cells. The global regulatory landscape of AD-MSCs was constructed: 25 TFs and 16 lncRNAs regulated 21 KEGG pathways through 27 mRNAs. Furthermore, we obtained a core stemness-related gene set consisting of ITGAV, MAD2L1, and PCNA. These genes were expressed at higher levels in Y1 cells than in E3 cells.
The multiomics global landscape of stemness-related gene clusters was determined for AD-MSCs, which may be helpful for selecting AD-MSCs with increased stemness.
When used for different therapeutic procedures aimed at treating various forms of tissue damage, the autologous fat graft is indeed a solution. Recent studies have shown that adipose-derived stromal vascular fractions (AD-SVFs) [1, 2] and adipose-derived stem cells (ASCs) [3, 4] can enhance tissue regeneration potential . ASCs are localized in AD-SVFs, which have a heterogeneous mesenchymal cell set  and can be obtained by 2 methods: enzymatic digestion and mechanical filtration [7, 8]. ASCs have the ability to differentiate into cells of mesenchymal origin in vitro; this includes osteoblasts, adipocytes, and chondrocytes. They can also enhance different tissues in vivo, including cartilage, muscle, bone, and fat . As shown in previous studies, adipose-derived human follicle stem cells (AD-HFSCs) are helpful for hair regrowth . Adipose-derived mesenchymal stem cells (AD-MSCs) contribute to wound healing and soft tissue defects . In summary, SVFs and AD-MSCs are ideal tools for generating medicine.
Some studies have shown that platelet-rich plasma (PRF) and hyaluronic acid (HA) dressings alone or in combination are also helpful to promote tissue repair [10,11,12]. PRF can promote tissue repair in two ways by providing a bridge for tissue incorporation and new blood vessel formation [13, 14] as well as producing growth factors (b-FGF, PDGF, VEGF, EGF, TGF-β, and IGF-1) . When AD-MSCs are injected into damaged parts of the human body, they can promote tissue regeneration through their own proliferation and differentiation. Furthermore, growth factors can accelerate this process by driving cell growth . Interestingly, AD-SVFs/AD-MSCs/fat grafts are used in combination with PRF to help improve wound healing and soft tissue defects [16, 17]. More importantly, some studies reported that the Wnt/TGF-β/β-catenin signaling pathway is necessary for the growth of cells in wound healing and soft tissue defects [9, 18, 19].
The stemness of AD-MSCs is mainly reflected in their ability to maintain self-renewal, cell differentiation, and proliferation . A previous study showed that the differentiation capacity is maintained with aging; however, AD-MSCs from younger donors may exhibit a higher proliferation rate . In addition, long-term cell culture with serial passaging may negatively affect the stemness of stem cells . These findings indicate that the stemness of AD-MSCs is influenced by the ages of donors and the number of cell passages. A previous study indicated that CD13, CD29, CD44, CD73, CD90, CD105, and CD106 in AD-MSCs are upregulated in expansion culture compared to those in the stromal vascular fraction (SVF) [23, 24]. In addition, the senescence-related proteins p53, p21, and p16 are strongly expressed . These studies show that AD-MSCs from different passages may present different molecular characteristics. However, the molecular characteristics of AD-MSCs from different passages from donors of different ages have not been well elucidated.
This study attempts to elucidate the multiomics and molecular characteristics of AD-MSCs derived from donors of different ages and from different cell passages. We found that the stemness of AD-MSCs from younger donors was increased. In addition, the stemness of early-passage (passage 3) AD-MSCs was the highest compared with that of intermediate-passage (passage 6) and late-passage (passage 10) AD-MSCs. These findings may provide guidance on how to select AD-MSCs with increased stemness, as well as whether AD-MSCs should be injected directly or post-expanded.
Sample acquisition and cell culture
The study was performed with the permission of the ethics committee of the First People’s Hospital of Nanning, Guangxi Zhuang Autonomous Region. We obtained written informed consent from donors. Human lipoaspirate from the abdominal subcutaneous tissue of healthy females (one young patient who was 27 years old and an elderly patient who was 63 years old) was stored for less than 48 h at 4 °C before processing. The adipose tissue (AT) was gently aspirated using 10-ml Luer-Lok syringes and then placed in a tube, which was centrifuged for 3 min at 1200g. Then, we added collagenase solution (type II; Worthington Biochemical Corp., Lakewood, N.J.) to the AT, which was shaken gently in a 37 °C water bath to digest the adipose tissue. The digested AT was centrifuged again for 10 min and then filtered to remove the debris at 180g. Next, erythrocyte lysis buffer was added to the cell pellet (stromal vascular part), which was resuspended and centrifuged for 10 min at 180g. The stromal vascular fraction was plated and cultured in tissue culture-treated flasks in adipose-derived stem cell plating medium (0.001% dexamethasone, Dulbecco’s Modified Eagle Medium/Nutrient Mixture F12, 10% fetal bovine serum, and 1% Fungizone (Bristol-Myers Squibb, New York, N.Y.) and 1% penicillin/streptomycin. After overnight incubation, the nonadherent cells were removed, and the medium was changed to fresh medium. After 20 to 48 h, the cells were expanded until they were nearly confluent to obtain the passage 0 AD-MSCs. The expansion rate and morphology of the AD-MSCs were observed under a microscope. When the AD-MSCs were at a confluence of 80% to 90%, they were passaged at a ratio of 1:3. The cell number of passage 0 AD-MSCs was 1 × 106 and 2 × 106 in young patient and the elderly patient, respectively. When amplifying passages 3, 6, and 10, the cell number was 1 × 108, 1 × 1010, and 2 × 1012 in the young patient. While in the elderly patient, the cell number was 1 × 108, 5 × 109, and 2 × 1011 when amplifying passage 3, 6, and 10. Subsequently, 1 × 108 cells from each kind AD-MSCs were collected for sequencing. AD-MSCs from the two patients were harvested at three time points during amplification for multiomics analysis to obtain the early-passage (passage 3), intermediate-passage (passage 6), and late-passage (passage 10) cells. A total of 6 types of AD-MSCs (E1, E2, E3, Y1, Y2, and Y3) were obtained; E represents the elderly patient, Y represents the young patient, and 1, 2, and 3 represent passages 3, 6, and 10. The workflow of this study is shown in Fig. 1.
RNA extraction and sequencing
All RNA was isolated from AD-MSCs by TRIzol purification (Invitrogen), and genomic DNA was removed with the help of gDNA eliminator columns from the RNeasy Mini Kit (Qiagen). In addition, spectrophotometry (NanoDrop, Thermo Fisher Scientific) was used to assess the RNA quality along with a Bioanalyzer (Agilent). The RNA was enriched for mRNA (Dynabeads mRNA purification kit, Invitrogen) by poly(A) isolation, and the mRNA was fragmented and subjected to first-strand cDNA synthesis. After second-strand cDNA synthesis, the double-stranded cDNA was quantified and purified. The library was constructed by using a protocol from Illumina with the NEBNext DNA Sample Prep Reagent kits (NEB). The double-stranded cDNA was end-repaired, and the adapters were ligated to the DNA fragments; poly(A) tracts RNA sequences were then added. After size selection (200~600 bp) and UDGase treatment to ensure strand specificity, the adapter-modified DNA fragments were amplified by PCR. Next, 76-base paired-end sequencing was performed on an Illumina instrument.
The strand-specific paired-end reads were screened for ribosomal RNA by alignment against known rRNA sequences (RefSeq) using Bowtie software . The rRNA-subtracted paired-end reads were aligned with TopHat . The maximum multibit was set to 1 m, and a micro-exon search was performed. In addition, a gene model was provided as a gene transfer format (GTF) file (Ensembl BDGP5.25.60). The aligned reads in valid pairs were subjected to FPKM estimation through Cuffilinks. Correction and bias detection were performed at this step, which only retained fragments compatible with Ensembl annotation (BDGP5.25.60). These reads were counted toward the number of mapped hits used in the FPKM denominator.
DNA methylation assay
LCM-DNAs from AD-MSCs were fragmented into 100–500 bp fragments. Then, the genomic DNA was extracted from the AD-MSCs using a high-salt procedure, and the HELP assay was performed as previously described. The assay involved interrogating the cytosine methylation status and performing comparative isoschizomer profiling on a genomic scale. The genomic DNA from AD-MSCs was digested by HpaII (a methylcytosine-sensitive enzyme) and MspI, and then the MspI- and HpaII-produced fragments were amplified through PCR mediated by ligation. Both of the amplified fractions were submitted to Roche-NimbleGen, Inc. (Madison, WI). Then, these fractions were labeled and hybridized with a human hg17 custom-designed oligonucleotide array (50-mers) covering 25,626 HpaII-amplifiable fragments (HAFs), which were located at gene promoters. HpaII-amplifiable fragments are defined as genomic sequences. Each fragment on the array was represented by 15 individual probes distributed randomly throughout the microarray slide. The microarray covered 50,000 CpGs corresponding to 14,000 gene promoters. The signal intensities were calculated as the robust mean of their component probe-level signal intensities after intensive quality control. The log2(HpaII/MspI) value represented methylation and was analyzed as a continuous variable. For most loci, each fragment was categorized as either hypomethylated if the log ratio was greater than zero or methylated if the centered log HpaII/MspI ratio was less than zero. The average methylation level for all CpG islands in a gene was considered the methylation level of the corresponding gene in the present study.
Each sample was homogenized with 5 volumes of glass sand. Then, 200 μl radioimmunoprecipitation assay (RIPA) buffer containing phenylmethylsulfonyl fluoride (PMSF) was added to lyse the cells. Then, the samples were sonicated. All soluble proteins were collected after centrifugation at 13,400g for 20 min. The concentration of the protein was measured by a quantitative protein kit (2D-quant-kit, GEHealthcare).
One hundred micrograms of protein was reduced, alkylated, and precipitated by the chloroform precipitation method to digest the protein. First, 55 mM Reducing Reagent (8-plex iTRAQ kit, AB Sciex, USA) was added to the proteinase and incubated at 60 °C for 60 min. Then, the protein was mixed with Cysteine Blocking Reagent (8-plex iTRAQ kit, AB Sciex, USA). After deionized water and 70% ethanol were added to the 10 KD ultrafiltration cartridge, the protein solution was poured into the ultrafiltration cartridge. The solution was centrifuged for 20 min at 13,400g. The ultrafiltration cartridge was centrifuged in 0.25 M TEAB (triethylammonium bicarbonate). The protein pellets were reconstituted in 0.25 M urea/50 mM TEAB and digested with 2% trypsin overnight (Promega).
A high-performance liquid chromatography (HPLC) system (Phenomenex columns; Gemini-NX 3u C18 110A; 150*2.00 mM) was used to perform strong cation exchange (SCX) fractionation chromatography. The isobaric tags for relative and absolute quantification (iTRAQ)-labeled peptides were separated with a linear gradient formed by mobile phase B (80% ACN, 20 mM HCOONH4, pH 10) and mobile phase A (20 mM HCOONH4, pH 10). The flow for peptide elution was set to a rate of 20 μl/min. The fractions were collected and acidified with trifluoroacetic acid (50%). The fractions were vacuum-dried prior to LC-MS/MS. The fractions were dissolved in buffer and pelleted. The supernatant was loaded onto analytical columns and identified with a Q Exactive system (Thermal Scientific). The components of the mobile phase used for LC-MS were formic acid (0.1%) and 80% CAN. The flow rate of the analytical columns was set at 350 nl/min. The peptides were analyzed by using a 3-step gradient for 65 min. The parameters of the first-grade MS were a maximum injection time of 40 ms and a resolution of 70,000. The scan range was from 350 to 1800 m/z. The second-grade MS spectra were acquired at a resolution of 17,500. The top 20 precursors were selected for each MS cycle.
The raw counts of the mapped reads were aggregated using featureCounts . The gene-level quantification was performed with a gene transfer format (GTF) file (UCSC Genome Bioinformatics: Frequently Asked Questions: Data File Forma. https://genome.ucsc.edu/FAQ/FAQformat.html#format4. Accessed on 12 January 2016.). Then, to obtain the mRNA/lncRNA expression profiles, the data were normalized using the “voom” function of the limma package  in R. The protein expression profile and the methylation profile were constructed according to the proteome sequencing data and the methylation sequencing data, respectively, without standardization.
Calculation of stemness
Stemness is considered to represent the potential for differentiation and self-renewal of the cell of origin, which possesses the ability to produce all cell types in the adult organism . The TCGAbiolinks package  in R was used to calculate the stemness for the six kinds of AD-MSCs based on mRNA expression using the messenger RNA (mRNA) expression profiles.
The correlation between stemness and mRNA/long noncoding RNA (lncRNA) expression/protein expression was used to identify mRNAs/lncRNAs/proteins that are highly related to stemness. The correlation between stemness and mRNA/lncRNA/DNA methylation was also explored. The Pearson correlation coefficient (r) was calculated using the Hmisc package  in R. P < 0.01 and r > 0.9 were considered to indicate a strong correlation. The correlation among coding gene (mRNA) DNA methylation, protein expression, and corresponding mRNA expression and the long noncoding gene (lncRNA) DNA methylation and corresponding lncRNA expression were also explored.
Short time-series expression miner (STEM) analysis
The mRNAs highly related to stemness were used to perform STEM analysis using STEM v1.3.8 . The mRNAs were organized into different profiles (clusters) based on expression patterns using STEM analysis. Donor age (young and elderly) and the AD-MSCs passage number (passages 3, 6, and 10) were considered time points.
Functional enrichment analysis and gene set variation analysis (GSVA)
To explore the biological functions of the stemness-related gene profiles, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analyses were performed using the clusterProfiler package  in R. Then, the most significant (ranked by P value) biological processes (BPs) and KEGG pathways were visualized with a bubble diagram. GSVA  was used to score the individual samples according to the BPs and KEGG pathways, and each sample received a GSVA score. In addition, KEGG and BP enrichment analysis were also performed with the ClueGO plug-in  in Cytoscape software . The GO and KEGG networks were used to determine the relationships between the GO terms according to the similarity of their related genes. A P value adjusted by the false discovery rate < 0.05 was considered significant.
Stemness-related transcriptional regulatory network
Interactions between stemness-related lncRNAs and their target genes were downloaded from the RNAInter database . The interactions between transcription factors (TFs) and stemness-related mRNAs were downloaded from the TRRUST v2 database . Correlation analysis between the lncRNA/TFs and their targets was performed to reduce the number of false positives and noise. Subsequently, in addition to the enrichment analysis, an AD-MSCs stemness-related lncRNA/TF-target-KEGG pathway network was constructed.
Multiomics atlas of genes related to AD-MSC stemness
A total of 19,189 mRNAs, 21,965 lncRNAs, and 3205 proteins were included in these analyses. The DNA methylation levels of 8070 mRNAs and 4308 lncRNAs were measured. Stemness was closely associated with the cell passage of the AD-MSCs in the following order of decreasing association: early passage > intermediate passage > late passage (Fig. 2a). A total of 1286 mRNAs, 90 proteins, and 804 lncRNAs were found to be highly related to stemness. The DNA methylation levels of 50 mRNAs and 29 lncRNAs were highly related to stemness (Fig. 2b). However, there were only 20 mRNAs and corresponding proteins that were highly related to stemness (Fig. 2c). In addition, the expression and DNA methylation level of the lncRNA STARD13 were both highly related to stemness. The expression of STARD13 was positively related to stemness, while the methylation level of STARD13 DNA was negatively related to stemness (Fig. 2d).
Clusters of stemness-related genes associated with AD-MSC passage or donor age
In the STEM analysis, although 15 mRNA profiles were identified, only 5 mRNA profiles including a total of 555 genes were considered significant (Table S1). Two mRNA profiles (0 and 3) showed significantly gradually upregulated expression patterns that occurred during cell proliferation and passaging. Three mRNA profiles (11, 12, and 15) showed significantly gradually downregulated expression patterns (Fig. 3a). An mRNA profile comprising 35 genes (Table S2) was identified as representing expression patterns that were significantly upregulated with age (Fig. 3b).
Stemness-related gene clusters involved in multiple functional pathways
The five AD-MSCs passage-related mRNA profiles, which included 555 genes, were used to perform functional enrichment analysis. The enrichment results showed that there were 2032 associated BPs (Table S3) and 56 associated KEGG pathways (Table S4). Unsurprisingly, multiple BPs, such as the PI3K-Akt signaling pathway , p53 signaling pathway , and apoptosis  (Fig. 4a), were mainly related to stem cells. The genes in profiles 11 and 15 were involved in some BPs related to the immune response and were negatively related to stemness. The five AD-MSCs passage-related mRNA profiles were involved in proliferation-, differentiation-, and migration-related pathways (Fig. 4b). In the ClueGO analysis, a total of 218 BPs were identified. The results of the ClueGO analysis also suggested that the stemness-related mRNA profiles were associated with cell proliferation-related BPs, such as DNA replication, the cell cycle phase transition, and regulation of the cell cycle (Fig. 4c). In addition, the stemness-related mRNA profiles were involved in DNA replication, necroptosis, the cell cycle, systemic lupus erythematosus, and homologous recombination pathways (Fig. 4d).
Global regulatory landscape of stemness-related gene clusters at the multiomics level
The interactions of stemness-related lncRNAs and mRNAs were extracted from the RNAInter database, and the interactions of TFs and stemness-related mRNAs were extracted from the TRRUST v2 database. By combining these interactions with the enrichment analysis, sixteen stemness-related lncRNAs were implicated in nine KEGG pathways by regulating 13 stemness-related mRNAs in profiles 0, 3, and 12 (Fig. 5a). Fifteen TFs that may regulate 15 stemness-related mRNAs in profiles 0, 3, and 11 were involved in 19 KEGG pathways (Fig. 5b). Subsequently, the AD-MSCs stemness-related lncRNA/TF-target-KEGG pathway networks were constructed and included 16 lncRNAs, 15 TFs, 27 mRNAs, and 21 KEGG pathways (Fig. 5c). Previous studies proposed that the pathways involving AD-MSCs are mainly involved in differentiation, proliferation, and migration . Therefore, we found 9 pathways of interest for further study, including DNA replication, cell cycle, homologous recombination, cellular senescence, apoptosis, PI3K-Akt signaling, ECM-receptor interaction, cytokine-cytokine receptor interaction, and p53 signaling pathways. In DNA replication, the MCM complex unwinds the double-strand parental DNA. Subsequently, with the help of DNA polymerase (POLA1 and POLD1), the DNA clamp (PCNA), and the clamp loader (EFC4), the new DNA is synthesized. We found that CDK1 could affect the cell cycle in the p53 signaling pathway. In addition, the specific interaction between cells and ECMs is mediated by integrins containing ITGAV via the cytokine-cytokine receptor interaction and extracellular matrix (ECM)-receptor interaction. In apoptosis, BIRC5 and TUBA1B could be affected, thereby affecting cell apoptosis. In homologous recombination, BLM is helpful to maintain genome integrity. In cellular senescence, the p21 complex containing CCNA2 and CCNB1 may be affected as well as its downstream regulator (E2F1). Eventually, the cell cycle is affected. In the cell cycle, the Maps1 complex containing MAD2L1 and the complex containing PCNA are affected by the DNA damage checkpoint; therefore, CCNA2 is affected (Fig. 5d).
Multiomics characteristics of core stemness-related genes and their lncRNA regulators
Forty genes were shown to be highly related to both mRNA and protein expression. Three genes were identified by combining pathways of interest, including ITGAV, MAD2L1, and PCNA. These genes were defined as the core stemness genes. In addition, the 3 genes with the highest similarity were ITGAV, MAD2L1, and PCNA (Fig. 6a), which further verified our findings. Among them, only MAD2L1 had a methylation profile. As shown in Fig. 2, the stemness of AD-MSCs was the highest for Y1 cells and the lowest for E1 cells. Therefore, Y1 and E1 cells were chosen to determine the expression of the core stemness genes. Compared to those of E1 cells, the mRNA and protein expression profiles of Y1 cells showed that MAD2L1 was upregulated, but MAD2L1 was downregulated according to the methylation profile. Compared to E1 cells, Y1 cells showed the upregulation of PCNA according to the protein and mRNA expression profiles. The expression profile of ITGAV was exactly the opposite of that of PCNA (Fig. 6b).
AD-MSCs, which have a high potential for differentiation into many kinds of mature cells in an appropriate environment, are considered to be an ideal treatment in regenerative therapies . In the current study, we found that stemness of AD-MSCs is closely associated with the cell passage in the decreasing order of early passage > intermediate passage > late passage. The potential biological functions of stemness-related mRNAs were highly related to proliferation, differentiation, and migration. In addition, we explored the regulatory network of AD-MSCs.
In this study, we found that the stemness of Y1, E1, and Y2 cells was higher than that of E2, Y3, and E3 cells. This is consistent with previous studies, which showed that as the age and the number of cell passages increased, the proliferation and differentiation potential of AD-MSCs decreased [44, 45]. In addition, 5 stemness gene clusters were obtained, which showed that the maintenance of stemness in AD-MSCs was the result of the interaction of multiple genes.
The comprehensive regulatory network indicated that lncRNAs and TFs could influence mRNA expression in AD-MSCs, thereby affecting their phenotype. Among the 9 pathways that were analyzed further, some had been reported in previous studies related to stem cells, such as the cell cycle , homologous recombination , apoptosis , PI3K-Akt signaling , and p53 signaling pathways . In addition, we found that DNA replication, cellular senescence, cytokine-cytokine receptor interactions, and ECM-receptor interactions were related to stem cells.
Furthermore, we found that the genes in the above pathways were widely involved in the regulation of the cell cycle and intercellular functions in cells in the current study. MCM5  and CCNA2  play roles in the regulation of the mammalian cell cycle. In addition, MCM2 plays an important role in cell division and DNA replication . Moreover, the protein encoded by PCNA plays a central role in recruiting and retaining many of the enzymes required for DNA replication and repair . POLA1 , POLD1 , PCNA , and RFC4  can participate in DNA replication by forming DNA polymerase. MAD2L1 is required for execution of the mitotic checkpoint . CDK1 is a catalytic subunit of a protein kinase complex that induces cell entry into mitosis . CCNB1 is expressed predominantly in the G2/M phase of cell division . CDC45 is an important component of the replication fork, participating in DNA unwinding . TUBA1B is involved in mitosis, cell movement, intracellular movement, and other biological processes . When activated in response to DNA damage, E2F1 can promote proliferation or apoptosis . BIRC5 is an inhibitor of apoptosis that can regulate apoptosis . BLM is a tumor suppressor that maintains genome integrity . ITGAV is a kind of α-V integrin that functions in cell surface adhesion . According to our study, some of the differentially expressed genes were related to proliferation, differentiation, and metastasis in AD-MSCs.
In this study, we identified 3 core stemness-related genes (ITGAV, MAD2L1, and PCNA) that comprise a core stemness-related gene set. As shown in previous studies, MAD2L1, PCNA, and ITGAV are related to the proliferation and differentiation of cells. ITGAV is a kind of mRNA, and its product belongs to the integrin α chain family, which serves as a major receptor for differentiation and cell proliferation . MAD2L1 mainly functions during mitosis to ensure that all chromosomes properly align at the metaphase plate . The protein encoded by PCNA is a cofactor of DNA polymerase δ and is involved in the processes of DNA replication and DNA repair . Therefore, the core stemness-related gene set identified in this study may be helpful to assess the stemness of AD-MSCs.
Our study presents several limitations. First, the study is based on bioinformatics prediction; therefore, further experimental verification is needed. In addition, the sample size involved in this study was small, and we need to verify the results in a larger data cohort later.
The multiomics global landscape of stemness-related gene clusters was constructed for AD-MSCs. Moreover, we identified a core stemness-related gene set comprised of ITGAV, PCNA, and MAD2L, which may be helpful for selecting AD-MSCs with increased stemness.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Adipose-derived mesenchymal stem cells
Kyoto Encyclopedia of Genes and Genomes
Stromal vascular fraction
Gene transfer format
High-performance liquid chromatography
Strong cation exchange
Long non-coding RNA
Gene set variation analysis
Adipose-derived stem cells
Adipose-derived human follicle stem cells
Isobaric tags for relative and absolute quantification
Gentile P, Kothari A, Casella D, et al. Fat graft enhanced with adipose-derived stem cells in aesthetic breast augmentation: clinical, histological, and instrumental evaluation. Aesthet Surg J. 2019.sjz292.
Gentile P, Scioli MG, Bielli A, et al. Comparing different nanofat procedures on scars: role of the stromal vascular fraction and its clinical implications. Regen Med. 2017;12(8):939–52.
Gentile P, Casella D, Palma E, et al. Engineered fat graft enhanced with adipose-derived stromal vascular fraction cells for regenerative medicine: clinical, histological and instrumental evaluation in breast reconstruction. J Clin Med. 2019;8(4):504.
Gentile P, Garcovich S. Concise review: adipose-derived stem cells (ASCs) and adipocyte-secreted exosomal microRNA (A-SE-miR) modulate cancer growth and promote wound repair. J Clin Med. 2019;8(6):855.
Gentile P, Piccinno M S, Calabrese C. Characteristics and potentiality of human adipose-derived stem cells (hASCs) obtained from enzymatic digestion of fat graft. Cells. 2019;8(3):282.
Mitchell JB, Mcintosh K, Zvonic S, et al. Immunophenotype of human adipose-derived cells: temporal changes in stromal-associated and stem cell-associated markers. Stem Cells. 2006;24(2):376–85.
Gentile P, Scioli MG, Orlandi A, et al. Breast reconstruction with enhanced stromal vascular fraction fat grafting: what is the best method? Plast Reconstr Surg Glob Open. 2015;3(6):e406.
Gentile P, Calabrese C, de Angelis B, et al. Impact of the different preparation methods to obtain human adipose-derived stromal vascular fraction cells (AD-SVFs) and human adipose-derived mesenchymal stem cells (AD-MSCs): enzymatic digestion versus mechanical centrifugation. Int J Mol Sci. 2019;20(21):5471.
Gentile P. Autologous cellular method using micrografts of human adipose tissue derived follicle stem cells in androgenic alopecia. Int J Mol Sci. 2019;20(14):3446.
Cervelli V, Lucarini L, Spallone D, et al. Use of platelet-rich plasma and hyaluronic acid in the loss of substance with bone exposure. Adv Skin Wound Care. 2011;24(4):176–81.
Nicoli F, Balzani A, Lazzeri D, et al. Severe hidradenitis suppurativa treatment using platelet-rich plasma gel and Hyalomatrix. Int Wound J. 2015;12(3):338–43.
Gentile P, Bottini DJ, Spallone D, et al. Application of platelet-rich plasma in maxillofacial surgery: clinical evaluation. J Craniofac Surg. 2010;21(3):900–4.
Gentile P, Colicchia GM, Nicoli F, et al. Complex abdominal wall repair using a porcine dermal matrix. Surg Innov. 2013;20(6):NP12–5.
Scioli MG, Bielli A, Gentile P, et al. Combined treatment with platelet-rich plasma and insulin favours chondrogenic and osteogenic differentiation of human adipose-derived stem cells in three-dimensional collagen scaffolds. J Tissue Eng Regen Med. 2017;11(8):2398–410.
Gentile P, Scioli MG, Bielli A, et al. Stem cells from human hair follicles: first mechanical isolation for immediate autologous clinical use in androgenetic alopecia and hair loss. Stem Cell Investig. 2017;4:58.
Cervelli V, Bocchini I, di Pasquali C, et al. P.R.L. platelet rich lipotransfert: our experience and current state of art in the combined use of fat and PRP. Biomed Res Int. 2013;2013:434191.
Cervelli V, Gentile P. Use of cell fat mixed with platelet gel in progressive hemifacial atrophy. Aesthet Plast Surg. 2009;33(1):22–7.
Gentile P, Scioli MG, Cervelli V, et al. Autologous micrografts from scalp tissue: trichoscopic and long-term clinical evaluation in male and female androgenetic alopecia. Biomed Res Int. 2020;2020:7397162.
Gentile P, Garcovich S. Advances in regenerative stem cell therapy in androgenic alopecia and hair loss: Wnt pathway, growth-factor, and mesenchymal stem cell signaling impact analysis on cell growth and hair follicle development. Cells. 2019;8(5):466.
el Atat O, Antonios D, Hilal G, et al. An evaluation of the stemness, paracrine, and tumorigenic characteristics of highly expanded, minimally passaged adipose-derived stem cells. PLoS One. 2016;11(9):e0162332.
Shi YY, Nacamuli RP, Salim A, et al. The osteogenic potential of adipose-derived mesenchymal cells is maintained with aging. Plast Reconstr Surg. 2005;116(6):1686–96.
Kumar A, Xu Y, Yang E, et al. Fidelity of long-term cryopreserved adipose-derived stem cells for differentiation into cells of ocular and other lineages. Exp Eye Res. 2019;189:107860.
Pachon-Pena G, Yu G, Tucker A, et al. Stromal stem cells from adipose tissue and bone marrow of age-matched female donors display distinct immunophenotypic profiles. J Cell Physiol. 2011;226(3):843–51.
Zhang X, Hirai M, Cantero S, et al. Isolation and characterization of mesenchymal stem cells from human umbilical cord blood: reevaluation of critical factors for successful isolation and high ability to proliferate and differentiate to chondrocytes as compared to mesenchymal stem cells from bone marrow and adipose tissue. J Cell Biochem. 2011;112(4):1206–18.
Jin HJ, Bae YK, Kim M, et al. Comparative analysis of human mesenchymal stem cells from bone marrow, adipose tissue, and umbilical cord blood as sources of cell therapy. Int J Mol Sci. 2013;14(9):17986–8001.
Langmead B, Trapnell C, Pop M, et al. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):R25.
Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25(9):1105–11.
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30.
Ritchie ME, Phipson B, Wu D, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
Malta TM, Sokolov A, Gentles AJ, et al. Machine learning identifies stemness features associated with oncogenic dedifferentiation. Cell. 2018;173(2):338–54 e15.
Mounir M, Lucchetta M, Silva TC, et al. New functionalities in the TCGAbiolinks package for the study and integration of cancer data from GDC and GTEx. PLoS Comput Biol. 2019;15(3):e1006701.
Koziol JA, Jia Z. The concordance index C and the Mann-Whitney parameter Pr(X>Y) with randomly censored data. Biom J. 2009;51(3):467–74.
Ernst J, Bar-Joseph Z. STEM: a tool for the analysis of short time series gene expression data. BMC Bioinformatics. 2006;7:191.
Yu G, Wang LG, Han Y, et al. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.
Hanzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics. 2013;14:7.
Bindea G, Mlecnik B, Hackl H, et al. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. 2009;25(8):1091–3.
Kohl M, Wiese S, Warscheid B. Cytoscape: software for visualization and analysis of biological networks. Methods Mol Biol. 2011;696:291-303.
Lin Y, Liu T, Cui T, et al. RNAInter in 2020: RNA interactome repository with increased coverage and annotation. Nucleic Acids Res. 2020;48(D1):D189-D97.
Han H, Cho JW, Lee S, et al. TRRUST v2: an expanded reference database of human and mouse transcriptional regulatory interactions. Nucleic Acids Res. 2018;46(D1):D380–D6.
Song F, Jiang D, Wang T, et al. Mechanical stress regulates osteogenesis and adipogenesis of rat mesenchymal stem cells through PI3K/Akt/GSK-3beta/beta-catenin signaling pathway. Biomed Res Int. 2017;2017:6027402.
Okuda A, Kurokawa S, Takehashi M, et al. Poly (ADP-ribose) polymerase inhibitors activate the p53 signaling pathway in neural stem/progenitor cells. BMC Neurosci. 2017;18(1):14.
Fulda S. Regulation of apoptosis pathways in cancer stem cells. Cancer Lett. 2013;338(1):168–73.
Richardson SM, Kalamegam G, Pushparaj PN, et al. Mesenchymal stem cells in regenerative medicine: focus on articular cartilage and intervertebral disc regeneration. Methods. 2016;99:69–80.
Dominici M, le Blanc K, Mueller I, et al. Minimal criteria for defining multipotent mesenchymal stromal cells. The International Society for Cellular Therapy position statement. Cytotherapy. 2006;8(4):315–7.
Stenderup K, Justesen J, Clausen C, et al. Aging is associated with decreased maximal life span and accelerated senescence of bone marrow stromal cells. Bone. 2003;33(6):919–26.
Mens MMJ, Ghanbari M. Cell cycle regulation of stem cells by microRNAs. Stem Cell Rev Rep. 2018;14(3):309–22.
Mujoo K, Pandita RK, Tiwari A, et al. Differentiation of human induced pluripotent or embryonic stem cells decreases the DNA damage repair by homologous recombination. Stem Cell Reports. 2017;9(5):1660–74.
Tsuruga H, Yabuta N, Hashizume K, et al. Expression, nuclear localization and interactions of human MCM/P1 proteins. Biochem Biophys Res Commun. 1997;236(1):118–25.
Murphy M, Stinnakre MG, Senamaud-Beaufort C, et al. Delayed early embryonic lethality following disruption of the murine cyclin A2 gene. Nat Genet. 1997;15(1):83–6.
Mincheva A, Todorov I, Werner D, et al. The human gene for nuclear protein BM28 (CDCL1), a new member of the early S-phase family of proteins, maps to chromosome band 3q21. Cytogenet Cell Genet. 1994;65(4):276–7.
Baple EL, Chambers H, Cross HE, et al. Hypomorphic PCNA mutation underlies a human DNA repair disorder. J Clin Invest. 2014;124(7):3137–46.
Toukoki C, Gryllos I. PolA1, a putative DNA polymerase I, is coexpressed with PerR and contributes to peroxide stress defenses of group A Streptococcus. J Bacteriol. 2013;195(4):717–25.
Palles C, Cazier JB, Howarth KM, et al. Germline mutations affecting the proofreading domains of POLE and POLD1 predispose to colorectal adenomas and carcinomas. Nat Genet. 2013;45(2):136–44.
Boehm EM, Gildenberg MS, Washington MT. The many roles of PCNA in eukaryotic DNA replication. Enzymes. 2016;39:231–54.
Xiang J, Fang L, Luo Y, et al. Levels of human replication factor C4, a clamp loader, correlate with tumor progression and predict the prognosis for colorectal cancer. J Transl Med. 2014;12:320.
Li Y, Benezra R. Identification of a human mitotic checkpoint gene: hsMAD2. Science. 1996;274(5285):246–8.
Al Jord A, Shihavuddin A, Servignat D'aout R, et al. Calibrated mitotic oscillator drives motile ciliogenesis. Science. 2017;358(6364):803–6.
Gong D, Ferrell JE Jr. The roles of cyclin A2, B1, and B2 in early and late mitotic events. Mol Biol Cell. 2010;21(18):3149–61.
Fenwick AL, Kliszczak M, Cooper F, et al. Mutations in CDC45, encoding an essential component of the pre-initiation complex, cause Meier-Gorlin syndrome and craniosynostosis. Am J Hum Genet. 2016;99(1):125–38.
Aiken J, Moore JK, Bates EA. TUBA1A mutations identified in lissencephaly patients dominantly disrupt neuronal migration and impair dynein activity. Hum Mol Genet. 2019;28(8):1227–43.
Morris EJ, Ji JY, Yang F, et al. E2F1 represses beta-catenin transcription and is antagonized by both pRB and CDK8. Nature. 2008;455(7212):552–6.
Caldas H, Holloway MP, Hall BM, et al. Survivin-directed RNA interference cocktail is a potent suppressor of tumour growth in vivo. J Med Genet. 2006;43(2):119–28.
Wang Y, Li S, Smith K, et al. Intrachromosomal recombination between highly diverged DNA sequences is enabled in human cells deficient in Bloom helicase. DNA Repair (Amst). 2016;41:73–84.
Kim H, Wrann CD, Jedrychowski M, et al. Irisin mediates effects on bone and fat via alphaV integrin receptors. Cell. 2018;175(7):1756–68 e17.
Wu Y, Tan L, Chen J, et al. MAD2 combined with mitotic spindle apparatus (MSA) and anticentromere antibody (ACA) for diagnosis of small cell lung cancer (SCLC). Med Sci Monit. 2018;24:7541–7.
Prestel A, Wichmann N, Martins JM, et al. The PCNA interaction motifs revisited: thinking outside the PIP-box. Cell Mol Life Sci. 2019;76(24):4923–43.
This work was financially supported by the National Nature Science Foundation of China (81760346, 81860341, 81771271), the National Key R&D Program of China (2018YFA0108304), the Guangxi Natural Science Foundation (2018GXNSFAA281148, 2019JJA140071), the Scientific Research & Technology Development Program of Nanning (20173021-2, 20183037-1, 20191034, 20193100), and the Yong River Program of innovation and entrepreneurship of Nanning (2018-01-07); The First Batch of High-level Talent Scientific Research Projects of the Affiliated Hospital of Youjiang Medical University for Nationalities in 2019 (Contract No. R20196307).
Ethics approval and consent to participate
The study was conducted with the approval of the ethics committee of First People’s Hospital of Nanning, Guangxi Zhuang Autonomous Region. We obtained written informed consent from donors.
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.
. The 555 genes belonging to the 5 stemness-related gene clusters associated with cell passages.
. The thirty-five genes belonging to the single stemness-related gene cluster associated with age.
. GO_BPs enriched in adipose-derived mesenchymal cells. GO: gene ontology, BPs: biological processes.
. KEGG pathways enriched in AD-MSCs. AD-MSCs: adipose-derived mesenchymal cells.
About this article
Cite this article
Lu, GM., Rong, YX., Liang, ZJ. et al. Multiomics global landscape of stemness-related gene clusters in adipose-derived mesenchymal stem cells. Stem Cell Res Ther 11, 310 (2020). https://doi.org/10.1186/s13287-020-01823-3
- Multiomics global landscape
- Adipose-derived mesenchymal stem cells
- Stemness gene clusters
- Short time-series miner analysis