Clinical Super-Resolution Computed Tomography of Bone Microstructure: Application in Musculoskeletal and Dental Imaging.

PURPOSE: Clinical cone-beam computed tomography (CBCT) devices are limited to imaging features of half a millimeter in size and cannot quantify the tissue microstructure. We demonstrate a robust deep-learning method for enhancing clinical CT images, only requiring a limited set of easy-to-acquire training data. METHODS: Knee tissue from five cadavers and six total knee replacement patients, and 14 teeth from eight patients were scanned using laboratory CT as training data for the developed super-resolution (SR) technique. The method was benchmarked against ex vivo test set, 52 osteochondral samples are imaged with clinical and laboratory CT. A quality assurance phantom was imaged with clinical CT to quantify the technical image quality. To visually assess the clinical image quality, musculoskeletal and maxillofacial CBCT studies were enhanced with SR and contrasted to interpolated images. A dental radiologist and surgeon reviewed the maxillofacial images. RESULTS: The SR models predicted the bone morphological parameters on the ex vivo test set more accurately than conventional image processing. The phantom analysis confirmed higher spatial resolution on the SR images than interpolation, but image grayscales were modified. Musculoskeletal and maxillofacial CBCT images showed more details on SR than interpolation; however, artifacts were observed near the crown of the teeth. The readers assessed mediocre overall scores for both SR and interpolation. The source code and pretrained networks are publicly available. CONCLUSION: Model training with laboratory modalities could push the resolution limit beyond state-of-the-art clinical musculoskeletal and dental CBCT. A larger maxillofacial training dataset is recommended for dental applications.

Cerebral tissue oxygenation response to brain irradiation measured during clinical radiotherapy.

Significance: Cancer therapy treatments produce extensive changes in the physiological and morphological properties of tissues, which are also individual dependent. Currently, a key challenge involves developing more tailored cancer therapy, and consequently, individual biological response measurement during therapy, such as tumor hypoxia, is of high interest. This is the first time human cerebral haemodynamics and cerebral tissue oxygenation index (TOI) changes were measured during the irradiation in clinical radiotherapy and functional near-infrared spectroscopy (fNIRS) technique was demonstrated as a feasible technique for clinical use in radiotherapy, based on 34 online patient measurements. Aim: Our aim is to develop predictive biomarkers and noninvasive real-time methods to establish the effect of radiotherapy during treatment as well as to optimize radiotherapy dose planning for individual patients. In particular, fNIRS-based technique could offer an effective and clinically feasible online technique for continuous monitoring of brain tissue hypoxia and responses to chemo- and radiotherapy, which involves modulating tumor oxygenation to increase or decrease tumor hypoxia. We aim to show that fNIRS is feasible for repeatability measuring in patient radiotherapy, the temporal alterations of tissue oxygenation induced by radiation. Approach: Fiber optics setup using multiwavelength fNIRS was built and combined with a medical linear accelerator to measure cerebral tissue oxygenation changes during the whole-brain radiotherapy treatment, where the radiation dose is given in whole brain area only preventing dosage to eyes. Correlation of temporal alterations in cerebral haemodynamics and TOI response to brain irradiation was quantified. Results: Online fNIRS patient measurement of cerebral haemodynamics during clinical brain radiotherapy is feasible in clinical environment, and results based on 34 patient measurements show strong temporal alterations in cerebral haemodynamics and decrease in TOI during brain irradiation and confirmed the repeatability. Our proof-of-concept study shows evidently that irradiation causes characteristic immediate changes in brain tissue oxygenation. Conclusions: In particular, TOI seems to be a sensitive parameter to observe the tissue effects of radiotherapy. Monitoring the real-time interactions between the subjected radiation dose and corresponding haemodynamic effects may provide important tool for the researchers and clinicians in the field of radiotherapy. Eventually, presented fNIRS technique could be used for improving dose planning and safety control for individual patients.

