Osteoclasts recycle via osteomorphs during RANKL-stimulated bone resorption.

Osteoclasts are large multinucleated bone-resorbing cells formed by the fusion of monocyte/macrophage-derived precursors that are thought to undergo apoptosis once resorption is complete. Here, by intravital imaging, we reveal that RANKL-stimulated osteoclasts have an alternative cell fate in which they fission into daughter cells called osteomorphs. Inhibiting RANKL blocked this cellular recycling and resulted in osteomorph accumulation. Single-cell RNA sequencing showed that osteomorphs are transcriptionally distinct from osteoclasts and macrophages and express a number of non-canonical osteoclast genes that are associated with structural and functional bone phenotypes when deleted in mice. Furthermore, genetic variation in human orthologs of osteomorph genes causes monogenic skeletal disorders and associates with bone mineral density, a polygenetic skeletal trait. Thus, osteoclasts recycle via osteomorphs, a cell type involved in the regulation of bone resorption that may be targeted for the treatment of skeletal diseases.

ENTRAIN: integrating trajectory inference and gene regulatory networks with spatial data to co-localize the receptor-ligand interactions that specify cell fate.

MOTIVATION: Cell fate is commonly studied by profiling the gene expression of single cells to infer developmental trajectories based on expression similarity, RNA velocity, or statistical mechanical properties. However, current approaches do not recover microenvironmental signals from the cellular niche that drive a differentiation trajectory. RESULTS: We resolve this with environment-aware trajectory inference (ENTRAIN), a computational method that integrates trajectory inference methods with ligand-receptor pair gene regulatory networks to identify extracellular signals and evaluate their relative contribution towards a differentiation trajectory. The output from ENTRAIN can be superimposed on spatial data to co-localize cells and molecules in space and time to map cell fate potentials to cell-cell interactions. We validate and benchmark our approach on single-cell bone marrow and spatially resolved embryonic neurogenesis datasets to identify known and novel environmental drivers of cellular differentiation. AVAILABILITY AND IMPLEMENTATION: ENTRAIN is available as a public package at https://github.com/theimagelab/entrain and can be used on both single-cell and spatially resolved datasets.

Single-Cell RNA Sequencing: Unravelling the Bone One Cell at a Time.

PURPOSE OF REVIEW: Bone is a complex tissue populated by a highly heterogeneous mix of cell types in different compartments. The endosteal compartment is a key site for bone remodelling and provides a supportive microenvironment to harbour haematopoietic and mesenchymal stem cells, as well as cancer cells that grow in bone. The purpose of this review is to summarize recent findings of studies in bone using single-cell RNA sequencing and emergent spatial RNA sequencing to describe different bone-resident cell types and their molecular programs. RECENT FINDINGS: Single-cell RNA sequencing identified novel and transcriptionally distinct cell clusters within different bone cell lineages, including MSCs, osteoblasts, chondrocytes, fibroblasts, osteoclasts and cells of the vasculature. Spatial transcriptomics methods provide information on the localization of the different cell populations. Single-cell transcriptomics provided valuable insights into long-standing knowledge gaps in the cellular heterogeneity of bone-resident cells in unprecedented detail, paving the way for studies to further investigate the different cell populations and to develop cell-based therapies for bone diseases.

A niche-dependent myeloid transcriptome signature defines dormant myeloma cells.

