Involvement of microglia-expressed MS4A6A in the onset of glioblastoma.

Glioblastoma multiforme (GBM) represents the deadliest form of brain tumour, characterized by its low survival rate and grim prognosis. Cytokines released from glioma-associated microglia/macrophages are involved in establishing the tumour microenvironment, thereby crucially promoting GBM progression. MS4A6A polymorphism was confirmed to be associated with neurodegenerative and polymorphism disease pathobiology, but whether it participates in the regulation of GBM and the underlying mechanisms is still not elucidated. Here, we found that MS4A6A was significantly upregulated in GBM patient samples. The results from the single-cell RNA-sequencing (scRNA-seq) database and immunostaining demonstrated the specific expression of MS4A6A in microglial cells. In vitro, microglial overexpression of MS4A6A stimulated the proliferation and migration of glioblastoma cells. Moreover, high MS4A6A mRNA expression was related to poor prognosis in GBM patients. Our study highlights the potential of MS4A6A as a promising biomarker for GBM, which may provide novel strategies for its prevention, diagnosis and treatment.

Targeting cytohesin-1 suppresses acute myeloid leukemia progression and overcomes resistance to ABT-199.

Adhesion molecules play essential roles in the homeostatic regulation and malignant transformation of hematopoietic cells. The dysregulated expression of adhesion molecules in leukemic cells accelerates disease progression and the development of drug resistance. Thus, targeting adhesion molecules represents an attractive anti-leukemic therapeutic strategy. In this study, we investigated the prognostic role and functional significance of cytohesin-1 (CYTH1) in acute myeloid leukemia (AML). Analysis of AML patient data from the GEPIA and BloodSpot databases revealed that CYTH1 was significantly overexpressed in AML and independently correlated with prognosis. Functional assays using AML cell lines and an AML xenograft mouse model confirmed that CYTH1 depletion significantly inhibited the adhesion, migration, homing, and engraftment of leukemic cells, delaying disease progression and prolonging animal survival. The CYTH1 inhibitor SecinH3 exerted in vitro and in vivo anti-leukemic effects by disrupting leukemic adhesion and survival programs. In line with the CYTH1 knockdown results, targeting CYTH1 by SecinH3 suppressed integrin-associated adhesion signaling by reducing ITGB2 expression. SecinH3 treatment efficiently induced the apoptosis and inhibited the growth of a panel of AML cell lines (MOLM-13, MV4-11 and THP-1) with mixed-lineage leukemia gene rearrangement, partly by reducing the expression of the anti-apoptotic protein MCL1. Moreover, we showed that SecinH3 synergized with the BCL2-selective inhibitor ABT-199 (venetoclax) to inhibit the proliferation and promote the apoptosis of ABT-199-resistant leukemic cells. Taken together, our results not only shed light on the role of CYTH1 in cell-adhesion-mediated leukemogenesis but also propose a novel combination treatment strategy for AML.

Proteomic characteristics of saliva in patients with different subgroups of IgG4-RD.

Background: Immunoglobulin G4-related disease (IgG4-RD) is a newly defined disease entity, with great heterogeneity among IgG4-RD subgroups with different organ involvement patterns. Identification of the proteomic characteristics of IgG4-RD subgroups will be critical for the understanding of the pathogenic mechanisms of IgG4-RD. Method: In this study, we performed proteomic analysis using Tandem Mass Tags (TMT) technology with "high field" mass analyzer with improved resolution and sequencing speed to investigate the proteomic profile of saliva and plasma samples from ten untreated IgG4-RD patients and five healthy controls (HCs). Differentially expressed proteins (DEPs) were identified by "t test" function in R package. Functional enrichment analysis was used to investigate pathways enriched in IgG4-RD samples. Results: Most salivary DEPs identified in IgG4-RD patients compared with HCs were mainly enriched in neutrophil mediated GO bioprocess. Within the comparisons between four IgG4-RD subgroups, more DEPs were identified in the comparison of Mikulicz group and Head and neck group. Among four subgroups of IgG4-RD, Head and neck group showed the most distinctive proteomic expression pattern when compared with HCs. Moreover, "Neutrophil mediated process" related GO bioprocess was commonly identified between comparisons of Mikulicz group and Head and neck group, Head and neck group and Retroperitoneal aorta group, Head and neck group and HCs, IgG4-RD patients with saliva gland involvement and those without saliva gland involvement. Key DEPs that involved in this GO bioprocess were identified. Besides, we performed proteomic analysis for plasma samples between ten IgG4-RD and five HCs and there were several DEPs identified overlapped in saliva and plasma. Conclusion: We identified multiple processes/factors and several signaling pathways in saliva that may be involved in the IgG4-RD pathogenesis.

