Osteogenesis depends on commissioning of a network of stem cell transcription factors that act as repressors of adipogenesis

In vitro determination of osteo-adipogenic lineage choice of bone marrow stromal/stem cells (BMSCs)

Bone marrow stromal/stem cells (BMSCs) are primitive and heterogeneous cells that can be differentiated into osteoblasts, adipocytes and other subsets. Their bone-fat lineage commitment is responsible for the homeostasis of bone marrow microenvironment. However, there are little effective methods and evidence to simultaneously visualise the lineage commitment of BMSCs. Here we provide a bivalent differentiation medium that can enable BMSCs differentiation into osteoblasts and adipocytes in vitro, and establish a method to simultaneously distinguish osteoblasts or adipocytes from the heterogeneous BMSCs based on Alizarin red S and Oil red O staining, which have been used for detection of specific mineralized nodules and lipid droplets, respectively. This assay provides a specifically simple but effective and low-cost method to evaluate the efficiency of osteo-adipogenic (OA) allocation of BMSCs.►Researchers can utilize the bivalent differentiation medium to evaluate the efficiency of osteogenic and adipogenic differentiation of BMSCs in vitro.

The E3 ubiquitin-protein ligase UHRF1 promotes adipogenesis and limits fibrosis by suppressing GPNMB-mediated TGF-β signaling

The E3 ubiquitin-ligase UHRF1 is an epigenetic regulator coordinating DNA methylation and histone modifications. However, little is known about how it regulates adipogenesis or metabolism. In this study, we discovered that UHRF1 is a key regulatory factor for adipogenesis, and we identified the altered molecular pathways that UHRF1 targets. Using CRISPR/Cas9-based knockout strategies, we discovered the whole transcriptomic changes upon UHRF1 deletion. Bioinformatics analyses revealed that key adipogenesis regulators such PPAR-γ and C/EBP-α were suppressed, whereas TGF-β signaling and fibrosis markers were upregulated in UHRF1-depleted differentiating adipocytes. Furthermore, UHRF1-depleted cells showed upregulated expression and secretion of TGF-β1, as well as the glycoprotein GPNMB. Treating differentiating preadipocytes with recombinant GPNMB led to an increase in TGF-β protein and secretion levels, which was accompanied by an increase in secretion of fibrosis markers such as MMP13 and a reduction in adipogenic conversion potential. Conversely, UHRF1 overexpression studies in human cells demonstrated downregulated levels of GPNMB and TGF-β, and enhanced adipogenic potential. In conclusion, our data show that UHRF1 positively regulates 3T3-L1 adipogenesis and limits fibrosis by suppressing GPNMB and TGF-β signaling cascade, highlighting the potential relevance of UHRF1 and its targets to the clinical management of obesity and linked metabolic disorders.

Semaglutide 6-months therapy of type 2 diabetes mellitus restores adipose progenitors potential to develop metabolically active adipocytes

BACKGROUND

Nowadays type 2 diabetes mellitus (T2DM) leads to population mortality growth. Today glucagon-like peptide type 1 receptor agonists (GLP-1 RA) are one of the most promising glucose-lowered drugs with anorexigenic and cardioprotective effects. The present study aims to determine the effects of GLP-1 RA semaglutide 6-month therapy on T2DM patient metabolic parameters and adipose progenitor cell health.

METHODS

T2DM patients (N = 8) underwent clinical characterization and subcutaneous fat biopsy at start point and after semaglutide 6-month therapy. Adipose-derived stem cells (ADSC) were isolated by enzymatic method. Cell proliferation analysis was performed by MTT and immunocytochemistry. White and beige adipogenesis was analyzed by BODIPY493/503 staining and confocal microscopy. Adipocyte's metabolic properties were estimated by 3H- and 14C-based metabolic assays. Thermogenesis analysis was performed by ERthermAC staining and confocal microscopy. Protein markers were assessed by Western blotting.

RESULTS

Semaglutide 6-month therapy demonstrated significant anorexigenic and glucose-lowering effects. However, insulin sensitivity (HOMA-IR and M-index) was unchanged after therapy. Semaglutide 6-month therapy increased ADSC proliferation and white and beige adipogenesis. Moreover, lipid droplets fragmentation was observed in beige adipocytes. Both white and beige adipocytes after semaglutide therapy demonstrated 2-3 fold growth of glucose uptake without changes in insulin sensitivity. Newly formed white adipocytes demonstrated glucose utilization for active ATP synthesis, whereas beige adipocytes for canonical thermogenesis.

CONCLUSIONS

Our study has revealed that semaglutide 6-month therapy has not only systemic anorexigenic effects, but can markedly improve adipose tissue health. We have demonstrated critical restoration of ADSC renewal functions, which potentially can be involved in semaglutide based weight loss.

Xenopus tropicalis osteoblast-specific open chromatin regions reveal promoters and enhancers involved in human skeletal phenotypes and shed light on early vertebrate evolution

While understanding the genetic underpinnings of osteogenesis has far-reaching implications for skeletal diseases and evolution, a comprehensive characterization of the osteoblastic regulatory landscape in non-mammalian vertebrates is still lacking. Here, we compared the ATAC-Seq profile of Xenopus tropicalis (Xt) osteoblasts to a variety of non mineralizing control tissues, and identified osteoblast-specific nucleosome free regions (NFRs) at 527 promoters and 6747 distal regions. Sequence analyses, Gene Ontology, RNA-Seq and ChIP-Seq against four key histone marks confirmed that the distal regions correspond to bona fide osteogenic transcriptional enhancers exhibiting a shared regulatory logic with mammals. We report 425 regulatory regions conserved with human and globally associated to skeletogenic genes. Of these, 35 regions have been shown to impact human skeletal phenotypes by GWAS, including one trps1 enhancer and the runx2 promoter, two genes which are respectively involved in trichorhinophalangeal syndrome type I and cleidocranial dysplasia. Intriguingly, 60 osteoblastic NFRs also align to the genome of the elephant shark, a species lacking osteoblasts and bone tissue. To tackle this paradox, we chose to focus on dlx5 because its conserved promoter, known to integrate regulatory inputs during mammalian osteogenesis, harbours an osteoblast-specific NFR in both frog and human. Hence, we show that dlx5 is expressed in Xt and elephant shark odontoblasts, supporting a common cellular and genetic origin of bone and dentine. Taken together, our work (i) unravels the Xt osteogenic regulatory landscape, (ii) illustrates how cross-species comparisons harvest data relevant to human biology and (iii) reveals that a set of genes including bnc2, dlx5, ebf3, mir199a, nfia, runx2 and zfhx4 drove the development of a primitive form of mineralized skeletal tissue deep in the vertebrate lineage.

Cytoskeletal rearrangement precedes nucleolar remodeling during adipogenesis

Differentiation of adipose progenitor cells into mature adipocytes entails a dramatic reorganization of the cellular architecture to accommodate lipid storage into cytoplasmic lipid droplets. Lipid droplets occupy most of the adipocyte volume, compressing the nucleus beneath the plasma membrane. How this cellular remodeling affects sub-nuclear structure, including size and number of nucleoli, remains unclear. We describe the morphological remodeling of the nucleus and the nucleolus during in vitro adipogenic differentiation of primary human adipose stem cells. We find that cell cycle arrest elicits a remodeling of nucleolar structure which correlates with a decrease in protein synthesis. Strikingly, triggering cytoskeletal rearrangements mimics the nucleolar remodeling observed during adipogenesis. Our results point to nucleolar remodeling as an active, mechano-regulated mechanism during adipogenic differentiation and demonstrate a key role of the actin cytoskeleton in defining nuclear and nucleolar architecture in differentiating human adipose stem cells.

SOD2 promotes the immunosuppressive function of mesenchymal stem cells at the expense of adipocyte differentiation

The potent immunomodulatory function of mesenchymal stem/stromal cells (MSCs) elicited by proinflammatory cytokines IFN-γ and TNF-α (IT) is critical to resolve inflammation and promote tissue repair. However, little is known about how the immunomodulatory capability of MSCs is related to their differentiation competency in the inflammatory microenvironment. In this study, we demonstrate that the adipocyte differentiation and immunomodulatory function of human adipose tissue-derived MSCs (MSC(AD)s) are mutually exclusive. Mitochondrial reactive oxygen species (mtROS), which promote adipocyte differentiation, were decreased in MSC(AD)s due to IT-induced upregulation of superoxide dismutase 2 (SOD2). Furthermore, knockdown of SOD2 led to enhanced adipogenic differentiation but reduced immunosuppression capability of MSC(AD)s. Interestingly, the adipogenic differentiation was associated with increased mitochondrial biogenesis and upregulation of peroxisome proliferator-activated receptor gamma coactivator-1 alpha (PPARGC1A/PGC-1α) expression. IT inhibited PGC-1α expression and decreased mitochondrial mass but promoted glycolysis in an SOD2-dependent manner. MSC(AD)s lacking SOD2 were compromised in their therapeutic efficacy in DSS-induced colitis in mice. Taken together, these findings indicate that the adipogenic differentiation and immunomodulation of MSC(AD)s may compete for resources in fulfilling the respective biosynthetic needs. Blocking of adipogenic differentiation by mitochondrial antioxidant may represent a novel strategy to enhance the immunosuppressive activity of MSCs in the inflammatory microenvironment.

