Automated high-resolution three-dimensional fluorescence imaging of large biological specimens

Spatial relationship between bone formation and mechanical stimulus within cortical bone: Combining 3D fluorochrome mapping and poroelastic finite element modelling

Bone is a dynamic tissue and adapts its architecture in response to biological and mechanical factors. Here we investigate how cortical bone formation is spatially controlled by the local mechanical environment in the murine tibia axial loading model (C57BL/6). We obtained 3D locations of new bone formation by performing ‘slice and view’ 3D fluorochrome mapping of the entire bone and compared these sites with the regions of high fluid velocity or strain energy density estimated using a finite element model, validated with ex-vivo bone surface strain map acquired ex-vivo using digital image correlation. For the comparison, 2D maps of the average bone formation and peak mechanical stimulus on the tibial endosteal and periosteal surface across the entire cortical surface were created. Results showed that bone formed on the periosteal and endosteal surface in regions of high fluid flow. Peak strain energy density predicted only the formation of bone periosteally. Understanding how the mechanical stimuli spatially relates with regions of cortical bone formation in response to loading will eventually guide loading regime therapies to maintain or restore bone mass in specific sites in skeletal pathologies.

•

3D spatial representation of new bone formation after loading is shown by fluorochrome mapping of the entire mouse tibia

•

Regions of new bone formation spatially associate with regions of high strain and fluid mechanical stimulus in a FE model

•

The FE model was validated with the strains on the bone surface determined ex-vivo using digital image correlation

•

Regions of new bone formation co-localize in sites of peak fluid flow, both endosteally and periosteally

•

Peak strain energy density was able to predict only periosteal bone formation

Application of autofluorescence robotic histology for quantitative evaluation of the 3-dimensional morphology of murine articular cartilage

Murine models of osteoarthritis (OA) are increasingly important for understating pathogenesis and for testing new therapeutic approaches. Their translational potential is, however, limited by the reduced size of mouse limbs which requires a much higher resolution to evaluate their articular cartilage compared to clinical imaging tools. In experimental models, this tissue has been predominantly assessed by time-consuming histopathology using standardized semi-quantitative scoring systems. This study aimed to develop a novel imaging method for 3-dimensional (3D) histology of mouse articular cartilage, using a robotic system-termed here "3D histocutter"-which automatically sections tissue samples and serially acquires fluorescence microscopy images of each section. Tibiae dissected from C57Bl/6 mice, either naïve or OA-induced by surgical destabilization of the medial meniscus (DMM), were imaged using the 3D histocutter by exploiting tissue autofluorescence. Accuracy of 3D imaging was validated by ex vivo contrast-enhanced micro-CT and sensitivity to lesion detection compared with conventional histology. Reconstructions of tibiae obtained from 3D histocutter serial sections showed an excellent agreement with contrast-enhanced micro-CT reconstructions. Furthermore, osteoarthritic features, including articular cartilage loss and osteophytes, were also visualized. An in-house developed software allowed to automatically evaluate articular cartilage morphology, eliminating the subjectivity associated to semi-quantitative scoring and considerably increasing analysis throughput. The novelty of this methodology is, not only the increased throughput in imaging and evaluating mouse articular cartilage morphology starting from conventionally embedded samples, but also the ability to add the third dimension to conventional histomorphometry which might be useful to improve disease assessment in the model.

Multiscale imaging of bone microdamage

Bone is a structural and hierarchical composite that exhibits remarkable ability to sustain complex mechanical loading and resist fracture. Bone quality encompasses various attributes of bone matrix from the quality of its material components (type-I collagen, mineral and non-collagenous matrix proteins) and cancellous microarchitecture, to the nature and extent of bone microdamage. Microdamage, produced during loading, manifests in multiple forms across the scales of hierarchy in bone and functions to dissipate energy and avert fracture. Microdamage formation is a key determinant of bone quality, and through a range of biological and physical mechanisms, accumulates with age and disease. Accumulated microdamage in bone decreases bone strength and increases bone's propensity to fracture. Thus, a thorough assessment of microdamage, across the hierarchical levels of bone, is crucial to better understand bone quality and bone fracture. This review article details multiple imaging modalities that have been used to study and characterize microdamage; from bulk staining techniques originally developed by Harold Frost to assess linear microcracks, to atomic force microscopy, a modality that revealed mechanistic insights into the formation diffuse damage at the ultrastructural level in bone. New automated techniques using imaging modalities, such as microcomputed tomography are also presented for a comprehensive overview.

