Characterization complex collagen fiber architecture in knee joint using high-resolution diffusion imaging

Explore the Pattern of Biomechanical Alterations in Vocal Fold Scar and Its Objective Quantitative Assessment Method

OBJECTIVES

This research aims to discern the evolving nature of the biomechanical properties of vocal fold scarring by calculating Young's modulus for the vocal fold cover layer, the body layer, and the structure as a whole. The study also investigates the potential of diffusion tensor imaging (DTI) for determining these biomechanical characteristics quantitatively.

METHODS

A total of six adult female Beagles were divided into two groups (A and B groups) for the creation of unilateral vocal fold scar models, each group containing three subjects. Five months postmodel creation, larynxes were excised and placed within a 9.4T BioSpec MRI system (Bruker, Germany) for scanning. Subsequently, the vocal folds were segregated from the larynx. In A group of Beagles, the vocal fold cover layer and body layers were separated, whereas in B group they remained intact. All samples were then subjected to cyclic tensile testing using an Instron MicroTester 5948, with Young's modulus computed for the vocal fold cover layer and body layers in the A group and for the intact vocal fold in the B group. Differences in the overall Young's modulus between the vocal fold scarred side and the healthy side were analyzed, and a Pearson correlation analysis was performed between DTI parameters and the outcomes of the stress-strain experiments.

RESULTS

A statistically significant discrepancy in the overall Young's modulus was identified between the scar and healthy sides of the vocal fold (P = 0.0401). The Young's modulus also displayed a significant difference between the scar and healthy sides of the vocal fold cover layer (P = 0.0241). No meaningful divergence was observed in the elastic modulus between the scar and healthy sides of the vocal fold body layer (P > 0.05). Postseparation, Young's modulus for both the cover and body layers of the scarred vocal fold were less than that of the same layers on the healthy side. However, Young's modulus of the entirety of the vocal fold on the scar side was greater than that of the whole vocal fold on the healthy side. The fractional anisotropy (FA) of the vocal fold cover layer had a significant correlation with the elastic modulus (r = 0.812, P = 0.050), as did the Tensor trace (r = -0.821, P = 0.045). The FA of the vocal fold body layer showed no significant correlation with the elastic modulus (r = -0.725, P = 0.103), while the Tensor trace demonstrated a significant correlation (r = 0.911, P = 0.012).

CONCLUSIONS

Biomechanical alterations in vocal fold scars demonstrate a closer association with adhesion bands, thus emphasizing the importance of adhesion band loosening for the restoration of vibratory function within vocal fold scarring. DTI emerges as a potent noninvasive quantitative instrument for assessing these biomechanical changes, as well as for quantitatively gauging the severity of vocal fold scarring.

Simultaneous multi-slice technique for reducing acquisition times in diffusion tensor imaging of the knee: a feasibility study

OBJECTIVES

To explore the feasibility of simultaneous multi-slice (SMS) technique for reducing acquisition times in readout-segmented echo planar imaging (RESOLVE) for diffusion tensor imaging (DTI) of the knee.

MATERIALS AND METHODS

A total of 30 healthy volunteers and 23 patients with knee acute injury (12 cases with anterior ligament (ACL) tears and 16 cases with patellar cartilage (PC) injury) were enrolled in this prospective study. Three DTI protocols were used: conventional RESOLVE-DTI with 12 directions (protocol 1), SMS-RESOLVE-DTI with 12 directions (protocol 2) and 20 directions (protocol 3). DTI parameters of gastrocnemius, ACL and posterior cruciate ligament (PCL), and PC from three protocols were quantitatively assessed.

RESULTS

For volunteers, protocol 2 significantly reduced acquisition time by 38.6% and 34.2% compared to protocols 1 and 3 while maintaining similar high-quality images and similar diffusive parameters, except for the fractional anisotropy (FA) and axial diffusivity (AD) of the PC between protocols 2 and 1 (P < 0.05). For injured ACL and PC, protocols 1 and 2 showed similar accurate diffusive parameters (except for AD, P = 0.025) and similar diagnostic efficacy, which demonstrated significantly lower FA and higher radial diffusivity (RD) in protocols 1 and 2 compared to volunteers (P < 0.05).

