Single cell RNA sequencing reveals emergent notochord-derived cell subpopulations in the postnatal nucleus pulposus.

Intervertebral disc degeneration is a leading cause of chronic low back pain. Cell-based strategies that seek to treat disc degeneration by regenerating the central nucleus pulposus (NP) hold significant promise, but key challenges remain. One of these is the inability of therapeutic cells to effectively mimic the performance of native NP cells, which are unique amongst skeletal cell types in that they arise from the embryonic notochord. In this study, we use single cell RNA sequencing to demonstrate emergent heterogeneity amongst notochord-derived NP cells in the postnatal mouse disc. Specifically, we established the existence of progenitor and mature NP cells, corresponding to notochordal and chondrocyte-like cells, respectively. Mature NP cells exhibited significantly higher expression levels of extracellular matrix (ECM) genes including aggrecan, and collagens II and VI, along with elevated transforming growth factor-beta and phosphoinositide 3 kinase-protein kinase B signaling. Additionally, we identified Cd9 as a novel surface marker of mature NP cells, and demonstrated that these cells were localized to the NP periphery, increased in numbers with increasing postnatal age, and co-localized with emerging glycosaminoglycan-rich matrix. Finally, we used a goat model to show that Cd9+ NP cell numbers decrease with moderate severity disc degeneration, suggesting that these cells are associated with maintenance of the healthy NP ECM. Improved understanding of the developmental mechanisms underlying regulation of ECM deposition in the postnatal NP may inform improved regenerative strategies for disc degeneration and associated low back pain.

Activin A marks a novel progenitor cell population during fracture healing and reveals a therapeutic strategy.

Insufficient bone fracture repair represents a major clinical and societal burden and novel strategies are needed to address it. Our data reveal that the transforming growth factor-ß superfamily member Activin A became very abundant during mouse and human bone fracture healing but was minimally detectable in intact bones. Single-cell RNA-sequencing revealed that the Activin A-encoding gene Inhba was highly expressed in a unique, highly proliferative progenitor cell (PPC) population with a myofibroblast character that quickly emerged after fracture and represented the center of a developmental trajectory bifurcation producing cartilage and bone cells within callus. Systemic administration of neutralizing Activin A antibody inhibited bone healing. In contrast, a single recombinant Activin A implantation at fracture site in young and aged mice boosted: PPC numbers; phosphorylated SMAD2 signaling levels; and bone repair and mechanical properties in endochondral and intramembranous healing models. Activin A directly stimulated myofibroblastic differentiation, chondrogenesis and osteogenesis in periosteal mesenchymal progenitor culture. Our data identify a distinct population of Activin A-expressing PPCs central to fracture healing and establish Activin A as a potential new therapeutic tool.

Single Cell RNA Sequencing Reveals Emergent Notochord-Derived Cell Subpopulations in the Postnatal Nucleus Pulposus.

Intervertebral disc degeneration is a leading cause of chronic low back pain. Cell-based strategies that seek to treat disc degeneration by regenerating the central nucleus pulposus hold significant promise, but key challenges remain. One of these is the inability of therapeutic cells to effectively mimic the performance of native nucleus pulposus cells, which are unique amongst skeletal cell types in that they arise from the embryonic notochord. In this study we use single cell RNA sequencing to demonstrate emergent heterogeneity amongst notochord-derived nucleus pulposus cells in the postnatal mouse disc. Specifically, we established the existence of early and late stage nucleus pulposus cells, corresponding to notochordal progenitor and mature cells, respectively. Late stage cells exhibited significantly higher expression levels of extracellular matrix genes including aggrecan, and collagens II and VI, along with elevated TGF-ß and PI3K-Akt signaling. Additionally, we identified Cd9 as a novel surface marker of late stage nucleus pulposus cells, and demonstrated that these cells were localized to the nucleus pulposus periphery, increased in numbers with increasing postnatal age, and co-localized with emerging glycosaminoglycan-rich matrix. Finally, we used a goat model to show the Cd9+ nucleus pulposus cell numbers decrease with moderate severity disc degeneration, suggesting that these cells are associated with maintenance of the healthy nucleus pulposus extracellular matrix. Improved understanding of the developmental mechanisms underlying regulation of ECM deposition in the postnatal NP may inform improved regenerative strategies for disc degeneration and associated low back pain.

