Analysis and Imaging of Osteocytes

Summary

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.

Abstract

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 (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 (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.

Introduction

Osteocytes are postulated to regulate bone mass in response to physical exercise¹. 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 signaling^2,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 osteocytes^5,6. These changes in bone architecture cause osteocytes to experience different levels of FFSS and strains^7,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 LCN^9,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 bone^13,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 3D^16,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 around¹⁸.

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 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 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 bone^5,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.

Protocol

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

1. 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 °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
2. Embed bones quickly in a fast polymerizing acrylic (Table of Materials) following the manufacturer'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.
3. Cut thick 300 µm thick transverse slices from a standardized site above the third trochanter using a diamond saw and store them at 4 °C in 70% ethanol prior to FITC staining.
4. Polish the sections using 600, 800, and then 1200 grit sandpaper to a final thickness of ~90-100 μm.
   NOTE: Attaining proper thickness was ensured using a digital caliper.
5. Rinse the sections in 70%, 95%, and 100% ethanol for 5 min each.
6. Stain them in 1% FITC in 100% ethanol for 4 h at room temperature (RT) in the dark with moderate shaking.
7. Wash the sections in 100% ethanol for 30 min with gentle shaking in the dark. Then, air-dry them overnight in the dark.
8. 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.

2. Confocal microscopy

1. Use a confocal microscope for imaging of FITC-stained bone slices.
2. Use a 100x 1.44NA oil objective with a digital zoom of 1.7 and a step size of 0.126 μm to collect detailed Z-stacks of 400 Z-planes at 1024 x 1024 pixels, 0.089 μm pixel resolution.
3. 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.
4. 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.
5. Use the 100x Z-stacks for the computer modeling of osteocytes.

3. Computer modeling

1. Import the collected 100x images, in TIFF format, to ImageJ software to build a sequence of images of the LCN in the Z-direction.
2. Import the Z-stacks into the 3D image-based processing software to construct a mask of the LCN after defining the image orientation.
3. 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.
4. 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 µm, 14 µm, and 19 µm.
5. 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.
6. Convert the lacunar-canalicular mask into an object using the Calculated Part operation in the 3D image-based processing software.
7. 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 µm and 0.08 µm, respectively^9,18.
8. Export the objects (STL format) as the last step in the 3D image-based processing software.
9. Import two layers of the LCN and osteocyte-dendritic membranes into the 3D data optimization software to generate a volume mesh.
10. 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.
11. Remove the inverted normal parts, intersecting triangles, and bad contours using the Auto-Fix operation in Fix Wizard.
12. Replace the overlapping triangles manually by defining new ones or automatically via the Fill Hole Normal operation.
13. Improve the mesh quality using operations including filtering sharp triangles, small edges, and small shells.
14. 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.
15. 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.

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.

1. Choose the young osteocyte as the base model and modify it to build other distinct osteocyte models by applying morphological alterations.
2. Generate an osteocyte model with different canalicular densities from the base model by changing the image thresholding in the 3D image-based processing software.
   1. 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.
3. 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.

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.

1. Create a Fluid Flow in the simulation software to prepare the models for the CFD analysis.
2. 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.
3. 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.
4. 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.
5. Click on mesh and select linear tetrahedral elements using an element size of 0.06 µ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.
6. 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.
7. Export the mesh (fluent file format) as it loads faster in the setup in the next step.
8. 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.
9. To imitate physiological conditions, exert a fluid inlet pressure of 300 Pa and 0 Pa on the inlets and outlets, respectively^19,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.
10. Treat the interstitial laminar fluid as water⁹, 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 fluid⁹.
11. 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.

6. CFD post processing

1. 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.
2. 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.
3. 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.

Results

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. Geo...

Discussion

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 a...

Materials

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