Cystatin M/E ameliorates bone resorption through increasing osteoclastic cell estrogen influx

Abstract In multiple myeloma (MM), increased osteoclast differentiation leads to the formation of osteolytic lesions in most MM patients. Bisphosphonates, such as zoledronic acid (ZA), are used to ameliorate bone resorption, but due to risk of serious side effects as well as the lack of repair of existing lesions, novel anti-bone resorption agents are required. Previously, the absence of osteolytic lesions in MM was strongly associated with elevated levels of cystatin M/E (CST6), a cysteine protease inhibitor, secreted by MM cells. In this study, both MM- and ovariectomy (OVX)-induced osteoporotic mouse models were used to compare the effects of recombinant mouse CST6 (rmCst6) and ZA on preventing bone loss. µCT showed that rmCst6 and ZA had similar effects on improving percent bone volume, and inhibited differentiation of non-adherent bone marrow cells into mature osteoclasts. Single-cell RNA sequencing showed that rmCst6 and not ZA treatment reduced bone marrow macrophage percentage in the MM mouse model compared to controls. Protein and mRNA arrays showed that both rmCst6 and ZA significantly inhibit OVX-induced expression of inflammatory cytokines. For OVX mice, ERα protein expression in bone was brought to sham surgery level by only rmCst6 treatments. rmCst6 significantly increased mRNA and protein levels of ERα and significantly increased total intracellular estrogen concentrations for ex vivo osteoclast precursor cell cultures. Based on these results, we conclude that CST6 improves MM or OVX bone loss models by increasing the expression of estrogen receptors as well as the intracellular estrogen concentration in osteoclast precursors, inhibiting their maturation.


Introduction
Multiple myeloma (MM) is a malignancy of terminally differentiated B-cells that is localized primarily in the bone marrow (BM) but also can be present in peripheral blood and tissue/organs.MM cells expand in the BM, produce extra and abnormal proteins and may crowd out healthy BM cells suppressing BM function 1 .MM symptoms include hypercalcemia, anemia, renal insufficiency, and osteolysis.Osteolysis, a hallmark of MM, is the cause of severe complications seen in nearly 80% of MM cases and is the result of interactions between MM cells and the BM microenvironment leading to increased osteoclast differentiation and suppressed osteoblast differentiation, resulting in increased bone resorption and the presence of osteolytic bone lesions 2,3 .These lesions frequently cause bone pain and lead to complications, including pathological fractures, vertebral collapse, spinal cord compression, hypercalcemia, and generalized osteoporosis, compromising MM patient quality of life, impairing survival odds, and increasing treatment costs for MM patients 4 .
Bisphosphonates (BPs) are extensively utilized in clinical practice for the treatment of diseases associated with high bone resorption, such as osteoporosis, Paget's disease, and cancerinduced bone disease [5][6][7] .BPs are pyrophosphate analogs that bind to exposed bone areas of hydroxyapatite crystals 8 .During bone remodeling, they are absorbed by osteoclasts, and through inhibition of intracellular farnesyl pyrophosphate synthase as well as the suppression of GTPase prenylation and interference with downstream pathways, BPs suppress formation of osteoclasts from precursors and induce apoptosis of mature osteoclasts 9,10.Based on their chemical structure, BPs are divided into two main groups: nitrogen and non-nitrogen containing 5 .Among BPs, zoledronic acid (ZA) is the most extensively used in treating cancer-induced bone disease 11 .In subsets of MM patients, clinical trials showed that ZA combined with other novel anti-MM agents reduced skeletal related events (SREs), prolonged the period between remission and recurrence, and improved the overall survival 12,13 .However, although ZA and other anti-resorption drugs are effective at inhibiting osteoclastic bone resorption, the inability of these compounds to repair existing osteolytic lesions and the potential adverse side effects associated with long-term use of such anti-bone resorption drugs, including renal impairment and osteonecrosis of jaw (ONJ), necessitate the development of novel agents 3,14,15 .
Other anti-resorption drugs, like the RANKL monoclonal antibody Denosumab were approved by the FDA to treat MM bone disease; however, the same adverse effects seen with ZA are also found in treated patients with Denosumab 15 .Recently, we combined PET-CT scanning with global gene expression profiling of BM CD138-selected plasma cells (PC) from 512 newly diagnosed MM patients to show that the absence of osteolytic lesions is linked to elevated expression of cystatin M/E (CST6), a cysteine protease inhibitor, secreted by MM cells.
Recombinant CST6 protein inhibits the activity of the osteoclast-specific protease cathepsin K, blocks osteoclast differentiation and function, and inhibits bone destruction in ex vivo and in vivo myeloma models.Suppression of cathepsin L blocks cleavage of p100 to p52 as well as degradation of TRAF3, suppressing the alternative NF-B pathway 16 .Furthermore, Li et al.     reported that CST6 and CST6 peptides (containing the conserved QLVAG residues) inhibit breast cancer bone metastasis by suppressing cathepsin B activity 17 .
In our current study, we utilized both MM-and ovariectomy (OVX)-induced osteoporosis mouse models to compare the effects of ZA and recombinant mouse CST6 (rmCst6) on bone resorption.Single-cell RNA-seq was used to show BM cell populations in murine MM mice following treatment with either rmCst6 or ZA.Previously, it was shown that in breast cancer, loss of CST6 led to a subsequent loss of estrogen receptor alpha (ERα) 18 .Estrogen is an important regulator of bone turnover, decreasing the rate of bone resorption and increasing the rate of bone formation, and its actions on bone cells are carried out through interactions with estrogen receptors such as ERα and ERβ.We thus decided to investigate CST6 effect on intracellular estrogen concentration in osteoclast precursors on protecting against bone deterioration in both MM osteolytic bone disease mouse models and OVX mouse models.The effect of CST6 on estrogen transport and estrogen related genes in bone tissue as well as osteoclast precursors were also investigated.It is our hope that this research will assist in the development of novel bone anti-resorption drugs for the treatment of MM osteolytic lesions as well as for other bone resorption disorders such as osteoporosis.