Csf1 from marrow adipogenic precursors is required for osteoclast formation and hematopoiesis in bone.

Colony-stimulating factor 1 (Csf1) is an essential growth factor for osteoclast progenitors and an important regulator for bone resorption. It remains elusive which mesenchymal cells synthesize Csf1 to stimulate osteoclastogenesis. We recently identified a novel mesenchymal cell population, marrow adipogenic lineage precursors (MALPs), in bone. Compared to other mesenchymal subpopulations, MALPs expressed Csf1 at a much higher level and this expression was further increased during aging. To investigate its role, we constructed MALP-deficient Csf1 CKO mice using AdipoqCre. These mice had increased femoral trabecular bone mass, but their cortical bone appeared normal. In comparison, depletion of Csf1 in the entire mesenchymal lineage using Prrx1Cre led to a more striking high bone mass phenotype, suggesting that additional mesenchymal subpopulations secrete Csf1. TRAP staining revealed diminished osteoclasts in the femoral secondary spongiosa region of Csf1 CKOAdipoq mice, but not at the chondral-osseous junction nor at the endosteal surface of cortical bone. Moreover, Csf1 CKOAdipoq mice were resistant to LPS-induced calvarial osteolysis. Bone marrow cellularity, hematopoietic progenitors, and macrophages were also reduced in these mice. Taken together, our studies demonstrate that MALPs synthesize Csf1 to control bone remodeling and hematopoiesis.

Study on Regional Strata Movement during Deep Mining of Erdos Coal Field and Its Control.

Disasters such as rock bursts and mine earthquakes became increasingly serious with the increase in mining depth in Erdos Coal Field and became serious problems that restrict high-strength continuous mining of coal mines. In this study, strata movement and energy polling distribution of ultrathick weak-bonding sandstone layers were controlled by the local filling−caving multi-faces coordinated mining technique, which was based on the analysis of subsidence and overlying structural characteristics in the Yingpanhao mining area. Moreover, the influencing factors and the control effect laws were investigated. Surface subsidence and energy polling distribution control effects of different mining modes were compared, which confirmed the superiority of local filling based on the main key stratum. According to the results, the maximum surface subsidence velocity of the first mining face was 1.24 mm/d, which indicates the presence of a logistic functional relationship between the mining degree and subsidence factors. When the mining degree was close to full mining, the practical surface subsidence was smaller than the corresponding logistic functional value. The largest influencing factor for the strata movement control effect of partial filling mining based on the main key stratum was the width of the caving face, followed by the filling ratio, section pillar width, and width of the filling face, successively. With respect to the influencing degree on the energy polling distribution of partial filling mining based on the main key stratum, the order followed as section pillar width > filling ratio > caving working face > width of backfilling working face. Additionally, the comparative analysis from the perspectives of control effect, resource utilization, and cost-effectiveness demonstrated that partial filling mining based on the main key stratum was one of the techniques with high cost-effectiveness in controlling strata movement and relieving rock bursts, mining earthquakes, and subsidence disasters.

Activating EGFR Signaling Attenuates Osteoarthritis Development Following Loading Injury in Mice.

Posttraumatic osteoarthritis (PTOA) results in joint pain, loss of joint function, and impaired quality of daily life in patients with limited treatment options. We previously demonstrated that epidermal growth factor receptor (EGFR) signaling is essential for maintaining chondroprogenitors during articular cartilage development and homeostasis. Here, we used a nonsurgical, loading-induced PTOA mouse model to investigate the protective action of EGFR signaling. A single bout of cyclic tibial loading at a peak force of 6 N injured cartilage at the posterior aspect of lateral femoral condyle. Similar loading at a peak force of 9 N ruptured the anterior cruciate ligament, causing additional cartilage damage at the medial compartment and ectopic cartilage formation in meniscus and synovium. Constitutively overexpression of an EGFR ligand, heparin binding EGF-like growth factor (HBEGF), in chondrocytes significantly reduced cartilage injury length, synovitis, and pain after 6 N loading and mitigated medial side cartilage damage and ectopic cartilage formation after 9 N loading. Mechanistically, overactivation of EGFR signaling protected chondrocytes from loading-induced apoptosis and loss of proliferative ability and lubricant synthesis. Overexpressing HBEGF in adult cartilage starting right before 6 N loading had similar beneficial effects. In contrast, inactivating EGFR in adult cartilage led to accelerated PTOA progression with elevated cartilage Mankin score and synovitis score and increased ectopic cartilage formation. As a therapeutic approach, we constructed a nanoparticle conjugated with the EGFR ligand TGFα. Intra-articular injections of this nanoconstruct once every 3 weeks for 12 weeks partially mitigated PTOA symptoms in cartilage and synovium after 6 N loading. Our findings demonstrate the anabolic actions of EGFR signaling in maintaining articular cartilage during PTOA development and shed light on developing a novel nanomedicine for PTOA. © 2022 American Society for Bone and Mineral Research (ASBMR).