CONCLUSIONS

The 12-direction SMS-RESOLVE-DTI demonstrated a favorable balance between acquisition time and image quality, making it a promising alternative to conventional DTI for evaluating ligament and cartilage injuries.

ADVANCES IN KNOWLEDGE

The SMS technique greatly reduces acquisition time while maintaining image quality, which signified the possibility of DTI's clinical application.

Quantitative MRI methods for the assessment of structure, composition, and function of musculoskeletal tissues in basic research and preclinical applications

Osteoarthritis (OA) is a disabling chronic disease involving the gradual degradation of joint structures causing pain and dysfunction. Magnetic resonance imaging (MRI) has been widely used as a non-invasive tool for assessing OA-related changes. While anatomical MRI is limited to the morphological assessment of the joint structures, quantitative MRI (qMRI) allows for the measurement of biophysical properties of the tissues at the molecular level. Quantitative MRI techniques have been employed to characterize tissues' structural integrity, biochemical content, and mechanical properties. Their applications extend to studying degenerative alterations, early OA detection, and evaluating therapeutic intervention. This article is a review of qMRI techniques for musculoskeletal tissue evaluation, with a particular emphasis on articular cartilage. The goal is to describe the underlying mechanism and primary limitations of the qMRI parameters, their association with the tissue physiological properties and their potential in detecting tissue degeneration leading to the development of OA with a primary focus on basic and preclinical research studies. Additionally, the review highlights some clinical applications of qMRI, discussing the role of texture-based radiomics and machine learning in advancing OA research.

High angular resolution diffusion imaging (HARDI) of porcine menisci: a comparison of diffusion tensor imaging and generalized q-sampling imaging

Background

Diffusion magnetic resonance imaging (MRI) allows for the quantification of water diffusion properties in soft tissues. The goal of this study was to characterize the 3D collagen fiber network in the porcine meniscus using high angular resolution diffusion imaging (HARDI) acquisition with both diffusion tensor imaging (DTI) and generalized q-sampling imaging (GQI).

Methods

Porcine menisci (n=7) were scanned ex vivo using a three-dimensional (3D) HARDI spin-echo pulse sequence with an isotropic resolution of 500 µm at 7.0 Tesla. Both DTI and GQI reconstruction techniques were used to quantify the collagen fiber alignment and visualize the complex collagen network of the meniscus. The MRI findings were validated with conventional histology.

Results

DTI and GQI exhibited distinct fiber orientation maps in the meniscus using the same HARDI acquisition. We found that crossing fibers were only resolved with GQI, demonstrating the advantage of GQI over DTI to visualize the complex collagen fiber orientation in the meniscus. Furthermore, the MRI findings were consistent with conventional histology.

Conclusions

HARDI acquisition with GQI reconstruction more accurately resolves the complex 3D collagen architecture of the meniscus compared to DTI reconstruction. In the future, these technologies have the potential to nondestructively assess both normal and abnormal meniscal structure.

Impact of Deep Learning Denoising Algorithm on Diffusion Tensor Imaging of the Growth Plate on Different Spatial Resolutions

