Detection of trabecular bone microdamage by micro-computed tomography

Microdamage is an important component of bone quality and affects bone remodeling. Improved techniques to assess microdamage without the need for histological sectioning would provide insight into the role of microdamage in trabecular bone strength by allowing the spatial distribution of damage within the trabecular microstructure to be measured. Nineteen cylindrical trabecular bone specimens were prepared and assigned to two groups. The specimens in group I were damaged to 3% compressive strain and labeled with BaSO(4). Group II was not loaded, but was labeled with BaSO(4). Micro-computed tomography (Micro-CT) images of the specimens were obtained at 10 microm resolution. The median intensity of the treated bone tissue was compared between groups. Thresholding was also used to measure the damaged area fraction in the micro-CT scans. The histologically measured damaged area fraction, the median CT intensity, and the micro-CT measured damaged area fraction were all higher in the loaded group than in the unloaded group, indicating that the micro-CT images could differentiate the damaged specimen group from the unloaded specimens. The histologically measured damaged area fraction was positively correlated with the micro-CT measured damaged area fraction and with the median CT intensity of the bone, indicating that the micro-CT images can detect microdamage in trabecular bone with sufficient accuracy to differentiate damage levels between samples. This technique provides a means to non-invasively assess the three-dimensional distribution of microdamage within trabecular bone test specimens and could be used to gain insight into the role of trabecular architecture in microdamage formation.     

Three-dimensional analysis of nonhuman primate trabecular architecture using micro-computed tomography

Until recently, detailed analyses of the architecture of nonhuman primate cancellous bone have not been possible due to a combination of methodological constraints, including poor resolution imaging or destructive protocols. The development of micro-computed tomography (microCT) and morphometric methods associated with this imaging modality offers anthropologists a new means to study the comparative architecture of cancellous bone. Specifically, microCT will allow anthropologists to investigate the relationship between locomotor behavior and trabecular structure. We conducted a preliminary study on the trabecular patterns in the proximal humerus and femur of Hylobates lar, Ateles paniscus, Macaca mulatta, and Papio anubis to investigate the quantitative differences in their trabecular architecture and evaluate the potential of microCT in anthropological inquiry. MicroCT allows the researcher to evaluate variables beyond simple two-dimensional orientations and radiographic densities. For example, this methodology facilitates the study of trabecular thickness and bone volume fraction using three-dimensional data. Results suggest that density-related parameters do not reliably differentiate suspensory-climbing species from quadrupedal species. However, preliminary results indicate that measurements of the degree of anisotropy, a measure of trabecular orientation uniformity, do distinguish suspensory-climbing taxa from more quadrupedal species. The microCT method is an advance over conventional radiography and medical CT because it can accurately resolve micron-sized struts that make up cancellous bone, and from these images a wide array of parameters that have been demonstrated to be related to cancellous bone mechanical properties can be measured. Methodological problems pertinent to any comparative microCT study of primate trabecular architecture are discussed.     

Determination of vertebral endplate deformation under load using micro-computed tomography

Endplate strength plays an important role in preventing vertebral failure of normal vertebrae and in cases where surgical intervention has replaced the disc with an implant or has altered the vertebral loading. We have developed a non-contact method based on the principles of image guided failure analysis, mechanical testing, and micro-computed tomography analysis, which allows for in vitro quantification of endplate deformation under axial load. The method allows for the implementation of a repeated measures experimental design, each specimen acting as its own control. Our methodology was validated using cadaveric functional spine units, loaded stepwise from 200 N to a maximum of 2000 N. The loading protocol was repeated over two days, allowing time for recovery of the disc mechanical properties. We found no meaningful difference in measured force, stiffness, and endplate deformation between day 1 and day 2. The mean fiducial registration error was less than 0.015 mm for all three axes. Endplate deformation could be reproducibly estimated. The root mean squared error was 0.03 mm, which is the effective precision of the method. Using this micro-CT based method, the effect of interbody implants, grafts, disc replacement strategies, and surgical procedures such as nucleotomy and vertebral cement augmentation on endplate mechanical behaviour can be ascertained.     