Evaluate rmCst6 protein and ZA inhibition in MM cell-induced bone resorption in vivo
The 5TGM1-KaLwRij murine MM model was utilized to compare the effect of CST6 protein and ZA in vivo on treating MM induced bone disease.One million 5TGM1 cells were inoculated into C57BL/KaLwRij mice via the tail vein and mice were treated with purified rmCst6, ZA or PBS (Figure 1A).Intravenous (i.v) injection of purified rmCst6 protein (200 μg/kg, twice per week) and subcutaneous (s.c) injection of ZA (100 μg/kg, twice per week) improved bone compared to PBS injection, while also significantly decreasing the number of osteolytic lesions in MM-bearing mice (Figure 1B and 1C).µCT reconstruction of mouse tibia showed that rmCst6 protein and ZA significantly increased trabecular bone volume over total volume (BV/TV), trabecular number (Tb.N) and bone mineral density (BMD).However, rmCst6 or ZA treatment of MM mouse models had no significant effect on trabecular thickness (Tb.Th)     or trabecular separation (Tb.Sp) (Figure 1B and 1C).Additional data for MM mouse model µCT results are listed in Table 1.Histomorphometric analyses demonstrated that rmCst6 and ZA administration significantly reduced osteoclast (OC) numbers as well as the proportion of bone surface occupied by osteoclasts in MM-bearing mice (Figure 1D and 1E).ZA treatment caused the appearance of bone to be somewhat osteopetrotic when compared to rmCst6 (Figure 1B).To determine whether CST6 or ZA influenced MM tumor burden, flow cytometry was performed to detect the bone marrow GFP+5TGM1 cells in the tibiae at the time of sacrifice.There was no difference between MM bearing mice and MM mice treated with rmCst6 or ZA (Figure 1F and 1G).Furthermore, ELISA measurements of tumor-specific M protein, IgG2b, in serum from MM-bearing mice with or without rmCst6 and ZA therapy after 25 days revealed no difference between the control and either treatment group (Figure 1H).Finally, ELISA analyses showed that levels of the C-terminal telopeptide of type-I collagen (CTX-1), which is a biomarker of the rate of bone turnover and osteoclast activity, were significantly reduced in mice treated with ZA and rmCst6 (Figure 1I).

Evaluate rmCst6 protein and ZA effects in an ovariectomized (OVX) mouse model
We next determined if CST6 could inhibit bone loss in a murine model of estrogen deprivation-induced osteoporosis similar to ZA. Six-month old C57/BL6 ovariectomized (OVX) mice were treated with PBS, ZA or rmCst6 for 6 weeks (Figure 2A).After 6 weeks, mice were sacrificed and tibias were analyzed by CT and histology.Compared with the sham group, tibias from the OVX mice exhibited significant bone loss, and treatment with either rmCst6 or ZA appeared to improve bone quality back to sham surgery levels or better (Figure 2B and 2C).
Quantitative analysis confirmed that bone parameters, including BV/TV, Tb.N, and BMD improved in the OVX + rmCst6 mice compared to OVX mice.BV/TV, Tb.N, Tb.Sp, and BMD bone parameters improved in OVX + ZA mice compared to just OVX mice.Tb.Th was not improved in OVX mice after treatment with ZA or rmCst6 (Figure 2C).Similar to the MM mouse model, CT demonstrated that a thick band of calcified trabeculae with distorted architecture somewhat osteopetrotic under the growth plate in OVX mice treated after ZA (Figure 2B).Data for OVX mouse model µCT results are listed in Table 2. Histomorphometric analyses demonstrated that rmCst6 administration significantly reduced the number of TRAPase positive cells in OVX mice.However, after 6 weeks of ZA treatment, quantitative statistical histomorphometry showed a significant increase in TRAPase positive cells in the trabecular bone region, especially on the surface of calcified cartilage that fills the tibia metaphysis (Figure 2D and 2E).The number of adipocyte-like cells in the tibia metaphysis was also measured; there was an initial rise after OVX that was reduced to sham levels after treatment with rmCst6 and to below sham levels after treatment with ZA (Figure 2F).Finally, ELISA analysis showed that CTX-1 and P1NP levels were significantly reduced in mice treated with ZA and rmCst6 protein following OVX compared to OVX mice alone (Figure 2G and 2H).

Cell composition of MM mouse bone marrow is altered by rmCst6 and ZA treatment
Single cell RNA sequencing (scRNA-seq) was used to examine the effect of ZA and rmCst6 on the bone marrow (BM) cell composition in MM mouse models (Figure 3A).Based on expression key genes cells were sorted into specific categories.Uniform Manifold Approximation and Projection (UMAP) plot of BM mononuclear cells indicated that MM cells induced an increase in the percentage of monocyte progenitors and a decrease of B cell percentage in BM (Figure 3B and 3C).When compared to other groups, only rmCst6 treated mice showed a decrease in BM macrophage percentage and an increase in BM monocyte percentage (Figure 3B and 3C).ZA treated mice alone showed an increase in BM mature neutrophil percentage (Figure 3B).
Since it is known that macrophages are osteoclast precursors, the change in percentage of previously identified macrophage subtypes was also investigated using UMAP plots of BM macrophages following MM mouse treatment with either PBS, rmCst6 or ZA 16 .rmCst6 and ZA treatment was found to decrease percentage of M0 and M4 macrophages, identified as early precursors of osteoclasts and tumor associated macrophages with high expression of osteoclast differentiation regulators (Jun and c-Fos) 16,19 .rmCst6 and ZA also decreased the percentage of M5 macrophages, while ZA alone decreased percentage of M3 macrophages (Supplemental Figure 1).rmCst6 treatment was unexpectedly shown to increase percentage of M7 macrophages while ZA treatment also increased percentage of M2 macrophages.M3 macrophages produce inflammatory cytokines and are thought to have tumor suppressing abilities 20 .M1, M2, M5 and M7 macrophages were classified as being involved in neurological disorders and viral infections but not osteoclastogenesis 16 .
Effects of rmCST6 and ZA on the viability and differentiation of osteoclasts, osteoblasts and chondrocytes rmCst6 has previously been shown to suppress MM induced osteolytic bone disease through interfering with osteoclast differentiation and function but not viability 16 .On the other hand, ZA was shown to suppress skeletal related events in MM induced osteolytic bone disease through inactivation and apoptosis of osteoclasts 1 .To compare the effects of rmCst6 and ZA on suppressing osteoclastogenesis, mouse BM monocytes were induced to differentiate into osteoclasts by addition of M-CSF (10 ng/ml) and RANKL (10 ng/ml) in the presence or absence of rmCst6 or ZA.TRAPase staining showed that rmCst6 and ZA significantly suppressed the formation of TRAPase-positive multinuclear osteoclasts in a dose-dependent manner (Figure 4A).However, different from rmCst6, ZA appears to be more effective at promoting cell death, as no cells (osteoclasts or precursors) are present at 5 µM ZA (Figure 4A and 4B).Perhaps the ZA specific effects on osteoclastic cells might explain ZA-related osteonecrosis.
While ZA was shown to effectively suppress osteoclastogenesis through cell death there is controversial evidence suggesting ZA can impact osteoblast function 5,10 .To evaluate the effects of rmCst6 and ZA on osteoblast cell viability and differentiation, mouse osteoblast progenitor MC3T3-E1 cells were exposed to different doses of ZA and 200 ng/ml rmCst6, the most effective dose for preventing osteoclast differentiation.Alkaline phosphase (ALP), a relative early marker of osteoblast differentiation, staining was performed to assess the effects of ZA and CST6 on the MC3T3-E1 differentiation.Following 14 days treatment, ALP staining showed that ZA, but not CST6, significantly suppress the osteoblast differentiation (Figure 4C).
The mineralization capacity of cultured osteoblasts treated with ZA and rmCst6 for 21 days were also evaluated using Alizarin Red assay.Again, ZA but not rmCst6, significantly inhibited the formation of mineralized nodules in a dose-dependent manner (Figure 4C).Because of the effect of ZA on pre-osteoblast differentiation and mineralization, cell viability of MC3T3-E1 cells treated with increasing concentrations of ZA was analyzed.CCK-8 cell viability assay showed that ZA had a significant dose-dependent decrease in cell viability starting at 1 µM (Figure 4D).This data suggests that compared with rmCst6, ZA decreased the osteoblast activity and function.
In addition to osteoblastogenesis, endochondrogenesis is also critical for bone formation.
To assess the effect of rmCst6 or ZA on chondrogenesis, teratocarcinoma stem cell line ATDC5 was utilized to determine the effect of both treatments on pre-chondrocyte differentiation, mineralization, and cell viability 21 .Following 14 days of treatment in chondrocyte differentiation medium, alcian blue staining demonstrated that the glycosaminoglycan (GAG)-rich extracellular matrix (ECM) found during chondrogenesis was suppressed by ZA treatment but not rmCst6 (Figure 4E).ATDC5 cells can mineralize surrounding ECM to produce mineral nodules.As such, alizarin red staining was performed and it was shown that ZA treatment and not rmCst6 inhibited ATDC5 mineralization (Figure 4E).Finally, CCK-8 cell viability assay showed that, similar to osteoblast precursors starting at 2 μM ZA, ATDC5 cell viability was significantly decreased in a dose-dependent manner (Figure 4F).