To assess the impact of a deep learning (DL) denoising reconstruction algorithm applied to identical patient scans acquired with two different voxel dimensions, representing distinct spatial resolutions, this IRB-approved prospective study was conducted at a tertiary pediatric center in compliance with the Health Insurance Portability and Accountability Act. A General Electric Signa Premier unit (GE Medical Systems, Milwaukee, WI) was employed to acquire two DTI (diffusion tensor imaging) sequences of the left knee on each child at 3T: an in-plane 2.0 × 2.0 mm2 with section thickness of 3.0 mm and a 2 mm3 isovolumetric voxel; neither had an intersection gap. For image acquisition, a multi-band DTI with a fat-suppressed single-shot spin-echo echo-planar sequence (20 non-collinear directions; b-values of 0 and 600 s/mm2) was utilized. The MR vendor-provided a commercially available DL model which was applied with 75% noise reduction settings to the same subject DTI sequences at different spatial resolutions. We compared DTI tract metrics from both DL-reconstructed scans and non-denoised scans for the femur and tibia at each spatial resolution. Differences were evaluated using Wilcoxon-signed ranked test and Bland-Altman plots. When comparing DL versus non-denoised diffusion metrics in femur and tibia using the 2 mm × 2 mm × 3 mm voxel dimension, there were no significant differences between tract count (p = 0.1, p = 0.14) tract volume (p = 0.1, p = 0.29) or tibial tract length (p = 0.16); femur tract length exhibited a significant difference (p < 0.01). All diffusion metrics (tract count, volume, length, and fractional anisotropy (FA)) derived from the DL-reconstructed scans, were significantly different from the non-denoised scan DTI metrics in both the femur and tibial physes using the 2 mm3 voxel size (p < 0.001). DL reconstruction resulted in a significant decrease in femorotibial FA for both voxel dimensions (p < 0.01). Leveraging denoising algorithms could address the drawbacks of lower signal-to-noise ratios (SNRs) associated with smaller voxel volumes and capitalize on their better spatial resolutions, allowing for more accurate quantification of diffusion metrics.

Applications of Diffusion-Weighted MRI to the Musculoskeletal System

Diffusion-weighted imaging (DWI) is an established MRI technique that can investigate tissue microstructure at the scale of a few micrometers. Musculoskeletal tissues typically have a highly ordered structure to fulfill their functions and therefore represent an optimal application of DWI. Even more since disruption of tissue organization affects its biomechanical properties and may indicate irreversible damage. The application of DWI to the musculoskeletal system faces application-specific challenges on data acquisition including susceptibility effects, the low T2 relaxation time of most musculoskeletal tissues (2-70 msec) and the need for sub-millimetric resolution. Thus, musculoskeletal applications have been an area of development of new DWI methods. In this review, we provide an overview of the technical aspects of DWI acquisition including diffusion-weighting, MRI pulse sequences and different diffusion regimes to study tissue microstructure. For each tissue type (growth plate, articular cartilage, muscle, bone marrow, intervertebral discs, ligaments, tendons, menisci, and synovium), the rationale for the use of DWI and clinical studies in support of its use as a biomarker are presented. The review describes studies showing that DTI of the growth plate has predictive value for child growth and that DTI of articular cartilage has potential to predict the radiographic progression of joint damage in early stages of osteoarthritis. DTI has been used extensively in skeletal muscle where it has shown potential to detect microstructural and functional changes in a wide range of muscle pathologies. DWI of bone marrow showed to be a valuable tool for the diagnosis of benign and malignant acute vertebral fractures and bone metastases. DTI and diffusion kurtosis have been investigated as markers of early intervertebral disc degeneration and lower back pain. Finally, promising new applications of DTI to anterior cruciate ligament grafts and synovium are presented. The review ends with an overview of the use of DWI in clinical routine. EVIDENCE LEVEL: 5 TECHNICAL EFFICACY: Stage 3.

Ultrashort echo time magnetic resonance imaging of the osteochondral junction

Osteoarthritis is a common chronic degenerative disease that causes pain and disability with increasing incidence worldwide. The osteochondral junction is a dynamic region of the joint that is associated with the early development and progression of osteoarthritis. Despite the substantial advances achieved in the imaging of cartilage and application to osteoarthritis in recent years, the osteochondral junction has received limited attention. This is primarily related to technical limitations encountered with conventional MR sequences that are relatively insensitive to short T2 tissues and the rapid signal decay that characterizes these tissues. MR sequences with ultrashort echo time (UTE) are of great interest because they can provide images of high resolution and contrast in this region. Here, we briefly review the anatomy and function of cartilage, focusing on the osteochondral junction. We also review basic concepts and recent applications of UTE MR sequences focusing on the osteochondral junction.

Using Diffusion Tensor Imaging to Explore the Changes in the Microstructure of Canine Vocal Fold Scar Tissue

