Quantification of the bone lacunocanalicular network from 3D X-ray phase nanotomography images

There is a growing interest in developing 3D microscopy for the exploration of thick biological tissues. Recently, 3D X-ray nanocomputerised tomography has proven to be a suitable technique for imaging the bone lacunocanalicular network. This interconnected structure is hosting the osteocytes which play a major role in maintaining bone quality through remodelling processes. 3D images have the potential to reveal the architecture of cellular networks, but their quantitative analysis remains a challenge due to the density and complexity of nanometre sized structures and the need to handle and process large datasets, for example, 2048³ voxels corresponding to 32 GB per individual image in our case. In this work, we propose an efficient image processing approach for the segmentation of the network and the extraction of characteristic parameters describing the 3D structure. These parameters include the density of lacunae, the porosity of lacunae and canaliculi, and morphological features of lacunae (volume, surface area, lengths, anisotropy etc.). We also introduce additional parameters describing the local environment of each lacuna and its canaliculi. The method is applied to analyse eight human femoral cortical bone samples imaged by magnified X-ray phase nanotomography with a voxel size of 120 nm, which was found to be a good compromise to resolve canaliculi while keeping a sufficiently large field of view of 246 μm in 3D. The analysis was performed on a total of 2077 lacunae showing an average length, width and depth of 17.1 μm × 9.2 μm × 4.4 μm, with an average number of 58.2 canaliculi per lacuna and a total lacuno-canalicular porosity of 1.12%. The reported descriptive parameters provide information on the 3D organisation of the lacuno-canalicular network in human bones.