BioMed Research International
Volume 2015 (2015), Article ID 376474, 14 pages
http://dx.doi.org/10.1155/2015/376474

Research Article

Strain Amplification Analysis of an Osteocyte under Static and Cyclic Loading: A Finite Element Study

Liping Wang,1,2 Jianghui Dong,³ and Cory J. Xian¹

¹Sansom Institute for Health Research, School of Pharmacy and Medical Sciences, University of South Australia, Adelaide, SA 5001, Australia
²College of Mechanical and Electronic Engineering, Shanghai Jianqiao University, Shanghai 201319, China
³School of Computer Science, Engineering and Mathematics, Flinders University, Adelaide, SA 5042, Australia

Received 23 June 2014; Revised 21 September 2014; Accepted 25 September 2014

Academic Editor: Despina Deligianni

Copyright © 2015 Liping Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Osteocytes, the major type of bone cells which reside in their lacunar and canalicular system within the bone matrix, function as biomechanosensors and biomechanotransducers of the bone. Although biomechanical behaviour of the osteocyte-lacunar-canalicular system has been investigated in previous studies mostly using computational 2-dimensional (2D) geometric models, only a few studies have used the 3-dimensional (3D) finite element (FE) model. In the current study, a 3D FE model was used to predict the responses of strain distributions of osteocyte-lacunar-canalicular system analyzed under static and cyclic loads. The strain amplification factor was calculated for all simulations. Effects on the strain of the osteocyte system were investigated under 500, 1500, 2000, and 3000 microstrain loading magnitudes and 1, 5, 10, 40, and 100 Hz loading frequencies. The maximum strain was found to change with loading magnitude and frequency. It was observed that maximum strain under 3000-microstrain loading was higher than those under 500, 1500, and 2000 microstrains. When the loading strain reached the maximum magnitude, the strain amplification factor of 100 Hz was higher than those of the other frequencies. Data from this 3D FE model study suggests that the strain amplification factor of the osteocyte-lacunar-canalicular system increases with loading frequency and loading strain increasing.

1. Introduction

While bone-forming cells osteoblasts and resorbing cells osteoclasts make up only small proportions (<5 and <1%, resp.), osteocytes make up ~90–95% of all bone cells in adult skeletons [1, 2]. Osteocytes are the longest living bone cells, which can live up to decades within their mineralized environment and have the potential to live as long as the organism itself [3]. The osteocytes are believed to be the cells that sense the mechanical loads of the bone. The major known function of osteocytes is to translate the mechanical strain signals into biochemical signals between osteocytes and cells in the bone [4–6]. Osteocytes are embedded within the bone matrix, can remodel the perilacunar matrix (PCM), and are connected to each other by slender cell processes located within the small bony tubes of the canaliculi [1, 4, 7].

Osteocytes are surrounded by the PCM, which is a thin hyaluronan-rich coating. The PCM plays an important role in many cell-surface phenomena; it transfers mechanical signals between the cell body and cell processes and between PCM and the surrounding extracellular matrix (ECM) [8–13]. The PCM represents the entire space between the ECM and the osteocyte cell [14]. Due to these functional relations, the osteocyte-lacunar-canalicular system should be investigated based on the cell-PCM-ECM geometric structure (Figure 1).

(a)

(b)

Figure 1: The osteocyte-lacunar-canalicular system. (a) Photomicrograph of osteocytes embedded in bone matrix (adopted from [34]); (b) schematic geometry of cubic volume of bone surrounding one osteocyte lacuna with canaliculi and processes.

With the known biomechanosensory function of the osteocyte-lacunar-canalicular system in bone, dynamic mechanical stimulus like cyclic loading has been frequently used for the microscale environment biomechanics investigation [15]. Hsieh and Turner studied the anabolic effect of mechanical loading on bone tissue by modulating the loading frequency in the range of 1–10 Hz [16]. Gururaja et al. developed a 2D poroelastic model to investigate load-induced fluid flow in the lacunar-canalicular system of bone subjected to harmonic bending excitations in the range of 1–10⁶ Hz [17]. Wang et al. investigated the net tracer transport in bone via the lacunar-canalicular porosity by using a mathematical model during cyclic mechanical loading at loading frequency 0.5–100 Hz [18].

