Name: Analysis and Imaging of Osteocytes (Video) | JoVE Journal
Uploaded: 2024-11-29T12:00:00.000+00:00
Duration: 10 min 19 s
Description: Learn computational fluid dynamics analysis using confocal microscopy and finite element modeling to evaluate fluid flow shear stress on osteocyte dendritic membranes in bone mechanosensation research.

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

The authors would like to acknowledge the National Science Foundation (NSF, award number NSF-CMMI-1662284 PI: T Ganesh), National Institute of Health (NIH - NIA P01 AG039355 PI: LF Bonewald) and (NIH/SIG S10OD021665 and S10RR027668 PI: SL Dallas), and the University of Missouri-Kansas City School of Graduate Studies Research Grant Program.

Animal experiments were carried out with the approval of the Institutional Animal Care and Use Committee at the University of Missouri, Kansas City (UMKC), and conformed to relevant federal guidelines.
1. Bone preparation process
<ol>
	<li>Collect femurs from 4-month-old and 22-month-old female C57BL6 mice and fix them in cold 4% paraformaldehyde in PBS for 24 h at 4 &#176;C with gentle rocking, then rinse them in PBS and store them in 70% ethanol before embedding. 
	NOTE: Fixative volume should be about 20x the tissue volume</li>
	<li>Embed bones quickly in a fast polymerizing acrylic (Table of Materials) following the manufacturer&#39;s instructions. 
	NOTE: It is important to use a fast polymerizing resin for this step (~10 min). The purpose is to support the bone tissue during sectioning using the diamond saw, but without the resin penetrating into the LCN, which would block the FITC stain from penetrating in.</li>
	<li>Cut thick 300 &#181;m thick transverse slices from a standardized site above the third trochanter using a diamond saw and store them at 4 &#176;C in 70% ethanol prior to FITC staining.</li>
	<li>Polish the sections using 600, 800, and then 1200 grit sandpaper to a final thickness of ~90-100 &#956;m. 
	NOTE: Attaining proper thickness was ensured using a digital caliper.</li>
	<li>Rinse the sections in 70%, 95%, and 100% ethanol for 5 min each.</li>
	<li>Stain them in 1% FITC in 100% ethanol for 4 h at room temperature (RT) in the dark with moderate shaking.</li>
	<li>Wash the sections in 100% ethanol for 30 min with gentle shaking in the dark. Then, air-dry them overnight in the dark.</li>
	<li>To mount, place the section into a drop of mounting media on a glass microscope slide. Use forceps to position the section as flat as possible against the slide, avoiding generating air bubbles and using the surrounding resin to manipulate the sample. Mount a coverslip on the specimen.</li>
</ol>
2. Confocal microscopy
<ol>
	<li>Use a confocal microscope for imaging of FITC-stained bone slices.</li>
	<li>Use a 100x 1.44NA oil objective with a digital zoom of 1.7 and a step size of 0.126 &#956;m to collect detailed Z-stacks of 400 Z-planes at 1024 x 1024 pixels, 0.089 &#956;m pixel resolution.</li>
	<li>Use a 488 nm laser for excitation, with an emission collection window of 496-596 nm. Collect the image stacks using compensation settings to correct for signal loss with increasing imaging depth.</li>
	<li>Increase the accuracy and resolution of the images by using image collection techniques such as oversampling and increased line averaging. Also, collect the images with a magnification of 5x, 20x, and 100x resolution of transverse sections of femur sections as shown in Figure 1. 
	NOTE: The low-resolution (5x) image in Figure 1 shows the complete cross-section area of the femur, with three regions selected for 100x imaging fields.</li>
	<li>Use the 100x Z-stacks for the computer modeling of osteocytes.</li>
</ol>
3. Computer modeling
<ol>
	<li>Import the collected 100x images, in TIFF&#160;format, to ImageJ software to build a sequence of images of the LCN in the Z-direction.</li>
	<li>Import the Z-stacks into the 3D image-based processing software to construct a mask of the LCN after defining the image orientation.</li>
	<li>Threshold the original image from the young and aged mice between 30,012-45,677 Hounsfield units and 15,000-46,701 Hounsfield units, respectively, to closely resemble the LCN. Adjust the Threshold in the section menu to change pixel intensity limits for inclusion in a mask.</li>
	<li>Crop one lacuna with its canaliculi as the region of interest (ROI) from the stack using the Crop Mask operation. Define the ROI so that it encircles the lacuna in the center of the cube, and all its connected canaliculi extend to the sides of the cube. Encase the lacuna in an imaginary bigger cube with side lengths of 21 &#181;m, 14 &#181;m, and 19 &#181;m.</li>
	<li>Since the model is created out of multiple parts, perform a region grow operation to select the connected pixel regions, remove noise and despeckle to generate a uniform LCN.</li>
	<li>Convert the lacunar-canalicular mask into an object using the Calculated Part operation in the 3D image-based processing software.</li>
	<li>Build the osteocyte and dendritic membranes by reducing the LCN volume using the Smoothing operation. Perform this operation several times to achieve a lacunar and canalicular space thickness of 0.75 &#181;m and 0.08 &#181;m, respectively9,18.</li>
	<li>Export the objects (STL format) as the last step in the 3D image-based processing software.</li>
	<li>Import two layers of the LCN and osteocyte-dendritic membranes into the 3D data optimization software to generate a volume mesh.</li>
	<li>Utilize the Fix Wizard tool in the software to identify mesh problems in each part. Check the mesh quality in the diagnostic section of Fix Wizard after each operation.</li>
	<li>Remove the inverted normal parts, intersecting triangles, and bad contours using the Auto-Fix operation in Fix Wizard.</li>
	<li>Replace the overlapping triangles manually by defining new ones or automatically via the Fill Hole Normal operation.</li>
	<li>Improve the mesh quality using operations including filtering sharp triangles, small edges, and small shells.</li>
	<li>After improving the mesh quality, combine two surfaces of the LCN and osteocyte-dendritic membranes into one surface (lacunar-canalicular fluidic space) that belongs to both parts using a non-manifold assembly.</li>
	<li>Create a volumetric model of the lacunar-canalicular space using the Remesh operation and then export it as an STL file. Adjust the object scale to be in micrometers in the export section.</li>