Three-dimensional characterization of resorption cavity size and location in human vertebral trabecular bone

The number and size of resorption cavities in cancellous bone are believed to influence rates of bone loss, local tissue stress and strain and potentially whole bone strength. Traditional two-dimensional approaches to measuring resorption cavities in cancellous bone report the percent of the bone surface covered by cavities or osteoclasts, but cannot measure cavity number or size. Here we use three-dimensional imaging (voxel size 0.7×0.7×5.0 μm) to characterize resorption cavity location, number and size in human vertebral cancellous bone from nine elderly donors (7 male, 2 female, ages 47-80 years). Cavities were 30.10 ± 8.56 μm in maximum depth, 80.60 ± 22.23∗10(3) μm(2) in surface area and 614.16 ± 311.93∗10(3) μm(3) in volume (mean ± SD). The average number of cavities per unit tissue volume (N.Cv/TV) was 1.25 ± 0.77 mm(-3). The ratio of maximum cavity depth to local trabecular thickness was 30.46 ± 7.03% and maximum cavity depth was greater on thicker trabeculae (p<0.05, r(2)=0.14). Half of the resorption cavities were located entirely on nodes (the intersection of two or more trabeculae) within the trabecular structure. Cavities that were not entirely on nodes were predominately on plate-like trabeculae oriented in the cranial-caudal (longitudinal) direction. Cavities on plate-like trabeculae were larger in maximum cavity depth, cavity surface area and cavity volume than cavities on rod-like trabeculae (p<0.05). We conclude from these findings that cavity size and location are related to local trabecular microarchitecture.

Three-dimensional dynamic bone histomorphometry

Dynamic bone histomorphometry is the standard method for measuring bone remodeling at the level of individual events. Although dynamic bone histomorphometry is an invaluable tool for understanding osteoporosis and other metabolic bone diseases, the technique's two-dimensional nature requires the use of stereology and prevents measures of individual remodeling event number and size. Here, we used a novel three-dimensional fluorescence imaging technique to achieve measures of individual resorption cavities and formation events. We performed this three-dimensional histomorphometry approach using a common model of postmenopausal osteoporosis, the ovariectomized rat. The three-dimensional images demonstrated the spatial relationship between resorption cavities and formation events consistent with the hemiosteonal model of cancellous bone remodeling. Established ovariectomy was associated with significant increases in the number of resorption cavities per unit bone surface (2.38 ± 0.24 mm⁻² sham surgery versus 3.86 ± 0.35 mm⁻² bilateral ovariectomy [OVX], mean ± SD, p < 0.05) and total volume occupied by cavities per unit bone volume (0.38% ± 0.06% sham versus 1.12% ± 0.18% OVX, p < 0.001), but there was no difference in surface area per resorption cavity, maximum cavity depth, or cavity volume. In addition, we found that established ovariectomy is associated with increased size of bone formation events because of the merging of formation events (23,700 ± 6,890 µm² sham verusus 33,300 ± 7,950 µm² OVX). No differences in mineral apposition rate (determined in 3D) were associated with established ovariectomy. That established estrogen depletion is associated with increased number of remodeling events with only subtle changes in remodeling event size suggests that circulating estrogens may have their primary effect on the origination of new basic multicellular units with relatively little effect on the progression and termination of active remodeling events.

Three-dimensional surface texture visualization of bone tissue through epifluorescence-based serial block face imaging