To date, some researchers have investigated the osteocyte-lacunar-canalicular system by experimental and computational models [14, 17, 19–22], and the unit of microstrain has been widely used in the loading and data analyses in the cell modelling. McCreadie and Hollister found the strain levels in and surrounding osteocytes were 1700 to 2700 microstrains [21]. Burr et al. measured the strain by surgically implanted strain gages and found that the strain was below 2000 microstrains even under conditions of strenuous activity [23]. Verbruggen et al. conducted a computational investigation of the strain amplification of osteocytes in vivo in response to mechanical stimulation under vigorous physiological activity (3000 με loading levels) [14]. In addition, researchers have introduced methodologies to reconstruct biorealistic osteocyte models. A number of methods have been used to study 3D reconstruction of the osteocyte-lacunar-canalicular system. The methods include scanning electron microscopy (SEM) [24–27], serial focused ion beam/scanning electron microscopy (FIB/SEM) [28], transmission electron microscopy (TEM) [29, 30], ultra high voltage electron microscopes (UHVEM) [31], light microscopy (LM) [32, 33], confocal laser scanning microscopy (CLSM) [14, 34–37], atomic force microscopy (AFM) [38–40], computed tomographic (CT) scanning (such as nano-CT and synchrotron radiation micro-CT) [36, 37, 41, 42], and X-ray (like X-ray phase nanotomography and X-ray microscope (TXM)) [43–45]. However, only a few of these models were used in finite element (FE) analysis of the osteocyte-lacunar-canalicular system. Verbruggen and his coworkers developed FE models of lacunar-canalicular network using confocal microscope imaging in rats [14, 46].

While most previous computational geometries were 2D models, several researchers have proposed 3D models (including biorealistic and idealized FE models) for the computational investigation [14, 21, 22, 47–49]. For example, a 3D poroelastic FE model of rat tibia was developed to study the mechanical loads in modulating local flow distributions and concentration gradients within bone tissue. However, due to the simplicity of the geometry, only a thin layer of elements on the endosteal and periosteal surfaces was used in defining the boundary conditions [49]. As results of FE studies depend on many conditions (e.g., geometry, mesh density, element type and order, and material properties), more credible results can be obtained if more accurate geometries are used in the 3D models, when the other parameters have been confirmed.

Thus, although the osteocyte-lacunar-canalicular system has been investigated through the pure experimental tool, mathematics theory, or 2D FE models in previous studies, only a few researchers have used the 3D FE model. In the current study, a 3D FE model was proposed to predict the biomechanical behaviour of the osteocyte-lacunar-canalicular system. Firstly, a 3D geometry of an osteocyte-lacunar-canalicular system incorporating cell body, PCM, and ECM, which is similar to the real dimension, was developed. Secondly, the biomechanical behaviour of the system was investigated by using the 3D model subjected to the static and cyclic loading of different frequencies and magnitudes.

2. Material and Methods

2.1. Geometric Modelling and Meshing

An idealized model of osteocyte-lacunar-canalicular system was developed (Figure 2). The triaxial lacunar osteocyte ellipsoid with the Cartesian equation for a general ellipsoid is shown as [50]

where , , and are the osteocyte lacunar semiaxes in the , , and direction, respectively. is the minor axis, the major axis, and the intermediate axis of the osteocyte lacuna in the local coordinate system. In this study, , , and are 4.6, 9.45, and 2.4, respectively, with the values obtained from human femurs [42]. In this system, the PCM and osteocyte are embedded in the ECM block, whilst the osteocyte is in the center of the ECM block. While the mean length of each direction of ECM is 43 μm for human [50], osteocyte cell bodies are surrounded by a pronounced ~0.5–1 μm thick layer of PCM [51]. The PCM thickness was assumed to be 0.75 μm in this study. Beno et al. [50] estimated the numbers of canaliculi emanating from osteocyte lacunae for different species (six species: chick, rabbit, bovine, horse, dog, and human) by using the slicing method and surface area method based on the data from Remaggi et al. and Ferretti et al. [52, 53]. Canaliculi were idealized as straight cylindrical channels and eighteen canaliculi were modelled as channels in the ECM [50] and using diameter of 0.25 μm [4, 54]. Furthermore, the average width of the canaliculi around the osteocyte process was ~80 nm [4, 55].

(a)

(b)

(c)

(d)

(e)

Figure 2: (a) The schematic diagram of an osteocyte lacuna, where is the intermediate axis, the major axis, and the minor axis of the lacuna in the local coordinate system. Osteocyte lacuna and canaliculi are shown schematically: (b) osteocyte cell body (green) with processes (blue), (c) - plane, (d) - plane, and (e) - plane.

The FE analysis software ABAQUS 6.12 (SIMULIA, Providence, RI, USA) was used for simulations assuming fully saturated media. Because the model is symmetrical, only 1/8 symmetry model was used in all the simulations (Figure 3). The reduced integration first order solid elements in Abaqus suffered from hourglassing when a mesh was coarse. In addition, the fully integration first order solid elements exhibited shear locking. To avoid the elements suffering from hourglassing and shear locking, the reduced integration with second order solid elements is recommended [56]. Twenty-node hexahedral reduced integrated elements (C3D20R) were used for all regions. The hexahedral FE mesh was mapped according to the high quality mesh geometries using twenty-node C3D20R elements in the whole FE analysis. Total number of nodes and elements were 453199 and 107868, respectively, in the 1/8 symmetry model (Table 1). Tie contact interfaces were used to ensure the PCM attached to the ECM and osteocyte and to prevent any relative movement during the simulation.