rmCst6 and ZA bring OVX induced inflammatory cytokine levels back to control
Inflammation is a key factor in osteoclastogenesis, as such the anti-inflammatory effects (and thus the anti-bone resorptive effects) of the treatments used in this study (rmCst6 and ZA) as well as the proinflammatory effects of OVX were measured using an inflammatory cytokine array Raybiotech.The 25 targets with the largest initial change in membrane intensity from the Sham surgery group to the OVX-PBS group are shown (Figure 5).The anti-inflammatory effects of rmCst6 and ZA appear to be equivalent.Data and statistics for 25 targets with the largest initial change in membrane intensity following OVX are listed in Table 3.The complete data set for inflammatory cytokine array is listed in supplemental file 1.

ZA and rmCst6 treatment upregulate different genomic pathways in bone of OVX mice
For the OVX model it was shown that both rmCst6 and ZA ameliorated OVX induced inflammation (and thus potentially bone resorption) back to sham surgery levels.To determine the mechanisms by which rmCst6 and ZA suppress bone resorption, sequencing analysis of RNA isolated from the tibia of OVX mice treated with either rmCst6 or ZA was performed (Figure 6A).Between the OVX-PBS and the OVX-rmCst6 treatment groups, there were 175 differentially expressed genes unique to the rmCst6 treatment; between the OVX-PBS and the OVX-ZA treatment groups, there were 12 differentially expressed genes unique to the ZA treatment and there are 3 differentially expressed genes shared between both OVX-PBS vs OVX-rmCst6 and OVX-PBS vs OVX-ZA (Figure 6B).The 10 genes with the largest fold compared to OVX-PBS are listed for both the rmCst6 and ZA treatment (Figure 6C and 6D).Of particular interest, for rmCst6 treatment, the solute carrier organic anion transporter family member 1a4 (Slco1a4), an organic anion transporter, had the largest fold increase (log2 fold = 3.60), while for ZA treatment, MMP9, which plays a role in apoptotic pathways, had the greatest fold increase (log2 fold = 3.33) 22,23 .Since ZA, is known for apoptotic effects on bone cells, the effect of rmCst6 on Slco1a4 gene expression was investigated further.Basal Slco1a4 RNA levels were found for different bone tissues and cells using PCR and agarose gel electrophoresis (Tissue: tibia, vertebrae, isolated bone marrow cells and empty femur; Cell types: Raw 264.7,ST2, OB6, ATDC5 and MLOY4) (Figure 6E and 6F).Real-time PCR for RNA isolated from OVX mice L3-L5 vertebrae matched RNA-seq results, mRNA levels were decreased following OVX and only rmCst6 treatment brought Slco1a4 mRNA back to sham surgery levels (Figure 6G).To determine cell types where rmCst6 treatment increases Slco1a4 gene expression realtime PCR was performed for macrophage Raw 264.7 cells (+/-20 ng/ml RANKL) and stromal ST2 cells (+/-1 mM A2P).For Raw 264.7 cells, only rmCst6 (+ 20 ng/ml RANKL) displayed a dose dependent increase in Slco1a4 mRNA levels up to 100 ng/ml rmCst6 (Figure 6H).rmCst6 treatment also increased Slco1a4 mRNA levels in osteoblast precursor ST2 cells (+/-1 mM A2P), however based on one-way ANOVA changes were not considered significant (Figure 6I and 6J).Normalized, log(2) gene expression for RNA-seq data is listed in supplemental file 2.

rmCst6 treatment increases intracellular estrogen concentration of osteoclast precursors
Due to structural and sequence similarities in the solute carrier organic anion transporter family, Slco1a4 is predicted to play a role in Na + independent transport of estrogen and its derivatives across the plasma membrane into the cell [24][25][26] .Previously it was found that loss of CST6 in breast cancer led to a loss of ERα expression 18 .As such, we investigated the effect of rmCst6 treatment on estrogen transport and its downstream effects.Protein expression of ERα in L3-L5 was investigated for mice subject to sham surgery or OVX and treated with ZA or rmCst6.When compared to sham surgery mice, ERα protein levels are decreased for OVX mice and brought back to near sham levels only by treatment with rmCst6 but not ZA (Figure 7A).
ERα mRNA and protein levels were analyzed at the cellular level using Raw 264.7 cells (+/-20 ng/ml RANKL).For Raw 264.7 cells incubated in the presence and absence of RANKL, treatment with rmCst6 significantly increased ERα mRNA levels compared to control.For Raw 264.7 cells incubated in the absence of RANKL, 1 µM ZA also increased ERα mRNA levels.
ERα protein levels were also increased for Raw 264.7 cells (+/-20 ng/ml RANKL) (Figure 7B and 7C).Protein band intensity data quantified by imageJ and normalized to β-actin intensity are listed in Table 4 and Table 5.Total intracellular estrogen (Estrone; E1, 17-β-estradiol; E2 or Estriol; E3) concentration was evaluated in non-adherent mouse bone marrow cells (pre-osteoclasts) isolated from the femur following 72 hr.treatment with either ZA or 100 ng/ml rmCST6.Intracellular total estrogen concentration was significantly increased following treatment with 100 ng/ml rmCst6.
Treatment with ZA did not significantly alter intracellular estrogen concentration compared to control (Figure 7D).Immunohistochemistry was used to investigate alterations in estrogen responsive genes in Raw 264.7 cells treated with 100 ng/ml rmCST6.At 24 hr.incubation, ERα expression was increased following rmCst6 treatment.PPARγ expression, which is suppressed by estrogen, was lowered following rmCst6 treatment (Figure 7E).Raw 264.7 cells were incubated for 7 days with 20 ng/ml RANKL to determine the rmCst6 affects the expression of estrogen responsive genes in macrophage-like cells differentiating into osteoclasts.Similar to the 24 hr.incubation period, treatment with 100 ng/ml rmCst6 caused an increase in the expression of ERα and a decrease in the expression of PPARγ (Figure 7F).We validated the remarkable influx of estrogen into Raw 264.7 cells after treatment with 100 ng/ml rmCst6 for 24 hours (Figure 7G).