Serial block face imaging is a microscopy technique in which the top of a specimen is cut or ground away and a mosaic of images is collected of the newly revealed cross-section. Images collected from each slice are then digitally stacked to achieve 3D images. The development of fully automated image acquisition devices has made serial block face imaging more attractive by greatly reducing labour requirements. The technique is particularly attractive for studies of biological activity within cancellous bone as it has the capability of achieving direct, automated measures of biological and morphological traits and their associations with one another. When used with fluorescence microscopy, serial block face imaging has the potential to achieve 3D images of tissue as well as fluorescent markers of biological activity. Epifluorescence-based serial block face imaging presents a number of unique challenges for visualizing bone specimens due to noise generated by sub-surface signal and local variations in tissue autofluorescence. Here we present techniques for processing serial block face images of trabecular bone using a combination of non-uniform illumination correction, precise tiling of the mosaic in each cross-section, cross-section alignment for vertical stacking, removal of sub-surface signal and segmentation. The resulting techniques allow examination of bone surface texture that will enable 3D quantitative measures of biological processes in cancellous bone biopsies.

Voxel size and measures of individual resorption cavities in three-dimensional images of cancellous bone

Cavities formed by osteoclasts on the surface of cancellous bone during bone remodeling (resorption cavities) are believed to act as stress risers and impair cancellous bone strength and stiffness. Although resorption cavities are readily detected as eroded surfaces in histology sections, identification of resorption cavities in three-dimensional images of cancellous bone has been rare. Here we use sub-micrometer resolution images of rat lumbar vertebral cancellous bone obtained through serial milling (n=5) to determine how measures of the number and surface area of resorption cavities are influenced by image resolution. Three-dimensional images of a 1 mm cube of cancellous bone were collected at 0.7x0.7x5.0 microm/voxel using fluorescence based serial milling and uniformly coarsened to four other resolutions ranging from 1.4x1.4x5.0 to 11.2x11.2x10 microm/voxel. Cavities were identified in the three-dimensional image as an indentation on the cancellous bone surface and were confirmed as eroded surfaces by viewing two-dimensional cross-sections (mimicking histology techniques). The number of cavities observed in the 0.7x0.7x5.0 microm/voxel images (22.0+/-1.43, mean+/-SD) was not significantly different from that in the 1.4x1.4x5.0 microm/voxel images (19.2+/-2.59) and an average of 79% of the cavities observed at both of these resolutions were coincident. However, at lower resolutions, cavity detection was confounded by low sensitivity (<20%) and high false positive rates (>40%). Our results demonstrate that when image voxel size exceeds 1.4x1.4x5.0 microm/voxel identification of resorption cavities by bone surface morphology is highly inaccurate. Experimental and computational studies of resorption cavities in three-dimensional images of cancellous bone may therefore require images to be collected at resolutions of 1.4 microm/pixel in-plane or better to ensure consistent identification of resorption cavities.

Validity of serial milling-based imaging system for microdamage quantification

Understanding the three-dimensional distribution of microdamage within trabecular bone may help provide a better understanding of the mechanisms of bone failure. Toward that end, a novel serial milling-based fluorescent imaging system was developed for quantifying microscopic damage in three dimensions throughout cores of trabecular bone. The overall goal for this study was to compare two-dimensional (2D), surface-based measures of microdamage extracted from this new imaging system against those from more conventional histological section analyses. Human vertebral trabecular cores were isolated, stained en bloc with a series of chelating fluorochromes, monotonically loaded, and underwent microdamage quantification via the two methods. Bone area fraction measured by the new system was significantly correlated to that measured by histological point counting (p<0.001, R(2)=0.80). Additionally, the new system produced statistically equivalent (p=0.021) measures of damage fraction (mean+/-SD), Dx.AF=0.047+/-0.021, to that obtained from stereological point counting, Dx.AF=0.048+/-0.017, at a 10% difference level. These results demonstrate that this serial milling-based fluorescent imaging system provides a destructive yet practical alternative to more conventional histologic section analysis in addition to its ability to provide a better understanding of the three-dimensional nature of microdamage.