Table 1: Material properties and element information applied in the FE model of osteocyte lacunae.

Components
Young’s Module, MPa
Poisson’s ratio
Type of element
Number of elements
Reference

ECM
16
0.38
C3D20R
77480
[14]

PCM (with canaliculi)
0.043
0.4
C3D20R
7183
[10]

Osteocyte cell body and processes
0.00447
0.3
C3D20R
23205
[58]

Figure 3: One-eighth symmetry FE mesh model includes ECM, PCM with canaliculi, and osteocyte cell body with processes.

2.2. Materials

All materials were assumed to be isotropic and linearly elastic. The PCM material properties are difficult to measure directly because of its size and connectivity to both cell body and ECM [10]. Meanwhile, it is hard to use the experimental data to define the material properties of the PCM surrounding the osteocyte [14]. An inverse FE approach was used to calculate the linear elastic properties of the PCM. The obtained PCM moduli ranged from 43 to 240 kPa based on this method [10]. Therefore, an elastic modulus of 43 kPa was assigned to the PCM and with Poisson’s ratio of 0.4 [9, 14, 57]. The properties of the cells in the lacuna space were assigned a Young’s modulus of 10 MPa and Poisson’s ratio of 0.4 [21]. A modulus of  Pa and Poisson’s ratio of 0.3 were attributed to the osteocyte cell body and processes [58]. ECM Young’s modulus of 16 GPa and Poisson’s ratio of 0.38 were applied. Young’s modulus of 43 kPa and Poisson’s ratio of 0.4 were used in PCM with canaliculi, and Young’s modulus of 4.47 kPa and Poisson’s ratio of 0.3 were used for osteocyte cell body and processes. The material was modelled as linear elastic isotropic (Table 1).

2.3. Boundary and Loading Conditions

To simulate the mechanical environment of the lacunocanalicular system, static and cyclic loads were applied and the responses of strain and stress distributions were analyzed. The ramped static compressive loads of 500, 1500, 2000, and 3000 microstrain were used in this study [14, 22]. Generally, for most animals, the peak principal compressive strains are about 2000 to 3000 microstrains. However, peak principal strain in the human tibia at the location measured is slightly less than these values. In some cases such as uphill and downhill zigzag running, the principal compressive and shear strains are the greatest, reaching nearly 2000 microstrains, about three times higher than that recorded during walking [23].

Furthermore, a sinusoidal loading was used for compression [59] (Figure 4). The maximal loading magnitudes were 500, 1500, 2000, and 3000 microstrains, respectively, in the simulations. The axial displacement boundary conditions were applied to the model faces to simulate the compression. The nodes on the opposing faces of applied displacement loading were constrained symmetrically to prevent rigid body motion [14, 22, 23]. To examine the influence of frequencies on strain, the loading frequencies of 1, 5, 10, 40, and 100 were investigated [16, 17, 59, 60].

Figure 4: Effects of load cycle on the model.

3. Results

3.1. Static Loading

To simulate the mechanical environment of the lacunocanalicular system, the global static compressive loads of 500, 1500, 2000, and 3000 microstrains were applied and the responses of strain and stress distributions were analyzed. The strain amplification factor was computed as the local maximum strain divided by the applied global strain as described [9, 22]. The contour plots of half symmetry FE model (Figures 5(a)–5(d)) show the predicted strain distribution of the lacunar-canalicular FE model under 500, 1500, 2000, and 3000 microstrains, respectively, and the maximum principal strains were about 1705, 5129, 6848, and 10300 microstrains, respectively. When subjected to a nominal continuum strain level approximately equal to that measured in humans in vivo during rigorous activity like 2,000 microstrains, the osteocyte level strains can be as high as 12,000 to 15,000 microstrains (1.2% to 1.5%) [61]. Obviously, the maximum principal strain is located in the canaliculi near to osteocyte cell body and the interface between canaliculi and osteocyte. For the osteocyte cell, the principal strains in different areas followed the following trend: innermost center > center > outermost section. For PCM, the principal strains in the top and bottom fractions near to osteocyte cell body were higher than the rest, whilst those in the left and right side sections were lower than the other parts in the ECM.

(a)

(b)

(c)

(d)

Figure 5: Strain distribution of half FE model in - plane (left), - plane (middle), and - plane (right), respectively, under different global loading, (a) 500-microstrain global loading, (b) 1500-microstrain global loading, (c) 2000-microstrain global loading, and (d) 3000-microstrain global loading.