Bright ultrashort echo time SWIFT MRI signal at the osteochondral junction is not located in the calcified cartilage.

In this study, we aimed to precisely localize the hyperintense signal that is generated at the osteochondral junction when using ultrashort echo time magnetic resonance imaging (MRI) and to investigate the osteochondral junction using sweep imaging with Fourier transformation (SWIFT) MRI. Furthermore, we seek to evaluate what compositional properties of the osteochondral junction are the sources of this signal. In the study, we obtained eight samples from a tibial plateau dissected from a 68-year-old male donor, and one additional osteochondral sample of bovine origin. The samples were imaged using high-resolution ultrashort echo time SWIFT MRI and microcomputed tomography (µCT) scans. Localization of the bright signal in the osteochondral junction was performed using coregistered data sets. Potential sources of the signal feature were examined by imaging the bovine specimen with variable receiver bandwidths and by performing variable flip angle T1 relaxation time mapping. The results of the study showed that the hyperintense signal was found to be located entirely in the deep noncalcified articular cartilage. The intensity of this signal at the interface varied between the specimens. Further tests with bovine specimens indicated that the imaging bandwidth and T1 relaxation affect the properties of the signal. Based on the present results, the calcified cartilage has low signal intensity even in SWIFT imaging. Concomitantly, it appears that the bright signal seen in ultrashort echo time imaging resides within the noncalcified cartilage. Furthermore, the most likely sources of this signal are the rapid T1 relaxation of the deep cartilage and the susceptibility-induced effects arising from the calcified tissues.

Raman microspectroscopic analysis of the tissue-specific composition of the human osteochondral junction in osteoarthritis: A pilot study.

This study investigates the influence of osteoarthritis (OA) disease severity on the bio-composition of the osteochondral junction at the human tibial plateau using Raman microspectroscopy. We specifically aim to analyze the spatial composition of mineralized osteochondral tissues, i.e., calcified cartilage (CC) and subchondral bone plate (SBP) from unfixed, hydrated specimens. We hypothesize that the mineralization of CC and SBP decreases in advanced OA. Twenty-eight cylindrical osteochondral samples (d = 4 mm) from tibial plateaus of seven cadaveric donors were harvested and sorted into three groups following histopathological grading: healthy (n = 5), early OA (n = 8), and advanced OA (n = 15). Raman spectra were subjected to multivariate cluster analyses to identify different tissues. Finally, the tissue-specific composition was analyzed, and the impact of OA was statistically evaluated with linear mixed models. Cluster analyses of Raman spectra successfully distinguished CC and SBP as well as a tidemark region and uncalcified cartilage. CC was found to be more mineralized and the mineral was more crystalline compared with SBP. Both tissues exhibited similar compositional changes as a function of histopathological OA severity. In early OA, the mineralization tends to increase, and the mineral contains fewer carbonate substitutions. Compared with early OA, mineral crystals are rich in carbonate while the overall mineralization decreases in advanced OA. This Raman spectroscopic study advances the methodology for investigating the complex osteochondral junction from native tissue. The developed methodology can be used to elucidate detailed tissue-specific changes in the chemical composition with advancing OA. STATEMENT OF SIGNIFICANCE: In this study, Raman microspectroscopy was utilized to investigate the influence of osteoarthritic degeneration on the tissue-specific biochemical composition of the human osteochondral junction. Multivariate cluster analyses allowed us to characterize subtle compositional changes in the calcified cartilage and subchondral bone plate as well as in the tidemark region. The compositional differences found between the calcified cartilage and subchondral bone plate in both organic and mineral phases will serve as critical benchmark parameters when designing biomaterials for osteochondral repair. We found tissue-specific changes in the mineralization and carbonate substitution as a function of histopathological OA severity. Our developed methodology can be used to investigate the metabolic changes in the osteochondral junction associated with osteoarthritis.

