Identification of co-expression network correlated with different periods of adipogenic and osteogenic differentiation of BMSCs by weighted gene co-expression network analysis (WGCNA).
Ontology highlight
ABSTRACT: Background
The differentiation of bone marrow mesenchymal stem cells is a complex and dynamic process. The gene expression pattern and mechanism of different periods of adipogenic and osteogenic differentiation remain unclear. Additionally, the interaction between these two lineage determination requires further exploration.Results
Five modules that were most significantly associated with osteogenic or adipogenic differentiation of BMSCs were selected for further investigation. Biological terms (e.g. ribosome biogenesis, TNF-α signalling pathway, glucose import and fatty acid metabolism) along with hub transcription factors (e.g. PPARG and YY1) and hub miRNAs (e.g. hsa-mir-26b-5p) were enriched in different modules. The expression pattern of 6 hub genes, ADIPOQ, FABP4, SLC7A5, SELPLG, BIRC3, and KLHL30 was validated by RT-qPCR. Finally, cell staining experiments extended the findings of bioinformatics analysis.Conclusion
This study identified the key genes, biological functions, and regulators of each time point of adipogenic and osteogenic differentiation of BMSCs and provided novel evidence and ideas for further research on the differentiation of BMSCs.
SUBMITTER: Liu Y
PROVIDER: S-EPMC8035768 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Osteosarcoma (OS), the most common malignant bone tumor, accounts for the heavy healthy threat in the period of children and adolescents. OS occurrence usually correlates with early metastasis and high death rate. This study aimed to better understand the mechanism of OS metastasis.Based on Gene Expression Omnibus (GEO) database, we downloaded 4 expression profile data sets associated with OS metastasis, and selected differential expressed genes. Weighted gene co-expression network analysis (WGCNA) approach allowed us to investigate the most OS metastasis-correlated module. Gene Ontology functional and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were used to give annotation of selected OS metastasis-associated genes.We select 897 differential expressed genes from OS metastasis and OS non-metastasis groups. Based on these selected genes, WGCNA further explored 142 genes included in the most OS metastasis-correlated module. Gene Ontology functional and KEGG pathway enrichment analyses showed that significantly OS metastasis-associated genes were involved in pathway correlated with insulin-like growth factor binding.Our research figured out several potential molecules participating in metastasis process and factors acting as biomarker. With this study, we could better explore the mechanism of OS metastasis and further discover more therapy targets.
| S-EPMC6023727 | biostudies-literature
Project description:Osteoporosis, or steroid-induced osteonecrosis of the hip, is accompanied by increased bone marrow adipogenesis. Such a disorder of adipogenic/osteogenic differentiation, affecting bone-marrow-derived mesenchymal stem cells (BMSCs), contributes to bone loss during aging. Here, we investigated the effects of extracellular vesicles (EVs) isolated from human (h)BMSCs during different stages of osteogenic differentiation on the osteogenic and adipogenic differentiation capacity of naïve (undifferentiated) hBMSCs. We observed that all EV groups increased viability and proliferation capacity and suppressed the apoptosis of naïve hBMSCs. In particular, EVs derived from hBMSCs at late-stage osteogenic differentiation promoted the osteogenic potential of naïve hBMSCs more effectively than EVs derived from naïve hBMSCs (naïve EVs), as indicated by the increased gene expression of COL1A1 and OPN. In contrast, the adipogenic differentiation capacity of naïve hBMSCs was inhibited by treatment with EVs from osteogenic differentiated hBMSCs. Proteomic analysis revealed that osteogenic EVs and naïve EVs contained distinct protein profiles, with pro-osteogenic and anti-adipogenic proteins encapsulated in osteogenic EVs. We speculate that osteogenic EVs could serve as an intercellular communication system between bone- and bone-marrow adipose tissue, for transporting osteogenic factors and thus favoring pro-osteogenic processes. Our data may support the theory of an endocrine circuit with the skeleton functioning as a ductless gland.
| S-EPMC9406723 | biostudies-literature
Project description:Weighted gene coexpression network analysis (WGCNA) is a novel approach that can quickly analyze the relationships between genes and traits. In this study, the milk yield, lactose, fat, and protein of Holstein dairy cows were detected in a lactation cycle. Meanwhile, a total of 18 gene expression profiles were detected using mammary glands from six lactation stages (day 7 to calving, -7 d; day 30 post-calving, 30 d; day 90 post-calving, 90 d; day 180 post-calving, 180 d; day 270 post-calving, 270 d; day 315 post-calving, 315 d). On the basis of the 18 profiles, WGCNA identified for the first time 10 significant modules that may be related to lactation stage, milk yield, and the main milk composition content. Genes in the 10 significant modules were examined with gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. The results revealed that the galactose metabolism pathway was a potential candidate for milk yield and milk lactose synthesis. In -7 d, ion transportation was more frequent and cell proliferation related terms became active. In late lactation, the suppressor of cytokine signaling 3 (SOCS3) might play a role in apoptosis. The sphingolipid signaling pathway was a potential candidate for milk fat synthesis. Dairy cows at 315 d were in a period of cell proliferation. Another notable phenomenon was that nonlactating dairy cows had a more regular circadian rhythm after a cycle of lactation. The results provide an important theoretical basis for the further molecular breeding of dairy cows.
| S-EPMC7911360 | biostudies-literature
Project description:Weighted gene co-expression network analysis (WGCNA) is used to detect clusters with highly correlated genes. Measurements of correlation most typically rely on linear relationships. However, a linear relationship does not always model pairwise functional-related dependence between genes. In this paper, we first compared 6 different correlation methods in their ability to capture complex dependence between genes in three different tissues. Next, we compared their gene-pairwise coefficient results and corresponding WGCNA results. Finally, we applied a recently proposed correlation method, Hellinger correlation, as a more sensitive correlation measurement in WGCNA. To test this method, we constructed gene networks containing co-expression gene modules from RNA-seq data of human frontal cortex from Alzheimer's disease patients. To test the generality, we also used a microarray data set from human frontal cortex, single cell RNA-seq data from human prefrontal cortex, RNA-seq data from human temporal cortex, and GTEx data from heart. The Hellinger correlation method captures essentially similar results as other linear correlations in WGCNA, but provides additional new functional relationships as exemplified by uncovering a link between inflammation and mitochondria function. We validated the network constructed with the microarray and single cell sequencing data sets and a RNA-seq dataset of temporal cortex. We observed that this new correlation method enables the detection of non-linear biologically meaningful relationships among genes robustly and provides a complementary new approach to WGCNA. Thus, the application of Hellinger correlation to WGCNA provides a more flexible correlation approach to modelling networks in gene expression analysis that uncovers novel network relationships.
| S-EPMC9307959 | biostudies-literature
Project description:Long noncoding RNAs (lncRNAs) are critical gene expression regulators and are involved in several bone diseases. To explore the potential roles of lncRNAs in osteonecrosis of the femoral head (ONFH), we investigated for the first time the lncRNA expression profile of bone marrow mesenchymal stem cells (BMSCs) from patients with steroid-induced ONFH (SONFH) with microarray and bioinformatics analysis. A total of 1878 lncRNAs and 838 mRNAs were significantly up-regulated while 1842 lncRNAs and 1937 mRNAs were statistically down-regulated in the SONFH group compared with control group. The results validated by qRT-PCR were consistent with the microarray profiling data, especially involved in upregulation and downregulation of critical lncRNAs as well as mRNAs expression related to adipogenic and osteogenic differentiation. Pathway analyses revealed 40 signaling pathways with significant differences, especially the signaling pathways to regulate stem cell pluripotency. The CNC and ceRNA network indicated that lncRNA RP1-193H18.2, MALAT1 and HOTAIR were associated with abnormal osteogenic and adipogenic differentiation of BMSCs in the patients with SONFH. Our results suggest the lncRNA expression profiles were closely associated with the abnormal adipogenic and osteogenic transdifferentiation of BMSCs during the development of SONFH and explore a new insight into the molecular mechanisms of SONFH.
| S-EPMC6002551 | biostudies-literature
Project description:Background: Long non-coding RNAs (lncRNAs), which are over 200 nt in length, have a key role in tumorigenesis and disease progression. To explore the role of prognostic lncRNAs in adult acute myeloid leukemia (AML), the expression profiles of lncRNAs and mRNAs in AML were analyzed. Methods: The RNAseq data of 167 adult AML patients and the corresponding clinical information were downloaded from The Cancer Genome Atlas (TCGA), which is a publicly available database. The RPKM values of the RNAseq data were subjected to weighted gene co-expression network analysis (WGCNA) in modularization. Results: We identified survival specific lncRNAs and mRNAs, which were divided into modules by coexpression analysis. The lncRNAs were mainly annotated into "Fc gamma R-mediated phagocytosis". The hub lncRNA and co-expressed mRNAs were further selected for analysis of risk stratification. LncRNA-LOC646762 may contribute to AML through the "endocytosis" signaling pathway. Finally, the expression levels of LOC646762 and co-expressed CCND3, CBR1, C10orf54, CD97 and BLOC1S1 in the adult AML patients and healthy volunteers were validated by qRT-PCR, and then their roles in prognosis and risk stratification were identified. Conclusions: Prognostic lncRNA-LOC646762, which may contribute to AML through the "endocytosis" signaling pathway, may act as a biomarker for predicting the survival of adult AML patients, as well as for risk stratification.
| S-EPMC6746144 | biostudies-literature
Project description:The present study aimed to identify the recurrence‑associated genes in colon cancer, which may provide theoretical evidence for the development of novel methods to prevent tumor recurrence. Colon cancer and matched normal samples microarray data (E‑GEOD‑39582) were downloaded from ArrayExpress. Genes with significant variation were identified, followed by the screening of differentially expressed genes (DEGs). Subsequently, the co‑expression network of DEGs was constructed using the weighted correlation network analysis (WGCNA) method, which was verified using the validation dataset. The significant modules associated with recurrence in the network were subsequently screened and verified in another independent dataset E‑GEOD‑33113. Function and pathway enrichment analyses were also conducted to determine the roles of selected genes. Survival analysis was performed to identify the association between these genes and survival. A total of 434 DEGs were identified in the colon samples, and stress‑associated endoplasmic reticulum protein family member 2 (SERP2) and long non‑coding RNA‑0219 (LINC0219) were determined to be the vital DEGs between all the three sub‑type groups with different clinical features. The brown module was identified to be the most significant module in the co‑expression network associated with the recurrence of colon cancer, which was verified in the E‑GEOD‑33113 dataset. Top 10 genes in the brown module, including EGF containing fibulin like extracellular matrix protein 2 (EFEMP2), fibrillin 1 (FBN1) and secreted protein acidic and cysteine rich (SPARC) were also associated with survival time of colon cancer patients. Further analysis revealed that the function of cell adhesion, biological adhesion, extracellular matrix (ECM) organization, pathways of ECM‑receptor interaction and focal adhesion were the significantly changed terms in colon cancer. In conclusion, SERP2, EFEMP2, FBN1, SPARC, and LINC0219 were revealed to be the recurrence‑associated molecular and prognostic indicators in colon cancer by WGCNA co‑expression network analysis.
| S-EPMC5865817 | biostudies-literature
Project description:Our previous studies shown that syndecan-1 (SDC1) may be a novel class of biomarkers for the diagnosis and treatment of glioma, but its specific roles and the in-depth molecular mechanism remain elusive. Here, we used Estimation of STromal and Immune cells in Malignant Tumor tissues using Expression data (ESTIMATE) algorithms and single-sample Gene Set Enrichment Analysis (ssGSEA) algorithms to evaluate the immune score of tumor samples and quantify the relative infiltration of immune cells in the tumor microenvironment (TME), respectively, in different data sets obtained from the Chinese Glioma Genome Atlas and The Cancer Gene Atlas. Next, we calculate the correlation of the immune score and immune cells with SDC1, respectively. To identify the specific process regulated by SDC1, the differentially expressed genes (DEGs) analysis between the high and low expression of SDC1 of glioma samples were used to discover the hub genes through Weighted Gene Coexpression Network Analysis (WGCNA). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed cardinal biological processes and pathways involved in genes and tumor grade correlation and survival analysis verified its significance in glioma. The results show that SDC1 is associated with the immune infiltration of glioma in the TME, especially activated CD4+T cells and CD8+T cells. The three data sets filter 8,887 DEGs, the genes in the blue modules were selected as hub genes in WGCNA. GO and KEGG analysis found eight genes in the blue modules involved in antigen processing and presentation in T cells in glioma. Kaplan-Meier estimator and log-rank test statistic determined that the introduced genes are associated with poor prognosis in glioma. Protein-protein network interaction analysis showed that SDC1 may regulate antigen processing and presentation through CTSL or CD4 in glioma. Finally, this study provided insights and clues for the next research direction of SDC1 and identified the key pathways and genes that might participate in the immune escape of glioma. These results might provide a new insight on the study of immune infiltration of glioma in the future.
| S-EPMC9165731 | biostudies-literature
Project description:BACKGROUND: Advances in "omics" technologies have revolutionized the collection of biological data. A matching revolution in our understanding of biological systems, however, will only be realized when similar advances are made in informatic analysis of the resulting "big data." Here, we compare the capabilities of three conventional and novel statistical approaches to summarize and decipher the tomato metabolome. METHODOLOGY: Principal component analysis (PCA), batch learning self-organizing maps (BL-SOM) and weighted gene co-expression network analysis (WGCNA) were applied to a multivariate NMR dataset collected from developmentally staged tomato fruits belonging to several genotypes. While PCA and BL-SOM are appropriate and commonly used methods, WGCNA holds several advantages in the analysis of highly multivariate, complex data. CONCLUSIONS: PCA separated the two major genetic backgrounds (AC and NC), but provided little further information. Both BL-SOM and WGCNA clustered metabolites by expression, but WGCNA additionally defined "modules" of co-expressed metabolites explicitly and provided additional network statistics that described the systems properties of the tomato metabolic network. Our first application of WGCNA to tomato metabolomics data identified three major modules of metabolites that were associated with ripening-related traits and genetic background.
| S-EPMC3198806 | biostudies-literature
Project description:Trichomes are attractive cells for terpenoid biosynthesis and accumulation in Artemisia annua. However, the molecular process underlying the trichome of A. annua is not yet fully elucidated. In this study, an analysis of multi-tissue transcriptome data was performed to examine trichome-specific expression patterns. A total of 6646 genes were screened and highly expressed in trichomes, including artemisinin biosynthetic genes such as amorpha-4,11-diene synthase (ADS) and cytochrome P450 monooxygenase (CYP71AV1). Mapman and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed that trichome-specific genes were mainly enriched in lipid metabolism and terpenoid metabolism. These trichome-specific genes were analyzed by a weighted gene co-expression network analysis (WGCNA), and the blue module linked to terpenoid backbone biosynthesis was determined. Hub genes correlated with the artemisinin biosynthetic genes were selected based on TOM value. ORA, Benzoate carboxyl methyltransferase (BAMT), Lysine histidine transporter-like 8 (AATL1), Ubiquitin-like protease 1 (Ulp1) and TUBBY were revealed as key hub genes induced by methyl jasmonate (MeJA) for regulating artemisinin biosynthesis. In summary, the identified trichome-specific genes, modules, pathways and hub genes provide clues and shed light on the potential regulatory mechanisms of artemisinin biosynthesis in trichomes in A. annua.
| S-EPMC10217836 | biostudies-literature