- Open Access
lncRNA-mRNA expression profiles and functional networks of mesenchymal stromal cells involved in monocyte regulation
Stem Cell Research & Therapyvolume 10, Article number: 207 (2019)
The goals of this study were to explore the expression profiles and functional networks of long non-coding RNAs (lncRNAs) and messenger RNAs (mRNAs) in mesenchymal stromal cells (MSCs) involved in regulating the function of monocytes and to clarify the mechanisms by which MSCs exert immunoregulatory effects on monocytes.
MSCs and CD14+ monocytes were separately isolated. The immunoregulatory effects of MSCs on monocytes were determined by flow cytometry. lncRNAs and mRNAs that were differentially expressed (DE) between the control group (MSCs only) and co-culture group (MSCs co-cultured with monocytes) were identified through high-throughput sequencing and bioinformatic analyses and were confirmed by qRT-PCR. Bioinformatic analyses were performed to identify the critical biological functions and signalling pathways involved in MSC-mediated monocyte regulation and to identify the functional networks formed between DE mRNAs and lncRNAs.
MSCs showed a strong ability to induce monocyte migration but inhibited monocyte differentiation into M1 macrophages. A total of 145 DE lncRNAs and 768 DE mRNAs were identified between the control and co-culture groups. Significant fold changes in lncRNAs and mRNAs were confirmed by qRT-PCR. GO analysis demonstrated that DE mRNAs and lncRNAs were highly associated with terms related to binding and biological regulation. KEGG analysis revealed 122 significantly regulated pathways, including the cytokine-cytokine receptor pathway and chemokine signalling pathway. Interaction and co-expression networks were constructed for DE mRNAs and lncRNAs, and several key microRNAs were identified in the competitive endogenous RNA (ceRNA) network. Target genes of the DE lncRNAs were analysed to predict critical mRNA-lncRNA axes involved in the immunoregulatory function of MSCs.
Our research describes the lncRNA and mRNA expression profiles and functional networks involved in MSC-mediated regulation of monocytes. These results provide possible molecular mechanisms for the immunoregulatory function of MSCs and may help to elucidate possible molecular therapeutic targets in MSCs for the treatment of autoimmune diseases.
Mesenchymal stromal cells (MSCs) are cells with multiple forms of differentiation potential and self-renewal ability . These cells are derived from the mesoderm and play a role in regulating many types of immune cells , including monocytes [3, 4], dendritic cells  and T lymphocytes . These powerful features contribute to their critical role in the clinical treatment of various diseases, such as systemic lupus erythaematosus (SLE) [7, 8] and graft-versus-host disease (GVDH) [9, 10]. In addition, according to recent studies, abnormal immunoregulation by MSCs could lead to several autoimmune diseases [11,12,13]. Therefore, more in-depth studies should be conducted to clarify the specific mechanisms by which MSCs exert immunoregulatory functions. A better understanding of these mechanisms may help to improve the curative effect of MSCs and illuminate the pathogenesis of autoimmune diseases.
Monocytes, which are derived from haematopoietic stem cells, are an important type of immune cell in vivo . These cells develop in the bone marrow, migrate to inflamed tissue and then differentiate into macrophages and dendritic cells . These processes are controlled by many factors that are critical for maintaining homeostasis of the immune system in vivo [16,17,18]. MSCs can efficiently regulate monocyte migration and differentiation by secreting various cytokines and chemokines [2, 4, 19]. However, the specific immunoregulatory mechanisms by which MSCs control monocytes must be addressed.
Long non-coding RNAs (lncRNAs) are a type of non-protein coding RNA greater than 200 nt in length . lncRNAs are important epigenetic regulators and thus play crucial roles in various cell biology behaviours . Specifically, lncRNAs are widely involved in the regulation of immune system homeostasis [22, 23]. However, it is still unclear whether lncRNAs control the immunoregulatory function of MSCs in immunocytes, particularly in monocytes.
The current study presents an integrative analysis of lncRNA-mRNA expression profiles and functional networks involved in MSC-mediated regulation of monocyte functions. These results improve our understanding of the role of lncRNAs in the immunoregulatory ability of MSCs and could indicate potential targets to improve the curative effect of MSCs.
Cell isolation and culture
MSCs were purified and isolated according to our previously reported methods . Donor information is listed in Additional file 1: Table S1. Briefly, the bone marrow was extracted with a sterile bone needle, centrifuged to obtain the upper layer by density gradient, centrifuged again and then cultured in Dulbecco’s modified Eagle’s medium (DMEM; Gibco, New York, USA) containing 10% foetal bovine serum (FBS; Gibco, New York, USA). The medium was replaced once every 3 days. In all experiments, MSCs at passages 3–5 were used. Molecular markers of MSCs were detected by flow cytometry. MSC osteogenic, adipogenic and chondrogenic differentiation was confirmed by previously reported methods . Peripheral blood mononuclear cells (PBMCs) were isolated by density gradient centrifugation. CD14+ monocytes were isolated from PBMCs using CD14 MicroBeads (Miltenyi Biotec, Bergisch Gladbach, Germany) according to the manufacturer’s protocol.
Co-culture of MSCs and CD14+ monocytes
MSCs were co-cultured with CD14+ monocytes using Polycarbonate Membrane Transwell® Inserts (Corning, New York, USA). The system includes six-well plates containing inserts with a 0.4-μm pore size. 1 × 105 MSCs were cultured in the base of the wells, and 1 × 106 CD14+ monocytes were seeded in the upper inserts. Both cell types were cultured in RPMI 1640 medium (Gibco, New York, USA) containing 10% FBS at 37 °C in a 5% CO2 atmosphere. The 5.0-μm pore Transwell system was used for the migration assay.
For phenotypic analyses, MSCs were centrifuged and then resuspended in phosphate-buffered saline (PBS). Cells were then incubated with human CD29-phycoerythrin (PE), CD34-allophycocyanin (APC), CD44-fluorescein isothiocyanate (FITC), CD45-FITC, CD105-FITC or HLA-DR-PE antibodies for 30 min. To assess the purity of CD14+ monocytes, we incubated cells with the CD14-FITC antibody for 30 min. For macrophage polarization assays, CD14+ monocytes were cultured with or without MSCs for 5 days, incubated with an anti-HLA-DR-PE antibody and then incubated with fixation medium (Invitrogen) for 15 min (Additional file 4). After three washes, the cells were incubated with a permeabilization medium (Invitrogen) plus an anti-CD68 antibody (BD Pharmingen) for 30 min. For migration assays, CD14+ monocytes that migrated into the lower chambers after 12 h of co-culture were collected and counted by flow cytometry. All labelled cells were detected using a BD Influx cell sorter (BD Biosciences). All of the antibodies used for flow cytometry were purchased from BD Biosciences (New York, USA).
Library construction and high-throughput sequencing
Five MSC samples cultured with CD14+ monocytes (co-culture group; samples B1–B5) and five samples cultured without CD14+ monocytes (control group; samples A1–A5) were separately treated with TRIzol (TAKARA). RNA was collected from each sample according to the manufacturer’s protocol. RNA integrity was evaluated using the Agilent 2200 TapeStation (Agilent Technologies, USA); each sample that had an RIN above 7.0. rRNA was removed using a ribosomal RNA depletion kit, and the RNA was fragmented (average fragment length was approximately 200 nt) and reverse-transcribed into single-stranded cDNA. Then, double-stranded cDNA was synthesized, purified and treated with terminal repair and ligation primers according to the instructions from the NEBNext® Ultra™ RNA Library Prep Kit for Illumina (NEB, USA). After PCR amplification and purification, libraries were paired-end sequenced (PE150, sequencing reads were 150 bp) at Guangzhou RiboBio Co., Ltd. (Guangzhou, China) using the IlluminaHiSeq 3000 platform.
Real-time quantitative reverse transcription-polymerase chain reaction (qRT-PCR)
Total RNA was isolated from MSCs cultured with or without CD14+ monocytes using TRIzol according to the manufacturer’s protocol. Complementary DNA was transcribed using the PrimeScript RT reagent kit (TaKaRa, Dalian, China). qRT-PCR was then performed, and the results were analysed using the 2−ΔΔCt method. A detailed method can be found in our previous study . The forward and reverse primers for each gene are listed in Additional file 2: Table S2.
The raw data were filtered to remove low-quality reads, evaluate sequencing quality and remove ribosomal RNA. The resulting high-quality data were normalized to the expected number of reads per kilobase of transcript sequence per million base pairs sequenced (RPKM) and analysed using the DESeq2 method based on the negative binomial generalized linear model. Genes displaying significant fold changes (log2FoldChange > 1; q value< 0.05) between sample groups were considered for additional investigation.
mRNA and lncRNA annotation and enrichment analysis
GO term enrichment analyses were performed to acquire annotation and enrichment information. Briefly, all genes of the species were selected as background genes, and GO was performed with KOBAS3.0 software, which provides label classification of gene function and gene product attributes (http://www.geneontology.org). P values were calculated using the hypergeometric distribution method. P < 0.05 was used as the significance threshold to obtain high-frequency annotations with statistical significance relative to controls.
KEGG pathway annotation was used to categorize DE genes by biological pathway. P values were calculated using Fisher’s exact test. P < 0.05 was used as the threshold for determining statistical significance of signal transduction and disease pathway enrichment. DE mRNAs and enriched pathways were mapped using KEGG pathway annotation with KOBAS3.0 software (http://www.genome.jp/kegg).
Signal network analysis
Based on the KEGG pathway annotation and mRNA and lncRNA DE results, R was used to generate pathway nets representing the interactions between the pathways associated with DE genes.
Interaction analysis and co-expression network analysis
We used the natural language processing (NLP) method to perform text mining from the PubMed summary database in order to obtain information about all genes or proteins related to the target object. Based on the above data, we obtained the overall construction of the interaction relationship and then built an mRNA-lncRNA co-expression network using the Cytoscape software. The co-expression network was constructed by calculating the Pearson correlation coefficients and P-values between multiple genes. The transcripts were filtered using a COR of > 0.85 and a P value of < 0.05.
ceRNA network analysis
The mRNAs and lncRNAs selected for the co-expression network analysis were the same as those used to predict miRNA targets using the miRbase. The miRNAs obtained through these predictions were screened using the miRanda and TargetScan programmes. lncRNAs and mRNAs possessing microRNA recognition elements (MREs) for the targeted miRNAs were predicted using RNA22. The competitive endogenous RNA (ceRNA) network was constructed and illustrated using Cytoscape (v3.4.0).
Statistical analyses were performed using SPSS 22.0 software (Chicago, IL, USA). All qPCR results are expressed as the mean ± standard deviation (SD). Spearman correlation was used to determine the relationship between lncRNAs and their target genes. P values < 0.05 were considered significant.
MSC-mediated regulation of CD14+ monocytes
MSCs were isolated and purified as previously described . Flow cytometric analyses demonstrated that all MSCs were positive for CD29, CD44 and CD105 and negative for CD34, CD45 and HLA-DR (Fig. 1a). In addition, MSCs were capable of differentiating into osteoblasts, adipoblasts and chondroblasts (Fig. 1b). The purity of isolated CD14+ monocytes was 99.0% (Fig. 1c). Through polarization assays, we demonstrated that the ratio of CD14+ monocytes that differentiated into M1 macrophages was significantly lower when these cells were co-cultured with MSCs compared with those cultured without MSCs (Fig. 1d). Additionally, migration assays showed that CD14+ monocyte migration was increased upon co-culture with MSCs (Fig. 1e). These results indicated that MSCs potently regulate CD14+ monocytes, promoting CD14+ monocyte migration but inhibiting M1 macrophage polarization.
Identification of differentially expressed lncRNAs and mRNAs
A total of 768 mRNAs were differentially expressed (DE) in MSCs co-cultured with CD14+ monocytes compared to MSCs cultured alone. Among these genes, 461 mRNAs were upregulated and 307 mRNAs were downregulated. The DE mRNAs are depicted using a clustergram (Fig. 2a) and volcano plots (Fig. 2c). The 20 mRNAs with the largest fold changes are shown in Table 1. Several cytokines are contained within this list, including CCL8, CXCL6, CCL20, CXCL5, CXCL8 and CXCL3, which are secreted by MSCs and are important for regulating the function of CD14+ monocytes. A total of 145 lncRNAs, including 102 upregulated and 43 downregulated lncRNAs, were DE in MSCs co-cultured with CD14+ monocytes compared to MSCs cultured alone. The DE lncRNAs are depicted in a clustergram (Fig. 2b) and volcano plots (Fig. 2d). The 10 lncRNAS with the largest fold changes are shown in Table 2.
Validation of DE mRNA and lncRNA expression levels
To confirm the results of RNA sequencing, several important DE mRNAs and lncRNAs were assessed by qPCR. As shown in Fig. 3a, CXCL2, CXCL3, CXCL6, CXCL8, IL-6, CCL2, MMP3, CFB, LIF, TNFAIP6 and SOCS3 were upregulated in MSCs co-cultured with CD14+ monocytes, while GREM2, DKK1, OSR1, CTGF and ROR1 were downregulated (Fig. 3b). Moreover, qPCR results demonstrated that LINC00473, ENSG00000231083.1, LOC101928674, ENSG00000237927.1, ENSG00000232949.1 and PDZNR3-AS1 were upregulated, while LINC01111, ENSG00000235513.1, PLCE1-AS1, ZFHX4-AS1 and LINC01279 were downregulated in MSCs co-cultured with CD14+ monocytes (Fig. 3c, d). All qPCR results were consistent with the RNA sequencing results, confirming the reliability of the sequencing data.
GO and KEGG analysis
We performed GO analysis on the DE mRNAs and lncRNAs. The top 10 GO terms related to biological processes, cellular components and molecular function are shown in Fig. 4a and Table 3. In the biological process domain, the top 5 GO terms associated with DE mRNAs were single-organism process, single-organism cellular process, cellular process, biological regulation and response to stimulus. In the cellular component domain, the top 5 GO terms associated with DE mRNAs were cell, cell part, intracellular, intracellular part and cytoplasm. In the molecular function domain, the top 5 GO terms associated with DE mRNAs were binding, protein binding, receptor binding, catalytic activity and carbohydrate derivative binding. KEGG analysis of the DE mRNAs determined that 122 pathways were significantly altered in MSCs co-cultured with CD14+ monocytes. The top 30 affected pathways are shown in Fig. 4b. The top 10 pathways and DE mRNAs associated with these pathways are shown in Table 4. The top pathways included cytokine-cytokine receptor interactions, TNF signalling pathway, chemokine signalling pathway and NF-kappa B signalling pathway, which contribute to the immunoregulatory function of MSCs, were identified.
Interaction and co-expression network analysis
The interactions between proteins coded by DE mRNAs are shown in Fig. 5a. CCL2, IL1B, COL7A1 and ICAM1 are the key genes that interacted with many other DE mRNAs in this network. In addition, a co-expression network was constructed for DE lncRNAs and mRNAs (Fig. 5b). Within these RNA networks, LOC101929122 has the maximum number of targets, including 28 DE mRNAs, and WNT5A has the maximum number of co-expressed lncRNAs. The top 10 co-expression pairs are shown in Table 5. lncRNAs act as ceRNAs to exert their biological function ; thus, we constructed a ceRNA network from all of the DE mRNAs and lncRNAs in order to identify possible regulatory miRNAs (Fig. 5c). MiR-939-5p, miR-940 and miR-8075 were enriched in the ceRNA network, indicating that they play critical roles in the DE mRNA-lncRNA functional network and in the immunoregulatory function of MSCs.
Analysis of lncRNA target mRNAs
Next, we analysed the possible target genes of the DE lncRNAs. A total of 1076 target genes of these DE lncRNAs were identified. Target genes with a combined score higher than 0.9 are shown in Fig. 6a. Figure 6a shows that HIF1A was the target gene of three DE lncRNAs, including NR_047116.1, NR_045406.1 and NR_144368.1. Moreover, IL6, an important cytokine, was the target of the DE lncRNA NR_131935.1. A Venn diagram analysis demonstrated that 56 mRNAs were found at the intersection of the 768 DE mRNAs and the 1076 DE lncRNA target genes (Fig. 6b). A list of these 56 mRNAs is shown in Additional file 3: Table S3. The list includes several known inflammatory-related molecules, such as IL6, SOCS3 and IDO1. We also performed KEGG pathway analysis on the DE lncRNA target genes to identify pathway clusters. The protein-protein interactions between these targets and the relationships between these pathways are shown in Fig. 6c, including the TNF signalling pathway, cytokine-cytokine receptor interactions and the HIF-1 signalling pathway.
In our present research, we utilized high-throughput sequencing followed by bioinformatic analysis to analyse the mRNA and lncRNA expression profiles and functional networks of MSCs co-cultured with CD14+ monocytes. These findings were then confirmed by qPCR. KEGG pathway analysis indicated that some key pathways, such as cytokine-cytokine receptor interactions, the TNF signalling pathway, the chemokine signalling pathway and the NF-kappa B signalling pathway, contribute to the immunoregulatory function of MSCs in monocytes. Other bioinformatic analyses revealed interactions between mRNAs, lncRNAs and miRNAs. Additionally, we identified target genes of the DE lncRNAs as well as the intersection of these genes with DE mRNAs and performed KEGG pathway analysis on the DE genes. Our results provide a model that can be used to explore the roles of lncRNAs and mRNAs in MSC-mediated immunoregulatory mechanisms.
The homeostasis of the immune system depends on the balance between immunocytes and their regulatory cells . MSCs are one of the most important immunoregulatory cell types [13, 28]. MSCs are pluripotent stem cells that regulate the functions of many immune cells, including T lymphocytes, B lymphocytes and dendritic cells [4, 19, 29]. Previous studies have demonstrated that MSCs can affect monocytes and macrophages in a paracrine manner in a co-culture system. For example, Ko et al.  found that MSC-primed monocytes/macrophages express high levels of MHC class II, B220, CD11b and IL-10 and exhibit T cell–suppressive activities in a TNF-α–stimulated gene/protein (TSG)-6–dependent manner. Similar findings are described in Melief et al. . This group showed that MSCs can promote monocyte survival and induce differentiation towards macrophage type 2 cells, which express CD206 and CD163 and secrete high levels of IL-10 and CCL-18. This powerful immunoregulatory ability allows MSCs to play wide-ranging roles in maintaining the dynamic equilibrium of the immune system in vivo. Monocytes are critical constituents of the immune system . Monocytes are active immunocytes that migrate to specific tissues and then differentiate into macrophages and dendritic cells upon stimulation by various factors . MSCs play an essential role in accelerating monocyte migration and inhibiting M1 macrophage differentiation; these observations were confirmed in this study [32,33,34]. On the one hand, MSCs exhibit a considerable therapeutic effect in several diseases through the secretion of cytokines that modulate monocyte function [9, 28, 35]. On the other hand, abnormal immunoregulation by MSCs could lead to several autoimmune diseases [11,12,13]. Therefore, it is important to elucidate a detailed mechanism of MSC-mediated immunoregulation of monocytes in order to improve their curative effects and to illuminate the pathogenesis of autoimmune diseases.
mRNA expression profiles reflect the biological behaviours and functions of cells. mRNA expression profiles are altered when cells are co-cultured with other cells. These changes in mRNA expression profiles are closely related to their functions in other cells . To explore the mechanism of MSC functions in monocytes, we analysed the mRNA profile by gene sequencing and bioinformatics analysis. This analysis identified 768 DE mRNAs. Among these DE mRNAs, CCL8, CXCL6, CCL20, CXCL5, CXCL8, CXCL3 and CXCL10 were the most significantly upregulated. These cytokines play important roles in the migration and polarization of monocytes [37,38,39,40]. These results indicate that MSCs affect monocyte migration and polarization mainly through these key cytokines. In addition, GO analysis demonstrated enrichment of specific molecular functions, such as cytokine activity, indicating the critical role of these cytokines in the immunoregulatory ability of MSCs. KEGG analysis identified 122 signalling pathways associated with the DE mRNAs. Cytokine-cytokine receptor interactions, the TNF signalling pathway, the chemokine signalling pathway and the NF-kappa B signalling pathway were the pathways with the largest significant differences. In a previous study, NF-kappa B-mediated activation of MSCs led to the secretion of TNF-α and other cytokines, which was related to MSC-based therapeutic efficacy . Therefore, we suggest that the NF-kappa B signalling pathway is activated in MSCs co-cultured with monocytes. This activation leads to the secretion of cytokines, including TNF-α, CCL and CXCL, which ultimately regulate the migration and polarization of monocytes. Notably, the SLE and rheumatoid arthritis pathways were also among the pathways identified as having the largest significant difference by KEGG analysis. These results are consistent with those of a previous study and emphasize the importance of MSC-mediated immunoregulation of monocytes in the pathogenesis and clinical treatment of autoimmune diseases [42, 43].
mRNA expression profiles are under the control of a series of epigenetic regulators, of which lncRNAs are an important component . A lncRNA is a transcript longer than 200 nt that does not have protein coding ability but exerts indispensable regulatory functions . Numerous studies have confirmed that lncRNAs regulate the differentiation of MSCs [45, 46]. For example, the lncRNA GAS5 negatively regulates MSC adipogenesis through the miR-18a/CTGF axis . However, the specific lncRNAs involved in the immunoregulatory capacity of MSCs remain unknown. In this study, we first detected lncRNA expression profiles involved in MSC immunoregulatory functions. We found that 145 lncRNAs were DE in MSCs upon co-culture with monocytes. These DE lncRNAs may contribute to the differential expression of mRNAs and thus may be involved in the immunoregulatory function of MSCs. One of the most critical roles of lncRNAs is their ability to act as ceRNAs for miRNAs . In this study, we constructed a ceRNA network between DE mRNAs and lncRNAs, which allowed us to identify several key miRNAs, including miR-939-5p, miR-940 and miR-8075. MiR-940 promotes the osteogenic differentiation of MSCs . Additionally, exosomal miR-940 induces M2 polarization of monocytes . Here, we suggest that a DE lncRNA-miR-940-mRNA network exists in MSCs and plays an active role in their immunoregulatory function.
To further study the mechanism of these lncRNAs, we identified 1076 target mRNAs of the DE lncRNAs. KEGG analysis identified a total of 33 enriched pathways, including the TNF signalling pathway and cytokine-cytokine receptor interaction pathway, further emphasizing their importance in our study. However, it is more important to determine whether these targets overlap with DE mRNAs. Using a Venn diagram analysis, we identified 56 genes that were included in both the DE mRNA cluster and the target mRNA of the DE lncRNA cluster. One of these 56 genes was IL6, which is a cytokine that participates in the immunoregulatory function of MSCs . Our results indicate that DE lncRNA NR_131935.1 and its target gene, DE mRNA IL6 (shown in Fig. 6a), could comprise another functional axis important for MSC-mediated immunoregulation of monocytes.
Ankylosing spondylitis (AS) is an autoimmune disease characterized by chronic inflammation and pathological osteogenesis . Abnormal MSC immunoregulation contributes to the pathogenesis of autoimmune diseases [51, 52]. Recently, we demonstrated that CCL2 is overexpressed in MSCs from AS patients, inducing the monocyte disorder and eventually leading to the chronic inflammation characteristic of AS . However, the mechanism and underlying cause of CCL2 overexpression remain unclear. CCL2, also called MCP1, is a secreted protein involved in immunoregulatory and inflammatory processes . In the present study, MSC CCL2 expression was significantly upregulated upon co-culture with monocytes, confirming that CCL2 is involved in MSC-mediated immunoregulation of monocytes, as previously reported . In addition, SOCS3, a regulator of CCL2 expression, was identified within the DE mRNA cluster and is the target gene of a DE lncRNA, LOC101928674. Therefore, we suggest that a LOC101928674-SOCS3-CCL2 axis may exist and may be involved in MSC-mediated immunoregulation of monocytes. We further suggest that abnormal LOC101928674 expression may contribute to CCL2 overexpression in MSCs from AS patients and that a LOC101928674-SOCS3-CCL2 axis may be involved in the pathogenesis of AS. However, further research should be performed to confirm this hypothesis.
In summary, our research presents lncRNA and mRNA expression profiles and functional networks of MSCs involved in the regulation of monocytes. Several candidate mechanisms were identified, including the LOC101928674-SOCS3-CCL2 axis. Our results provide possible molecular mechanisms for the immunoregulatory function of MSCs, which may help to elucidate the pathogenesis of autoimmune diseases and improve the clinical use of MSCs. However, the specific lncRNAs that function in MSCs and their associated mechanisms are still unknown. Further studies should address these questions.
Our research describes the MSC lncRNA and mRNA expression profiles and functional networks involved in the regulation of monocytes. Our results provide possible molecular mechanisms for the immunoregulatory function of MSCs. Several candidate mechanisms were discussed, including the LOC101928674-SOCS3-CCL2 axis. These proposed mechanisms may help to elucidate possible molecular therapeutic targets for MSC-based autoimmune diseases.
Availability of data and materials
Dulbecco’s modified Eagle medium
Foetal bovine serum
Long non-coding RNA
MicroRNA recognition elements
Mesenchymal stromal cell
Natural language processing
Peripheral blood mononuclear cells
Real-time quantitative reverse transcription-polymerase chain reaction
Systemic lupus erythaematosus
Dominici M, Le Blanc K, Mueller I, Slaper-Cortenbach I, Marini F, Krause D, et al. Minimal criteria for defining multipotent mesenchymal stromal cells. The International Society for Cellular Therapy position statement. Cytotherapy. 2006;8:315–7.
Bernardo ME, Fibbe WE. Mesenchymal stromal cells: sensors and switchers of inflammation. Cell Stem Cell. 2013;13:392–402.
Francois M, Romieu-Mourez R, Li M, Galipeau J. Human MSC suppression correlates with cytokine induction of indoleamine 2,3-dioxygenase and bystander M2 macrophage differentiation. Mol Ther. 2012;20:187–95.
Amouzegar A, Mittal SK, Sahu A, Sahu SK, Chauhan SK. Mesenchymal stem cells modulate differentiation of myeloid progenitor cells during inflammation. Stem Cells. 2017;35:1532–41.
Chiesa S, Morbelli S, Morando S, Massollo M, Marini C, Bertoni A, et al. Mesenchymal stem cells impair in vivo T-cell priming by dendritic cells. Proc Natl Acad Sci U S A. 2011;108:17384–9.
Rashedi I, Gomez-Aristizabal A, Wang XH, Viswanathan S, Keating A. TLR3 or TLR4 activation enhances mesenchymal stromal cell-mediated treg induction via notch signaling. Stem Cells. 2017;35:265–75.
Wang D, Zhang H, Liang J, Wang H, Hua B, Feng X, et al. A long-term follow-up study of allogeneic mesenchymal stem/stromal cell transplantation in patients with drug-resistant systemic lupus erythematosus. Stem Cell Rep. 2018;10:933–41.
Zhu Y, Feng X. Genetic contribution to mesenchymal stem cell dysfunction in systemic lupus erythematosus. Stem Cell Res Ther. 2018;9:149.
Fujii S, Miura Y, Fujishiro A, Shindo T, Shimazu Y, Hirai H, et al. Graft-versus-host disease amelioration by human bone marrow mesenchymal stromal/stem cell-derived extracellular vesicles is associated with peripheral preservation of naive T cell populations. Stem Cells. 2018;36:434–45.
Dotoli GM, De Santis GC, Orellana MD, de Lima Prata K, Caruso SR, Fernandes TR, et al. Mesenchymal stromal cell infusion to treat steroid-refractory acute GvHD III/IV after hematopoietic stem cell transplantation. Bone Marrow Transplant. 2017;52:859–62.
He Y, Xu LL, Feng FE, Wang QM, Zhu XL, Wang CC, et al. Mesenchymal stem cell deficiency influences megakaryocytopoiesis through the TNFAIP3/NF-kappaB/SMAD pathway in patients with immune thrombocytopenia. Br J Haematol. 2018;180:395–411.
Geng L, Tang X, Zhou K, Wang D, Wang S, Yao G, et al. MicroRNA-663 induces immune dysregulation by inhibiting TGF-beta1 production in bone marrow-derived mesenchymal stem cells in patients with systemic lupus erythematosus. Cell Mol Immunol. 2018;16:260–74.
Pistoia V, Raffaghello L. Mesenchymal stromal cells and autoimmunity. Int Immunol. 2017;29:49–58.
Shi C, Pamer EG. Monocyte recruitment during infection and inflammation. Nat Rev Immunol. 2011;11:762–74.
Jakubzick CV, Randolph GJ, Henson PM. Monocyte differentiation and antigen-presenting functions. Nat Rev Immunol. 2017;17:349–62.
Tabas I, Lichtman AH. Monocyte-macrophages and T cells in atherosclerosis. Immunity. 2017;47:621–34.
Nahrendorf M. Myeloid cell contributions to cardiovascular health and disease. Nat Med. 2018;24:711–20.
Murray PJ. Immune regulation by monocytes. Semin Immunol. 2018;35:12–8.
Ma S, Xie N, Li W, Yuan B, Shi Y, Wang Y. Immunobiology of mesenchymal stem cells. Cell Death Differ. 2014;21:216–25.
Ulitsky I, Bartel DP. lincRNAs: genomics, evolution, and mechanisms. Cell. 2013;154:26–46.
Kopp F, Mendell JT. Functional classification and experimental dissection of long noncoding RNAs. Cell. 2018;172:393–407.
Wang P, Xue Y, Han Y, Lin L, Wu C, Xu S, et al. The STAT3-binding long noncoding RNA lnc-DC controls human dendritic cell differentiation. Science. 2014;344:310–3.
Du M, Yuan L, Tan X, Huang D, Wang X, Zheng Z, et al. The LPS-inducible lncRNA Mirt2 is a negative regulator of inflammation. Nat Commun. 2017;8:2049.
Wang P, Li Y, Huang L, Yang J, Yang R, Deng W, et al. Effects and safety of allogenic mesenchymal stem cell intravenous infusion in active ankylosing spondylitis patients who failed NSAIDs: a 20-week clinical trial. Cell Transplant. 2014;23:1293–303.
Li M, Xie Z, Wang P, Li J, Liu W, Tang S, et al. The long noncoding RNA GAS5 negatively regulates the adipogenic differentiation of MSCs by modulating the miR-18a/CTGF axis as a ceRNA. Cell Death Dis. 2018;9:554.
Xie Z, Sa T, Ye G, Wang P, Li J, Liu W, et al. Interleukin-6/interleukin-6 receptor complex promotes osteogenic differentiation of bone marrow-derived mesenchymal stem cells. Stem Cell Res Ther. 2018;9:13.
Thomson DW, Dinger ME. Endogenous microRNA sponges: evidence and controversy. Nat Rev Genet. 2016;17:272–83.
Galipeau J, Sensebe L. Mesenchymal stromal cells: clinical challenges and therapeutic opportunities. Cell Stem Cell. 2018;22:824–33.
Shi Y, Wang Y, Li Q, Liu K, Hou J, Shao C, et al. Immunoregulatory mechanisms of mesenchymal stem and stromal cells in inflammatory diseases. Nat Rev Nephrol. 2018;14:493–507.
Ko JH, Lee HJ, Jeong HJ, Kim MK, Wee WR, Yoon SO, et al. Mesenchymal stem/stromal cells precondition lung monocytes/macrophages to produce tolerance against allo- and autoimmunity in the eye. Proc Natl Acad Sci U S A. 2016;113:158–63.
Melief SM, Schrama E, Brugman MH, Tiemessen MM, Hoogduijn MJ, Fibbe WE, et al. Multipotent stromal cells induce human regulatory T cells through a novel pathway involving skewing of monocytes toward anti-inflammatory macrophages. Stem Cells. 2013;31:1980–91.
Jiang XX, Zhang Y, Liu B, Zhang SX, Wu Y, Yu XD, et al. Human mesenchymal stem cells inhibit differentiation and function of monocyte-derived dendritic cells. Blood. 2005;105:4120–6.
Spaggiari GM, Abdelrazik H, Becchetti F, Moretta L. MSCs inhibit monocyte-derived DC maturation and function by selectively interfering with the generation of immature DCs: central role of MSC-derived prostaglandin E2. Blood. 2009;113:6576–83.
Morrison TJ, Jackson MV, Cunningham EK, Kissenpfennig A, McAuley DF, O'Kane CM, et al. Mesenchymal stromal cells modulate macrophages in clinically relevant lung injury models by extracellular vesicle mitochondrial transfer. Am J Respir Crit Care Med. 2017;196:1275–86.
Panes J, Garcia-Olmo D, Van Assche G, Colombel JF, Reinisch W, Baumgart DC, et al. Expanded allogeneic adipose-derived mesenchymal stem cells (Cx601) for complex perianal fistulas in Crohn's disease: a phase 3 randomised, double-blind controlled trial. Lancet. 2016;388:1281–90.
Han Z, Wang J, Gao L, Wang Q, Wu J. Aberrantly expressed messenger RNAs and long noncoding RNAs in degenerative nucleus pulposus cells co-cultured with adipose-derived mesenchymal stem cells. Arthritis Res Ther. 2018;20:182.
Asano K, Takahashi N, Ushiki M, Monya M, Aihara F, Kuboki E, et al. Intestinal CD169(+) macrophages initiate mucosal inflammation by secreting CCL8 that recruits inflammatory monocytes. Nat Commun. 2015;6:7802.
Metcalf TU, Wilkinson PA, Cameron MJ, Ghneim K, Chiang C, Wertheimer AM, et al. Human monocyte subsets are transcriptionally and functionally altered in aging in response to pattern recognition receptor agonists. J Immunol. 2017;199:1405–17.
Yoshida K, Korchynskyi O, Tak PP, Isozaki T, Ruth JH, Campbell PL, et al. Citrullination of epithelial neutrophil-activating peptide 78/CXCL5 results in conversion from a non-monocyte-recruiting chemokine to a monocyte-recruiting chemokine. Arthritis Rheumatol. 2014;66:2716–27.
O'Kane CM, Boyle JJ, Horncastle DE, Elkington PT, Friedland JS. Monocyte-dependent fibroblast CXCL8 secretion occurs in tuberculosis and limits survival of mycobacteria within macrophages. J Immunol. 2007;178:3767–76.
Zhang Y, Chiu S, Liang X, Gao F, Zhang Z, Liao S, et al. Rap1-mediated nuclear factor-kappaB (NF-κB) activity regulates the paracrine capacity of mesenchymal stem cells in heart repair following infarction. Cell Death Discov. 2015;1:15007.
Liu Y, Mu R, Wang S, Long L, Liu X, Li R, et al. Therapeutic potential of human umbilical cord mesenchymal stem cells in the treatment of rheumatoid arthritis. Arthritis Res Ther. 2010;12:R210.
Barbado J, Tabera S, Sanchez A, Garcia-Sancho J. Therapeutic potential of allogeneic mesenchymal stromal cells transplantation for lupus nephritis. Lupus. 2018;27:2161–5.
Atianand MK, Caffrey DR, Fitzgerald KA. Immunobiology of long noncoding RNAs. Annu Rev Immunol. 2017;35:177–98.
Wang L, Wu F, Song Y, Li X, Wu Q, Duan Y, et al. Long noncoding RNA related to periodontitis interacts with miR-182 to upregulate osteogenic differentiation in periodontal mesenchymal stem cells of periodontitis patients. Cell Death Dis. 2016;7:e2327.
Huang Y, Zheng Y, Jia L, Li W. Long noncoding RNA H19 promotes osteoblast differentiation via TGF-beta1/Smad3/HDAC signaling pathway by deriving miR-675. Stem Cells. 2015;33:3481–92.
Hashimoto K, Ochi H, Sunamura S, Kosaka N, Mabuchi Y, Fukuda T, et al. Cancer-secreted hsa-miR-940 induces an osteoblastic phenotype in the bone metastatic microenvironment via targeting ARHGAP1 and FAM134A. Proc Natl Acad Sci U S A. 2018;115:2204–9.
Chen X, Ying X, Wang X, Wu X, Zhu Q, Wang X. Exosomes derived from hypoxic epithelial ovarian cancer deliver microRNA-940 to induce macrophage M2 polarization. Oncol Rep. 2017;38:522–8.
Deng Y, Zhang Y, Ye L, Zhang T, Cheng J, Chen G, et al. Umbilical cord-derived mesenchymal stem cells instruct monocytes towards an IL10-producing phenotype by secreting IL6 and HGF. Sci Rep. 2016;6:37566.
Taurog JD, Chhabra A, Colbert RA. Ankylosing spondylitis and axial spondyloarthritis. N Engl J Med. 2016;374:2563–74.
Tang Y, Ma X, Zhang H, Gu Z, Hou Y, Gilkeson GS, et al. Gene expression profile reveals abnormalities of multiple signaling pathways in mesenchymal stem cell derived from patients with systemic lupus erythematosus. Clin Dev Immunol. 2012;2012:826182.
Bocelli-Tyndall C, Bracci L, Spagnoli G, Braccini A, Bouchenaki M, Ceredig R, et al. Bone marrow mesenchymal stromal cells (BM-MSCs) from healthy donors and auto-immune disease patients reduce the proliferation of autologous- and allogeneic-stimulated lymphocytes in vitro. Rheumatology (Oxford). 2007;46:403–8.
Xie Z, Wang P, Li J, Li Y, Wang S, Wu X, et al. MCP1 triggers monocyte dysfunctions during abnormal osteogenic differentiation of mesenchymal stem cells in ankylosing spondylitis. J Mol Med (Berl). 2017;95:143–54.
Zhang K, Luo J. Role of MCP-1 and CCR2 in alcohol neurotoxicity. Pharmacol Res. 2018;139:360–6.
Howe CL, LaFrance-Corey RG, Goddery EN, Johnson RK, Mirchia K. Neuronal CCL2 expression drives inflammatory monocyte infiltration into the brain during acute virus infection. J Neuroinflammation. 2017;14:238.
The authors thank Yuxiao Ye and American Journal Experts for providing English language editing of the manuscript.
This work was supported by grants from the National Natural Science Foundation of China (81702120, 81672128, 81672097) and the Science and Technology Project of Guangdong Province (2015B020228001).
Ethics approval and consent to participate
Sample collection was approved by the Ethics Committee of Sun Yat-sen Memorial Hospital, Sun Yat-Sen University. The bone marrow MSCs used in this study were obtained from the Center for Biotherapy, Sun Yat-sen Memorial Hospital, Sun Yat-Sen University.
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.
Table S1. Donor information. (DOCX 14 kb)
Table S2. Primers used for qRT-PCR. (DOCX 17 kb)
Table S3. The intersection of DE mRNAs and DE lncRNA target genes. (DOCX 13 kb)
Figure S1. The proportion of M2 macrophages in 0 day and 5 days after spontaneous differentiation and the flow diagram of TNFα and IL1β by MSCs when co-culture or not with monocytes. (JPG 309 kb)