Discussion
In this study, we show that rmCst6 treatment recovers bone mass in both MM and OVX mouse models without the increased trabecular thickness comparable to the commonly prescribed bisphosphonate ZA. rmCst6 was shown to be as effective as ZA at suppressing osteoclast cell proliferation and viability without cytotoxic effects towards osteoblast and chondrocyte cell lines.When compared to ZA, rmCst6 was shown to significantly increase intracellular estrogen concentration in mouse bone marrow cells.In both male and female bone, estrogen plays an important role in regulating bone turnover, impacting both osteoclast and osteoblast function through a variety of different pathways [27][28][29][30] .As such, biomolecules with the ability to regulate estrogen influx in bone cells would be highly prized as novel agents for combating diseases associated with increased bone resorption.
Based on UMAP results for MM mice, estrogen may attenuate monocyte polarization into macrophage cell type through suppressing expression of key inflammatory factors, explaining the rmCst6 induced increase in monocyte cells and decrease in macrophages 16,[31][32][33] .A depletion in macrophage cell percentage would in turn lead to a depletion in osteoclasts and thus suppression of bone resorption.Both rmCst6 and ZA treatments were also shown to affect the different macrophage subtypes present in the bone marrow.M0, M4 and M5 macrophage subtypes were significantly decreased in both rmCst6 and ZA treated MM mouse models.M0 macrophages are undifferentiated macrophages with potential to polarize into different macrophage sub-types, including osteoclasts 19 .The specific role of M4 macrophages is currently unknown, however this subtype has previously been shown to have increased expression of osteoclast differentiation genes 16 .As with the decrease in total BM macrophage cell percentage the decrease in M0 and M4 macrophage subtypes following rmCst6 treatment can potentially be explained as estrogen preventing monocyte polarization to macrophage cell types and macrophage polarization into osteoclasts 16,[31][32][33][34] .For ZA treatment, decrease in M0 and M4 macrophages may be from known apoptotic effects of ZA 22 .In the context of MM osteolysis as well as other bone resorption disorders, suppression of macrophage cell percentage, as well as a decrease in M0 and M4 macrophage subtypes would explain the bone protective effects of both rmCst6 and ZA.
Currently the decrease in M5 and M3 macrophage following ZA treatment can be explained in the context of MM.For ZA, apoptosis may again explain the decrease M5 and M3 macrophages.Potential apoptosis of anti-tumor acting M3 macrophages hint that ZA alone may not be most effective at treating MM osteolytic bone disease 20,22 .For rmCst6 treatment suppression of M5 macrophages and increased expression of M7 macrophages are more difficult to explain.As M5 and M7 macrophages are thought to play a role in the immune system, Cst6 may impact M5 and M7 macrophage levels through inhibition of key cysteine proteases involved in immune system regulation 35 .However, will be needed to clarify the precise mechanism explaining rmCst6 induced increase in M5 and M7 macrophages.M2 macrophage involvement in oxidative stress repair pathways may explain their increase following ZA treatment due to increased ROS 22,36 .
To further investigate the ability of rmCst6 as an inhibitor of bone resorption, rmCst6's impact on osteoclastogensis on mouse primary non-adherent BM cells was compared to ZA. rmCst6 was found to have a dose dependent effect on osteoclast number while ZA, even at low concentrations significantly depleted number of osteoclast and precursor cells.In addition to osteoclasts, ZA was also found to have significant suppressive effects on cell vitality of osteoblast (MC3T3-E1) and chondrocyte (ATDC5) precursors hinting towards toxicity from long term ZA use 37,38 .In contrast to ZA, rmCst6 at the highest concentration available (200 ng/ml) had no cytotoxic effects on osteoblast and chondrocyte cell vitality, further hinting towards its potential as an alternative to ZA.However, the efficacy, bioavailability and long-term safety of rmCst6 compared to ZA will need to be addressed in future studies.
In the OVX model, both rmCst6 and ZA were shown to decrease OVX induced increase in inflammatory cytokines back to sham surgery levels.While inflammatory cytokine array was performed for only OVX model, in the context of sex steroid deficiency as well as MM induced bone loss, suppression of inflammatory factors would suppress osteoclastogenesis and thus bone resorption.Sex steroid deficiency is known to cause an increase in adipocyte-like cells present in the bone marrow, explaining increase in adipocyte like cells from OVX mouse model histology results.This increase in adipose tissue would in turn cause an increase in inflammatory factors and thus bone resorption 39 .CST6's ability to promote estrogen efflux may explain this decrease in adipose tissue 40 .In addition, estrogen has been shown to suppress expression of inflammatory cytokines such as TNFα, BLC (CXCL13) and SCF [41][42][43] .ZA was also shown to suppress adipocyte-like cell formation and expression inflammatory factors.In the context of the OVX model, ZA induced effect on inflammatory factors may also be through suppression of adipose tissue formation.Currently the mechanism by which ZA reduces adipogenesis is unknown.However, based on previous results for osteoblast and chondrocyte cell lines, high concentrations of ZA may also promote apoptosis of adipocyte-like cells 22 .ZA was also shown to significantly increase the number of TRAP positive cells in BM of OVX mice tibia.Increased number of TRAP positive cells and mRNA levels of MMP9 may be due to ZA causing an increase in the number of "non-attached" osteoclasts undergoing prolonged apoptosis 22,44,45 .In the future the impact of rmCst6 as well as ZA on the expression of inflammatory cytokines in MM mice will need be studied.
The mechanism by which rmCst6 promotes estrogen influx in macrophage precursors is currently unknown.Based on OVX model RNA-seq data, expression of the organic anion transport protein Slco1a4 is increased following treatment with rmCst6.In humans, the Slco1a4 gene analog includes both estrone-3-sulfate and 17β-estradiol-glucuronide as substrates [24][25][26] .
While the OVX mouse model is "sex steroid deficient" it has been shown that levels of estrone (E1), are higher in post-menopausal women, potentially explaining the bone protective ability of rmCst6 in this model 46 .Based on real-time PCR results from cell lines, upregulation of Slco1a4 following rmCst6 treatment may partially explain the increase in intracellular estrogen concentration in bone marrow primary cell lines.While Raw 264.7 and ST2 cell lines can differentiate into mature osteoclasts and osteoblasts, respectively, they are not primary precursor cells, and as such, the expression of Slco1a4 will need to be determined in primary osteoclast or osteoblast precursor cells isolated from bone marrow.It will also need to be determined if Slco1a4 is upregulated in MM mice following rmCst6 treatment.To establish Slco1a4 as a factor in rmCst6 induced estrogen influx, further research is needed.
Cystatins are not thought to make the cell membrane more permeable to biomolecules such as estrogen.To explain the increase in intracellular estrogen levels in osteoclast precursors following rmCst6 treatment, estrogen transport would need to be increased from upregulated expression of plasma membrane localized ERα (and potentially Slco1a4).However, a mechanism explaining how rmCst6 increases ERα expression is currently unknown.Previously it was proposed that ERβ may cause an upregulation of cystatins in triple negative breast cancer 47 .These cystatins may suppress canonical TGFβ signaling leading to decreased Smad2/3 phosphorylation 48 .Inactivated Smad2/3 decreases transport of Smad4, a known inhibitor of ERα transcriptional activity into the nucleus 49 .Increased ERα activity may lead to a mechanism whereby ERα regulates its own expression, leading to increased estrogen influx.
We hope that the work performed in this study will assist in the development of novel treatments focused on using recombinant CST6 to ameliorate bone resorption in MM osteolytic bone disease as well as other condition induced increased bone resorption.Here we found that increased intracellular estrogen influx in pre-osteoclastic primary BM cells treated with rmCst6 support a mechanism by which Cst6 mediates increased intracellular uptake of estrogen, leading to decreased bone resorption and increased bone formation [27][28][29][30] .Previously, loss of Cst6 was shown to negatively impact ERα gene expression in breast cancer.Yet, the mechanism explaining how Cst6 promotes estrogen influx is currently unknown.It is our belief that this research will lead to further studies focused on understanding the mechanism by which Cst6 is able to promote estrogen influx and suppress bone resorption.Our future efforts to understand the mechanism by which CST6 suppresses bone resorption will be investigated through use of an ERα conditional knockout mouse model and by investigating TGFβ/Smad signaling pathways that CST6 may potentially regulate.
In ex vivo, mouse bone marrow cells from flushed femur were first cultured in αminimum essential media (αMEM) supplemented with 10% FBS, 1% P/S, and 4 mM glutamine for 24 hrs.Following this, non-adherent bone marrow cells were cultured onto tissue culture plates for analysis of intracellular estrogen concentration.