Epigenetic regulators controlling osteogenic lineage commitment and bone formation

Bone formation and homeostasis are controlled by environmental factors and endocrine regulatory cues that initiate intracellular signaling pathways capable of modulating gene expression in the nucleus. Bone-related gene expression is controlled by nucleosome-based chromatin architecture that limits the accessibility of lineage-specific gene regulatory DNA sequences and sequence-specific transcription factors. From a developmental perspective, bone-specific gene expression must be suppressed during the early stages of embryogenesis to prevent the premature mineralization of skeletal elements during fetal growth in utero. Hence, bone formation is initially inhibited by gene suppressive epigenetic regulators, while other epigenetic regulators actively support osteoblast differentiation. Prominent epigenetic regulators that stimulate or attenuate osteogenesis include lysine methyl transferases (e.g., EZH2, SMYD2, SUV420H2), lysine deacetylases (e.g., HDAC1, HDAC3, HDAC4, HDAC7, SIRT1, SIRT3), arginine methyl transferases (e.g., PRMT1, PRMT4/CARM1, PRMT5), dioxygenases (e.g., TET2), bromodomain proteins (e.g., BRD2, BRD4) and chromodomain proteins (e.g., CBX1, CBX2, CBX5). This narrative review provides a broad overview of the covalent modifications of DNA and histone proteins that involve hundreds of enzymes that add, read, or delete these epigenetic modifications that are relevant for self-renewal and differentiation of mesenchymal stem cells, skeletal stem cells and osteoblasts during osteogenesis.

Transcriptomic landscape of human induced pluripotent stem cell-derived osteogenic differentiation identifies a regulatory role of KLF16

Osteogenic differentiation is essential for bone development and metabolism, but the underlying gene regulatory networks have not been well investigated. We differentiated mesenchymal stem cells, derived from 20 human induced pluripotent stem cell lines, into preosteoblasts and osteoblasts, and performed systematic RNA-seq analyses of 60 samples for differential gene expression. We noted a highly significant correlation in expression patterns and genomic proximity among transcription factor (TF) and long noncoding RNA (lncRNA) genes. We identified TF-TF regulatory networks, regulatory roles of lncRNAs on their neighboring coding genes for TFs and splicing factors, and differential splicing of TF, lncRNA, and splicing factor genes. TF-TF regulatory and gene co-expression network analyses suggested an inhibitory role of TF KLF16 in osteogenic differentiation. We demonstrate that in vitro overexpression of human KLF16 inhibits osteogenic differentiation and mineralization, and in vivo Klf16+/- mice exhibit increased bone mineral density, trabecular number, and cortical bone area. Thus, our model system highlights the regulatory complexity of osteogenic differentiation and identifies novel osteogenic genes.

Unveiling Polysomal Long Non-Coding RNA Expression on the First Day of Adipogenesis and Osteogenesis in Human Adipose-Derived Stem Cells

Understanding the intricate molecular mechanisms governing the fate of human adipose-derived stem cells (hASCs) is essential for elucidating the delicate balance between adipogenic and osteogenic differentiation in both healthy and pathological conditions. Long non-coding RNAs (lncRNAs) have emerged as key regulators involved in lineage commitment and differentiation of stem cells, operating at various levels of gene regulation, including transcriptional, post-transcriptional, and post-translational processes. To gain deeper insights into the role of lncRNAs' in hASCs' differentiation, we conducted a comprehensive analysis of the lncRNA transcriptome (RNA-seq) and translatome (polysomal-RNA-seq) during a 24 h period of adipogenesis and osteogenesis. Our findings revealed distinct expression patterns between the transcriptome and translatome during both differentiation processes, highlighting 90 lncRNAs that are exclusively regulated in the polysomal fraction. These findings underscore the significance of investigating lncRNAs associated with ribosomes, considering their unique expression patterns and potential mechanisms of action, such as translational regulation and potential coding capacity for microproteins. Additionally, we identified specific lncRNA gene expression programs associated with adipogenesis and osteogenesis during the early stages of cell differentiation. By shedding light on the expression and potential functions of these polysome-associated lncRNAs, we aim to deepen our understanding of their involvement in the regulation of adipogenic and osteogenic differentiation, ultimately paving the way for novel therapeutic strategies and insights into regenerative medicine.

DNA methylation restricts coordinated germline and neural fates in embryonic stem cell differentiation

As embryonic stem cells (ESCs) transition from naive to primed pluripotency during early mammalian development, they acquire high DNA methylation levels. During this transition, the germline is specified and undergoes genome-wide DNA demethylation, while emergence of the three somatic germ layers is preceded by acquisition of somatic DNA methylation levels in the primed epiblast. DNA methylation is essential for embryogenesis, but the point at which it becomes critical during differentiation and whether all lineages equally depend on it is unclear. Here, using culture modeling of cellular transitions, we found that DNA methylation-free mouse ESCs with triple DNA methyltransferase knockout (TKO) progressed through the continuum of pluripotency states but demonstrated skewed differentiation abilities toward neural versus other somatic lineages. More saliently, TKO ESCs were fully competent for establishing primordial germ cell-like cells, even showing temporally extended and self-sustained capacity for the germline fate. By mapping chromatin states, we found that neural and germline lineages are linked by a similar enhancer dynamic upon exit from the naive state, defined by common sets of transcription factors, including methyl-sensitive ones, that fail to be decommissioned in the absence of DNA methylation. We propose that DNA methylation controls the temporality of a coordinated neural-germline axis of the preferred differentiation route during early development.

Multiomics analysis identifies novel facilitators of human dopaminergic neuron differentiation

Midbrain dopaminergic neurons (mDANs) control voluntary movement, cognition, and reward behavior under physiological conditions and are implicated in human diseases such as Parkinson's disease (PD). Many transcription factors (TFs) controlling human mDAN differentiation during development have been described, but much of the regulatory landscape remains undefined. Using a tyrosine hydroxylase (TH) human iPSC reporter line, we here generate time series transcriptomic and epigenomic profiles of purified mDANs during differentiation. Integrative analysis predicts novel regulators of mDAN differentiation and super-enhancers are used to identify key TFs. We find LBX1, NHLH1 and NR2F1/2 to promote mDAN differentiation and show that overexpression of either LBX1 or NHLH1 can also improve mDAN specification. A more detailed investigation of TF targets reveals that NHLH1 promotes the induction of neuronal miR-124, LBX1 regulates cholesterol biosynthesis, and NR2F1/2 controls neuronal activity.

The unique role of bone marrow adipose tissue in ovariectomy-induced bone loss in mice

BACKGROUND

Accumulation of bone marrow adipose tissue (BMAT) is always seen in osteoporosis induced by estrogen deficiency. Herein, we aimed to investigate the mechanisms and consequences of this phenomenon by establishing a mouse model of osteoporosis caused by ovariectomy (OVX)-mimicked estrogen deficiency.

METHODS

Micro-CT, osmium tetroxide staining, and histological analyses were performed to examine the changes in bone microstructure, BMAT and white adipose tissue (WAT) in OVX mice compared to sham mice. The osteogenesis and adipogenesis of primary bone marrow stromal cells (BMSCs) isolated from sham and OVX mice were compared in vitro. The molecular phenotypes of BMAT and WAT were determined and compared by quantitative PCR (qPCR). Bone marrow adipocyte-conditioned medium (BMA CM) was prepared from sham or OVX mice for coculture assays, and BMSCs or bone marrow monocytes/macrophages (BMMs) were isolated and subjected to osteoblast and osteoclast differentiation, respectively. Cell staining and qPCR were used to assess the effects of BMAT on bone metabolism.

RESULTS

OVX-induced estrogen deficiency induced reductions in both cortical and trabecular bone mass along with an expansion of BMAT volume. At the cellular level, loss of estrogen inhibited BMSC osteogenesis and promoted BMSC adipogenesis, whereas addition of estradiol exerted the opposite effects. In response to estrogen deficiency, despite the common proinflammatory molecular phenotype observed in both fat depots, BMAT, unlike WAT, unexpectedly exhibited an increase in adipocyte differentiation and lipolytic activity as well as the maintenance of insulin sensitivity. Importantly, BMAT, but not WAT, presented increased mRNA levels of both BMP receptor inhibitors (Grem1, Chrdl1) and Rankl following OVX. In addition, treatment with BMA CM, especially from OVX mice, suppressed the osteoblast differentiation of BMSCs while favoring the osteoclast differentiation of BMMs.