The landscape of T and B lymphocytes interaction and synergistic effects of Th1 and Th2 type response in the involved tissue of IgG4-RD revealed by single cell transcriptome analysis.

OBJECTIVES: To investigate the landscape of T-B cell interaction, immune receptor profiles and effects of different types of immune responses in the involved tissues of IgG4-RD. METHODS: Single cell RNA sequencing, bulk sample RNA sequencing, immune receptor repertoire analysis (both BCR and TCR), multi-color flow cytometry, and in-vitro assays with model cells (e.g. EBV-immortalized B cells from IgG4-RD patient) and histologic methods were applied to investigate the immunopathological features of IgG4-RD from multiple aspects. RESULTS: Ectopic germinal center formation was observed in IgG4-RD patients at advanced disease stage, and a large part of B cells in involved tissue were germinal center B cell-like. Germinal center reaction in IgG4-RD led to the irregularities of both TCR and BCR clones in the involved tissues, and limited clonal overlaps among different samples. Enhanced Th1- and Th2-type responses were observed in involved tissues of IgG4-RD and patients with both increased Th1- and Th2-type response related cell subsets possessed more severe inflammatory indices. Analyses to the origin of IGHG4 transcripts in IgG4-RD indicated that IgG4 could be switched from IgM directly, or from other IgG subclasses. In vitro assays with EBV-immortalized B cells, fibroblasts and epithelial cells revealed the effects of Th1-type and Th2-type responses on germinal center reaction, ectopic expression of MHC-II molecules, and formation of tertiary lymphoid structures. CONCLUSIONS: Synergistic effects of Th1- and Th2-type responses were involved in the pathogenesis of IgG4-RD via their influences on both acute inflammatory processes and the chronicity and complexity of IgG4-RD.

Involvement of MST1/mTORC1/STAT1 activity in the regulation of B-cell receptor signalling by chemokine receptor 2.

BACKGROUND: CCR2 is involved in maintaining immune homeostasis and regulating immune function. This study aims to elucidate the mechanism by which CCR2 regulates B-cell signalling. METHODS: In Ccr2-knockout mice, the development and differentiation of B cells, BCR proximal signals, actin movement and B-cell immune response were determined. Besides, the level of CCR2 in PBMC of SLE patients was analysed by bioinformatics. RESULTS: CCR2 deficiency reduces the proportion and number of follicular B cells, upregulates BCR proximal signalling and enhances the oxidative phosphorylation of B cells. Meanwhile, increased actin filaments aggregation and its associated early-activation events of B cells are also induced by CCR2 deficiency. The MST1/mTORC1/STAT1 axis in B cells is responsible for the regulation of actin remodelling, metabolic activities and transcriptional signalling, specific MST1, mTORC1 or STAT1 inhibitor can rescue the upregulated BCR signalling. Glomerular IgG deposition is obvious in CCR2-deficient mice, accompanied by increased anti-dsDNA IgG level. Additionally, the CCR2 expression in peripheral B cells of SLE patients is decreased than that of healthy controls. CONCLUSIONS: CCR2 can utilise MST1/mTORC1/STAT1 axis to regulate BCR signalling. The interaction between CCR2 and BCR may contribute to exploring the mechanism of autoimmune diseases.

Gut Microbiota Aberration in Patients of Systemic Sclerosis and Bleomycin-Induced Mice Model.

Systemic sclerosis (SSc) is an immune-mediated systemic autoimmune disease with unknown etiology, which has high morbidity and mortality. Current treatments to dispose of this disorder are limited. And there are still no ideal animal models that can fully replicate the four basic pathophysiological features of SSc, including vascular lesions, fibrosis, inflammation, and autoimmunity, let alone animal models specifically designed to study gastrointestinal lesions. It's essential to seek and establish appropriate animal models to explore the role of gut microbiota in the pathogenesis of SSc. In this study, we found similar gut microbiota aberration in patients of SSc and bleomycin (BLM)-induced mice model through 16S rRNA gene sequencing. In terms of phylum-level differences, the relative abundance of Bacteroidetes was significantly decreased and Firmicutes increased in the SSc patients and the mice. Notably, the genera of Lactobacillus, commonly used as a probiotic additive, was also elevated in SSc patients and BLM mice, which was consistent with a few of studies. Therefore, the model can likely mimic the pathological changes of gut microbiota in patients with SSc, which may offer an important potential platform for the in-depth understanding of gut microbiota aberration in patients with SSc and to devise potential disease-modifying treatments.