The predicted results of strain amplification factors of the FE model were calculated for different loads (Figure 6), which were compared with the data in the literature [14, 22]. The strain amplification factors were 3.41, 3.419, 3.24, and 3.43 for 500, 1500, 2000, and 3000 microstrain global loads, respectively. The strain amplification factors for loads of 500 and 3000 microstrains were comparable with the values of the idealized model and confocal image-derived model from Verbruggen et al. [14], while that for the loading of 1500 microstrains was comparable with the values of the idealized model from Verbruggen et al. [14] and Nicolella et al. [62]. Finally, the strain amplification factor for 2000-microstrain loading was comparable with those from the idealized model of Rath Bonivtch et al. [22] and Nicolella et al. [62].

Figure 6: The strain amplification factor of FE mesh model compared with reference data from the literature.

3.2. Cyclic Loading

The mechanical environment response simulations of the osteocyte-lacunocanalicular system were also investigated under cyclic loads with the loading frequencies of 1, 5, 10, 40, and 100 Hz and the loads of 500, 1500, 2000, 3000, and 3000 microstrains. The results of strains were normalized at time increments of , where is the relevant period of loading. Figures 7, 8, 9, 10, and 11 show curves of strain amplification factors versus    and strain distribution under loading frequencies of 1, 5, 10, 40, and 100 Hz and strain loads of 500, 1500, 2000, 3000, and 3000 microstrain, respectively.

(a)

(b)

(c)

(d)

Figure 7: Strain amplification factors and strain distribution by application of a sinusoidal global strain load of 500 microstrains at the frequency of 1 Hz, (a) strain amplification factors versus one loading period, (b) when , , , , and , strain distribution in - plane, (c) when , , , , and , strain distribution in - plane, and (d) when , , , , and , strain distribution in - plane.

(a)

(b)

(c)

(d)

Figure 8: Strain amplification factors and strain distribution by application of a sinusoidal global strain load of 1500 microstrains at the frequency of 5 Hz, (a) strain amplification factors versus one loading period, (b) when , , , , and , strain distribution in - plane, (c) when , , , , and , strain distribution in - plane, and (d) when , , , , and , strain distribution in - plane.

(a)

(b)

(c)

(d)

Figure 9: Strain amplification factors and strain distribution by application of a sinusoidal global strain load of 2000 microstrains at the frequency of 10 Hz, (a) strain amplification factors versus one loading period, (b) when , , , , and , strain distribution in - plane, (c) when , , , , and , strain distribution in - plane, and (d) when , , , , and , strain distribution in - plane.

(a)

(b)

(c)

(d)

Figure 10: Strain amplification factors and strain distribution by application of a sinusoidal global strain load of 3000 microstrains at the frequency of 40 Hz, (a) strain amplification factors versus one loading period, (b) when , , , , and , strain distribution in - plane, (c) when , , , , and , strain distribution in - plane, and (d) when , , , , and , strain distribution in - plane.

(a)

(b)

(c)

(d)

Figure 11: Strain amplification factors and strain distribution by application of a sinusoidal global strain load of 3000 microstrains at the frequency of 100 Hz, (a) strain amplification factors versus one loading period, (b) when , , , , and , strain distribution in - plane, (c) when , , , , and , strain distribution in - plane, and (d) when , , , , and , strain distribution in - plane.

The influence of loading frequency and loading magnitude on strain amplification was investigated, and the trends of strain amplification factors at different frequencies (1, 5, 10, 40, and 100 Hz) and loading strains (500, 1500, 2000, and 3000 microstrains) were shown in Figure 12 for .

Figure 12: Comparison of strain amplification factors at different frequencies and loading strains all at .

4. Discussion

In this work, a 3D FE model of osteocyte-lacunar-canalicular was developed to investigate the effects of cyclic loading at different frequencies on strain responses of the osteocyte-lacunar-canalicular system. The compressive loads of 500, 1500, 2000, and 3000 microstrains were applied in both static simulation and cyclic simulation. The strain amplification factors of all simulations were also calculated.

Some previous computational studies have examined idealized models of the osteocyte mechanical environment in 3D [21, 22] and have studied both idealized models and biorealistic models [14, 46]; these previous studies represented the first development of accurate 3D FE geometries of the lacuna-osteocyte-canaliculus-ECM system to predict osteocyte mechanobiology. In the previous studies, the shapes of idealized lacunar osteocyte were modelled as revolved ellipsoid with the major and minor axes for idealized models. However, in our study, one triaxial lacunar osteocyte ellipsoid FE model was developed to predict the stain of osteocyte-lacunar-canalicular system for human.