CONCLUSION

Our study illustrates that OVX-induced estrogen deficiency results in bone loss and BMAT expansion by triggering imbalance between the osteogenesis and adipogenesis of BMSCs. Furthermore, expanded BMAT, unlike typical WAT, may negatively regulate bone homeostasis through paracrine inhibition of osteoblast-mediated bone formation and promotion of osteoclast-mediated bone resorption.

Single-cell transcriptome analysis reveals the effectiveness of cytokine priming irrespective of heterogeneity in mesenchymal stromal cells

BACKGROUND AIMS

Mesenchymal stromal cells (MSCs) are recognized as a potential cell-based therapy for regenerative medicine. Short-term inflammatory cytokine pre-stimulation (cytokine priming) is a promising approach to enhance regenerative efficacy of MSCs. However, it is unclear whether their intrinsic heterogenic nature causes an unequal response to cytokine priming, which might blunt the accessibility of clinical applications.

METHODS

In this study, by analyzing the single-cell transcriptomic landscape of human bone marrow MSCs from a naïve to cytokine-primed state, we elucidated the potential mechanism of superior therapeutic potential in cytokine-primed MSCs.

RESULTS

We found that cytokine-primed MSCs had a distinct transcriptome landscape. Although substantial heterogeneity was identified within the population in both naïve and primed states, cytokine priming enhanced the several characteristics of MSCs associated with therapeutic efficacy irrespective of heterogeneity. After cytokine-priming, all sub-clusters of MSCs possessed high levels of immunoregulatory molecules, trophic factors, stemness-related genes, anti-apoptosis markers and low levels of multi-lineage and senescence signatures, which are critical for their therapeutic potency.

CONCLUSIONS

In conclusion, our results provide new insights into MSC heterogeneity under cytokine stimulation and suggest that cytokine priming reprogrammed MSCs independent of heterogeneity.

Metformin potentiates immunosuppressant activity and adipogenic differentiation of human umbilical cord-mesenchymal stem cells

Metformin, a first-line drug for type-2 diabetes, displays pleiotropic effects on inflammation, aging, and cancer. Obesity triggers a low-grade chronic inflammation leading to insulin resistance, characterized by increased pro-inflammatory cytokines produced by adipocytes and infiltrated immune cells, which contributes to metabolic syndrome. We investigated metformin's differentiation and immunoregulatory properties of human umbilical cord-mesenchymal stem cells (UC-MSC), as cellular basis of its beneficial role in metabolic dysfunctions. Isolation, characterization and multilineage differentiation of UC-MSC were performed using standard protocols and flow-cytometry. Metformin effects on UC-MSC growth was assessed by colony formation and MTT assay, gene and protein expression by qRT-PCR, and western blot analysis. Proliferation of peripheral blood mononuclear cells (PBMCs) co-cultured with metformin-treated UC-MSC-conditioned media was evaluated by dye dilution assay. We show that metformin decreases proliferation and colony formation of UC-MSCs and enhances their adipogenic lineage commitment. Metformin (3 mM) increases PPARγ and downregulates FABP4 mRNA both in basal and in adipogenic culture conditions; however, the modulation of PPARγ expression is unrelated to the antiproliferative effects. Moreover, metformin inhibits UC-MSC inflammatory activity reducing the expression of IL-6, MCP-1, and COX-2. Conditioned media, collected from metformin-treated UC-MSCs, down-regulate CD3+ T lymphocyte growth in stimulated PBMCs and, in particular, reduce the CD8+ T cell population. These results indicate that metformin may favor new adipocyte formation and potentiate immune suppressive properties of UC-MSCs. Thus, adipose tissue regeneration and anti-inflammatory activity may represent possible mechanisms by which metformin exerts its positive effect on lipid metabolism.

Nutrient regulation of bone marrow adipose tissue: skeletal implications of weight loss

Adipose tissue is a dynamic component of the bone marrow, regulating skeletal remodelling and secreting paracrine and endocrine factors that can affect haematopoiesis, as well as potentially nourishing the bone marrow during periods of stress. Bone marrow adipose tissue is regulated by multiple factors, but particularly nutrient status. In this Review, we examine how bone marrow adipocytes originate, their function in normal and pathological states and how bone marrow adipose tissue modulates whole-body homoeostasis through actions on bone cells, haematopoietic stem cells and extra-medullary adipocytes during nutritional challenges. We focus on both rodent models and human studies to help understand the unique marrow adipocyte, its response to the external nutrient environment and its effects on the skeleton. We finish by addressing some critical questions that to date remain unanswered.

Rare variant association analysis in 51,256 type 2 diabetes cases and 370,487 controls informs the spectrum of pathogenicity of monogenic diabetes genes

We meta-analyzed array data imputed with the TOPMed reference panel and whole-genome sequence (WGS) datasets and performed the largest, rare variant (minor allele frequency as low as 5×10-5) GWAS meta-analysis of type 2 diabetes (T2D) comprising 51,256 cases and 370,487 controls. We identified 52 novel variants at genome-wide significance (p<5 × 10-8), including 8 novel variants that were either rare or ancestry-specific. Among them, we identified a rare missense variant in HNF4A p.Arg114Trp (OR=8.2, 95% confidence interval [CI]=4.6-14.0, p = 1.08×10-13), previously reported as a variant implicated in Maturity Onset Diabetes of the Young (MODY) with incomplete penetrance. We demonstrated that the diabetes risk in carriers of this variant was modulated by a T2D common variant polygenic risk score (cvPRS) (carriers in the top PRS tertile [OR=18.3, 95%CI=7.2-46.9, p=1.2×10-9] vs carriers in the bottom PRS tertile [OR=2.6, 95% CI=0.97-7.09, p = 0.06]. Association results identified eight variants of intermediate penetrance (OR>5) in monogenic diabetes (MD), which in aggregate as a rare variant PRS were associated with T2D in an independent WGS dataset (OR=4.7, 95% CI=1.86-11.77], p = 0.001). Our data also provided support evidence for 21% of the variants reported in ClinVar in these MD genes as benign based on lack of association with T2D. Our work provides a framework for using rare variant imputation and WGS analyses in large-scale population-based association studies to identify large-effect rare variants and provide evidence for informing variant pathogenicity.

Chinese Ecliptae herba (Eclipta prostrata (L.) L.) extract and its component wedelolactone enhances osteoblastogenesis of bone marrow mesenchymal stem cells via targeting METTL3-mediated m6A RNA methylation

ETHNOPHARMACOLOGICAL RELEVANCE

Chinese Ecliptae herba (Eclipta prostrata (L.) L.) is an ethnomedicinal herb, which is used mainly to nourish kidney and thus strengthen bones according to traditional Chinese medicine theory. Pharmacological studies have supported the ethnomedicine use, showing that Ecliptae herba extract has an anti-osteoporotic effect in vivo and promoted osteoblast proliferation and activity in vitro. However, the molecular mechanism of Ecliptae herba on osteoblast differentiation from bone marrow mesenchymal stem cells (BMSC), the progenitors of osteoblasts, is still unclear.

AIM OF THE STUDY

N6-methyladenosine (m6A) mRNA epigenetic modification may play a key role in promoting osteoblastic differentiation, and thus treating osteoporosis. This study sought to assess the mechanism through which Eclipate herba and its component wedelolactone influence m6A modification during the process of osteoblastogenesis from BMSC.

MATERIAL AND METHODS

The alkaline phosphatase (ALP) and Alizarin red S (ARS) staining were applied to determine osteoblastogenesis from BMSC. Western blot and quantitative real-time PCR were performed. RNA sequencing analysis was used to determine the characteristics of m6A methylation. Stable knocking down of METTL3 using lentiviral-based shRNA was performed.

RESULTS

Upon 9 d treatment of BMSC with ethyl acetate extract of Ecliptae herba (MHL), ALP activity and ossification level increased in comparison with osteogenic medium (OS)-treated control. The expression of methyltransferase METTL3 and METTL14 was significantly increased, but WTAP expression had no change in response to MHL treatment. Knocking down of METTL3 resulted in a decrease in MHL-induced ALP activity, ossification level as well as mRNA expression of Osterix and Osteocalcin, two bone formation-related markers. The level of m6A increased when BMSC was treated with MHL for 9 d. RNA sequencing analysis indicated that MHL treatment altered mRNA m6A modification of genes associated with osteoblastogenesis. By kyoto encyclopedia of genes and genomes (KEGG) pathway analysis, HIF-1α, PI3K/Akt, and Hippo signaling pathways were enriched and associated with m6A modification. The expression of m6A-modified genes including HIF-1α, VEGF-A, and RASSF1, was upregulated by MHL, but the upregulation was reversed after METTL3 knockdown. Additionally, the enhanced expression of METTL3 was also observed after treatment with wedelolactone, a component from MHL.