The era of targeted therapies has seen significant improvements in depth of response, progression-free survival, and overall survival for patients with multiple myeloma. Despite these improvements in clinical outcome, patients inevitably relapse and require further treatment. Drug-resistant dormant myeloma cells that reside in specific niches within the skeleton are considered a basis of disease relapse but remain elusive and difficult to study. Here, we developed a method to sequence the transcriptome of individual dormant myeloma cells from the bones of tumor-bearing mice. Our analyses show that dormant myeloma cells express a distinct transcriptome signature enriched for immune genes and, unexpectedly, genes associated with myeloid cell differentiation. These genes were switched on by coculture with osteoblastic cells. Targeting AXL, a gene highly expressed by dormant cells, using small-molecule inhibitors released cells from dormancy and promoted their proliferation. Analysis of the expression of AXL and coregulated genes in human cohorts showed that healthy human controls and patients with monoclonal gammopathy of uncertain significance expressed higher levels of the dormancy signature genes than patients with multiple myeloma. Furthermore, in patients with multiple myeloma, the expression of this myeloid transcriptome signature translated into a twofold increase in overall survival, indicating that this dormancy signature may be a marker of disease progression. Thus, engagement of myeloma cells with the osteoblastic niche induces expression of a suite of myeloid genes that predicts disease progression and that comprises potential drug targets to eradicate dormant myeloma cells.

Mevalonate kinase deficiency leads to decreased prenylation of Rab GTPases.

Mevalonate kinase deficiency (MKD) is caused by mutations in a key enzyme of the mevalonate-cholesterol biosynthesis pathway, leading to recurrent autoinflammatory disease characterised by enhanced release of interleukin-1ß (IL-1ß). It is currently believed that the inflammatory phenotype of MKD is triggered by temperature-sensitive loss of mevalonate kinase activity and reduced biosynthesis of isoprenoid lipids required for the prenylation of small GTPase proteins. However, previous studies have not clearly shown any change in protein prenylation in patient cells under normal conditions. With lymphoblast cell lines from two compound heterozygous MKD patients, we used a highly sensitive in vitro prenylation assay, together with quantitative mass spectrometry, to reveal a subtle accumulation of unprenylated Rab GTPases in cells cultured for 3 days or more at 40 °C compared with 37 °C. This included a 200% increase in unprenylated Rab7A, Rab14 and Rab1A. Inhibition of sterol regulatory element-binding protein (SREBP) activation by fatostatin led to more pronounced accumulation of unprenylated Rab proteins in MKD cells but not parent cells, suggesting that cultured MKD cells may partially overcome the loss of isoprenoid lipids by SREBP-mediated upregulation of enzymes required for isoprenoid biosynthesis. Furthermore, while inhibition of Rho/Rac/Rap prenylation promoted the release of IL-1ß, specific inhibition of Rab prenylation by NE10790 had no effect in human peripheral blood mononuclear cells or human THP-1 monocytic cells. These studies demonstrate for the first time that mutations in mevalonate kinase can lead to a mild, temperature-induced defect in the prenylation of small GTPases, but that loss of prenylated Rab GTPases is not the cause of enhanced IL-1ß release in MKD.

A pilot study to compare the detection of HPV-16 biomarkers in salivary oral rinses with tumour p16(INK4a) expression in head and neck squamous cell carcinoma patients.

BACKGROUND: Human papilloma virus-16 (HPV-16) infection is a major risk factor for a subset of head and neck squamous cell carcinoma (HNSCC), in particular oropharyngeal squamous cell carcinoma (OPSCC). Current techniques for assessing the HPV-16 status in HNSCC include the detection of HPV-16 DNA and p16(INK4a) expression in tumor tissues. When tumors originate from hidden anatomical sites, this method can be challenging. A non-invasive and cost-effective alternative to biopsy is therefore desirable for HPV-16 detection especially within a community setting to screen at-risk individuals. METHODS: The present study compared detection of HPV-16 DNA and RNA in salivary oral rinses with tumor p16(INK4a) status, in 82 HNSCC patients using end-point and quantitative polymerase chain reaction (PCR). RESULTS: Of 42 patients with p16(INK4a)-positive tumours, 39 (sensitivity = 92.9 %, PPV = 100 % and NPV = 93 %) had oral rinse samples with detectable HPV-16 DNA, using end-point and quantitative PCR. No HPV-16 DNA was detected in oral rinse samples from 40 patients with p16(INK4a) negative tumours, yielding a test specificity of 100 %. For patients with p16(INK4a) positive tumours, HPV-16 mRNA was detected using end-point reverse transcription PCR (RT-PCR) in 24/40 (sensitivity = 60 %, PPV = 100 % and NPV = 71 %), and using quantitative RT-PCR in 22/40 (sensitivity = 55 %, PPV = 100 % and NPV = 69 %). No HPV-16 mRNA was detected in oral rinse samples from the p16(INK4a)-negative patients, yielding a specificity of 100 %. CONCLUSIONS: We demonstrate that the detection of HPV-16 DNA in salivary oral rinse is indicative of HPV status in HNSCC patients and can potentially be used as a diagnostic tool in addition to the current methods.