Quantifying Complex Micro-Topography of Degenerated Articular Cartilage Surface by Contrast-Enhanced Micro-Computed Tomography and Parametric Analyses.

One of the earliest changes in osteoarthritis (OA) is a surface discontinuity of the articular cartilage (AC), and these surface changes become gradually more complex with OA progression. We recently developed a contrast enhanced micro-computed tomography (µCT) method for visualizing AC surface in detail. The present study aims to introduce a µCT analysis technique to parameterize these complex AC surface features and to demonstrate the feasibility of using these parameters to quantify degenerated AC surface. Osteochondral plugs (n = 35) extracted from 19 patients undergoing joint surgery were stained with phosphotungstic acid and imaged using µCT. The surface micro-topography of AC was analyzed with developed method. Standard root mean square roughness (Rq ) was calculated as a reference, and the Area Under Curve (AUC) for receiver operating characteristic analysis was used to compare the acquired quantitative parameters with semi-quantitative visual grading of µCT image stacks. The parameters quantifying the complex micro-topography of AC surface exhibited good sensitivity and specificity in identifying surface continuity (AUC: 0.93, [0.80 0.99]), fissures (AUC: 0.94, [0.83 0.99]) and fibrillation (AUC: 0.98, [0.88 1.0]). Standard Rq was significantly smaller compared with the complex roughness (CRq ) already with mild surface changes with all surface reference parameters - continuity, fibrillation, and fissure sum. Furthermore, only CRq showed a significant difference when comparing the intact surface with lowest fissure sum score. These results indicate that the presented method for evaluating complex AC surfaces exhibit potential to identify early OA changes in superficial AC and is dynamic throughout OA progression. © 2019 The Authors. Journal of Orthopaedic Research® Published by Wiley Periodicals, Inc. on behalf of the Orthopaedic Research Society. Society. 9999:1-12, 2019.

An Automatic Regularization Method: An Application for 3-D X-Ray Micro-CT Reconstruction Using Sparse Data.

X-ray tomography is a reliable tool for determining the inner structure of 3-D object with penetrating X-rays. However, traditional reconstruction methods, such as Feldkamp-Davis-Kress (FDK), require dense angular sampling in the data acquisition phase leading to long measurement times, especially in X-ray micro-tomography to obtain high-resolution scans. Acquiring less data using greater angular steps is an obvious way for speeding up the process and avoiding the need to save huge data sets. However, computing 3-D reconstruction from such a sparsely sampled data set is difficult because the measurement data are usually contaminated by errors, and linear measurement models do not contain sufficient information to solve the problem in practice. An automatic regularization method is proposed for robust reconstruction, based on enforcing sparsity in the 3-D shearlet transform domain. The inputs of the algorithm are the projection data and a priori known expected degree of sparsity, denoted as . The number Cpr can be calibrated from a few dense-angle reconstructions and fixed. Human subchondral bone samples were tested, and morphometric parameters of the bone reconstructions were then analyzed using standard metrics. The proposed method is shown to outperform the baseline algorithm (FDK) in the case of sparsely collected data. The number of X-ray projections can be reduced up to 10% of the total amount 300 projections over 180° with uniform angular step while retaining the quality of the reconstruction images and of the morphometric parameters.

Iterative and discrete reconstruction in the evaluation of the rabbit model of osteoarthritis.

Micro-computed tomography (µCT) is a standard method for bone morphometric evaluation. However, the scan time can be long and the radiation dose during the scan may have adverse effects on test subjects, therefore both of them should be minimized. This could be achieved by applying iterative reconstruction (IR) on sparse projection data, as IR is capable of producing reconstructions of sufficient image quality with less projection data than the traditional algorithm requires. In this work, the performance of three IR algorithms was assessed for quantitative bone imaging from low-resolution data in the evaluation of the rabbit model of osteoarthritis. Subchondral bone images were reconstructed with a conjugate gradient least squares algorithm, a total variation regularization scheme, and a discrete algebraic reconstruction technique to obtain quantitative bone morphometry, and the results obtained in this manner were compared with those obtained from the reference reconstruction. Our approaches were sufficient to identify changes in bone structure in early osteoarthritis, and these changes were preserved even when minimal data were provided for the reconstruction. Thus, our results suggest that IR algorithms give reliable performance with sparse projection data, thereby recommending them for use in µCT studies where time and radiation exposure are preferably minimized.