Multiple Processes May Involve in the IgG4-RD Pathogenesis: An Integrative Study via Proteomic and Transcriptomic Analysis.

Immunoglobulin G4-related disease (IgG4-RD) is a newly defined disease entity, while the exact pathogenesis is still not clear. Identifying the characters of IgG4-RD in proteomic and transcriptomic aspects will be critical to investigate the potential pathogenic mechanisms of IgG4-RD. We performed proteomic analysis realized with iTRAQ technique for serum samples from eight treatment-naive IgG4-RD patients and eight healthy volunteers, and tissue samples from two IgG4-RD patients and two non-IgG4-RD patients. Transcriptomic data (GSE40568 and GSE66465) was obtained from the GEO Dataset for validation. The weighted correlation network analysis (WGCNA) was applied to detect the gene modules correlated with IgG4-RD. KEGG pathway analysis was used to investigate pathways enriched in IgG4-RD samples. As a result, a total of 980 differentially expressed proteins (DEPs) in tissue and 94 DEPs in serum were identified between IgG4-RD and control groups. Three hundred fifty-four and two hundred forty-seven genes that most correlated with IgG4-RD were detected by WGCNA analysis in tissue and PBMC, respectively. We also found that DEPs in IgG4-RD samples were enriched in several immune-related activities including bacterial/viral infections and platelet activation as well as some immune related signaling pathways. In conclusion, we identified multiple processes/factors and several signaling pathways that may involve in the IgG4-RD pathogenesis, and found out some potential therapeutic targets for IgG4-RD.

The expression and regulation of HOX genes and membrane proteins among different cytogenetic groups of acute myeloid leukemia.

BACKGROUND: The cytogenetic aberrations were considered as markers for diagnosis and prognosis in acute myeloid leukemia (AML), while the expression and regulation under different cytogenetic groups remain to be fully elucidated. METHODS: In this paper, for favorable, poor, and cytogenetically normal groups of AML patients, we performed comprehensive bioinformatics analyses including identifying differentially expressed genes (DEGs) and microRNAs (miRNAs) among them, functional enrichment and regulatory networks. RESULTS: We found that DEGs were enriched in membrane-related processes. Eleven genes and two miRNAs were significantly differentially expressed among these three AML groups. In survival analysis, membrane-related genes and several miRNAs were significant on prognostic outcome. Notably, six HOXA and three HOXB genes were significantly in low expression and high methylation in AML with favorable cytogenetics. Meanwhile, the miRNA-HOX gene co-regulatory networks revealed that HOXA5 was a hub node and regulated an AML oncogene SPARC. CONCLUSION: Our work may provide novel insights to the molecular characteristics and classification between AML with different cytogenetics.

A 6-Membrane Protein Gene score for prognostic prediction of cytogenetically normal acute myeloid leukemia in multiple cohorts.

Background: Cytogenetically normal acute myeloid leukemia (CN-AML) is a large proportion of AMLs with diverse prognostic outcomes. Identifying membrane protein genes as prognostic factors to stratify CN-AML patients will be critical to improve their outcomes. Purpose: This study aims to identify prognostic factors to stratify CN-AML patients to choose better treatments and improve their outcomes. Methods: CN-AML data were from TCGA cohort (n = 79) and four GEO datasets. We identified independent prognostic genes by Cox regression and Kaplan-Meier methods, and constructed linear regression model using LASSO algorithm. The prediction error curve was calculated using R package "pec". Results: Based on independent prognostic membrane genes, we constructed a regression model for CN-AML prognosis prediction: score = (0.0492 * CD52) - (0.0018 * CD96) + (0.0131 * EMP1) + (0.2058 * TSPAN2) + (0.0234 * STAB1) - (0.3658 * MBTPS1), which was named as MPG6 (6-Membrane Protein Gene) score. Tested in multiple CN-AML datasets, consistent results showed that CN-AML patients with high MPG6 score had poor survival, higher WBC count and shorter EFS. Comparing with other reported scoring models, the benchmark result of MPG6 achieved better association with survival in multiple cohorts. Moreover, by combining with other clinical indicators in CN-AML, MPG6 could improve the performance of survival prediction and serve as a robust prognostic factor. Conclusions: We identified the MPG6 score as a stable indicator with great potential for clinical application in risk stratification and outcome prediction in CN-AML.