Quantification of lung fibrosis and emphysema in mice using automated micro-computed tomography

Background

In vivo high-resolution micro-computed tomography allows for longitudinal image-based measurements in animal models of lung disease. The combination of repetitive high resolution imaging with fully automated quantitative image analysis in mouse models of lung fibrosis lung benefits preclinical research. This study aimed to develop and validate such an automated micro-computed tomography analysis algorithm for quantification of aerated lung volume in mice; an indicator of pulmonary fibrosis and emphysema severity.

Methodology

Mice received an intratracheal instillation of bleomycin (n = 8), elastase (0.25U elastase n = 9, 0.5U elastase n = 8) or saline control (n = 6 for fibrosis, n = 5 for emphysema). A subset of mice was scanned without intervention, to evaluate potential radiation-induced toxicity (n = 4). Some bleomycin-instilled mice were treated with imatinib for proof of concept (n = 8). Mice were scanned weekly, until four weeks after induction, when they underwent pulmonary function testing, lung histology and collagen quantification. Aerated lung volumes were calculated with our automated algorithm.

Principal Findings

Our automated image-based aerated lung volume quantification method is reproducible with low intra-subject variability. Bleomycin-treated mice had significantly lower scan-derived aerated lung volumes, compared to controls. Aerated lung volume correlated with the histopathological fibrosis score and total lung collagen content. Inversely, a dose-dependent increase in lung volume was observed in elastase-treated mice. Serial scanning of individual mice is feasible and visualized dynamic disease progression. No radiation-induced toxicity was observed. Three-dimensional images provided critical topographical information.

Conclusions

We report on a high resolution in vivo micro-computed tomography image analysis algorithm that runs fully automated and allows quantification of aerated lung volume in mice. This method is reproducible with low inherent measurement variability. We show that it is a reliable quantitative tool to investigate experimental lung fibrosis and emphysema in mice. Its non-invasive nature has the unique benefit to allow dynamic 4D evaluation of disease processes and therapeutic interventions.     

Visualization and quantification of whole rat heart laminar structure using high-spatial resolution contrast-enhanced MRI

It has been shown by histology that cardiac myocytes are organized into laminae and this structure is important in function, both influencing the spread of electrical activation and enabling myocardial thickening in systole by laminar sliding. We have carried out high-spatial resolution three-dimensional MRI of the ventricular myolaminae of the entire volume of the isolated rat heart after contrast perfusion [dimeglumine gadopentate (Gd-DTPA)]. Four ex vivo rat hearts were perfused with Gd-DTPA and fixative and high-spatial resolution MRI was performed on a 9.4T MRI system. After MRI, cryosectioning followed by histology was performed. Images from MRI and histology were aligned, described, and quantitatively compared. In the three-dimensional MR images we directly show the presence of laminae and demonstrate that these are highly branching and are absent from much of the subepicardium. We visualized these MRI volumes to demonstrate laminar architecture and quantitatively demonstrated that the structural features observed are similar to those imaged in histology. We showed qualitatively and quantitatively that laminar architecture is similar in the four hearts. MRI can be used to image the laminar architecture of ex vivo hearts in three dimensions, and the images produced are qualitatively and quantitatively comparable with histology. We have demonstrated in the rat that: 1) laminar architecture is consistent between hearts; 2) myolaminae are absent from much of the subepicardium; and 3) although localized orthotropy is present throughout the myocardium, tracked myolaminae are branching structures and do not have a discrete identity.     

In vivo characterization of lung morphology and function in anesthetized free-breathing mice using micro-computed tomography.