CONCLUSIONS

These results suggested a previously uncharacterized mechanism of MHL and wedelolactone on osteoblastogenesis, by which METTL3-mediated m6A methylation is involved and thus contributes to the enhancement of osteoblastogenesis.

Increased body mass index is linked to systemic inflammation through altered chromatin co-accessibility in human preadipocytes

Obesity-induced adipose tissue dysfunction can cause low-grade inflammation and downstream obesity comorbidities. Although preadipocytes may contribute to this pro-inflammatory environment, the underlying mechanisms are unclear. We used human primary preadipocytes from body mass index (BMI) -discordant monozygotic (MZ) twin pairs to generate epigenetic (ATAC-sequence) and transcriptomic (RNA-sequence) data for testing whether increased BMI alters the subnuclear compartmentalization of open chromatin in the twins' preadipocytes, causing downstream inflammation. Here we show that the co-accessibility of open chromatin, i.e. compartmentalization of chromatin activity, is altered in the higher vs lower BMI MZ siblings for a large subset ( ~ 88.5 Mb) of the active subnuclear compartments. Using the UK Biobank we show that variants within these regions contribute to systemic inflammation through interactions with BMI on C-reactive protein. In summary, open chromatin co-accessibility in human preadipocytes is disrupted among the higher BMI siblings, suggesting a mechanism how obesity may lead to inflammation via gene-environment interactions.

Single-Cell RNA-Sequencing Reveals the Skeletal Cellular Dynamics in Bone Repair and Osteoporosis

The bone is an important organ that performs various functions, and the bone marrow inside the skeleton is composed of a complex intermix of hematopoietic, vascular, and skeletal cells. Current single-cell RNA sequencing (scRNA-seq) technology has revealed heterogeneity and sketchy differential hierarchy of skeletal cells. Skeletal stem and progenitor cells (SSPCs) are located upstream of the hierarchy and differentiate into chondrocytes, osteoblasts, osteocytes, and bone marrow adipocytes. In the bone marrow, multiple types of bone marrow stromal cells (BMSCs), which have the potential of SSPCs, are spatiotemporally located in distinct areas, and SSPCs' potential shift of BMSCs may occur with the advancement of age. These BMSCs contribute to bone regeneration and bone diseases, such as osteoporosis. In vivo lineage-tracing technologies show that various types of skeletal lineage cells concomitantly gather and contribute to bone regeneration. In contrast, these cells differentiate into adipocytes with aging, leading to senile osteoporosis. scRNA-seq analysis has revealed that alteration in the cell-type composition is a major cause of tissue aging. In this review, we discuss the cellular dynamics of skeletal cell populations in bone homeostasis, regeneration, and osteoporosis.

Adipogenic and SWAT cells separate from a common progenitor in human brown and white adipose depots

Adipocyte function is a major determinant of metabolic disease, warranting investigations of regulating mechanisms. We show at single-cell resolution that progenitor cells from four human brown and white adipose depots separate into two main cell fates, an adipogenic and a structural branch, developing from a common progenitor. The adipogenic gene signature contains mitochondrial activity genes, and associates with genome-wide association study traits for fat distribution. Based on an extracellular matrix and developmental gene signature, we name the structural branch of cells structural Wnt-regulated adipose tissue-resident (SWAT) cells. When stripped from adipogenic cells, SWAT cells display a multipotent phenotype by reverting towards progenitor state or differentiating into new adipogenic cells, dependent on media. Label transfer algorithms recapitulate the cell types in human adipose tissue datasets. In conclusion, we provide a differentiation map of human adipocytes and define the multipotent SWAT cell, providing a new perspective on adipose tissue regulation.

Systematic dissection of coordinated stromal remodeling identifies Sox10+ glial cells as a therapeutic target in myelofibrosis

Remodeling of the tissue niche is often evident in diseases, yet, the stromal alterations and their contribution to pathogenesis are poorly characterized. Bone marrow fibrosis is a maladaptive feature of primary myelofibrosis (PMF). We performed lineage tracing and found that most collagen-expressing myofibroblasts were derived from leptin-receptor-positive (LepR+) mesenchymal cells, whereas a minority were from Gli1-lineage cells. Deletion of Gli1 did not impact PMF. Unbiased single-cell RNA sequencing (scRNA-seq) confirmed that virtually all myofibroblasts originated from LepR-lineage cells, with reduced expression of hematopoietic niche factors and increased expression of fibrogenic factors. Concurrently, endothelial cells upregulated arteriolar-signature genes. Pericytes and Sox10+ glial cells expanded drastically with heightened cell-cell signaling, suggesting important functional roles in PMF. Chemical or genetic ablation of bone marrow glial cells ameliorated fibrosis and improved other pathology in PMF. Thus, PMF involves complex remodeling of the bone marrow microenvironment, and glial cells represent a promising therapeutic target.

Recent advances of osteoimmunology research in rheumatoid arthritis: From single-cell omics approach

ABSTRACT

Cellular immune responses as well as generalized and periarticular bone loss are the key pathogenic features of rheumatoid arthritis (RA). Under the pathological conditions of RA, dysregulated inflammation and immune processes tightly interact with skeletal system, resulting in pathological bone damage via inhibition of bone formation or induction of bone resorption. Single-cell omics technologies are revolutionary tools in the field of modern biological research.They enable the display of the state and function of cells in various environments from a single-cell resolution, thus making it conducive to identify the dysregulated molecular mechanisms of bone destruction in RA as well as the discovery of potential therapeutic targets and biomarkers. Here, we summarize the latest findings of single-cell omics technologies in osteoimmunology research in RA. These results suggest that single-cell omics have made significant contributions to transcriptomics and dynamics of specific cells involved in bone remodeling, providing a new direction for our understanding of cellular heterogeneity in the study of osteoimmunology in RA.

Assessing Biomaterial-Induced Stem Cell Lineage Fate by Machine Learning-Based Artificial Intelligence

Current functional assessment of biomaterial-induced stem cell lineage fate in vitro mainly relies on biomarker-dependent methods with limited accuracy and efficiency. Here a "Mesenchymal stem cell Differentiation Prediction (MeD-P)" framework for biomaterial-induced cell lineage fate prediction is reported. MeD-P contains a cell-type-specific gene expression profile as a reference by integrating public RNA-seq data related to tri-lineage differentiation (osteogenesis, chondrogenesis, and adipogenesis) of human mesenchymal stem cells (hMSCs) and a predictive model for classifying hMSCs differentiation lineages using the k-nearest neighbors (kNN) strategy. It is shown that MeD-P exhibits an overall accuracy of 90.63% on testing datasets, which is significantly higher than the model constructed based on canonical marker genes (80.21%). Moreover, evaluations of multiple biomaterials show that MeD-P provides accurate prediction of lineage fate on different types of biomaterials as early as the first week of hMSCs culture. In summary, it is demonstrated that MeD-P is an efficient and accurate strategy for stem cell lineage fate prediction and preliminary biomaterial functional evaluation.

The enhancer landscape predetermines the skeletal regeneration capacity of stromal cells

Multipotent stromal cells are considered attractive sources for cell therapy and tissue engineering. Despite numerous experimental and clinical studies, broad application of stromal cell therapeutics is not yet emerging. A major challenge is the functional diversity of available cell sources. Here, we investigated the regenerative potential of clinically relevant human stromal cells from bone marrow (BMSCs), white adipose tissue, and umbilical cord compared with mature chondrocytes and skin fibroblasts in vitro and in vivo. Although all stromal cell types could express transcription factors related to endochondral ossification, only BMSCs formed cartilage discs in vitro that fully regenerated critical-size femoral defects after transplantation into mice. We identified cell type-specific epigenetic landscapes as the underlying molecular mechanism controlling transcriptional stromal differentiation networks. Binding sites of commonly expressed transcription factors in the enhancer and promoter regions of ossification-related genes, including Runt and bZIP families, were accessible only in BMSCs but not in extraskeletal stromal cells. This suggests an epigenetically predetermined differentiation potential depending on cell origin that allows common transcription factors to trigger distinct organ-specific transcriptional programs, facilitating forward selection of regeneration-competent cell sources. Last, we demonstrate that viable human BMSCs initiated defect healing through the secretion of osteopontin and contributed to transient mineralized bone hard callus formation after transplantation into immunodeficient mice, which was eventually replaced by murine recipient bone during final tissue remodeling.

Gene silencing dynamics are modulated by transiently active regulatory elements

