TY - JOUR AU - Jeanne,, Sylvie AB - Summary Objectives Mechanobiology phenomena constitute a major element of the cellular and tissue response during orthodontic treatment and the implantation of a biomaterial. Better understanding these phenomena will improve the effectiveness of our treatments. The objective of this work is to validate a model of three-dimensional (3D) culture of osteoblasts to study mechanobiology. Materials and methods The hFOB 1.19 cell line was cultured either traditionally on a flat surface or in aggregates called spheroids. They were embedded in 0.8% low-melting agarose type VII and placed in a polyethylene terephthalate transwell insert. Compressive forces of 1 and 4 g/cm2 were applied with an adjustable weight. Proliferation was evaluated by measuring diameters, monitoring glucose levels, and conducting Hoechst/propidium iodide staining. Enzyme-linked immunosorbent assays focusing on the pro-inflammatory mediators interleukin (IL)-6 and IL-8 and bone remodelling factor osteoprotegerin were performed to evaluate soluble factor synthesis. quantitative reverse transcription-polymerase chain reaction was performed to evaluate bone marker transcription. Results The 3D model shows good cell viability and permits IL dosing. Additionally, three gene expression profiles are analysable. Limitations The model allows analysis of conventional markers; larger exploration is needed for better understanding osteoblast mechanobiology. However, it only allows an analysis over 3 days. Conclusion The results obtained by applying constant compressive forces to 3D osteoblastic cultures validate this model system for exploring biomolecule release and analysing gene transcription. In particular, it highlights a disturbance in the expression of markers of osteogenesis. Introduction Orthodontic tooth movement (OTM) is initiated by the interaction of cellular and molecular components of the periodontal ligament and alveolar bone (1). Conventional knowledge dictates that during OTM, the applied force leads to the creation of a compression side and a tension side. By convention, the tension side is marked by bone formation, whereas the compression side is marked by bone resorption. The applied orthodontic forces are transmitted via the stress of the extracellular matrix to the cells in alveolar bone and periodontal ligament. The transmission of this mechanical signal alters the vascularization of the periodontal ligament (2) and creates microcracks in bone (3). This leads to the local release of various molecules, such as neurotransmitters, cytokines, growth factors, colony-stimulating factors, and arachidonic acid metabolites (4). The liberation of pro-inflammatory, angiogenic, and osteogenic mediators triggers ligament and bone remodelling prior to OTM (5, 6). The first phase during remodelling comprises the differentiation and activation of the precursors of osteoblasts and osteoclasts and is mediated by the release of interleukins IL-1, IL-6, IL-10, and IL-33; tumour necrosis factor alpha or the interferon γ (7). During the second phase, the osteoclast differentiation is completed. It is controlled by the binding of receptor activator of nuclear factor kappa-B ligand (RANKL) on the receptor activator of nuclear factor-kappa (RANK) membrane receptor of pre-osteoclasts. In parallel, the decoy receptor to RANKL, called osteoprotegerin (OPG), shows decreased expression. In fact, compressive forces upregulate the RANKL/OPG ratio and promote osteoclast differentiation (8). During bone remodelling, osteoblasts play a key role by regulating both resorption and formation (9). Osteoblasts transduce mechanical stimuli into biochemical signals through multiple mechanoreceptors, such as integrins, pericellular sensors, focal adhesions, ion channels, cadherins, connexins, and lipidic rafts (10). These proteins activate intracellular signalling cascades such as the MAPK, PI3K/Akt, GTPases, Wnt, and calcium pathways to modulate osteoblastic differentiation by up- or downregulating gene expression (11). There is a need to better understand how all these mechanisms function. Osteoblasts also secrete bone matrix components before mineralization. The synthesized proteins reflect the stage of osteoblast differentiation, thus serving as phenotypic markers. Alkaline phosphatase (ALP) and collagen I (Col1) are commonly used as markers of phenotypes under mechanical stress (12). Col1 is a major structural protein of the bone matrix, whereas ALP is believed to reflect the first step of bone mineralization. In addition, the sequential development of bone depends on the transcription factor Runt-related transcription factor 2 (RunX2). Compressive forces have been shown to increase RunX2 expression (13). Col1 and ALP expression are observed at an early stage of cell differentiation (14). As the cells progress to the mineralization phase, peak levels of non-collagenous proteins are reached: osteopontin (OP) and osteocalcin (OCN) are expressed (15). To better understand the changes taking place in bone during OTM, it is important to elucidate the cellular and molecular interactions within the bone tissue. The complexity of in vivo studies makes analysing these communication channels difficult. As a result, models using cells isolated from alveolar bone have been established to simulate the application of orthodontic forces (16). Currently, cell culture is widely used in biomedical research to better elucidate the physiology, pathology, and therapy of many processes and diseases. Much of our understanding of biological mechanisms such as migration, differentiation, and force sensing has been acquired from two-dimensional (2D; also called monolayer) cell culture (17). Although conventional 2D cell culture has predominated research fields, recent advancements in tissue engineering have allowed a shift towards culture using three-dimensional (3D) structures (18). Regardless, the generic term ‘3D cell culture’ actually covers a profound number of models, all of which have the characteristic of culturing cells within a matrix comprising poly-l-lactide acid, poly-lactic-co-glycolic acid, collagen, agarose, etc. The majority of these constructs are in fact isolated cells ‘standing up’ in the scaffold (19). These authors suggest going further by integrating micro-tissues in the form of so-called spheroid cellular aggregates in an agarose matrix. The use of agarose matrices limits cellular attachment to the physiological extracellular matrix and therefore allows a strict evaluation of the role of mechanical compression without interference from the microenvironment. In the field of orthodontics, the originally described 2D compression model (20) is still considered the gold standard. It is a simple method for the application of a static, compressive, unidirectional force on a cell layer. Conversely, 3D models are used to answer questions such as how do cells detect forces (mechano-sensing), how do cells convert mechanical stress into molecular signals (mechano-transduction), and how do these molecular signals influence the specific cellular response or to a specific force. Despite these possibilities, this 3D culture model requires strong validation before being included as a reference in the field. The objective of this work was to validate a 3D cell culture weight-based model to elucidate the complex cascade of secondary events to the application of various intensities of orthodontic forces. Materials and methods Cell culture The human foetal osteoblastic cell line hFOB 1.19 (ATCC® CCL-11372™) was seeded at a density of 6 × 103 cells/cm2 in Dulbecco’s modified Eagle medium (Lonza, Belgium) mixed with Hams F12 (Lonza, Belgium) at a ratio of 1:1 and supplemented with 10% foetal bovine serum (Gibco, Life Technologies, UK), 2 mmol/l glutamine (Lonza, Belgium), and antibiotics (100 U/ml penicillin, 100 µg/ml streptomycin; Lonza, Belgium). The cells were cultured in a humidified atmosphere containing 5% CO2 at a temperature of 34°C, during which time they were conditionally immortalized with a gene encoding a temperature-sensitive mutant (tsA58) of SV40 large T antigen (21). Upon reaching confluence, cells were passaged by trypsin/ethylenediaminetetraacetic acid treatment (Lonza, Belgium). 2D model Cells were re-suspended in 0.8% low-melting agarose solution at a density of 4 × 105 cells/ml for gene assays or 2 × 105 cells/ml for protein assays. One millilitre of the solution was added into a 6-well transwell polyethylene terephthalate (PET) insert with a pore size 3 µm and a pore density of 8 × 105 pores/cm2 (Corning, Lowell, Massachusetts); the total gel thickness was 2 mm. The mixture was allowed to solidify for 20 minutes at room temperature. Then, 1.5 and 2.5 ml of medium were added to the upper and lower chambers of the well, respectively. Finally, a thin glass slide was placed on top of the agarose gel. Spheroids Flat-bottom 96-well microplates were treated with 50 µL of 1% agarose (Sigma-Aldrich, Saint-Louis, Missouri, USA) prepared in phosphate-buffered saline (Lonza, Belgium) to form a thin non-adhesive film to prevent cell adhesion (22). Cells were suspended at 1 × 104 cells/ml and aliquots of 200 µl i.e. 2000 cells/well were seeded into individual wells to initiate spheroid formation, and the microplates were incubated at 34°C in a 5% CO2 atmosphere. Cells were cultured in complete medium during spheroid growth that was not replaced. After 72 hours of spheroid formation, each aggregate was collected by aspiration with a 100 µl pipette. After centrifugation at 800 rpm for 3 minutes, 96 or 192 spheroids were re-suspended in 1 ml of 0.8% agarose low-melting solution type VII (Sigma-Aldrich, Saint Louis, Missouri, USA). Then, 1 ml of the solution was added into a 6-well transwell PET insert with a pore size of 3 µm and a pore density of 8 × 105 pores/cm2 in the membrane (Corning, Lowell, Massachusetts); the total gel thickness was 2 mm. The spheroid–agarose mixture was allowed to solidify for 20 minutes at room temperature. As previously described, 1.5 and 2.5 ml of medium were added to the upper and lower chambers of the well, respectively. Finally, a thin 22 mm diameter glass slide was placed on top of the spheroid–gel construct. Densities of 96 spheroids/ml for protein assays or 192 spheroids/ml for gene assays were used. Static compressive force In each of the two previous constructs, a static and unidirectional compressive force was applied. A plastic test tube cap was filled with lead pellets to adjust the weights. Based on the diameter of the glass slide that distributes the force over the agarose, the lead weight was adjusted to deliver forces of 1 or 4 g/cm2 (Figures 1 and 2). Figure 1. Open in new tabDownload slide Illustrations of the cell culture models for mechanical stress applications. Figure 1. Open in new tabDownload slide Illustrations of the cell culture models for mechanical stress applications. Figure 2. Open in new tabDownload slide Clinical photo of the setup. General overview (A) and focus on the weight resting on the agarose matrix placed at the bottom of the insert (B). Figure 2. Open in new tabDownload slide Clinical photo of the setup. General overview (A) and focus on the weight resting on the agarose matrix placed at the bottom of the insert (B). Cell survival and proliferation Spheroid diameters The increase in spheroid diameter was monitored by obtaining photographs via optical microscopy. The spheroids embedded in agarose were followed and compared with non-embedded spheroids constituting the control group. Collection of images at the standard ×4 magnification allowed the use of a measurement layer in Photoshop CC software. Spheroids appeared round or ovoid; in the latter case, the largest diameter was retained. Glucose consumption The growth and vitality of the hFOB 1.19 cells embedded in an agarose matrix were assessed by the indirect measurement of metabolic activity (23). Each day, glucose uptake of the cells was measured in the medium with a cobas c 111 analyser (Roche Diagnostics GmbH, Allemagne). Based on the final volume of the solution, the cumulative amount of glucose consumed was calculated. Fresh culture medium without cells was used as a reference. Each test was done in triplicate. Propidium iodide and Hoechst staining Staining was performed using 5 µg/ml Hoechst 33258 and 1 µg/ml propidium iodide. The ratio between viable and dead cells was established as the ratio of both channel intensities. For the 2D model, cells were incubated for 2 hours at 34°C and analysed by confocal microscopy. A total of five pictures of each condition was analysed. For the 3D model, cells were incubated for 20 hours incubation at 34°C with Hoechst, treated with propidium iodide 2 hours before reading, and analysed by light sheet fluorescence microscopy. For each condition, between three and six spheroids were imaged and analysed. RNA extraction and real-time polymerase chain reaction Total RNA was extracted from hFOB 1.19 osteoblasts at 24 and 72 hours using a Nucleospin RNA extraction kit (Macherey Nagel, Germany), with some adaptations of the manufacturer’s instructions with respect to the context of spheroids embedded in agarose matrix. The volume of reagent RA1 was multiplied by 2, i.e. increased to 700 µl. Quantitative real-time polymerase chain reaction was carried out with a SYBR Green PCR kit (Applied Biosystems, Foster City, California, USA) in an ABI Prism 7900 HT Fast system (Applied Biosystems) under the following cycling conditions: 2 minutes at 50°C; 10 minutes at 95°C; 40 cycles of 15 seconds at 95°C and 1 minutes at 60°C; and a final dissociation step. The primer sequences used in this experiment are listed in Table 1. The transcripts of the 18S and HPRT housekeeping gene were used normalize the mRNA levels (24, 25). The results were normalized using the geometric mean Ct of HPRT and 18S. Quantitative results were analysed with SDS v2.3 software (Applied Biosystems). Relative gene expression was calculated with the 2−ΔΔCt method (26). The assay for each gene was carried out in triplicate. The validation of the primers, the purity and integrity of the RNA, and the amplification efficiency have been verified according to the recommendations of Taylor and Mrkusich (27). Table 1. Description and characteristics of primers used for real-time PCR. Gene . Sequence 5′–3′ . Exon position . Product size (bp) . Primer efficiency, Ep (%) . Coefficient of determination, R2 . ALP F: AGAACCCCAAAGGCTTCTTC F: E7 74 98.8 0.999 R: CTTGGCTTTTCCTTCATGGT R: E8 Col1A1 F: CGAAGACATCCCACCAATCAC F: E1/2 98 97.2 0.999 R: TGTCGCAGACGCAGAT R: E2 RunX2 F: ACCCAGAAGGCACAGACAGAAG F: E4/5 82 97.2 0.999 R: AGGAATGCGCCCTAAATCACT R: E5 OCN F: GCAGCGAGGTAGTGAAGAGA F: E3 137 92.7 0.996 R: GATGTGGTCAGCCAACTCGT R: E4 OP F: TCACCTGTGCCATACCAGTTAAA F: E1/2 85 98.0 0.991 R: GCCACAGCATCTGGGTATTTG R: E3 18S F: ATTAAGGGTGTGGGCCGAAG F: E1/2 111 96.1 0.999 R: GGTGATCACACGTTCCACCT R: E2/3 HPRT F: AGCTTGCTGGTGAAAAGG F: E6/7 107 97.6 0.993 R: TCATTATAGTCAAGGGCATATC R: E8 Gene . Sequence 5′–3′ . Exon position . Product size (bp) . Primer efficiency, Ep (%) . Coefficient of determination, R2 . ALP F: AGAACCCCAAAGGCTTCTTC F: E7 74 98.8 0.999 R: CTTGGCTTTTCCTTCATGGT R: E8 Col1A1 F: CGAAGACATCCCACCAATCAC F: E1/2 98 97.2 0.999 R: TGTCGCAGACGCAGAT R: E2 RunX2 F: ACCCAGAAGGCACAGACAGAAG F: E4/5 82 97.2 0.999 R: AGGAATGCGCCCTAAATCACT R: E5 OCN F: GCAGCGAGGTAGTGAAGAGA F: E3 137 92.7 0.996 R: GATGTGGTCAGCCAACTCGT R: E4 OP F: TCACCTGTGCCATACCAGTTAAA F: E1/2 85 98.0 0.991 R: GCCACAGCATCTGGGTATTTG R: E3 18S F: ATTAAGGGTGTGGGCCGAAG F: E1/2 111 96.1 0.999 R: GGTGATCACACGTTCCACCT R: E2/3 HPRT F: AGCTTGCTGGTGAAAAGG F: E6/7 107 97.6 0.993 R: TCATTATAGTCAAGGGCATATC R: E8 Open in new tab Table 1. Description and characteristics of primers used for real-time PCR. Gene . Sequence 5′–3′ . Exon position . Product size (bp) . Primer efficiency, Ep (%) . Coefficient of determination, R2 . ALP F: AGAACCCCAAAGGCTTCTTC F: E7 74 98.8 0.999 R: CTTGGCTTTTCCTTCATGGT R: E8 Col1A1 F: CGAAGACATCCCACCAATCAC F: E1/2 98 97.2 0.999 R: TGTCGCAGACGCAGAT R: E2 RunX2 F: ACCCAGAAGGCACAGACAGAAG F: E4/5 82 97.2 0.999 R: AGGAATGCGCCCTAAATCACT R: E5 OCN F: GCAGCGAGGTAGTGAAGAGA F: E3 137 92.7 0.996 R: GATGTGGTCAGCCAACTCGT R: E4 OP F: TCACCTGTGCCATACCAGTTAAA F: E1/2 85 98.0 0.991 R: GCCACAGCATCTGGGTATTTG R: E3 18S F: ATTAAGGGTGTGGGCCGAAG F: E1/2 111 96.1 0.999 R: GGTGATCACACGTTCCACCT R: E2/3 HPRT F: AGCTTGCTGGTGAAAAGG F: E6/7 107 97.6 0.993 R: TCATTATAGTCAAGGGCATATC R: E8 Gene . Sequence 5′–3′ . Exon position . Product size (bp) . Primer efficiency, Ep (%) . Coefficient of determination, R2 . ALP F: AGAACCCCAAAGGCTTCTTC F: E7 74 98.8 0.999 R: CTTGGCTTTTCCTTCATGGT R: E8 Col1A1 F: CGAAGACATCCCACCAATCAC F: E1/2 98 97.2 0.999 R: TGTCGCAGACGCAGAT R: E2 RunX2 F: ACCCAGAAGGCACAGACAGAAG F: E4/5 82 97.2 0.999 R: AGGAATGCGCCCTAAATCACT R: E5 OCN F: GCAGCGAGGTAGTGAAGAGA F: E3 137 92.7 0.996 R: GATGTGGTCAGCCAACTCGT R: E4 OP F: TCACCTGTGCCATACCAGTTAAA F: E1/2 85 98.0 0.991 R: GCCACAGCATCTGGGTATTTG R: E3 18S F: ATTAAGGGTGTGGGCCGAAG F: E1/2 111 96.1 0.999 R: GGTGATCACACGTTCCACCT R: E2/3 HPRT F: AGCTTGCTGGTGAAAAGG F: E6/7 107 97.6 0.993 R: TCATTATAGTCAAGGGCATATC R: E8 Open in new tab Enzyme-linked immunosorbent assay To determine the IL-6, IL-8, and OPG concentrations produced by hFOB 1.19 cells in response to mechanical compression, supernatants from 2D and 3D cell cultures were harvested, and the concentrations were measured by enzyme-linked immunosorbent assay (DuoSet ELISA Development System Human IL-6, OPG, and CXCL8, R&D Systems, Minneapolis, Minnesota, USA) according to the manufacturer’s instructions. Supernatants were collected at 24 and 72 hours, and their absorbance at 450 nm was read with a spectrophotometer. Concentrations were calculated from the standard curves. Three independent experiments were performed for each condition. Statistical analysis Data analysis was performed with GraphPad Prism v6.0 (GraphPad Software, San Diego, California, USA). Two-way analysis of variance and Tukey’s post hoc tests were used to determine significant differences between groups. Data are expressed as the means ± standard deviation, and differences were considered significant when P < 0.05. Results Cell viability Cell viability was checked in both models using two approaches. First, a metabolic- (glucose consumption over time) and optic-based (spheroid diameters) method was used. Over a period of 4 days, glucose consumption was maintained in both the 2D and 3D cultures indicating cellular viability; furthermore, the spheroid diameters stabilized over time as there was no statistically significant difference between the culture times for each stress intensity (Figure 3). Figure 3. Open in new tabDownload slide Cumulative glucose uptake of hFOB 1.19 cells cultivated for 96 hours under various compressive weights in a 2D culture (A) and a 3D spheroid culture (B). Monitoring of spheroid diameters under compressive strain (C). Figure 3. Open in new tabDownload slide Cumulative glucose uptake of hFOB 1.19 cells cultivated for 96 hours under various compressive weights in a 2D culture (A) and a 3D spheroid culture (B). Monitoring of spheroid diameters under compressive strain (C). The second approach consisted of the measurement of Hoechst 33258 staining, which labels viable cells, and propidium iodide staining, which labels dead cells. The Hoechst/IP ratio in the 2D model demonstrated constant viability at 24 hours regardless of the compressive force exerted and a slightly reduced viability of 15% at 72 hours at the highest pressure applied. In the spheroid model, at 24 hours, the viability was significantly increased in response to the strongest compressive force and was not significantly decreased at 72 hours in response to a pressure of 1 g/cm2 (Figure 4) (Supplemental data). Figure 4. Open in new tabDownload slide (A) Hoechst/PI ratio based on immunohistochemical staining as recorded by confocal microscopy for 2D culture and (C) images of cells in the 2D model after 72 hours of compressive stress. (B) Light sheet fluorescence microscopy for 3D culture and (D) images of spheroids after 72 hours of compressive stress. Figure 4. Open in new tabDownload slide (A) Hoechst/PI ratio based on immunohistochemical staining as recorded by confocal microscopy for 2D culture and (C) images of cells in the 2D model after 72 hours of compressive stress. (B) Light sheet fluorescence microscopy for 3D culture and (D) images of spheroids after 72 hours of compressive stress. Expression of osteogenic markers In the 2D model, ALP gene expression dropped significantly between 24 and 72 hours (−50%) (Figure 5A), whereas in the 3D model, the opposite phenomenon was observed with a significant increase (+250%) in gene expression over time (Figure 5B). At 72 hours, ALP increased by a factor of 4–10 between the 2D and 3D models depending on the compressive force. The influence of the pressure exerted was more significant at 72 hours in both cases. For the 2D model, the increase was correlated with the applied pressure (+40% for 1 g/cm2 and +47% for 4 g/cm2), but this increase was different for spheroids in which ALP expression peaked at 1 g/cm2 (+44 and +20% for 24 and 72 hours, respectively) but was negligible at 4 g/cm2. Figure 5. Open in new tabDownload slide Quantitative reverse transcription-polymerase chain reaction demonstrating the effect of different magnitudes of compressive force for 24 or 72 hours on the mRNA expression of ALP (A and B), Col1 (C and D), OCN (E and F), OP (G and H), and Runx2 (I and J) relative to that of the endogenous control (HPRT). Numbers are presented as the mean and SEM (*P < 0.05 compared with 0 g/cm2, #P < 0.05 compared with 1 g/cm2). Figure 5. Open in new tabDownload slide Quantitative reverse transcription-polymerase chain reaction demonstrating the effect of different magnitudes of compressive force for 24 or 72 hours on the mRNA expression of ALP (A and B), Col1 (C and D), OCN (E and F), OP (G and H), and Runx2 (I and J) relative to that of the endogenous control (HPRT). Numbers are presented as the mean and SEM (*P < 0.05 compared with 0 g/cm2, #P < 0.05 compared with 1 g/cm2). In the 2D model, at 24 hours, COL1 expression gradually decreased with increased pressure exerted, whereas at 72 hours, only the 4 g/cm2 pressure caused a decrease of nearly 55% (Figure 5C). In the 3D model, the expression remained constant regardless of the time under pressure or the pressure exerted (Figure 5D). OCN expression was almost constant at 24 and 72 hours in the 2D model (Figure 5E). In contrast, it was significantly increased in the 3D model at 24 and 72 hours for the strongest pressure (≥80% compared with the control in both cases) (Figure 5F). OP expression decreased with pressure and time in the 2D model (Figure 5G), with the largest decrease observed after the application of the largest force for 3 days (greater than 50% decrease). This reduction in expression depending of the pressure was less marked in the 3D model since the largest decrease recorded was only 23% at 24 hours and no difference was detected at 72 hours (Figure 5H). In the 2D model, RunX2 expression was significantly modified in response to the strongest pressure applied (greater than 150% compared with the control) (Figure 5I). In the spheroid model, RunX2 expression was significantly increased with the applied pressure at 24 hours (+36 and +47%, respectively) and was significantly increased for the highest pressure at 72 hours (+50%) (Figure 5J). Quantification of soluble factors The release of IL-6 in the 2D model was low, with a peak value of 500 ng/ml, while in the 3D model, the release was 40–80 times higher. In the 2D model, the measured quantities were independent of applied pressure at 24 hours and had doubled at 72 hours in a pressure-dependent manner. In the 3D cultures, the IL-6 concentrations were very high and were proportional to the pressure applied twofold at 24 hours, with larger quantities released at 72 hours (+200%) (Figure 6A and 6B). Figure 6. Open in new tabDownload slide Effect of varying magnitudes of compressive force on the protein expression of IL-6 (A and B), IL-8 (C and D), and OPG (E and F) after 24 or 72 hours of compressive force. The experiments were performed in triplicate. Numbers are presented as the mean and SEM. IL, interleukin; OPG, osteoprotegerin (*P < 0.05 compared with 0 g/cm2). Figure 6. Open in new tabDownload slide Effect of varying magnitudes of compressive force on the protein expression of IL-6 (A and B), IL-8 (C and D), and OPG (E and F) after 24 or 72 hours of compressive force. The experiments were performed in triplicate. Numbers are presented as the mean and SEM. IL, interleukin; OPG, osteoprotegerin (*P < 0.05 compared with 0 g/cm2). Regarding the release of IL-8 in the 2D model, the measured quantities were low (approximately 400 ng/ml). The only significant value was at 72 hours with a pressure of 4 g/cm2, where the amount of IL-8 was 10 times greater than that in the control. For spheroids, the values were significantly higher than those in the 2D cultures (from 11- to 40-fold). Unlike that observed with IL-6, there was no pressure-dependent effect, as no significant difference was highlighted. There was a time-dependent effect since the amount of IL-8 increased 250% at 72 hours compared with the levels at 24 hours (Figure 6C and 6D). For OPG, the difference was not marked as previously described in the 2D model and did not exceed a factor of 5 for each condition tested with amounts ranging from 133 to 1042 ng/ml. The quantity fluctuated little in response to pressure. Similarly, for the 3D culture, no significant difference was observed in response to the forces applied for 24 and 72 hours. The amount of protein was greater (+700%) after 3 days of pressure (Figure 6E and 6F). Discussion Many cell lines are available to create culture models. However, we wanted to use a human line, closer to the clinical context. The hFOB 1.19, a human osteoblastic cell line was immortalized by transfected temperature-sensitive SV40 T-antigen. This makes it possible to free of the inter-individual variability of the primo-cultures. This line is closer to the healthy cell than the tumoural lines and allows suitable assays and comparison to others studies. In addition, the selected line allows the formation of spheroids, their embedding in agarose and mechanical stress application. We have control that the aggregates did not dissociate during the mechanical stressing. The hFOB 1.19 line appeared to be the most appropriate since its spheroids have a high rate of internal cohesion. Mechanical stress applied here was from 1 to 4 g/cm2. This is in the range of the most commonly applied forces in orthodontic weight-loaded approach models (16). Micro-tissue constructs such as spheroids raise challenging questions such as control of the spatio-temporal distribution of nutrients and oxygen, the primary basis for cell viability and proliferation (18). Cell aggregates formed in a 3D culture create a gradient of nutrient access. The first validation step of the model was to control the viability of cells embedded in the agarose matrix. Conventional colorimetric methods based on the reduction of resazurin (Alamar blue assay) or tetrazolium ([3-{4,5-dimethylthiazol-2-yl}-2,5-diphenyltetrazolium bromide] tetrazolium assay) are indeed frequently used to evaluate viability in 2D cell culture. However, they are not applicable to spheroids. Tight cell–cell junctions can affect the absorption and diffusion kinetics of the dyes, thus modifying the reading of the test and making the results more difficult to interpret (28, 29). To avoid the problem, we evaluated viability by monitoring diameters of the spheroids and indirectly by measuring the glucose uptake from the medium. In one hand, the diameter of the spheroids was not changed for the lowest pressures. On the other hand, for the most important pressure, the diameter increased strongly during the first hours then stabilized. This is explained by the loss of the spherical aspect of the aggregates in this condition in favour of an ogival aspect whose large diameter we have measured. At the end of the 72 hours compression period, glucose was still metabolized in both the 2D and 3D models. This result leads to the conclusion that hFOB 1.19 cells are still proliferating at this time regardless of the model used. This behaviour also indicates that the pore sizes of the insert membrane and matrix mesh are suitable for the exchange of nutrients and waste products between the two compartments. The surface of the spheroid aggregates exhibits the highest levels of proliferation, whereas the interior possesses the highest number of quiescent or necrotic cells (30). These data were also supported on our micro-tissues by light sheet fluorescence microscopy after labelling cells that had lost their membrane integrity with propidium iodide and living cells with Hoechst. Globally, glucose consumption and the use of fluorochromes confirmed preserved viability in spheroids compared with that of isolated cells embedded in agarose. These results are in agreement with those of previous studies (23, 31). A major limitation of culture in a matrix is the limited cell survival to 3 or 4 days. Only the poly-l-lactide matrix allows a longer follow-up, although limited to 2 weeks (32). This is why we did not carry out a follow-up on a longer term. In addition, our work focuses on the induction of tooth movement which take place during the first hours of application of orthodontic forces (3). As several studies demonstrate their variation according to the mechanical stress imposed on cells, markers of osteoblastic activity were analysed by qRT-PCR (16). We used two reference genes according to the work of Vazquez and Rauh: 18S and HPRT (24, 25). This is a different choice from Kirschneck’s study which compares 13 reference genes for the study of the consequences of dental movement on periodontal ligament fibroblasts (33). This is justified by a different culture model (3D versus 2D) and different cell populations (osteoblasts versus fibroblasts). The comparison of the gene expression between the 2D and 3D cell culture models revealed three situations: 1) There was no variation in the gene expression profile between the two models. 2) The expression of the markers was increased in the 3D model. 3) The expression was decreased in the 2D model but was maintained in the 3D model. For the first observation, there was no difference in the expression profile between the two models with regard to RunX2, a master gene that transcriptionally promotes osteoblast differentiation. Its expression is essential for the homeostasis of bone tissue. This transcription factor has previously been shown to regulate the expression of various osteoblastic genes (34). In both model systems, its expression was increased in relation to the intensity and the duration of the pressure exerted. This confirms earlier works, such as Shen’s study, which used murine osteoblasts embedded in a collagen matrix and reported that compressive stress ranging from 0 to 5 g/cm2 enhances osteoblastic differentiation (35). In the case of ALP and OCN, the gene expression of both products was greatly increased in the spheroids compared with the 2D model. ALP is a very early marker of osteoblast differentiation and contributes to bone mineralization. Its expression was almost constant in the 2D cultures (no significant variations as a function of time or compressive force). In the 3D model, there was a difference between 24 and 72 hours of incubation with a maximum (×5 compared with the 24 hours control) value observed at 1 g/cm2. These results are consistent with those of the Elashry study, which demonstrated that 10 days of mechano-stimulation were needed to elicit a significant increase in ALP expression (30). Other studies confirm the increase in the expression of ALP in the spheroid model (36). This increase is also observed in saliva during orthodontic treatment (37). This expression is due to both mechanical and genetic factors: the interaction of osteoblasts in 3D space and tension within the spheroids. It is well known that the OCN gene, as a target of β-catenin, plays a decisive role in osteogenesis. Therefore, the expression of genes involved in the Wnt/β-catenin pathway contributes to spheroid formation (38, 39). The expression of the OCN gene is characteristic of the late stages of osteoblastic differentiation. OCN, the most abundant non-collagenous protein in bone, can be considered a signal that recruits fuses osteoclastic progenitors until they mature (40). Hoang demonstrated the ability of this protein to bind hydroxyapatite crystals and to inhibit their growth (41). On the one hand, its expression is almost constant in the 2D model regardless of the time and the compressive force applied. On the other hand, its expression was clearly increased with time in the 3D model and peaked in response to the highest compressive force. This is in agreement with results in the literature reporting that OCN is increased during orthodontic treatment in vivo (7) as well as in vitro (42). In the last experiment, Col1 and OP gene expression levels decreased as a function of time and pressure applied in the 2D model but were constant in the 3D model. Col1 is a marker of the later stages of osteoblastic differentiation. Its expression decreased proportionally to the weight applied in the 2D model at 24 hours and remained low at 72 hours. Regarding spheroids, neither time nor force applied had an impact on Col1 expression. Many studies have demonstrated the induction of collagen expression by mechanical forces (12, 43). However, Min reported that the mechano-response of fibroblasts varies according to the force applied: gene induction with weak compression and deactivation with stronger forces (44). In spheroids, the abundance of intercellular junctions allows the distribution of the force exerted to all the cells, unlike the cells isolated in agarose matrix (39). This could explain the decrease in collagen expression in the 2D model and the consistent gene expression in the 3D model. OP gene expression was significantly diminished as a function of time and pressure exerted in the 2D model, while its expression was almost constant in the 3D model. As described by Singh, OP is a molecule released by osteoblasts and acts as a signal for the recruitment, differentiation, and attachment of osteoclasts (45). This phenomenon allows for the degradation of the non-mineralized matrix and characterizes the later stages of osteoblastic differentiation. OP has been shown to have a regulatory effect on hydroxyapatite crystal growth and potently inhibits the mineralization of osteoblast cultures in a phosphate-dependent manner. Bone remodelling is vital for OTM (46); thus, it is crucial that OP expression be maintained regardless of the forces applied and the time of application (45). Some authors showed the increase in its concentration on the compression side from the initiation of the OTM in vivo with a peak reached on the second day after initiating pressure. Following the increase in the expression of osteopontin, the number of osteoclasts increased 17-fold compared with the control, and numerous resorption depressions were observed on the pressure side of the alveolar bone (47). Regarding the release of IL-6 and IL-8, in this study, a temporal effect between 24 and 72 hours was observed in both models. These cytokines are characteristic of the inflammatory response and osteoblastic activity. Moreover, these markers of osteoblast activity are usually found in the crevicular fluid of teeth subjected to OTM (4). IL-6, which is produced by various cells in the medullary microenvironment (including monocytes, fibroblasts, and T cells), is an important regulator of osteoclast formation in vivo (48). It is also produced by osteoblasts and induces bone resorption both alone and in concert with other bone resorption agents (49). Although one study found that compressive forces have little impact on alveolar bone cells in 2D culture models (50), other studies demonstrated that static compressive forces promote IL-6 release in similar models (51, 52). In our study, the quantities of IL-6 collected in the 3D model were 40 times higher than those collected in the 2D model. The increase in IL-6 release in the 3D model could be explained by a model effect: spheroidal culture leads to increased transcription of ILs (48), and IL-6 expression is greatly enhanced by mechanical stimulation (53). IL-8 stimulates both osteoclastogenesis and bone resorption in human osteoclasts derived from peripheral blood mononuclear cells. It also regulates the expression of the essential osteoclastogenic factor RANKL (54). Similar to IL-6, IL-8 was released in greater quantities in the 3D model. In the 2D model, its release was very low at 24 hours but had increased at 72 hours. These results agree with those of another study showing that IL-8 expression was only modulated by the application of compressive forces ranging from 0.5 to 3 g/cm2 during the first 24 hours (55). Finally, it has been shown that OPG plays a key role in RANK/RANKL coupling. OPG is a soluble cytokine that inhibits osteoclastogenesis via its endogenous antagonist activity that neutralizes the biological effects of RANKL and therefore acts as an inhibitor of bone resorption. The biological effects of OPG on bone cells include inhibition of the terminal stages of osteoclast differentiation, suppression of mature osteoclast activation, and induction of apoptosis. It is obvious that the RANKL/OPG interaction plays a crucial role during tooth movement. An elevated RANKL/OPG ratio in the gingival fluid of displaced teeth has been demonstrated (8). In this study, OPG expression in the 2D model was almost non-existent at 24 hours and was weak at 72 hours. In the 3D model, the released quantities at 72 hours were seven times higher than those at 24 hours. The results showed that the quantities of OPG released were only weakly dependent on the compressive forces applied in the two models and that in both models, the strongest pressure induced a decrease in OPG release. These results are consistent with those conducted on bone cell lines (56). However, other authors have shown that OPG expression is clearly upregulated at both the mRNA and protein levels when cells were exposed to mechanical stimuli (57). However, it should be noted that these studies used mesenchymal stem cells and not osteoblasts. Conclusion In this study, a 3D model was developed to measure the impact of mechanical stress on osteoblasts. The results showed that cell viability and metabolism were conserved; however, differences in responses to mechanical stimuli were observed between 2D and 3D cultures. Indeed, there is a proportionality in the responses in the spheroids that is a time- and pressure-dependent in terms of the expression of osteogenesis markers and the release of cytokines. Spheroids are comparable to micro-tissues, and the results obtained in this study are more reminiscent to in vivo physiology than 2D models. It is therefore a seductive and more realistic approach to studying bone mechanobiology. Supplementary material Supplementary data are available at European Journal of Orthodontics online. Supplemental data. Exploration by light sheet microscopy of a spheroid grown in an agarose matrix and compressed by a force of 4 g/cm2 for 72 hours. The staining performed here is Hoechst 33258 (blue) and propidium iodide (red). Acknowledgements The authors acknowledge the members of the Microscopy Rennes Imaging Center (BIOSIT, Université de Rennes 1), in particular, Drs X. Pinson and S. Dutertre. They also acknowledge the Biochemistry-Toxicology Laboratory of the University Hospital of Rennes, Pr C. Bendavid and Dr M. Bourgeois. Funding This work was supported by the French Society of Orthodontics and Dento Facial Orthopedics (SFODF); the National Institute for Medical Research (INSERM) and the French Federation of Orthodontics (FFO). Ethical approval human and animal experiments Human commercial cells were used according to the Laboratory accreditation DC-2012-1573. Conflicts of interest None to declare. References 1. Isola , G. , Matarese , G., Cordasco , G., Perillo , L. and Ramaglia , L. ( 2016 ) Mechanobiology of the tooth movement during the orthodontic treatment: a literature review . Minerva Stomatologica , 65, 299–327. OpenURL Placeholder Text WorldCat 2. Schwarz , A.M . ( 1932 ) Tissue changes incidental to orthodontic tooth movement . International Journal of Orthodontia, Oral Surgery and Radiography , 18 , 331 – 352 . Google Scholar Crossref Search ADS WorldCat 3. Henneman , S. , Von den Hoff , J.W. and Maltha , J.C. ( 2008 ) Mechanobiology of tooth movement . European Journal of Orthodontics , 30 , 299 – 306 . Google Scholar Crossref Search ADS PubMed WorldCat 4. Krishnan , V. and Davidovitch , Z. ( 2006 ) Cellular, molecular, and tissue-level reactions to orthodontic force . American Journal of Orthodontics and Dentofacial Orthopedics , 129 , 469.e1 – 469.32 . Google Scholar Crossref Search ADS WorldCat 5. Feller , L. , Khammissa , R.A., Schechter , I., Thomadakis , G., Fourie , J. and Lemmer , J. ( 2015 ) Biological events in periodontal ligament and alveolar bone associated with application of orthodontic forces . The Scientific World Journal , 2015 , 876509 . Google Scholar PubMed OpenURL Placeholder Text WorldCat 6. Masella , R.S. and Meister , M. ( 2006 ) Current concepts in the biology of orthodontic tooth movement . American Journal of Orthodontics and Dentofacial Orthopedics , 129 , 458 – 468 . Google Scholar Crossref Search ADS PubMed WorldCat 7. Vansant , L. , Cadenas De Llano-Pérula , M., Verdonck , A. and Willems , G. ( 2018 ) Expression of biological mediators during orthodontic tooth movement: a systematic review . Archives of Oral Biology , 95 , 170 – 186 . Google Scholar Crossref Search ADS PubMed WorldCat 8. Yamaguchi , M . ( 2009 ) RANK/RANKL/OPG during orthodontic tooth movement . Orthodontics & Craniofacial Research , 12 , 113 – 119 . Google Scholar Crossref Search ADS PubMed WorldCat 9. Tripuwabhrut , P. , Mustafa , M., Gjerde , C.G., Brudvik , P. and Mustafa , K. ( 2013 ) Effect of compressive force on human osteoblast-like cells and bone remodelling: an in vitro study . Archives of Oral Biology , 58 , 826 – 836 . Google Scholar Crossref Search ADS PubMed WorldCat 10. Brézulier , D. , Pellen-Mussi , P., Sorel , O. and Jeanne , S. ( 2018 ) La mécanobiologie osseuse, un domaine émergeant : revue de littérature . Orthodontie Française , 89 , 343 – 353 . Google Scholar Crossref Search ADS PubMed WorldCat 11. Miyashita , S. , Ahmed , N.E., Murakami , M., Iohara , K., Yamamoto , T., Horibe , H., Kurita , K., Takano-Yamamoto , T. and Nakashima , M. ( 2017 ) Mechanical forces induce odontoblastic differentiation of mesenchymal stem cells on three-dimensional biomimetic scaffolds . Journal of Tissue Engineering and Regenerative Medicine , 11 , 434 – 446 . Google Scholar Crossref Search ADS PubMed WorldCat 12. Pavlin , D. , Dove , S.B., Zadro , R. and Gluhak-Heinrich , J. ( 2000 ) Mechanical loading stimulates differentiation of periodontal osteoblasts in a mouse osteoinduction model: effect on type I collagen and alkaline phosphatase genes . Calcified Tissue International , 67 , 163 – 172 . Google Scholar Crossref Search ADS PubMed WorldCat 13. Baumert , U. , Golan , I., Becker , B., Hrala , B.P., Redlich , M., Roos , H.A., Palmon , A., Reichenberg , E. and Müssig , D. ( 2004 ) Pressure simulation of orthodontic force in osteoblasts: a pilot study . Orthodontics & Craniofacial Research , 7 , 3 – 9 . Google Scholar Crossref Search ADS PubMed WorldCat 14. Wlodarski , K.H. and Reddi , A.H. ( 1986 ) Alkaline phosphatase as a marker of osteoinductive cells . Calcified Tissue International , 39 , 382 – 385 . Google Scholar Crossref Search ADS PubMed WorldCat 15. Lian , J.B. and Stein , G.S. ( 1995 ) Development of the osteoblast phenotype: molecular mechanisms mediating osteoblast growth and differentiation . The Iowa Orthopaedic Journal , 15 , 118 – 140 . Google Scholar PubMed OpenURL Placeholder Text WorldCat 16. Janjic , M. , Docheva , D., Trickovic Janjic , O., Wichelhaus , A. and Baumert , U. ( 2018 ) In vitro weight-loaded cell models for understanding mechanodependent molecular pathways involved in orthodontic tooth movement: a systematic review . Stem Cells International . https://www.hindawi.com/journals/sci/2018/3208285/ (cited 30 October 2018). OpenURL Placeholder Text WorldCat 17. Baker , B.M. and Chen , C.S. ( 2012 ) Deconstructing the third dimension: how 3D culture microenvironments alter cellular cues . Journal of Cell Science , 125 , 3015 – 3024 . Google Scholar Crossref Search ADS PubMed WorldCat 18. Duval , K. , Grover , H., Han , L.H., Mou , Y., Pegoraro , A.F., Fredberg , J. and Chen , Z. ( 2017 ) Modeling physiological events in 2D vs. 3D cell culture . Physiology (Bethesda, Md.) , 32 , 266 – 277 . Google Scholar PubMed OpenURL Placeholder Text WorldCat 19. Yang , L. , Yang , Y., Wang , S., Li , Y. and Zhao , Z. ( 2015 ) In vitro mechanical loading models for periodontal ligament cells: from two-dimensional to three-dimensional models . Archives of Oral Biology , 60 , 416 – 424 . Google Scholar Crossref Search ADS PubMed WorldCat 20. Kanai , K . ( 1992 ) Initial effects of continuously applied compressive stress to human periodontal ligament fibroblasts . Journal of Japan Orthodontic Society , 51 , 153 – 163 . OpenURL Placeholder Text WorldCat 21. Harris , S.A. , Enger , R.J., Riggs , B.L. and Spelsberg , T.C. ( 1995 ) Development and characterization of a conditionally immortalized human fetal osteoblastic cell line . Journal of Bone and Mineral Research , 10 , 178 – 186 . Google Scholar Crossref Search ADS PubMed WorldCat 22. Friedrich , J. , Seidel , C., Ebner , R. and Kunz-Schughart , L.A. ( 2009 ) Spheroid-based drug screen: considerations and practical approach . Nature Protocols , 4 , 309 – 324 . Google Scholar Crossref Search ADS PubMed WorldCat 23. Detsch , R. , Alles , S., Hum , J., Westenberger , P., Sieker , F., Heusinger , D., Kasper , C. and Boccaccini , A.R. ( 2015 ) Osteogenic differentiation of umbilical cord and adipose derived stem cells onto highly porous 45S5 Bioglass®-based scaffolds . Journal of Biomedical Materials Research. Part A , 103 , 1029 – 1037 . Google Scholar Crossref Search ADS PubMed WorldCat 24. Rauh , J. , Jacobi , A. and Stiehler , M. ( 2015 ) Identification of stable reference genes for gene expression analysis of three-dimensional cultivated human bone marrow-derived mesenchymal stromal cells for bone tissue engineering . Tissue Engineering. Part C, Methods , 21 , 192 – 206 . Google Scholar Crossref Search ADS PubMed WorldCat 25. Vazquez , M. , Evans , B.A.J., Riccardi , D., Evans , S.L, Ralphs , J.R., Dillingham , C.M., et al. ( 2014 ) A new method to investigate how mechanical loading of osteocytes controls osteoblasts . Frontiers in Endocrinology , 5 . https://www.frontiersin.org/articles/10.3389/fendo.2014.00208/full (cited 19 December 2019). OpenURL Placeholder Text WorldCat 26. Livak , K.J. and Schmittgen , T.D. ( 2001 ) Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method . Methods (San Diego, California) , 25 , 402 – 408 . Google Scholar Crossref Search ADS WorldCat 27. Taylor , S.C. and Mrkusich , E.M. ( 2014 ) The state of RT-quantitative PCR: firsthand observations of implementation of minimum information for the publication of quantitative real-time PCR experiments (MIQE) . Journal of Molecular Microbiology and Biotechnology , 24 , 46 – 52 . Google Scholar Crossref Search ADS PubMed WorldCat 28. Kijanska, M., Kelm, J. (2016) In vitro 3D Spheroids and microtissues: ATP-based cell viability and toxicity assays. In: Sittampalam, G. S., Coussens, N. P., Brimacombe, K., et al. (eds.), Assay Guidance Manual [Internet]. Eli Lilly & Company and the National Center for Advancing Translational Sciences, Bethesda, MD. 29. Walzl , A. , Unger , C., Kramer , N., Unterleuthner , D., Scherzer , M., Hengstschläger , M., Schwanzer-Pfeiffer , D. and Dolznig , H. ( 2014 ) The resazurin reduction assay can distinguish cytotoxic from cytostatic compounds in spheroid screening assays . Journal of Biomolecular Screening , 19 , 1047 – 1059 . Google Scholar Crossref Search ADS PubMed WorldCat 30. Edmondson , R. , Broglie , J.J., Adcock , A.F. and Yang , L. ( 2014 ) Three-dimensional cell culture systems and their applications in drug discovery and cell-based biosensors . Assay and Drug Development Technologies , 12 , 207 – 218 . Google Scholar Crossref Search ADS PubMed WorldCat 31. Cheng , G. , Tse , J., Jain , R.K. and Munn , L.L. ( 2009 ) Micro-environmental mechanical stress controls tumor spheroid size and morphology by suppressing proliferation and inducing apoptosis in cancer cells . PLoS ONE , 4 , e4632 . Google Scholar Crossref Search ADS PubMed WorldCat 32. Liao , W. , Okada , M., Sakamoto , F., Okita , N., Inami , K., Nishiura , A., Hashimoto , Y. and Matsumoto , N. ( 2013 ) In vitro human periodontal ligament-like tissue formation with porous poly-l-lactide matrix . Materials Science & Engineering. C, Materials for Biological Applications , 33 , 3273 – 3280 . Google Scholar Crossref Search ADS PubMed WorldCat 33. Kirschneck , C. , Batschkus , S., Proff , P., Köstler , J., Spanier , G. and Schröder , A. ( 2017 ) Valid gene expression normalization by RT-qPCR in studies on hPDL fibroblasts with focus on orthodontic tooth movement and periodontitis . Scientific Reports , 7 , 14751 . Google Scholar Crossref Search ADS PubMed WorldCat 34. Ducy , P. , Zhang , R., Geoffroy , V., Ridall , A.L. and Karsenty , G. ( 1997 ) Osf2/Cbfa1: a transcriptional activator of osteoblast differentiation . Cell , 89 , 747 – 754 . Google Scholar Crossref Search ADS PubMed WorldCat 35. Shen , X.Q. , Geng , Y.M., Liu , P., Huang , X.Y., Li , S.Y., Liu , C.D., Zhou , Z. and Xu , P.P. ( 2017 ) Magnitude-dependent response of osteoblasts regulated by compressive stress . Scientific Reports , 7 , 44925 . Google Scholar Crossref Search ADS PubMed WorldCat 36. Moritani , Y. , et al. ( 2018 ) Spheroid culture enhances osteogenic potential of periodontal ligament mesenchymal stem cells . Journal of Periodontal Research , 53 , 870 – 882 . Google Scholar Crossref Search ADS PubMed WorldCat 37. Al-Khatieeb , M.M. , Rafeeq , R.A. and Saleem , A.I. ( 2018 ) Relationship between orthodontic force applied by monoblock and salivary levels of alkaline phosphatase and lactate dehydrogenase enzymes . The Journal of Contemporary Dental Practice , 19 , 1346 – 1351 . Google Scholar Crossref Search ADS PubMed WorldCat 38. Vincan , E , Schwab , R.H.M., Flanagan , D.J., Moselen , J.M., Tran , B.M., Barker , N., et al. ( 2018 ) Chapter ten—The central role of wnt signaling and organoid technology in personalizing anticancer therapy . In: Larraín , J., Olivares , G. (eds.), Progress in Molecular Biology and Translational Science. WNT Signaling in Health and Disease . Academic Press , vol. 153 , pp. 299 – 319 . http://www.sciencedirect.com/science/article/pii/S1877117317301874 (cited 4 June 2019). Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC 39. Marie , P.J. and Haÿ , E. ( 2013 ) Cadherins and Wnt signalling: a functional link controlling bone formation . BoneKEy Reports , 2 . https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3722765/ (cited 14 January 2019). OpenURL Placeholder Text WorldCat 40. Neve , A. , Corrado , A. and Cantatore , F.P. ( 2013 ) Osteocalcin: skeletal and extra-skeletal effects . Journal of Cellular Physiology , 228 , 1149 – 1153 . Google Scholar Crossref Search ADS PubMed WorldCat 41. Hoang , Q.Q. , Sicheri , F., Howard , A.J. and Yang , D.S. ( 2003 ) Bone recognition mechanism of porcine osteocalcin from crystal structure . Nature , 425 , 977 – 980 . Google Scholar Crossref Search ADS PubMed WorldCat 42. Shunzhi , Y. , Zhonghai , L. and Ning , Y. ( 2017 ) Mechanical stress affects the osteogenic differentiation of human ligamentum flavum cells via the BMP‑Smad1 signaling pathway . Molecular Medicine Reports , 16 , 7692 – 7698 . Google Scholar Crossref Search ADS PubMed WorldCat 43. Redlich , M. , Palmon , A., Zaks , B., Geremi , E., Rayzman , S. and Shoshan , S. ( 1998 ) The effect of centrifugal force on the transcription levels of collagen type I and collagenase in cultured canine gingival fibroblasts . Archives of Oral Biology , 43 , 313 – 316 . Google Scholar Crossref Search ADS PubMed WorldCat 44. Min , J. , Li , B., Liu , C., Guo , W., Hong , S., Tang , J. and Hong , L. ( 2017 ) Extracellular matrix metabolism disorder induced by mechanical strain on human parametrial ligament fibroblasts . Molecular Medicine Reports , 15 , 3278 – 3284 . Google Scholar Crossref Search ADS PubMed WorldCat 45. Singh , A. , Gill , G., Kaur , H., Amhmed , M. and Jakhu , H. ( 2018 ) Role of osteopontin in bone remodeling and orthodontic tooth movement: a review . Progress in Orthodontics , 19 , 18 . Google Scholar Crossref Search ADS PubMed WorldCat 46. Nomura , S. and Takano-Yamamoto , T. ( 2000 ) Molecular events caused by mechanical stress in bone . Matrix Biology , 19 , 91 – 96 . Google Scholar Crossref Search ADS PubMed WorldCat 47. Terai , K. , Takano-Yamamoto , T., Ohba , Y., Hiura , K., Sugimoto , M., Sato , M., Kawahata , H., Inaguma , N., Kitamura , Y. and Nomura , S. ( 1999 ) Role of osteopontin in bone remodeling caused by mechanical stress . Journal of Bone and Mineral Research , 14 , 839 – 849 . Google Scholar Crossref Search ADS PubMed WorldCat 48. Kurihara , N. , Bertolini , D., Suda , T., Akiyama , Y. and Roodman , G.D. ( 1990 ) IL-6 stimulates osteoclast-like multinucleated cell formation in long term human marrow cultures by inducing IL-1 release . Journal of Immunology (Baltimore, Md.: 1950) , 144 , 4226 – 4230 . Google Scholar PubMed OpenURL Placeholder Text WorldCat 49. Ishimi , Y. , Miyaura , C., Jin , C.H., Akatsu , T., Abe , E., Nakamura , Y., Yamaguchi , A., Yoshiki , S., Matsuda , T. and Hirano , T. ( 1990 ) IL-6 is produced by osteoblasts and induces bone resorption . Journal of Immunology (Baltimore, Md.: 1950) , 145 , 3297 – 3303 . Google Scholar PubMed OpenURL Placeholder Text WorldCat 50. Shi , J. , Baumert , U., Folwaczny , M. and Wichelhaus , A. ( 2018 ) Influence of static forces on the expression of selected parameters of inflammation in periodontal ligament cells and alveolar bone cells in a co-culture in vitro model . Clinical Oral Investigations . doi:10.1007/s00784-018-2697-2 (cited 27 October 2018). OpenURL Placeholder Text WorldCat 51. Tripuwabhrut , P. , Mustafa , K., Brudvik , P. and Mustafa , M. ( 2012 ) Initial responses of osteoblasts derived from human alveolar bone to various compressive forces . European Journal of Oral Sciences , 120 , 311 – 318 . Google Scholar Crossref Search ADS PubMed WorldCat 52. Asano , M. , Yamaguchi , M., Nakajima , R., Fujita , S., Utsunomiya , T., Yamamoto , H. and Kasai , K. ( 2011 ) IL-8 and MCP-1 induced by excessive orthodontic force mediates odontoclastogenesis in periodontal tissues . Oral Diseases , 17 , 489 – 498 . Google Scholar Crossref Search ADS PubMed WorldCat 53. Wang , X. , Fan , J., Zhang , M., Sun , Z. and Xu , G. ( 2013 ) Gene expression changes under cyclic mechanical stretching in rat retinal glial (Müller) cells . PloS One , 8 , e63467 . Google Scholar Crossref Search ADS PubMed WorldCat 54. Bendre , M.S. , Montague , D.C., Peery , T., Akel , N.S., Gaddy , D. and Suva , L.J. ( 2003 ) Interleukin-8 stimulation of osteoclastogenesis and bone resorption is a mechanism for the increased osteolysis of metastatic bone disease . Bone , 33 , 28 – 37 . Google Scholar Crossref Search ADS PubMed WorldCat 55. Koyama , Y. , Mitsui , N., Suzuki , N., Yanagisawa , M., Sanuki , R., Isokawa , K., Shimizu , N. and Maeno , M. ( 2008 ) Effect of compressive force on the expression of inflammatory cytokines and their receptors in osteoblastic Saos-2 cells . Archives of Oral Biology , 53 , 488 – 496 . Google Scholar Crossref Search ADS PubMed WorldCat 56. Grimm , S. , Walter , C., Pabst , A., Goldschmitt , J., Wehrbein , H. and Jacobs , C. ( 2015 ) Effect of compressive loading and incubation with clodronate on the RANKL/OPG system of human osteoblasts . Journal of Orofacial Orthopedics , 76 , 531 – 542 . Google Scholar Crossref Search ADS PubMed WorldCat 57. Zhang , L. , Liu , W., Zhao , J., Ma , X., Shen , L., Zhang , Y., Jin , F. and Jin , Y. ( 2016 ) Mechanical stress regulates osteogenic differentiation and RANKL/OPG ratio in periodontal ligament stem cells by the Wnt/β-catenin pathway . Biochimica et Biophysica Acta , 1860 , 2211 – 2219 . Google Scholar Crossref Search ADS PubMed WorldCat © The Author(s) 2020. Published by Oxford University Press on behalf of the European Orthodontic Society. All rights reserved. For permissions, please email: journals.permissions@oup.com This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic.oup.com/journals/pages/open_access/funder_policies/chorus/standard_publication_model) TI - Development of a 3D human osteoblast cell culture model for studying mechanobiology in orthodontics JF - The European Journal of Orthodontics DO - 10.1093/ejo/cjaa017 DA - 2020-09-11 UR - https://www.deepdyve.com/lp/oxford-university-press/development-of-a-3d-human-osteoblast-cell-culture-model-for-studying-mdg2WTebs0 SP - 387 EP - 395 VL - 42 IS - 4 DP - DeepDyve ER -