Lung morphology and function in human subjects can be monitored with computed tomography (CT). Because many human respiratory diseases are routinely modeled in rodents, a means of monitoring the changes in the structure and function of the rodent lung is desired. High-resolution images of the rodent lung can be attained with specialized micro-CT equipment, which provides a means of monitoring rodent models of lung disease noninvasively with a clinically relevant method. Previous studies have shown respiratory-gated images of intubated and respirated mice. Although the image quality and resolution are sufficient in these studies to make quantitative measurements, these measurements of lung structure will depend on the settings of the ventilator and not on the respiratory mechanics of the individual animals. In addition, intubation and ventilation can have unnatural effects on the respiratory dynamics of the animal, because the airway pressure, tidal volume, and respiratory rate are selected by the operator. In these experiments, important information about the symptoms of the respiratory disease being studied may be missed because the respiration is forced to conform to the ventilator settings. In this study, we implement a method of respiratory-gated micro-CT for use with anesthetized free-breathing rodents. From the micro-CT images, quantitative analysis of the structure of the lungs of healthy unconscious mice was performed to obtain airway diameters, lung and airway volumes, and CT densities at end expiration and during inspiration. Because the animals were free breathing, we were able to calculate tidal volume (0.09 +/- 0.03 ml) and functional residual capacity (0.16 +/- 0.03 ml).     

Computed tomography and automated image analysis of prehistoric femora

Non-invasive characterization of limb bone cross-sectional geometry would be useful for biomechanical analyses of skeletal collections. Computed tomography (CT) is potentially the method of choice. Additionally, CT images are suitable for automated analysis. CT is here shown to be both accurate and precise in the analysis of cross-sectional geometry of prehistoric femora. Beam hardening artifacts can be reduced by using a water bath. As the availability of CT for research increases, both bone density and geometry could be determined simultaneously with this method.     

Quantification of porosity in Acropora pulchra (Brook 1891) using X-ray micro-computed tomography techniques

Skeletal density and porosity characteristics are key parameters for investigations into scleractinian coral growth and for assessing the effects of a range of anthropogenic influences on coral reefs. Typically, skeletal density is measured by using planar X-rays of thin slabs cut from mound-shaped colonies or, for branching forms, by using methods based on Archimedean principles. This paper describes a novel non-destructive technique based on micro-computed tomography (micro-CT) to measure porosity of branching coral skeleton. This approach incorporates methods for segmenting internal and external portions of branch and for distinguishing between skeleton and air, whilst accounting for the effects of beam hardening. Measurements were obtained from colonies of branching Acropora pulchra collected across a reef-flat transect at King Reef, central nearshore Great Barrier Reef. The results show significant variation in porosity within and among branches sampled from individual colonies, but not within a reef-flat transect. Micro-CT techniques yield comparable results to traditional methods based on Archimedean principles, but offer advantages in their suitability for a wider range of coral specimens because of the non-destructive nature of the technique and in their more rigorous control of model parameters that can bias results.     

Modeling the onset and propagation of trabecular bone microdamage during low-cycle fatigue

Relatively small amounts of microdamage have been suggested to have a major effect on the mechanical properties of bone. A significant reduction in mechanical properties (e.g. modulus) can occur even before the appearance of microcracks. This study uses a novel non-linear microdamaging finite-element (FE) algorithm to simulate the low-cycle fatigue behavior of high-density trabecular bone. We aimed to investigate if diffuse microdamage accumulation and concomitant modulus reduction, without the need for complete trabecular strut fracture, may be an underlining mechanism for low-cycle fatigue failure (defined as a 30% reduction in apparent modulus). A microCT constructed FE model was subjected to a single cycle monotonic compression test, and constant and variable amplitude loading scenarios to study the initiation and accumulation of low-cycle fatigue microdamage. Microcrack initiation was simulated using four damage criteria: 30%, 40%, 50% and 60% reduction in bone element modulus (el-MR). Evaluation of structural (apparent) damage using the four different tissue level damage criteria resulted in specimen fatigue failure at 72, 316, 969 and 1518 cycles for the 30%, 40%, 50% and 60% el-MR models, respectively. Simulations based on the 50% el-MR model were consistent with previously published experimental findings. A strong, significant non-linear, power law relationship was found between cycles to failure (N) and effective strain (Deltasigma/E(0)): N=1.394x10(-25)(Deltasigma/E(0))(-12.17), r(2)=0.97, p<0.0001. The results suggest that microdamage and microcrack propagation, without the need for complete trabecular strut fracture, are mechanisms for high-density trabecular bone failure. Furthermore, the model is consistent with previous numerical fatigue simulations indicating that microdamage to a small number of trabeculae results in relatively large specimen modulus reductions and rapid failure.     