Transcriptional enhancers have been extensively characterized, but cis-regulatory elements involved in acute gene repression have received less attention. Transcription factor GATA1 promotes erythroid differentiation by activating and repressing distinct gene sets. Here, we study the mechanism by which GATA1 silences the proliferative gene Kit during murine erythroid cell maturation and define stages from initial loss of activation to heterochromatinization. We find that GATA1 inactivates a potent upstream enhancer but concomitantly creates a discrete intronic regulatory region marked by H3K27ac, short noncoding RNAs, and de novo chromatin looping. This enhancer-like element forms transiently and serves to delay Kit silencing. The element is ultimately erased via the FOG1/NuRD deacetylase complex, as revealed by the study of a disease-associated GATA1 variant. Hence, regulatory sites can be self-limiting by dynamic co-factor usage. Genome-wide analyses across cell types and species uncover transiently active elements at numerous genes during repression, suggesting that modulation of silencing kinetics is widespread.

O-GlcNAc glycosylation orchestrates fate decision and niche function of bone marrow stromal progenitors

In mammals, interactions between the bone marrow (BM) stroma and hematopoietic progenitors contribute to bone-BM homeostasis. Perinatal bone growth and ossification provide a microenvironment for the transition to definitive hematopoiesis; however, mechanisms and interactions orchestrating the development of skeletal and hematopoietic systems remain largely unknown. Here, we establish intracellular O-linked β-N-acetylglucosamine (O-GlcNAc) modification as a posttranslational switch that dictates the differentiation fate and niche function of early BM stromal cells (BMSCs). By modifying and activating RUNX2, O-GlcNAcylation promotes osteogenic differentiation of BMSCs and stromal IL-7 expression to support lymphopoiesis. In contrast, C/EBPβ-dependent marrow adipogenesis and expression of myelopoietic stem cell factor (SCF) is inhibited by O-GlcNAcylation. Ablating O-GlcNAc transferase (OGT) in BMSCs leads to impaired bone formation, increased marrow adiposity, as well as defective B-cell lymphopoiesis and myeloid overproduction in mice. Thus, the balance of osteogenic and adipogenic differentiation of BMSCs is determined by reciprocal O-GlcNAc regulation of transcription factors, which simultaneously shapes the hematopoietic niche.

Emerging RUNX2-Mediated Gene Regulatory Mechanisms Consisting of Multi-Layered Regulatory Networks in Skeletal Development

Skeletal development is tightly coordinated by chondrocytes and osteoblasts, which are derived from skeletal progenitors, and distinct cell-type gene regulatory programs underlie the specification and differentiation of cells. Runt-related transcription factor 2 (Runx2) is essential to chondrocyte hypertrophy and osteoblast differentiation. Genetic studies have revealed the biological functions of Runx2 and its involvement in skeletal genetic diseases. Meanwhile, molecular biology has provided a framework for our understanding of RUNX2-mediated transactivation at a limited number of cis-regulatory elements. Furthermore, studies using next-generation sequencing (NGS) have provided information on RUNX2-mediated gene regulation at the genome level and novel insights into the multiple layers of gene regulatory mechanisms, including the modes of action of RUNX2, chromatin accessibility, the concept of pioneer factors and phase separation, and three-dimensional chromatin organization. In this review, I summarize the emerging RUNX2-mediated regulatory mechanism from a multi-layer perspective and discuss future perspectives for applications in the treatment of skeletal diseases.

Molecular Dissection of Somatic Skeletal Disease in Neurofibromatosis Type 1

Neurofibromatosis type 1 (NF1) is a tumor predisposition syndrome caused by heterozygous NF1 gene mutations. Patients with NF1 present with pleiotropic somatic secondary manifestations, including development of bone pseudarthrosis after fracture. Somatic NF1 gene mutations were reproducibly identified in patient-derived pseudarthrosis specimens, suggesting a local mosaic cell population including somatic pathologic cells. The somatic cellular pathogenesis of NF1 pseudarthroses remains unclear, though defects in osteogenesis have been posited. Here, we applied time-series single-cell RNA-sequencing (scRNA-seq) to patient-matched control and pseudarthrosis-derived primary bone stromal cells (BSCs). We show that osteogenic specification to an osteoblast progenitor cell population was evident for control bone-derived cells and haploinsufficient pseudarthrosis-derived cells. Similar results were observed for somatic patient fracture-derived NF1^-/- cells; however, expression of genetic pathways associated with skeletal mineralization were significantly reduced in NF1^-/- cells compared with fracture-derived NF1^+/- cells. In mice, we show that Nf1 expressed in bone marrow osteoprogenitors is required for the maintenance of the adult skeleton. Results from our study implicate impaired Clec11a-Itga11-Wnt signaling in the pathogenesis of NF1-associated skeletal disease. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Screening for genes, miRNAs and transcription factors of adipogenic differentiation and dedifferentiation of mesenchymal stem cells

BACKGROUND

The purpose of present study was to reveal the molecular mechanisms responsible for both adipogenic differentiation and dedifferentiation of mesenchymal stem cells (MSCs).

METHODS

Microarray data GSE36923 were obtained from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) between adipogenically differentiated cells vs undifferentiated bone marrow-derived MSCs, adipogenically differentiated cells vs dedifferentiated cells samples at day 7 and adipogenically differentiated cells vs dedifferentiated cells samples at day 35 were screened, and overlapped DEGs across the three groups were analyzed. The underlying functions of the upregulated and downregulated DEGs were investigated by Gene ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analysis. The protein-protein interaction network was constructed, and hub genes were obtained subsequently. Hub genes were verified with GSE113253 dataset, and then miRNA-gene network and TF-gene network were constructed.

RESULTS

A total of 284 upregulated DEGs and 376 downregulated DEGs overlapped across the three groups. PPAR signaling pathway, AMPK signaling pathway, insulin signaling pathway, carbon metabolism, pyruvate metabolism, fatty acid metabolism, regulation of lipolysis in adipocytes, biosynthesis of amino acids, citrate cycle (TCA cycle) and 2-Oxocarboxylic acid metabolism were the top 10 pathways involving in the upregulated DEGs, and graft-versus-host disease, allograft rejection, viral myocarditis, cell adhesion molecules, phagosome, type I diabetes mellitus, antigen processing and presentation, autoimmune thyroid disease, intestinal immune network for IgA production and rheumatoid arthritis were the top 10 pathways in downregulated DEGs. After validation, the 8 hub genes were IL6, PPARG, CCL2, FASN, CEBPA, ADIPOQ, FABP4 and LIPE. Ten key miRNAs were hsa-mir-27a-3p, hsa-mir-182-5p, hsa-mir-7-5p, hsa-mir-16-5p, hsa-mir-1-3p, hsa-mir-155-5p, hsa-mir-21-3p, hsa-mir-34a-5p, hsa-mir-27a-5p and hsa-mir-30c-5p, and 10 key TFs were TFDP1, GTF2A2, ZNF584, NRF1, ZNF512, NFRKB, CEBPG, KLF16, GLIS2 and MXD4.

CONCLUSION

Our study constructed miRNA-gene network and TF-gene network involved in both adipogenic differentiation and dedifferentiation of MSCs, contributing to enhancing the efficiency of MSCs transplantation in soft tissue defect repair and developing more potent remedies for adipogenesis-related skeletal disorders.

Identifying Biomarkers for Osteogenic Potency Assay Development

There has been extensive exploration of how cells may serve as advanced therapy medicinal products to treat skeletal pathologies. Osteoblast progenitors responsible for production of extracellular matrix that is subsequently mineralized during bone formation have been characterised as a rare bone marrow subpopulation of cell culture plastic adherent cells. Conveniently, they proliferate to form single-cell derived colonies of fibroblastoid cells, termed colony forming unit fibroblasts that can subsequently differentiate to aggregates resembling small areas of cartilage or bone. However, donor heterogeneity and loss of osteogenic differentiation capacity during extended cell culture have made the discovery of reliable potency assay biomarkers difficult. Nonetheless, functional osteoblast models derived from telomerised human bone marrow stromal cells have allowed extensive comparative analysis of gene expression, microRNA, morphological phenotypes and secreted proteins. This chapter highlights numerous insights into the molecular mechanisms underpinning osteogenic differentiation of multipotent stromal cells and bone formation, discussing aspects involved in the choice of useful biomarkers for functional attributes that can be quantitively measured in osteogenic potency assays.

Improved Protocol to Study Osteoblast and Adipocyte Differentiation Balance