Correlations of low-field NMR and variable-field NMR parameters with osteoarthritis in human articular cartilage under load.

NMR experiments carried out at magnetic fields below 1 T provide new relaxation parameters unavailable with conventional clinical scanners. Contrast of T1 generally becomes larger towards low fields, as slow molecular reorientation processes dominate relaxation at the corresponding Larmor frequencies. This advantage has to be considered in the context of lower sensitivity and frequently reduced spatial resolution. The layered structure of cartilage is one example where a particularly strong variation of T1 across the tissue occurs, being affected by degenerative diseases such as osteoarthritis (OA). Furthermore, the presence of 1 H-14  N cross-relaxation, leading to so-called quadrupolar dips in the 1 H relaxation time dispersion, provide insight into the concentration and mobility of proteoglycans and collagen in cartilage, both being affected by OA. In this study, low-field imaging and variable-field NMR relaxometry were combined for the first time for tissue samples, employing unidirectional load to probe the mechanical properties. 20 human knee cartilage samples were placed in a compression cell, and studied by determining relaxation profiles without and with applied pressure (0.6 MPa) at 50 µm in-plane resolution, and comparing with volume-averaged T1 dispersion. Samples were subsequently stored in formalin, prepared for histology and graded according to the Mankin score system. Quadrupolar dips and thickness change under load showed the strongest correlation with Mankin grade. Average T1 and change of maximum T1 under load, as well as its position, correlate with thickness and thickness change. Furthermore, T1 (ω) above 25 mT was found to correlate with thickness change. While volume-averaged T1 is not a suitable indicator for OA, its change due to mechanical load and its extreme values are suggested as biomarkers available in low-field MRI systems. The shape of the dispersion T1 (ω) represents a promising access to understanding and quantifying molecular dynamics in tissue, pointing toward future in vivo tissue studies.

Effects of Articular Cartilage Constituents on Phosphotungstic Acid Enhanced Micro-Computed Tomography.

Contrast-enhanced micro-computed tomography (CEµCT) with phosphotungstic acid (PTA) has shown potential for detecting collagen distribution of articular cartilage. However, the selectivity of the PTA staining to articular cartilage constituents remains to be elucidated. The aim of this study was to investigate the dependence of PTA for the collagen content in bovine articular cartilage. Adjacent bovine articular cartilage samples were treated with chondroitinase ABC and collagenase to degrade the proteoglycan and the collagen constituents in articular cartilage, respectively. Enzymatically degraded samples were compared to the untreated samples using CEµCT and reference methods, such as Fourier-transform infrared imaging. Decrease in the X-ray attenuation of PTA in articular cartilage and collagen content was observed in cartilage depth of 0-13% and deeper in tissue after collagen degradation. Increase in the X-ray attenuation of PTA was observed in the cartilage depth of 13-39% after proteoglycan degradation. The X-ray attenuation of PTA-labelled articular cartilage in CEµCT is associated mainly with collagen content but the proteoglycans have a minor effect on the X-ray attenuation of the PTA-labelled articular cartilage. In conclusion, the PTA labeling provides a feasible CEµCT method for 3D characterization of articular cartilage.

Micro-CT Analysis of Bone Healing in Rabbit Calvarial Critical-Sized Defects with Solid Bioactive Glass, Tricalcium Phosphate Granules or Autogenous Bone.