Detection of microspheres in vivo using multispectral optoacoustic tomography

We introduce a new approach to detect individual microparticles that contain NIR fluorescent dye by multispectral optoacoustic tomography in the context of the hemoglobin-rich environment within murine liver. We encapsulated a near infrared (NIR) fluorescent dye within polystyrene microspheres, then injected them into the ileocolic vein, which drains to the liver. NIR absorption was determined using multispectral optoacoustic tomography. To quantitate the minimum diameter of microspheres, we used both colorimetric and spatial information to segment the regions in which the microspheres appear. Regional diameter was estimated by doubling the maximum regional distance. We found that the minimum microsphere size threshold for detection by multispectral optoacoustic tomography images is 78.9 µm.     

Quantitative measures of femoral fracture repair in rats derived by micro-computed tomography

Although fracture healing is frequently studied in pre-clinical models of long bone fractures using rodents, there is a dearth of objective quantitative techniques to assess successful healing. Biomechanical testing is possibly the most quantitative and relevant to a successful clinical outcome, but it is a destructive technique providing little insight into the cellular mechanisms associated with healing. The advent of X-ray computed tomography (CT) has provided the opportunity to quantitatively and non-destructively assess bone structure and density, but it is unknown how measurements derived using this technology relate to successful healing. To examine possible relationships, we used a pre-clinical model to test for statistically significant correlations between quantitative characteristics of the callus by micro-CT (μCT) and the bending strength, stiffness, and energy-to-failure of the callus as assessed by three-point bending of excised bones. A closed, transverse fracture was generated in the mid-shaft of rat femurs by impact loading. Shortly thereafter, the rats received a one-time, local injection of either the vehicle or one of four doses of lovastatin. Following sacrifice after 4 weeks of healing, fractured femurs were extracted for μCT analysis and then three-point bending. Setting the region of interest to be 3.2 mm above and below the fracture line, we acquired standard and new μCT-derived measurements. The mineralized callus volume and the mineral density of the callus correlated positively with callus strength (r_xy=?0.315, p=0.016 and r_xy=0.444, p<0.0005, respectively) and stiffness (r_xy=?0.271, p=0.040 and r_xy=0.325, p=0.013, respectively), but the fraction of the callus that mineralized and the moment of inertia of the callus did not. This fraction did correlate with energy-to-failure (r_xy=?0.343, p=0.0085). Of the μCT-derived measurements, quantifying defects within the outer bridging cortices of the callus produced the strongest correlation with both callus strength (r_xy=0.557, p<0.0001) and stiffness (r_xy=0.468, p=0.0002). By both reducing structural defects and increasing mineralization, lovastatin appears to increase the callus strength.     

Image interpolation allows accurate quantitative bone morphometry in registered micro-computed tomography scans

Time-lapsed in vivo micro-computed tomography is a powerful tool to analyse longitudinal changes in the bone micro-architecture. Registration can overcome problems associated with spatial misalignment between scans; however, it requires image interpolation which might affect the outcome of a subsequent bone morphometric analysis. The impact of the interpolation error itself, though, has not been quantified to date. Therefore, the purpose of this ex vivo study was to elaborate the effect of different interpolator schemes [nearest neighbour, tri-linear and B-spline (BSP)] on bone morphometric indices. None of the interpolator schemes led to significant differences between interpolated and non-interpolated images, with the lowest interpolation error found for BSPs (1.4%). Furthermore, depending on the interpolator, the processing order of registration, Gaussian filtration and binarisation played a role. Independent from the interpolator, the present findings suggest that the evaluation of bone morphometry should be done with images registered using greyscale information.     

A micro-computed tomography study of the trabecular bone structure in the femoral head