</ol>
4. Geometry modification technique in the 3D image-based processing software and 3D data optimization software
NOTE: The geometry modification technique is used to model changes in osteocyte morphology, such as canalicular density and diameter and lacunar-canalicular thickness owing to aging or bone disease.
<ol>
	<li>Choose the young osteocyte as the base model and modify it to build other distinct osteocyte models by applying morphological alterations.</li>
	<li>Generate an osteocyte model with different canalicular densities from the base model by changing the image thresholding in the 3D image-based processing software.
	<ol>
		<li>Select a lower threshold to reduce the light intensity of the image and obtain a lacuna with fewer canaliculi. The advantage of the thresholding technique is that the lacuna shape and size remain the same, and only the influence of canalicular density is studied. Figure 2 shows the simulated aged model generated from the young osteocyte using the geometry modification technique.</li>
	</ol>
	</li>
	<li>Develop osteocyte models with different lacunar-canalicular space thicknesses or dendrite/canalicular diameters in the 3D image-based processing software and 3D data optimization software. Build larger or smaller osteocyte models by Wrapping or Smoothing operations, respectively. Figure 3 shows six osteocyte models with altered geometry developed from the young osteocyte.</li>
</ol>
5. CFD analysis
NOTE: After generating the volumetric osteocyte models, several steps, including geometry, mesh, and setup, are conducted in the CFX module of the simulation software.
<ol>
	<li>Create a Fluid Flow in the simulation software to prepare the models for the CFD analysis.</li>
	<li>Import the developed confocal image-based geometries into the geometry section of CFX, known as ANSYS SpaceClaim (3D modeling tool). Set the unit dimensions to nanometers in the setting.</li>
	<li>Geometry appears as two facets of the LCN and osteocyte-dendritic processes. Click on Facet on the top menu and remove the geometrical errors such as intersections, sharp or overconnected edges, and vertices, openings, or holes for each facet.</li>
	<li>Click on Subtract on the Facet menu to reduce the smaller facet and osteocyte-dendritic processes from the larger facet, the LCN, to achieve a single body of lacunar-canalicular space. Then, right-click on the generated facet and convert it from facets to a solid domain without merging faces. Figure 4 depicts the cross-sectional area of the young osteocyte model, which represents the lacunar-canalicular space.</li>
	<li>Click on mesh and select linear tetrahedral elements using an element size of 0.06 &#181;m. Refine the mesh with a mesh convergence study to have enough elements in the tiny dendritic system to ensure the results are independent of the mesh size.</li>
	<li>Select the surface and choose the canaliculi on the top side of the imaginary cube as fluid inlets. Select the canaliculi on the other five faces as fluid outlets using Box select.</li>
	<li>Export the mesh (fluent file format) as it loads faster in the setup in the next step.</li>
	<li>Create another Fluid Flow in the simulation software and import the fluent mesh into the setup section of CFX. Define two boundary conditions of inlets and outlets for the faces preselected as inlets/outlets using the Insert Boundary option.</li>
	<li>To imitate physiological conditions, exert a fluid inlet pressure of 300 Pa and 0 Pa on the inlets and outlets, respectively19,20. Treat the remaining surfaces as walls with a no-slip condition in that fluid which has zero velocity at the interface of walls. Fluid flows from the inlets around the dendrites and osteocyte cell body and exits from the other canaliculi assigned as outlets.</li>
	<li>Treat the interstitial laminar fluid as water9, chosen from the material library. Set the Heat Transfer, Combustion, and Thermal Radiation sections to None, as no heat transfer is defined in the problem. Select the Turbulence mode as the fluid characteristic in the LCN, which is a laminar fluid9.</li>
	<li>Run the software using Double Precision and Direct Start as the submission type. Monitor mass and momentum until residuals drop and become constant. After solution convergence, measure FFSS data using the CFD-post section of the CFD software.</li>
</ol>
6. CFD post processing
<ol>
	<li>To depict FFSS experienced by osteocytes and their dendrites, insert a new contour in the results section of the CFD software. Create an FFSS contour by choosing the wall shear on the osteocyte-dendritic membranes as the variable at the domain. 
	NOTE: To better show high FFSS on dendritic membranes, the range of FFSS is set to User Specified to modify min/max values of FFSS.</li>
	<li>Insert a velocity streamlines contour inside the lacunar-canalicular domain starting from the inlets. Set the sampling to equally spaced and select the number of points as 2500. An animation section in the CFD software accurately displays in 3D how the fluid particles flow inside the lacunar-canalicular spaces using the velocity streamline graph.</li>
	<li>Use the Function Calculator tool in the CFD software for analyzing the magnitude of FFSS or velocity based on geometrical parameters, especially where there are various osteocyte models (i.e., young vs. aged). Measure the volume and surface area of the lacunar-canalicular space as geometrical parameters along with maximum, minimum, or average FFSS values.</li>
</ol>

3D Modeling of Osteocytes for CFD Analysis to Study Mechanosensation in Bone