OBJECTIVES: The purpose of the present study was to evaluate bone healing in rabbit critical-sized calvarial defects using two different synthetic scaffold materials, solid biodegradable bioactive glass and tricalcium phosphate granules alongside solid and particulated autogenous bone grafts. MATERIAL AND METHODS: Bilateral full thickness critical-sized calvarial defects were created in 15 New Zealand white adult male rabbits. Ten defects were filled with solid scaffolds made of bioactive glass or with porous tricalcium phosphate granules. The healing of the biomaterial-filled defects was compared at the 6 week time point to the healing of autologous bone grafted defects filled with a solid cranial bone block in 5 defects and with particulated bone combined with fibrin glue in 10 defects. In 5 animals one defect was left unfilled as a negative control. Micro-computed tomography (micro-CT) was used to analyze healing of the defects. RESULTS: Micro-CT analysis revealed that defects filled with tricalcium phosphate granules showed new bone formation in the order of 3.89 (SD 1.17)% whereas defects treated with solid bioactive glass scaffolds showed 0.21 (SD 0.16)%, new bone formation. In the empty negative control defects there was an average new bone formation of 21.8 (SD 23.7)%. CONCLUSIONS: According to findings in this study, tricalcium phosphate granules have osteogenic potential superior to bioactive glass, though both particulated bone with fibrin glue and solid bone block were superior defect filling materials.

Bone healing in rabbit calvarial critical-sized defects filled with stem cells and growth factors combined with granular or solid scaffolds.

PURPOSE: In pediatric neurosurgery, decompressive craniectomy and correction of congenital cranial anomalies can result in major cranial defects. Corrective cranioplasty for the repair of these critical-sized defects is not only a cosmetic issue. The limited availability of suitable autogenous bone and the morbidity of donor site harvesting have driven the search for new approaches with biodegradable and bioactive materials. This study aimed to assess the healing of rabbit calvarial critical-sized defects filled with osteogenic material, either with bioactive glass scaffolds or tricalcium phosphate granules in various combinations with adipose stem cells or bone marrow stem cells, BMP-2, BMP-7, or VEGF to enhance osteogenesis. METHODS: Eighty-two bicortical full thickness critical-sized calvarial defects were operated. Five defects were left empty as negative control defects. The remaining 77 defects were filled with solid bioactive glass scaffolds or tricalcium phosphate granules seeded with adipose or bone marrow derived stem cells in combination with BMP-2, BMP-7, or VEGF. The defects were allowed to heal for 6 weeks before histologic and micro-CT analyses. RESULTS: Micro-CT examination at the 6-week post-operative time point revealed that defects filled with stem cell-seeded tricalcium phosphate granules resulted in new bone formation of 6.0 %, whereas defects with bioactive glass scaffolds with stem cells showed new bone formation of 0.5 to 1.7 %, depending on the growth factor used. CONCLUSIONS: This study suggests that tricalcium phosphate granules combined with stem cells have osteogenic potential superior to solid bioactive glass scaffolds with stem cells and growth factors.

Imaging of subchondral bone by optical coherence tomography upon optical clearing of articular cartilage.

Optical clearing is an effective method to reduce light scattering of biological tissues that provides significant enhancement of light penetration into the biological tissues making non-invasive diagnosis more feasible. In current report Optical Coherence Tomography (OCT) in conjunction with optical clearing is applied for assessment of deep cartilage layers and cartilage-bone interface. The solution of Iohexol in water has been used as an optical clearing agent. The cartilage-bone boundary becomes visible after 15 min of optical clearing that enabling non-invasive estimation of its roughness: Sa = 10 ± 1 µm. The results show that for 0.9 mm thick cartilage optical clearing is stopped after 50 min with an increase of refractive index from 1.386 ± 0.008 to 1.510 ± 0.009. Current approach enables more reliable detection of arthroscopically inaccessible regions, including cartilage-bone boundary and subchondral bone, and potentially improves accuracy of the osteoarthritis diagnosis.