Molecular stress-inducing compounds increase osteoclast formation in a heat shock factor 1 protein-dependent manner.

Many anticancer therapeutic agents cause bone loss, which increases the risk of fractures that severely reduce quality of life. Thus, in drug development, it is critical to identify and understand such effects. Anticancer therapeutic and HSP90 inhibitor 17-(allylamino)-17-demethoxygeldanamycin (17-AAG) causes bone loss by increasing osteoclast formation, but the mechanism underlying this is not understood. 17-AAG activates heat shock factor 1 (Hsf1), the master transcriptional regulator of heat shock/cell stress responses, which may be involved in this negative action of 17-AAG upon bone. Using mouse bone marrow and RAW264.7 osteoclast differentiation models we found that HSP90 inhibitors that induced a heat shock response also enhanced osteoclast formation, whereas HSP90 inhibitors that did not (including coumermycin A1 and novobiocin) did not affect osteoclast formation. Pharmacological inhibition or shRNAmir knockdown of Hsf1 in RAW264.7 cells as well as the use of Hsf1 null mouse bone marrow cells demonstrated that 17-AAG-enhanced osteoclast formation was Hsf1-dependent. Moreover, ectopic overexpression of Hsf1 enhanced 17-AAG effects upon osteoclast formation. Consistent with these findings, protein levels of the essential osteoclast transcription factor microphthalmia-associated transcription factor were increased by 17-AAG in an Hsf1-dependent manner. In addition to HSP90 inhibitors, we also identified that other agents that induced cellular stress, such as ethanol, doxorubicin, and methotrexate, also directly increased osteoclast formation, potentially in an Hsf1-dependent manner. These results, therefore, indicate that cellular stress can enhance osteoclast differentiation via Hsf1-dependent mechanisms and may significantly contribute to pathological and therapeutic related bone loss.

Heat-shock factor 1 both positively and negatively affects cellular clonogenic growth depending on p53 status.

HSF1 (heat-shock factor 1) is the master regulator of the heat-shock response; however, it is also activated by cancer-associated stresses and supports cellular transformation and cancer progression. We examined the role of HSF1 in relation to cancer cell clonogenicity, an important attribute of cancer cells. Ectopic expression or HSF1 knockdown demonstrated that HSF1 positively regulated cancer cell clonogenic growth. Furthermore, knockdown of mutant p53 indicated that HSF1 actions were mediated via a mutant p53-dependent mechanism. To examine this relationship more specifically, we ectopically co-expressed mutant p53(R273H) and HSF1 in the human mammary epithelial cell line MCF10A. Surprisingly, within this cellular context, HSF1 inhibited clonogenicity. However, upon specific knockdown of endogenous wild-type p53, leaving mutant p53(R273H) expression intact, HSF1 was observed to greatly enhance clonogenic growth of the cells, indicating that HSF1 suppressed clonogenicity via wild-type p53. To confirm this we ectopically expressed HSF1 in non-transformed and H-Ras(V12)-transformed MCF10A cells. As expected, HSF1 significantly reduced clonogenicity, altering wild-type p53 target gene expression levels consistent with a role of HSF1 increasing wild-type p53 activity. In support of this finding, knockdown of wild-type p53 negated the inhibitory effects of HSF1 expression. We thus show that HSF1 can affect clonogenic growth in a p53 context-dependent manner, and can act via both mutant and wild-type p53 to bring about divergent effects upon clonogenicity. These findings have important implications for our understanding of HSF1's divergent roles in cancer cell growth and survival as well as its disparate effect on mutant and wild-type p53.