Transient expansion and myofibroblast conversion of adipogenic lineage precursors mediate bone marrow repair after radiation.

Radiation causes a collapse of bone marrow cells and elimination of microvasculature. To understand how bone marrow recovers after radiation, we focused on mesenchymal lineage cells that provide a supportive microenvironment for hematopoiesis and angiogenesis in bone. We recently discovered a nonproliferative subpopulation of marrow adipogenic lineage precursors (MALPs) that express adipogenic markers with no lipid accumulation. Single-cell transcriptomic analysis revealed that MALPs acquire proliferation and myofibroblast features shortly after radiation. Using an adipocyte-specific Adipoq-Cre, we validated that MALPs rapidly and transiently expanded at day 3 after radiation, coinciding with marrow vessel dilation and diminished marrow cellularity. Concurrently, MALPs lost most of their cell processes, became more elongated, and highly expressed myofibroblast-related genes. Radiation activated mTOR signaling in MALPs that is essential for their myofibroblast conversion and subsequent bone marrow recovery at day 14. Ablation of MALPs blocked the recovery of bone marrow vasculature and cellularity, including hematopoietic stem and progenitors. Moreover, VEGFa deficiency in MALPs delayed bone marrow recovery after radiation. Taken together, our research demonstrates a critical role of MALPs in mediating bone marrow repair after radiation injury and sheds light on a cellular target for treating marrow suppression after radiotherapy.

Superoxide dismutase-loaded porous polymersomes as highly efficient antioxidant nanoparticles targeting synovium for osteoarthritis therapy.

Oxidative stress and the reactive oxygen species (ROS) have important roles in osteoarthritis (OA) development and progression. Scavenging ROS by exogenous antioxidant enzymes could be a promising approach for OA treatment. However, the direct use of antioxidant enzymes, such as superoxide dismutase (SOD), is challenging due to a lack of effective drug delivery system to knee joints. This study utilized a highly efficient antioxidative nanoparticle based on SOD-loaded porous polymersome nanoparticles (SOD-NPs) for delivery of SOD to mouse knee joints. The resultant SOD-NPs had prolonged mouse joint retention time with predominant accumulation in synovium but not in articular cartilage. Examining human synovial explants revealed that SOD-NPs minimize oxidative damages induced by OA-like insults. Intra-articular injections of SOD-NPs in mice receiving OA surgery were effective in attenuating OA initiation and preventing its further progression. Mechanistically, SOD-NPs reduced ROS production and the synthesis of catabolic proteases in both articular cartilage and synovium. Hence, our work demonstrates the therapeutic potential of SOD-NPs and indicate that targeting synovium holds a great promise for OA therapy.

EGFR Signaling Is Required for Maintaining Adult Cartilage Homeostasis and Attenuating Osteoarthritis Progression.

