Murine Oncostatin M Acts via Leukemia Inhibitory Factor Receptor to Phosphorylate Signal Transducer and Activator of Transcription 3 (STAT3) but Not STAT1, an Effect That Protects Bone Mass*

Oncostatin M (OSM) and leukemia inhibitory factor (LIF) are IL-6 family members with a wide range of biological functions. Human OSM (hOSM) and murine LIF (mLIF) act in mouse cells via a LIF receptor (LIFR)-glycoprotein 130 (gp130) heterodimer. In contrast, murine OSM (mOSM) signals mainly via an OSM receptor (OSMR)-gp130 heterodimer and binds with only very low affinity to mLIFR. hOSM and mLIF stimulate bone remodeling by both reducing osteocytic sclerostin and up-regulating the pro-osteoclastic factor receptor activator of NF-κB ligand (RANKL) in osteoblasts. In the absence of OSMR, mOSM still strongly suppressed sclerostin and stimulated bone formation but did not induce RANKL, suggesting that intracellular signaling activated by the low affinity interaction of mOSM with mLIFR is different from the downstream effects when mLIF or hOSM interacts with the same receptor. Both STAT1 and STAT3 were activated by mOSM in wild type cells or by mLIF/hOSM in wild type and Osmr−/− cells. In contrast, in Osmr−/− primary osteocyte-like cells stimulated with mOSM (therefore acting through mLIFR), microarray expression profiling and Western blotting analysis identified preferential phosphorylation of STAT3 and induction of its target genes but not of STAT1 and its target genes; this correlated with reduced phosphorylation of both gp130 and LIFR. In a mouse model of spontaneous osteopenia caused by hyperactivation of STAT1/3 signaling downstream of gp130 (gp130Y757F/Y757F), STAT1 deletion rescued the osteopenic phenotype, indicating a beneficial effect of promoting STAT3 signaling over STAT1 downstream of gp130 in this low bone mass condition, and this may have therapeutic value.


Oncostatin M (OSM) and leukemia inhibitory factor (LIF) are IL-6 family members with a wide range of biological functions. Human OSM (hOSM) and murine LIF (mLIF) act in mouse cells via a LIF receptor (LIFR)-glycoprotein 130 (gp130) heterodimer. In contrast, murine OSM (mOSM) signals mainly via an OSM receptor (OSMR)-gp130 heterodimer and binds with only very low affinity to mLIFR. hOSM and mLIF stimulate bone remodeling by both reducing osteocytic sclerostin and up-regulating the pro-osteoclastic factor receptor activator of NF-B ligand (RANKL) in osteoblasts.
In the absence of OSMR, mOSM still strongly suppressed sclerostin and stimulated bone formation but did not induce RANKL, suggesting that intracellular signaling activated by the low affinity interaction of mOSM with mLIFR is different from the downstream effects when mLIF or hOSM interacts with the same receptor. Both STAT1 and STAT3 were activated by mOSM in wild type cells or by mLIF/ hOSM in wild type and Osmr ؊/؊ cells. In contrast, in Osmr ؊/؊ primary osteocyte-like cells stimulated with mOSM (therefore acting through mLIFR), microarray expression profiling and Western blotting analysis identified preferential phosphorylation of STAT3 and induction of its target genes but not of STAT1 and its target genes; this correlated with reduced phosphorylation of both gp130 and LIFR. In a mouse model of spontaneous osteopenia caused by hyperactivation of STAT1/3 signaling downstream of gp130 (gp130 Y757F/Y757F ), STAT1 deletion rescued the osteopenic phenotype, indicating a beneficial effect of promoting STAT3 signaling over STAT1 downstream of gp130 in this low bone mass condition, and this may have therapeutic value.
Oncostatin M (OSM) 2 is a member of the IL-6 cytokine superfamily with a diverse range of activities (1). It was first identified as a secreted product of macrophage-like cells that inhibited proliferation of melanoma-, neuroblastoma-, and lung cancer-derived cell lines (2), and although it suppresses breast cancer cell proliferation (3), it has been suggested to support metastasis of breast cancer in the skeleton due to proosteoclastic activities (4) and to stimulate Kaposi's sarcoma (5). Studies in knock-out mice have shown that OSM supports hematopoiesis (6) and bone formation at the expense of adipogenesis (7), but it has also been implicated in pulmonary tissue fibrosis (8), cardiac disease and repair (9), prostate cancer (10), asthma (11), periodontal disease (12), and both rheumatoid and osteoarthritis (13).
The diverse roles of OSM may be explained by its ability to act through two distinct signaling complexes. Human OSM (hOSM) signals by binding first to the ubiquitously expressed glycoprotein 130 co-receptor ␣ subunit (gp130) and then recruiting, with equal affinity, either the hOSM receptor (hOSMR) or the human leukemia inhibitory factor receptor (hLIFR) ␤ subunits (14), which may then activate distinct intracellular signaling pathways (1). Murine OSM (mOSM) behaves differently with the mOSM-gp130 complex binding with high affinity to mOSMR and very low affinity to mLIFR (15).
In the skeleton, OSM is a paracrine factor expressed by osteoblasts (7), osteocytes (7), and macrophages (16 -19). OSM stimulates the formation of bone-resorbing osteoclasts indirectly by acting on osteoblasts (20) to promote their expression of receptor activator of NF-B ligand (RANKL) (21). OSM also promotes osteoblast differentiation at the expense of adipogenesis (7,22) and stimulates bone formation in vivo (7) at least in part by acting on osteocytes (23) (matrix-embedded osteoblasts) and inhibiting their production of the Wnt antagonist sclerostin (7,23). These actions mediated by OSMR are required for normal levels of bone formation and resorption in juvenile and adult skeletons (7) and for normal skeletal response to parathyroid hormone (24), a therapeutic agent used to reduce fracture risk in osteoporosis. Surprisingly, in the absence of OSMR, mOSM can signal through LIFR to suppress sclerostin expression and stimulate bone formation (7). Unlike LIF and hOSM, when mOSM acted through mLIFR, it did not significantly stimulate RANKL production (7). This suggested a novel pathway that could increase bone mass by stimulating bone formation without promoting bone resorption.
The mechanism by which mOSM-mLIFR activates downstream signals different from those activated by canonical murine LIF (mLIF)-mLIFR signaling is unknown. This study sought to (a) identify unique downstream signaling pathways of mOSM-mLIFR compared with hOSM-mLIFR and mLIF-mLIFR, (b) determine mechanisms by which such distinct signaling pathways could be activated, and (c) ascertain the benefit of promoting the mOSM-mLIFR-specific pathway in low bone mass.