In the static simulation, the strain amplification factors were compared with the existing data in literature (Figure 6). The results of the simulation were found to be in good agreement with the results in literature [14, 22, 62]. The strain amplification factors of simulation were ~3.41, 3.419, 3.424, and 3.43 under 500, 1500, 2000, and 3000 microstrains, respectively, whilst Verbruggen et al. found that these values were 2.2, 2.2, and 4.2 and 7.8, 8.7, and 23.3 for idealized model and confocal image-derived model, respectively [14]. In addition, the strain amplification factor of simulation was about 3.58 under 2000 microstrains, and the result of Rath Bonivtch et al. was 2.45 for idealized model [22]. Compared with the average strain amplification factors (varying from 1.1 to 3.8) from another study [62], the results of our FE analysis were consistent with these measurements. Thus, the simulation of static loading with our system was validated by the reference data in the literature, suggesting that our FE model of osteocyte-lacunar-canalicular is reliable.

For the cyclic loading simulations, the loading frequency and loading magnitude were investigated in this study, and we performed a total of 20 simulations with different loading strain ad 5 different frequencies (1, 5, 10, 40, and 100 Hz) which are within the physiologically meaningful frequency range of 1–100 Hz [17, 59, 63]. The maximum strain increased with increasing loading frequency and increasing loading magnitude (loading strain). The maximum principal strain occurred at , which matched the maximum value of loading strain. Because the sinusoidal loading was applied in the cyclic loading simulations, the principal strain varied with time. When the loading strain reached the maximum magnitude value, the strain amplification factor of 100 Hz was higher than the other frequencies. When , the strain amplification factor increased with loading frequency and with loading strain. However, the strain amplification factor increased only slightly with the increasing loading magnitude for the same loading frequency and increased slightly with the increasing loading frequency for the same loading magnitude.

There has been a lack of experimental data regarding the strain amplification factors for the osteocyte-lacunar-canalicular system. Due to the complex microstructural organization of osteocyte-lacunar-canalicular system, the osteocytes cannot be reliably and easily estimated from global strain measurements [36]. The results of static loading in our study are comparable with FE analysis estimates [14, 22] and digital image measurements [36]. We have to point out that the cyclic loading data of our model is yet to be validated in future studies, as it is beyond the scope of present work. Furthermore, since pore fluid pressure is an important load-induced phenomenon of the osteocyte-lacunar-canalicular system as it in large part dictates the shear stress and degree of chemotransport experienced by the cells, future studies will be required to investigate the pore pressure, flow velocity, and poroelasticity in the 3D poroelastic FE model. Such data will facilitate the understanding of the biomechanical responses to loading and chemotransport around the osteocyte cell body in the osteocyte-lacunar-canalicular system.

Taken together, to understand the changes in strain of osteocytes with different loading frequencies and loading magnitudes, the current study has developed a 3D FE model of osteocyte-lacunar-canalicular system which was used to predict the responses of strain distributions of the osteocyte system under various static and cyclic loads. The study found that the strain amplification factor increased with loading frequency and increasing loading strain.

Conflict of Interests

The authors declared that they have no conflict of interests to this work.

Acknowledgments

Jianghui Dong was supported by Flinders University Ph.D. scholarship and Cory J. Xian is supported by Australian NHMRC Senior Research Fellowship.

References

L. F. Bonewald, “The amazing osteocyte,” Journal of Bone and Mineral Research, vol. 26, no. 2, pp. 229–238, 2011.

T. A. Franz-Odendaal, B. K. Hall, and P. E. Witten, “Buried alive: how osteoblasts become osteocytes,” Developmental Dynamics, vol. 235, no. 1, pp. 176–190, 2006.

M. L. Knothe Tate, J. R. Adamson, A. E. Tami, and T. W. Bauer, “The osteocyte,” International Journal of Biochemistry and Cell Biology, vol. 36, no. 1, pp. 1–8, 2004.

L.-D. You, S. Weinbaum, S. C. Cowin, and M. B. Schaffler, “Ultrastructure of the osteocyte process and its pericellular matrix,” Anatomical Record Part A: Discoveries in Molecular, Cellular, and Evolutionary Biology, vol. 278, no. 2, pp. 505–513, 2004.

K. Zhang, C. Barragan-Adjemian, L. Ye et al., “E11/gp38 selective expression in osteocytes: regulation by mechanical strain and role in dendrite elongation,” Molecular and Cellular Biology, vol. 26, no. 12, pp. 4539–4552, 2006.

L. F. Bonewald, “Mechanosensation and transduction in osteocytes,” BoneKEy Osteovision, vol. 3, no. 10, pp. 7–15, 2006.

J. Klein-Nulend, A. D. Bakker, R. G. Bacabac, A. Vatsa, and S. Weinbaum, “Mechanosensation and transduction in osteocytes,” Bone, vol. 54, no. 2, pp. 182–190, 2013.

J. K. Terzis and K. J. Smith, “Repair of severed peripheral nerves: comparison of the 'de Medinaceli' and standard microsuture methods,” Experimental Neurology, vol. 96, no. 3, pp. 672–680, 1987.