The uppermost superficial zone of articular cartilage is the first line of defense against the initiation of osteoarthritis (OA). We previously used Col2-Cre to demonstrate that epidermal growth factor receptor (EGFR), a tyrosine kinase receptor, plays an essential role in maintaining superficial chondrocytes during articular cartilage development. Here, we showed that EGFR activity in the articular cartilage decreased as mice age. In mouse and human OA samples, EGFR activity was initially reduced at the superficial layer and then resurged in cell clusters within the middle and deep zone in late OA. To investigate the role of EGFR signaling in postnatal and adult cartilage, we constructed an inducible mouse model with cartilage-specific EGFR inactivation (Aggrecan-CreER EgfrWa5/flox , Egfr iCKO). EdU incorporation revealed that postnatal Egfr iCKO mice contained fewer slow-cycling cells than controls. EGFR deficiency induced at 3 months of age reduced cartilage thickness and diminished superficial chondrocytes, in parallel to alterations in lubricin production, cell proliferation, and survival. Furthermore, male Egfr iCKO mice developed much more severe OA phenotypes, including cartilage erosion, subchondral bone plate thickening, cartilage degeneration at the lateral site, and mechanical allodynia, after receiving destabilization of the medial meniscus (DMM) surgery. Similar OA phenotypes were also observed in female iCKO mice. Moreover, tamoxifen injections of iCKO mice at 1 month post-surgery accelerated OA development 2 months later. In summary, our data demonstrated that chondrogenic EGFR signaling maintains postnatal slow-cycling cells and plays a critical role in adult cartilage homeostasis and OA progression. © 2022 American Society for Bone and Mineral Research (ASBMR).

An ECHO of Cartilage: In Silico Prediction of Combinatorial Treatments to Switch Between Transient and Permanent Cartilage Phenotypes With Ex Vivo Validation.

A fundamental question in cartilage biology is: what determines the switch between permanent cartilage found in the articular joints and transient hypertrophic cartilage that functions as a template for bone? This switch is observed both in a subset of OA patients that develop osteophytes, as well as in cell-based tissue engineering strategies for joint repair. A thorough understanding of the mechanisms regulating cell fate provides opportunities for treatment of cartilage disease and tissue engineering strategies. The objective of this study was to understand the mechanisms that regulate the switch between permanent and transient cartilage using a computational model of chondrocytes, ECHO. To investigate large signaling networks that regulate cell fate decisions, we developed the software tool ANIMO, Analysis of Networks with interactive Modeling. In ANIMO, we generated an activity network integrating 7 signal transduction pathways resulting in a network containing over 50 proteins with 200 interactions. We called this model ECHO, for executable chondrocyte. Previously, we showed that ECHO could be used to characterize mechanisms of cell fate decisions. ECHO was first developed based on a Boolean model of growth plate. Here, we show how the growth plate Boolean model was translated to ANIMO and how we adapted the topology and parameters to generate an articular cartilage model. In ANIMO, many combinations of overactivation/knockout were tested that result in a switch between permanent cartilage (SOX9+) and transient, hypertrophic cartilage (RUNX2+). We used model checking to prioritize combination treatments for wet-lab validation. Three combinatorial treatments were chosen and tested on metatarsals from 1-day old rat pups that were treated for 6 days. We found that a combination of IGF1 with inhibition of ERK1/2 had a positive effect on cartilage formation and growth, whereas activation of DLX5 combined with inhibition of PKA had a negative effect on cartilage formation and growth and resulted in increased cartilage hypertrophy. We show that our model describes cartilage formation, and that model checking can aid in choosing and prioritizing combinatorial treatments that interfere with normal cartilage development. Here we show that combinatorial treatments induce changes in the zonal distribution of cartilage, indication possible switches in cell fate. This indicates that simulations in ECHO aid in describing pathologies in which switches between cell fates are observed, such as OA.

The critical role of Hedgehog-responsive mesenchymal progenitors in meniscus development and injury repair.

Meniscal tears are associated with a high risk of osteoarthritis but currently have no disease-modifying therapies. Using a Gli1 reporter line, we found that Gli1+ cells contribute to the development of meniscus horns from 2 weeks of age. In adult mice, Gli1+ cells resided at the superficial layer of meniscus and expressed known mesenchymal progenitor markers. In culture, meniscal Gli1+ cells possessed high progenitor activities under the control of Hh signal. Meniscus injury at the anterior horn induced a quick expansion of Gli1-lineage cells. Normally, meniscal tissue healed slowly, leading to cartilage degeneration. Ablation of Gli1+ cells further hindered this repair process. Strikingly, intra-articular injection of Gli1+ meniscal cells or an Hh agonist right after injury accelerated the bridging of the interrupted ends and attenuated signs of osteoarthritis. Taken together, our work identified a novel progenitor population in meniscus and proposes a new treatment for repairing injured meniscus and preventing osteoarthritis.