OBJECTIVE

To apply diffusion tensor imaging (DTI) in measurement of the diffusion characteristics of water molecules in vocal fold scar tissue, combined with the analysis of textural characteristics of collagen fibers in the cover layer of the vocal folds to explore the feasibility of DTI in the qualitative and quantitative diagnosis of vocal fold scars and the evaluation of microstructural changes of vocal fold scar tissue.

METHODS

A unilateral injury was created using micro-cup forceps in the left vocal fold of six beagles. The contralateral normal vocal fold was used as a self-control. Five months postinjury, the larynges were excised and placed into a magnetic resonance imaging (MRI) system (9.4T BioSpec MRI, Bruker, German) for scanning and extraction of the diffusion parameters, fractional anisotropy (FA) and tensor trace in the anterior, middle, and posterior portions of the vocal fold cover layer. These parameters were then analyzed for statistical significance between the scarred vocal fold and the normal vocal fold. After MRI scanning, the tissue of the vocal folds was divided into anterior, middle, and posterior parts for sectioning and staining with hematoxylin and eosin, and samples were subsequently digitally scanned for texture analysis. The irregularity parameters, energy, contrast, correlation, and homogeneity, of collagen fibers of the vocal folds and the mean gray value of collagen fibers were calculated by the gray-level co-occurrence matrix (GLCM) texture analysis method. The differences in the mean value of the two sides of the vocal fold were compared. In addition, Pearson correlation analysis was performed between DTI parameters and irregularity parameters.

RESULTS

The FA of the left vocal fold cover layer was significantly lower compared to the self-control group (P = 0.0366), and the tensor trace value on the left vocal fold cover layer was significantly higher compared to the self-control group (P = 0.0353). The FA was significantly higher in the anterior part of the right vocal fold cover layer compared to the middle and posterior parts of the same side (P = 0.0352), and the tensor trace was significantly lower in the anterior part of the right vocal fold cover layer compared to the middle and posterior parts of the same side (P = 0.0298). There were no significant differences in FA and tensor trace between the middle and posterior parts of the vocal fold cover layer. The mean gray value of the left vocal folds cover layer was significantly smaller than the right vocal fold cover layer (P = 0.0219), the energy of the left vocal fold cover layer was significantly smaller than that of the right vocal fold cover layer (P < 0.0001), the contrast of the left vocal folds cover layer was significantly larger than that of the right vocal fold cover layer (P = 0.0002), the correlation of the left vocal folds cover layer was significantly smaller than the right vocal fold cover layer (P = 0.0002), and the homogeneity of the left vocal folds cover layer was significantly smaller than the right vocal fold cover layer (P = 0.0003). Pearson correlation analysis yielded values of r = 0.926, P = 0.000 between the FA and mean gray value; r = -0.918, P = 0.000 between FA and energy; r = -0.924, P = 0.000 between the FA and homogeneity, r = -0.949, P = 0.000 between tensor trace and mean gray value; r = 0.893, P = 0.000 between the tensor trace and energy; and r = 0.929, P = 0.000 between the tensor trace and homogeneity.

CONCLUSION

FA and tensor trace can be used as effective parameters to reflect microstructural changes in vocal fold scars. DTI is an objective and quantitative method of analyzing vocal fold scarring, and it noninvasively evaluates the microstructure of vocal fold collagen fibers.

Magic angle effect on diffusion tensor imaging in ligament and brain

PURPOSE

To evaluate the magic angle effect on diffusion tensor imaging (DTI) measurements in rat ligaments and mouse brains.

METHODS

Three rat knee joints and three mouse brains were scanned at 9.4 T using a modified 3D diffusion-weighted spin echo pulse sequence with the isotropic spatial resolution of 45 μm. The b value was 1000 s/mm2 for rat knee and 4000 s/mm2 for mouse brain. DTI model was used to investigate the quantitative metrics at different orientations with respect to the main magnetic field. The collagen fiber structure of the ligament was validated with polarized light microscopy (PLM) imaging.

RESULTS