Adipogenesis-osteoblastogenesis balance-rupture is relevant in multiple diseases. Current human mesenchymal stem cells (hMSCs) in vitro differentiation models are expensive, and are hardly reproducible. Their scarcity and variability make an affordable and reliable method to study adipocyte-osteoblast-equilibrium difficult. Moreover, media composition has been inconstant throughout the literature. Our aims were to compare improved differentiation lab-made media with consensus/commercial media, and to identify a cell-line to simultaneously evaluate both MSCs differentiations. Lab-made media were compared with consensus and commercial media in C3H10T1/2 and hMSC, respectively. Lab-made media were tested on aged women primary pre-osteoblast-like cells. To determine the optimum cell line, C3H10T1/2 and hMSC-TERT cells were differentiated to both cell fates. Differentiation processes were evaluated by adipocytic and osteoblastic gene-markers expression and staining. Lab-made media significantly increased consensus medium induction and overcame commercial media in hMSCs differentiation to adipocytes and osteoblasts. Pre-osteoblast-like cells only properly differentiate to adipocyte. Lab-made media promoted adipocyte gene-markers expression in C3H10T1/2 and hMSC-TERT, and osteoblast gene-markers in C3H10T1/2. Oil Red O and Alizarin Red staining supported these findings. Optimized lab-made media were better at differentiating MSCs compared to consensus/commercial media, and evidenced the adipogenic commitment of pre-osteoblast-like cells from aged-women. C3H10T1/2 is an optimum MSC line by which to study adipocyte-osteoblast differentiation balance.

Lineage-specific rearrangement of chromatin loops and epigenomic features during adipocytes and osteoblasts commitment

Human mesenchymal stem cells (hMSCs) can be differentiated into adipocytes and osteoblasts. The processes are driven by the rewiring of chromatin architectures and transcriptomic/epigenomic changes. Here, we induced hMSCs to adipogenic and osteogenic differentiation, and performed 2 kb resolution Hi-C experiments for chromatin loops detection. We also generated matched RNA-seq, ChIP-seq and ATAC-seq data for integrative analysis. After comprehensively comparing adipogenesis and osteogenesis, we quantitatively identified lineage-specific loops and screened out lineage-specific enhancers and open chromatin. We reveal that lineage-specific loops can activate gene expression and facilitate cell commitment through combining enhancers and accessible chromatin in a lineage-specific manner. We finally proposed loop-mediated regulatory networks and identified the controlling factors for adipocytes and osteoblasts determination. Functional experiments validated the lineage-specific regulation networks towards IRS2 and RUNX2 that are associated with adipogenesis and osteogenesis, respectively. These results are expected to help better understand the chromatin conformation determinants of hMSCs fate commitment.

Long-range chromosomal interactions increase and mark repressed gene expression during adipogenesis

Obesity perturbs central functions of human adipose tissue, centred on differentiation of preadipocytes to adipocytes, i.e., adipogenesis. The large environmental component of obesity makes it important to elucidate epigenetic regulatory factors impacting adipogenesis. Promoter Capture Hi-C (pCHi-C) has been used to identify chromosomal interactions between promoters and associated regulatory elements. However, long range interactions (LRIs) greater than 1 Mb are often filtered out of pCHi-C datasets, due to technical challenges and their low prevalence. To elucidate the unknown role of LRIs in adipogenesis, we investigated preadipocyte differentiation to adipocytes using pCHi-C and bulk and single nucleus RNA-seq data. We first show that LRIs are reproducible between biological replicates, and they increase >2-fold in frequency across adipogenesis. We further demonstrate that genomic loci containing LRIs are more epigenetically repressed than regions without LRIs, corresponding to lower gene expression in the LRI regions. Accordingly, as preadipocytes differentiate into adipocytes, LRI regions are more likely to contain repressed preadipocyte marker genes; whereas these same LRI regions are depleted of actively expressed adipocyte marker genes. Finally, we show that LRIs can be used to restrict multiple testing of the long-range cis-eQTL analysis to identify variants that regulate genes via LRIs. We exemplify this by identifying a putative long range cis regulatory mechanism at the LYPLAL1/TGFB2 obesity locus. In summary, we identify LRIs that mark repressed regions of the genome, and these interactions increase across adipogenesis, pinpointing developmental regions that need to be repressed in a cell-type specific way for adipogenesis to proceed.

Comprehensive comparison of gene expression diversity among a variety of human stem cells

Several factors, including tissue origins and culture conditions, affect the gene expression of undifferentiated stem cells. However, understanding the basic identity across different stem cells has not been pursued well despite its importance in stem cell biology. Thus, we aimed to rank the relative importance of multiple factors to gene expression profile among undifferentiated human stem cells by analyzing publicly available RNA-seq datasets. We first conducted batch effect correction to avoid undefined variance in the dataset as possible. Then, we highlighted the relative impact of biological and technical factors among undifferentiated stem cell types: a more influence on tissue origins in induced pluripotent stem cells than in other stem cell types; a stronger impact of culture condition in embryonic stem cells and somatic stem cell types, including mesenchymal stem cells and hematopoietic stem cells. In addition, we found that a characteristic gene module, enriched in histones, exhibits higher expression across different stem cell types that were annotated by specific culture conditions. This tendency was also observed in mouse stem cell RNA-seq data. Our findings would help to obtain general insights into stem cell quality, such as the balance of differentiation potentials that undifferentiated stem cells possess.

Insights into skeletal stem cells

The tissue-resident skeletal stem cells (SSCs), which are self-renewal and multipotent, continuously provide cells (including chondrocytes, bone cells, marrow adipocytes, and stromal cells) for the development and homeostasis of the skeletal system. In recent decade, utilizing fluorescence-activated cell sorting, lineage tracing, and single-cell sequencing, studies have identified various types of SSCs, plotted the lineage commitment trajectory, and partially revealed their properties under physiological and pathological conditions. In this review, we retrospect to SSCs identification and functional studies. We discuss the principles and approaches to identify bona fide SSCs, highlighting pioneering findings that plot the lineage atlas of SSCs. The roles of SSCs and progenitors in long bone, craniofacial tissues, and periosteum are systematically discussed. We further focus on disputes and challenges in SSC research.

Lineage-selective super enhancers mediate core regulatory circuitry during adipogenic and osteogenic differentiation of human mesenchymal stem cells

Human mesenchymal stem cells (hMSCs) can be differentiated into osteoblasts and adipocytes. During these processes, super enhancers (SEs) play important roles. Here, we performed comprehensive characterization of the SEs changes associated with adipogenic and osteogenic differentiation of hMSCs, and revealed that SEs changed more dramatically compared with typical enhancers. We identified a set of lineage-selective SEs, whose target genes were enriched with cell type-specific functions. Functional experiments in lineage-selective SEs demonstrated their specific roles in directed differentiation of hMSCs. We also found that some key transcription factors regulated by lineage-selective SEs could form core regulatory circuitry (CRC) to regulate each other's expression and control the hMSCs fate determination. In addition, we found that GWAS SNPs of osteoporosis and obesity were significantly enriched in osteoblasts-selective SEs or adipocytes-selective SEs, respectively. Taken together, our studies unveiled important roles of lineage-selective SEs in hMSCs differentiation into osteoblasts and adipocytes.

Vitamin C epigenetically controls osteogenesis and bone mineralization

Vitamin C deficiency disrupts the integrity of connective tissues including bone. For decades this function has been primarily attributed to Vitamin C as a cofactor for collagen maturation. Here, we demonstrate that Vitamin C epigenetically orchestrates osteogenic differentiation and function by modulating chromatin accessibility and priming transcriptional activity. Vitamin C regulates histone demethylation (H3K9me3 and H3K27me3) and promotes TET-mediated 5hmC DNA hydroxymethylation at promoters, enhancers and super-enhancers near bone-specific genes. This epigenetic circuit licenses osteoblastogenesis by permitting the expression of all major pro-osteogenic genes. Osteogenic cell differentiation is strictly and continuously dependent on Vitamin C, whereas Vitamin C is dispensable for adipogenesis. Importantly, deletion of 5hmC-writers, Tet1 and Tet2, in Vitamin C-sufficient murine bone causes severe skeletal defects which mimic bone phenotypes of Vitamin C-insufficient Gulo knockout mice, a model of Vitamin C deficiency and scurvy. Thus, Vitamin C's epigenetic functions are central to osteoblastogenesis and bone formation and may be leveraged to prevent common bone-degenerating conditions.

High-fat diet decreases H3K27ac in mice adipose-derived stromal cells

OBJECTIVE

The study goal was to analyze the effects of a high-fat diet (HFD) on the histone 3 lysine 27 (H3K27) posttranscriptional modifications and the expression of histone-modifying enzymes in adipose-derived stromal cells (ASCs) from white adipose tissue (WAT).

METHODS

Male C57BL/6J mice received control or HFD for 12 weeks. The ASCs were isolated from subcutaneous and visceral (epididymal) WAT, cultivated, and evaluated for expression of H3K27 trimethylation (H3K27me3) and H3K27 acetylation (H3K27ac) by Western blot. The transcription of histone-modifying enzymes was analyzed by real-time polymerase chain reaction.

RESULTS