HSP90 inhibitors enhance differentiation and MITF (microphthalmia transcription factor) activity in osteoclast progenitors.

The HSP90 (heat-shock protein 90) inhibitor 17-AAG (17-allylamino-demethoxygeldanamycin) increases osteoclast formation both in vitro and in vivo, an action that can enhance cancer invasion and growth in the bone microenvironment. The cellular mechanisms through which 17-AAG exerts this action are not understood. Thus we sought to clarify the actions of 17-AAG on osteoclasts and determine whether other HSP90 inhibitors had similar properties. We determined that 17-AAG and the structurally unrelated HSP90 inhibitors CCT018159 and NVP-AUY922 dose-dependently increased RANKL [receptor activator of NF-κB (nuclear factor κB) ligand]-stimulated osteoclastogenesis in mouse bone marrow and pre-osteoclastic RAW264.7 cell cultures. Moreover, 17-AAG also enhanced RANKL- and TNF (tumour necrosis factor)-elicited osteoclastogenesis, but did not affect RANKL-induced osteoclast survival, suggesting that only differentiation mechanisms are targeted. 17-AAG affected the later stages of progenitor maturation (after 3 days of incubation), whereas the osteoclast formation enhancer TGFß (transforming growth factor ß) acted prior to this, suggesting different mechanisms of action. In studies of RANKL-elicited intracellular signalling, 17-AAG treatment did not increase c-Fos or NFAT (nuclear factor of activated T-cells) c1 protein levels nor did 17-AAG increase activity in luciferase-based NF-κB- and NFAT-response assays. In contrast, 17-AAG treatment (and RANKL treatment) increased both MITF (microphthalmia-associated transcription factor) protein levels and MITF-dependent vATPase-d2 (V-type proton ATPase subunit d2) gene promoter activity. These results indicate that HSP90 inhibitors enhance osteoclast differentiation in an NFATc1-independent manner that involves elevated MITF levels and activity.

Temporal patterns of osteoclast formation and activity following withdrawal of RANKL inhibition.

Rebound bone loss following denosumab discontinuation is an important clinical challenge. Current treatment strategies to prevent this fail to suppress the rise and overshoot in osteoclast-mediated bone resorption. In this study, we use a murine model of denosumab treatment and discontinuation to show the temporal changes in osteoclast formation and activity during RANKL inhibition and withdrawal. We show that the cellular processes that drive the formation of osteoclasts and subsequent bone resorption following withdrawal of RANKL inhibition precede the rebound bone loss. Furthermore, a rise in serum TRAP and RANKL levels is detected before markers of bone turnover used in current clinical practice. These mechanistic advances may provide insight into a more defined window of opportunity to intervene with sequential therapy following denosumab discontinuation.

Stopping denosumab, a medication commonly used to improve bone mass by blocking formation of bone resorbing osteoclasts, leads to a rebound loss in the bone which was gained during treatment. Current strategies to prevent this bone loss fail in most cases as they are unable to prevent the rise and overshoot in bone resorption by osteoclasts. Thie stems from an incomplete understanding of how osteoclasts behave during denosumab treatment and after treatment is discontinued. We use a mouse model of this phenomenon to show how osteoclast formation and activity changes throughout this process. We show that increases in the processes that drive the formation of osteoclasts can be detected in the circulation before bone loss occurs. These findings could therefore provide insight into a targeted 'window of opportunity' to intervene and prevent the rebound bone loss following stopping denosumab in patients.

Minimally invasive longitudinal intravital imaging of cellular dynamics in intact long bone.