<ol><li>Li, M. C. M., Chow, S. K. H., Wong, R. M. Y., Qin, L., Cheung, W. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+role+of+osteocytes-specific+molecular+mechanism+in+regulation+of+mechanotransduction+-+A+systematic+review%5BTitle%5D)+AND+%22J+Orthop+Translat%22%5BJournal%5D)">The role of osteocytes-specific molecular mechanism in regulation of mechanotransduction - A systematic review.</a> J Orthop Translat. 29, 1-9 (2021).</li>
<li>Fritton, S. P., Weinbaum, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Fluid+and+solute+transport+in+bone%3A+Flow-induced+mechanotransduction%5BTitle%5D)+AND+%22Annu+Rev+Fluid+Mech%22%5BJournal%5D)">Fluid and solute transport in bone: Flow-induced mechanotransduction.</a> Annu Rev Fluid Mech. 41, 347-374 (2009).</li>
<li>Klein-Nulend, J., Bacabac, R. G., Bakker, A. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Mechanical+loading+and+how+it+affects+bone+cells%3A+the+role+of+the+osteocyte+cytoskeleton+in+maintaining+our+skeleton%5BTitle%5D)+AND+%22Eur+Cell+Mater%22%5BJournal%5D)">Mechanical loading and how it affects bone cells: the role of the osteocyte cytoskeleton in maintaining our skeleton.</a> Eur Cell Mater. 24, 278-291 (2012).</li>
<li>Knothe Tate, M. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((%26%2334%3BWhither+flows+the+fluid+in+bone%3F%26%2334%3B+An+osteocyte%27s+perspective%5BTitle%5D)+AND+%22J+Biomech%22%5BJournal%5D)">&#34;Whither flows the fluid in bone?&#34; An osteocyte's perspective.</a> J Biomech. 36 (10), 1409-1424 (2003).</li>
<li>Tiede-Lewis, L. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Degeneration+of+the+osteocyte+network+in+the+C57BL%2F6+mouse+model+of+aging%5BTitle%5D)+AND+%22Aging+%28Albany+NY%29%22%5BJournal%5D)">Degeneration of the osteocyte network in the C57BL/6 mouse model of aging.</a> Aging (Albany NY). 9 (10), 2190-2208 (2017).</li>
<li>Lai, X., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+dependences+of+osteocyte+network+on+bone+compartment%2C+age%2C+and+disease%5BTitle%5D)+AND+%22Bone+Res%22%5BJournal%5D)">The dependences of osteocyte network on bone compartment, age, and disease.</a> Bone Res. 3 (1), 15009(2015).</li>
<li>Schurman, C. A., Verbruggen, S. W., Alliston, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Disrupted+osteocyte+connectivity+and+pericellular+fluid+flow+in+bone+with+aging+and+defective+TGF-%CE%B2+signaling%5BTitle%5D)+AND+%22Proc+Natl+Acad+Sci+U+S+A%22%5BJournal%5D)">Disrupted osteocyte connectivity and pericellular fluid flow in bone with aging and defective TGF-β signaling.</a> Proc Natl Acad Sci U S A. 118 (25), e2023999118(2021).</li>
<li>van Tol, A. F., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+mechanoresponse+of+bone+is+closely+related+to+the+osteocyte+lacunocanalicular+network+architecture%5BTitle%5D)+AND+%22Proc+Natl+Acad+Sci+U+S+A%22%5BJournal%5D)">The mechanoresponse of bone is closely related to the osteocyte lacunocanalicular network architecture.</a> Proc Natl Acad Sci U S A. 117 (51), 32251-32259 (2020).</li>
<li>Verbruggen, S. W., Vaughan, T. J., McNamara, L. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Fluid+flow+in+the+osteocyte+mechanical+environment%3A+a+fluid-structure+interaction+approach%5BTitle%5D)+AND+%22Biomech+Model+Mechanobiol%22%5BJournal%5D)">Fluid flow in the osteocyte mechanical environment: a fluid-structure interaction approach.</a> Biomech Model Mechanobiol. 13 (1), 85-97 (2014).</li>
<li>Ganesh, T., Laughrey, L. E., Niroobakhsh, M., Lara-Castillo, N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Multiscale+finite+element+modeling+of+mechanical+strains+and+fluid+flow+in+osteocyte+lacunocanalicular+system%5BTitle%5D)+AND+%22Bone%22%5BJournal%5D)">Multiscale finite element modeling of mechanical strains and fluid flow in osteocyte lacunocanalicular system.</a> Bone. 137, 115328(2020).</li>
<li>Schneider, P., Meier, M., Wepf, R., Müller, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Serial+FIB%2FSEM+imaging+for+quantitative+3D+assessment+of+the+osteocyte+lacuno-canalicular+network%5BTitle%5D)+AND+%22Bone%22%5BJournal%5D)">Serial FIB/SEM imaging for quantitative 3D assessment of the osteocyte lacuno-canalicular network.</a> Bone. 49 (2), 304-311 (2011).</li>
<li>Kamioka, H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+method+for+observing+silver-stained+osteocytes+in+situ+in+3-microm+sections+using+ultra-high+voltage+electron+microscopy+tomography%5BTitle%5D)+AND+%22Microsc+Microanal%22%5BJournal%5D)">A method for observing silver-stained osteocytes in situ in 3-microm sections using ultra-high voltage electron microscopy tomography.</a> Microsc Microanal. 15 (5), 377-383 (2009).</li>
<li>Ciani, C., Doty, S. B., Fritton, S. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((An+effective+histological+staining+process+to+visualize+bone+interstitial+fluid+space+using+confocal+microscopy%5BTitle%5D)+AND+%22Bone%22%5BJournal%5D)">An effective histological staining process to visualize bone interstitial fluid space using confocal microscopy.</a> Bone. 44 (5), 1015-1017 (2009).</li>
<li>Sharma, D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Alterations+in+the+osteocyte+lacunar-canalicular+microenvironment+due+to+estrogen+deficiency%5BTitle%5D)+AND+%22Bone%22%5BJournal%5D)">Alterations in the osteocyte lacunar-canalicular microenvironment due to estrogen deficiency.</a> Bone. 51 (3), 488-497 (2012).</li>
<li>Verbruggen, S. W., Vaughan, T. J., McNamara, L. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Strain+amplification+in+bone+mechanobiology%3A+a+computational+investigation+of+the+in+vivo+mechanics+of+osteocytes%5BTitle%5D)+AND+%22J+R+Soc+Interface%22%5BJournal%5D)">Strain amplification in bone mechanobiology: a computational investigation of the in vivo mechanics of osteocytes.</a> J R Soc Interface. 9 (75), 2735-2744 (2012).</li>
<li>Goggin, P. M., Zygalakis, K. C., Oreffo, R. O., Schneider, P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((High-resolution+3D+imaging+of+osteocytes+and+computational+modelling+in+mechanobiology%3A+insights+on+bone+development%2C+ageing%2C+health+and+disease%5BTitle%5D)+AND+%22Eur+Cell+Mater%22%5BJournal%5D)">High-resolution 3D imaging of osteocytes and computational modelling in mechanobiology: insights on bone development, ageing, health and disease.</a> Eur Cell Mater. 31, 264-295 (2016).</li>
<li>Kamioka, H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Microscale+fluid+flow+analysis+in+a+human+osteocyte+canaliculus+using+a+realistic+high-resolution+image-based+three-dimensional+model%5BTitle%5D)+AND+%22Integr+Biol%22%5BJournal%5D)">Microscale fluid flow analysis in a human osteocyte canaliculus using a realistic high-resolution image-based three-dimensional model.</a> Integr Biol. 4 (10), 1198-1206 (2012).</li>
<li>Wang, L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((In+situ+measurement+of+solute+transport+in+the+bone+lacunar-canalicular+system%5BTitle%5D)+AND+%22Proc+Natl+Acad+Sci+U+S+A%22%5BJournal%5D)">In situ measurement of solute transport in the bone lacunar-canalicular system.</a> Proc Natl Acad Sci U S A. 102 (33), 11911-11916 (2005).</li>
<li>Manfredini, P., Cocchetti, G., Maier, G., Redaelli, A., Montevecchi, F. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Poroelastic+finite+element+analysis+of+a+bone+specimen+under+cyclic+loading%5BTitle%5D)+AND+%22J+Biomech%22%5BJournal%5D)">Poroelastic finite element analysis of a bone specimen under cyclic loading.</a> J Biomech. 32 (2), 135-144 (1999).</li>
<li>Steck, R., Niederer, P., Knothe Tate, M. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+finite+element+analysis+for+the+prediction+of+load-induced+fluid+flow+and+mechanochemical+transduction+in+bone%5BTitle%5D)+AND+%22J+Theor+Biol%22%5BJournal%5D)">A finite element analysis for the prediction of load-induced fluid flow and mechanochemical transduction in bone.</a> J Theor Biol. 220 (2), 249-259 (2003).</li>
<li>You, L., Cowin, S. C., Schaffler, M. B., Weinbaum, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+model+for+strain+amplification+in+the+actin+cytoskeleton+of+osteocytes+due+to+fluid+drag+on+pericellular+matrix%5BTitle%5D)+AND+%22J+Biomech%22%5BJournal%5D)">A model for strain amplification in the actin cytoskeleton of osteocytes due to fluid drag on pericellular matrix.</a> J Biomech. 34 (11), 1375-1386 (2001).</li>
<li>Weinbaum, S., Cowin, S. C., Zeng, Y. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+model+for+the+excitation+of+osteocytes+by+mechanical+loading-induced+bone+fluid+shear+stresses%5BTitle%5D)+AND+%22J+Biomech%22%5BJournal%5D)">A model for the excitation of osteocytes by mechanical loading-induced bone fluid shear stresses.</a> J Biomech. 27 (3), 339-360 (1994).</li>
<li>Niroobakhsh, M., Laughrey, L. E., Dallas, S. L., Johnson, M. L., Ganesh, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Computational+modeling+based+on+confocal+imaging+predicts+changes+in+osteocyte+and+dendrite+shear+stress+due+to+canalicular+loss+with+aging%5BTitle%5D)+AND+%22Biomech+Model+Mechanobiol%22%5BJournal%5D)">Computational modeling based on confocal imaging predicts changes in osteocyte and dendrite shear stress due to canalicular loss with aging.</a> Biomech Model Mechanobiol. 23 (1), 129-143 (2024).</li>
<li>Dallas, S. L., Moore, D. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Using+confocal+imaging+approaches+to+understand+the+structure+and+function+of+osteocytes+and+the+lacunocanalicular+network%5BTitle%5D)+AND+%22Bone%22%5BJournal%5D)">Using confocal imaging approaches to understand the structure and function of osteocytes and the lacunocanalicular network.</a> Bone. 138, 115463(2020).</li>
<li>Boyde, A., Wolfe, L. A., Maly, M., Jones, S. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Vital+confocal+microscopy+in+bone%5BTitle%5D)+AND+%22Scanning%22%5BJournal%5D)">Vital confocal microscopy in bone.</a> Scanning. 17 (2), 72-85 (1995).</li>
</ol>