Landscape of cancer diagnostic biomarkers from specifically expressed genes.

Although there has been great progress in cancer treatment, cancer remains a serious health threat to humans because of the lack of biomarkers for diagnosis, especially for early-stage diagnosis. In this study, we comprehensively surveyed the specifically expressed genes (SEGs) using the SEGtool based on the big data of gene expression from the The Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression (GTEx) projects. In 15 solid tumors, we identified 233 cancer-specific SEGs (cSEGs), which were specifically expressed in only one cancer and showed great potential to be diagnostic biomarkers. Among them, three cSEGs (OGDH, MUDENG and ACO2) had a sample frequency >80% in kidney cancer, suggesting their high sensitivity. Furthermore, we identified 254 cSEGs as early-stage diagnostic biomarkers across 17 cancers. A two-gene combination strategy was applied to improve the sensitivity of diagnostic biomarkers, and hundreds of two-gene combinations were identified with high frequency. We also observed that 13 SEGs were targets of various drugs and nearly half of these drugs may be repurposed to treat cancers with SEGs as their targets. Several SEGs were regulated by specific transcription factors in the corresponding cancer, and 39 cSEGs were prognosis-related genes in 7 cancers. This work provides a survey of cancer biomarkers for diagnosis and early diagnosis and new insights to drug repurposing. These biomarkers may have great potential in cancer research and application.

Identification of STAB1 in Multiple Datasets as a Prognostic Factor for Cytogenetically Normal AML: Mechanism and Drug Indications.

Cytogenetically normal acute myeloid leukemia (CN-AML) presents with diverse outcomes in different patients and is categorized as an intermediate prognosis group. It is important to identify prognostic factors for CN-AML risk stratification. In this study, using the TCGA CN-AML dataset, we found that the scavenger receptor stabilin-1 (STAB1) is a prognostic factor for poor outcomes and validated it in three other independent CN-AML datasets. The high STAB1 expression (STAB1high) group had shorter event-free survival compared with the low STAB1 expression (STAB1low) group in both the TCGA dataset (n = 79; p = 0.0478) and GEO: GSE6891 dataset (n = 187; p = 0.0354). Differential expression analysis between the STAB1high and STAB1low groups revealed that upregulated genes in the STAB1high group were enriched in pathways related to cell adhesion and migration and immune responses. We confirmed that STAB1 suppression inhibits cell growth in KG1a and NB4 leukemia cells. Expression correlation analyses between STAB1 and cancer drug targets suggested that patients in the STAB1low group are more sensitive to the BCL2 inhibitor venetoclax, and we confirmed it in cell lines. In conclusion, we identified STAB1 as a prognostic factor for CN-AML in multiple datasets, explored its underlying mechanism, and provided potential therapeutic indications.

Similar Transition Processes in Synovial Fibroblasts from Rheumatoid Arthritis and Osteoarthritis: A Single-Cell Study.

Rheumatoid arthritis (RA) and osteoarthritis (OA) are common rheumatic disorders that primarily involve joints. The inflammation of the synovium can be observed in both of the two diseases. Synovial fibroblasts (SFs) play an important role in the inflammatory process of the synovium. The functional states of synovial fibroblasts are heterogeneous, and the detailed transition process of their functional states is still unclear. By using transcriptomic data of SFs at a single-cell level, we found a similar transition process for SFs in RA and OA. We also identified the potential regulatory effects of the WNT signaling pathway, the TGF-ß signaling pathway, the FcεRI signaling pathway, and the ERBB signaling pathway on modifying the SFs' functional state. These findings indicate potentially overlapped pathogenic mechanisms in these two diseases, which may help uncover new therapeutic targets to ameliorate disease progression.

High IL2RA mRNA expression is an independent adverse prognostic biomarker in core binding factor and intermediate-risk acute myeloid leukemia.