When compared with control, HFD ASCs showed a decrease in H3K27ac enrichment in subcutaneous and visceral WAT and ATP-citrate lyase expression in subcutaneous WAT. Curiously, the expression of CREB-binding protein was increased in visceral ASCs from HFD-fed mice.

CONCLUSIONS

These results show that an HFD significantly reduces acetylation of H3K27 in ASCs and the expression of ATP-citrate lyase in subcutaneous ASCs, suggesting that, in this fat depot, the H3K27ac reduction could be partly due to lower acetyl-coenzyme A availability. H3K27ac is an epigenetic mark responsible for increasing the transcription rate and its reduction can have an important impact on ASC proliferation and differentiation potential.

TCF7L2 plays a complex role in human adipose progenitor biology, which might contribute to genetic susceptibility to type 2 diabetes

INTRODUCTION

Non-coding genetic variation at TCF7L2 is the strongest genetic determinant of type 2 diabetes (T2D) risk in humans. TCF7L2 encodes a transcription factor mediating the nuclear effects of WNT signaling in adipose tissue (AT). In vivo studies in transgenic mice have highlighted important roles for TCF7L2 in adipose tissue biology and systemic metabolism.

OBJECTIVE

To map the expression of TCF7L2 in human AT, examine its role in human adipose cell biology in vitro, and investigate the effects of the fine-mapped T2D-risk allele at rs7903146 on AT morphology and TCF7L2 expression.

METHODS

Ex vivo gene expression studies of TCF7L2 in whole and fractionated human AT. In vitro TCF7L2 gain- and/or loss-of-function studies in primary and immortalized human adipose progenitor cells (APCs) and mature adipocytes (mADs). AT phenotyping of rs7903146 T2D-risk variant carriers and matched controls.

RESULTS

Adipose progenitors (APs) exhibited the highest TCF7L2 mRNA abundance compared to mature adipocytes and adipose-derived endothelial cells. Obesity was associated with reduced TCF7L2 transcript levels in whole subcutaneous abdominal AT but paradoxically increased expression in APs. In functional studies, TCF7L2 knockdown (KD) in abdominal APs led to dose-dependent activation of WNT/β-catenin signaling, impaired proliferation and dose-dependent effects on adipogenesis. Whilst partial KD enhanced adipocyte differentiation, near-total KD impaired lipid accumulation and adipogenic gene expression. Over-expression of TCF7L2 accelerated adipogenesis. In contrast, TCF7L2-KD in gluteal APs dose-dependently enhanced lipid accumulation. Transcriptome-wide profiling revealed that TCF7L2 might modulate multiple aspects of AP biology including extracellular matrix secretion, immune signaling and apoptosis. The T2D-risk allele at rs7903146 was associated with reduced AP TCF7L2 expression and enhanced AT insulin sensitivity.

CONCLUSIONS

TCF7L2 plays a complex role in AP biology and has both dose- and depot-dependent effects on adipogenesis. In addition to regulating pancreatic insulin secretion, genetic variation at TCF7L2 might also influence T2D risk by modulating AP function.

Lipolysis regulates major transcriptional programs in brown adipocytes

β-Adrenergic signaling is a core regulator of brown adipocyte function stimulating both lipolysis and transcription of thermogenic genes, thereby expanding the capacity for oxidative metabolism. We have used pharmacological inhibitors and a direct activator of lipolysis to acutely modulate the activity of lipases, thereby enabling us to uncover lipolysis-dependent signaling pathways downstream of β-adrenergic signaling in cultured brown adipocytes. Here we show that induction of lipolysis leads to acute induction of several gene programs and is required for transcriptional regulation by β-adrenergic signals. Using machine-learning algorithms to infer causal transcription factors, we show that PPARs are key mediators of lipolysis-induced activation of genes involved in lipid metabolism and thermogenesis. Importantly, however, lipolysis also activates the unfolded protein response and regulates the core circadian transcriptional machinery independently of PPARs. Our results demonstrate that lipolysis generates important metabolic signals that exert profound pleiotropic effects on transcription and function of cultured brown adipocytes.

LncMIR181A1HG is a novel chromatin-bound epigenetic suppressor of early stage osteogenic lineage commitment

Bone formation requires osteogenic differentiation of multipotent mesenchymal stromal cells (MSCs) and lineage progression of committed osteoblast precursors. Osteogenic phenotype commitment is epigenetically controlled by genomic (chromatin) and non-genomic (non-coding RNA) mechanisms. Control of osteogenesis by long non-coding RNAs remains a largely unexplored molecular frontier. Here, we performed comprehensive transcriptome analysis at early stages of osteogenic cell fate determination in human MSCs, focusing on expression of lncRNAs. We identified a chromatin-bound lncRNA (MIR181A1HG) that is highly expressed in self-renewing MSCs. MIR181A1HG is down-regulated when MSCs become osteogenic lineage committed and is retained during adipogenic differentiation, suggesting lineage-related molecular functions. Consistent with a key role in human MSC proliferation and survival, we demonstrate that knockdown of MIR181A1HG in the absence of osteogenic stimuli impedes cell cycle progression. Loss of MIR181A1HG enhances differentiation into osteo-chondroprogenitors that produce multiple extracellular matrix proteins. RNA-seq analysis shows that loss of chromatin-bound MIR181A1HG alters expression and BMP2 responsiveness of skeletal gene networks (e.g., SOX5 and DLX5). We propose that MIR181A1HG is a novel epigenetic regulator of early stages of mesenchymal lineage commitment towards osteo-chondroprogenitors. This discovery permits consideration of MIR181A1HG and its associated regulatory pathways as targets for promoting new bone formation in skeletal disorders.

Type 2 Diabetes Mellitus Facilitates Shift of Adipose-Derived Stem Cells Ex Vivo Differentiation toward Osteogenesis among Patients with Obesity

Objective: Sedentary behavior with overnutrition provokes the development of obesity, insulin resistance, and type 2 diabetes mellitus (T2DM). The main progenitor cells of adipose tissue are adipose-derived stem cells (ADSCs) which can change differentiation, metabolic, and secretory phenotypes under obesity conditions. The purpose of this study was to evaluate ADSC osteogenesis activity among patients with obesity in normal glucose tolerance (NGT) and T2DM conditions. Methods: In the study, ADSCs from donors with obesity were used. After clinical characterization, all patients underwent bariatric surgery and ADSCs were isolated from subcutaneous fat biopsies. ADSCs were subjected to osteogenic differentiation, stained with Alizarin Red S, and harvested for real-time PCR and Western blotting. Cell senescence was evaluated with a β-galactosidase-activity-based assay. Results: Our results demonstrated the significantly increased calcification of ADSC on day 28 of osteogenesis in the T2DM group. These data were confirmed by the statistically significant enhancement of RUNX2 gene expression, which is a master regulator of osteogenesis. Protein expression analysis showed the increased expression of syndecan 1 and collagen I before and during osteogenesis, respectively. Moreover, T2DM ADSCs demonstrated an increased level of cellular senescence. Conclusion: We suggest that T2DM-associated cellular senescence can cause ADSC differentiation to shift toward osteogenesis, the impaired formation of new fat depots in adipose tissue, and the development of insulin resistance. The balance between ADSC adipo- and osteogenesis commitment is crucial for the determination of the metabolic fate of patients and their adipose tissue.

Osteogenic capacity of mesenchymal stem cells from patients with osteoporotic hip fractures in vivo

PURPOSE

Human mesenchymal stem cells (MSCs) have the ability to differentiate into bone-forming osteoblasts. The aim of this study was to elucidate if MSCs from patients with OP show a senescent phenotype and explore their bone-forming ability in vivo.

MATERIALS AND METHODS

MSCs from patients with OP and controls with osteoarthritis (OA) were implanted into the subcutaneous tissue of immunodeficient mice for histological analysis and expression of human genes by RT-PCR. The expression of senescence-associated phenotype (SASP) genes, as well as p16, p21, and galactosidase, was studied in cultures of MSCs.

RESULTS

In vivo bone formation was evaluated in 103 implants (47 OP, 56 OA). New bone was observed in 45% of the implants with OP cells and 46% of those with OA cells (p = 0.99). The expression of several bone-related genes (collagen, osteocalcin, alkaline phosphatase, sialoprotein) was also similar in both groups. There were no differences between groups in SASP gene expression, p16, and p21 expression, or in senescence-associated galactosidase activity.

CONCLUSION

Senescence markers and the osteogenic capacity in vivo of MSCs from patients with OP are not inferior to that of cells from controls of similar age with OA. This supports the interest of future studies to evaluate the potential use of autologous MSCs from OP patients in bone regeneration procedures.

Osteogenic Commitment of Human Periodontal Ligament Cells Is Predetermined by Methylation, Chromatin Accessibility and Expression of Key Transcription Factors