5TGM1/KaLwRij MM mouse model
Six to eight-week-old female C57BL/KaLwRij mice (Harlan Mice, Netherlands) were injected with 1 × 10 6 5TGM1-GFP cells intravenously via the tail vein and randomized into 3 groups (MM mice, n = 5 per group).Five days after injection of tumor cells, mice were treated with either PBS (100 µl), rmCst6 (200 μg/kg) via intravenous (i.v) injection twice a week, or Zoledronic Acid (ZA) (100 μg/kg) via subcutaneous (s.c) injection twice a week.On day 25 post-tumor cell inoculation, when most mice had started to develop paraplegia, the experiment was terminated, and mice were sacrificed.Blood samples were collected every week.All animal procedures adhered to a protocol approved by the local Institutional Animal Care and Use Committee (IACUC) at the University of Arkansas for Medical Sciences.

Micro-computed tomography (µCT)
µCT analysis of mice tibia was performed as previously described 50 .Tibias from both MM and OVX models were dissected and fixed in 10% neutral-buffered formalin for 2 days.For µCT, region of interest (ROI) was selected to include the entire epiphysis and metaphysis of one end of the tibia as it contains trabecular bone.Scans were acquired at 60kV and 166uA; Al 0.5mm filter; 10uM Pixel size.After scanning, tibia images were reconstructed using the Skyscan NRecon program with a beam hardening correction of 40.Trabecular and cortical bone microarchitecture were analyzed using the Skyscan CT Analyzer program.The following bone parameters were calculated using µCT analysis: bone volume fraction (BV/TV%), trabecular number (Tb N, 1/mm), trabecular thickness (Tb Th, mm), trabecular separation (Tb Sp, mm), and bone mineral density (BMD).

Bone histomorphometry
Following µCT, tibia from both MM and OVX mouse models were decalcified in 5% EDTA solution (pH 7.0) for 7 days at room temperature and embedded in paraffin.Bone sections (5 μm thickness) were stained for H&E and tartrate-resistant acid phosphatase (TRAPase) staining using a Leukocyte Acid Phosphatase Kit (Sigma-Aldrich, St. Louis, MO).For the MM mouse models, histomorphometric analyses was performed using the OsteoMeasure software (OsteoMetrics,) with a Zeiss Axioskop2 microscope.For the OVX mouse models, TRAPase positive cells and adipocyte-like cells were counted using Nikon NIS elements platform and normalized to total area of ROI at 20X magnification (285474.816µm 2 ).At 10X and 20X magnifications, regions were selected to include areas with minimal damage.In addition to this, regions were also chosen to include all cell and tissue types (adipocyte-like cells, TRAPase positive cells and bone matrix).For both MM and OVX mouse models, histomorphometric parameters were averaged for each treatment and graphed using GraphPad prism 9.0.

Evaluation of intracellular bone turnover markers and estrogen levels
The serum levels of CTX-1 and P1NP were examined using a CTX-1 ELISA kit and P1NP ELISA kit (Novus) according to the manufacturer's instructions for both OVX mice.Intracellular total estrogen concentration for non-adherent BMC treated with rmCst6 (100 ng/ml) or ZA (1 µM) was examined using a mouse estrogen ELISA kit (Mybiosource: MBS766177).
Cells were cultured in αMEM in triplicate with treatments for 72 hrs. at 37° C.After 72 hrs., cells were lysed and assayed according to the manufacturer's instructions.Total intracellular estrogen concentration for non-adherent BMC was measured using triplicate analysis and experiment was repeated three times to give N = 9 for each treatment (Control, 1 µM ZA, 100 ng/ml rmCst6).Results for estrogen concentration (pg/ml) were normalized by dividing by intracellular protein concentration (µg/ml).To reduce variability between plates, data was normalized for each plate by dividing estrogen concentration (pg/µg) for each sample by average estrogen concentration of control treatment for said plate, then multiplying by the average estrogen concentration for all samples across three plates.Following normalization, results were plotted using GraphPad Prism 9.0.