BACKGROUND: Elevated protein expressions of CD markers such as IL2RA/CD25, CXCR4/CD184, CD34 and CD56 are associated with adverse prognosis in acute myeloid leukemia (AML). However, the prognostic value of mRNA expressions of these CD markers in AML remains unclear. Through our pilot evaluation, IL2RA mRNA expression appeared to be the best candidate as a prognostic biomarker. Therefore, the aim of this study is to characterize the prognostic value of IL2RA mRNA expression and evaluate its potential to refine prognostification in AML. METHODS: In a cohort of 239 newly diagnosed AML patients, IL2RA mRNA expression were measured by TaqMan realtime quantitative PCR. Morphological, cytogenetics and mutational analyses were also performed. In an intermediate-risk AML cohort with 66 patients, the mRNA expression of prognostic biomarkers (BAALC, CDKN1B, ERG, MECOM/EVI1, FLT3, ID1, IL2RA, MN1 and WT1) were quantified by NanoString technology. A TCGA cohort was analyzed to validate the prognostic value of IL2RA. For statistical analysis, Mann-Whitney U test, Fisher exact test, logistic regression, Kaplan-Meier and Cox regression analyses were used. RESULTS: In AML cohort of 239 patients, high IL2RA mRNA expression independently predicted shorter relapse free survival (RFS, p < 0.001) and overall survival (OS, p < 0.001) irrespective of age, cytogenetics, FLT3-ITD or c-KIT D816V mutational status. In core binding factor (CBF) AML, high IL2RA mRNA expression correlated with FLT3-ITD status (p = 0.023). Multivariable analyses revealed that high IL2RA expression (p = 0.002), along with c-KIT D816V status (p = 0.013) significantly predicted shorter RFS, whereas only high IL2RA mRNA expression (p = 0.014) significantly predicted shorter OS in CBF AML. In intermediate-risk AML in which multiple gene expression markers were tested by NanoString, IL2RA significantly correlated with ID1 (p = 0.006), FLT3 (p = 0.007), CDKN1B (p = 0.033) and ERG (p = 0.030) expressions. IL2RA (p < 0.001) and FLT3 (p = 0.008) expressions remained significant in predicting shorter RFS, whereas ERG (p = 0.008) and IL2RA (p = 0.044) remained significant in predicting shorter OS. Similar analyses in TCGA intermediate-risk AML showed the independent prognostic role of IL2RA in predicting event free survival (p < 0.001) and OS (p < 0.001). CONCLUSIONS: High IL2RA mRNA expression is an independent and adverse prognostic factor in AML and specifically stratifies patients to worse prognosis in both CBF and intermediate-risk AML.

Enhanced Response of Acute Monocytic Leukemia Cells to Low-dose Cytarabine by 1,25-dihydroxyvitamin D3.

Low-dose cytarabine combined with differentiating or DNA hypomethylating agents, such as vitamin D compounds, is a potential regimen to treat acute myeloid leukemia (AML) patients who are unfit for high-intensity chemotherapy. The present study aimed to determine which subset of AML would be most responsive to low-dose cytarabine with the differentiating agent 1,25-dihydroxyvitamin D3 (1,25-D3). Here, firstly, cBioPortal database was used and we found out that vitamin D receptor (VDR) was highly expressed in acute monocytic leukemia (M5) and high VDR expression was associated with a poor survival of AML patients. Then, we confirmed that 1,25-D3 at clinical available concentration could induce more significant differentiation in acute monocytic leukemia cell lines (U937, MOLM-13, THP-1) and blasts from M5 patients than in non-monocytic cell lines (KGla and K562) and blasts from M2 patient. Finally, it was shown that the combination of 1,25-D3 and low-dose cytarabine further increased the differentiating rate, growth inhibition and G0/G1 arrest, while mild changes were found in the apoptosis in acute monocytic leukemia cell lines. Our study demonstrates that the enhanced response of acute monocytic leukemia cells to low-dose cytarabine by 1,25-D3 might indicate a novel therapeutic direction for patients with acute monocytic leukemia, especially for elderly and frail ones.

SEGtool: a specifically expressed gene detection tool and applications in human tissue and single-cell sequencing data.