Marrow adipogenic lineage precursor: A new cellular component of marrow adipose tissue.

Bone marrow mesenchymal stromal cells are a highly heterogenic cell population containing mesenchymal stem cells as well as other cell types. With the advance of single cell transcriptome analysis, several recent reports identified a prominent subpopulation of mesenchymal stromal cells that specifically express adipocyte markers but do not contain lipid droplets. We name this cell type marrow adipogenic lineage precursor, MALP, and consider it as a major cellular component of marrow adipose tissue. Here, we review the discovery of MALPs and summarize their unique features and regulatory roles in bone. We further discuss how these findings advance our understanding of bone remodeling, mesenchymal niche regulation of hematopoiesis, and marrow vasculature maintenance.

Gli1 Defines a Subset of Fibro-adipogenic Progenitors that Promote Skeletal Muscle Regeneration With Less Fat Accumulation.

Skeletal muscle has remarkable regenerative ability after injury. Mesenchymal fibro-adipogenic progenitors (FAPs) are necessary, active participants during this repair process, but the molecular signatures of these cells and their functional relevance remain largely unexplored. Here, using a lineage tracing mouse model (Gli1-CreER Tomato), we demonstrate that Gli1 marks a small subset of muscle-resident FAPs with elevated Hedgehog (Hh) signaling. Upon notexin muscle injury, these cells preferentially and rapidly expanded within FAPs. Ablation of Gli1+ cells using a DTA mouse model drastically reduced fibroblastic colony-forming unit (CFU-F) colonies generated by muscle cells and impaired muscle repair at 28 days. Pharmacologic manipulation revealed that Gli1+ FAPs rely on Hh signaling to increase the size of regenerating myofiber. Sorted Gli1+ FAPs displayed superior clonogenicity and reduced adipogenic differentiation ability in culture compared to sorted Gli1- FAPs. In a glycerol injury model, Gli1+ FAPs were less likely to give rise to muscle adipocytes compared to other FAPs. Further cell ablation and Hh activator/inhibitor treatments demonstrated their dual actions in enhancing myogenesis and reducing adipogenesis after injury. Examining single-cell RNA-sequencing dataset of FAPs from normal mice indicated that Gli1+ FAPs with increased Hh signaling provide trophic signals to myogenic cells while restrict their own adipogenic differentiation. Collectively, our findings identified a subpopulation of FAPs that play an essential role in skeletal muscle repair. © 2021 American Society for Bone and Mineral Research (ASBMR).

Targeting cartilage EGFR pathway for osteoarthritis treatment.

Osteoarthritis (OA) is a widespread joint disease for which there are no disease-modifying treatments. Previously, we found that mice with cartilage-specific epidermal growth factor receptor (EGFR) deficiency developed accelerated knee OA. To test whether the EGFR pathway can be targeted as a potential OA therapy, we constructed two cartilage-specific EGFR overactivation models in mice by overexpressing heparin binding EGF-like growth factor (HBEGF), an EGFR ligand. Compared to wild type, Col2-Cre HBEGF-overexpressing mice had persistently enlarged articular cartilage from adolescence, due to an expanded pool of chondroprogenitors with elevated proliferation ability, survival rate, and lubricant production. Adult Col2-Cre HBEGF-overexpressing mice and Aggrecan-CreER HBEGF-overexpressing mice were resistant to cartilage degeneration and other signs of OA after surgical destabilization of the medial meniscus (DMM). Treating mice with gefitinib, an EGFR inhibitor, abolished the protective action against OA in HBEGF-overexpressing mice. Polymeric micellar nanoparticles (NPs) conjugated with transforming growth factor-α (TGFα), a potent EGFR ligand, were stable and nontoxic and had long joint retention, high cartilage uptake, and penetration capabilities. Intra-articular delivery of TGFα-NPs effectively attenuated surgery-induced OA cartilage degeneration, subchondral bone plate sclerosis, and joint pain. Genetic or pharmacologic activation of EGFR revealed no obvious side effects in knee joints and major vital organs in mice. Together, our studies demonstrate the feasibility of using nanotechnology to target EGFR signaling for OA treatment.

The Expressions of Dickkopf-Related Protein 1 and Frizzled-Related Protein Are Negatively Correlated to Local Inflammation and Osteoarthritis Severity.