The goal of this study was to characterize the trabecular microarchitecture of the femoral head using micro-computed tomography (ICT). Femoral head specimens were obtained from subjects following total hip replacement. Cylindrical cores from the specimens were scanned to obtain 3-D images with an isotropic resolution of 26 Im. Bone structural parameters were evaluated on a per millimeter basis: relative bone volume (BV/TV), trabecular number (Tb.N), thickness (Tb.Th) and separation (Tb.Sp), structure model index (SMI), and connectivity (Conn.D). The ICT data show that the first two millimeters, starting at the joint surface, are characterized by more plate-like trabeculae, and are significantly denser than the underlying trabecular bone. Regional differences in the trabecular architecture reveal that the superior pole has significantly higher BV/TV, Tb.N and Tb.Th values, with lower Tb.Sp compared to the inferior and side poles. Because subchondral bone is essential in the load attenuation of joints, the difference in bone structure between the subchondral and trabecular bone might arise from the different functions each have within joint-forming bones. The denser trabecular structure of the superior pole as compared to the inferior pole can be interpreted as a functional adaptation to higher loading in this area.     

X-ray tomography of whole cells

X-ray tomography has been shown to provide insights into the internal structure of whole cells that can't be obtained by any other means. With recent advances in instrumentation and the advent of automated cryogenic sample stages, it has become possible to collect isotropic tomographic data from radiation-sensitive cells and to compute reconstructions with a high degree of fidelity to a resolution of 50 nm. The new generation of X-ray optics will extend this resolution limit to 15 nm or better. The development of X-ray tomography of whole cells generates opportunities to study cells, and cellular processes, in a completely new way.     

Mapping molecular agents distributions in whole mice hearts using born-normalized optical projection tomography

To date there is a lack of tools to map the spatio-temporal dynamics of diverse cells in experimental heart models. Conventional histology is labor intensive with limited coverage, whereas many imaging techniques do not have sufficiently high enough spatial resolution to map cell distributions. We have designed and built a high resolution, dual channel Born-normalized near-infrared fluorescence optical projection tomography system to quantitatively and spatially resolve molecular agents distribution within whole murine heart. We validated the use of the system in a mouse model of monocytes/macrophages recruitment during myocardial infarction. While acquired, data were processed and reconstructed in real time. Tomographic analysis and visualization of the key inflammatory components were obtained via a mathematical formalism based on left ventricular modeling. We observed extensive monocyte recruitment within and around the infarcted areas and discovered that monocytes were also extensively recruited into non-ischemic myocardium, beyond that of injured tissue, such as the septum.     

Development of a fluorescent light technique for evaluating microdamage in bone subjected to fatigue loading.