The authors have nothing to disclose.

This protocol outlines a confocal imaging technique for visualization and computational modeling of the osteocytes. Before confocal imaging, the bone preparation process for sectioning and staining bone samples is performed. Confocal images of 100x magnification are imported into various software to develop computer models of osteocytes and the lacunar-canalicular space. A CFD analysis is conducted lastly on the confocal image-based models to model FFSS surrounding the osteocytes and dendritic membranes due to physical activity, which is technically challenging in vivo due to the inaccessibility of osteocytes in bone matrix. The provided protocol is also helpful in determining other stress or strains using a finite element or fluid-structure interaction (FSI) approach in multiphysics modeling.
Confocal imaging has been used to visualize the LCN7,9,24,25 owing to the high-resolution images that can be achieved using FITC staining of the bone. Although the diameter of the canaliculi is close to the resolution limit of the microscope, reliable confocal images of the canaliculi can be obtained, which can capture most of the canaliculi. Due to the smaller diameter of the osteocyte dendrites, it is technically more challenging to accurately capture images of the dendrites. A helpful technique was therefore developed to model the osteocyte and dendrites from confocal images of the LCN by reducing their volume using a smoothing operation in the software. Some other published studies used a similar approach to offset the LCN to generate the osteocyte-dendritic membrane9. The other benefit of the smoothing technique is in improving some imaging errors in CLSM, such as the elongation of canalicular diameter that occurs in the Z-axis due to the lower imaging resolution in the Z-axis, which is a feature of optical imaging. Deconvolution is also considered an alternative method to decrease the optical stretching in the Z-direction since it improves the signal-to-noise in the original data set. The lacunar-canalicular space is modeled from LCN and osteocyte-dendritic membranes to analyze the fluid profile inside the LCN.
Notably, a new technique was introduced with respect to existing methods to develop geometry-modified osteocyte models without sacrificing another mouse. The geometry modification technique also helps to simulate the alteration of bone microstructure due to aging or bone diseases such as osteoporosis, diabetes and conditions such as perlecan deficiency, etc. In this technique, a confocal image-based osteocyte model of a young mouse was used to model the morphology of aged or other bone-diseased osteocytes by modifying parameters such as canalicular density, lacunar-canalicular space, and dendrite diameter). This approach was validated by the similarity of the two aged osteocyte models, one from confocal image-based modeling of an actual osteocyte from a 22-month-old mouse and the other from geometry modification modeling of an osteocyte originally imaged from a 4-month-old mouse by reducing its canalicular density. These two aged osteocyte models showed a similar FFSS magnitude, which was the lowest among all osteocyte models. The geometry modification technique could have many applications in different biological&#160;fields as it can predict comparable results to osteocyte models generated using real osteocytes from an animal. Some of the limitations of this work include losing some of the detail of surface roughness due to the smoothing technique and the fact that the canalicular diameters are at the limit of resolution of the confocal microscopes, which could potentially lead to some inaccuracy for canaliculi with smaller diameters, and also fluid flow idealization instead of modeling using oscillatory flow.
To demonstrate the technique&#39;s functionality, osteocyte FFSS values were measured for distinct osteocyte models generated from a young and an aged mouse or using the geometry modification technique. One of the noteworthy findings using these computer models was that the volume of lacunar-canalicular space is directly proportional to the average FFSS experienced by the osteocytes. This new finding, which had not been investigated before, may have various applications in determining the osteocyte mechano-responsiveness as it can correlate the osteocyte network morphology to detected FFSS by the osteocytes and their processes.