OBJECTIVE: To investigate the presence of WNT antagonists Dickkopf-related protein 1 (DKK1), Frizzled-related protein (FRZB) and BMP antagonist Gremlin 1 (GREM1) in synovial fluid (SF) and serum, respectively, from end-stage knee osteoarthritis (OA) patients, and correlate their expression with other markers of OA. DESIGN: In a cross-sectional study, SF and serum were collected from OA patients (n = 132). The concentrations of DKK1, FRZB and GREM1 in SF and serum were determined using immunoassays. Correlation measurements were performed between groups and previously assessed disease markers, such as synovium nitric oxide (NO), inerleukin-1ß (IL1ß), tumor necrosis factor-α (TNFα), and prostaglandin E2 (PGE2). RESULTS: The OA patients with the celecoxib treatment till surgery have higher median SF FRZB values compared with the control (no treatment); the celecoxib 3-days before surgery stopped treatment group has higher median serum FRZB values than the control and the naproxen treatment group. The combinational analysis of SF DKK1 and SF FRZB negatively correlated with macroscopic cartilage scores and histological synovium scores in OA patients. The expression of DKK1 and FRZB in SF showed the same expression trend as their expression in serum. Furthermore, the SF concentration of DKK1 was positively correlated with FRZB in both SF and serum. In contrast, it was negatively correlated with synovium NO and IL1ß. SF FRZB was negatively correlated with synovium NO, IL1ß, cartilage PGE2, and age. CONCLUSIONS: Our findings suggest DKK1 and FRZB were negatively correlated with OA severity and multiple pro-inflammatory cytokines. Our data indicate that DKK1 and FRZB can be joint disease-specific biomarkers.

A Novel Enzymatic Digestion Approach for Isolation and Culture of Rodent Bone Marrow Mesenchymal Progenitors.

Bone marrow mesenchymal stem cells (MSCs) are promising therapeutic tools for tissue repair and treatment of a number of human diseases. As a result, there is substantial interest in characterizing and expanding these cells to uncover their therapeutic potential. Bone marrow mesenchymal progenitors, containing both MSCs and their proliferative progeny, are commonly isolated from the central region of rodent long bones. However, challenges exist in expanding these central mesenchymal progenitors in culture. We have designed an enzymatic digestion protocol to isolate mesenchymal progenitors within rodent long bones that resides close to the bone surface, which we termed endosteal mesenchymal progenitors. These cells are more metabolically active and more responsive to external stimuli compared to central mesenchymal progenitors. Therefore, they represent a biologically important target for MSC research. This chapter describes the approach in detail how to isolate and culture endosteal mesenchymal progenitors as well as their central counterparts from rodent long bones.

Bone marrow adipogenic lineage precursors promote osteoclastogenesis in bone remodeling and pathologic bone loss.

Bone is maintained by coupled activities of bone-forming osteoblasts/osteocytes and bone-resorbing osteoclasts. Alterations in this relationship can lead to pathologic bone loss such as osteoporosis. It is well known that osteogenic cells support osteoclastogenesis via production of RANKL. Interestingly, our recently identified bone marrow mesenchymal cell population-marrow adipogenic lineage precursors (MALPs) that form a multidimensional cell network in bone-was computationally demonstrated to be the most interactive with monocyte-macrophage lineage cells through high and specific expression of several osteoclast regulatory factors, including RANKL. Using an adipocyte-specific Adipoq-Cre to label MALPs, we demonstrated that mice with RANKL deficiency in MALPs have a drastic increase in trabecular bone mass in long bones and vertebrae starting from 1 month of age, while their cortical bone appears normal. This phenotype was accompanied by diminished osteoclast number and attenuated bone formation at the trabecular bone surface. Reduced RANKL signaling in calvarial MALPs abolished osteolytic lesions after LPS injections. Furthermore, in ovariectomized mice, elevated bone resorption was partially attenuated by RANKL deficiency in MALPs. In summary, our studies identified MALPs as a critical player in controlling bone remodeling during normal bone metabolism and pathological bone loss in a RANKL-dependent fashion.

Single cell transcriptomics identifies a unique adipose lineage cell population that regulates bone marrow environment.