L. G. Alexopoulos, L. A. Setton, and F. Guilak, “The biomechanical role of the chondrocyte pericellular matrix in articular cartilage,” Acta Biomaterialia, vol. 1, no. 3, pp. 317–325, 2005.

A. J. Michalek and J. C. Iatridis, “A numerical study to determine pericellular matrix modulus and evaluate its effects on the micromechanical environment of chondrocytes,” Journal of Biomechanics, vol. 40, no. 6, pp. 1405–1409, 2007.

C. Fotia, G. M. L. Messina, G. Marletta, N. Baldini, and G. Ciapetti, “Hyaluronan-based pericellular matrix: substrate electrostatic charges and early cell adhesion events,” European Cells and Materials, vol. 26, pp. 133–149, 2013.

B. Wang, X. Lai, C. Price et al., “Perlecan-containing pericellular matrix regulates solute transport and mechanosensing within the osteocyte lacunar-canalicular system,” Journal of Bone and Mineral Research, vol. 29, no. 4, pp. 878–891, 2014.

I. Prasadam, S. Farnaghi, J. Q. Feng et al., “Impact of extracellular matrix derived from osteoarthritis subchondral bone osteoblasts on osteocytes: role of integrinβ1 and focal adhesion kinase signaling cues,” Arthritis Research & Therapy, vol. 15, no. 5, article R150, 2013.

S. W. Verbruggen, T. J. Vaughan, and L. M. McNamara, “Strain amplification in bone mechanobiology: a computational investigation of the in vivo mechanics of osteocytes,” Journal of the Royal Society Interface, vol. 9, no. 75, pp. 2735–2744, 2012.

S. C. Cowin, L. Moss-Salentijn, and M. L. Moss, “Candidates for the mechanosensory system in bone,” Journal of Biomechanical Engineering, vol. 113, no. 2, pp. 191–197, 1991.

Y.-F. Hsieh and C. H. Turner, “Effects of loading frequency on mechanically induced bone formation,” Journal of Bone and Mineral Research, vol. 16, no. 5, pp. 918–924, 2001.

S. Gururaja, H. J. Kim, C. C. Swan, R. A. Brand, and R. S. Lakes, “Modeling deformation-induced fluid flow in cortical bone's canalicular-lacunar system,” Annals of Biomedical Engineering, vol. 33, no. 1, pp. 7–25, 2005.

L. Wang, S. C. Cowin, S. Weinbaum, and S. P. Fritton, “Modeling tracer transport in an osteon under cyclic loading,” Annals of Biomedical Engineering, vol. 28, no. 10, pp. 1200–1209, 2000.

G. P. Watchmaker, C. A. Gumucio, R. E. Crandall, M. A. Vannier, and P. M. Weeks, “Fascicular topography of the median nerve: a computer based study to identify branching patterns,” The Journal of Hand Surgery, vol. 16, no. 1, pp. 53–59, 1991.

K. Sun, J. Zhang, T. Chen et al., “Three-dimensional reconstruction and visualization of the median nerve from serial tissue sections,” Microsurgery, vol. 29, no. 7, pp. 573–577, 2009.

B. R. McCreadie and S. J. Hollister, “Strain concentrations surrounding an ellipsoid model of lacunae and osteocytes,” Computer Methods in Biomechanics and Biomedical Engineering, vol. 1, no. 1, pp. 61–68, 1997.

A. Rath Bonivtch, L. F. Bonewald, and D. P. Nicolella, “Tissue strain amplification at the osteocyte lacuna: a microstructural finite element analysis,” Journal of Biomechanics, vol. 40, no. 10, pp. 2199–2206, 2007.

D. B. Burr, C. Milgrom, D. Fyhrie et al., “In vivo measurement of human tibial strains during vigorous activity,” Bone, vol. 18, no. 5, pp. 405–410, 1996.

U. E. Pazzaglia, T. Congiu, V. Sibilia, and D. Quacci, “Osteoblast-osteocyte transformation. A SEM densitometric analysis of endosteal apposition in rabbit femur,” Journal of Anatomy, vol. 224, no. 2, pp. 132–141, 2014.

U. E. Pazzaglia, T. Congiu, P. C. Brunelli, L. Magnano, and A. Benetti, “The long bone deformity of osteogenesis imperfecta III: analysis of structural changes carried out with scanning electron microscopic morphometry,” Calcified Tissue International, vol. 93, no. 5, pp. 453–461, 2013.

U. E. Pazzaglia, T. Congiu, A. Pienazza, M. Zakaria, M. Gnecchi, and C. Dell'Orbo, “Morphometric analysis of osteonal architecture in bones from healthy young human male subjects using scanning electron microscopy,” Journal of Anatomy, vol. 223, no. 3, pp. 242–254, 2013.