Osteocytes are postulated to regulate bone mass in response to physical exercise1. Membrane deformation of osteocytes and their dendritic processes due to mechanical loading, subjects them to FFSS, which is detected by the osteocytes and triggers intracellular signaling2,3,4. Bone microstructure undergoes through deterioration or alterations in its lacunar-canalicular morphology due to aging or bone diseases such as osteoporosis and diabetes and in conditions such as perlecan deficiency that causes impaired mechano-responsiveness of osteocytes5,6. These changes in bone architecture cause osteocytes to experience different levels of FFSS and strains7,8. Importantly, FFSS experienced by osteocytes in response to mechanical loading is difficult to quantify in vivo because they are embedded in the calcified bone matrix.
Confocal image-based modeling is a powerful technique to overcome the limitations of studying inaccessible osteocytes in their natural environment by replicating computer models of the LCN9,10. Processing and modeling the interconnected network of LCN in 3D has been challenging. There are several imaging techniques, such as Transmission electron microscopy (TEM), scanning electron microscopy (SEM), serial block face sectioning, and serial focused ion beam scanning electron microscopy (FIB/SEM)2,11,12. A valuable technique was developed to visualize bone13,14,15 and generate 3D osteocyte models via confocal laser scanning microscopy (CLSM). CLSM was chosen here for computational modeling rather than other imaging techniques due to its ability to image all of the lacuna volume and most of the canaliculi in 3D16,17. The LCN geometry can be generated using CLSM for osteocyte Finite Element Analysis (FEA) to predict bone strains. However, fluid analysis to predict FFSS experienced by osteocytes is more complicated as it requires modeling of the cell membrane of the osteocyte and its dendrites within the LCN to enable modeling of the narrow lacunar-canalicular space, in which the interstitial fluid moves around18.
In this protocol, fluorescein isothiocyanate (FITC) dye is applied to undecalcified thick bone sections before confocal microscopy to label the LCN inside the bone, and osteocyte-dendritic membranes are modeled based on imaging data from the LCN. The lacunar-canalicular space is simulated using computational modeling, and physiological loading due to physical activity is modeled using a CFD approach. The osteocytes are subjected to a fluid pressure gradient in the CFD software to analyze the fluid profile inside the LCN and measure FFSS on the osteocyte and dendritic membranes. Furthermore, an FEA&#160;approach can measure osteocyte strains or stresses by applying compressive mechanical loading.
A geometry modification technique was also developed to modify the microstructures derived from images of young, healthy bone in order to simulate the altered lacunar-canalicular morphology in aged animals or those&#160;with bone disease. Alterations of the bone microstructure included reducing the number of canaliculi with aging, reducing the lacunar-canalicular space area to model what happens in perlecan deficiency and increasing it to model aging effects, and reducing the canalicular and dendritic wall area to model diabetic bone5,6. The geometry modification technique allows us to compare FFSS experienced by osteocytes in bone with different microstructures, such as young versus aged or bones in healthy versus diseased animals.
Overall, confocal image-based modeling is a valuable tool for simulating the morphology of osteocytes in healthy bone as well as in aging/disease-associated changes in osteocyte morphology. Furthermore, osteocyte morphological parameters, such as surface area and volume of the lacunar-canalicular space, can be measured and compared in various bones to predict cellular responses to mechanical strain.