A new method using fluorescent light microscopy has been developed to visualize and evaluate bone microdamage. We report the findings of two different experiments with a common aim of comparing the fluorescent light technique to the brightfield method for quantifying microdamage in bone. In Experiment 1, 36 canine femurs were tested in four-point cyclic bending until they had lost between 5 and 43% of their stiffness. The loaded portion of the bone was stained en bloc with basic fuchsin for the presence of damage. Standard point counting techniques were used to calculate fractional damaged area (Dm.Ar = Cr.Ar/B.Ar, mm2/mm2) under brightfield and fluorescent microscopy. In Experiment 2, bone microdamage adjacent to endosseous implants, subjected to fatigue loading (150,000 cycles, 2 Hz and 37 degrees C) ex vivo was examined. The bone around the implant was either allowed to heal (adapted specimen) for 12 weeks after placement in dog mid-femoral diaphyses prior to testing or was loaded immediately to simulate non-healed bone surrounding endosseous implants (non-adapted). Crack numerical density (Cr.Dn = Cr.N/B.Ar, #/mm2), crack surface density (Cr.S.Dn = Tt.Cr.Le/B.Ar, mm/mm2) and fractional damaged area were calculated separately by both techniques in the adapted and non-adapted specimens. In both Experiments 1 and 2, significantly more microdamage was detected by the fluorescent technique than by the brightfield method. Also, there was a trend towards higher intraobserver repeatability when using the fluorescent method. These results suggest that the brightfield technique underestimates microdamage accumulation and that the fluorescent technique better represents the actual amounts of microdamage present. The results demonstrate that the fluorescent method provides an accurate and precise approach for bone microdamage evaluation, and that it improves the prediction of stiffness loss from damage accumulation.     

Synchrotron X-ray micro-computed tomography as a tool for in situ elucidation of insect bacteriomes

Obligate bacterial endosymbionts are common, influential associates of arthropods, and are often found in specific organs termed bacteriomes. Three dimensional images of bacteriomes of the leafhopper Orosius albicinctus (Hemiptera: Cicadellidae) were reconstructed from synchrotron-based X-ray micro-computed tomography (CT). Results show that bilateral bacteriomes are located between the first and second abdominal tergites, are mushroom-shaped and consist two different types of tissue. Fluorescence in situ hybridization reveals that the primary bacterial symbiont Sulcia muelleri is in the ‘cap’ part of the of organ. The technique allows a noninvasive, in situ, means of visualizing bacteriomes and will facilitate understanding their form and function.     

A comparison of the femoral head and neck trabecular architecture of Galago and Perodicticus using micro-computed tomography (microCT)

Innovations in micro-computed tomography (microCT) in the medical field have resulted in the development of techniques that allow the precise quantification of bone density and fabric related parameters of trabecular bone. For the purpose of this study, the technique was applied to a small sample of Perodicticus potto and Galago senegalensis femora to see if differences in loading environment elicit the predicted effects on trabecular structure. While the overall bone volume was approximately three times larger in the potto, there was no significant difference in the apparent volume density in the two taxa. When regional differences in the proximal femur were examined, the cancellous bone of the femoral head of Perodicticus potto and Galago senegalensis, while not differing in volume density, showed differences in trabecular orientation, with the potto having more randomly oriented trabeculae than the bushbaby. This was as hypothesized, given that the bushbaby submits its femora to more stereotypical loading environments than the potto. In the femoral neck, the cancellous bone was not only more randomly oriented, it was also denser in the potto compared with the bushbaby. This suggests that trabecular morphology may be extremely sensitive to certain differences in the loading environment and that this information, combined with information on cortical bone structure and external geometry, will result in a more complete understanding of how bone shape and composition correspond to loading and locomotor patterns. Ultimately, a synthesis of these different lines of evidence may have considerable applications in paleontological studies that attempt to reconstruct bone use from morphology.     

Predicting the permeability of trabecular bone by micro-computed tomography and finite element modeling

The intrinsic permeability of bone plays an important role in the transport of nutrients and minerals within the tissue, and affects the mechanical stimuli that are related to the fate of the stem cells. The objective of this study was to establish a method to assess trabecular bone permeability using experimental and finite element (FE) modeling approaches based on micro computed tomography (µCT) images. Human cadaveric tibia cube specimens (N=23) were scanned with µCT. The permeability was measured experimentally using a custom-developed constant-head permeameter, and computationally by a poroelastic formulation to simulate the fluid flow within the discretized bone matrix and pore phase. The average of the experimentally measured permeability was 4.84×10⁻¹⁰ m² with a standard deviation of 3.70×10⁻¹⁰ m². A regression model of the µCT determined that the maximum bone area to total area ratio (maxBA/TA) for all slices that are perpendicular to the direction of fluid flow explained 84% of the variability of the natural logarithm of the experimentally measured permeability. The 2D measure of maxBA/TA performed better than 3D measures in general, although some parameters were reasonably well associated with permeability such as bone volume ratio (BV/TV, r=−0.71), the bone surface/bone volume (BS/BV, r=0.73), and the trabecular thickness (TbTh, r=−0.71). The correlation between the permeability predicted with FE models and experimentally measured permeability was reasonable (r=0.69), but the FE approach did not accurately represent the wide variability of permeability measured experimentally. The results of this study suggest that the changes in the trabecular bone microarchitecture have an exponential relationship with permeability, and the use of µCT-based 2D measurement of maxBA/TA performs well at predicting permeability, thus providing a convenient approach to measure this important aspect affecting biomechanical functions in the tissue.