Intravital two-photon microscopy enables deep-tissue imaging at high temporospatial resolution in live animals. However, the endosteal bone compartment and underlying bone marrow pose unique challenges to optical imaging as light is absorbed, scattered and dispersed by thick mineralized bone matrix and the adipose-rich bone marrow. Early bone intravital imaging methods exploited gaps in the cranial sutures to bypass the need to penetrate through cortical bone. More recently, investigators have developed invasive methods to thin the cortical bone or implant imaging windows to image cellular dynamics in weight-bearing long bones. Here, we provide a step-by-step procedure for the preparation of animals for minimally invasive, nondestructive, longitudinal intravital imaging of the murine tibia. This method involves the use of mixed bone marrow radiation chimeras to unambiguously double-label osteoclasts and osteomorphs. The tibia is exposed by a simple skin incision and an imaging chamber constructed using thermoconductive T-putty. Imaging sessions up to 12 h long can be repeated over multiple timepoints to provide a longitudinal time window into the endosteal and marrow niches. The approach can be used to investigate cellular dynamics in bone remodeling, cancer cell life cycle and hematopoiesis, as well as long-lived humoral and cellular immunity. The procedure requires an hour to complete and is suitable for users with minimal prior expertise in small animal surgery.

Visualisation of tumour cells in bone in vivo at single-cell resolution.

The skeleton is a common site for the establishment of distant metastases. Once cancers occupy bone, the prognosis is poor as disease recurrence and visceral spread is imminent. Understanding the pathways and cellular interactions, which regulate tumour cell seeding, dormancy and growth in bone, is pertinent to improving outcomes for patients with advanced cancers. Advances in imaging techniques have facilitated the development of the concept that the behavior of bone marrow resident cells dictates the fate of tumour cells upon arrival in bone. This review summarises recent findings achieved through intravital imaging. It highlights the importance of developing both longitudinal static and acute dynamic data to develop our understanding of tumour cell engraftment, dormancy, activation and the subsequent establishment of metastases. We also describe how imaging techniques have developed our knowledge of the elements that make up the complex bone microenvironment which tumour cells interact with to survive and grow. We also discuss how through combining these imaging insights with single cell RNA sequencing data, we are entering a new era of research which has the power to define the cell-cell interactions which control tumour cell growth in bone.

Osteocyte transcriptome mapping identifies a molecular landscape controlling skeletal homeostasis and susceptibility to skeletal disease.

Osteocytes are master regulators of the skeleton. We mapped the transcriptome of osteocytes from different skeletal sites, across age and sexes in mice to reveal genes and molecular programs that control this complex cellular-network. We define an osteocyte transcriptome signature of 1239 genes that distinguishes osteocytes from other cells. 77% have no previously known role in the skeleton and are enriched for genes regulating neuronal network formation, suggesting this programme is important in osteocyte communication. We evaluated 19 skeletal parameters in 733 knockout mouse lines and reveal 26 osteocyte transcriptome signature genes that control bone structure and function. We showed osteocyte transcriptome signature genes are enriched for human orthologs that cause monogenic skeletal disorders (P = 2.4 × 10-22) and are associated with the polygenic diseases osteoporosis (P = 1.8 × 10-13) and osteoarthritis (P = 1.6 × 10-7). Thus, we reveal the molecular landscape that regulates osteocyte network formation and function and establish the importance of osteocytes in human skeletal disease.

An atlas of genetic influences on osteoporosis in humans and mice.

Osteoporosis is a common aging-related disease diagnosed primarily using bone mineral density (BMD). We assessed genetic determinants of BMD as estimated by heel quantitative ultrasound in 426,824 individuals, identifying 518 genome-wide significant loci (301 novel), explaining 20% of its variance. We identified 13 bone fracture loci, all associated with estimated BMD (eBMD), in ~1.2 million individuals. We then identified target genes enriched for genes known to influence bone density and strength (maximum odds ratio (OR) = 58, P = 1 × 10-75) from cell-specific features, including chromatin conformation and accessible chromatin sites. We next performed rapid-throughput skeletal phenotyping of 126 knockout mice with disruptions in predicted target genes and found an increased abnormal skeletal phenotype frequency compared to 526 unselected lines (P < 0.0001). In-depth analysis of one gene, DAAM2, showed a disproportionate decrease in bone strength relative to mineralization. This genetic atlas provides evidence linking associated SNPs to causal genes, offers new insight into osteoporosis pathophysiology, and highlights opportunities for drug development.