The signal intensity, signal-to-noise ratio (SNR), and DTI metrics in the ligament were strongly dependent on the collagen fiber orientation with respect to the main magnetic field from both simulation and actual MRI scans. The variation of fractional anisotropy (FA) was about ~32%, and the variation of mean diffusivity (MD) was ~11%. These findings were further validated with the numerical simulation at different SNRs (~10.0 to 86.0). Compared to the ligament, the DTI metrics showed little orientation dependence in mouse brains.

CONCLUSION

Magic angle effect plays an important role in DTI measurements in the highly ordered collagen-rich tissues, while MD showed less orientation dependence than FA.

Single shot quantitative polarized light imaging system for rapid planar biaxial testing of soft tissues

Quantitative Polarized Light Imaging (QPLI) is an established technique used to compute the orientation of collagen fibers based on their birefringence. QPLI systems typically require rotating linear polarizers to obtain sufficient data to estimate orientation, which limits acquisition speeds and is not ideal for its application to mechanical testing. In this paper, we present a QPLI system designed with no moving parts; a single shot technique which is ideal to characterize collagen fiber orientation and kinematics during mechanical testing. Our single shot QPLI system (ssQPLI) sorts polarized light into four linear polarization states that are collected simultaneously by four cameras. The ssQPLI system was validated using samples with known orientation and retardation, and we demonstrate its use with planar biaxial testing of mouse skin. The ssQPLI system was accurate with a mean orientation error of 1.35° ± 1.58°. Skin samples were tested with multiple loading protocols and in each case the mean orientation of the collagen network reoriented to align in the direction of primary loading as expected. In summary, the ssQPLI system is effective at quantifying collagen fiber organization, and, when combined with mechanical testing, can rapidly provide pixel-wise measures of fiber orientation during biaxial loading.

Tractography of Porcine Meniscus Microstructure Using High-Resolution Diffusion Magnetic Resonance Imaging

To noninvasively evaluate the three-dimensional collagen fiber architecture of porcine meniscus using diffusion MRI, meniscal specimens were scanned using a 3D diffusion-weighted spin-echo pulse sequence at 7.0 T. The collagen fiber alignment was revealed in each voxel and the complex 3D collagen network was visualized for the entire meniscus using tractography. The proposed automatic segmentation methods divided the whole meniscus to different zones (Red-Red, Red-White, and White-White) and different parts (anterior, body, and posterior). The diffusion tensor imaging (DTI) metrics were quantified based on the segmentation results. The heatmap was generated to investigate the connections among different regions of meniscus. Strong zonal-dependent diffusion properties were demonstrated by DTI metrics. The fractional anisotropy (FA) value increased from 0.13 (White-White zone) to 0.26 (Red-Red zone) and the radial diffusivity (RD) value changed from 1.0 × 10^-3 mm²/s (White-White zone) to 0.7 × 10^-3 mm²/s (Red-Red zone). Coexistence of both radial and circumferential collagen fibers in the meniscus was evident by diffusion tractography. Weak connections were found between White-White zone and Red-Red zone in each part of the meniscus. The anterior part and posterior part were less connected, while the body part showed high connections to both anterior part and posterior part. The tractography based on diffusion MRI may provide a complementary method to study the integrity of meniscus and nondestructively visualize the 3D collagen fiber architecture.

Effects of Angular Resolution and b Value on Diffusion Tensor Imaging in Knee Joint

OBJECTIVE

To investigate the influences of the diffusion gradient directions (angular resolution) and the strength of the diffusion gradient (b value) on diffusion tensor imaging (DTI) metrics and tractography of various connective tissues in knee joint.

DESIGN

Two rat knee joints were scanned on a preclinical 9.4-T system using a 3-dimensional diffusion-weighted spin echo pulse sequence. One protocol with b value of 500, 1500, and 2500 s/mm² were acquired separately using 43 diffusion gradient directions. The other protocol with b value of 1000 s/mm² was performed using 147 diffusion gradient directions. The in-plane resolution was 45 µm isotropic. Fractional anisotropy (FA) and mean diffusivity (MD) were compared at different angular resolution. Tractography was quantitatively evaluated at different b values and angular resolutions in cartilage, ligament, meniscus, and growth plate.