Identification of Downstream Gene Signatures by Microarray-
Microarray analysis indicated that the expression changes caused by mLIF and hOSM in Osmr Ϫ/Ϫ primary calvarial osteoblasts (defined as mLIF-mLIFR effects and hOSM-mLIFR effects, respectively) were very highly correlated (Fig. 1A), supporting previous observations that hOSM acts entirely through LIFR in murine cells (25). A similar correlation between mLIF and hOSM responses was observed in WT (data not shown). primary calvarial osteoblasts treated with mLIF or hOSM for 6 h. The genas biological correlation between the log -fold change profiles is 99.7%. B, scatterplot of log 2 -fold expression changes at 1 h in 17-day-differentiated Osmr Ϫ/Ϫ primary calvarial osteoblasts stimulated by mOSM versus control (y axis) and by the mean of hOSM and mLIF versus control (x axis). Genes highlighted in cyan are highly up-regulated by both mOSM and hOSM/mLIF. The black dashed diagonal line shows the least square regression through the origin of y on x with slope 0.367. The cyan dashed line is chosen to have slope 0.275, which is 75% of the least square slope. Gene symbols are shown for selected genes with large -fold changes. C, scatterplot of log 2 -fold expression changes at 6 h in Osmr Ϫ/Ϫ osteoblasts. Genes highlighted in cyan are highly up-regulated by both mOSM and hOSM/mLIF. The black dashed diagonal line shows the least square regression through the origin of y on x with slope 0.4. The cyan dashed line is chosen to have slope 0.3, which is 75% of the least square slope. Gene symbols are shown for selected genes with large -fold changes. D and E, real time PCR validation on three independent cultures of wild type and Osmr Ϫ/Ϫ cells treated with mOSM, hOSM, or mLIF showing Socs1 (at 1 h), Oasl2 (at 1 h), and Gbp2 (at 6 h) regulated by hOSM, mLIF, and mOSM in wild type cells (black bars), but not by mOSM in Osmr Ϫ/Ϫ cells (white bars) (D). Also shown are regulation of Socs3 (at 1 h), Bcl3 (at 6 h), and Sbno2 (at 6 h) by hOSM, mLIFR, and mOSM in Osmr Ϫ/Ϫ cells (white bars) (E). Values are means from three independent cultures, each performed in triplicate; error bars represent S.E. F, Venn diagram showing number of genes (bold text) and selected targets differentially expressed between mLIF/hOSM-treated and control in wild type cells (mLIF-mLIFR/hOSM-mLIFR) versus those regulated by mOSM treatment in wild type cells but not regulated by mOSM in Osmr Ϫ/Ϫ cells (mOSM-mOSMR). The subset of (hOSM/mLIF)-LIFR genes highly up-regulated by mOSM in Osmr Ϫ/Ϫ cells (mOSM-mLIFR) identified by the analysis in B and C are highlighted in the circle outlined with a dashed line. Ctrl, control.
For this reason, mLIF-mLIFR and hOSM-mLIFR responses were averaged for all subsequent analyses.
In contrast, in the top 10 genes most highly regulated by the canonical (hOSM/mLIF)-LIFR pathway that were not regulated by mOSM-LIFR were known STAT1/IFN␥-responsive genes (Fig. 1, B and C, and Table 2). At 1 h, the top nine known genes regulated were all STAT1/IFN␥ targets: Mx2 (37), Irgm1  (49), and Oas1g (49). A subset of genes was confirmed by quantitative real time PCR (qPCR) in a further three independent sets of cultures from Osmr Ϫ/Ϫ and WT cells (Fig. 1, D and E). This indicates that the interaction of mOSM with LIFR, unlike canonical mLIF/LIFR signaling, did not induce STAT1-responsive genes but retained the ability to induce STAT3 target genes, albeit at a lower level.
Tnfsf11 (the gene for RANKL) was significantly stimulated downstream of mOSM-LIFR in contrast to our earlier PCRbased analysis. As we reported previously, Tnfsf11 was not upregulated by mOSM in Osmr Ϫ/Ϫ cells at 1 h after treatment (7), but at 6 h after treatment Tnfsf11 was up-regulated by mOSM 2.8-fold. This response was 9.6-fold lower than the response observed to mOSM in WT cells (p ϭ 2.43 E-14) and less than the response in Osmr Ϫ/Ϫ cells treated with hOSM/mLIF (Tnfsf11 was increased 6.1-fold for both hOSM and mLIF).
Our hypothesis was that osteoblastic genes promoting osteoblast differentiation would be regulated through mOSM-LIFR, whereas genes supporting osteoclast formation (catabolic) would be regulated by mOSM through OSMR. To identify those genes regulated by only mOSM through OSMR (mOSM-OSMR) and not by mOSM through LIFR, we determined genes differentially expressed in WT cells treated with mOSM versus control that were not differentially expressed in Osmr Ϫ/Ϫ cells treated with mOSM versus control; mOSM-OSMR targets numbered 2669 genes (supplemental Table 1). This gene set was compared with those differentially expressed in mLIF/
Confirmation of LIFR-dependent STAT3 Induction-Western blotting analysis showed that mOSM, hOSM, and mLIF treatment of WT primary calvarial osteoblasts stimulated STAT1 and STAT3 phosphorylation (Fig. 2, A and B). Although Osmr Ϫ/Ϫ cells responded to hOSM and mLIF with phosphorylation of the lower molecular weight form of STAT1, no STAT1 phosphorylation was detected in response to mOSM ( Fig. 2A). STAT3 phosphorylation at tyrosine 705 was robustly stimulated by mOSM, hOSM, and mLIF in both WT cells and Osmr Ϫ/Ϫ cells at 15 and 30 min (15 min shown in Fig. 2B), indicating an intact STAT3 phosphorylation response in the absence of OSMR. To confirm that the mOSM-induced increase in STAT3 phosphorylation was mediated by LIFR, we tested the response when cells were pretreated with a specific LIFR antagonist (LA) (55). This blocked the induction of STAT3 phosphorylation by mOSM, hOSM, and mLIF in Osmr Ϫ/Ϫ cells (Fig. 2C), showing that STAT3 phosphorylation in the absence of OSMR was mediated by LIFR. In wild type cells, LA blocked the effects of hOSM and mLIF but not mOSM, indicating that mOSM acts preferentially through OSMR.
Use of a specific STAT3 inhibitor (ML116) (56) validated that the Socs3 gene response to mOSM-LIFR was STAT3-dependent (Fig. 3, A and B). In WT cells, Socs3 mRNA levels induced by mOSM and mLIF were inhibited only by the highest dose of ML116 (10 g/ml), whereas the increased Socs3 mRNA levels in Osmr Ϫ/Ϫ cells treated with mOSM were completely blocked by even the lowest dose of ML116 (Fig. 3A). Phosphorylation of STAT3 in response to mOSM in Osmr Ϫ/Ϫ cells was reduced by ML116 treatment as was the response to mLIF in both wild type and Osmr Ϫ/Ϫ cells (Fig. 3, B and C).