<table><tbody><tr><td>1,200 Grit sandpaper</td><td>Buehler</td><td>30-5170-012-100</td><td></td></tr><tr><td>3-Matic software</td><td>Materialise</td><td>https://www.materialise.com/en/industrial/software/3-matic</td><td>3D data optimization software</td></tr><tr><td>600 grit sandpaper</td><td>Buehler</td><td>30-5118-600-100</td><td></td></tr><tr><td>800 Grit sandpaper</td><td>Buehler</td><td>30-5170-800-100</td><td></td></tr><tr><td>ANSYS software</td><td>ANSYS</td><td>https://www.ansys.com/</td><td>simulation software</td></tr><tr><td>Fluorescein Isothiocyanate (FITC)</td><td>Sigma-Aldrich</td><td>F7250</td><td></td></tr><tr><td>ImageJ software</td><td></td><td>https://imagej.net/ij/</td><td></td></tr><tr><td>Immersion Oil for Microscopes</td><td>Leica Microsystems</td><td>195371-10-9</td><td></td></tr><tr><td>Leica TCS Sp5 II confocal microscope&amp;nbsp;</td><td>Leica Microsystems</td><td>TCS Sp5 II&amp;nbsp;</td><td></td></tr><tr><td>Leitz 1600 inner hole diamond saw</td><td>Leica&amp;nbsp;</td><td></td><td></td></tr><tr><td>MIMICS Innovation Suite software</td><td>Materialise</td><td>https://www.materialise.com/en/healthcare/mimics-innovation-suite</td><td>3D image-based processing software</td></tr><tr><td>Permount mount medium</td><td>Fisher scientific</td><td>SP15-500</td><td></td></tr><tr><td>Sampl-Kwick Fast Cure Acrylic Kit</td><td>Buehler</td><td>20-3560</td><td></td></tr><tr><td>Single Platform Laboratory Shaker</td><td>Reliable scientific INC</td><td>Model 55S</td><td></td></tr></tbody></table>

analysis and imaging of osteocytes

Analysis and Imaging of Osteocytes

Osteocytes are the bone cells that are thought to respond to mechanical strains and fluid flow shear stress (FFSS) by activating various biological pathways in a process known as mechanotransduction. Confocal image-derived models of osteocyte networks are a valuable tool for conducting Computational Fluid Dynamics&#160;(CFD) analysis to evaluate shear stresses on the osteocyte membrane, which cannot be determined by direct measurement. Computational modeling using these high-resolution images of the microstructural architecture of bone was used to numerically simulate the mechanical loading exerted on bone and understand the load-induced stimulation of osteocytes.
This study elaborates on the methods to develop 3D single osteocyte models using confocal microscope images of the Lacunar-Canalicular Network&#160;(LCN) to perform CFD analysis utilizing various computational modeling software. Prior to confocal microscopy, the mouse bones are sectioned and stained with Fluorescein isothiocyanate (FITC) dye to label the LCN. At 100x resolution, Z-stack images are collected using a confocal microscope and imported into MIMICS software (3D image-based processing software) to construct a surface model of the LCN and osteocyte-dendritic processes.
These surfaces are then subtracted using a Boolean operation in 3-Matic software (3D data optimization software) to model the lacunar fluidic space around the osteocyte cell body and canalicular space around the dendrites containing lacunocanalicular fluid. 3D volumetric fluid geometry is imported into ANSYS software (simulation software) for CFD analysis. ANSYS CFX (CFD software) is used to apply physiological loading on the bone as fluid pressure, and the wall shear stresses on the osteocytes and dendritic processes are determined. The morphology of the LCN affects the shear stress values sensed by the osteocyte cell membrane and cell processes. Therefore, the details of how confocal image-based models are developed can be valuable in understanding osteocyte mechanosensation and can lay the groundwork for future studies in this area.