Different tissues and diseases have distinct transcriptional profilings with specifically expressed genes (SEGs). So, the identification of SEGs is an important issue in the studies of gene function, biological development, disease mechanism and biomarker discovery. However, few accurate and easy-to-use tools are available for RNA sequencing (RNA-seq) data to detect SEGs. Here, we presented SEGtool, a tool based on fuzzy c-means, Jaccard index and greedy annealing method for SEG detection automatically and self-adaptively ignoring data distribution. Testing result showed that our SEGtool outperforms the existing tools, which was mainly developed for microarray data. By applying SEGtool to Genotype-Tissue Expression (GTEx) human tissue data set, we detected 3181 SEGs with tissue-related functions. Regulatory networks reveal tissue-specific transcription factors regulating many SEGs, such as ETV2 in testis, HNF4A in liver and NEUROD1 in brain. Applied to a case study of single-cell sequencing (SCS) data from embryo cells, we identified many SEGs in specific stages of human embryogenesis. Notably, SEGtool is suitable for RNA-seq data and even SCS data with high specificity and accuracy. An implementation of SEGtool R package is freely available at http://bioinfo.life.hust.edu.cn/SEGtool/.

Gene expression, regulation of DEN and HBx induced HCC mice models and comparisons of tumor, para-tumor and normal tissues.

BACKGROUND: Hepatocellular carcinoma (HCC) is the leading cause of cancer mortality. Chemical and virus induction are two major risk factors, however, the potential molecular mechanisms of their differences remain elusive. In this study, to identify the similarities and differences between chemical and virus induced HCC models, we compared the gene expression profiles between DEN and HBx mice models, as well as the differences among tumor, para-tumor and normal tissues. METHODS: We sequenced both gene and microRNA (miRNA) expression for HCC tumor tissues, para-tumor and normal liver tissues from DEN model mice (30-week-old) and downloaded the corresponding microarray expression data of HBx model from GEO database. Then differentially expressed genes (DEGs), miRNAs and transcription factors (TFs) were detected by R packages and performed functional enrichment analysis. To explore the gene regulatory network in HCC models, miRNA and TF regulatory networks were constructed by target prediction. RESULTS: For model comparison, although DEGs between tumor and normal tissues in DEN and HBx models only had a small part of overlapping, they shared common pathways including lipid metabolism, oxidation-reduction process and immune process. For tissue comparisons in each model, genes in oxidation-reduction process were down-regulated in tumor tissues and genes in inflammatory response showed the highest expression level in para-tumor tissues. Genes highly expressed in both tumor and para-tumor tissues in two models mainly participated in immune and inflammatory response. Genes expressed in HBx model were also involved in cell proliferation and cell migration etc. Network analysis revealed that several miRNAs such as miR-381-3p, miR-142a-3p, miR-214-3p and TFs such as Egr1, Atf3 and Klf4 were the core regulators in HCC. CONCLUSIONS: Through the comparative analyses, we found that para-tumor tissue is a highly inflammatory tissue while the tumor tissue is specific with both inflammatory and cancer signaling pathways. The DEN and HBx mice models have different gene expression pattern but shared pathways. This work will help to elucidate the molecular mechanisms underlying different HCC models.

Interferon-γ alters the immune-related miRNA expression of microvesicles derived from mesenchymal stem cells.

Increasing studies have demonstrated that interferon gamma (IFN-γ), which serves as a critical inflammatory cytokine, is essential to induce the immunosuppressive effects of mesenchymal stem cells (MSCs). However, the mechanisms underlying the enhanced immunosuppressive effects of IFN-γ-stimulated MSCs (γMSCs) are not fully understood. MSC-derived microvesicles (MSC-MVs) have been viewed as potential pivotal mediators of the immunosuppressive effects of MSCs. Moreover, microRNAs (miRNAs) are important regulators of immunological processes and can be shuttled from cell to cell by MVs. The aim of our study was to analyze the the miRNA expression signature of MVs derived from γMSCs (γMSC-MVs), which may provide better understanding of the immunosuppressive property of their parent cells. Through miRNA microarray and bioinformatics analysis, we found 62 significantly differentially expressed miRNAs (DEMs) in γMSC-MVs compared with MSC-MVs. And the potential target genes and signaling pathways regulated by DEMs were predicted and analyzed. Interestingly, many DEMs and predicted signaling pathways had been demonstrated to be involved in immunoregulation. Furthermore, the network between immunoregulation-related pathways and relevant DEMs was constructed. Collectively, our research on the miRNA repertoires of γMSC-MVs not only provides new perspectives into the mechanisms underlying the enhanced immunosuppressive property of γMSCs, but also paves the way to clinical application of these potent organelles in the future.