probes regulated by (hOSM/mLIF)-LIFR but lost by mOSM-LIFR
The top probes ordered by mLIF -fold change lying below the cyan dashed lines in Fig. 1, C and D, respectively, are shown. At 1 h, the selected genes have log 2 FC response to mOSM Ͻ27.5% of that in hOSM/mLIF. At 6 h, the selected genes have log 2 FC response to mOSM Ͻ30% of that in hOSM/mLIF. were also comparable with those induced by mLIF and hOSM in Osmr Ϫ/Ϫ and WT primary calvarial osteoblasts (Fig. 4). This suggests that the difference in genes transcribed in response to mOSM-LIFR versus mLIF-LIFR does not relate to a difference in STAT3 phosphorylation patterns.
To determine whether differences in gp130 and LIFR phosphorylation could be responsible for the bias toward STAT3 signaling of mOSM-LIFR versus mLIF-LIFR, we assessed phosphorylation of these receptor subunits in response to each ligand (Fig. 5, A-C). In WT cells, mOSM treatment resulted in robust OSMR and gp130 phosphorylation (Fig. 5, A and B). OSMR was not phosphorylated in response to mLIF (Fig. 5B) and was not detected in Osmr Ϫ/Ϫ cells as expected (Fig. 5B). LIFR phosphorylation was robustly detected in response to mLIF in WT and Osmr Ϫ/Ϫ cells (Fig. 5C). Although mLIF and hOSM both induced phosphorylation of gp130 in Osmr Ϫ/Ϫ cells, gp130 phosphorylation was barely detected in response to mOSM and in some experiments could not be detected at all (Fig. 5A). In addition, LIFR phosphorylation in response to mOSM in Osmr Ϫ/Ϫ cells was also only barely detected (Fig. 5C). This suggests that, in Osmr Ϫ/Ϫ cells, although mLIF induces phosphorylation of both gp130 and LIFR, resulting in phosphorylation of both STAT3 and STAT1, the ability of mOSM to activate STAT3 without STAT1 signaling may relate to reduced phosphorylation of gp130 and LIFR (Fig. 5D).
Assessing Possible Therapeutic Benefit in the Skeleton-To assess whether shifting the balance of signaling toward STAT3 rather than STAT1 might promote bone formation, we analyzed the bone phenotype of 12-week-old male global Stat1 Ϫ/Ϫ mice but observed no significant difference in trabecular bone volume (Fig. 6A), trabecular thickness (Fig. 6B), trabecular number (Fig. 6C), or trabecular separation (Fig. 6D). A significant increase in periosteal circumference without any change in femoral length was observed, suggesting elevated periosteal bone formation when STAT3 signaling is greater than STAT1 (mean periosteal circumference ϮS.E.: WT, 7.50 Ϯ 0.14; To further explore the implications of biased STAT3 signaling in a disease setting, we used a previously described mouse model of osteopenia caused by hyperactivation of STAT1/3 signaling downstream of gp130 by deletion of the binding site for the negative regulators SOCS3 and SHP2 (gp130 Y757F/Y757F ) (58). When these mice were crossed with Stat1 Ϫ/Ϫ mice so that only STAT3 signaling downstream of gp130 was hyperactivated (36), the gp130 Y757F/Y757F osteopenic phenotype was rescued; trabecular bone volume, trabecular number, trabecular thickness, and trabecular separation all returned to normal (Fig.  6, E-I). This was not reproduced when gp130 Y757F/Y757F mice were crossed with Stat3 ϩ/Ϫ mice, a cross that normalizes STAT3 signaling downstream of gp130 (36). Osteoclast and

Discussion
In this study, we report a full set of gene responses to mOSM in the absence of OSMR indicating that mOSM, which we have previously shown to promote bone formation independently of OSMR (7), regulates multiple gene targets downstream of mLIFR. Unlike canonical signaling of mLIF and hOSM where both STAT1 and STAT3 are phosphorylated downstream of mLIFR, mOSM induces phosphorylation of STAT3 but not STAT1 and stimulates a STAT3-responsive subset of genes downstream of mLIFR. This may be explained by formation of a signaling complex with very low levels of gp130 and LIFR phosphorylation (Fig. 5B). Supporting an anabolic action of this pathway, favoring STAT3 over STAT1 signaling downstream of gp130 protected against osteopenia in a mouse with hyperresponsiveness to IL-6 family cytokines (gp130 Y575F/Y757F ). This suggests such a pathway could be exploited in skeletal pathologies where IL-6 family signaling is elevated, such as inflammatory or metastatic bone loss. Reduced trabecular bone volume due to hyperactivation of STAT1/3 signaling downstream of gp130 in the gp130 Y575F/Y757F mice has been reported previously to result from a high level of bone remodeling (58). Although both bone resorption and bone formation are elevated, the level of resorption outstrips that of formation, resulting in an osteopenic phenotype, a similar etiology to bone loss associated with estrogen deficiency (60). Elevated signaling by gp130-dependent cytokines has been postulated not only to cause osteoporosis after menopause (61) but also to play a role in focal and systemic bone loss that results from colitis (62), inflammatory arthritis (63, 64), multiple myeloma (65), Gorham-Stout disease (66), Paget's disease (67), and breast cancer metastasis to bone (68). The osteopenic phenotype of gp130 Y575F/Y757F mice was rescued by blocking STAT1 signaling. Although bone remodeling remained high, the balance was shifted to favor bone formation, and bone mass was protected. STAT1 signaling blockade may therefore provide a therapeutic benefit in patients with bone loss related to elevated IL-6 family cytokine signaling. The greater level of periosteal bone formation in Stat1 Ϫ/Ϫ mice may relate to the loss of STAT1-mediated suppression of periosteal development; STAT1 is expressed at high levels in the developing perichondrium (69). The limited benefit of STAT1 deletion without hyperactivation of STAT3 (i.e. only periosteal bone formation was increased in Stat1 Ϫ/Ϫ mice) suggests that this approach might not benefit all forms of bone fragility.
When IL-6 family cytokines form a receptor-ligand complex, it is phosphorylation of the intracellular domain of the receptor subunits, both gp130 (70,71) and LIFR (72), that provides docking sites for STAT proteins. These are subsequently phosphorylated by JAKs on specific residues allowing homo-or heterodimerization and translocation to the nucleus to modify transcription. This study identifies that the low affinity interaction of mOSM with mLIFR results in transcription of a subset of genes downstream of canonical (mLIF/hOSM)-mLIFR signaling. These are known STAT3-responsive genes, including   OCTOBER 7, 2016 • VOLUME 291 • NUMBER 41 than STAT1, target gene transcription appeared to result from an ability of mOSM-mLIFR to phosphorylate STAT3 at both Ser-727 and Tyr-705 but not to phosphorylate STAT1 in response to mOSM in Osmr Ϫ/Ϫ cells.

Defining mOSM Action through mLIFR
The only gene that was not a classical STAT3 target in the top 10 genes was Cish, which is usually regarded as a STAT5-responsive gene (34). Previously we did not detect STAT5 phosphorylation in response to mOSM in Osmr Ϫ/Ϫ osteoblasts (7), and other STAT5 target genes, such as Pim1 and Bcl2 (73), were not in the mOSM-LIFR gene set, suggesting that Cish transcription may occur in response to a broader range of STAT proteins.
The question of how mOSM can act through gp130-LIFR to regulate only a subset of those genes regulated by mLIF/hOSM through the same receptor complex is common to all receptors with multiple ligands that produce different biological effects within a single cell type. These include not only LIFR and OSMR (1,74,75) but also type I interferons (76). In the case of mOSM, our data suggest that impaired gp130 and LIFR phosphorylation may be responsible for the bias toward STAT3 signaling, but how receptor phosphorylation level is controlled remains unknown. It may relate to altered binding conformation, binding affinity, or binding stability compared with mLIF/ hOSM. This could be better understood by structural studies as achieved with the binding of mLIF to LIFR (75). No structural data yet have defined the binding of mOSM or hOSM to either OSMR or LIFR; it remains structurally unclear why hOSM can bind to both LIFR and OSMR, whereas mOSM binds with higher affinity to OSMR (1). One possible model is that mOSM engages LIFR and gp130 in a manner that leads to incomplete transactivation of JAKs so that the STAT3 docking site on LIFR is phosphorylated, but STAT1 and SOCS3 binding sites are not. This would lead to selective STAT3 signaling and possibly its enhancement by lack of SOCS3 feedback even at low affinity interactions.
STAT3 docking sites on LIFR␤ are phosphorylated at Tyr-765, Tyr-812, Tyr-904, and Tyr-914 in mouse (Tyr-981, Tyr-1001, and Tyr-1028 in human), whereas STAT3 docking sites on gp130 are phosphorylated at Tyr-767, Tyr-814, Tyr-905, and Tyr-915 (75). Any one or a combination of these may be modified. In addition, the single SOCS3 binding site on gp130 is phosphorylated at Tyr-757 in mouse (Tyr-759 in human), and a similar SOCS3 binding site on LIFR␤ may also be phosphorylated (75). Identifying the specific differences in phosphorylation will require an extensive analysis of mutant receptors.
The focus of mOSM-LIFR action on STAT3 signaling is a similar effect on gene expression as that described previously for hOSM-hLIFR in human breast cancer cell lines (3). In that work, a LIF antagonist was used to distinguish shared and distinct signaling pathways of hOSM-hOSMR and hOSM-hLIFR. They reported that hOSM action through both hOSMR and hLIFR resulted in changes in genes downstream of STAT3 and STAT1 (e.g. Socs3 and Irf1), but regulation of STAT1 genes was lost when hOSM interacted with hLIFR. The similarities we observe suggest that mOSM in mouse cells may interact with mLIFR in a manner similar to hOSM-hLIFR in human cells. This suggests that human analogues that mimic the anabolic action of mOSM-STAT3 signaling could be developed when more information about the structure of the interactions of mOSM-mLIFR and hOSM-hLIFR becomes available. This study was conducted in primary calvarial osteoblasts because genetic deletion of OSMR in these cells allowed us to identify mOSM-dependent regulation of sclerostin (7), an inhibitor of bone formation, in the absence of OSMR. Whether mOSM is capable of initiating STAT3 phosphorylation in other OSMR-deficient cell types is not known. We could detect no STAT3 phosphorylation response to mOSM in Ba/F3 cells expressing mgp130 and mLIFR. This suggests that accessory factors that promote mOSM-LIFR complex formation may be present in calvarial osteoblasts. One possibility is sortilin, which has been shown to promote ciliary neurotrophic factor and LIF signaling (77) and which we have detected by qPCR in primary calvarial osteoblasts (data not shown). However, addition of sortilin to mgp130/mLIFR-positive Ba/F3 cells did not promote STAT3 phosphorylation.
This study began with the aim of identifying the pathway through which mOSM-mLIFR is able to promote bone formation without stimulating bone resorption. Although in our ear-lier work we observed that mOSM could promote bone formation in the absence of OSMR, and its effect on promoting RANKL (Tnfsf11) transcription was not detectable in OSMRdeficient cells, the present study detected a delayed, but still significant, increase in Tnfsf11 mRNA in Osmr Ϫ/Ϫ cells in response to mOSM. Although our hypothesis that mOSM-LIFR might promote an exclusively "anabolic" gene signature whereas mOSM-OSMR promoted an exclusively "catabolic" gene signature was not upheld, the rescue of the gp130 Y757F/Y757F osteopenic phenotype suggests that an elevated ratio of STAT3: STAT1 responses may be more anabolic than a balanced ratio or one in favor of STAT1. This may be particularly relevant when STAT3 is hyperactivated as in inflammation-or cancer-induced bone loss. This is supported by an earlier report that STAT1 inhibition promoted bone formation in fracture healing (78), the greater periosteal circumference we observed in Stat1 Ϫ/Ϫ mice, and high alkaline phosphatase activity in cultured osteoblasts with constitutive STAT3 activation (16).
To conclude, this study identifies that mOSM acting through the mLIFR phosphorylates STAT3 but not STAT1, resulting in specific regulation of STAT3-responsive genes, thereby activating a more specific intracellular signaling pathway than those induced by hOSM or hLIF acting through the same receptor. Our data suggest that this specificity relates to reduced gp130 Defining mOSM Action through mLIFR OCTOBER 7, 2016 • VOLUME 291 • NUMBER 41 and LIFR phosphorylation by mOSM-mLIFR. Targeted activation of STAT3 signaling downstream of gp130 without promoting STAT1 signaling provides a protective effect in the skeleton, suggesting that mimics of the mOSM-mLIFR interaction or of this pattern of STAT3-specific activation could provide benefit for skeletal fragility.

Experimental Procedures
Cell Culture and Microarray-Osmr Ϫ/Ϫ mice (6) were obtained from Prof. Atsushi Miyajima (The University of Tokyo), backcrossed onto C57BL/6 for six generations, and maintained as a heterozygous breeding colony at our institute for Ͼ5 years. Primary calvarial osteoblasts were generated as described previously (79) from homozygous litters of Osmr Ϫ/Ϫ and wild type cousin-bred mice. On three occasions, WT and Osmr Ϫ/Ϫ cells were defrosted, expanded, and differentiated for 17 days in osteoblast differentiation medium. At 17 days, when these cells express elevated levels of sclerostin, they were serum-starved overnight and then treated with 10 ng/ml hOSM (R&D Systems), mOSM (R&D Systems), or mLIF (Chemicon). Cell digests were collected at 1 and 6 h after treatment commenced. mRNA was prepared as described previously (7) and hybridized to Illumina MouseWG-6 V2 BeadChips. Raw probe intensities for both regular and control probes were exported from Illumina's Genome Studio software.
Bioinformatics analysis was conducted using the limma software package (80). The raw intensity values were background-corrected and normalized using the neqc function, which performs "normexp" background correction and quantile normalization using parameters estimated from the control probes (81). Probes were filtered as not expressed if they failed to achieve a detection p value of 0.05 for at least three samples. The Illumina manifest file (MouseWG-6_V2_0_R1_11278593_A) was used for probe annotation. Gene symbol aliases were converted to current official gene symbols using the Bioconductor organism package org.Mm.eg.db. Linear models were used to test for expression differences between treatments while adjusting for batch effects associated with the date on which the cell cultures were performed. Differential expression between treatments was assessed using empirical Bayes moderated t statistics (82), allowing for an intensity-dependent trend in the standard errors. The mLIF and hOSM responses were combined and compared with others using appropriate linear model contrasts. The false discovery rate was controlled at less than 0.05 using the method of Benjamini and Hochberg (59). Biological correlations between log -fold change profiles were estimated by limma's genas function, which estimates genuine associations between profiles, adjusting for any technical component of the correlation. Lists of STAT3 target genes were obtained from the Molecular Signatures Database v5.1 (83), and enrichment was tested using limma's kegga function. The microarray data are available as Gene Expression Omnibus (GEO) series GSE83418.
Phosflow and Western Blotting-Samples for Phosflow were generated by culturing mouse primary calvarial osteoblasts from C57BL/6 mice for 22-24 days in differentiation medium on three separate occasions. Cells were seeded and grown to 80 -90% confluence in culture in complete medium and serum-starved overnight in ␣MEM supplemented with 2% heat-inactivated fetal bovine serum (FBS) ϩ 50 g/ml ascorbate. The following day cells were stimulated with 50 ng/ml mOSM, hOSM, or mLIF for 15 min in low serum conditions. Cells were rinsed twice with PBS, fixed in 1% formalin at 37°C for 10 min, and permeabilized in 90% methanol for 30 min on ice. After overnight storage at Ϫ80°C, cells were washed twice in Stain Buffer (PBS plus 1% FBS, 0.09% sodium azide, 0.5 mM EDTA) and stained with the following antibodies: AF647 pSTAT3 Ser-727, AF488 pSTAT3 Tyr-705, and phycoerythrin total STAT3 (all BD Biosciences). Cells were washed, resuspended in Stain Buffer, and analyzed on an LSR Fortessa (BD Biosciences). Mean fluorescence intensity (MFI) was determined for all three biological replicates using FlowJo software. Gating was based on fluorescent signal above background such that all signals were included for cytokine induction of pSTAT3 Tyr-705 and total STAT3; signals in the secondary peak were included for cytokine induction of pSTAT3 Ser-727.
Whole cell lysates were prepared for Western blotting analysis from wild type or Osmr Ϫ/Ϫ primary calvarial osteoblasts maintained in normal growth medium (␣MEM ϩ 10% FBS) (JRH Biosciences batch numbers 074, 414, and 379) or under differentiating conditions (␣MEM ϩ 15% FBS (JRH Biosciences batch numbers 074, 414, and 379; Assay Matrix catalogue number A137A11) and 50 g/ml ascorbate (Sigma)) for 21 days in 10-cm 2 dishes. Cultures were serum-starved in 2% FBS overnight prior to cytokine treatment. Cells were treated with mOSM (R&D Systems catalogue number 495-MO, lot number TX0313011), hOSM (R&D Systems catalogue number 295-OM, lot number DY0808111), or mLIF (Merck Millipore catalogue number LIF2010, lot number 2561112) at 50 ng/ml for 15 or 30 min. LA (55) was used at 2.5 g/ml for 30 min prior to the addition of cytokines; STAT3 inhibitor (ML116) (56) was used for 1 h prior to the addition of cytokines at 0.1, 1, or 10 g/ml. Cells were then washed twice with ice-cold PBS before cell lysates were prepared by shaking at 4°C with modified radioimmune precipitation assay buffer (84), proteinase inhibitor (Sigma), and phosphatase inhibitor (Dako) before lysates were centrifuged. Protein concentration was determined according to the manufacturer's protocol (Pierce). 20 -50 g of protein were loaded onto 4 -12% gradient gels (Invitrogen) under reducing conditions for electrophoresis before being transferred to nitrocellulose membranes (iBlot, Invitrogen) and probed with pSTAT1(Tyr-701) ( Immunoprecipitation was performed by dividing the volume of lysate collected equally into three, one lysate per antibody. Immunoprecipitation was carried out using 30 l of Protein A-Sepharose fast flow beads (GE Healthcare catalogue number 17-1279-01, lot number 10168456) and 2 g of either gp130 (M-20) antibody (Santa Cruz Biotechnology catalogue number sc-656, lot number D1515), OSMR␤ (C-7) antibody (Santa Cruz Biotechnology catalogue number sc-376380, lot number L2711), or LIFR (C-19) antibody (Santa Cruz Biotechnology catalogue number sc-659, lot number F12115), rotating overnight at 4°C. Postimmunoprecipitation the beads were washed three times in PBS and once in ice-cold modified radioimmune precipitation assay buffer. Loading buffer was added, and the beads were heated at 60°C for 10 min prior to being loaded onto 4 -12% gradient gels (Invitrogen) under reducing conditions for electrophoresis before being transferred to nitrocellulose membranes (iBlot) and probed with pan-phosphorylation antibody 4G10 (Merck Millipore catalogue number 05-321, lot number 2593884). Protein bands were detected with ECL chemiluminescence and film according to the manufacturers' instructions.
Quantitative Real Time PCR-cDNA synthesis of 1 g of DNase-treated (Ambion) RNA was performed using Affini-tyScript (Agilent Technologies catalogue number 600559) according to the manufacturer's instructions. Stock cDNA was diluted 1:5-1:10, and qPCR was performed using either an inhouse master mix or 10ϫ AmpliTaq Gold with Brilliant II SYBR Green Mastermix (Agilent Technologies catalogue number 600828). Primers used (Table 3) were designed using NCBI Primer Blast or PrimerQuest. Samples were dispensed into an optically clear 96-well plate (4Titude, Thermo Scientific catalogue number 4ti-0750) and run on a Stratagene Mx3000P (Agilent Technologies) with two-step cycling conditions (95°C for 10 min followed by 95°C for 30 s and 60°C for 1 min) for 40 cycles followed by a dissociation step (95°C for 1 min, 55°C for 30 s, and 95°C for 30 s). Postrun samples were analyzed using Stratagene software MxPro and reported using linear ⌬ threshold cycle (⌬Ct) values normalized to ␤ 2 -microglobulin (B2m) or hypoxanthine phosphoribosyltransferase 1 (Hprt1). mRNA levels of these two housekeeping genes were not modified by the cytokine or antagonist treatments used.
Breeding and Analysis of Stat1-null and gp130 Y757F/Y757F Mice -Mice homozygous for the gp130 Y757F/Y757F knock-in mutation and their osteopenic phenotype have been described previously (58) as have the earlier crosses of these mice with the Stat1 Ϫ/Ϫ and Stat3 ϩ/Ϫ mice (36). Stat1 Ϫ/Ϫ mice (85) were derived from the same colony, bred from heterozygous parents, and compared with littermate wild type mice. For the cross with gp130 Y757F/Y757F mice, mice were from mixed litters including gp130 Y757F/Y757F /Stat1 Ϫ/Ϫ and gp130 Y757F/Y757F littermates and gp130 Y757F/Y757F /Stat3 ϩ/Ϫ and gp130 Y757F/Y757F littermates; wild type mice were from the same breeding colony. All mice were bred on a mixed 129/C57BL/6 background. All mice were bred and maintained in the same facility as their controls; tissues were collected for analysis at 12-14 weeks of age. Micro-CT analysis of Stat1 Ϫ/Ϫ mice was performed using a SkyScan 1076 system (Bruker micro-CT, Kontich, Belgium) on distal femoral specimens from 12-week-old mice as described previously (86). Briefly, images were acquired using the following settings: 9-m voxel resolution; 0.5-mm aluminum filter; 50-kV voltage; 100-A current; 2600-ms exposure time; rotation, 0.5°; frame averaging, 1. Images were reconstructed and analyzed using SkyScan software NRecon (version 1.6.8.0), DataViewer (version 1.4.4), and CT Analyzer (version 1.11.8.0). The femoral trabecular analysis region of interest was determined by identifying the intracondylar notch at the distal end of the femur and calculating 15% of the total femur length toward the femoral midshaft where we then analyzed a region of interest of 15% of the total femur length; a threshold of 43-255 was used to define trabecular bone. Tibial samples from gp130 Y757F/Y757F mice were analyzed by histomorphometry in the proximal secondary spongiosa using the OsteoMeasure system (Osteometrics, Decatur GA) after embedding in methylmethacrylate resin as described previously (87). Microscope images were collected on a Leica DMRB microscope coupled to an Olympus DP72 camera and CellSens software under a 4.0 Plan objective.

Gene
Primer sequence GenBank TM accession no.