Syeda Masooma Naqvi†
Juan Alberto Panadero Pérez†
Vatsal Kumar
Anneke S. K. Verbruggen
Laoise M. McNamara*- Mechanobiology and Medical Device Research Group, Department of Biomedical Engineering, College of Engineering and Informatics, National University of Ireland Galway, Galway, Ireland
Recent in vitro studies have revealed that the mechanobiological responses of osteoblasts and osteocytes are fundamentally impaired during estrogen deficiency. However, these two-dimensional (2D) cell culture studies do not account for in vivo biophysical cues. Thus, the objectives of this study are to (1) develop a three-dimensional (3D) osteoblast and osteocyte model integrated into a bioreactor and (2) apply this model to investigate whether estrogen deficiency leads to changes in osteoblast to osteocyte transition, mechanosensation, mineralization, and paracrine signaling associated with bone resorption by osteoclasts. MC3T3-E1s were expanded in media supplemented with estrogen (17β-estradiol). These cells were encapsulated in gelatin-mtgase before culture in (1) continued estrogen (E) or (2) no further estrogen supplementation. Constructs were placed in gas permeable and water impermeable cell culture bags and maintained at 5% CO2 and 37°C. These bags were either mechanically stimulated in a custom hydrostatic pressure (HP) bioreactor or maintained under static conditions (control). We report that osteocyte differentiation, characterized by the presence of dendrites and staining for osteocyte marker dentin matrix acidic phosphoprotein 1 (DMP1), was significantly greater under estrogen withdrawal (EW) compared to under continuous estrogen treatment (day 21). Mineralization [bone sialoprotein (BSP), osteopontin (OPN), alkaline phosphatase (ALP), calcium] and gene expression associated with paracrine signaling for osteoclastogenesis [receptor activator of nuclear factor kappa-β ligand (RANKL)/osteoprotegerin OPG ratio] were significantly increased in estrogen deficient and mechanically stimulated cells. Interestingly, BSP and DMP-1 were also increased at day 1 and day 21, respectively, which play a role in regulation of biomineralization. Furthermore, the increase in pro-osteoclastogenic signaling may be explained by altered mechanoresponsiveness of osteoblasts or osteocytes during EW. These findings highlight the impact of estrogen deficiency on bone cell function and provide a novel in vitro model to investigate the mechanisms underpinning changes in bone cells after estrogen deficiency.
Introduction
Osteoporosis is a debilitating bone disease, in which severe bone loss occurs leading to fractures of the hip, wrist, or vertebrae (Balena et al., 1993). The disease is most commonly manifested in postmenopausal women when estrogen production is deficient, and affects approximately 30% of postmenopausal women (Melton et al., 2005). Bisphosphonate drugs target bone resorbing osteoclasts but only reduce osteoporosis fractures by 50% (Randell et al., 2002). However, osteoporosis is not simply a disease of bone loss and recent studies have demonstrated that bone tissue composition is fundamentally altered at the tissue level (Gadeleta et al., 2000; Harris et al., 2000; Zioupos and Aspden, 2000; Efstathiadou et al., 2001; Ciarelli et al., 2003; Mansell and Bailey, 2003; McCreadie et al., 2006; McNamara et al., 2006; Zhu et al., 2008). In particular, changes in bone tissue mineral occur, which are anatomically distinct, do not occur ubiquitously throughout the proximal femur, and may be associated with alterations in the cell mechanical environment (Riggs et al., 2002; McNamara et al., 2005, 2006; Brennan et al., 2011a, b, 2012b, 2014a; Vaughan et al., 2015; Verbruggen et al., 2015, 2016; O’Sullivan et al., 2019; Parle et al., 2020).
Throughout life, maintenance of normal adult bone relies on biophysical stimulation, whereby osteocytes continually receive mechanical cues from daily skeletal loading (Knothe Tate, 2003; Coughlin and Niebur, 2012; Birmingham et al., 2013) and elicit biochemical signaling to osteoblasts and osteoclasts (Klein-Nulend et al., 1996; McAllister, 2000; Nauman et al., 2001; Kapur et al., 2002; Malone et al., 2007; Wang et al., 2008; Birmingham et al., 2012; Li et al., 2012; Lee et al., 2014). Cellular mechanosensors have been identified in osteoblasts and osteocytes, which transduce mechanical signals into the cell through its interactions with the actin cytoskeleton (Ziambaras et al., 1998; Rubin et al., 2006; Malone et al., 2007; Litzenberger et al., 2010; Hoey et al., 2012). Osteocyte cell processes interact with the surrounding extracellular matrix by means of αvβ3 integrins (McNamara et al., 2009; Cabahug-Zuckerman et al., 2018), which connect the matrix to the intracellular cytoskeleton in combination with other proteins in complexes known as focal adhesions (Hughes et al., 1993; Zaidel-Bar et al., 2004). A variety of intracellular signaling cascades (MAP kinase, Rho–ROCK, and Wnt/β-catenin) are activated upon exposure to mechanical stimulation (Rubin et al., 2006; Bonewald and Johnson, 2008; Arnsdorf et al., 2009; Jacobs et al., 2010). Osteocytes produce biochemicals and proteins that activate osteoblasts and osteoclasts to remodel bone (Poole et al., 2005; Robling et al., 2008; Kramer et al., 2010; Lewiecki, 2011; Wijenayaka et al., 2011; Tu et al., 2012; Baron and Kneissel, 2013; Delgado-Calle et al., 2017). In particular, WNT/β-catenin signaling promotes bone formation by osteoblasts through inhibition of sclerostin, a protein produced by osteocytes (Papanicolaou et al., 2009; Wijenayaka et al., 2011; Moustafa et al., 2012; Spatz et al., 2015; Thompson et al., 2015; Hemmatian et al., 2018). Mechanical loading decreases sclerostin production by osteocytes (Holdsworth et al., 2019) enabling bone formation. Osteocytes also regulate osteoclastogenesis and bone resorption through the upregulation or downregulation of receptor activator of NFκB ligand (RANKL) and osteoprotegerin (OPG). RANKL is a ligand for receptor activator of NFκB (RANK) present in the cell membrane of osteoclast precursors. The binding between RANKL and RANK triggers signaling pathways, which promote proliferation and differentiation of mature osteoclasts (Wada et al., 2006; Nakashima et al., 2011), which are responsible for bone resorption. OPG is a soluble glycoprotein that acts as a decoy receptor for RANKL, by binding to RANKL and prevents RANK-RANKL association, thus antagonizing the signaling that promotes osteoclastogenesis (Han et al., 2018).
We previously developed an in vitro two-dimensional (2D) postmenopausal model to study changes in bone biology during estrogen deficiency. We reported that estrogen withdrawal (EW) altered the osteoblast actin cytoskeleton, PGE2 release and expression of RUNX2, COX2, and OPN genes compared to cells that continuously received estrogen (Brennan et al., 2012a, c, 2014b). However, these and other studies of bone mechanobiology and pathophysiology have largely relied on 2D in vitro cell culture and applied mechanical stimulation in the form of flow in fluid shear stress (Alford et al., 2003; Litzenberger et al., 2010; Li et al., 2012; Deepak et al., 2017). Yet such 2D approaches do not represent the complexity of in vivo biophysical cues, wherein osteocytes are embedded within a three-dimensional (3D) matrix and are simultaneously influenced by mechanical cues arising due to daily physical activity. Indeed, a study of osteocytes cultured in collagen-coated 3D scaffolds revealed a significant increase in expression of important genes (SOST, RANKL) compared to the same cells cultured in 2D tissue culture plastic (Spatz et al., 2015). Moreover, osteocyte-selective gene expression (E11/GP38) by differentiated human osteocytes cultured in 3D (with biphasic calcium phosphate particles) was comparable to that of mature human cortical osteocytes, whereas gene expression by human osteoblasts cultured on 2D tissue culture plastic was limited (Boukhechba et al., 2009). Specific matrix properties have been identified that govern osteocyte differentiation and mineralization by osteoblasts in vitro (Mullen et al., 2013; Mc Garrigle et al., 2016). In particular, it has been shown that matrix stiffness (0.58 kPa) controls the phenotypic shift from osteoblasts to osteocytes in a 3D environment in terms of osteocyte dendrite formation. Such strategies provide a means of investigating osteocyte biology in a 3D environment.
To more faithfully represent the complexity of bone mechanobiology in vivo, a 3D culture and bioreactor approach is required. In vitro bioreactors provide a physical growth environment for cells and tissues by subjecting them to various types of mechanical forces, such as hydrostatic pressure (HP) (Klein-Nulend et al., 1986; Pörtner et al., 2005; Chen and Hu, 2006; Henstock et al., 2013; Reinwald et al., 2015). We have previously developed a HP bioreactor that can enhance the in vitro mineralization potential of human mesenchymal stem cells (Freeman et al., 2017). HP has further been shown to enhance mineralization and bone formation when applied to osteoblasts, osteocytes, cell−seeded constructs and ex vivo cultures of chick femurs (Roelofsen et al., 1995; Liu et al., 2010; Henstock et al., 2013; Reinwald et al., 2015). However, no such 3D bioreactor approach has been applied to investigate the interaction between estrogen deficiency, osteoblast-osteocyte differentiation, and mechanical loading.
In this study, we tested the hypothesis that estrogen deficiency alters mineralization and pro-osteoclastogenic potential of osteocytes in 3D constructs under cyclic loading. Specifically, we (1) developed a simplified 3D osteocyte model by combining matrix-based and bioreactor strategies and (2) applied this model to investigate whether estrogen deficiency leads to changes in osteoblast to osteocyte transition, mechanosensation, mineralization, and paracrine biochemical signaling associated with bone resorption by osteoclasts.
Materials and Methods
Osteoblast Cell Expansion and Estrogen Pre-treatment
To closely mimic the in vivo pre-menopausal physiological environment, and to accustom the cells to the estrogen environment prior to our estrogen deficiency experiments, we first expanded MC3T3-E1 osteoblastic cells (ATCC, United States) in standard αMEM supplemented with 10% fetal bovine serum (FBS), 100 U/mL penicillin, 100 μg/mL streptomycin, 2 mM L-glutamine, and 1 × 10–8 M 17β-estradiol (a naturally occurring estrogen derived from cholesterol) and maintained the cells at 37°C in a humidified 5% CO2 environment. This estrogen dosage is within the normal physiological circulating serum levels in mice (Porter et al., 2002; Wood et al., 2007; Haisenleder et al., 2011). Cell culture media was replenished every 3–4 days. The pre-treatment duration (14 days) was determined following preliminary experiments and was deemed to be an appropriate duration to allow cells to become accustomed to estrogen before subsequent withdrawal (Brennan et al., 2014b). All reagents used in this study were from Sigma–Aldrich unless otherwise stated.
Gelatin-Mtgase Hydrogel Encapsulation
The estrogen pre-treated MC3T3-E1 cells were suspended in sterile microbial transglutaminase (mtgase) solution containing Activa WM (Aginomoto foods) at 4 × 106 cells/mL. This mtgase/cell suspension was then enzymatically crosslinked with an equal volume of 6% gelatin (type A, 175 Bloom) (ratio of 1:1), to yield a final mtgase concentration of 0.3% w/w, a cell density of 2 × 106 cells/mL and a gelatin concentration of 3%. This concentration for gelatin-mtgase was established based on previously published work from our lab that demonstrated a matrix stiffness of 0.58 kPa controls the phenotypic shift from osteoblasts to osteocytes in a 3D environment in terms of osteocyte dendrite formation (Mc Garrigle et al., 2016). To produce rectangular constructs (3 mm × 4 mm × 13 mm), the gelatin-mtgase cell suspension was pipetted into custom made polydimethylsiloxane (PDMS) wells and allowed to cool at 4°C for 6 min.
Estrogen Deficiency
These estrogen pre-treated MC3T3-E1 cells encapsulated in gelatin were then cultured under the following experimental conditions: (1) standard growth medium with estrogen supplementation (E) (1 × 10–8 M 17β-estradiol) or (2) standard growth medium with no estrogen supplementation, to mimic the onset of estrogen deficiency in postmenopausal osteoporosis, referred to as EW. For the estrogen supplemented group, estrogen was supplemented every 3 days for the duration of the HP and static experiments described below.
Hydrostatic Pressure (HP) Experiments
After 24 h of incubation under the different estrogen conditions [continued estrogen supplementation (E), EW], constructs from each experimental group were transferred from well plates to sterile, heat-sealed, gas permeable, and water impermeable, EVO cell culture bags (Quest Biomedical, West Midlands, United Kingdom) containing 30 mL of their respective media, as described above. Ten constructs were placed within each EVO cell culture bag and air was removed. For mechanical stimulation, the bags were placed in a custom developed HP bioreactor and were subjected to an intermittent HP regime of 270 kPa and a frequency of 1 Hz for 1 h per day, 5 days per week for the duration of the study (21 days). A pressure of 270 kPa replicates the physiological pressures experienced by osteocytes in the lacunar–canalicular network of load-bearing bones (Zhang et al., 1998; Reinwald et al., 2015; Reinwald and El Haj, 2018). This pressure was found to promote bone growth and mineralization in a developmental model of organotypically cultured ex vivo chick fetal femurs (Henstock et al., 2013). Static control groups were also investigated for all experimental conditions (E, EW), and were similarly contained in heat-sealed EVO cell culture bags and maintained in an open water bath (37°C) for the same duration as the mechanical stimulation periods. After each loading bout, mechanically stimulated and static bags were transferred to a humidified incubator and maintained at 5% CO2 and 37°C. Constructs were assessed at days 1, 7, and 21 by immunofluorescent staining, micro-CT, biochemical analyses, and PCR, as are described further below.
DNA Assay
In order to assess the effects of EW on cell number, DNA content was measured after 1, 7, and 21 day(s) of culture. At each time point, constructs were washed twice with PBS, frozen, and stored at −80°C. Constructs were digested with 3.88 U/mL papain in 0.1 M sodium acetate, 5 mM L-cysteine-HCl, 0.05 M EDTA, pH 6.0 at 60°C under constant rotation for 18 h. Cell number was evaluated using the Hoechst 33258 DNA assay as previously described (Haugh et al., 2011).
Immunofluorescent Staining
Constructs from each experimental group were stained for Actin, DAPI, dentin matrix acidic phosphoprotein 1 (DMP1), αvβ3, and vinculin. Constructs were fixed using 4% paraformaldehyde after 1, 7, and 21 day(s) of culture at 4°C overnight. Cells within the constructs were permeabilized with 0.5% Triton-X in PBS for 10 min at 4°C under agitation. After three washing steps in 1% BSA, the constructs were immersed in 1% BSA blocking solution for 1 h under agitation.
For DMP1 staining, constructs were incubated in mouse monoclonal anti-DMP1 antibody (1:100) at 4°C overnight (Clone 8G10.3, Merck) and then treated with a DylightTM 549 conjugate goat anti-mouse secondary antibody (1:100) (Jackson Immunoresearch) for 1 h under agitation at room temperature. Constructs were further counterstained with phalloidin-FITC at 1.25 μg/mL (1:400) to stain the actin cytoskeleton and DAPI dilactate (1:2000) to stain the nucleus. Negative controls were previously performed on identical MC3T3-E1 encapsulated gelatin constructs, by omitting the primary antibody incubation step and it was found that no DMP1 staining was observed (Mc Garrigle et al., 2016). In this study, Day 1 cell encapsulated gelatin constructs acted as negative controls, since there is no expected DMP1 expression at this early time point and negative staining for DMP1 would confirm that there is no aspecific staining due to the gelatin constructs or in cells prior to differentiation.
We performed αvβ3 and vinculin staining to investigate whether the integrin αvβ3, previously shown to play an important role in osteocyte mechanotransduction under EW in 2D (Geoghegan et al., 2019), is also implicated under EW in 3D. For integrin αvβ3 and associated membrane protein vinculin staining, constructs were incubated in mouse anti-vinculin primary antibody (1:100) at 4°C overnight and then incubated in a DylightTM 549 conjugate goat anti-mouse secondary antibody (1:100) (Jackson Immunoresearch) for 1 h under agitation at room temperature. Constructs were incubated in an anti-αvβ3 antibody directly conjugated to Alexa Fluor® 488 (1:100) (Santa Cruz) overnight at 4°C and then incubated in DAPI dilactate (1:200) to stain the nucleus.
To quantify all immunofluorescent staining, Z-stack imaging was carried out using a Fluoview FV1000 confocal laser scanning microscope system (Olympus) at a magnification of either 20x (air) or 60x (oil immersion) with a step size of 5 or 1 μm, respectively. All stacks were obtained at the same intensity setting between groups. NIH ImageJ software and FV10-ASW 2.0 Viewer software were used to analyze the images. The images used for analysis were created by combining the stacks of images at maximum intensity projections.
Quantification of Dendritic Cells and Intensity of DMP1 and Integrin αvβ3
Confocal images of actin/DAPI stained constructs were used for quantification of dendritic cells. Cell processes were defined as cellular features composed of actin, located at the cell membrane, which extended for a distance of at least 10 μm from the cell body, as previously described (Mullen et al., 2014). Using this method, cell morphologies were classified as follows: (1) “not connected dendritic” cells exhibited a small cell body and long thin cell processes of at least 10 μm, without forming interconnections with neighboring cells within the constructs, (2) “not dendritic” cells exhibited a spherical or cuboidal morphology and no cell processes, and (3) “interconnected dendritic” cells exhibited a small cell body and long thin cell processes of at least 10 μm, forming interconnections with neighboring cells within the constructs. The percentage of each type of cell classification was quantified using NIH ImageJ software particle analysis, combining the stacks of images at maximum intensity projections. Osteocytes are defined as interconnected dendritic cells.
Cell area and overall fluorescence intensity were measured using the images stained. Cell area was measured by thresholding the images to remove background fluorescence and then using the “wand tool” to select the region of interest around each cell. When the software could not detect the cell perimeter, the region of interest was drawn manually with the “freehand tool.” The cell area thresholds were then used to determine the intensity of the staining. All fluorescent intensities were measured using the integrated density of each cell with their corresponding background integrated density subtracted.
Micro-Computed Tomography (μCT)
Like constructs for immunofluorescence staining, constructs for micro-computed tomography (μCT) were fixed using 4% paraformaldehyde at 4°C overnight. A construct from each static and mechanically stimulated group was analyzed for mineralization using μCT scanning methods (Scanco Medical μCT100) with the following parameters: 9 mm holder, 45 kVp tube voltage, 200 μA intensity, Al filter, 5 μm voxel size, 300 ms integration time, and a frame averaging value of 2. Samples were scanned in air and paper strips were used to minimize movement within the chamber. A density threshold of 370 mg HA/cm3 (3964 Hounsfield units) was established to exclude hydrogel material and where only mineralized regions remained. This threshold level is higher than the standard range for isolating mineral (90–120 mg HA/cm3) because hydrogel and collagen materials have higher density ranges (Vetsch et al., 2015, 2017; Manju et al., 2018; Raina et al., 2018) (Supplementary Figure S1). μCT scans required a Gaussian filter (0.8 sigma, 1.0 voxel support) to reduce noise. A volume of interest (VOI) of 3.95 mm3 was 3D reconstructed using Scanco Medical software (Switzerland) and evaluated for mineralization properties such as bone volume fraction (BV/TV), bone mineral density distribution (BMDD), and full width at half max (FWHM) to assess the quantity and distribution of mineralized material within the constructs (Vetsch et al., 2015; Melke et al., 2016).
Extracellular ALP Activity
Cell culture medium was sampled for each experimental group and subsequently analyzed for ALP secreted to the media. ALP activity was determined using a colorimetric assay of enzyme activity (SIGMAFAST p-NPP Kit), which uses p-nitrophenyl phosphate (pNPP) (nmol) as a phosphatase substrate, with ALP enzyme as a standard. After 1, 7, and 21 day(s) EVO cell culture bags containing constructs were maintained in an incubator for 2 h after the last bout of mechanical stimulation before media was removed, frozen and stored at −80°C. After thawing, 40 μL of the medium was added to a 96-well plate in triplicate with a 50 μL of pNPP solution. The plate was shielded from direct light at room temperature for one hour. The plate was read at 405 nm in a Synergy HT Multi-mode microplate reader. Relative colorimetric readings were converted to ALP activity by creating a standard ALP curve. ALP production was normalized to DNA content to establish whether changes in mineralization were related to (1) the overall mineral capacity of individual cells or (2) changes in the cell population available to produce that mineral.
Mineralization
Calcium content within constructs was determined using a Calcium Liquicolour kit (Stanbio Laboratories, Syntec, Ireland) according to the manufacturer’s protocol. After 1, 7, and 21 day(s) of incubation, constructs were washed twice with PBS, frozen, and stored at −80°C. Constructs were then thawed and digested by adding 1 mL of 1 M hydrochloric acid (HCL) to each construct and maintaining the constructs at 60°C under constant rotation for 18 h. 10 μL each of the digested constructs and assay standards was added to a 96-well plate and 200 μL of the working solution was added. The plate was read on a synergy HT Multi-mode microplate reader at an absorbance of 550 nm, as previously described (Freeman et al., 2015).
Quantitative Real-Time PCR
Total RNA was isolated with a custom CTAB method (Koster et al., 2016). Briefly, each construct was digested with 250 μL of RNAse-free proteinase K solution 5 mg/mL, containing 5 mM CaCl2. The digested constructs were mixed with 500 μL of CTAB solution (2% [w/v] Cetyl trimethylammonium bromide, 1.4 M NaCl, 20 mM EDTA, and 100 mM Tris, pH 8) containing 1% (v/v) β-mercaptoethanol. Then 500 μL of chloroform were added to each construct and centrifuged at 14,000 g for 2 min at room temperature. The upper/aqueous phase was transferred to a fresh tube and 800 μL of isopropanol (Fisher) were added to precipitate total RNA by centrifugation (14,000 g, 15 min, room temperature). The pellet was washed with 600 μL 70% (v/v) ethanol and dissolved in 20 μL low TE buffer (Thermo Fisher Scientific) for 15 min at 65°C. RNA purity and yield were assessed using a spectrophotometer (DS-11 FX, DeNovix), with 260/280 and 260/230 ratios over 1.8 for all constructs. 1 μg of RNA from each construct was transcribed into cDNA using Qiagen Quantinova reverse transcription kits and thermal cycler (5PRIMEG/O2, Prime). Relative gene expression was studied by quantitative real-time polymerase chain reaction (qRT-PCR). The genes of interest included ALP (alpl- alkaline phosphatase) for mineralization, RANKL (TNFSF11—Receptor Activator for Nuclear Factor κ B Ligand) and OPG (TNFRSF11B—osteoprotegerin) for pro-osteoclastogenic response, bone sialoprotein (BSP), OCN (bglap—osteocalcin), and OPN (spp1—osteopontin) for osteoblast phenotype and DMP1 a marker of differentiation to osteocyte phenotype. Rpl13A was used as reference gene (Table 1). qRT-PCR was carried out with 12.5 ng of cDNA per construct, using Qiagen Quantinova SYBR Green PCR kit and a StepOne Plus PCR machine (Thermo Fisher Scientific). The PCR was conducted with an enzyme activation step at 95°C for 10 min, and then for 40 cycles with the following steps: 10 s extension at 95°C, specific annealing temperatures for each primer pair (Table 1) for 10 s, and extension 72°C for 20 s. Analysis of the results was done using the Pfaffl method (Pfaffl, 2001), and the results are expressed as relative quantitative changes to constructs from day 1 of static culture under E.
Table 1. List of primers employed for qRT-PCR, including sequences (5′–3′) and annealing temperatures (Tm).
Statistical Analysis
Statistical analyses were performed using GraphPad Prism (version 5) software. Two−way ANOVA was used for analysis of variance with Bonferroni’s post hoc tests to compare between groups. The results are displayed as mean ± standard deviation. Significance was accepted at a level of p ≤ 0.05. The entire experiment was repeated three times with three to four replicates analyzed for each experimental group per repeat.
Results
3D Osteocyte Model Under Continued Estrogen Treatment
We first sought to establish whether this model may be considered an appropriate 3D osteocyte model under normal estrogen conditions (E) and compared static and mechanically stimulated constructs. Actin staining revealed that the MC3T3-E1 cells began to form dendritic cell processes by day 7 under both static and mechanically stimulated conditions (Figures 1A,B). By day 21, interconnected dendritic cells were present (Figures 1C–F). Simultaneously, DMP1 staining was positive at day 21 (Figures 1C,D,G,H) indicative of phenotypic differentiation toward osteocytes. Interestingly, the intensity of DMP1 staining decreased significantly with mechanical stimulation (p < 0.05) (Figure 1I). Day 1 MC3T3 encapsulated gelatin constructs stained negatively for the DMP1 marker, which confirms that there was no aspecific staining due to the gelatin constructs or in cells prior to differentiation (Supplementary Figure S2). Next, we performed immunostaining for integrin αvβ3 on day 21 constructs. All constructs stained positively for integrin αvβ3 and vinculin (Figures 1J–O). Semi-quantitative analysis revealed a significant decrease in αvβ3 intensity under mechanical stimulation compared to static culture (p < 0.005, Figure 1P).
Figure 1. Constructs under continued estrogen treatment under static conditions and mechanical stimulation. Actin (Green) staining at day 7 (A,B). Red arrows indicate dendritic cells. Actin (green)/DMP1 (red) staining at day 21 (C–H). Merged channels for Actin/DMP1 staining (C,D). Green channel for Actin staining only (E,F). Red channel for DMP1 staining only (G,H). DMP1 intensity (I). αvβ3 (green)/vinculin (red) staining at day 21 (J–O). Merged channels for αvβ3/vinculin staining (J,K). Green channel for αvβ3 staining only (L,M). Red channel for vinculin staining only (N,O). αvβ3 intensity (P). Images were taken from the construct surface to a depth of approximately 65 μm. Scale bar = 50 μm. *p < 0.05, **p < 0.01.
The time of cell culture had a significant effect on gene expression under all the conditions. For example, the expression of DMP1 (Figure 2A) was upregulated at day 21 (p-value < 0.0005) compared to earlier time points of days 7 and 1. Our DMP1 gene expression data correlates with the DMP1 staining, which revealed positive staining at day 21 only. Interestingly, the difference in DMP1 intensity observed between static and mechanically stimulated conditions was not apparent in gene expression data. Expressions of OCN (Figure 2B) and ALP (Figure 2C) peaked at day 7 while by day 21, they were downregulated.
Figure 2. Constructs under continued estrogen treatment under static conditions and mechanical stimulation. Gene expression of DMP1 (A), OCN (B), and ALP (C). All constructs normalized to static constructs at day 1, &significance compared to earlier/later time point(s) for the same experimental group; p < 0.005.
To confirm that we had mineralization as well as osteocyte differentiation, we performed μCT on our 3D models. Analysis of constructs at day 21 confirmed the presence of clustered regions of density values at and above levels assigned to bone mineral (Vetsch et al., 2015). Figures 3A,B present BMDD analysis of these constructs with corresponding results in Table 2. BV/TV reveal the mechanically stimulated E construct contains ∼4.61% bone-like mineral. The percentage volumes of high mineral density (above the 75th percentile marker) are approximately 16% in the static construct but lower in the mechanically stimulated group (∼13%). Consequentially, the low mineral density (below the 25th percentile marker) is lesser in the static (∼6%), compared to the mechanically stimulated (∼15%) constructs. However, the volume of most frequent mineral value (Mmode) was highest for the mechanically stimulated E sample (597.55 mg HA/ccm), and lowest for its static counterpart (533.6 mg HA/ccm), indicating the highest density recorded occupies a larger volume of the 3D construct following mechanical stimulation. FWHM, an indicator of mineral heterogeneity, is greater in the mechanically stimulated constructs than the static constructs. 3D reconstructions of the thresholded regions (Figures 3C–F) demonstrate large spaces between mineral clusters. Figures 3G,H reflect this observation, where views of the constructs in grayscale cross-sections (distributions at 0.3 mm intervals) display bone-like mineral along the external layer of the construct, and some highly clustered regions are formed within.
Figure 3. Constructs under continued estrogen treatment under static conditions and mechanical stimulation. Graphs representing BMDD analysis, featuring full width half max bars (FWHM) to quantify heterogeneity, the most frequent mineral density (Mmode), and markers indicating mineral density at the 25th (Mlow) and 75th (Mhigh) percentiles (A,B). 3D reconstructions of the VOI (C–F). Orthogonal view (C,E). Top view (D,F). Scale bars = 100 μm. Grayscale views of cross-section of constructs, alongside four sub-images of example slices, 0.3 mm apart, throughout the constructs (G,H). All scale bars = 1.0 mm.
Table 2. BMDD results for static and mechanically stimulated estrogen and estrogen withdrawn samples, including BV/TV and volume percentages of mineral in the lower, medium, and upper quadrants of distributed density.
3D Osteocyte Estrogen Withdrawal Model
We next investigated osteocyte differentiation under conditions that mimic the onset of estrogen deficiency in postmenopausal osteoporosis. Cells under EW began to form processes early in the culture period (day 7, Figures 4A,B) and by day 21 most of the dendritic cells were interconnected (Figures 4C–F). Furthermore, DMP1 staining revealed increased DMP1 under EW when compared to E (Figures 4C,D,G,H). Interestingly, after quantification of the mean intensity per cell, the DMP1 intensity increased significantly with mechanical stimulation (p < 0.05) compared to static controls, which was contrary to that observed under estrogen conditions (Figure 4I). Furthermore, there was a significant increase in DMP1 intensity under EW (93,458 ± 82,594) compared to E (28,778 ± 21,376) for mechanically stimulated constructs only (p < 0.0001). Next, we performed immunostaining of integrin αvβ3 on day 21 constructs. EW constructs stained positively for integrin αvβ3 and vinculin at day 21 (Figures 4J–O), and semi-quantitative analysis revealed that there was a significant decrease in mechanically stimulated groups compared to static groups (p < 0.001, Figure 4P). This is similar to results obtained for the 3D osteocyte model under continued estrogen treatment.
Figure 4. Constructs under estrogen withdrawal (EW) under static conditions and mechanical stimulation. Actin (green) staining at day 7 (A,B). Red arrows indicate dendritic cells. Actin (green)/DMP1 (red) staining at day 21 (C–H). Merged channels for actin/DMP1 staining (C,D). Green channel for actin staining only (E,F). Red channel for DMP1 staining only (G,H). DMP1 intensity (I). αvβ3 (green)/vinculin (red) staining at day 21 (J–O). Merged channels for αvβ3/vinculin staining (J,K). Green channel for αvβ3 staining only (L,M). Red channel for vinculin staining only (N,O). αvβ3 intensity (P). Images were taken from the construct surface to a depth of approximately 65 μm. Scale bar = 50 μm. *p < 0.05, ***p < 0.005.
DMP1 gene expression (Figure 5A) was upregulated at day 21 (p-value < 0.0005). Our DMP1 gene expression data correlate with the DMP1 staining, which revealed positive staining at day 21 only. Interestingly, the difference in DMP1 intensity observed between estrogen conditions was not apparent in gene expression data. Expressions of OCN (Figure 5B) and ALP (Figure 5C) peaked at day 7 while at day 21, they were downregulated.
Figure 5. Constructs under estrogen withdrawal (EW) under static conditions and mechanical stimulation. Gene expression of DMP1 (A), OCN (B), and ALP (C). All constructs normalized to static constructs at day 1, &significance compared to earlier/later time point(s) for the same experimental group; p < 0.005.
Similar to the osteocyte model under continued estrogen treatment, μCT analysis of estrogen deficient constructs at day 21 confirmed the presence of clustered regions of density values at and above levels assigned to bone mineral (Vetsch et al., 2015). Figures 6A,B present BMDD analysis of these four constructs with corresponding results in Table 2. Like the osteocyte model under continued estrogen treatment, the percentage volume of high mineral density is approximately 16% in the static constructs but notably lower in the mechanically stimulated group. Consequentially, the low mineral density range (below the 25th percentile marker) is greater in the static, compared to the mechanically stimulated constructs. In contrast to the osteocyte model under continued estrogen treatment, the volume of most frequent mineral value (Mmode) was highest for the static EW sample (577.93 mg HA/ccm), and lowest for its mechanically stimulated counterpart (560.43 mg HA/ccm), indicating the highest density recorded occupies a larger space within the 3D construct under static conditions. Like the osteocyte model under continued estrogen treatment, FWHM is greater in the mechanically stimulated constructs than the static constructs. 3D reconstructions of the thresholded regions (Figures 6C–F) demonstrate that the mechanically stimulated constructs formed clustered bone-like mineral regions throughout the cross-sections. Figures 6G,H display these distributions in grayscale cross-sections at 0.3 mm intervals.
Figure 6. Constructs under estrogen withdrawal (EW) under static conditions and mechanical stimulation. Graphs representing BMDD analysis, featuring full width half max bars (FWHM) to quantify heterogeneity, the most frequent mineral density (Mmode), and markers indicating mineral density at the 25th (Mlow) and 75th (Mhigh) percentiles (A,B). 3D reconstructions of the VOI (C–F). Orthogonal view (C,E). Top view (D,F). Scale bars = 100 μm. Grayscale views of cross-section of constructs, alongside four sub-images of example slices, 0.3 mm apart, throughout the constructs (G,H). All scale bars = 1.0 mm.
Estrogen Withdrawal Increased Cell Interconnectivity
There was significantly greater DNA content for constructs under E conditions, compared to EW, for both static and mechanically stimulated conditions at day 1. The difference between E and EW was also evident for mechanically stimulated conditions only at day 7 (Figure 7A). DNA content at day 21 was significantly lower than days 1 and 7. The percentage of cells with interconnected processes was significantly higher for the constructs maintained under EW (77.0 ± 21.8%) compared to those maintained under E (29.0 ± 6.1%) for static conditions only (Figure 7B). Under mechanical stimulation, this same trend was observed between EW and E (76.7 ± 21.9 and 50.7 ± 18.6%, respectively), although the difference was not significant. Figures 7C–J illustrate representative actin/DAPI images at day 21, which revealed that interconnected dendritic cells were present under both static and mechanically stimulated conditions. These images were included in the quantification of cell interconnectivity.
Figure 7. Constructs under both continued estrogen treatment (E) and estrogen withdrawal (EW) under static conditions and mechanical stimulation. Total DNA (ng/mg) content normalized to wet weight at days 1, 7, and 21 (A). Percentage cells after 21 days of culture (B). Actin (green) and Dapi (blue) staining after 21 days of culture illustrating cell morphology (C–F); scale bar = 50 μm. Magnifications of actin/DAPI stained images illustrating interconnections after 21 days of culture (G–J); scale bar = 10 μm. Images were taken from the construct surface to a depth of approximately 65 μm; ***significance compared to E for the same experimental group and time point, &significance compared to earlier time point(s) for the same experimental group, p < 0.001.
Estrogen Withdrawal Increased Osteoblast Mineralization
There was a significant increase in ALP activity under EW compared to E for all time points under mechanical stimulation and for days 1 and 7 under static conditions (Figure 8A). Similarly, there was a significant increase in calcium production for constructs maintained in EW under mechanical stimulation compared to E for all time points (Figure 8B). There were no differences in mineral deposition between static and mechanically stimulated groups for E. However, for EW groups, there was significantly more ALP activity and calcium content under mechanical stimulation at day 7 compared to those maintained under static conditions. Of note, both ALP activity and calcium content increased significantly from days 1 and 7 to day 21. Immunofluorescence of collagen type I illustrated pericellular deposition at day 21 (Supplementary Figure S3) for all groups.
Figure 8. Constructs under both continued estrogen treatment and estrogen withdrawal (EW) under static conditions and mechanical stimulation. ALP activity normalized to DNA at days 1, 7, and 21 (A). Calcium normalized to DNA at days 1, 7, and 21 (B); #significance compared to static for the same experimental group and time point, &significance compared to earlier time point(s) for the same experimental group, ***significance compared to E for the same experimental group and time point; p < 0.001.
Estrogen Withdrawal Upregulates Early Expression of Pro-mineralization Markers and Late Expression of RANKL
Expression of BSP (Figure 9A) was upregulated only in EW under mechanical stimulation at day 1 (p < 0.001). By day 21, its expression significantly increased, relative to earlier time points, for all groups regardless of culture condition. Like BSP at day 1, expression of OPN (Figure 9B) was upregulated in EW under mechanical stimulation (p < 0.001). However, at days 7 and 21, its expression decreased to levels comparable to constructs under E.
Figure 9. Constructs under both continued estrogen treatment and estrogen withdrawal (EW) under static conditions and mechanical stimulation. Gene expression of BSP (A), OPN (B), RANKL (C), and OPG (D). Ratio of expression between RANKL and OPG (E). All constructs normalized to constructs of day 1 and E, #significance compared to static for the same experimental group and time point, &significance compared to earlier time point(s) for the same experimental group, p < 0.05 ***significance compared to E for the same experimental group and time point; p < 0.001, *significance compared to E for the same experimental group and time point; p < 0.05.
After 21 days of culture, gene expression of RANKL was significantly increased under EW for both static and mechanically stimulated conditions (Figure 9C). OPG expression on the other hand was not significantly different among all groups (Figure 9D). Notably, there was a significantly higher RANKL/OPG ratio (p < 0.005) for constructs under EW at day 21 (Figure 9E) for both static and mechanically stimulated conditions (2.04 ± 0.27 and 2.47 ± 0.39 fold increase, respectively).
Discussion
In this paper, we developed a 3D osteocyte model combining matrix-based and bioreactor strategies to more faithfully represent in vivo biomechanical cues. We demonstrated that the strategy successfully enabled osteoblast to osteocyte differentiation, as evidenced by interconnected DMP1 positive dendritic cells that also stained positive for αvβ3 and vinculin. Importantly, we report here for the first time that withdrawal of estrogen altered the osteocyte response compared to continued estrogen supplementation. Specifically osteocytes cultured under EW conditions demonstrated increased DMP1 intensity, cell interconnectivity, mineralization, and paracrine signaling for osteoclastogenesis (RANKL/OPG). Here, the interaction between mechanical stimulation and EW is of importance, whereby significant changes due to EW were only apparent in the presence of HP, in terms of increases in DMP1 intensity (day 21), ALP activity (days 7 and 21), calcium production (days 7 and 21), expression of markers BSP and OPN (day 1), and ratio of RANKL:OPG (days 7 and 21) under EW.
It must be noted that there are several limitations associated with this study. First, the MC3T3-E1 cell line may not represent the behavior of primary osteoblasts. However, these cells are an accepted, commercially available model of primary osteoblast function, because they express ALP, produce mineral, and can differentiate into osteocyte-like cells in vitro (Sudo et al., 1983). Moreover, the use of cell lines eliminates donor variability associated with primary cell culture and thus MC3T3-E1 cells have been studied extensively (Luppen et al., 2003; Balkan et al., 2005; Nakano et al., 2007; Kyono et al., 2012; Wang et al., 2012; Hoac et al., 2013). Second, the withdrawal of estrogen in the current study is abrupt and thus unlike in vivo estrogen depletion during menopause in females. However, ovariectomized animal models also undergo an immediate reduction in circulating estrogen levels and these have been shown to be a good representation of osteoporosis (Brennan et al., 2009; Kennedy et al., 2009). Another limitation was that there was only one sample from each group for micro-CT and so our conclusions regarding mineralization are based on the biochemical analysis. Nonetheless, micro-CT data provided a preliminary spatial understanding of mineral deposition, which indicated an increase in low-density mineral and the most frequent mineral density under EW and mechanical stimulation. Our determination of osteocyte-like cells was based on DMP1 positive staining and the presence of dendritic cell processes, following our previous approach (Mullen et al., 2013; Mc Garrigle et al., 2016). However, it should be noted that DMP1 is expressed by both late osteoblasts and osteocytes (Feng et al., 2002) and so these cells may not be mature osteocytes. Future studies could analyze later markers of osteocyte differentiation, in particular the osteocytic marker Sost/sclerostin. Finally, the cell culture media used for these experiments contains phenol red, which has reported estrogenic activity (Berthois et al., 1986; Welshons et al., 1988). However, the concentration of phenol red in αMEM culture media (∼11 mg/L) used in this study is within the range of DMEM/F12 standard culture media (∼9 mg/L), which was shown to be insufficient to cause estrogenic effects (Moreno-Cuevas and Sirbasku, 2000). Moreover, we used the same batch of media for both experimental groups and further supplemented this batch with estradiol (1 × 10–8 M) for the estrogen groups only and so we do not have a concern that phenol red confounded the results reported here. Indeed, we anticipate that the use of indicator free media would exacerbate the differences between the EW group and the continued estrogen supplemented group.
This study reported changes in osteoblast-osteocyte transition with 3D encapsulation inducing the differentiation of osteoblast-like cells without the addition of osteogenic growth factors. Specifically, under continued estrogen supplementation, cell processes began to form early in the culture period and by end of culture, we observed a largely interconnected network of cells within the 3D construct. This network was characterized by long osteocyte-like interconnected dendrites. Indeed, osteocytes form a complex interconnected dendritic network in vivo allowing for communication with their neighbors and with osteoblasts and osteoclasts on bone surfaces (Bonewald, 2011; Divieti Pajevic and Krause, 2019) and thus for regulation of bone formation and bone resorption. Furthermore, it is established that as an osteoblast differentiates into an osteocyte, expression of the osteoblast marker enzyme ALP is greatly reduced (Mikuni-Takagaki et al., 1995; Nakano et al., 2004), as well as reduced expression of OCN (Bonewald, 2011) and increased expression of DMP1 (Narayanan et al., 2003; Rios et al., 2005). In this study, gene expression of ALP and OCN was upregulated early in the culture period and downregulated by the end of culture for all groups. In contrast, DMP1 was upregulated by the end of culture and immunostaining of DMP1 revealed positive staining also. Of note, the DMP1 staining appeared to be nuclear for both experimental conditions. In addition to its role in the MAP kinase signaling cascade (Kulkarni et al., 2000), DMP1 may function as a transcription factor in the nucleus (Narayanan et al., 2003). One study examined nuclear localization of DMP1 in MC3T3-E1 preosteoblast cells (Siyam et al., 2012) and demonstrated two subpopulations with either nuclear or cytoplasmic localization of DMP1. In addition, we found positive integrin αvβ3 staining by the end of culture. These results suggest the 3D encapsulation is committing the MC3T3 cells to ultimate differentiation into osteocytes. Indeed, it has previously been demonstrated that 3D encapsulation of osteoblast-like cells induces differentiation (Atkins et al., 2009; Boukhechba et al., 2009; Woo et al., 2011; Uchihashi et al., 2013). However, in the current study, we show that osteocyte differentiation and the formation of interconnections are governed by a 3D matrix without the addition of osteogenic growth factors.
Osteoblasts and osteocytes express estrogen receptor α and β (ERα and ERβ), which play a role in regulating both cell survival and mechanosensation (Lee et al., 2003; Saxon et al., 2012; Castillo et al., 2014). In osteoblasts, estrogen has been shown to play a protective role by inhibiting osteoblast apoptosis and increasing its lifespan. Estrogen exerts these effects due to activation of the Src/Shc/ERK signaling pathway and downregulating JNK, which alter the activity of transcription factors such as CREB and c-Jun/cFos (Kousteni et al., 2001, 2003). Supplementation of estrogen in osteoblasts has also been shown to increase Opg expression (Tomkinson et al., 1997, 1998; Emerton et al., 2010; Florencio-Silva et al., 2018) and augment Cox-2 (via β1 integrins and ERs) response to fluid shear stress (Sterck et al., 1998; Bakker et al., 2005, 2006; Voisin and McNamara, 2015; Deepak et al., 2017; Geoghegan et al., 2019), and decrease RANKL and SOST expression (Brennan et al., 2012c, 2014b). Estrogen supplementation has been shown to have a protective role against osteocyte apoptosis (Plotkin et al., 2005; Marathe et al., 2012), increase connexion 43 gap junction expression and mechanosensitivity (Ren et al., 2013), and increase osteogenic signaling by MLO-Y4 osteocytes (Deepak et al., 2017). Conversely, estrogen deficiency alters osteoblast and osteocyte responses (Sterck et al., 1998; Bakker et al., 2005, 2006; Voisin and McNamara, 2015; Deepak et al., 2017; Geoghegan et al., 2019), and an in particular is associated with an increase in osteocyte apoptosis (Brennan et al., 2012c, 2014b). Ovariectomized rodents and osteoporotic humans display significant increases in osteocyte apoptotic markers (Tomkinson et al., 1997, 1998; Emerton et al., 2010; Florencio-Silva et al., 2018).
Importantly, we observed the notable changes in osteoblast-osteocyte differentiation under postmenopausal conditions, stimulated by EW, as evidenced by actin staining, associated quantification of interconnected dendrites and DMP1 staining. Comparing EW to continued estrogen supplementation, we found the number of interconnected cells identified late in the culture period was significantly greater under EW. DMP1 was also identified within the matrix by immunofluorescence late in the culture period, and its intensity was increased under EW and mechanical stimulation, when compared to continued estrogen supplementation. Of note, DMP1 gene expression was upregulated over time with no significant difference between groups. Like continued estrogen supplementation, we found positive integrin αvβ3 staining by the end of culture. The quantification from fluorescence staining of αvβ3 in this study revealed no significant differences between estrogen conditions. Importantly, our results show for the first time the effects of estrogen deficiency in a 3D construct.
It is interesting that early in the culture period, MC3T3 cell number was lower under EW when compared to cells that received continued estrogen supplementation. This reduction in cell number could be explained by cell death or a reduction in proliferation and concomitant differentiation (Ruijtenberg and van den Heuvel, 2016). Notably, we also observed increased expression of osteogenic markers BSP and OPN mRNA under mechanical stimulation at early timepoints. Both proteins increase bone cell adhesion in the mineralized collagen matrix during bone formation (McKee and Nanci, 1996; Giachelli and Steitz, 2000; Malaval et al., 2008). BSP is produced by osteoblasts, while OPN is produced by both osteoblasts and osteoclasts (Denhardt and Noda, 1998; Lenton et al., 2015; Icer and Gezmen-Karadag, 2018). Although BSP acts as an initiator of hydroxyapatite crystal formation in the bone matrix, OPN has an antagonistic effect and is a strong inhibitor of hydroxyapatite crystal growth (Boskey et al., 2002; Harmey et al., 2006). The upregulation of both of these genes under EW may be a reflection of the increase in both mineralization and bone resorption observed during early osteoporosis (McNamara, 2010). Taken together the early reduced cell number and BSP and OPN mRNA expression might indicate that under EW the MC3T3s transitioned toward osteocytes as soon as estrogen was withdrawn.
With regard to mineralization, groups maintained under EW resulted in greater ALP activity early in the culture period and this difference was maintained with time in culture. EW caused no effect during the time frame of the study at the transcription level (ALP gene expression). Of note, there was a difference in gene expression and protein activity at day 21, with reduced ALP gene expression and increased ALP protein activity compared to day 7. It should be noted that gene expression is transient, and precedes protein expression, and as such our gene expression time points (days 1, 7, and 21) may not have captured the upregulation of ALP gene expression, but our results do reflect ultimate changes in ALP proteins that arise in estrogen deficiency. However, cells maintained under EW demonstrated increased calcium content as evidenced by biochemical analyses. This correlates with our previous 2D in vitro study which reported that EW caused significantly increased osteoblast mineralization compared to continued estrogen supplementation (Brennan et al., 2014b). Moreover, higher mineral content has previously been reported in osteoporotic bone compared to healthy bone from human patients (Dickenson et al., 1981). In another study, analogs of gonadotrophin releasing hormone (GnRH) induced short-term estrogen suppression which resulted in bone accumulation with a higher mineralization density (Boyde et al., 1998). In vivo studies have shown similar results; estrogen deficiency caused by ovariectomy increased plasma levels of systemic factors that stimulate proliferation and differentiation of osteoblasts in rats and thus increased bone formation compared sham-operated rats (Yokose et al., 1996). In the same study, osteoblasts from the ovariectomized rats produced more ALP and mineralized bone nodules compared to osteoblasts from the sham-operated rats. It is important to note that the increased calcium content observed under EW in our study suggests that EW itself may be able to enhance the deposition of mineral in the extracellular matrix even in the absence of osteogenic growth factors. The increased expression of BSP and OPN with mechanical loading and EW very early in the culture period, as well as the increased intensity of DMP1 by the end of culture under EW, may provide one explanation for the increase in calcium content. These three phosphoproteins of the SIBLING family have been widely implicated in the regulation of biomineralization due to an abundance of phosphorylated serine and threonine acid-rich domains in their composition (He et al., 2019). These domains can bind to calcium and therefore act as nucleators, inhibitors, anchoring molecules, growth modifiers, or as substrates for mineral deposition. Furthermore, BSP and DMP1 have been shown to promote mineralization (Padovano et al., 2015). Gelatin is the only component of our scaffolds and there is no exogenous source of calcium phosphate, therefore all mineral nucleation must be related to the encapsulated cells. It is possible that mineral formation arises either due to some biological mechanism or nucleation on dead cells/cell debris. As cell number was lower under EW and mechanical stimulation, cell death may be a contributor to the higher mineralization reported in that group. However, these same conditions corresponded to higher gene expression of BSP and OPN at the early time point of day 1, which suggests early activation of biological mineralization processes. Interestingly, preliminary μCT analysis revealed that mineral depositions were thinly formed along the construct periphery and clustered toward the center, a pattern evident of bone matrix in similarly cultured scaffolds with hMSCs in a silk fibroin scaffold (Hagenmüller et al., 2007; Vetsch et al., 2015). Moreover, deposition of mineral was more heterogeneous within mechanically stimulated constructs compared to the static constructs for both estrogen treated and EW conditions, but that withdrawal of estrogen induces greater homogeneity of bone-like mineral compared to continued estrogen supplementation. It should be noted that the micro-CT mineralization results are strictly indicative as no replicas were available for duplicate analysis and therefore could not be statistically analyzed.
Here we report an early upregulation of genes BSP and OPN, which encode proteins important for osteoclast adherence (Ross et al., 1993) and act as substrates for the ECM degrading protein TRAP (Ek-Rylander et al., 1994). We also report a late upregulation of the gene RANKL, which encodes the protein RANKL important for bone resorption by osteoclasts (Buckley and Fraser, 2002). These findings are consistent with previously published work in 2D which report exacerbated osteoblast-induced osteoclastogenesis (Allison and McNamara, 2019) and increased RANKL mRNA expression by MLOY4 osteocytes under EW (Geoghegan et al., 2019). During estrogen deficiency, osteoblasts alter paracrine regulation of osteoclast differentiation (TRAP, CTSK, NFATc1) and resorption (Allison and McNamara, 2019). In osteocytes, EW alters integrin-based mechanosensors in MLYO4 osteocyte-like cells, resulting in defective COX2 responses to oscillatory fluid flow, and an increased RANKL/OPG ratio (Geoghegan et al., 2019). Together, these findings suggest that both osteoblasts and osteocytes play a role in enhancing osteoclastogenesis during estrogen deficiency, which may compound the direct effects of estrogen deficiency on osteoclasts. Interestingly, in this study, RANKL was upregulated under EW regardless of whether or not the constructs were mechanically stimulated. This is consistent with the increased RANKL mRNA expression by osteocytes mentioned earlier, which was demonstrated in static conditions (Geoghegan et al., 2019). Also, we report in this study an increase of dendritic cells under EW in both static and mechanically stimulated conditions. Thus, it is suggested that EW may itself induce these effects. We propose that there may be an increase in available osteocyte membrane-bound RANKL under EW, due to both the increased RANKL mRNA expression and increased number of dendrites we have observed in this study. It is established that estrogen regulates bone resorption, in part, by modulating the expression of membrane-bound RANKL on dendrites of osteocytes, which bind with RANK on osteoclast precursors and in this way induce osteoclastogenesis (Figures 10A,B). Based on our results, we hypothesize that EW may enhance osteoclastogenesis as a result of the interaction of RANKL with RANK (Honma et al., 2014) (Figure 10B). Interestingly, it has previously been proposed that estrogen deficiency results in an increase in PGE2 production in bone marrow stromal cells, which then induces the expression of RANKL leading to accelerated osteoclastogenesis through interaction with RANK (Kanematsu et al., 2000). Here, we propose that estrogen deficiency results in an increase in number of dendrites, which then would make more RANKL available to bind with RANK and ultimately lead to accelerated osteoclastogenesis.
Figure 10. Schematic representation of the proposed mechanism for accelerated osteoclastogenesis in estrogen withdrawal. (A) RANK on osteoclast precursor binds with RANKL on dendrites of osteocytes under normal estrogen conditions. (B) Estrogen deficiency results in an increase in the number of dendrites which suggests increased expression of RANKL leading to accelerated osteoclastogenesis as a result of interaction of RANKL with RANK on osteoclast precursors.
The current study reveals a link between EW, osteocyte mechanosensitivity, differentiation, mineralization, and expression of pro-osteoclastogenic signaling in osteocyte-like cells embedded in 3D hydrogels subjected to mechanical stimulation without growth factors or mineralization media. Importantly, we report for the first time that osteoblast-osteocyte differentiation increased under estrogen deficiency, as confirmed by actin staining, quantification of interconnected dendrites, and DMP1 staining. These findings highlight the impact of estrogen deficiency on bone cell function and provide a novel in vitro osteocyte model to investigate the mechanisms underpinning changes in bone cells after estrogen deficiency.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.
Author Contributions
SN and LM conceived the experimental design. SN, JP, and LM drafted the manuscript. SN and JP performed the cell culture. SN and VK performed biochemical analyses, histological analyses, and immunofluorescence. SN and VK performed phase contrast and confocal imaging. JP performed real-time qPCR and quantification of confocal images. AV performed micro-CT scanning. LM supervised the project. All authors read and approved the final manuscript.
Funding
This study has emanated from research conducted with the financial support of Science Foundation Ireland (SFI) co-funded under the European Regional Development fund (14/IA/2884) and the Irish Research Council (IRC) Consolidator grant (MEMETic).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors acknowledge the facilities and scientific and technical assistance of the Centre for Microscopy and Imaging at the National University of Ireland Galway (www.imaging.nuigalway.ie). The authors would also like to acknowledge Servier Medical Art (www.servier.com) for their image bank used to produce Figure 10.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2020.00601/full#supplementary-material
References
Alford, A. I., Jacobs, C. R., and Donahue, H. J. (2003). Oscillating fluid flow regulates gap junction communication in osteocytic MLO-Y4 cells by an ERK1/2 MAP kinase-dependent mechanism. Bone 33, 64–70. doi: 10.1016/s8756-3282(03)00167-4
Allison, H., and McNamara, L. M. (2019). Inhibition of osteoclastogenesis by mechanically stimulated osteoblasts is attenuated during estrogen deficiency. Am. J. Physiol. Cell. Physiol. 317, C969–C982.
Arnsdorf, E. J., Tummala, P., Kwon, R. Y., and Jacobs, C. R. (2009). Mechanically induced osteogenic differentiation–the role of RhoA, ROCKII and cytoskeletal dynamics. J. Cell. Sci. 122(Pt 4), 546–553. doi: 10.1242/jcs.036293
Atkins, G. J., Welldon, K. J., Wijenayaka, A. R., Bonewald, L. F., and Findlay, D. M. (2009). Vitamin K promotes mineralization, osteoblast-to-osteocyte transition, and an anticatabolic phenotype by {gamma}-carboxylation-dependent and -independent mechanisms. Am. J. Physiol. Cell. Physiol. 297, C1358–C1367.
Bakker, A. D., Klein-Nulend, J., Tanck, E., Heyligers, I. C., Albers, G. H., Lips, P., et al. (2006). Different responsiveness to mechanical stress of bone cells from osteoporotic versus osteoarthritic donors. Osteoporos. Int. 17, 827–833. doi: 10.1007/s00198-006-0072-7
Bakker, A. D., Klein-Nulend, J., Tanck, E., Albers, G. H., Lips, P., and Burger, E. H. (2005). Additive effects of estrogen and mechanical stress on nitric oxide and prostaglandin E2 production by bone cells from osteoporotic donors. Osteoporos. Int. 16, 983–989. doi: 10.1007/s00198-004-1785-0
Balena, R., Toolan, B. C., Shea, M., Markatos, A., Myers, E. R., Lee, S. C., et al. (1993). The effects of 2-year treatment with the aminobisphosphonate alendronate on bone metabolism, bone histomorphometry, and bone strength in ovariectomized nonhuman primates. J. Clin. Invest. 92, 2577–2586. doi: 10.1172/jci116872
Balkan, W., Burnstein, K. L., Schiller, P. C., Perez-Stable, C., D’Ippolito, G., Howard, G. A., et al. (2005). Androgen-induced mineralization by MC3T3-E1 osteoblastic cells reveals a critical window of hormone responsiveness. Biochem. Biophys. Res. Commun. 328, 783–789. doi: 10.1016/j.bbrc.2004.12.090
Baron, R., and Kneissel, M. (2013). WNT signaling in bone homeostasis and disease: from human mutations to treatments. Nat. Med. 19, 179–192. doi: 10.1038/nm.3074
Berthois, Y., Katzenellenbogen, J. A., and Katzenellenbogen, B. S. (1986). Phenol red in tissue culture media is a weak estrogen: implications concerning the study of estrogen-responsive cells in culture. Proc. Natl. Acad. Sci. U.S.A. 83, 2496–2500. doi: 10.1073/pnas.83.8.2496
Birmingham, E., Grogan, J. A., Niebur, G. L., McNamara, L. M., and McHugh, P. E. (2013). Computational modelling of the mechanics of trabecular bone and marrow using fluid structure interaction techniques. Ann. Biomed. Eng. 41, 814–826. doi: 10.1007/s10439-012-0714-1
Birmingham, E., Niebur, G. L., McHugh, P. E., Shaw, G., Barry, F. P., and McNamara, L. M. (2012). Osteogenic differentiation of mesenchymal stem cells is regulated by osteocyte and osteoblast cells in a simplified bone niche. Eur. Cell. Mater. 23, 13–27. doi: 10.22203/ecm.v023a02
Bonewald, L. F. (2011). The amazing osteocyte. J. Bone. Miner. Res. 26, 229–238. doi: 10.1002/jbmr.320
Bonewald, L. F., and Johnson, M. L. (2008). Osteocytes, mechanosensing and Wnt signaling. Bone 42, 606–615. doi: 10.1016/j.bone.2007.12.224
Boskey, A. L., Spevak, L., Paschalis, E., Doty, S. B., and McKee, M. D. (2002). Osteopontin deficiency increases mineral content and mineral crystallinity in mouse bone. Calcif. Tissue Int. 71, 145–154. doi: 10.1007/s00223-001-1121-z
Boukhechba, F., Balaguer, T., Michiels, J. F., Ackermann, K., Quincey, D., Bouler, J. M., et al. (2009). Human primary osteocyte differentiation in a 3D culture system. J. Bone Miner. Res. 24, 1927–1935. doi: 10.1359/jbmr.090517
Boyde, A., Compston, J. E., Reeve, J., Bell, K. L., Noble, B. S., Jones, S. J., et al. (1998). Effect of estrogen suppression on the mineralization density of iliac crest biopsies in young women as assessed by backscattered electron imaging. Bone 22, 241–250. doi: 10.1016/s8756-3282(97)00275-5
Brennan, M. A., Gleeson, J. P., Browne, M., O’Brien, F. J., Thurner, P. J., and McNamara, L. M. (2011a). Site specific increase in heterogeneity of trabecular bone tissue mineral during oestrogen deficiency. Eur. Cell. Mater. 21, 396–406. doi: 10.22203/ecm.v021a30
Brennan, M. A., Gleeson, J. P., O’Brien, F. J., and McNamara, L. M. (2014a). Effects of ageing, prolonged estrogen deficiency and zoledronate on bone tissue mineral distribution. J. Mech. Behav. Biomed. Mater. 29, 161–170. doi: 10.1016/j.jmbbm.2013.08.029
Brennan, M. A., Haugh, M. G., O’Brien, F. J., and McNamara, L. M. (2014b). Estrogen Withdrawal from osteoblasts and osteocytes causes increased mineralization and apoptosis. Horm. Metab. Res. 46, 537–545. doi: 10.1055/s-0033-1363265
Brennan, M. A., McDermott, A., and McNamara, L. M. (2012a). The Impact of Shear Stress during Estrogen Withdrawal on Cell Viability and Mineral Production of Osteoblasts. San Francisco, USA: Orthopaedic Research Society.
Brennan, O., Kennedy, O. D., Lee, T. C., Rackard, S. M., and O’Brien, F. J. (2009). Biomechanical properties across trabeculae from the proximal femur of normal and ovariectomised sheep. J. Biomech. 42, 498–503. doi: 10.1016/j.jbiomech.2008.11.032
Brennan, O., Kennedy, O. D., Lee, T. C., Rackard, S. M., O’Brien, F. J., and McNamara, L. M. (2011b). The effects of estrogen deficiency and bisphosphonate treatment on tissue mineralisation and stiffness in an ovine model of osteoporosis. J. Biomech. 44, 386–390. doi: 10.1016/j.jbiomech.2010.10.023
Brennan, O., Kuliwaba, J., Lee, T., Parkinson, I., Fazzalari, N., McNamara, L., et al. (2012b). Temporal changes in bone composition, architecture and strength following estrogen deficiency in osteoporosis. Calcif. Tissue Int. 91, 440–449. doi: 10.1007/s00223-012-9657-7
Brennan, O., O’Brien, F., and McNamara, L. (2012c). Estrogen plus estrogen receptor antagonists alter mineral production by osteoblasts in vitro. Horm. Metab. Res. 44, 47–53. doi: 10.1055/s-0031-1291358
Buckley, K. A., and Fraser, W. D. (2002). Receptor activator for nuclear factor kappaB ligand and osteoprotegerin: regulators of bone physiology and immune responses/potential therapeutic agents and biochemical markers. Ann. Clin. Biochem. 39(Pt 6), 551–556. doi: 10.1177/000456320203900602
Cabahug-Zuckerman, P., Stout, R. F. Jr., Majeska, R. J., Thi, M. M., Spray, D. C., Weinbaum, S., et al. (2018). Potential role for a specialized beta3 integrin-based structure on osteocyte processes in bone mechanosensation. J. Orthop. Res. 36, 642–652.
Castillo, A. B., Triplett, J. W., Pavalko, F. M., and Turner, C. H. (2014). Estrogen receptor-beta regulates mechanical signaling in primary osteoblasts. Am. J. Physiol. Endocrinol. Metab. 306, E937–E944.
Chen, H.-C., and Hu, Y.-C. (2006). Bioreactors for tissue engineering. Biotechnol. Lett. 28, 1415–1423.
Ciarelli, T. E., Fyhrie, D. P., and Parfitt, A. M. (2003). Effects of vertebral bone fragility and bone formation rate on the mineralization levels of cancellous bone from white females. Bone 32, 311–315. doi: 10.1016/s8756-3282(02)00975-4
Coughlin, T. R., and Niebur, G. L. (2012). Fluid shear stress in trabecular bone marrow due to low-magnitude high-frequency vibration. J Biomech. 45, 2222–2229. doi: 10.1016/j.jbiomech.2012.06.020
Deepak, V., Kayastha, P., and McNamara, L. M. (2017). Estrogen deficiency attenuates fluid flow-induced [Ca2+]i oscillations and mechanoresponsiveness of MLO-Y4 osteocytes. FASEB J. 31, 3027–3039. doi: 10.1096/fj.201601280r
Delgado-Calle, J., Sato, A. Y., and Bellido, T. (2017). Role and mechanism of action of sclerostin in bone. Bone 96, 29–37. doi: 10.1016/j.bone.2016.10.007
Denhardt, D. T., and Noda, M. (1998). Osteopontin expression and function: role in bone remodeling. J. Cell Biochem. Suppl. 3, 92–102. doi: 10.1002/(sici)1097-4644(1998)72:30/31
Dickenson, R. P., Hutton, W. C., and Stott, J. R. (1981). The mechanical properties of bone in osteoporosis. J. Bone Joint Surg. Br. 63-B, 233–238.
Divieti Pajevic, P., and Krause, D. S. (2019). Osteocyte regulation of bone and blood. Bone 119, 13–18. doi: 10.1016/j.bone.2018.02.012
Efstathiadou, Z., Tsatsoulis, A., and Ioannidis, J. P. (2001). Association of collagen Ialpha 1 Sp1 polymorphism with the risk of prevalent fractures: a meta-analysis. J. Bone Miner. Res. 16, 1586–1592. doi: 10.1359/jbmr.2001.16.9.1586
Ek-Rylander, B., Flores, M., Wendel, M., Heinegard, D., and Andersson, G. (1994). Dephosphorylation of osteopontin and bone sialoprotein by osteoclastic tartrate-resistant acid phosphatase. Modulation of osteoclast adhesion in vitro. J. Biol. Chem. 269, 14853–14856.
Emerton, K. B., Hu, B., Woo, A. A., Sinofsky, A., Hernandez, C., Majeska, R. J., et al. (2010). Osteocyte apoptosis and control of bone resorption following ovariectomy in mice. Bone 46, 577–583. doi: 10.1016/j.bone.2009.11.006
Feng, J. Q., Zhang, J., Dallas, S. L., Lu, Y., Chen, S., Tan, X., et al. (2002). “Dentin matrix protein 1, a target molecule for cbfa1 in bone, is a unique bone marker gene. J. Bone Miner. Res. 17, 1822–1831. doi: 10.1359/jbmr.2002.17.10.1822
Florencio-Silva, R., Sasso, G. R. S., Sasso-Cerri, E., Simões, M. J., and Cerri, P. S. (2018). Effects of estrogen status in osteocyte autophagy and its relation to osteocyte viability in alveolar process of ovariectomized rats. Biomed. Pharmacother. 98, 406–415. doi: 10.1016/j.biopha.2017.12.089
Freeman, F. E., Haugh, M. G., and McNamara, L. M. (2015). An in vitro bone tissue regeneration strategy combining chondrogenic and vascular priming enhances the mineralization potential of mesenchymal stem cells in vitro while also allowing for vessel formation. Tissue Eng. Part A. 21, 1320–1332. doi: 10.1089/ten.tea.2014.0249
Freeman, F. E., Schiavi, J., Brennan, M. A., Owens, P., Layrolle, P., and McNamara, L. M. (2017). (∗) Mimicking the biochemical and mechanical extracellular environment of the endochondral ossification process to enhance the in vitro mineralization potential of human mesenchymal stem cells. Tissue Eng. Part A. 23, 1466–1478. doi: 10.1089/ten.tea.2017.0052
Gadeleta, S. J., Boskey, A. L., Paschalis, E., Carlson, C., Menschik, F., Baldini, T., et al. (2000). A physical, chemical, and mechanical study of lumbar vertebrae from normal, ovariectomized, and nandrolone decanoate-treated cynomolgus monkeys (Macaca fascicularis). Bone 27, 541–550. doi: 10.1016/s8756-3282(00)00362-8
Geoghegan, I. P., Hoey, D. A., and McNamara, L. M. (2019). Estrogen deficiency impairs integrin αvβ3-mediated mechanosensation by osteocytes and alters osteoclastogenic paracrine signalling. Sci. Rep. 9:4654.
Giachelli, C. M., and Steitz, S. (2000). Osteopontin: a versatile regulator of inflammation and biomineralization. Matrix Biol. 19, 615–622. doi: 10.1016/s0945-053x(00)00108-6
Hagenmüller, H., Hofmann, S., Kohler, T., Merkle, H. P., Kaplan, D. L., Vunjak-Novakovic, G., et al. (2007). Non-invasive time-lapsed monitoring and quantification of engineered bone-like tissue. Ann. Biomed. Eng. 35, 1657–1667. doi: 10.1007/s10439-007-9338-2
Haisenleder, D. J., Schoenfelder, A. H., Marcinko, E. S., Geddis, L. M., and Marshall, J. C. (2011). Estimation of estradiol in mouse serum samples: evaluation of commercial estradiol immunoassays. Endocrinology 152, 4443–4447. doi: 10.1210/en.2011-1501
Han, Y., You, X., Xing, W., Zhang, Z., and Zou, W. (2018). Paracrine and endocrine actions of bone—the functions of secretory proteins from osteoblasts, osteocytes, and osteoclasts. Bone Res. 6:16.
Harmey, D., Johnson, K. A., Zelken, J., Camacho, N. P., Hoylaerts, M. F., Noda, M., et al. (2006). Elevated skeletal osteopontin levels contribute to the hypophosphatasia phenotype in Akp2(-/-) mice. J. Bone Miner. Res. 21, 1377–1386. doi: 10.1359/jbmr.060619
Harris, S. S., Patel, M. S., Cole, D. E., and Dawson-Hughes, B. (2000). Associations of the collagen type Ialpha1 Sp1 polymorphism with five-year rates of bone loss in older adults. Calcif. Tissue Int. 66, 268–271. doi: 10.1007/pl00005842
Haugh, M. G., Meyer, E. G., Thorpe, S. D., Vinardell, T., Duffy, G. P., and Kelly, D. J. (2011). Temporal and spatial changes in cartilage-matrix-specific gene expression in mesenchymal stem cells in response to dynamic compression. Tissue Eng. Part A. 17, 3085–3093. doi: 10.1089/ten.tea.2011.0198
He, L., Hao, Y., Zhen, L., Liu, H., Shao, M., Xu, X., et al. (2019). Biomineralization of dentin. J. Struc. Biol. 207, 115–122.
Hemmatian, H., Jalali, R., Semeins, C. M., Hogervorst, J. M. A., van Lenthe, G. H., Klein-Nulend, J., et al. (2018). mechanical loading differentially affects osteocytes in fibulae from lactating mice compared to osteocytes in virgin mice: possible role for lacuna size. Calcif. Tissue Int. 103, 675–685. doi: 10.1007/s00223-018-0463-8
Henstock, J. R., Rotherham, M., Rose, J. B., and El Haj, A. J. (2013). Cyclic hydrostatic pressure stimulates enhanced bone development in the foetal chick femur in vitro. Bone 53, 468–477. doi: 10.1016/j.bone.2013.01.010
Hoac, B., Kiffer-Moreira, T., Millán, J. L., and McKee, M. D. (2013). Polyphosphates inhibit extracellular matrix mineralization in MC3T3-E1 osteoblast cultures. Bone 53, 478–486. doi: 10.1016/j.bone.2013.01.020
Hoey, D. A., Tormey, S., Ramcharan, S., O’Brien, F. J., and Jacobs, C. R. (2012). Primary cilia-mediated mechanotransduction in human mesenchymal stem cells. Stem. Cells 30, 2561–2570. doi: 10.1002/stem.1235
Holdsworth, G., Roberts, S. J., and Ke, H. Z. (2019). Novel actions of sclerostin on bone. J. Mol. Endocrinol. 62, R167–R185.
Honma, M., Ikebuchi, Y., Kariya, Y., and Suzuki, H. (2014). Regulatory mechanisms of RANKL presentation to osteoclast precursors. Curr. Osteoporos. Rep. 12, 115–120. doi: 10.1007/s11914-014-0189-0
Hughes, D. E., Salter, D. M., Dedhar, S., and Simpson, R. (1993). Integrin expression in human bone. J. Bone Miner. Res. 8, 527–533. doi: 10.1002/jbmr.5650080503
Icer, M. A., and Gezmen-Karadag, M. (2018). The multiple functions and mechanisms of osteopontin. Clin. Biochem. 59, 17–24. doi: 10.1016/j.clinbiochem.2018.07.003
Jacobs, C. R., Temiyasathit, S., and Castillo, A. B. (2010). Osteocyte mechanobiology and pericellular mechanics. Annu. Rev. Biomed. Eng. 12, 369–400. doi: 10.1146/annurev-bioeng-070909-105302
Kanematsu, M., Sato, T., Takai, H., Watanabe, K., Ikeda, K., and Yamada, Y. (2000). Prostaglandin E2 induces expression of receptor activator of nuclear factor-kappa B ligand/osteoprotegrin ligand on pre-B cells: implications for accelerated osteoclastogenesis in estrogen deficiency. J. Bone Miner. Res. 15, 1321–1329. doi: 10.1359/jbmr.2000.15.7.1321
Kapur, S., Baylink, D. J., and Lau, K. H. W. (2002). Fluid flow shear stress stimulates human osteoblast proliferation and differentiation through multiple interacting and competing signal transduction pathways. Bone 32, 241–251. doi: 10.1016/s8756-3282(02)00979-1
Kennedy, O. D., Brennan, O., Rackard, S. M., Staines, A., O’Brien, F. J., Taylor, D., et al. (2009). Effects of ovariectomy on bone turnover, porosity, and biomechanical properties in ovine compact bone 12 months postsurgery. J. Orthop. Res. 27, 303–309. doi: 10.1002/jor.20750
Klein-Nulend, J., Semeins, C. M., and Burger, E. H. (1996). Prostaglandin mediated modulation of transforming growth factor-beta metabolism in primary mouse osteoblastic cells in vitro. J. Cell. Physiol. 168, 1–7. doi: 10.1002/(sici)1097-4652(199607)168:1<1::aid-jcp1>3.0.co;2-t
Klein-Nulend, J., Veldhuijzen, J. P., and Burger, E. H. (1986). Increased calcification of growth plate cartilage as a result of compressive force in vitro. Arthritis Rheum. 29, 1002–1009. doi: 10.1002/art.1780290809
Knothe Tate, M. L. (2003). “Whither flows the fluid in bone?” An osteocyte’s perspective. J. Biomech. 36, 1409–1424. doi: 10.1016/s0021-9290(03)00123-4
Koster, N., Schmiermund, A., Grubelnig, S., Leber, J., Ehlicke, F., Czermak, P., et al. (2016). Single-step RNA extraction from different hydrogel-embedded mesenchymal stem cells for quantitative reverse transcription-polymerase chain reaction analysis. Tissue Eng. Part C. Methods 22, 552–560. doi: 10.1089/ten.tec.2015.0362
Kousteni, S., Bellido, T., Plotkin, L. I., O’Brien, C. A., Bodenner, D. L., Han, L., et al. (2001). Nongenotropic, sex-nonspecific signaling through the estrogen or androgen receptors: dissociation from transcriptional activity. Cell 104, 719–730. doi: 10.1016/s0092-8674(02)08100-x
Kousteni, S., Han, L., Chen, J.-R., Almeida, M., Plotkin, L. I., Bellido, T., et al. (2003). Kinase-mediated regulation of common transcription factors accounts for the bone-protective effects of sex steroids. J. Clin. Invest. 111, 1651–1664. doi: 10.1172/jci200317261
Kramer, I., Halleux, C., Keller, H., Pegurri, M., Gooi, J. H., Weber, P. B., et al. (2010). Osteocyte Wnt/beta-catenin signaling is required for normal bone homeostasis. Mol. Cell. Biol. 30, 3071–3085. doi: 10.1128/mcb.01428-09
Kulkarni, G. V., Chen, B., Malone, J. P., Narayanan, A. S., and George, A. (2000). Promotion of selective cell attachment by the RGD sequence in dentine matrix protein 1. Arch. Oral Biol. 45, 475–484. doi: 10.1016/s0003-9969(00)00010-8
Kyono, A., Avishai, N., Ouyang, Z., Landreth, G. E., and Murakami, S. (2012). FGF and ERK signaling coordinately regulate mineralization-related genes and play essential roles in osteocyte differentiation. J. Bone Miner. Metabol. 30, 19–30. doi: 10.1007/s00774-011-0288-2
Lee, K., Jessop, H., Suswillo, R., Zaman, G., and Lanyon, L. (2003). Bone adaptation requires oestrogen receptor-α. Nature 424, 389–389. doi: 10.1038/424389a
Lee, K. L., Hoey, D. A., Spasic, M., Tang, T., Hammond, H. K., and Jacobs, C. R. (2014). Adenylyl cyclase 6 mediates loading-induced bone adaptation in vivo. FASEB J. 28, 1157–1165. doi: 10.1096/fj.13-240432
Lenton, S., Seydel, T., Nylander, T., Holt, C., Härtlein, M., Teixeira, S., et al. (2015). Dynamic footprint of sequestration in the molecular fluctuations of osteopontin. J. R. Soc. Interface 12:20150506. doi: 10.1098/rsif.2015.0506
Lewiecki, E. M. (2011). New targets for intervention in the treatment of postmenopausal osteoporosis. Nat. Rev. Rheumatol. 7, 631–638. doi: 10.1038/nrrheum.2011.130
Li, J., Rose, E., Frances, D., Sun, Y., and You, L. (2012). Effect of oscillating fluid flow stimulation on osteocyte mRNA expression. J. Biomech. 45, 247–251. doi: 10.1016/j.jbiomech.2011.10.037
Litzenberger, J., Kim, J., Tummala, P., and Jacobs, C. (2010). Beta1 integrins mediate mechanosensitive signaling pathways in osteocytes. Calcif. Tissue Int. 86, 325–332. doi: 10.1007/s00223-010-9343-6
Liu, C., Zhao, Y., Cheung, W.-Y., Gandhi, R., Wang, L., and You, L. (2010). Effects of cyclic hydraulic pressure on osteocytes. Bone 46, 1449–1456. doi: 10.1016/j.bone.2010.02.006
Luppen, C. A., Smith, E., Spevak, L., Boskey, A. L., and Frenkel, B. (2003). Bone morphogenetic protein-2 restores mineralization in glucocorticoid-inhibited MC3T3-E1 osteoblast cultures. J. Bone Miner. Res. 18, 1186–1197. doi: 10.1359/jbmr.2003.18.7.1186
Malaval, L., Wade-Guéye, N. M., Boudiffa, M., Fei, J., Zirngibl, R., Chen, F., et al. (2008). Bone sialoprotein plays a functional role in bone formation and osteoclastogenesis. J. Exp. Med. 205, 1145–1153. doi: 10.1084/jem.20071294
Malone, A. M., Anderson, C. T., Tummala, P., Kwon, R. Y., Johnston, T. R., Stearns, T., et al. (2007). Primary cilia mediate mechanosensing in bone cells by a calcium-independent mechanism. Proc. Natl. Acad. Sci. U.S.A. 104, 13325–13330. doi: 10.1073/pnas.0700636104
Manju, V., Anitha, A., Menon, D., Iyer, S., Nair, S. V., and Nair, M. B. (2018). Nanofibrous yarn reinforced HA-gelatin composite scaffolds promote bone formation in critical sized alveolar defects in rabbit model. Biomed. Mater. 13:065011. doi: 10.1088/1748-605x/aadf99
Mansell, J. P., and Bailey, A. J. (2003). Increased metabolism of bone collagen in post-menopausal female osteoporotic femoral heads. Int. J. Biochem. Cell Biol. 35, 522–529. doi: 10.1016/s1357-2725(02)00312-6
Marathe, N., Rangaswami, H., Zhuang, S., Boss, G. R., and Pilz, R. B. (2012). Pro-survival effects of 17beta-estradiol on osteocytes are mediated by nitric oxide/cGMP via differential actions of cGMP-dependent protein kinases I and II. J. Biol. Chem. 287, 978–988. doi: 10.1074/jbc.m111.294959
Mc Garrigle, M. J., Mullen, C. A., Haugh, M. G., Voisin, M. C., and McNamara, L. M. (2016). Osteocyte differentiation and the formation of an interconnected cellular network in vitro. Eur. Cell. Mater. 31, 323–340. doi: 10.22203/ecm.v031a21
McAllister, T. (2000). Fluid shear stress stimulates prostaglandin and nitric oxide release in bone marrow-derived preosteoclast-like cells. Biochem. Biophys. Res. Commun. 270, 643–648. doi: 10.1006/bbrc.2000.2467
McCreadie, B. R., Morris, M. D., Chen, T. C., Sudhaker Rao, D., Finney, W. F., Widjaja, E., et al. (2006). Bone tissue compositional differences in women with and without osteoporotic fracture. Bone 39, 1190–1195. doi: 10.1016/j.bone.2006.06.008
McKee, M. D., and Nanci, A. (1996). Osteopontin at mineralized tissue interfaces in bone, teeth, and osseointegrated implants: ultrastructural distribution and implications for mineralized tissue formation, turnover, and repair. Microsc. Res. Tech. 33, 141–164. doi: 10.1002/(sici)1097-0029(19960201)33:2<141::aid-jemt5>3.0.co;2-w
McNamara, L. M. (2010). Perspective on post-menopausal osteoporosis: establishing an interdisciplinary understanding of the sequence of events from the molecular level to whole bone fractures. J. R. Soc. Interface 7, 353–372. doi: 10.1098/rsif.2009.0282
McNamara, L. M., Ederveen, A. G., Lyons, C. G., Price, C., Schaffler, M. B., Weinans, H., et al. (2006). Strength of cancellous bone trabecular tissue from normal, ovariectomized and drug-treated rats over the course of ageing. Bone 39, 392–400. doi: 10.1016/j.bone.2006.02.070
McNamara, L. M., Majeska, R. J., Weinbaum, S., Friedrich, V., and Schaffler, M. B. (2009). Attachment of osteocyte cell processes to the bone matrix. Anat. Rec. 292, 355–363. doi: 10.1002/ar.20869
McNamara, L. M., Prendergast, P. J., and Schaffler, M. B. (2005). Bone tissue material properties are altered during osteoporosis. J. Musculoskelet. Neuronal Interact. 5, 342–343.
Melke, J., Midha, S., Ghosh, S., Ito, K., and Hofmann, S. (2016). Silk fibroin as biomaterial for bone tissue engineering. Acta Biomater. 31, 1–16. doi: 10.1016/j.actbio.2015.09.005
Melton, L. J. III, Chrischilles, E. A., Cooper, C., Lane, A. W., and Riggs, B. L. (2005). How many women have osteoporosis? JBMR anniversary classic. JBMR, volume 7, Number 9, 1992. J. Bone Miner. Res. 20, 886–892.
Mikuni-Takagaki, Y., Kakai, Y., Satoyoshi, M., Kawano, E., Suzuki, Y., Kawase, T., et al. (1995). Matrix mineralization and the differentiation of osteocyte-like cells in culture. J. Bone Miner. Res. 10, 231–242. doi: 10.1002/jbmr.5650100209
Moreno-Cuevas, J. E., and Sirbasku, D. A. (2000). Estrogen mitogenic action. III. Is phenol red a “red herring”? In Vitro Cell. Dev. Biol. Anim. 36, 447–464.
Moustafa, A., Sugiyama, T., Prasad, J., Zaman, G., Gross, T. S., Lanyon, L. E., et al. (2012). Mechanical loading-related changes in osteocyte sclerostin expression in mice are more closely associated with the subsequent osteogenic response than the peak strains engendered. Osteoporos. Int. 23, 1225–1234. doi: 10.1007/s00198-011-1656-4
Mullen, C. A., Haugh, M. G., Schaffler, M. B., Majeska, R. J., and McNamara, L. M. (2013). Osteocyte differentiation is regulated by extracellular matrix stiffness and intercellular separation. J. Mech. Behav. Biomed. Mater. 28, 183–194. doi: 10.1016/j.jmbbm.2013.06.013
Mullen, C. A., Vaughan, T. J., Voisin, M. C., Brennan, M. A., Layrolle, P., and McNamara, L. M. (2014). Cell morphology and focal adhesion location alters internal cell stress. J. R. Soc. Interface 11:20140885. doi: 10.1098/rsif.2014.0885
Nakano, Y., Addison, W. N., and Kaartinen, M. T. (2007). ATP-mediated mineralization of MC3T3-E1 osteoblast cultures. Bone 41, 549–561. doi: 10.1016/j.bone.2007.06.011
Nakano, Y., Beertsen, W., van den Bos, T., Kawamoto, T., Oda, K., and Takano, Y. (2004). Site-specific localization of two distinct phosphatases along the osteoblast plasma membrane: tissue non-specific alkaline phosphatase and plasma membrane calcium ATPase. Bone 35, 1077–1085. doi: 10.1016/j.bone.2004.07.009
Nakashima, T., Hayashi, M., Fukunaga, T., Kurata, K., Oh-hora, M., Feng, J. Q., et al. (2011). Evidence for osteocyte regulation of bone homeostasis through RANKL expression. Nat. Med. 17, 1231–1234. doi: 10.1038/nm.2452
Narayanan, K., Ramachandran, A., Hao, J., He, G., Park, K. W., Cho, M., et al. (2003). Dual functional roles of dentin matrix protein 1. Implications in biomineralization and gene transcription by activation of intracellular Ca2+ store. J. Biol. Chem. 278, 17500–17508. doi: 10.1074/jbc.m212700200
Nauman, E. A., Satcher, R. L., Keaveny, T. M., Halloran, B. P., and Bikle, D. D. (2001). Osteoblasts respond to pulsatile fluid flow with short- term increases in PGE2 but no change in mineralization. J. Appl. Physiol. 90, 1849–1854. doi: 10.1152/jappl.2001.90.5.1849
O’Sullivan, L. M., Allison, H., Parle, E. E., Schiavi, J., and McNamara, L. M. (2019). Secondary alterations in bone mineralisation and trabecular thickening occur after long-term estrogen deficiency in ovariectomised rat tibiae, which do not coincide with initial rapid bone loss. Osteoporos. Int. 31, 587–599. doi: 10.1007/s00198-019-05239-5
Padovano, J. D., Ravindran, S., Snee, P. T., Ramachandran, A., Bedran-Russo, A. K., and George, A. (2015). DMP1-derived peptides promote remineralization of human dentin. J. Dent. Res. 94, 608–614. doi: 10.1177/0022034515572441
Papanicolaou, S. E., Phipps, R. J., Fyhrie, D. P., and Genetos, D. C. (2009). Modulation of sclerostin expression by mechanical loading and bone morphogenetic proteins in osteogenic cells. Biorheology 46, 389–399. doi: 10.3233/bir-2009-0550
Parle, E., Tio, S., Behre, A., Carey, J. J., Murphy, C. G., O’Brien, T. F., et al. (2020). Bone mineral is more heterogeneously distributed in the femoral heads of osteoporotic and diabetic patients: a pilot study. JBMR Plus 4:e10253.
Pfaffl, M. W. (2001). A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 29:e45.
Plotkin, L. I., Aguirre, J. I., Kousteni, S., Manolagas, S. C., and Bellido, T. (2005). Bisphosphonates and estrogens inhibit osteocyte apoptosis via distinct molecular mechanisms downstream of extracellular signal-regulated kinase activation. J. Biol. Chem. 280, 7317–7325. doi: 10.1074/jbc.m412817200
Poole, K. E. S., van Bezooijen, R. L., Loveridge, N., Hamersma, H., Papapoulos, S. E., Lowik, C. W., et al. (2005). Sclerostin is a delayed secreted product of osteocytes that inhibits bone formation. Faseb. J. 19, 1842–1844. doi: 10.1096/fj.05-4221fje
Porter, K. L., Chanda, S., Wang, H. Q., Gaido, K. W., Smart, R. C., and Robinette, C. L. (2002). 17beta-estradiol is a hormonal regulator of mirex tumor promotion sensitivity in mice. Toxicol. Sci. 69, 42–48. doi: 10.1093/toxsci/69.1.42
Pörtner, R., Nagel-Heyer, S., Goepfert, C., Adamietz, P., and Meenen, N. M. (2005). Bioreactor design for tissue engineering. J. Biosci. Bioeng. 100, 235–245.
Raina, D. B., Larsson, D., Mrkonjic, F., Isaksson, H., Kumar, A., Lidgren, L., et al. (2018). Gelatin- hydroxyapatite- calcium sulphate based biomaterial for long term sustained delivery of bone morphogenic protein-2 and zoledronic acid for increased bone formation: in-vitro and in-vivo carrier properties. J. Control Release 272, 83–96. doi: 10.1016/j.jconrel.2018.01.006
Randell, K. M., Honkanen, R. J., Kroger, H., and Saarikoski, S. (2002). Does hormone-replacement therapy prevent fractures in early postmenopausal women? J. Bone Miner. Res. 17, 528–533. doi: 10.1359/jbmr.2002.17.3.528
Reinwald, Y., and El Haj, A. J. (2018). Hydrostatic pressure in combination with topographical cues affects the fate of bone marrow-derived human mesenchymal stem cells for bone tissue regeneration. J. Biomed. Mater. Res. Part A 106, 629–640. doi: 10.1002/jbm.a.36267
Reinwald, Y., Leonard, K. H. L., Henstock, J. R., Whiteley, J. P., Osborne, J. M., Waters, S. L., et al. (2015). Evaluation of the growth environment of a hydrostatic force bioreactor for preconditioning of tissue-engineered constructs. Tissue Eng. Part C. Methods 21, 1–14. doi: 10.1089/ten.tec.2013.0476
Ren, J., Wang, X. H., Wang, G. C., and Wu, J. H. (2013). 17beta estradiol regulation of connexin 43-based gap junction and mechanosensitivity through classical estrogen receptor pathway in osteocyte-like MLO-Y4 cells. Bone 53, 587–596. doi: 10.1016/j.bone.2012.12.004
Riggs, B. L., Khosla, S., and Melton, L. J. III (2002). Sex steroids and the construction and conservation of the adult skeleton. Endocr. Rev. 23, 279–302. doi: 10.1210/edrv.23.3.0465
Rios, H. F., Ye, L., Dusevich, V., Eick, D., Bonewald, L. F., and Feng, J. Q. (2005). DMP1 is essential for osteocyte formation and function. J. Musculoskelet Neuronal Interact 5, 325–327.
Robling, A. G., Niziolek, P. J., Baldridge, L. A., Condon, K. W., Allen, M. R., Alam, I., et al. (2008). Mechanical stimulation of bone in vivo reduces osteocyte expression of Sost/sclerostin. J. Biol. Chem. 283, 5866–5875. doi: 10.1074/jbc.m705092200
Roelofsen, J., Klein-Nulend, J., and Burger, E. H. (1995). Mechanical stimulation by intermittent hydrostatic compression promotes bone-specific gene expression in vitro. J. Biomech. 28, 1493–1503. doi: 10.1016/0021-9290(95)00097-6
Ross, F. P., Chappel, J., Alvarez, J. I., Sander, D., Butler, W. T., Farach-Carson, M. C., et al. (1993). Interactions between the bone matrix proteins osteopontin and bone sialoprotein and the osteoclast integrin alpha v beta 3 potentiate bone resorption. J. Biol. Chem. 268, 9901–9907.
Rubin, J., Rubin, C., and Jacobs, C. R. (2006). Molecular pathways mediating mechanical signaling in bone. Gene 367, 1–16. doi: 10.1016/j.gene.2005.10.028
Ruijtenberg, S., and van den Heuvel, S. (2016). Coordinating cell proliferation and differentiation: Antagonism between cell cycle regulators and cell type-specific gene expression. Cell Cycle 15, 196–212. doi: 10.1080/15384101.2015.1120925
Saxon, L. K., Galea, G., Meakin, L., Price, J., and Lanyon, L. E. (2012). Estrogen receptors α and β have different gender-dependent effects on the adaptive responses to load bearing in cancellous and cortical bone. Endocrinology 153, 2254–2266. doi: 10.1210/en.2011-1977
Siyam, A., Wang, S., Qin, C., Mues, G., Stevens, R., D’Souza, R. N., et al. (2012). Nuclear localization of DMP1 proteins suggests a role in intracellular signaling. Biochem. Biophys. Res. Commun. 424, 641–646. doi: 10.1016/j.bbrc.2012.07.037
Spatz, J. M., Wein, M. N., Gooi, J. H., Qu, Y., Garr, J. L., Liu, S., et al. (2015). The Wnt inhibitor sclerostin is up-regulated by mechanical unloading in osteocytes in vitro. J. Biol. Chem. 290, 16744–16758. doi: 10.1074/jbc.m114.628313
Sterck, J. G., Klein-Nulend, J., Lips, P., and Burger, E. H. (1998). Response of normal and osteoporotic human bone cells to mechanical stress in vitro. Am. J. Physiol. 274(6 Pt 1), E1113–E1120.
Sudo, H., Kodama, H. A., Amagai, Y., Yamamoto, S., and Kasai, S. (1983). In vitro differentiation and calcification in a new clonal osteogenic cell line derived from newborn mouse calvaria. J. Cell Biol. 96, 191–198. doi: 10.1083/jcb.96.1.191
Thompson, W. R., Uzer, G., Brobst, K. E., Xie, Z., Sen, B., Yen, S. S., et al. (2015). Osteocyte specific responses to soluble and mechanical stimuli in a stem cell derived culture model. Sci. Rep. 5:11049.
Tomkinson, A., Gevers, E. F., Wit, J. M., Reeve, J., and Noble, B. S. (1998). The role of estrogen in the control of rat osteocyte apoptosis. J. Bone Miner. Res. 13, 1243–1250. doi: 10.1359/jbmr.1998.13.8.1243
Tomkinson, A., Reeve, J., Shaw, R. W., and Noble, B. S. (1997). The Death of Osteocytes via Apoptosis Accompanies Estrogen Withdrawal in Human Bone∗. J. Clin. Endocrinol. Metabol. 82, 3128–3135. doi: 10.1210/jc.82.9.3128
Tu, X. L., Rhee, Y., Condon, K. W., Bivi, N., Allen, M. R., Dwyer, D., et al. (2012). Sost downregulation and local Wnt signaling are required for the osteogenic response to mechanical loading. Bone 50, 209–217. doi: 10.1016/j.bone.2011.10.025
Uchihashi, K., Aoki, S., Matsunobu, A., and Toda, S. (2013). Osteoblast migration into type I collagen gel and differentiation to osteocyte-like cells within a self-produced mineralized matrix: A novel system for analyzing differentiation from osteoblast to osteocyte. Bone 52, 102–110. doi: 10.1016/j.bone.2012.09.001
Vaughan, T. J., Voisin, M., Niebur, G. L., and McNamara, L. M. (2015). Multiscale modeling of trabecular bone marrow: understanding the micromechanical environment of mesenchymal stem cells during osteoporosis. J. Biomech. Eng. 137. doi: 10.1115/1.4028986
Verbruggen, S. W., Mc Garrigle, M. J., Haugh, M. G., Voisin, M. C., and McNamara, L. M. (2015). Altered mechanical environment of bone cells in an animal model of short- and long-term osteoporosis. Biophys. J. 108, 1587–1598. doi: 10.1016/j.bpj.2015.02.031
Verbruggen, S. W., Vaughan, T. J., and McNamara, L. M. (2016). Mechanisms of osteocyte stimulation in osteoporosis. J. Mech. Behav. Biomed. Mater. 62, 158–168. doi: 10.1016/j.jmbbm.2016.05.004
Vetsch, J. R., Betts, D. C., Muller, R., and Hofmann, S. (2017). Flow velocity-driven differentiation of human mesenchymal stromal cells in silk fibroin scaffolds: a combined experimental and computational approach. PLoS One 12:e0180781. doi: 10.1371/journal.pone.0180781
Vetsch, J. R., Paulsen, S. J., Muller, R., and Hofmann, S. (2015). Effect of fetal bovine serum on mineralization in silk fibroin scaffolds. Acta Biomater. 13, 277–285. doi: 10.1016/j.actbio.2014.11.025
Voisin, M., and McNamara, L. M. (2015). Differential beta and beta integrin expression in bone marrow and cortical bone of estrogen deficient rats. Anat. Rec. 298, 1548–1559. doi: 10.1002/ar.23173
Wada, T., Nakashima, T., Hiroshi, N., and Penninger, J. M. (2006). RANKL-RANK signaling in osteoclastogenesis and bone disease. Trends Mol. Med. 12, 17–25. doi: 10.1016/j.molmed.2005.11.007
Wang, W., Li, F., Wang, K., Cheng, B., and Guo, X. (2012). PAPSS2 promotes alkaline phosphates activity and mineralization of osteoblastic MC3T3-E1 cells by crosstalk and Smads signal pathways. PloS One 7:e43475. doi: 10.1371/journal.pone.0043475
Wang, Y., McNamara, L. M., Schaffler, M. B., and Weinbaum, S. (2008). Strain amplification and integrin based signaling in osteocytes. J. Musculoskelet Neuronal Interact 8, 332–334.
Welshons, W. V., Wolf, M. F., Murphy, C. S., and Jordan, V. C. (1988). Estrogenic activity of phenol red. Mol. Cell. Endocrinol. 57, 169–178. doi: 10.1016/0303-7207(88)90072-x
Wijenayaka, A. R., Kogawa, M., Lim, H. P., Bonewald, L. F., Findlay, D. M., and Atkins, G. J. (2011). Sclerostin stimulates osteocyte support of osteoclast activity by a RANKL-dependent pathway. PLoS One 6:e25900. doi: 10.1371/journal.pone.0180781
Woo, S. M., Rosser, J., Dusevich, V., Kalajzic, I., and Bonewald, L. F. (2011). Cell line IDG-SW3 replicates osteoblast-to-late-osteocyte differentiation in vitro and accelerates bone formation in vivo. J. Bone Miner. Res. 26, 2634–2646. doi: 10.1002/jbmr.465
Wood, G. A., Fata, J. E., Watson, K. L., and Khokha, R. (2007). Circulating hormones and estrous stage predict cellular and stromal remodeling in murine uterus. Reproduction 133, 1035–1044. doi: 10.1530/rep-06-0302
Yokose, S., Ishizuya, T., Ikeda, T., Nakamura, T., Tsurukami, H., Kawasaki, K., et al. (1996). An estrogen deficiency caused by ovariectomy increases plasma levels of systemic factors that stimulate proliferation and differentiation of osteoblasts in rats. Endocrinology 137, 469–478. doi: 10.1210/endo.137.2.8593791
Zaidel-Bar, R., Cohen, M., Addadi, L., and Geiger, B. (2004). Hierarchical assembly of cell–matrix adhesion complexes. Biochem. Soci. Transac. 32, 416–420. doi: 10.1042/bst0320416
Zhang, D., Weinbaum, S., and Cowin, S. C. (1998). On the calculation of bone pore water pressure due to mechanical loading. Int. J. Sol. Struc. 35, 4981–4997. doi: 10.1016/s0020-7683(98)00105-x
Zhu, X., Luo, J., Chen, X., Wang, J., Wang, G., Li, H., et al. (2008). Expression characteristic and significance of interleukin-6, nuclear factor kappa beta, and bone formation markers in rat models of osteoporosis. Transl. Res. 152, 18–23. doi: 10.1016/j.trsl.2008.05.003
Ziambaras, K., Lecanda, F., Steinberg, T. H., and Civitelli, R. (1998). Cyclic stretch enhances gap junctional communication between osteoblastic cells. J. Bone Miner. Res. 13, 218–228. doi: 10.1359/jbmr.1998.13.2.218
Keywords: estrogen deficiency, osteoblasts, osteocytes, mineralization, 3D culture, mechanical stimulation
Citation: Naqvi SM, Panadero Pérez JA, Kumar V, Verbruggen ASK and McNamara LM (2020) A Novel 3D Osteoblast and Osteocyte Model Revealing Changes in Mineralization and Pro-osteoclastogenic Paracrine Signaling During Estrogen Deficiency. Front. Bioeng. Biotechnol. 8:601. doi: 10.3389/fbioe.2020.00601
Received: 30 March 2020; Accepted: 18 May 2020;
Published: 10 June 2020.
Edited by:
Bruce Alan Bunnell, Tulane University, United StatesReviewed by:
Lilian Irene Plotkin, Indiana University Bloomington, United StatesAstrid Diana Bakker, VU University Amsterdam, Netherlands
Copyright © 2020 Naqvi, Panadero Pérez, Kumar, Verbruggen and McNamara. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Laoise M. McNamara, bGFvaXNlLm1jbmFtYXJhQG51aWdhbHdheS5pZQ==
†These authors have contributed equally to this work