Periodontal ligament stem cells (PDLCs) can be used as a valuable source in cell therapies to regenerate bone tissue. However, the potential therapeutic outcomes are unpredictable due to PDLCs’ heterogeneity regarding the capacity for osteoblast differentiation and mineral nodules production. Here, we identify epigenetic (DNA (hydroxy)methylation), chromatin (ATAC-seq) and transcriptional (RNA-seq) differences between PDLCs presenting with low (l) and high (h) osteogenic potential. The primary cell populations were investigated at basal state (cultured in DMEM) and after 10 days of osteogenic stimulation (OM). At a basal state, the expression of transcription factors (TFs) and the presence of gene regulatory regions related to osteogenesis were detected in h-PDLCs in contrast to neuronal differentiation prevalent in l-PDLCs. These differences were also observed under stimulated conditions, with genes and biological processes associated with osteoblast phenotype activated more in h-PDLCs. Importantly, even after the induction, l-PDLCs showed hypermethylation and low expression of genes related to bone development. Furthermore, the analysis of TFs motifs combined with TFs expression suggested the relevance of SP1, SP7 and DLX4 regulation in h-PDLCs, while motifs for SIX and OLIG2 TFs were uniquely enriched in l-PDLCs. Additional analysis including a second l-PDLC population indicated that the high expression of OCT4, SIX3 and PPARG TFs could be predictive of low osteogenic commitment. In summary, several biological processes related to osteoblast commitment were activated in h-PDLCs from the onset, while l-PDLCs showed delay in the activation of the osteoblastic program, restricted by the persistent methylation of gene related to bone development. These processes are pre-determined by distinguishable epigenetic and transcriptional patterns, the recognition of which could help in selection of PDLCs with pre-osteoblastic phenotype.

Periosteal Skeletal Stem Cells and Their Response to Bone Injury

Bone exhibits remarkable self-repair ability without fibrous scars. It is believed that the robust regenerative capacity comes from tissue-resident stem cells, such as skeletal stem cells (SSCs). Roughly, SSC has two niches: bone marrow (BM) and periosteum. BM-SSCs have been extensively studied for years. In contrast, our knowledge about periosteal SSCs (P-SSCs) is quite limited. There is abundant clinical evidence for the presence of stem cell populations within the periosteum. Researchers have even successfully cultured “stem-like” cells from the periosteum in vitro. However, due to the lack of effective markers, it is difficult to evaluate the stemness of real P-SSCs in vivo. Recently, several research teams have developed strategies for the successful identification of P-SSCs. For the first time, we can assess the stemness of P-SSCs from visual evidence. BM-SSCs and P-SSCs not only have much in common but also share distinct properties. Here, we provide an updated review of P-SSCs and their particular responses to bone injury.

Osteogenesis-Related Long Noncoding RNA GAS5 as a Novel Biomarker for Osteonecrosis of Femoral Head

Background: The lack of effective biomarkers makes it difficult to achieve early diagnosis and intervention for osteonecrosis of the femoral head (ONFH). Hence, we aimed to identify novel long noncoding RNA (lncRNA) biomarkers for ONFH.

Methods: High-throughput RNA sequencing was performed to detect lncRNA and mRNA expression levels in subchondral bone samples from three patients with ONFH and three patients with femoral neck fractures. Integrated bioinformatics analyses were conducted to identify lncRNAs associated with ONFH development and their potential functions and signaling pathways. A co-expression network was constructed based on the gene time-series expression data in GSE113253. After selecting lncRNA GAS5 as a novel biomarker for ONFH, bone marrow mesenchymal stem cell (BMSC) osteogenic differentiation assays were performed to verify the association between lncRNA GAS5 and osteogenic differentiation. Alkaline phosphatase (ALP) staining and quantitative reverse transcription polymerase chain reaction (RT-qPCR) were used to measure the osteogenic phenotype and lncRNA GAS5 expression. Finally, for further validation, ONFH rat models were established, and lncRNA GAS5 expression in subchondral bone was detected by RT-qPCR.

Results: We identified 126 and 959 differentially expressed lncRNAs and genes, respectively. lncRNA GAS5 expression level was significantly downregulated in patients with ONFH compared to the control group patients. The BMSC osteogenic differentiation assays showed that ALP activity increased gradually from days 3 to 7, while the lncRNA GAS5 expression level was significantly upregulated in the osteogenic differentiation induction groups. Furthermore, in vivo experiments suggested that the bone volume/tissue volume value and trabecular thickness significantly decreased in the ONFH rat model group compared to the control group, whereas the trabecular space significantly increased in the ONFH group compared to the control group. In addition, the lncRNA GAS5 expression level significantly decreased in the ONFH rat model group.

Conclusion: The lncRNA GAS5 expression level was highly associated with BMSC osteogenic differentiation and was significantly downregulated in both the subchondral trabecular bone tissue of ONFH patients and ONFH rat models. Therefore, lncRNA GAS5 can serve as an ONFH osteogenic biomarker to provide an effective target for early diagnosis and molecular therapy of ONFH.

Integrative analysis reveals a lineage-specific circular RNA landscape for adipo-osteogenesis of human mesenchymal stem cells

Background

The balance between osteogenesis and adipogenesis of mesenchymal stem cells (MSCs) is critical to skeletal development and diseases. As a research hotspot, circular RNAs (circRNAs) have expanded our understanding of a hidden layer of the transcriptome. Yet, their roles during adipo-osteogenesis remain poorly described.

Methods

The identity of human MSCs derived from bone marrow and adipose were first determined by flow cytometry, cellular staining, and quantitative polymerase chain reaction (qPCR). Multi-strategic RNA-sequencing was performed using Poly A, RiboMinus and RiboMinus/RNase R methods. Integrative analysis was performed to identify lineage-specific expressed circRNAs. The structural and expressional characteristics were identified by Sanger sequencing and qPCR, respectively. The regulatory effects of adipogenesis-specific circ-CRLF1 were confirmed using siRNA transcfection and qPCR.

Results

We generated a whole transcriptome map during adipo-osteogenesis based on 10 Poly A, 20 RiboMinus and 20 RiboMinus/ RNase R datasets. A total of 31,326 circRNAs were identified and quantified from ~ 3.4 billion paired-end reads. Furthermore, the integrative analysis revealed that 1166 circRNA genes exhibited strong lineage-specific expression patterns. Their host genes were enriched in distinct biological functions, such as cell adhesion, cytokine signaling, and cell division. We randomly selected and validated the back-spliced junction sites and expression patterns of 12 lineage-specific circRNAs. Functional analysis indicated that circ-CRLF1 negatively regulated adipogenesis.

Conclusions

Our integrative analysis reveals an accurate and generally applicable lineage-specific circRNA landscape for adipo-osteogenesis of MSCs and provides a potential therapeutic target, circ-CRLF1, for the treatment of skeleton-related disease.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13287-022-02792-5.

Increased BMSC exosomal miR-140-3p alleviates bone degradation and promotes bone restoration by targeting Plxnb1 in diabetic rats

Background

Diabetes mellitus (DM) is considered to be an important factor for bone degeneration disorders such as bone defect nonunion, which is characterized by physical disability and tremendous economy cost to families and society. Exosomal miRNAs of BMSCs have been reported to participate in osteoblastogenesis and modulating bone formation. However, their impacts on the development of bone degeneration in DM are not yet known. The role of miRNAs in BMSCs exosomes on regulating hyperglycemia bone degeneration was investigated in the present study.

Results

The osteogenic potential in bone defect repair of exosomes derived from diabetes mellitus BMSCs derived exosomes (DM-Exos) were revealed to be lower than that in normal BMSCs derived exosomes (N-Exos) in vitro and in vivo. Here, we demonstrate that miR-140-3p level was significantly altered in exosomes derived from BMSCs, ADSCs and serum from DM rats. In in vitro experiments, upregulated miR-140-3p exosomes promoted DM BMSCs differentiation into osteoblasts. The effects were exerted by miR-140-3p targeting plxnb1, plexin B1 is the receptor of semaphoring 4D(Sema4D) that inhibited osteocytes differentiation, thereby promoting bone formation. In DM rats with bone defect, miR-140-3p upregulated exosomes were transplanted into injured bone and accelerated bone regeneration. Besides, miR-140-3p in the exosomes was transferred into BMSCs and osteoblasts and promoted bone regeneration by targeting the plexin B1/RohA/ROCK signaling pathway.

Conclusions

Normal-Exos and miR-140-3p overexpressed-Exos accelerated diabetic wound healing by promoting the osteoblastogenesis function of BMSCs through inhibition plexin B1 expression which is the receptor of Sema4D and the plexin B1/RhoA/ROCK pathway compared with diabetes mellitus-Exos. This offers a new insight and a new therapy for treating diabetic bone unhealing.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s12951-022-01267-2.