Single cell RNA sequencing (scRNA-seq)
scRNA-seq was performed as previously described 16 .Briefly, BM mononuclear cells from 5TGM1/KaLwRij (MM) mice were isolated were isolated at 25 days post tumor cell inoculation.Cells were isolated from healthy control mice (n = 2), MM mice (n = 3), MM mice + rmCst6 (n = 3) or MM mice + ZA (n = 2) treated mice sorted using fluorescence-activated cell sorting (FACS).Sorted GFP negative cells with a purity greater than 95% and viability higher than 95% were used for 10 X genomics scRNA-seq, resuspended in Mg-and Ca-free PBS (+0.04% BSA) and counted using a light microscope under 10X magnification.Single-cell emulsions were generated (Chromium Controller, Chromium Next GEM Chip G, Chromium NextGEM Single Cell 3' v3.1 kit; 10X Genomics, Pleasanton, CA, USA).Libraries were assessed for mass concentration, fragment size and validated.Initial sequencing was performed on an Illumina NovaSeq SP 100-cycle flow cell and data was assessed using the Cell Ranger Count output.

Bioinformatic analysis of scRNA-seq
Bioinformatics analysis of scRNA-seq data was performed as previously described 16 .
To visualize genes simultaneously in kernel joint density estimation, the Nebulosa package was used.Based on the kernel joint density of Adgr1 and Fcgr3, we sub-selected cells that have a high value of the kernel joint density for subclustering analysis to study the cellular heterogeneity of macrophage cells.Gene-set enrichment analysis of marker genes was performed on Gene Ontology annotation using piano package.

In ex vivo osteoclast differentiation
Primary mouse BM cells were flushed out from the femur and tibia 6-8-week C57BL/6J mice, BM erythrocytes were removed with the ACK Lysing Buffer (KD Medical, catalog # RGF3015) for 3 min.The remaining cells were cultured in α-MEM (Gibco, catalog # 12571048) supplemented with 10% FBS and 10 ng/ml M-CSF (PeproTech, catalog # 315-02) overnight in 10 cm tissue culture dish. 4 x 10 4 non-adherent cells were seeded into 96-well plate with α-MEM containing 10% FBS and 10 ng/ml M-CSF for 3 days to recruit macrophage and then osteoclast differentiation was induced by addition of 10 ng/ml RANKL (R&D, catalog # 462-TEC-010/CF) with or without rmCst6 protein or ZA for 3-5 days.Media changes were carried out every 2 days.After 3 days the cells were evaluated for TRAP staining.
On day 14, the activity of alkaline phosphatase (ALP) was evaluated using histological staining.
Cells were washed with PBS and fixed with 4% paraformaldehyde for 20 min at room temperature, and then stained with 1-step NBT-BCIP substrate solution (Thermo Scientific, catalog # 34042) for 20 minutes.After staining, cells were washed twice with water and observed under a light microscope.On day 21, cells were fixed and stained with 2% Alizarin Red (Sigma-Aldrich, catalog # A5533) at pH 4.2 to evaluate the mineralization using previously described methods 51 .

Cell counting kit (CCK)-8 assay
2,000 cells (either ATDC5 or MC3T3-E1 cell line) seeded in 96 well plate in 100 μl culture medium were allowed to adhere overnight and then were treated with different doses of Zoledronic Acid (0-20 μM) for 7 days in a humidified atmosphere of 5% CO2 at 37˚C.Following treatment, 10 μl CCK-8 solution (Apexbio Technology, catalog # K10181) was added to each well and the plate was incubated for 2 hours in the incubator.The absorbance was measured at a wavelength of 450 nm using a microplate reader (BioTek, USA).Results were averaged and plotted using GraphPad Prism 9.0.

Cytokine array
RayBio C-Series Mouse Cytokine antibody array C1000 (RayBiotech) was used according to the manufacturer's protocol with protein lysate from ovariectomized mice L3-L5 vertebrae.Protein was isolated from vertebrae by homogenization with RIPA buffer (Solarbio) followed by sonication.Protein concentrations were found for lysates using BCA assay and for each treatment samples were pooled.Final concentration for each pooled sample was 1.94 mg/ml.Samples were added to array membrane in triplicate according to the manufacturer's protocol.A heat map (GraphPad Prism 9.0) was created for array targets where CST6 or ZA treatments were shown to reduce cytokine levels elevated by OVX to near sham levels.

RNA isolation, real-time reverse transcription-polymerase chain reaction
For measurement of Slco1a4 mRNA levels, Raw 264.7 cells (20 ng/ml RANKL) or ST2 cells were cultured into 12 well tissue culture plates (2.5 x 10 5 cells/ml, 1 ml total volume).Cells were treated with either a PBS control solution or increasing concentrations of rmCst6 (50 ng/ml, 100 ng/ml or 200 ng/ml).Cells were incubated at 37° C for 72 hr.For ST2 cells incubated with ascorbate-2-phosphate (A2P), cells were initially cultured at 2000 cells/well and incubated at 37°C for 10 days, with the same treatments applied.Media and treatments were replaced every 48 hr.For measurement of ERα mRNA levels Raw 264.7 cells were cultured into 12 well tissue culture plates (2.5 x 10 5 cells/ml, 1 ml total volume), treated with PBS control solution, 1 µM ZA, or 100 ng/ml rmCst6 (+/-20 ng/ml RANKL) and incubated at 37° C for 48 hr.
For each treatment, triplicate samples were cultured.RNA was isolated using RNeasy plus Mini Kit from Qiagen.RNA concentration and purity (A260/A280) for RNA samples were determined using a Polarstar Omega plate reader.Reverse transcription was carried out using an iScript cDNA Synthesis Kit.Real-time polymerase chain reaction (RT-PCR) experiment was carried out using SYBR Green and with the QuantStudio 6 Flex real-time PCR system from Applied Biosystems.Samples were performed in triplicate, and gene expression was normalized to cyclophillin A and averaged for each treatment concentration.ERα, Slco1a4 and cyclophillin A primers for RT-PCR analysis were designed using Primer Express Software 2.0.0.
For analysis of basal Slco1a4 tissue and cell mRNA levels RNA was isolated from the following tissues using the trizol method: tibia, L3-L5 vertebrae, bone marrow cells, and empty femur as well as from the following cell types Raw 264.7,ST2, OB6, ATDC5 and MLOY4 and purified using RNeasy plus Mini Kit.RNA concentration and purity (A260/A280) for RNA was determined as previously described and cDNA was synthesized.QuantStudio 6 Flex real-time PCR system was used for determining basal expression in tissues and cell culture.Following real-time PCR, gene expression was then visualized using agarose gel electrophoresis.
For in vivo analysis of Slco1a4 and ERα mRNA levels total RNA was isolated from mouse L3-L5 vertebrae for all samples from each treatment group (Sham, OVX-PBS, OVX-CST6, OVX-ZA).Tissue was first homogenized using metal beads and a homogenizer (6500 rpm, 20 seconds, twice).After homogenization, total RNA was isolated using trizol and purified using the RNeasy plus Mini Kit.For real time PCR, RNA concentration and purity was determined as with previous RNA samples and cDNA was synthesized using using iScript cDNA Synthesis Kit.RT-PCR experiment was carried out using SYBR Green as described in the above paragraphs with gene expression being normalized to cyclophillin A. For all RT-PCR experiments, results were averaged and graphed using GraphPad Prism 9.0.

Bone tissue RNA sequencing analysis
For RNA sequencing experiments, RNA was isolated from in vivo (tibia; OVX mouse model; Sham, OVX-PBS, OVX-CST6, OVX-ZA; three samples per group) samples as described earlier.R12NA integrity was checked using Agilent Technologies 4200 Tapestation with RNA screen tape.Samples with RIN e greater than 8 were used for further analysis.RNA sequencing was performed in the UAMS developmental genomics core.Briefly, Sequence ready libraries were constructed using the Illumina Stranded mRNA Prep Ligation kit following the manufacturer's protocol, then sequenced on an Illumina NextSeq 500.Briefly, poly-A containing mRNA was isolated from total RNA using oligo-dt attached magnetic beads and then converted to cDNA.After fragmentation and end repair, cDNA 3' ends were adenylated and then ligated with index anchors.The anchor ligated fragments were indexed and amplified, then normalized, pooled and sequenced.

Analysis of RNA sequencing results was performed at the Louisiana Cancer Research
Center.Briefly, FASTQ files were uploaded to Partek Flow and contaminants were removed with Bowtie 2 (v2.2.5).Reads were aligned to STAR v2.7.8a using the mm10 version of the mouse genome.Aligned reads were quantified using Ensembl Transcripts release 102.Features with a maximum reads ≤ 5 were filtered out from the analysis.Normalization was done by TMM and log2 transformation.Pathway analysis (KEGG) was performed in Partek Flow.Heat maps comparing LS mean for treated samples (OVX-CST6, OVX-ZA) to OVX samples were created using GraphPad Prism 9.0.

Western blots
To analyze ERα protein expression, total protein from each sample was first isolated from L3-L5 vertebrae by homogenization with RIPA buffer (Solarbio) followed by sonication.Protein concentrations were found for lysates using BCA assay.For ERα protein expression in Raw 264.7 cells (+/-20 ng/ml RANKL) cells were cultured onto 6 well tissue culture plates with either control, 1 µM ZA or 100 ng/ml rmCst6 treatments and incubated at 37° C for 48 hr.(2.5 x 10 5 cells/ml, 2 ml total volume).After 48 hr. cells were lysed using RIPA buffer and protein concentration was found using BCA assay.Both protein tissue and protein cell samples were prepared for SDS PAGE, equal amount of protein were loaded for each sample and then ran on 10% Acrylamide/Bis-acrylamide gels (N = 6 for bone tissue, N = 3 for Raw 264.7 cells).
Following membrane transfer, mouse anti h/m/r ERα (MAB57151, R&D Systems) was used as the primary antibody and HRP conjugated anti-mouse IgG (R&D systems, HAF018) was used as the secondary antibody.Bands of interest were visualized and imaged under chemiluminescent detection using the Amersham Imager 600 System.Amido black staining (for tissue lysate samples) and mouse anti β-actin antibody (Sigma, A1978) (for cell lysate samples) were used as loading controls for western blots.Densitometry analysis of western blot bands was performed using ImageJ.

Immunohistochemistry
Raw 264.7 cells were grown to 10,000 cells/well in Nunc.Lab Tec Chamber slides (Thermo Fisher, 177399) in DMEM supplemented with 10% FBS and 1% P/S.Cells were treated with either control solution or 100 ng/ml rmCst6 and incubated at 37° C for either 24 hr.or 7 days.7-day incubation period was in the presence of 20 ng/ml RANKL.Following incubation, immunostaining was performed as follows.Media was aspirated and cells were rinsed with cold 4% paraformaldehyde in 1X PBS and incubated at room temperature for 1 min.
Cells were then fixed with 4% paraformaldehyde in 1X PBS at room temperature for 20 min.
Final cells were covered with DAPI-Fluoromount-G and observed using Nikon Eclipse T/2 epifluorescent microscope.

Statistical analysis
Statistical analysis was performed with GraphPad Prism 9.0 (GraphPad Software, Inc., San Diego, Ca, USA).Numerical variables were expressed as mean ± SD (Standard Deviation); n equals to the number of samples/group.For all experiments differences within groups were           Protein level significance was analyzed using one-tailed t-test (Table 4).(C) Raw cells treated with 1 µM ZA or 100 ng/ml of rmCst6 for 48 hr.(20 ng/ml RANKL) had increased mRNA and protein levels of ERα when compared to control treated samples (N = 3).Protein level significance was analyzed using one-tailed t-test (Table 5).

Figure 1 :
Figure 1: rmCst6 protein and ZA inhibit bone destruction in 5TGM1-C57BL/KaLwRij MM mice.(A) A schematic model for the MM mouse study.5TGM1 murine MM cells were injected into 8-week-old C57BL/KaLwRij female mice via tail vein.Recombinant mouse Cst6 (rmCst6) protein and ZA were administered on day 5 post tumor inoculation.At day 25, mice were sacrificed and samples were collected.(B) Reconstructed µCT images of tibia sagittal sections taken from MM mice show trabecular architecture change in bone for MM mice and MM mice treated with either rmCst6 or ZA.(C) Bar-plots present the number of bone lytic lesions on the right medial tibia surface and the trabecular bone parameters: BV/TV, Tb.N, Tb.Th, Tb.Sp, and BMD (N=5).(D) Best representative images of TRAPase staining shows TRAP positive cells (indicated with arrows) in tibiae derived from C57BL/KaLwRij mice injected with 5TGM1 MM cells with or without rmCst6 and ZA treatment.(E) Bar-plots represent histomorphometric analyses of TRAPase-stained number of osteoclast per bone perimeter (N.Oc/B.Pm) and osteoclast surface per bone surface (Oc.S/BS) (N=5).(F) Representative flow cytometry plots presented the GFP + 5TGM1 cells in BM of control and rmCst6 and ZA treated mice.(G) Bar-plots represent the percentage of GFP+ 5TGM1 cells in BM (N=5).(H) Tumor burden was assessed by measuring serum levels of IgG2b (mg/ml) by ELISA (N=5).(I) Bar-plots show serum levels of the bone resorption marker CTX-1 detected by ELISA (N=5).For all measurements data is represented as mean ± SD and was analyzed by oneway ANOVA with Tukey's multiple-comparisons. Denoted markings are considered significant *P < 0.05, **P < 0.01, ***P < 0.001.

Figure 2 :
Figure 2: rmCst6 protein and ZA suppress bone resorption in sex steroid deficient ovariectomized (OVX) mouse model.(A) A schematic model for the OVX mice study.Sixmonth female mice OVX were administered rmCst6 protein or ZA for 6 weeks by i.v or s.c injection twice a week.After 6 weeks, mice were sacrificed and samples were collected.(B) Reconstructed µCT images of tibia from OVX mice show trabecular architecture change in bone following sex steroid depletion and treatment with either rmCst6 or ZA.(C) Bar-plots show the trabecular bone parameters: BV/TV, Tb N, Tb Th, Tb Sp, and BMD (N=7).(D) Best representative images of TRAPase staining shows TRAPase positive cells (indicated with arrows) and adipocyte-like cells present in tibia from Sham, OVX, OVX + ZA and OVX + rmCst6 mice.(E) Bar-plots represents the histomorphometric analyses of the number of

Figure 3 :
Figure 3: scRNA-seq reveals BM microenvironment alterations after rmCst6 and ZA treatment.(A) Experimental workflow for scRNA-seq on BM mononuclear cells.5TGM1-GFP+ MM cells were injected into 8-week-old C57BL/KaLwRij female mice via tail intravenously.Hind limbs were extracted, and BM mononuclear cells from individual mice were sorted out by depleting 5TGM1-GFP+ MM cells.(B) The Uniform Manifold Approximation and Projection (UMAP) plot of BM mononuclear cells derived from control mice (n = 2), MM mice treated with PBS (n = 3), rmCst6 protein (n = 3) or ZA (n = 2).(C) Bar-views show the proportion of various cell types in BM mononuclear cells of control or MM-bearing mice treated with PBS, rmCst6 or ZA protein.

Figure 4 :
Figure 4: The effects of CST6 and ZA on the viability and differentiation of osteoclast, osteoblast and chondrocyte.(A) Mouse osteoclast precursors were differentiated into osteoclasts by addition of M-CSF and RANKL.Different concentrations of rmCst6 protein and ZA were added into the culture media for 4 days.TRAPase staining shows osteoclasts containing multiple nuclei.(B) Bar-plots show the quantification of TRAP + osteoclasts (n=4).(C) MC3T3-E1 cells were differentiated to osteoblasts with 50 μg/ml ascorbic acid and β-glycerophosphate.Alkaline phosphatase staining on day 14 (upper panel) and Alizarin red staining on day 21 (lower panel) showed ZA but not CST6 inhibited osteoblast differentiation and mineralization.(D) MC3T3-E1 cells were incubated with different doses of ZA for 7 days and cell viability was detected by CCK-8 assay (n=4).(E) ATDC5 cells were differentiated to chondrocytes with ITS medium for 14 days, Alcian blue staining (upper panel) showed the glycosaminoglycan (GAG) deposition was suppressed by ZA; ATDC5 cells were induced with 1X ITS medium plus ZA or CST6 for 21days.Endochondral ossification was detected by Alizarin red staining (lower panel).(F) ATDC5 cells were incubated with different doses of ZA for 7 days and cell viability was

Figure 5 :
Figure 5: OVX-induced increase in inflammatory cytokine expression is brought back to sham levels by treatment with either rmCst6 or ZA.Inflammatory cytokine array for Protein samples from L3-L5 mice vertebrae.The 25 targets with the largest initial change in membrane intensity from the sham surgery group to the OVX-PBS group are shown.Ovariectomy causes the upregulation of multiple inflammatory cytokines, which in turn can cause increased bone resorption.Treatment with either rmCst6 or ZA brought levels back to sham.For each group (Sham, OVX-PBS, OVX-CST6, OVX-ZA) L3-L5 vertebrae protein lysate from each sample in the group was pooled.For pooled samples triplicate analysis was performed.

Figure 6 :
Figure 6: RNA-seq analysis of tibia from OVX-mouse model reveals potential different pathways for CST6 and ZA treatment.(A) Experimental workflow for RNA-seq experiment on tibia isolated from OVX mice.(B) 175 unique differentially expressed genes were found for OVX-CST6 treated tibia compared to the 12 differentially expressed genes for OVX-ZA treated tibia.(C) The top 10 genes with largest fold increase in the tibia following rmCst6 treatment of OVX mice are listed in a heat map and sorted by LS Mean.Of interest Slco1a4, an organic anion transporter that may be involved in bringing estrogen into the cell had the largest fold increase (log2 fold = 3.60).(D) The top 10 genes with the greatest fold change (increase or decrease) following ZA treatment of OVX tibia are listed in a heat map sorted by LS Mean.(E-F) Basal tissue and cell level gene expression of Slco1a4 determined by PCR.(G) RNA-seq results for tibia were confirmed using real-time PCR on RNA isolated from L3-L5 vertebrae (n = 7).(H-J) Slco1a4 mRNA levels increases for both pre-osteoclastic Raw 264.7 cells treated with RANKL (20 ng/ml) as well as for pre-osteoblastic ST2 cells (+/-1 mM A2P) (n = 3).Only the increase for Raw 264.7 cells was considered significant.Data represented as mean +/-SD analyzed by one-way ANOVA with Tukey's multiple-comparisons. Denoted markings are considered significant *P < 0.05, **P < 0.01, ***P < 0.001.

Figure 7 :
Figure 7: ERα expression in bone tissue and cells and intracellular total estrogen levels in osteoclast precursors are increased following treatment with rmCST6.(A) ERα protein levels are decreased in L3-L5 vertebrae of OVX mice and expression is brought back to sham surgery levels only by rmCst6 treatment (N = 6).Amido black staining was used as a loading (D) Total intracellular estrogen levels (E1, E2, E3) for cultured osteoclast precursor cells isolated from mouse femur was increased for only 100 ng/ml rmCst6 treatment.Estrogen concentration was normalized to total protein lysate concentration.Graph is the result of triplicate experiments (N = 9).Three ELISA plates were scaled to each other by dividing estrogen concentration (pg/µg) for each sample by average estrogen concentration of control treatment for said plate, then multiplying by the average estrogen concentration for all samples across three plates.(E) Immunohistology staining of Raw 264.7 cells incubated for 24 hr.treated with PBS control solution or 100 ng/ml rmCST6.Cells treated with rmCst6 showed increased expression of ERα, whose transcription known to be upregulated by estrogen, and decreased expression of PPARγ, whose transcription is known to be downregulated by estrogen.(F) Immunohistology staining of Raw 264.7 cells incubated for 7 days treated with PBS control solution or 100 ng/ml rmCst6 (20 ng/ml RANKL).Cells treated with rmCst6 showed increased expression of ERα and decreased expression of PPARγ.(G)Estrogen influx is increased in Raw 264.7 cells by 24 hr.rmCst6 treatment.Raw 264.7 cells were plated onto tissue culture treated slides (Lab-Tek Chamber Slide System, 177437) at a concentration of 2.5 x 10 5 cells/ml and treated with PBS (left panel) or 100 ng/ml rmCst6 (right panel) and 1 pm Estradiol glow (E2 Glow, Jena Bioscience, PR-958S) for 24 hr.After 24 hr.media was aspirated and cells were fixed with cold 4% Paraformaldehyde for 20 min and rinsed with PBS.Following this, cells were covered with DAPI-Fluoromount-G and observed using a Nikon Eclipse T/2 epifluorescent microscope.DAPI ex wavelength: 350 nm, em wavelength: 470 nm; E2 Glow ex wavelength: 467 nm, em wavelength: 618 nm.Data represented as mean +/-SD analyzed by one-way ANOVA with Tukey's multiple-comparisons. Denoted markings are considered significant *P < 0.05, **P < 0.01, ***P < 0.001.