This protocol describes how to develop confocal-derived osteocyte models to investigate the amount of fluid flow shear stress an osteocyte and its dendritic processes are subjected to due to mechanical loading. An aged and a young C57BL6 mouse were selected to build young and aged confocal image-based osteocyte models. Six other simulated osteocyte models were generated from the same young osteocyte model using the geometry modification technique to study the alteration of LCN morphology due to aging or bone disease. Geometry-modified parameters included the lacunar-canalicular thickness, canalicular density, and diameter. Modified osteocyte models were developed by changing one parameter at a time while holding the other parameters the same. Therefore, the effect of each of the morphological parameters on FFSS is examined in isolation. The results revealed that the amount of FFSS experienced by an osteocyte was different in these osteocyte models, showing the dependency of FFSS on LCN morphology. FFSS values were predicted on various parts of the osteocyte model, such as the osteocyte cell body or dendritic membranes. Dendrites in all models experienced greater FFSS values than the osteocyte cell body membrane.
Figure 5 shows the FFSS contours for young vs. aged osteocyte models generated directly from confocal FITC image stacks and six other osteocyte models developed from Model 1 using the geometry modification technique. The average FFSS on osteocyte and dendritic membranes in the young osteocyte model (Model 1) was 0.42 Pa, which is more significant than the average FFSS of 0.13 Pa in the aged osteocyte model (Model 2). The obtained data show an increase in average FFSS in osteocytes with larger lacunar-canalicular space (Model 7 and Model 8) and a decrease in FFSS in osteocytes with lower canalicular density (Model 3 and Model 4), while FFSS did not significantly change with changing dendrite diameter in Model 5 and Model 6. Model 2 and Model 4 had the lowest FFSS values of 0.19 Pa and 0.13 Pa, respectively. The magnitude of FFSS in Model 8 is greater than in the other osteocyte models. Not only are most of the dendrites in red color in Model 8 but the osteocyte cell body also experiences higher FFSS than in other models.
FFSS experienced by osteocytes was also presented based on the lacunar-canalicular space volume and surface area, which were obtained by post-processing in ANSYS CFX to find out the correlation between FFSS and the altered structure of the LCN. It was also possible to measure the surface area of each model exposed to FFSS greater than 0.8 Pa, which is thought to be the level of FFSS that causes bone formation21,22. The osteocyte modeling showed a correlation between the average FFSS and the lacunar-canalicular space volume or surface area for all osteocyte models, as shown in Table 1. The aged osteocyte (Model 2) and simulated middle-aged osteocyte (Model 4) with the lowest canalicular density had the lowest lacunar-canalicular surface area of 1372.7 µm2 and 1122.6 µm2. Model 2 and Model 4 also had the lowest lacunar-canalicular space volume of 229.4 µm3 and 203.5 µm3, respectively, which can be correlated to their low amounts of FFSS. Model 6, with the largest dendrite diameter, had the largest lacunar-canalicular surface area. Model 8 had the largest lacunar-canalicular space volume and the greatest FFSS values.
<img alt="Fluorescence microscopy images at 5x, 20x, 100x magnification; cellular structure analysis." class="xfigimg" src="/files/ftp_upload/64699/64699fig01.jpg" data-altupdatedat="1788304762000" data-alt="Figure 1"> 
Figure 1: Single Z-plane confocal images of a transverse section of a FITC-stained femur taken with 5x, 20x, and 100x oil objectives. The Lacunar-canalicular network is visible in the image taken with the 100x objective. <a href="https://www.jove.com/files/ftp_upload/64699/64699fig01large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Osteocyte morphology diagram, young vs aged comparison, showing cell structure differences." class="xfigimg" src="/files/ftp_upload/64699/64699fig02.jpg" data-altupdatedat="1788304762000" data-alt="Figure 2"> 
Figure 2: Confocal image-based models of a young (Model 1) and an aged (Model 2) osteocyte generated from the FITC-stained femur of 4-month-old and 22-month-old C57BL6 mouse. The magnified image shows the lacunar-canalicular (LC) thickness for young osteocyte and the same for the aged osteocyte. Model 1 image has been modified with permission from Niroobakhsh et al.23. <a href="https://www.jove.com/files/ftp_upload/64699/64699fig02large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Simulated osteocyte models, dendrite diameter and lacuna-canalicular thickness variations, diagram." class="xfigimg" src="/files/ftp_upload/64699/64699fig03.jpg" data-altupdatedat="1788304762000" data-alt="Figure 3"> 
Figure 3: Six osteocyte models developed from Model 1 using the geometry modification technique. Canalicular density was reduced to generate Model 3 (simulated middle-aged) and Model 4 (simulated aged), while dendrite diameter and LC thickness were the same as Model 1. Model 5 and Model 6 were developed with the same canalicular density and LC thickness as Model 1, whereas their dendrite diameter was two and four times greater than Model 1. Model 7 and Model 8 had two times and four times greater lacunar-canalicular thickness, respectively, while canalicular density and dendrite diameter remained the same as in Model 1. Model 1 image has been modified with permission from Niroobakhsh et al.23. <a href="https://www.jove.com/files/ftp_upload/64699/64699fig03large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Cell structure diagram; membrane visualization; microscopy image; cellular biology model." class="xfigimg" src="/files/ftp_upload/64699/64699fig04.jpg" data-altupdatedat="1788304762000" data-alt="Figure 4"> 
Figure 4: 3D view confocal image-based young osteocyte model. 3D view of the cross-sectional area of the confocal image-based young osteocyte model depicting the modeling of the lacunar-canalicular space, which is used for fluid analysis around the osteocyte and dendrites. <a href="https://www.jove.com/files/ftp_upload/64699/64699fig04large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Osteocyte models diagram; confocal imaging, geometry-modified simulation; stress distribution." class="xfigimg" src="/files/ftp_upload/64699/64699fig05.jpg" data-altupdatedat="1788304762000" data-alt="Figure 5"> 
Figure 5: FFSS distribution for all osteocyte models shows the dependence of FFSS magnitude on osteocyte morphology. Red regions show the highest FFSS values, especially on the dendritic membranes rather than the osteocyte cell body. This image has been modified with permission from Niroobakhsh et al.23. <a href="https://www.jove.com/files/ftp_upload/64699/64699fig05large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td></td>
			<td>Models</td>
			<td>Lacunar-canalicular surface area (µm2)</td>
			<td>Lacunar-canalicular space volume (µm3)</td>
			<td>Ave FFSS (Pa)</td>
		</tr>
		<tr>
			<td rowspan="4">Confocal image-based osteocyte models</td>
			<td>Model 1</td>
			<td rowspan="2">2652.9</td>
			<td rowspan="2">261</td>
			<td rowspan="2">0.42</td>
		</tr>
		<tr>
			<td>Young (4-month-old)</td>
		</tr>
		<tr>
			<td>Model 2</td>
			<td rowspan="2">1372.7</td>
			<td rowspan="2">229.4</td>
			<td rowspan="2">0.19</td>
		</tr>
		<tr>
			<td>Aged (22-month-old)</td>
		</tr>
		<tr>
			<td rowspan="12">Simulated osteocyte models generated using geometry modification technique</td>
			<td>Model 3</td>
			<td rowspan="2">1530.4</td>
			<td rowspan="2">242.6</td>
			<td rowspan="2">0.27</td>
		</tr>
		<tr>
			<td>(Simulated middle-aged)</td>
		</tr>
		<tr>
			<td>Model 4</td>
			<td rowspan="2">1122.6</td>
			<td rowspan="2">203.5</td>
			<td rowspan="2">0.13</td>
		</tr>
		<tr>
			<td>(Simulated aged)</td>
		</tr>
		<tr>
			<td>Model 5</td>
			<td rowspan="2">3724.9</td>
			<td rowspan="2">302.9</td>
			<td rowspan="2">0.49</td>
		</tr>
		<tr>
			<td>(2 x dendrite diameter)</td>
		</tr>
		<tr>
			<td>Model 6</td>
			<td rowspan="2">4791</td>
			<td rowspan="2">337</td>
			<td rowspan="2">0.57</td>
		</tr>
		<tr>
			<td>(4 x dendrite diameter)</td>
		</tr>
		<tr>
			<td>Model 7</td>
			<td rowspan="2">3116.3</td>
			<td rowspan="2">440.6</td>
			<td rowspan="2">1.17</td>
		</tr>
		<tr>
			<td>(2 x LC thickness)</td>
		</tr>
		<tr>
			<td>Model 8</td>
			<td rowspan="2">3542.2</td>
			<td rowspan="2">658.2</td>
			<td rowspan="2">1.89</td>
		</tr>
		<tr>
			<td>(4 x LC thickness)</td>
		</tr>
	</tbody>