Author Correction: An atlas of genetic influences on osteoporosis in humans and mice.

In the version of this article initially published, in Fig. 5a, the data in the right column of 'DAAM2 gRNA1' were incorrectly plotted as circles indicating 'untreated' rather than as squares indicating 'treated'. The error has been corrected in the HTML and PDF versions of the article.

Presence of residual beta cells and co-existing islet autoimmunity in the NOD mouse during longstanding diabetes: a combined histochemical and immunohistochemical study.

During type 1 diabetes, most beta cells die by immune processes. However, the precise fate and characteristics of beta cells and islet autoimmunity after onset are unclear. Here, the extent of beta cell survival was determined in the non-obese diabetic (NOD) mouse during increasing duration of disease and correlated with insulitis. Pancreata from female NOD mice at diagnosis and at 1, 2, 3 and 4 weeks thereafter were analysed immunohistochemically for insulin, glucagon and somatostatin cells and glucose transporter-2 (glut2) and correlated with the degree of insulitis and islet immune cell phenotypes. Insulitis, although variable, persisted after diabetes and declined with increasing duration of disease. During this period, beta cells also declined sharply whereas glucagon and somatostatin cells increased, with occasional islet cells co-expressing insulin and glucagon. Glut2 was absent in insulin-containing cells from 1 week onwards. CD4 and CD8 T cells and macrophages persisted until 4 weeks, in islets with residual beta cells or extensive insulitis. We conclude that after diabetes onset, some beta cells survive for extended periods, with continuing autoimmunity and expansion of glucagon and somatostatin cells. The absence of glut2 in several insulin-positive cells suggests that some beta cells may be unresponsive to glucose.

Histone deacetylase activity mediates acquired resistance towards structurally diverse HSP90 inhibitors.

Heat shock protein 90 (HSP90) regulates multiple signalling pathways critical for tumour growth. As such, HSP90 inhibitors have been shown to act as effective anticancer agents in preclinical studies but, for a number of reasons, the same effect has not been observed in the clinical trials to date. One potential reason for this may be the presence of de novo or acquired resistance within the tumours. To investigate mechanisms of resistance, we generated resistant cell lines through gradual dose escalation of the HSP90 inhibitor 17-allylamino-17-demethoxygeldanamycin (17-AAG). The resultant resistant cell lines maintained their respective levels of resistance (7-240×) in the absence of 17-AAG and were also cross-resistant with other benzoquinone ansamycin HSP90 inhibitors. Expression of members of the histone deacetylase family (HDAC 1, 5, 6) was altered in the resistant cells. To determine whether HDAC activity contributed to resistance, pan-HDAC inhibitors (TSA and LBH589) and the class II HDAC-specific inhibitor SNDX275 were found to resensitize resistant cells towards 17-AAG and 17-dimethylaminoethylamino-17-demethoxygeldanamycin. Most significantly, resistant cells were also identified as cross-resistant towards structurally distinct HSP90 inhibitors such as radicicol and the second-generation HSP90 inhibitors CCT018159, VER50589 and AUY922. HDAC inhibition also resensitized resistant cells towards these classes of HSP90 inhibitors. In conclusion, we report that prolonged 17-AAG treatment results in acquired resistance of cancer cells towards not just 17-AAG but also to a spectrum of structurally distinct HSP90 inhibitors. This acquired resistance can be inhibited using clinically relevant HDAC inhibitors. This work supports the potential benefit of using HSP90 and HDAC inhibitors in combination within the clinical setting.