G. Marotti, M. A. Muglia, and D. Zaffe, “A SEM study of osteocyte orientation in alternately structured osteons,” Bone, vol. 6, no. 5, pp. 331–334, 1985.

P. Schneider, M. Meier, R. Wepf, and R. Müller, “Serial FIB/SEM imaging for quantitative 3D assessment of the osteocyte lacuno-canalicular network,” Bone, vol. 49, no. 2, pp. 304–311, 2011.

D. Sharma, C. Ciani, P. A. R. Marin, J. D. Levy, S. B. Doty, and S. P. Fritton, “Alterations in the osteocyte lacunar-canalicular microenvironment due to estrogen deficiency,” Bone, vol. 51, no. 3, pp. 488–497, 2012.

F. A. Shah, B. R. Johansson, P. Thomsen, and A. Palmquist, “Ultrastructural evaluation of shrinkage artefacts induced by fixatives and embedding resins on osteocyte processes and pericellular space dimensions,” Journal of Biomedical Materials Research Part A, 2014.

H. Kamioka, Y. Kameo, Y. Imai et al., “Microscale fluid flow analysis in a human osteocyte canaliculus using a realistic high-resolution image-based three-dimensional model,” Integrative Biology, vol. 4, no. 10, pp. 1198–1206, 2012.

U. E. Pazzaglia, T. Congiu, M. Marchese, G. Zarattini, and C. Dell'Orbo, “The canalicular system and the osteoblast domain in human secondary osteons,” Anatomia, Histologia, Embryologia, vol. 41, no. 6, pp. 410–418, 2012.

L. du Plessis and E. van Wilpe, “Thrombocyte morphology and morphometric observations in two vulture species,” Veterinary Clinical Pathology, vol. 38, no. 3, pp. 316–320, 2009.

Y. Sugawara, H. Kamioka, T. Honjo, K.-I. Tezuka, and T. Takano-Yamamoto, “Three-dimensional reconstruction of chick calvarial osteocytes and their cell processes using confocal microscopy,” Bone, vol. 36, no. 5, pp. 877–883, 2005.

H. Kamioka, T. Honjo, and T. Takano-Yamamoto, “A three-dimensional distribution of osteocyte processes revealed by the combination of confocal laser scanning microscopy and differential interference contrast microscopy,” Bone, vol. 28, no. 2, pp. 145–149, 2001.

R. P. van Hove, P. A. Nolte, A. Vatsa et al., “Osteocyte morphology in human tibiae of different bone pathologies with different bone mineral density—is there a role for mechanosensing?” Bone, vol. 45, no. 2, pp. 321–329, 2009.

A. Vatsa, R. G. Breuls, C. M. Semeins, P. L. Salmon, T. H. Smit, and J. Klein-Nulend, “Osteocyte morphology in fibula and calvaria—is there a role for mechanosensing?” Bone, vol. 43, no. 3, pp. 452–458, 2008.

Y. Lin and S. Xu, “AFM analysis of the lacunar-canalicular network in demineralized compact bone,” Journal of Microscopy, vol. 241, no. 3, pp. 291–302, 2011.

K. Ghosh, S. Gangodkar, P. Jain et al., “Imaging the interaction between dengue 2 virus and human blood platelets using atomic force and electron microscopy,” Journal of Electron Microscopy, vol. 57, no. 3, pp. 113–118, 2008.

G. C. Reilly, H. F. Knapp, A. Stemmer, P. Niederer, and M. L. Knothe Tate, “Investigation of the morphology of the lacunocanalicular system of cortical bone using atomic force microscopy,” Annals of Biomedical Engineering, vol. 29, no. 12, pp. 1074–1081, 2001.

P. Dong, A. Pacureanu, M. A. Zuluaga et al., “A new quantitative approach for estimating bone cell connections from nano-CT images,” in Proceedings of the 35th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC '13), pp. 3694–3697, Osaka, Japan, July 2013.

P. Dong, S. Haupert, B. Hesse et al., “3D osteocyte lacunar morphometric properties and distributions in human femoral cortical bone using synchrotron radiation micro-CT images,” Bone, vol. 60, pp. 172–185, 2014.

N. Nango, S. Kubota, A. Takeuchi et al., “Talbot-defocus multiscan tomography using the synchrotron X-ray microscope to study the lacuno-canalicular network in mouse bone,” Biomedical Optics Express, vol. 4, no. 6, pp. 917–923, 2013.

J. C. Andrews, E. Almeida, M. C. van der Meulen et al., “Nanoscale X-ray microscopic imaging of mammalian mineralized tissue,” Microscopy & Microanalysis, vol. 16, no. 3, pp. 327–336, 2010.

M. Langer, A. Pacureanu, H. Suhonen, Q. Grimal, P. Cloetens, and F. Peyrin, “X-ray phase nanotomography resolves the 3D human bone ultrastructure,” PLoS ONE, vol. 7, no. 8, Article ID e35691, 2012.