</table>
Table 1: Comparison of osteocyte models. Comparison of average FFSS, lacunar-canalicular space volume, and surface area in eight distinct osteocyte models, including two confocal image-based models and six simulated models generated from Model 1 using the geometry modification technique.

This study outlines the method to visualize and develop three-dimensional (3D) models of osteocytes within the lacunar-canalicular network (LCN) for computational fluid dynamics (CFD) analysis. The generated models using this method help to understand osteocyte mechanosensation in healthy or diseased bones.

Osteocytes are the bone cells that are thought to respond to mechanical strains and fluid flow shear stress (FFSS) by activating various biological ...

Research

JoVE Journal

Bioengineering

Treatment of primary mouse osteoblasts with retinoic acid produces a homogeneous population of ramified cells bearing morphological and molecular features of osteocytes. The method overcomes the difficulty of obtaining and maintaining primary osteocytes in culture, and can be advantageous to study cells derived from transgenic...

Application of Retinoic Acid to Obtain Osteocytes Cultures from Primary Mouse Osteoblasts

Osteoclasts are the principal bone-resorbing cell in the body. An ability to isolate osteoclasts in large numbers has resulted in significant advances in the understanding of osteoclast biology. In this protocol, we describe a method for isolation, cultivating and quantifying osteoclast activity in vitro.

Osteoclast Derivation from Mouse Bone Marrow

We describe a protocol for hybrid imaging, combining fluorescence-mediated tomography (FMT) with micro computed tomography (µCT). After fusion and reconstruction, we perform interactive organ segmentation to extract quantitative measurements of the fluorescence distribution.

Hybrid &#181;CT-FMT imaging and image analysis

Implantation of autologous and allogeneic bone grafts constitute accepted approaches to treat major craniofacial bone loss. Yet the effect of graft composition on the interplay among neovascularization, cell differentiation and bone formation is unclear. We present a multimodal imaging protocol aimed to elucidate the angiogenesis-osteogenesis interdependence at the graft proximity.

Computed Tomography and Optical Imaging of Osteogenesis-angiogenesis Coupling to Assess Integration of Cranial Bone Autografts and Allografts

Uranium is known to affect bone metabolism. Here, we present a protocol aimed at investigating the effect of natural uranium exposure on the viability, the differentiation, and the function of osteoclasts, the cells in charge of bone resorption.

Methods for Analyzing the Impacts of Natural Uranium on In Vitro Osteoclastogenesis

We present a protocol to compare the state of minerals in vesicles released by two human bone cell lines: hFOB 1.19 and Saos-2. Their mineralization profiles were analyzed by Alizarin Red-S (AR-S) staining, ultraviolet (UV) light visualization, transmission electron microscopy (TEM) imaging and energy dispersive X-ray microanalysis (EDX).

Analysis of Minerals Produced by hFOB 1.19 and Saos-2 Cells Using Transmission Electron Microscopy with Energy Dispersive X-ray Microanalysis

We employed a geological (coring) sampling protocol to procure cortical bone specimens of uniform size for SRµCT experiments from the anterior aspect of human femora. This method is minimally destructive, efficient, results in cylindrical specimens that minimize imaging artifacts from irregular sample shapes and improves microarchitectural visualization and analysis.

A Sectioning, Coring, and Image Processing Guide for High-Throughput Cortical Bone Sample Procurement and Analysis for Synchrotron Micro-CT

This protocol describes the dissection of neonatal dmp1-topaz mouse calvaria and isolation of osteocytes expressing the green fluorescent protein through cell digestion and fractionation, in addition to osteocyte preparation for fluorescence activated cell sorting (FACS).

Obtaining Primary Osteocytes Through Murine Calvarial Fractionation of GFP-Expressing Osteocytes

The current article describes performing an intravital imaging approach to observe mechanically induced calcium signaling of embedded osteocytes in vivo in real-time in response to tissue-level mechanical loading of the mouse third metatarsal.

A Fluorescent Intravital Imaging Approach to Study Load-Induced Calcium Signaling Dynamics in Mouse Osteocytes

This protocol presents an economical and efficient method for quantitative evaluation of bone microarchitecture in a mouse model of osteoporosis by combining Hematoxylin-Eosin (HE) staining and micro-computed tomography (Micro-CT) techniques.

Trabecular Bone Microarchitecture Evaluation in an Osteoporosis Mouse Model

We propose a protocol for assessing bone regeneration at various hierarchical levels. This multimodal approach allows for structural analysis across different length scales, effectively addressing challenges in sample preparation. It also provides a reliable workflow for evaluating the effectiveness of biomaterials in promoting bone regeneration.

Multimodal Approach to Assess Bone Regeneration and Scaffold Performance

Learn computational fluid dynamics analysis using confocal microscopy and finite element modeling to evaluate fluid flow shear stress on osteocyte dendritic membranes in bone mechanosensation research.

Analysis and Imaging of Osteocytes

Explore More Videos