Bone marrow mesenchymal lineage cells are a heterogeneous cell population involved in bone homeostasis and diseases such as osteoporosis. While it is long postulated that they originate from mesenchymal stem cells, the true identity of progenitors and their in vivo bifurcated differentiation routes into osteoblasts and adipocytes remain poorly understood. Here, by employing large scale single cell transcriptome analysis, we computationally defined mesenchymal progenitors at different stages and delineated their bi-lineage differentiation paths in young, adult and aging mice. One identified subpopulation is a unique cell type that expresses adipocyte markers but contains no lipid droplets. As non-proliferative precursors for adipocytes, they exist abundantly as pericytes and stromal cells that form a ubiquitous 3D network inside the marrow cavity. Functionally they play critical roles in maintaining marrow vasculature and suppressing bone formation. Therefore, we name them marrow adipogenic lineage precursors (MALPs) and conclude that they are a newly identified component of marrow adipose tissue.

ECHO, the executable CHOndrocyte: A computational model to study articular chondrocytes in health and disease.

Computational modeling can be used to investigate complex signaling networks in biology. However, most modeling tools are not suitable for molecular cell biologists with little background in mathematics. We have built a visual-based modeling tool for the investigation of dynamic networks. Here, we describe the development of computational models of cartilage development and osteoarthritis, in which a panel of relevant signaling pathways are integrated. In silico experiments give insight in the role of each of the pathway components and reveal which perturbations may deregulate the basal healthy state of cells and tissues. We used a previously developed computational modeling tool Analysis of Networks with Interactive Modeling (ANIMO) to generate an activity network integrating 7 signal transduction pathways resulting in a network containing over 50 nodes and 200 interactions. We performed in silico experiments to characterize molecular mechanisms of cell fate decisions. The model was used to mimic biological scenarios during cell differentiation using RNA-sequencing data of a variety of stem cell sources as input. In a case-study, we wet-lab-tested the model-derived hypothesis that expression of DKK1 (Dickkopf-1) and FRZB (Frizzled related protein, WNT antagonists) and GREM1 (gremlin 1, BMP antagonist) prevents IL1ß (Interleukin 1 beta)-induced MMP (matrix metalloproteinase) expression, thereby preventing cartilage degeneration, at least in the short term. We found that a combination of DKK1, FRZB and GREM1 may play a role in modulating the effects of IL1ß induced inflammation in human primary chondrocytes.

Periarticular Mesenchymal Progenitors Initiate and Contribute to Secondary Ossification Center Formation During Mouse Long Bone Development.

Long bone development involves the embryonic formation of a primary ossification center (POC) in the incipient diaphysis followed by postnatal development of a secondary ossification center (SOC) at each epiphysis. Studies have elucidated major basic mechanisms of POC development, but relatively little is known about SOC development. To gain insights into SOC formation, we used Col2-Cre Rosa-tdTomato (Col2/Tomato) reporter mice and found that their periarticular region contained numerous Tomato-positive lineage cells expressing much higher Tomato fluorescence (termed TomatoH ) than underlying epiphyseal chondrocytes (termed TomatoL ). With time, the TomatoH cells became evident at the SOC invagination site and cartilage canal, increased in number in the expanding SOC, and were present as mesenchymal lineage cells in the subchondral bone. These data were verified in two mouse lineage tracing models, Col2-CreER Rosa-tdTomato and Gli1-CreER Rosa-tdTomato. In vitro tests showed that the periarticular TomatoH cells from Col2/Tomato mice contained mesenchymal progenitors with multidifferentiation abilities. During canal initiation, the cells expressed vascular endothelial growth factor (VEGF) and migrated into epiphyseal cartilage ahead of individual or clusters of endothelial cells, suggesting a unique role in promoting vasculogenesis. Later during SOC expansion, chondrocytes in epiphyseal cartilage expressed VEGF, and angiogenic blood vessels preceded TomatoH cells. Gene expression analyses of microdissected samples revealed upregulation of MMPs in periarticular cells at the invagination site and suggested potential roles for novel kinase and growth factor signaling pathways in regulating SOC canal initiation. In summary, our data indicate that the periarticular region surrounding epiphyseal cartilage contains mesenchymal progenitors that initiate SOC development and form subchondral bone. Stem Cells 2019;37:677-689.