Biphasic regulation of osteoblast development via the ERK MAPK–mTOR pathway
Abstract
Emerging evidence supports that osteogenic differentiation of skeletal progenitors is a key determinant of overall bone formation and bone mass. Despite extensive studies showing the function of mitogen-activated protein kinases (MAPKs) in osteoblast differentiation, none of these studies show in vivo evidence of a role for MAPKs in osteoblast maturation subsequent to lineage commitment. Here, we describe how the extracellular signal-regulated kinase (ERK) pathway in osteoblasts controls bone formation by suppressing the mechanistic target of rapamycin (mTOR) pathway. We also show that, while ERK inhibition blocks the differentiation of osteogenic precursors when initiated at an early stage, ERK inhibition surprisingly promotes the later stages of osteoblast differentiation. Accordingly, inhibition of the ERK pathway using a small compound inhibitor or conditional deletion of the MAP2Ks Map2k1 (MEK1) and Map2k2 (MEK2), in mature osteoblasts and osteocytes, markedly increased bone formation due to augmented osteoblast differentiation. Mice with inducible deletion of the ERK pathway in mature osteoblasts also displayed similar phenotypes, demonstrating that this phenotype reflects continuous postnatal inhibition of late-stage osteoblast maturation. Mechanistically, ERK inhibition increases mitochondrial function and SGK1 phosphorylation via mTOR2 activation, which leads to osteoblast differentiation and production of angiogenic and osteogenic factors to promote bone formation. This phenotype was partially reversed by inhibiting mTOR. Our study uncovers a surprising dichotomy of ERK pathway functions in osteoblasts, whereby ERK activation promotes the early differentiation of osteoblast precursors, but inhibits the subsequent differentiation of committed osteoblasts via mTOR-mediated regulation of mitochondrial function and SGK1.
Editor's evaluation
This work provides a novel insight into regulation of osteogenesis by ERK-mTOR pathway. The authors proposed that the effect of Erk pathway would be mediated mTOR2-SGK1. The mitochondrial metabolisms appears to be involved in this regulation. This study is well performed, and the manuscript is clearly written.
https://doi.org/10.7554/eLife.78069.sa0Introduction
Signal transduction and transcription programs are coordinated to determine the commitment, proliferation, and differentiation of skeletal progenitors to osteoblasts. Once skeletal progenitors are committed to osteoprogenitors, they become preosteoblasts, proliferate, and differentiate into mature osteoblasts producing mineral and extracellular matrix proteins. Mature osteoblasts lastly terminally differentiate into osteocytes and embed themselves within the bone matrix (Salhotra et al., 2020). Proliferation is dominant at an early stage of osteogenesis, while proliferation rates decrease and extracellular mineralization increases during osteoblast differentiation (Infante and Rodríguez, 2018; Long, 2011). The balance between proliferation and differentiation during osteogenesis is tightly regulated by the sequential activation of signaling cascades, such as the mitogen-activated protein kinase (MAPK) and the mammalian/mechanistic target of rapamycin (mTOR) pathways.
The extracellular signal-regulated kinase (ERK) MAPK pathway transduces activation signals from extracellular cues, such as bone morphogenic proteins (BMPs) and wingless-related integration site (WNT) ligands, determining the specification of skeletal progenitors to osteoprogenitors (Ge et al., 2007; Greenblatt et al., 2013; Kim et al., 2019; Matsushita et al., 2009). Extensive studies have demonstrated the importance of the ERK MAPK pathway in promoting early commitment and differentiation of skeletal progenitors to the osteoblast lineage and skeletal mineralization (Greenblatt et al., 2022). For example, conditional deletion of Mapk1 (ERK2) and Mapk2 (ERK1) in the mesenchyme (Prrx1Cre) causes severe limb deformity and bone defects due to impaired osteogenesis (Matsushita et al., 2009). Similarly, mice lacking the MAP2Ks in the ERK pathway, Map2k1 (MAP2K1, MEK1) and Map2k2 (MAP2K2, MEK2), in osteoprogenitors (Sp7Cre) showed low bone mass and a cleidocranial dysplasia-like phenotype (CCD), similar to that seen in mice and humans with runt-related transcription factor 2 (RUNX2) haploinsufficiency (Ge et al., 2007; Ge et al., 2009; Kim et al., 2019). Thus, ERK-mediated activation of RUNX2 is required for the commitment of skeletal progenitors to osteoprogenitors and the subsequent proliferation of osteoprogenitors at early differentiation stages. Additionally, the ERK MAPK pathway mediates phosphorylation of ribosomal s6 kinase 2 (RSK2; Yang et al., 2004). Rsk2-deficient mice and humans are osteopenic due to impaired production of type I collagen, the major organic component of bone. RSK2 mutations have been reported to be associated with Coffin–Lowry syndrome showing prominent skeletal dysplasia features (David et al., 2005; Hanauer and Young, 2002). As neurofibromin 1 (NF1) functions as a negative regulator of the ERK MAPK pathway, loss of function mutations in NF1 cause neurofibromatosis type 1 (NF1) syndrome, characterized by a complex set of bone phenotypes that can include osteopenia and impaired fracture healing as well as cutaneous neurofibromas (Bok et al., 2020; Crawford and Schorry, 1999; de la Croix Ndong et al., 2014). Similar to NF1, gain of function mutations in protein-tyrosine phosphatase nonreceptor-type 11 (PTPN11, SHP2) enhance ERK signaling in osteoblasts, causing Noonan syndrome associated with impaired skeletal mineralization and short stature (Binder, 2009; Choudhry et al., 2012; Roberts et al., 2007). It remains unclear how modulation of ERK signaling can produce these diverse and clinically important phenotypes. These considerations motivate consideration of how the role of ERK may vary based on context. In particular, the role of ERK in regulating the stages of osteoblast differentiation and function after lineage commitment remains unresolved.
In addition to RUNX2, the ERK MAPK pathway controls the mTOR pathway, a signaling modulator of mitochondrial biogenesis that regulates cellular energy metabolism (Morita et al., 2013). Cross-regulation between these pathways in various cancer cells has been also reported (Mendoza et al., 2011). ERK inhibition by an MEK inhibitor activates the mTORC2–AKT signaling axis downstream of epidermal growth factor (EGF; Yu et al., 2002), and a similar increase in AKT activation by ERK inhibition was reported in anaplastic lymphoma kinase-addicted neuroblastoma (Umapathy et al., 2017). Mechanistically, impaired ERK activation by the MEK inhibitor trametinib treatment enhanced phosphorylation of the mTORC2 complex subunits, Rictor or stress-activated protein kinase (SAPK)-interacting protein 1 (SIN1), via activation of induced integrin-linked kinase, resulting in AKT activation (Pairo-Castineira et al., 2021; Umapathy et al., 2017). Thus, determining the function and the regulatory mechanism of ERK in osteoblast differentiation and bone formation is clinically significant as the ERK pathway is the most frequently mutated pathway in cancer and, accordingly, ERK pathway components, including rat sarcoma (RAS), rapidly accelerated fibrosarcoma, and MEK isoforms, are all well-established targets to treat a wide variety of ERK pathway-activated cancers (Hong et al., 2020; Hyman et al., 2015; Robert et al., 2015; Tiacci et al., 2021). Given that many of these patients are at an increased risk for fracture, it is important to weigh the risks and benefits when determining the role of the ERK pathway in the regulation of bone mass and the potential skeletal impact of ERK pathway-directed therapies.
Here, we establish dynamic roles for the ERK MAPK pathway during osteogenesis, with the ERK–mTOR signaling axis acting in osteoblasts to regulate mitochondria-dependent energy metabolism and SGK1-mediated production of pro-angiogenic and osteogenic factors. ERK inhibition in mature osteoblasts promotes bone formation via mTOR-mediated activation of mitochondria and SGK1. Accordingly, inhibition of the mTOR pathway using rapamycin attenuated bone accrual in MEK1/2-deficient mice. Our study uncovers a surprising dichotomy of ERK pathway functions in osteoblasts, with the ERK pathway promoting bone formation in osteoblast precursors but inhibiting bone formation in committed osteoblasts via mTOR pathway.
Results
A biphasic role for the ERK MAPK pathway in osteogenesis
To investigate the role of the ERK MPAK pathway in osteoporosis, the MEK inhibitor trametinib was treated to mice with estrogen deficiency-induced bone loss. Twelve-week-old mice were subjected to the ovariectomy (OVX) model of postmenopausal osteoporosis or sham surgery and were subsequently treated with trametinib via oral gavage for 8 weeks. While OVX surgery induced a significant bone loss in the femur of vehicle-treated mice, femoral bone loss was prevented by treatment with trametinib, as shown by a greater trabecular bone volume (Tb. BV) and trabecular thickness (Tb. Th; Figure 1A). Remarkably, dynamic histomorphometry analysis demonstrated an increase in bone formation rate per bone surface (BFR/BS) and osteoblast surface on the trabecular bone area (Ob.S/BS) of trametinib-treated OVX mice (Figure 1B, C). Thus, ERK inhibition prevents estrogen deficiency-induced bone loss by promoting bone formation due to augmented osteoblast activity.
To gain insights into the mechanism underlying this finding, we examined the effects of ERK inhibition on the proliferation and differentiation of osteoblast precursors using human bone marrow-derived mesenchymal stromal cells (BMSCs). Expression of key osteogenic commitment factors, RUNX2 and SP7 (Osterix), was markedly upregulated 8 days after osteogenic culture as expected, establishing that the BMSCs at this time point represent committed osteoblast-lineage cells (Figure 1—figure supplement 1). ERK inhibition starting at day 8 of culture promoted osteoblast differentiation, whereas ERK inhibition from day 0 of culture conversely inhibited proliferation and osteogenic potential of BMSCs (Figure 1D, E). Similar findings were observed with committed osteoblasts isolated from mouse calvarium, which also showed increased osteogenic differentiation in the presence of trametinib. This was evident in ERK inhibitor induced increases in alkaline phosphatase (ALP), mineralization, and expression of osteogenic genes, including integrin-binding sialoprotein (Ibsp) and osteocalcin (Bglap; Figure 1F–H). Notably, cell proliferation rates were largely unaffected (Figure 1I). These results suggest that while the ERK MAPK pathway plays a positive role in the early stages of osteogenesis, it acts as a negative regulator of the later stages of osteoblast differentiation. Thus, ERK inhibition in late-stage osteoblasts by trametinib treatment is likely to be responsible for the increased osteoblast activity observed in ERK inhibitor-treated OVX mice.
Impaired ERK signaling in mature osteoblasts promotes bone formation
To confirm these findings using genetic approaches in vivo, Map2k1 (MEK1) and Map2k2 (MEK2), which are upstream of ERK1/2, were conditionally deleted in mature osteoblasts and osteocytes by crossing with Dmp1Cre (Dentin matrix protein1) mice, which targets conditional gene deletion to mature osteoblasts and osteocytes (Lu et al., 2007). Deletion of MEK1/2 was confirmed in RNA from the tibiae of WT and Map2k1fl/fl;Dmp1Cre;Map2k2−/− (double knockout, dKODmp1) mice (Figure 2A), and resulted in a significant increase in femoral bone mass, as evidenced by greater trabecular bone volume (Tb. BV), thickness (Tb. Th), number (Tb. N), and cortical thickness (Ct. Th; Figure 2B, C). Likewise, skull thickness and L4 vertebrae trabecular bone volume/tissue volume (Tb. BV/TV) were markedly increased (Figure 2D, E). These skeletal phenotypes require the deletion of both MEK1 and MEK2 alleles as deletion of MEK1 (Map2k1 fl/fl;Dmp1Cre) or MEK2 (Map2k2−/−) alone did not show an increase in bone mass compared to WT femurs (Figure 2C), suggesting a redundant role for MEK1 and MEK2 in osteoblasts. Consistent with the increased bone mass observed, a three-point bending test of dKODmp1 mice revealed an increase in maximum bending moment and load at failure, demonstrating enhanced bone mechanical strength (Figure 2F). Thus, the bone formed by ERK pathway ablation in late-stage osteoblasts has the biomechanical properties of mature, physiologic bone and not that of fracture-prone pathologic bone. Finally, dynamic histomorphometry analysis demonstrated a significant increase in BFR/BS and mineral apposition rate (MAR) of dKODmp1 femurs, indicating enhanced in vivo osteoblast activity in the absence of MEK1/2 (Figure 2G). Therefore, deletion of MEK1/2 in late-stage osteoblast-lineage cells augmented their bone-forming activity and bone accrual and strength. Importantly, in vivo osteoclast activity and numbers were markedly elevated in dKODmp1 mice, as shown by an increase in tartrate-resistant acid phosphatase (TRAP)-positive osteoclasts, bone erosion surface, and serum levels of the bone resorption marker C-terminal telopeptide type I collagen (CTx-I; Figure 2—figure supplement 1A, B). This was accompanied with reduced expression of osteoprotegerin (OPG, Tnfrsf11b) without any alteration in receptor activator of nuclear factor kappa-Β ligand (RANKL, Tnfsf11) expression (Figure 2—figure supplement 1C), suggesting that the upregulation of RANKL/OPG ratio may be responsible for the enhanced osteoclast activity observed in these mice.
To confirm that deletion of MEK1/2 in mature osteoblasts is responsible for the development of these phenotypes, we generated inducible, osteoblast-specific MEK1/2-knockout mice by crossing Map2k1fl/fl;Map2k2−/− mice with BglapCreERT (osteocalcin) mice expressing a tamoxifen-induced Cre recombinase in mature osteoblasts (Maes et al., 2010; Map2k1fl/fl; BglapCreErt;Map2k2−/−, dKOBglap). Treatment of mice with tamoxifen (dKOBglap) significantly increased trabecular bone mass of the femur (Figure 2H,I), demonstrating that the inhibitory role of the ERK MAPK pathway is specific to mature osteoblasts and reflects the continuous function of ERK in these cells and is not secondary to a developmental abnormality.
Given previous studies showing the importance of the ERK MAPK pathway in osteocyte survival (Plotkin et al., 2005; Ru and Wang, 2020), terminal deoxynucleotidyl transferase dUTP nick-end labeling (TUNEL) staining was performed in dKODmp1 femurs to assess osteocyte cell death (Figure 2J, top). Osteocyte apoptosis was substantially increased in the cortical bone of dKODmp1 mice relative to control mice, which corresponded to decreased osteocyte numbers and increased empty lacunae (a cavity within the bone; Figure 2J, bottom, K). Likewise, transmission electron microscopy (TEM) analysis of dKODmp1 femurs demonstrated that terminally differentiated osteocytes embedded in the bone matrix undergo apoptosis (Figure 2 L3) while osteocytes at early (Figure 2 L1) and intermediate (Figure 2 L2) differentiation stages were grossly normal in the absence of MEK1/2 (Figure 2L). This is consistent with RT-PCR analysis showing that mRNA levels of early osteocyte markers, including dentin matrix protein1 (Dmp1) and phosphate regulating endopeptidase homolog X-linked (Phex) were comparable between control and dKODmp1 tibias. However, expression of the late osteocyte marker, sclerostin (Sost), was significantly decreased (Figure 2—figure supplement 1C). These results suggest that ERK activation is required for the survival of late/fully mature osteocytes, but not the transition of mature osteoblasts to early osteocytes. Taken together, these results demonstrate that ERK inhibition in mature osteoblasts promotes bone formation, while impaired ERK signaling in late osteocytes results in cell death.
ERK inhibition increases osteogenesis and production of angiogenic and osteogenic factors
As seen in dKODmp1 mice, Cre-induced deletion of MEK1 and MEK2 MAP2Ks in committed osteoblasts increased ALP activity, mineralization, and osteogenic gene expression. Deletion of MEK1 or MEK2 alone did not impact osteogenic differentiation (Figure 3A–C, Figure 3—figure supplement 1A). While proliferation rates of control osteoblasts were markedly increased under osteogenic conditions, little to no increase was observed in the absence of MEK1/2 (Figure 3D), suggesting that ERK inhibition may control a shift toward differentiation at the expense of proliferation. To gain insights into the mechanisms by which ERK inhibition promotes osteoblast differentiation, we performed whole transcriptome analysis of dKO osteoblasts 6 days after cell growth (GM) or osteogenic (OIM) culture conditions. Gene ontology analysis demonstrated that, under cell growth conditions, upregulated genes were highly enriched in the pathways involved in ossification, skeletal system morphogenesis, and connective tissue development (Figure 3E), consistent with an increase in ALP activity and expression of osteogenic genes in undifferentiated dKO osteoblasts (Figure 3—figure supplement 1B–F). These results suggest that ERK inhibition induces spontaneous osteogenic differentiation by upregulating transcription of genes favorable to osteogenic differentiation. Under osteogenic conditions, gene ontology analysis of differentiated dKO osteoblasts revealed upregulation of genes associated with extracellular matrix organization, while genes related to cell proliferation, including nuclear division and cell cycle, were downregulated (Figure 3F, Figure 3—figure supplement 2A–D). These results suggest that impaired ERK signaling in later stage committed osteoblasts upregulates transcriptional programs associated with osteogenic differentiation while suppressing the induction of genes associated with cell proliferation. Intriguingly, the genes associated with the pro-osteogenic WNT/β-catenin and TGF-β signaling pathways were highly upregulated in dKO osteoblasts under osteogenic, but not undifferentiated, conditions (Figure 3—figure supplement 3A–C). This is consistent with luciferase assays showing a significant increase in activities of key transcription factors of WNT and TGF-β signaling, β-catenin and Smads, in the absence of MEK1/2 under osteogenic, but not undifferentiated, conditions (Figure 3—figure supplement 4A). Likewise, protein levels of β-catenin and Smad3 in dKO cells cultured under osteogenic conditions were markedly increased (Figure 3—figure supplement 4B), suggesting that ERK inhibition at later stages of osteogenic differentiation enhances the pro-osteogenic WNT and TGF-β signaling pathways. Accordingly, pharmacological inhibition of these signaling pathways reversed enhanced osteogenesis of dKO cells, as shown by reduced ALP activity and osteogenic gene expression in a dose-dependent manner (Figure 3—figure supplement 4C–F). These results suggest that enhanced WNT and TGF-β signaling at later stages of osteogenic differentiation of dKO osteoblasts might be major contributors to drive osteogenic differentiation program.
While the expression and transcriptional activity of RUNX2 and SP7 were not strongly impacted in mature dKO osteoblasts (Figure 3—figure supplement 5A–C), ERK inhibition in mature osteoblasts markedly increased the expression and transcription activity of ATF4 (OSE1-luc; Figure 3—figure supplement 5D, E). These results suggest that during the later stages of osteoblast differentiation, the ERK MAPK pathway is required for post-translational modifications of ATF4, but not RUNX2 and SP7. This finding provides a clear mechanistic distinction between the role of the ERK pathway in early osteogenic commitment, which relates to regulation of RUNX2 and SP7, from the role of the ERK pathway in later stage osteoblasts which, by contrast, includes regulation of signaling pathways related to ATF4, WNT/β-catenin, and TGF-β.
Whole transcriptome analysis of dKO osteoblasts also revealed significant upregulation of secreted cytokines and/or growth factors involved in ossification and angiogenesis/vasculature development (Figure 3F), such as fibroblast growth factor 9 (Fgf9; Behr et al., 2010), ectonucleotide pyrophosphatase/phosphodiesterase 2 (Enpp2; Cholia et al., 2015), transforming growth factor-beta 2 (Tgfb2; Wu et al., 2016), insulin-like growth factors 1, 2 (Igf1, 2; Yakar et al., 2018), bone morphogenic proteins 1, 3 (Bmp1, 3; Wu et al., 2016), C-C Motif Chemokine Ligand 5 (Ccl5; Suffee et al., 2012), and Ccl11 (Salcedo et al., 2001; Figure 3G). Accordingly, conditioned medium (CM), collected from dKO osteoblasts, markedly increased osteogenic potentials of wildtype BMSCs (Figure 3H) and capillary tube formation of bone marrow-derived endothelial progenitor outgrowth cells (EPOCs, Figure 3I). Likewise, an immunofluorescence analysis of dKODmp1 femurs showed elevated levels of CD31- and endomucin-positive skeletal vasculature supporting osteoblast development (Xu et al., 2018; Figure 3J). These results suggest that pro-angiogenic and/or -osteogenic factors secreted from dKO osteoblasts enhance angiogenesis and osteogenesis, promoting bone formation. Taken together, MEK1/2 deficiency does not only increase osteoblast differentiation by upregulating osteogenic transcriptional programs, but also improves bone-forming environment by producing pro-angiogenic and/or -osteogenic factors.
ERK inhibition enhances mitochondrial energy metabolism in osteoblasts
It has been well established that the ERK MAPK pathway controls cellular metabolism in cancer cells and highly proliferating stem cells (Lee et al., 2017; Papa et al., 2019). Mature osteoblasts produce collagen for extracellular matrix formation using adenosine triphosphate (ATP) as a major energy resource (Lee et al., 2017). Moreover, glucose and glutamine, used for mitochondrial ATP production, are important for cell proliferation at early stages of osteoblast differentiation and matrix mineralization at later differentiation stages (Karner et al., 2015; Wei et al., 2015). Thus, we hypothesized that ERK-mediated regulation of glucose/glutamine metabolism may be involved in the determination of osteogenic potential and differentiation. Remarkably, the differentiation of dKO osteoblasts was substantially reduced by glutaminase inhibition but only modestly impacted by glucose deprivation (Figure 4A, B), demonstrating the importance of glutamine metabolism in ERK-mediated regulation of osteogenesis. Specifically, the responsiveness of dKO osteoblasts to the glutaminase inhibitor (bis-2-[5-phenylacetamido-1,3,4-thiadiazol-2-yl]ethyl sulfide or BPTES) was greater than that of control osteoblasts, as shown by a significant reduction in ALP activity, cell proliferation, and osteogenic gene expression in the presence of BPTES (Figure 4B, C). Likewise, these cells showed elevated levels of glutaminase (Gls), a key enzyme that produces glutamate from glutamine (Figure 4D), and increased glutamate production and release (Figure 4E), suggesting enhanced glutamine metabolism in dKO osteoblasts. Of note, in contrast to previous studies showing critical roles of glucose metabolism (Wei et al., 2015), effect of glucose depravation on osteoblast maturation was intact in the absence of MEK1/2 (Figure 4A), suggesting that the ERK MAPK pathway is not required for glucose metabolism during the later stages of osteoblast differentiation. This finding provides a clear mechanistic distinction between the role of the ERK pathway in early osteogenic commitment which relates to regulation of glucose metabolism, from the role of the ERK pathway in later stage osteoblasts which, by contrast, includes regulation of glutamine metabolism.
Given the importance of glutamine metabolism in mitochondria-mediated energy production, we next examined the intracellular ATP content in dKO osteoblasts and found increased ATP levels (Figure 4F). Additionally, mitochondrial numbers and membrane potentials and a ratio of mitochondrial DNA (mtDNA) levels to nuclear DNA (nDNA) levels in dKO or trametinib-treated osteoblasts were all markedly increased compared to control or vehicle-treated osteoblasts (Figure 4G–I, Figure 4—figure supplement 1). This corresponds to enhanced oxygen consumption rate (OCR), an indicator of cellular respiration, in dKO osteoblasts (Figure 4J), demonstrating that ERK inhibition promotes mitochondrial function and metabolic activity in committed osteoblasts. Finally, levels of reactive oxidative species (ROS; Figure 4K) and mRNA levels of peroxisome proliferator-activated receptor gamma coactivator 1-alpha (Pgc1a), dynamin-related protein 1 (Drp1), and mitofusin 1 (Mfn1; Figure 4L), genes that are important for mitochondrial biogenesis, fission, and fusion, were also markedly elevated in dKO osteoblasts. Taken together, these results suggest that ERK inhibition increases osteoblast differentiation via augmented glutamine metabolism and mitochondrial function.
ERK inhibition promotes bone formation via the mTOR pathway in osteoblasts
As a metabolic hub in response to nutrients and growth factors, the mTOR pathway plays a crucial role in mitochondria-mediated cell metabolism during skeletal development and homeostasis (Chen and Long, 2018). Two main mTOR signaling complexes, mTORC1 and mTORC2, are composed of discrete protein-binding partners to regulate growth, motility, and metabolism in osteoblasts. mTORC1 activation mainly promotes osteoblast proliferation and differentiation in response to bone anabolic signals and their regulatory processes (Chen et al., 2014; Lim et al., 2016). On the other hand, mTORC2 activation not only mediates osteogenesis and skeletal mineralization in response to growth factors or mechanical loading, but it also plays a role in the restoration of metabolic homeostasis in response to nutrient fluctuations (Szwed et al., 2021). Additionally, mTORC2 has been reported to be associated with age-related osteoporosis since impaired mTORC2 activation in aged mice switches the fate of skeletal progenitors to adipocytes from osteoblasts (Lai et al., 2016). Given that ERK functions upstream of the mTOR pathway in regulating mitochondria-mediated energy metabolism (Morita et al., 2013), we tested the ability of mTOR inhibition to reverse the enhanced osteogenic differentiation seen after MEK1/2 deletion (Figure 5A). Rapamycin was used to inhibit both mTORC1 and mTORC2 activation (Lamming et al., 2012) in dKO osteoblasts, demonstrating that rapamycin treatment attenuated ALP activity (Figure 5B) and expression of osteogenic genes (Figure 5C) in the absence of MEK1/2. Of note, in contrast to many contexts where rapamycin solely inhibits the mTORC1 pathway, including p70S6 kinase (p70S6K), rapamycin treatment was effective at inhibiting the phosphorylation of the mTORC2 downstream molecules, protein kinase B (AKT) and serum/glucocorticoid regulated kinase 1 (SGK1), in dKO osteoblasts (Figure 5—figure supplement 1). These results suggest that rapamycin treatment decreases the differentiation of dKO osteoblasts by inhibiting both mTORC1 and mTOCR2 pathways. To test this hypothesis in vivo, 2-week-old dKODmp1 mice were treated with rapamycin for 6 weeks and femoral bone mass was assessed in 8-week-old mice using microCT. While little effect of rapamycin on femoral bone mass was seen in control mice, trabecular bone mass and cortical bone thickness were both decreased in dKODmp1 mice when treated with rapamycin relative to vehicle, as shown by a significant reduction in trabecular bone volume (Tb. BV), number (Tb. N), thickness (Tb. Th), and cortical thickness (Ct.Th) (Figure 5D, E). These results suggest, as observed in dKO osteoblasts, mTOR inhibition can reverse the enhanced bone-forming activity by MEK1/2-deletion in mice. A histologic analysis of rapamycin-treated dKODmp1 femurs revealed that mTOR inhibition was also effective in preventing the osteocyte apoptosis occurring in the absence of MEK1/2 (Figure 5F). Additionally, osteoclast surface (Oc.S) in the trabecular bone was slightly decreased in these femurs, suggesting that rapamycin treatment may directly suppress osteoclast differentiation and activity (Glantschnig et al., 2003; Kneissel et al., 2004; Figure 5G). Taken together, impaired ERK signaling is likely to promote osteogenic differentiation and bone formation via activation of the mTOR pathway, whose inhibition can reverse dKODmp1 phenotypes in vitro and in vivo.
ERK inhibition promotes osteoblast differentiation via mTORC2/SGK1 activation
To understand how ERK controls mTOR activation in osteoblasts, we examined the effects of MEK1/2 deficiency on mTORC1 and mTORC2 signaling. While mTORC1 primarily promotes cell growth via activation of p70S6K and inactivation of eIF4E-binding protein 1 (4EBP1), mTORC2 regulates cell survival and proliferation via AKT phosphorylation as well as cellular metabolism via activation of SGK1 and N-Myc downstream regulated 1 (NDRG1) (Saxton and Sabatini, 2017). In comparison to control osteoblasts, dKO osteoblasts showed a significant increase in phosphorylation levels of mTORC2 downstream molecules, including AKT, SGK1, and NDRG1, whereas little to no alteration in the phosphorylation levels of the mTORC1 downstream molecules p70S6K and 4EBP1 were detected in the absence of MEK1/2 (Figure 6A). Together with previous studies (Brown et al., 2020; Umapathy et al., 2017), these results suggest that ERK inhibition in osteoblasts activates mTORC2- but not mTORC1-mediated signal transduction via phosphorylation of mTORC2 subunits, Rictor or SIN1.
In a previous study, SGK1, that functions downstream of mTORC2, played an important role in osteogenic trans-differentiation and calcification of vascular muscle cells (Voelkl et al., 2018). Since the genes associated with the ossification pathway in dKO osteoblasts were highly upregulated (Figure 3E, F), we hypothesized that ERK inhibition may upregulate osteoblast differentiation by activating the mTORC2–SGK1 signaling axis. Similar to rapamycin-mediated mTOR inhibition (Figure 5B, C), Sgk1 knockdown markedly attenuated ALP activity (Figure 6B), expression of osteogenic genes (Figure 6C), and mineralization (Figure 6D) in dKO osteoblasts, not in control osteoblasts, suggesting that SGK1 functions downstream of ERK during osteoblast differentiation. Notably, the enhanced expression and transcription activity of ATF4 in mature dKO osteoblasts were reversed by Sgk1 knockdown or mitochondrial inhibitor treatment, suggesting that mTOR–SGK1 and mitochondrial function control ERK-mediated regulation of ATF4 in mature osteoblasts (Figure 3—figure supplement 4D, E). Additionally, the ability of CM collected from dKO osteoblasts to promote angiogenesis of EPOCs and osteogenesis of wildtype BMSCs was ablated by Sgk1 knockdown, as shown by a significant decrease in migration of EPOCs (Figure 6E), and ALP activity and Col1a1 mRNA levels of BMSCs (Figure 6F). Notably, little to no alteration of mitochondrial numbers and membrane potentials (Figure 6G) and OCR (Figure 6H) was observed in trametinib-treated Sgk1-deficient osteoblasts or Sgk1-deficient dKO osteoblasts, suggesting that SGK1 is dispensable for ERK-mediated mitochondrial function. Taken together, the ERK–mTORC2 signaling axis mediates osteogenesis and production of pro-angio/osteogenic factors via SGK1 activation, while controlling osteogenesis via mitochondria-dependent energy metabolism (Figure 6I).
Discussion
This study demonstrates essential roles for the ERK–mTOR pathway in osteoblast development and function. While ERK activation is crucial for cell proliferation and osteogenic potential at the early stages of osteogenic differentiation, ERK-mediated regulation of the mTOR pathway determines osteogenic differentiation in later stages. In addition, ERK activation dampens the secretion of pro-angiogenic and osteogenic factors required to create an optimal microenvironment for bone formation. Mechanistically, impaired ERK signaling in committed osteoblasts enhances osteoblast differentiation via an increase in mitochondria-mediated metabolic demands, activation of mTORC2–SGK1 signaling axis, thereby promoting bone formation (Figure 6I). This corresponds to previous reports showing that mTORC2 activation is required for mitochondrial function and homeostasis (Wang et al., 2021) and that both glycolysis and oxidative phosphorylation by mitochondria are required to meet the energy demands needed for bone formation by osteoblasts (Lee et al., 2017). However, despite previous studies showing the importance of SGK1 in osteogenic trans-differentiation of vascular smooth muscle cells (Voelkl et al., 2018) and energy metabolism in different cellular types and functions (Mason et al., 2021), our findings suggest that while ERK–mTORC2–SGK1 signaling axis is required for osteogenesis and production of pro-angiogenic and -osteogenic factors, it is dispensable for mitochondrial function. Thus, our study identified two distinct pathways downstream of the ERK–mTORC2 signaling axis controlling bone formation, mitochondria- and SGK1-mediated osteogenesis and production of pro-angiogenic and -osteogenic factors. Notably, while we demonstrated the mTORC2–SGK1 signaling axis as a critical downstream pathway for ERK during late-stage differentiation, the ERK MAPK pathway is likely to have additional substrates that contribute to its regulation of osteogenesis. Since WNT ligands can activate both the ERK and mTOR pathways independently to regulate skeletal development and homeostasis (Chen and Long, 2018; Karner and Long, 2017; Shim et al., 2013), the ERK and mTORC2–SGK1 signaling axis may also function in parallel at different stages of osteogenesis. Additionally, further investigation will be necessary to define the mechanism by which the ERK MAPK pathway suppresses mTORC2 activation in committed osteoblasts via phosphorylation of mTORC2 subunits, Rictor or SIN1 (Brown et al., 2020; Umapathy et al., 2017).
ERK activation is dynamically changed in response to various extracellular stimuli during tissue development, homeostasis, and regeneration (Lavoie et al., 2020; Patel and Shvartsman, 2018). The ERK MAPK pathway plays a critical role in the regulation of cell quiescence/proliferation, differentiation, and commitment of adult stem cells and progenitors, as ERK activation promotes cell proliferation while suppressing cell differentiation (Krishnan et al., 2021). However, very few in vivo studies have reported dynamics of ERK signaling in the skeletal system. This study demonstrates distinct roles of the ERK MAPK pathway at different stages of osteogenesis. In particular, during late-stage osteogenic differentiation, ERK activation is important for mitochondria-mediated glutamine metabolism rather than glucose metabolism (Figure 4). Previous studies have shown that osteogenic precursors not only consume glutamine to generate ATP but also synthesize metabolites to promote osteoblast viability and matrix mineralization (Yu et al., 2019). Our findings reveal that impaired ERK signaling enhances mitochondria-mediated glutamine metabolism, thereby promoting osteogenic differentiation. Remarkably, this can be reversed by the mTOR inhibitor rapamycin in vitro and in vivo. Since systemic administration of rapamycin increases bone anabolic activity in both normal bone mass and osteopenic conditions (Bateman et al., 2019; Chagin, 2016; Izawa et al., 2017), rapamycin-induced attenuation of osteogenic differentiation and bone anabolic activity in MEK1/2-deficient mice and osteoblasts is unlikely to result from nonspecific or general effects of rapamycin unrelated to the specific defects in MEK1/2-deficient mice. Together with previous studies showing the association of the mTOR pathway with energy metabolism during skeletal development and homeostasis (Chen and Long, 2018), our data suggest that the ERK MAPK pathway controls energy metabolism in committed osteoblasts via the mTOR pathway.
This study also justifies an examination of the clinical impact of ERK pathway inhibitors on bone mass and fracture risk, as the biphasic complex contribution of the ERK pathway to bone formation makes empiric study of this question important. This is especially true given that the ERK pathway is the most frequently mutated pathway in cancer (Guo et al., 2020), and thus ERK inhibitors are important tools for targeted therapy in a wide range of tumor types. This study also shows that dual treatment with ERK and mTOR inhibitors may be associated with a strong decrease in bone mass, as mTOR inhibition would ablate the anabolic effect of ERK inhibition in late-stage osteoblasts, leaving only the negative effects of ERK on early osteogenic differentiation and commitment. In mice, trametinib treatment increased bone mass in osteoporotic bone (Figure 1A–C), indicating that the pro-anabolic effect of ERK inhibition in mature osteoblasts is dominant over the anti-anabolic effects on early-stage osteoblasts. However, the relative dominance of these opposing early versus late-stage effects may depend highly on dose and duration of therapy, in addition to clinical context and patient demographics. Thus, empirical study of the clinical effects of trametinib and other ERK pathway inhibitors will be important, especially given that many cancer patients being treated with these agents will be at increased risk for skeletal fractures due to the effects of either cancer, chemotherapy, or patient demographics impacting bone mass. Given the emerging role of osteoblasts in promoting and maintaining skeletal metastases, our finding that the ERK pathway has complex effects on early- and late-stage osteoblast-lineage cells also suggests that further examination of the impact of ERK pathway inhibition of skeletal metastases will be important to understand the risks and benefits of ERK pathway inhibition (Bado et al., 2021; Swami et al., 2017; Zhang et al., 2021).
Materials and methods
Cell culture, antibodies, and reagents
Request a detailed protocolHuman BMSCs were purchased from Cyagen Biosciences and were maintained in GM (HUXMX-90011) and cultured in osteogenic medium (GUXMX-90021). Primary osteoblasts were isolated from the calvaria of 5-day-old mouse pups using collagenase type II (50 mg/ml, Worthington, LS004176)/dispase II (100 mg/ml, Roche, 10165859001). MEK1/2-deficient osteoblasts were obtained from calvaria of 5-day-old Map2k1fl/fl;Map2k2−/− neonates, transduced with lentiviruses expressing vector control or Cre recombinase, and selected with puromycin. Osteoblasts were maintained in α-minimal essential medium (Gibco) containing 10% fetal bovine serum (Corning), 2 mM L-glutamine (Corning), 1% penicillin/streptomycin (Corning), and 1% nonessential amino acids (Corning) and differentiated with ascorbic acid (200 μM, Sigma, A8960) and β-glycerophosphate (10 mM, Sigma, G9422).
Antibodies specific to MEK1 (12671), P-ERK1/2 (4376), P-p70S6K (9205), P-4EBP1 (2855), P-AKT (4060), P-NDRG1 (5482), β-catenin (8480), and RUNX2 (12556) were purchased from Cell Signaling. Antibodies specific to P-SGK1 (44–1264 G) and GAPDH (CB1001) were purchased from Thermo Scientific and EMD Millipore, respectively. Antibodies specific to ERK1/2 (A0228), p70S6K (A2190), 4EBP (A1248), AKT (A11016), NDRG1 (A2142), SGK1 (A1025), and SMAD3 (A199115) were purchased from ABclonal. Trametinib (S2673) and BPTES (5301) were purchased from Selleck Chemicals and Tocris Bioscience, respectively. Rapamycin (sc-3054) for cell culture and antimycin A (sc-20246) were purchased from Santa Cruz Biotechnology. WNT inhibitor (IWR-1-endo, 72562) and TGF-β inhibitor (SB431542, 72232) were purchased from STEMCELL Technologies. All constructs encoding shRNAs were purchased from Sigma.
Mice
Map2k1fl/fl mice and Map2k2−/− mice were generated as previously reported, respectively (Bélanger et al., 2003; Bissonauth et al., 2006) and maintained on a 129/SvEv background. To generate mature osteoblast/osteocyte-specific double knockout (Map2k1fl/fl;Dmp1Cre;Map2k2−/−, dKODmp1), Map2k1fl/fl;Map2k2−/− mice were crossed with Dmp1Cre mice (Lu et al., 2007). Sex- and age-matched littermates were used as controls for all skeletal analyses. To generate mice that harbored inducible deletion of Map2k1 (MEK1) in mature osteoblasts, Map2k1fl/fl;Map2k2−/− mice were crossed with BglapCreERT (osteocalcin) mice (Park et al., 2012) (Map2k1fl/fl;BglapCreERT;Map2k2−/−, dKOBglap). For postnatal activation of CreERT, 75 mg/kg tamoxifen (Sigma, T5648) in corn oil (Sigma) was injected intraperitoneally into 8-week-old mice once a day, for 5 consecutive days. Rapamycin (R-5000, LC laboratories) was dissolved in ethanol and diluted in filter-sterilized vehicle (5% Tween-80, 5% PEG-40, 0.9% NaCl) immediately prior to intraperitoneal injection. Vehicle or 4 mg/kg rapamycin was administered every other day for 6 weeks. All animals were used under the NIH Guide for the Care and Use of Laboratory Animals and were handled according to the animal protocol approved by the University of Massachusetts Medical School Institutional Animal Care and Use Committee (IACUC).
OVX-induced bone loss
Request a detailed protocolTo induce postmenopausal osteoporosis, 12-week-old female mice were anesthetized and bilaterally ovariectomized or subjected to a sham operation. All mice were then randomly assigned to treatment with vehicle or trametinib (0.6 mg/kg; Hu-Lieskovan et al., 2015). Inhibitor or vehicle treatment was given daily by oral gavage. After 8 weeks, mice were euthanized and subjected to bone analysis.
Bone mechanical testing
Request a detailed protocolTo determine bone mechanical properties, 8-week-old male Map2k2−/− (Ctrl) and Map2k1 fl/fl;Dmp1Cre;Map2k2−/− (dKODmp1) femora were tested, as previously described (Melville et al., 2014). Right femora were loaded to failure in three‐point bending in the anterior‐posterior direction with a span length of 6 mm at a rate of 0.1 mm/s (858 Mini Bionix: MTS, Eden Prairie, MN, USA). Bending strength was calculated based on the force and displacement data from the tests, and midshaft geometry was measured with microCT.
MicroCT
Request a detailed protocolMicroCT analysis was performed as previously described (Kim et al., 2020). Sex-matched and control littermates were used, and analysis was performed by an investigator blinded to the genotypes of the animals under analysis. Femurs excised from the indicated mice were scanned using a microCT 35 (Scanco Medical) with a spatial resolution of 7 μm. For trabecular bone analysis of the distal femur, an upper 2.1 mm region beginning 280 μm proximal to the growth plate was contoured. For cortical bone analysis of femur and tibia, a midshaft region of 0.6 mm in length was used. MicroCT scans of skulls and lumbar vertebrae (L4) were performed using isotropic voxel sizes of 12 μm. 3D reconstruction images were obtained from contoured 2D images by methods based on distance transformation of the binarized images. All images presented are representative of the respective genotypes (n > 5).
Histology, dynamic histomorphometry, and immunofluorescence
Request a detailed protocolFor histological analysis, hindlimbs were dissected from the mice, fixed in 10% neutral buffered formalin for 2 days, and decalcified by daily changes of 0.5 M tetrasodium Ethylenediaminetetraacetic acid (EDTA, pH 7.4) for 3–4 weeks. Tissues were dehydrated by passage through an ethanol series, cleared twice in xylene, embedded in paraffin, and sectioned at 7 μm thickness along the coronal plate from anterior to posterior. Decalcified femoral sections were stained with hematoxylin and eosin (H&E) or TRAP. To detect apoptotic osteocytes in bone tissue, terminal deoxynucleotidyl transferase (TdT)-mediated dUTP-digoxigenin nick-end labeling (TUNEL) staining was performed using TUNEL assay kit (Abcam, ab206386) according to the manufacturer’s instructions.
For dynamic histomorphometry analysis, 25 mg/kg calcein (Sigma, C0875) and 50 mg/kg alizarin-3-methyliminodiacetic acid (Sigma, A3882) dissolved in 2% sodium bicarbonate solution were subcutaneously injected into mice at 5-day intervals. After 2-day fixation in 10% neutral buffered formalin, undecalcified femur samples were embedded in methylmethacrylate (Fukuda et al., 2013). A region of interest was defined and bone formation rate/bone surface (BFR/BS), MAR, osteoblast surface/bone surface (Ob.S/BS), osteoclast number/bone parameter (N.Oc/B.Pm), osteoclast surface/bone surface (Oc.S/BS), erosion surface/bone surface (ES/BS), osteocyte number/bone area (N.Ot/B.Ar), and empty lacunae were quantitated using a semiautomatic analysis system (OsteoMetrics, Atlanta, GA, USA). Measurements were taken on two sections/sample (separated by ~25 μm) and summed prior to normalization to obtain a single measure/sample in accordance with ASBMR standards (Dempster et al., 2013; Parfitt et al., 1987). This methodology underwent extensive quality control and validation and the results were assessed by a research specialist in a blinded fashion.
For immunofluorescence, decalcified and cryo-sectioned samples were incubated with antibodies for CD31 (1:100, BD Pharmingen, 553370) and endomucin (1:100; Santa Cruz, sc-65495) were used as primary antibodies and Alexa Fluor 488 (1:400, Thermo, A21206) and Alexa Fluor 594 (1:400, Thermo, A11032) were used as secondary antibodies according to the manufacturer’s instructions.
CTx-I measurement
Request a detailed protocolSerum level of CTx-I was measured according to the manufacturer’s instructions (Immunodiagnostic Systems, AC-06F1).
Transmission electron microscopy
Request a detailed protocolTEM images were obtained from JEOL JEM 1400. For TEM analysis of femoral bone, whole femur was fixed with 4% paraformaldehyde in 0.08 M Sorenson’s phosphate buffer at 4°C overnight and decalcified with 10% EDTA (pH 7.2) for 10 days. Decalcified bones were embedded in epoxy resin after dehydration through an ethanol series and transferred to a transitional solvent, propylene oxide. Ultrathin sections of metaphysis area were examined for osteocyte analysis (Cheville and Stasko, 2014).
Transcriptome analysis
Request a detailed protocolEight RNA-seq samples (two GM control, two GM dKO, two OIM control, and two OIM dKO samples) were mapped to the mouse reference genome (Mus_musculus.GRCm38.80) with STAR aligner (v.2.6.1b; Dobin et al., 2013; Dobin and Gingeras, 2015). After mapping, read counts were generated by using HTSeq-count (v.0.11.3; Anders et al., 2015). The read counts were used for a differential expression analysis between control and dKO groups using DESeq2 (v.1.28.1; Love et al., 2014) with the ashr shrinkage estimator (v.2.2.47; Stephens, 2017). Statistically significantly expressed genes were determined as having absolute log-fold change larger than 1.5 and having a p value less than 0.005. For upregulated genes (log-fold change [LFC] >0; n = 152 for GM and n = 2048 for OIM) and downregulated genes (LFC <0; n = 415 for GM and n = 1,576 for OIM) in dKO samples, gene ontology (GO) analysis was performed with enrichGO function in clusterProfiler package (v.3.18.1; Yu et al., 2012) with the genome-wide annotation for mouse (org.Mm.eg.db; v.3.12) in Bioconductor. Category net plot is used for visualization using the cnet function in clusterProfiler. For gene set enrichment analysis (GSEA), differentially expressed genes in dKO RNA samples were analyzed using the GSEA software (v.4.1.0) with molecular signatures database (v.7.4) as annotation. Significant GO terms as adjusted p value <0.01 and q value <0.005 were identified for GO enrichment analysis. Nominal p value <0.01 and q value <0.25 were considered as a significantly enriched pathway with default settings for GSEA.
Analysis of cell proliferation
Cell proliferation was assessed by alamarBlue staining (Thermo Scientific, DAL1100) or Bromodeoxyuridine (BrdU) incorporation assay (Abcam, ab126556).
Collection of CM and osteoblast and endothelial cell functional assay
Request a detailed protocolCM was collected from control or dKO osteoblasts at day 6 of osteogenic culture. Briefly, cells were washed with PBS twice and incubated with serum free medium for 24 hr. Mouse EPOCs were obtained from BioChain (7030031) and cultured in GM (BioChain Z7030035). For capillary tube formation assay, EPOCs (30,000 cells/well) were seeded in a 96-well plate precoated with Matrigel (BD) and incubated with CM or FGF2. After 5 hr of incubation at 37°C, the number of tube branches in each well was quantified by counting four random fields per well as previously described (Xu et al., 2018). Alternatively, CM was added to osteogenic culture of wildtype BMSCs and osteogenic differentiation was assessed after 6 days of culture.
Measurement of glutamate and ATP production
Request a detailed protocolExtracellular glutamate release was determined according to the manufacturer’s protocols (Promega, J8021). Alternatively, intracellular ATP levels were determined by reaction of ATP with recombinant firefly luciferase and its substrate D-luciferin according to the manufacturer’s instructions (Thermo Scientific, A22066).
Analysis of mDNA copy number
The mitochondrial copy number was determined by mtDNA/nDNA using real-time PCR for mt-ND1/Hk2 or mt-16srRNA/Hk2, as previously reported (Quiros et al., 2017). Sequences are provided in Table 1.
Assessment of mitochondrial number and membrane potential
Request a detailed protocolCellular mitochondria were stained with Mito Tracker (Biotium, 70075) and detected by confocal microscopy (Leica SP8). Alternatively, flow cytometery was used to assess mitochondrial number and membrane potential of the cells treated with Mito Tracker (Biotium, 70075) and TMRE (tetramethylrhodamine ethyl ester, Biotium, 70005), respectively, according to the manufacturers’ protocols. Flow cytometry analysis was performed on an LSR II (BD Biosciences).
Measurement of OCR and intracellular ROS levels
Intracellular OCR (Cayman, 600800) and cellular reactive oxidative species (ROS) levels (C10422, Invitrogen) were measured according to the manufacturers’ protocols.
Luciferase reporter assay
Request a detailed protocolMouse calvarial osteoblasts were transiently transfected with RUNX2- (OG2-luc), β-catenin- (Topflash-luc), Smad- (3TP-luc), SP7- (SP7-luc), or ATF4-responsive (OSE1-luc) reporter gene along with Renilla luciferase vector (Promega) using the Effectene transfection reagent (Qiagen). After 48 hr, dual luciferase assay was performed according to the manufacturer’s protocol (Promega) and luciferase activities were normalized to Renilla.
RT-PCR and immunoblotting
Request a detailed protocolTo prepare bone RNA samples from mouse limbs, the hindlimbs were dissected and skin/muscle tissues were removed. The remaining tibias were chopped and homogenized. Total RNA from cells or tissues were extracted using QIAzol (QIAGEN) and cDNA was synthesized using the High-Capacity cDNA Reverse Transcription Kit from Applied Biosystems. Quantitative real-time PCR was performed using SYBR Green PCR Master Mix (Bio-Rad, Hercules, CA) with Bio-Rad CFX Connect Real-Time PCR detection system. mRNA levels were normalized to the housekeeping gene Ribosomal protein, large, P0 (Rplp0). The primers used for PCR are described in Table 1.
For immunoblotting, cell lysates were prepared in lysis buffer (50 mM Tris–HCl [pH 7.8], 150 mM NaCl, 1% Triton X-100, 1 mM dithiothreitol, 0.2% sarkosyl acid and protease inhibitor cocktail [Sigma, P8340]). Protein samples were subjected to sodium dodecyl sulfate polyacrylamide gel electrophoresis and transferred to Immobilon-P membranes (Millipore). Membranes were immunoblotted with the indicated antibodies and developed with ECL (Thermo Scientific). Immunoblotting with antibody specific to GAPDH was used as a loading control.
Statistics and reproducibility
Request a detailed protocolAll experiments were performed a minimum of two to three times. For histological staining, flow cytometry, and immunoblotting, representative images are shown. All data are shown as the mean ± standard deviation (SD). We first performed the Shapiro–Wilk normality test for checking normal distributions of the groups. For comparisons between two groups, a two-tailed unpaired Student’s t-test was used if normality tests passed and a Mann–Whitney test was used if normality tests failed. For the comparisons of three to six groups, we used one-way analysis of variance if normality tests passed, followed by Tukey’s multiple comparison test for all pairs of groups. GraphPad PRISM software (ver.9.0.2, La Jolla, CA) was used for statistical analysis. p < 0.05 was considered statistically significant.
Data availability
All data generated or analyzed during this study are included in the manuscript and supporting files. Source data files have been provided for all the figures.
References
-
Effect of rapamycin on bone mass and strength in the α2(I)-G610C mouse model of osteogenesis imperfectaJournal of Cellular and Molecular Medicine 23:1735–1745.https://doi.org/10.1111/jcmm.14072
-
Mek2 is dispensable for mouse growth and developmentMolecular and Cellular Biology 23:4778–4787.https://doi.org/10.1128/MCB.23.14.4778-4787.2003
-
Noonan syndrome, the ras-MAPK signalling pathway and short statureHormone Research 71 Suppl 2:64–70.https://doi.org/10.1159/000192439
-
MEKK2 mediates aberrant ERK activation in neurofibromatosis type INature Communications 11:5704.https://doi.org/10.1038/s41467-020-19555-6
-
Effectors of mtor-autophagy pathway: targeting cancer, affecting the skeletonCurrent Opinion in Pharmacology 28:1–7.https://doi.org/10.1016/j.coph.2016.02.004
-
WNT7B promotes bone formation in part through mtorc1PLOS Genetics 10:e1004145.https://doi.org/10.1371/journal.pgen.1004145
-
Techniques in electron microscopy of animal tissueVeterinary Pathology 51:28–41.https://doi.org/10.1177/0300985813505114
-
Decreased bone mineralization in children with noonan syndrome: another consequence of dysregulated RAS mapkinase pathway?Molecular Genetics and Metabolism 106:237–240.https://doi.org/10.1016/j.ymgme.2012.04.003
-
Neurofibromatosis in children: the role of the orthopaedistThe Journal of the American Academy of Orthopaedic Surgeons 7:217–230.https://doi.org/10.5435/00124635-199907000-00002
-
Essential role of RSK2 in c-fos-dependent osteosarcoma developmentThe Journal of Clinical Investigation 115:664–672.https://doi.org/10.1172/JCI22877
-
STAR: ultrafast universal RNA-seq alignerBioinformatics 29:15–21.https://doi.org/10.1093/bioinformatics/bts635
-
Mapping RNA-seq reads with STARCurrent Protocols in Bioinformatics 51:11.https://doi.org/10.1002/0471250953.bi1114s51
-
Identification and functional characterization of ERK/MAPK phosphorylation sites in the runx2 transcription factorThe Journal of Biological Chemistry 284:32533–32543.https://doi.org/10.1074/jbc.M109.040980
-
M-CSF, tnfalpha and RANK ligand promote osteoclast survival by signaling through mtor/S6 kinaseCell Death and Differentiation 10:1165–1177.https://doi.org/10.1038/sj.cdd.4401285
-
Mitogen-activated protein kinase pathways in osteoblastsAnnual Review of Cell and Developmental Biology 29:63–79.https://doi.org/10.1146/annurev-cellbio-101512-122347
-
ERK/MAPK signalling pathway and tumorigenesisExperimental and Therapeutic Medicine 19:1997–2007.https://doi.org/10.3892/etm.2020.8454
-
Coffin-lowry syndrome: clinical and molecular featuresJournal of Medical Genetics 39:705–713.https://doi.org/10.1136/jmg.39.10.705
-
KRAS(G12C) inhibition with sotorasib in advanced solid tumorsThe New England Journal of Medicine 383:1207–1217.https://doi.org/10.1056/NEJMoa1917239
-
Improved antitumor activity of immunotherapy with BRAF and MEK inhibitors in BRAF(V600E) melanomaScience Translational Medicine 7:279ra41.https://doi.org/10.1126/scitranslmed.aaa4691
-
Vemurafenib in multiple nonmelanoma cancers with BRAF V600 mutationsThe New England Journal of Medicine 373:726–736.https://doi.org/10.1056/NEJMoa1502309
-
Osteogenesis and aging: lessons from mesenchymal stem cellsStem Cell Research & Therapy 9:244.https://doi.org/10.1186/s13287-018-0995-x
-
Inherited CD70 deficiency in humans reveals a critical role for the CD70-CD27 pathway in immunity to epstein-barr virus infectionThe Journal of Experimental Medicine 214:73–89.https://doi.org/10.1084/jem.20160784
-
Increased glutamine catabolism mediates bone anabolism in response to WNT signalingThe Journal of Clinical Investigation 125:551–562.https://doi.org/10.1172/JCI78470
-
Wnt signaling and cellular metabolism in osteoblastsCellular and Molecular Life Sciences 74:1649–1657.https://doi.org/10.1007/s00018-016-2425-5
-
The ERK MAPK pathway is essential for skeletal development and homeostasisInternational Journal of Molecular Sciences 20:E1803.https://doi.org/10.3390/ijms20081803
-
Loss of rictor with aging in osteoblasts promotes age-related bone lossCell Death & Disease 7:e2408.https://doi.org/10.1038/cddis.2016.249
-
ERK signalling: a master regulator of cell behaviour, life and fateNature Reviews. Molecular Cell Biology 21:607–632.https://doi.org/10.1038/s41580-020-0255-7
-
Energy metabolism of the osteoblast: implications for osteoporosisEndocrine Reviews 38:255–266.https://doi.org/10.1210/er.2017-00064
-
Building strong bones: molecular regulation of the osteoblast lineageNature Reviews. Molecular Cell Biology 13:27–38.https://doi.org/10.1038/nrm3254
-
DMP1-targeted cre expression in odontoblasts and osteocytesJournal of Dental Research 86:320–325.https://doi.org/10.1177/154405910708600404
-
Female mice lacking estrogen receptor-alpha in osteoblasts have compromised bone mass and strengthJournal of Bone and Mineral Research 29:370–379.https://doi.org/10.1002/jbmr.2082
-
The ras-ERK and PI3K-mtor pathways: cross-talk and compensationTrends in Biochemical Sciences 36:320–328.https://doi.org/10.1016/j.tibs.2011.03.006
-
Bone histomorphometry: standardization of nomenclature, symbols, and unitsReport of the ASBMR Histomorphometry Nomenclature Committee. J Bone Miner Res 2:595–610.https://doi.org/10.1002/jbmr.5650020617
-
Mechanical stimulation prevents osteocyte apoptosis: requirement of integrins, src kinases, and erksAmerican Journal of Physiology. Cell Physiology 289:C633–C643.https://doi.org/10.1152/ajpcell.00278.2004
-
Analysis of mtdna/ndna ratio in miceCurrent Protocols in Mouse Biology 7:47–54.https://doi.org/10.1002/cpmo.21
-
Improved overall survival in melanoma with combined dabrafenib and trametinibThe New England Journal of Medicine 372:30–39.https://doi.org/10.1056/NEJMoa1412690
-
Germline gain-of-function mutations in SOS1 cause noonan syndromeNature Genetics 39:70–74.https://doi.org/10.1038/ng1926
-
Eotaxin (CCL11) induces in vivo angiogenic responses by human CCR3+ endothelial cellsJournal of Immunology Baltimore, Md 166:7571–7578.https://doi.org/10.4049/jimmunol.166.12.7571
-
Mechanisms of bone development and repairNature Reviews. Molecular Cell Biology 21:696–711.https://doi.org/10.1038/s41580-020-00279-w
-
Schnurri-3 regulates ERK downstream of WNT signaling in osteoblastsThe Journal of Clinical Investigation 123:4010–4022.https://doi.org/10.1172/JCI69443
-
False discovery rates: a new dealBiostatistics 18:275–294.https://doi.org/10.1093/biostatistics/kxw041
-
Regulation and metabolic functions of mtorc1 and mtorc2Physiological Reviews 101:1371–1426.https://doi.org/10.1152/physrev.00026.2020
-
Vemurafenib plus rituximab in refractory or relapsed hairy-cell leukemiaThe New England Journal of Medicine 384:1810–1823.https://doi.org/10.1056/NEJMoa2031298
-
SGK1 induces vascular smooth muscle cell calcification through NF-κb signalingThe Journal of Clinical Investigation 128:3024–3040.https://doi.org/10.1172/JCI96477
-
Rictor/mtorc2 involves mitochondrial function in ES cells derived cardiomyocytes via mitochondrial connexin 43Acta Pharmacologica Sinica 42:1790–1797.https://doi.org/10.1038/s41401-020-00591-3
-
Targeting skeletal endothelium to ameliorate bone lossNature Medicine 24:823–833.https://doi.org/10.1038/s41591-018-0020-z
-
Insulin-like growth factors: actions on the skeletonJournal of Molecular Endocrinology 61:T115–T137.https://doi.org/10.1530/JME-17-0298
-
ERK negatively regulates the epidermal growth factor-mediated interaction of gab1 and the phosphatidylinositol 3-kinaseThe Journal of Biological Chemistry 277:19382–19388.https://doi.org/10.1074/jbc.M200732200
-
ClusterProfiler: an R package for comparing biological themes among gene clustersOMICS : A Journal of Integrative Biology 16:284–287.https://doi.org/10.1089/omi.2011.0118
Article and author information
Author details
Funding
NIH-NIAMS (R21AR077557)
- Jae-hyuck Shim
AAVAA Therapeutics
- Jae-hyuck Shim
Burroughs Wellcome Fund
- Matthew B Greenblatt
NIH-NIAMS (R01AR075585)
- Matthew B Greenblatt
Novartis Institutes for Biomedical Research Global Scholars Award
- Matthew B Greenblatt
Pershing Square Sohn Cancer Research Alliance award
- Matthew B Greenblatt
The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.
Acknowledgements
We would like to thank Dr. Hwanhee Oh for help with the skeletal analysis. We also thank Oksun Lee and Jihea Kim for experimental support and the many individuals who provided valuable reagents. This project was supported by NIH-NIAMS R21AR077557 and AAVAA Therapeutics. M.B.G. holds a Career Award for Medical Scientists from the Burroughs Wellcome Fund, NIH support under award R01AR075585, a Novartis Institutes for Biomedical Research Global Scholars Award, and a Pershing Square Sohn Cancer Research Alliance award.
Ethics
This study was performed in accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All of the animals were handled according to approved Institutional Animal Care and Use Committee (IACUC) protocols (202200036) of the University of Massachusetts Chan Medical School. The protocol was approved by the Committee on the Ethics of Animal Experiments of the University of Massachusetts Chan Medical School (202200036).
Copyright
© 2022, Kim et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,302
- views
-
- 360
- downloads
-
- 16
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Medicine
Sudden death after myocardial infarction (MI) is associated with electrophysiological heterogeneities and ionic current remodelling. Low ejection fraction (EF) is used in risk stratification, but its mechanistic links with pro-arrhythmic heterogeneities are unknown. We aim to provide mechanistic explanations of clinical phenotypes in acute and chronic MI, from ionic current remodelling to ECG and EF, using human electromechanical modelling and simulation to augment experimental and clinical investigations. A human ventricular electromechanical modelling and simulation framework is constructed and validated with rich experimental and clinical datasets, incorporating varying degrees of ionic current remodelling as reported in literature. In acute MI, T-wave inversion and Brugada phenocopy were explained by conduction abnormality and local action potential prolongation in the border zone. In chronic MI, upright tall T-waves highlight large repolarisation dispersion between the border and remote zones, which promoted ectopic propagation at fast pacing. Post-MI EF at resting heart rate was not sensitive to the extent of repolarisation heterogeneity and the risk of repolarisation abnormalities at fast pacing. T-wave and QT abnormalities are better indicators of repolarisation heterogeneities than EF in post-MI.
-
- Medicine
- Microbiology and Infectious Disease
Background:
Under which conditions antibiotic combination therapy decelerates rather than accelerates resistance evolution is not well understood. We examined the effect of combining antibiotics on within-patient resistance development across various bacterial pathogens and antibiotics.
Methods:
We searched CENTRAL, EMBASE, and PubMed for (quasi)-randomised controlled trials (RCTs) published from database inception to 24 November 2022. Trials comparing antibiotic treatments with different numbers of antibiotics were included. Patients were considered to have acquired resistance if, at the follow-up culture, a resistant bacterium (as defined by the study authors) was detected that had not been present in the baseline culture. We combined results using a random effects model and performed meta-regression and stratified analyses. The trials’ risk of bias was assessed with the Cochrane tool.
Results:
42 trials were eligible and 29, including 5054 patients, qualified for statistical analysis. In most trials, resistance development was not the primary outcome and studies lacked power. The combined odds ratio for the acquisition of resistance comparing the group with the higher number of antibiotics with the comparison group was 1.23 (95% CI 0.68–2.25), with substantial between-study heterogeneity (I2=77%). We identified tentative evidence for potential beneficial or detrimental effects of antibiotic combination therapy for specific pathogens or medical conditions.
Conclusions:
The evidence for combining a higher number of antibiotics compared to fewer from RCTs is scarce and overall compatible with both benefit or harm. Trials powered to detect differences in resistance development or well-designed observational studies are required to clarify the impact of combination therapy on resistance.
Funding:
Support from the Swiss National Science Foundation (grant 310030B_176401 (SB, BS, CW), grant 32FP30-174281 (ME), grant 324730_207957 (RDK)) and from the National Institute of Allergy and Infectious Diseases (NIAID, cooperative agreement AI069924 (ME)) is gratefully acknowledged.