S. W. Verbruggen, T. J. Vaughan, and L. M. McNamara, “Fluid flow in the osteocyte mechanical environment: a fluid-structure interaction approach,” Biomechanics and Modeling in Mechanobiology, vol. 13, no. 1, pp. 85–97, 2014.

D. D. Deligianni and C. A. Apostolopoulos, “Multilevel finite element modeling for the prediction of local cellular deformation in bone,” Biomechanics and Modeling in Mechanobiology, vol. 7, no. 2, pp. 151–159, 2008.

M. E. Jabaley, W. H. Wallace, and F. R. Heckler, “Internal topography of major nerves of the forearm and hand: a current view,” Journal of Hand Surgery, vol. 5, no. 1, pp. 1–18, 1980.

R. Steck, P. Niederer, and M. L. Knothe Tate, “A finite element analysis for the prediction of load-induced fluid flow and mechanochemical transduction in bone,” Journal of Theoretical Biology, vol. 220, no. 2, pp. 249–259, 2003.

T. Beno, Y.-J. Yoon, S. C. Cowin, and S. P. Fritton, “Estimation of bone permeability using accurate microstructural measurements,” Journal of Biomechanics, vol. 39, no. 13, pp. 2378–2387, 2006.

L. M. McNamara, R. J. Majeska, S. Weinbaum, V. Friedrich, and M. B. Schaffler, “Attachment of osteocyte cell processes to the bone matrix,” Anatomical Record, vol. 292, no. 3, pp. 355–363, 2009.

F. Remaggi, V. Canè, C. Palumbo, and M. Ferretti, “Histomorphometric study on the osteocyte lacuno-canalicular network in animals of different species. I. Woven-fibered and parallel-fibered bones,” Italian Journal of Anatomy and Embryology, vol. 103, no. 4, pp. 145–155, 1998.

M. Ferretti, M. A. Muglia, F. Remaggi, V. Canè, and C. Palumbo, “Histomorphometric study on the osteocyte lacuno-canalicular network in animals of different species. II. Parallel-fibered and lamellar bones,” Italian Journal of Anatomy and Embryology, vol. 104, no. 3, pp. 121–131, 1999.

U. E. Pazzaglia and T. Congiu, “The cast imaging of the osteon lacunar-canalicular system and the implications with functional models of intracanalicular flow,” Journal of Anatomy, vol. 222, no. 2, pp. 193–202, 2013.

L. Wang, Y. Wang, Y. Han et al., “In situ measurement of solute transport in the bone lacunar-canalicular system,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 33, pp. 11911–11916, 2005.

A. Tsouknidas, S. Savvakis, Y. Asaniotis, K. Anagnostidis, A. Lontos, and N. Michailidis, “The effect of kyphoplasty parameters on the dynamic load transfer within the lumbar spine considering the response of a bio-realistic spine segment,” Clinical Biomechanics, vol. 28, no. 9-10, pp. 949–955, 2013.

L. G. Alexopoulos, M. A. Haider, T. P. Vail, and F. Guilak, “Alterations in the mechanical properties of the human chondrocyte pericellular matrix with osteoarthritis,” Journal of Biomechanical Engineering, vol. 125, no. 3, pp. 323–333, 2003.

Y. Sugawara, R. Ando, H. Kamioka et al., “The alteration of a mechanical property of bone cells during the process of changing from osteoblasts to osteocytes,” Bone, vol. 43, no. 1, pp. 19–24, 2008.

G. C. Goulet, D. Coombe, R. J. Martinuzzi, and R. F. Zernicke, “Poroelastic evaluation of fluid movement through the lacunocanalicular system,” Annals of Biomedical Engineering, vol. 37, no. 7, pp. 1390–1402, 2009.

C. H. Turner, M. R. Forwood, and M. W. Otter, “Mechanotransduction in bone: Do bone cells act as sensors of fluid flow?” FASEB Journal, vol. 8, no. 11, pp. 875–878, 1994.

D. P. Nicolella and J. Lankford, “Microstructural strain near osteocyte lacuna in cortical bone in vitro,” Journal of Musculoskeletal Neuronal Interactions, vol. 2, no. 3, pp. 261–263, 2002.

D. P. Nicolella, D. E. Moravits, A. M. Gale, L. F. Bonewald, and J. Lankford, “Osteocyte lacunae tissue strain in cortical bone,” Journal of Biomechanics, vol. 39, no. 9, pp. 1735–1743, 2006.

S. P. Fritton, K. J. McLeod, and C. T. Rubin, “Quantifying the strain history of bone: spatial uniformity and self-similarity of low-magnitude strains,” Journal of Biomechanics, vol. 33, no. 3, pp. 317–325, 2000.