RESULTS

The ligament showed higher FA value compared with growth plate and cartilage. The FA values were largely overestimated at the angular resolution of 6. Compared with FA, MD showed less sensitivity to the angular resolution. The fiber tracking was failed at low angular resolution (6 diffusion gradient directions) or high b value (2500 s/mm²). The quantitative measurements of tract length and track volume were strongly dependent on angular resolution and b value.

CONCLUSIONS

To obtain consistent DTI outputs and tractography in knee joint, the scan may require a proper b value (ranging from 500 to 1500 s/mm²) and sufficient angular resolution (>14) with signal-to-noise ratio >10.

Characterization of anisotropic T2W signals from human knee femoral cartilage: The magic angle effect on a spherical surface

The aim of the current study was to propose a generalized magic angle effect (gMAE) function for characterizing anisotropic T2W signals of human knee femoral cartilage with a spherical surface in clinical studies. A gMAE model function f(α, ε) was formulated for an orientation-dependent (ε) transverse T₂ (i.e., 1/R₂ ) relaxation in cartilage assuming an axially symmetric distribution (α) of collagen fibers. T2W sagittal images were acquired on an adult volunteer's healthy knee at 3 T, and ROI-based average signals S(ε) were extracted from angularly and radially segmented femoral cartilage. Compared with the standard MAE (sMAE) functions in the deep (DZ, α = 0°) and in the superficial (SZ, α = 90°) zones, a general form of R₂ orientation-dependent function f(α, ε) was fitted to S(ε), including an isotropic R₂ contribution (internal reference [REF]). Goodness of fit was evaluated by root-mean-square deviations (RMSDs). An F-test and a paired t-test were respectively used to assess significant differences between the observed variances and means, with statistical significance set to p less than .05. As a symmetric orientation-dependence function with a varying dynamic range, the proposed gMAE model outperformed the previous sMAE functions manifested by significantly reduced RMSDs in the DZ (0.239 ± 0.122 vs. 0.267 ± 0.097, p = .014) and in the SZ (0.183 ± 0.081 vs. 0.254 ± 0.085, p < .001). The fitted average angle α (38.5 ± 34.6° vs. 45.1 ± 30.1°, p < .43) and REF (5.092 ± 0.369 vs. 5.305 ± 0.440, p < .001) were smaller in the DZ than those in SZ, in good agreement with the reported collagen fibril microstructural configurations and the nonbound water contribution to R₂ in articular cartilage. In conclusion, a general form of the magic angle effect function was proposed and demonstrated for better characterizing anisotropic T2W signals from human knee femoral cartilage at 3 T in clinical studies.

Effect of b Value on Imaging Quality for Diffusion Tensor Imaging of the Spinal Cord at Ultrahigh Field Strength

Objective

To explore the optimal b value setting for diffusion tensor imaging of rats' spinal cord at ultrahigh field strength (7 T).

Methods

Spinal cord diffusion tensor imaging data were collected from 14 rats (5 healthy, 9 spinal cord injured) with a series of b values (200, 300, 400, 500, 600, 700, 800, 900, and 1000 s/mm²) under the condition that other scanning parameters were consistent. The image quality (including image signal-to-noise ratio and image distortion degree) and data quality (i.e., the stability and consistency of the DTI-derived parameters, referred to as data stability and data consistency) were quantitatively evaluated. The min-max normalization method was used to process the calculation results of the four indicators. Finally, the image and data quality under each b value were synthesized to determine the optimal b value.

Results

b = 200 s/mm² and b = 900 s/mm² ranked in the top two of the comprehensive evaluation, with the best image quality at b = 200 s/mm² and the best data quality at b = 900 s/mm².

Conclusion

Considering the shortcomings of the ability of low b values to reflect the microstructure, b = 900 s/mm² can be used as the optimal b value for 7 T spinal cord diffusion tensor scanning.