1
|
Lockard CA, Hooijmans MT, Zhou X, Coolbaugh C, Damon BM. The impact of diffusion tensor imaging tractography settings on muscle fascicle architecture and diffusion parameter estimates: Tract length, completion, and curvature are most sensitive to tractography settings. NMR IN BIOMEDICINE 2024:e5205. [PMID: 38967274 DOI: 10.1002/nbm.5205] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/16/2023] [Revised: 05/17/2024] [Accepted: 06/05/2024] [Indexed: 07/06/2024]
Abstract
Diffusion-tensor (DT)-MRI tractography provides information about properties relevant to muscle health and function, including estimates of architectural properties such as fascicle length, pennation angle, and curvature and diffusion properties such as mean diffusivity (MD) and fractional anisotropy (FA). Tractography settings, including integration algorithms, thresholds for early tract termination, and tract smoothing approaches, impact the accuracy of the muscle property estimates. However, muscle DT-MRI tractography is performed using a variety of these settings, complicating comparisons between different studies. The effects of different tractography settings on muscle architecture estimates have not been fully explored, and optimized settings for muscle tractography have not yet been determined. We examined the influence of integration algorithm and termination check settings combined with a range of step sizes, termination criteria, and smoothing polynomial orders on tract characteristics, completion/reason for termination, and goodness of fit between fiber tracts and smoothing polynomials using 3-T DT-MR images of the lower leg muscles of seven healthy adults. We found that tract length and completion were highly sensitive to strict FA and intersegment angle thresholds (25%-69% reduction in complete fiber tracts from lowest to highest minimum FA threshold and 11%-36% reduction from highest to lowest intersegment angle threshold). Higher order polynomials (third and fourth order vs. second order) better fit the muscle fiber trajectories, but curvature estimates were highly sensitive to smoothing polynomial order (3.9-6.6 m-1 increase for second- vs. fourth-order fitting polynomials). Step size impacted curvature estimates, albeit to a lesser degree. Integration algorithm had little impact, and mean pennation angle, and tract-based FA and MD, were relatively insensitive to all parameters. The results demonstrate which muscle diffusion measures and architectural estimates are most sensitive to varying tractography settings and support the need for consistent reporting of tractography details to aid interpretation and comparison of results between studies.
Collapse
Affiliation(s)
- Carly A Lockard
- Carle Clinical Imaging Research Program, Stephens Family Clinical Research Institute, Carle Health, Urbana, Illinois, USA
| | - Melissa T Hooijmans
- Carle Clinical Imaging Research Program, Stephens Family Clinical Research Institute, Carle Health, Urbana, Illinois, USA
- Department of Radiology and Nuclear Medicine, Amsterdam UMC, University of Amsterdam, Amsterdam, The Netherlands
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, Tennessee, USA
| | - Xingyu Zhou
- Carle Clinical Imaging Research Program, Stephens Family Clinical Research Institute, Carle Health, Urbana, Illinois, USA
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, Tennessee, USA
- Department of Biomedical Engineering, Vanderbilt University, Nashville, Tennessee, USA
| | - Crystal Coolbaugh
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, Tennessee, USA
| | - Bruce M Damon
- Carle Clinical Imaging Research Program, Stephens Family Clinical Research Institute, Carle Health, Urbana, Illinois, USA
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, Tennessee, USA
- Department of Biomedical Engineering, Vanderbilt University, Nashville, Tennessee, USA
- Department of Radiology and Radiological Sciences, Vanderbilt University, Nashville, Tennessee, USA
- Department of Bioengineering, University of Illinois at Urbana-Champaign, Urbana, Illinois, USA
- Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, Illinois, USA
| |
Collapse
|
2
|
Hooijmans MT, Lockard CA, Zhou X, Coolbaugh C, Pineda Guzman R, Kersh ME, Damon BM. A registration strategy to characterize DTI-observed changes in skeletal muscle architecture due to passive shortening. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.04.11.589123. [PMID: 38645028 PMCID: PMC11030449 DOI: 10.1101/2024.04.11.589123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/23/2024]
Abstract
Skeletal muscle architecture is a key determinant of muscle function. Architectural properties such as fascicle length, pennation angle, and curvature can be characterized using Diffusion Tensor Imaging (DTI), but acquiring these data during a contraction is not currently feasible. However, an image registration-based strategy may be able to convert muscle architectural properties observed at rest to their contracted state. As an initial step toward this long-term objective, the aim of this study was to determine if an image registration strategy could be used to convert the whole-muscle average architectural properties observed in the extended joint position to those of a flexed position, following passive rotation. DTI and high-resolution fat/water scans were acquired in the lower leg of seven healthy participants on a 3T MR system in +20° (plantarflexion) and -10° (dorsiflexion) foot positions. The diffusion and anatomical images from the two positions were used to propagate DTI fiber-tracts from seed points along a mesh representation of the aponeurosis of fiber insertion. The -10° and +20° anatomical images were registered and the displacement fields were used to transform the mesh and fiber-tracts from the +20° to the -10° position. Student's paired t-tests were used to compare the mean architectural parameters between the original and transformed fiber-tracts. The whole-muscle average fiber-tract length, pennation angle, curvature, and physiological cross-sectional areas estimates did not differ significantly. DTI fiber-tracts in plantarflexion can be transformed to dorsiflexion position without significantly affecting the average architectural characteristics of the fiber-tracts. In the future, a similar approach could be used to evaluate muscle architecture in a contracted state.
Collapse
Affiliation(s)
- Melissa T. Hooijmans
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, United States of America
- Carle Clinical Imaging Research Program, Stephens Family Clinical Research Institute, Carle Health, Urbana, IL, United States of America
- Department of Radiology and Nuclear Medicine, Amsterdam UMC, University of Amsterdam, Amsterdam, The Netherlands
| | - Carly A. Lockard
- Carle Clinical Imaging Research Program, Stephens Family Clinical Research Institute, Carle Health, Urbana, IL, United States of America
| | - Xingyu Zhou
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, United States of America
- Carle Clinical Imaging Research Program, Stephens Family Clinical Research Institute, Carle Health, Urbana, IL, United States of America
- Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, United States of America
| | - Crystal Coolbaugh
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, United States of America
| | - Roberto Pineda Guzman
- Carle Clinical Imaging Research Program, Stephens Family Clinical Research Institute, Carle Health, Urbana, IL, United States of America
| | - Mariana E. Kersh
- Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, United States of America
- Department of Biomedical and Translational Sciences, Carle-Illinois College of Medicine, University of Illinois at Urbana-Champaign, Urbana, IL, United States of America
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States of America
| | - Bruce M. Damon
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, United States of America
- Carle Clinical Imaging Research Program, Stephens Family Clinical Research Institute, Carle Health, Urbana, IL, United States of America
- Department of Biomedical and Translational Sciences, Carle-Illinois College of Medicine, University of Illinois at Urbana-Champaign, Urbana, IL, United States of America
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States of America
- Department of Radiology and Radiological Sciences, Vanderbilt University, Nashville, TN, United States of America
- Department of Bioengineering, University of Illinois at Urbana-Champaign, Urbana, IL, United States of America
| |
Collapse
|
3
|
Zhou X, Lockard CA, Hooijmans MT, Damon BM. Predicted effects of image acquisition and analysis conditions on DTMRI tractography-based muscle architecture estimates. Magn Reson Med 2024; 91:1337-1353. [PMID: 38044800 DOI: 10.1002/mrm.29910] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2023] [Revised: 10/15/2023] [Accepted: 10/16/2023] [Indexed: 12/05/2023]
Abstract
PURPOSE To quantify the effects of the intrinsic signal pattern, image acquisition conditions, and data analysis conditions on diffusion-tensor MRI (DTMRI) tractography-based muscle architecture estimates using a sampling-reconstruction assessment framework. METHODS Numerical models of muscles were constructed with realistic architectural properties. DTMRI signals were computed at signal-to-noise ratio (SNR) of 24-96 and common voxel sizes. Fiber tracking was performed, and the results were compared with the known architectural properties. RESULTS SNR exerted the most significant impact on the outcome. The outcome variables approached asymptotes at SNR ≈ 54. Large in-plane voxel dimensions reduced the similarity between reconstructed fibers and the known architectural properties. Higher order polynomials helped reconstruct fibers with more complicated geometry but overfit noise for less complex geometries. The intrinsic fiber curvature also affected the robustness of polynomial smoothing to SNR. Other conditions, such as the fiber dimensionality, voxel aspect ratio, and slice thickness, did not affect the outcomes. CONCLUSION SNR ≥ 54 is recommended for accurate muscle architecture characterization using DTMRI. Averaged across all simulated conditions, the greatest percent errors under SNR = 54 were -5.6% and -4.0% for the pennation angle and fiber-tract length estimates, respectively. For fiber tracts with intermediate intrinsic curvature, the greatest percent error for the curvature estimate was 9.8% for SNR = 54. Smaller in-plane voxel size (≤1.5 mm) is preferred to minimize the estimation error in architectural properties. If necessary, slice thickness may be adjusted within typical ranges to achieve sufficient SNR when slices are aligned near the fiber direction. Third-order polynomial fitting is appropriate for smoothing fiber tracts.
Collapse
Affiliation(s)
- Xingyu Zhou
- Carle Clinical Imaging Research Program, Stephens Family Clinical Research Institute, Carle Health, Urbana, Illinois, USA
- Department of Biomedical Engineering, Vanderbilt University, Nashville, Tennessee, USA
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, Tennessee, USA
| | - Carly A Lockard
- Carle Clinical Imaging Research Program, Stephens Family Clinical Research Institute, Carle Health, Urbana, Illinois, USA
| | - Melissa T Hooijmans
- Department of Radiology and Nuclear Medicine, Amsterdam UMC, University of Amsterdam, Amsterdam Movement Sciences, Amsterdam, The Netherlands
| | - Bruce M Damon
- Carle Clinical Imaging Research Program, Stephens Family Clinical Research Institute, Carle Health, Urbana, Illinois, USA
- Department of Biomedical Engineering, Vanderbilt University, Nashville, Tennessee, USA
- Radiology and Radiological Sciences, Vanderbilt University, Nashville, Tennessee, USA
- Department of Bioengineering, Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, Illinois, USA
- Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, Illinois, USA
| |
Collapse
|
4
|
Damon BM, Ding Z, Hooijmans MT, Anderson AW, Zhou X, Coolbaugh CL, George MK, Landman BA. A MATLAB toolbox for muscle diffusion-tensor MRI tractography. J Biomech 2021; 124:110540. [PMID: 34171675 DOI: 10.1016/j.jbiomech.2021.110540] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2021] [Revised: 05/17/2021] [Accepted: 05/31/2021] [Indexed: 11/26/2022]
Abstract
Diffusion-tensor MRI fiber tractography has been used to reconstruct skeletal muscle architecture, but remains a specialized technique using custom-written data processing routines. In this work, we describe the public release of a software toolbox having the following design objectives: accomplish the pre-processing tasks of file input, image registration, denoising, and diffusion-tensor calculation; allow muscle-specific methods for defining seed points; make fiber-tract architectural measurements referenced to tendinous structures; visualize fiber tracts and other muscle structures of interest; analyze the goodness of outcomes; and provide a programming structure that allows the addition of new capabilities in future versions. The proper function of the code was verified using simulated datasets. The toolbox capabilities for characterizing human muscle structure in vivo were demonstrated in a case study. These capabilities included measurements of muscle morphology; contractile and non-contractile tissue volumes; fiber-tract length, pennation angle, curvature; and the physiological cross-sectional area,. The free public release of this software is a first step in creating of a community of users who use these tools in studies of muscle physiology and biomechanics. Users may further contribute to code development. Along with simulated and actual datasets for benchmarking, these tools will further create mechanisms for enhancing scientific rigor and developing and validating new code features. Planned future developments include additional options for image pre-processing, development of a graphical user interface, analysis of architectural patterns during muscle contraction, and integration of functional imaging data.
Collapse
Affiliation(s)
- Bruce M Damon
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, USA; Departments of Radiology and Radiological Sciences, Vanderbilt University, Nashville, TN 37232, USA; Departments of Biomedical Engineering, Vanderbilt University, Nashville, TN 37232, USA; Departments of Molecular Physiology and Biophysics, Vanderbilt University, Nashville, TN 37232, USA.
| | - Zhaohua Ding
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, USA; Departments of Radiology and Radiological Sciences, Vanderbilt University, Nashville, TN 37232, USA; Departments of Biomedical Engineering, Vanderbilt University, Nashville, TN 37232, USA; Departments of Electrical Engineering and Computer Science, Vanderbilt University, Nashville, TN 37232, USA
| | - Melissa T Hooijmans
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, USA
| | - Adam W Anderson
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, USA; Departments of Radiology and Radiological Sciences, Vanderbilt University, Nashville, TN 37232, USA; Departments of Biomedical Engineering, Vanderbilt University, Nashville, TN 37232, USA; Departments of Physics and Astronomy, Vanderbilt University, Nashville, TN 37232, USA
| | - Xingyu Zhou
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, USA; Departments of Biomedical Engineering, Vanderbilt University, Nashville, TN 37232, USA
| | - Crystal L Coolbaugh
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, USA
| | - Mark K George
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, USA
| | - Bennett A Landman
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, USA; Departments of Radiology and Radiological Sciences, Vanderbilt University, Nashville, TN 37232, USA; Departments of Biomedical Engineering, Vanderbilt University, Nashville, TN 37232, USA; Departments of Electrical Engineering and Computer Science, Vanderbilt University, Nashville, TN 37232, USA
| |
Collapse
|
5
|
Sahrmann AS, Stott NS, Besier TF, Fernandez JW, Handsfield GG. Soleus muscle weakness in cerebral palsy: Muscle architecture revealed with Diffusion Tensor Imaging. PLoS One 2019; 14:e0205944. [PMID: 30802250 PMCID: PMC6388915 DOI: 10.1371/journal.pone.0205944] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2018] [Accepted: 01/24/2019] [Indexed: 11/28/2022] Open
Abstract
Cerebral palsy (CP) is associated with movement disorders and reduced muscle size. This latter phenomenon has been observed by computing muscle volumes from conventional MRI, with most studies reporting significantly reduced volumes in leg muscles. This indicates impaired muscle growth, but without knowing muscle fiber orientation, it is not clear whether muscle growth in CP is impaired in the along-fiber direction (indicating shortened muscles and limited range of motion) or the cross-fiber direction (indicating weak muscles and impaired strength). Using Diffusion Tensor Imaging (DTI) we can determine muscle fiber orientation and construct 3D muscle architectures which can be used to examine both along-fiber length and cross-sectional area. Such an approach has not been undertaken in CP. Here, we use advanced DTI sequences with fast imaging times to capture fiber orientations in the soleus muscle of children with CP and age-matched, able-bodied controls. Cross sectional areas perpendicular to the muscle fiber direction were reduced (37 ± 11%) in children with CP compared to controls, indicating impaired muscle strength. Along-fiber muscle lengths were not different between groups. This study is the first to demonstrate along-fiber and cross-fiber muscle architecture in CP using DTI and implicates impaired cross-sectional muscle growth in children with cerebral palsy.
Collapse
Affiliation(s)
- Annika S. Sahrmann
- Auckland Bioengineering Institute, University of Auckland, Auckland, New Zealand
| | - Ngaire Susan Stott
- Department of Orthopaedic Surgery, Faculty of Medical and Health Sciences, University of Auckland, Auckland, New Zealand
| | - Thor F. Besier
- Auckland Bioengineering Institute, University of Auckland, Auckland, New Zealand
- Department of Engineering Science, Faculty of Engineering, University of Auckland, Auckland, New Zealand
| | - Justin W. Fernandez
- Auckland Bioengineering Institute, University of Auckland, Auckland, New Zealand
- Department of Engineering Science, Faculty of Engineering, University of Auckland, Auckland, New Zealand
| | - Geoffrey G. Handsfield
- Auckland Bioengineering Institute, University of Auckland, Auckland, New Zealand
- * E-mail:
| |
Collapse
|
6
|
Damon BM, Froeling M, Buck AKW, Oudeman J, Ding Z, Nederveen AJ, Bush EC, Strijkers GJ. Skeletal muscle diffusion tensor-MRI fiber tracking: rationale, data acquisition and analysis methods, applications and future directions. NMR IN BIOMEDICINE 2017; 30:10.1002/nbm.3563. [PMID: 27257975 PMCID: PMC5136336 DOI: 10.1002/nbm.3563] [Citation(s) in RCA: 62] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/10/2015] [Revised: 03/19/2016] [Accepted: 04/27/2016] [Indexed: 05/21/2023]
Abstract
The mechanical functions of muscles involve the generation of force and the actuation of movement by shortening or lengthening under load. These functions are influenced, in part, by the internal arrangement of muscle fibers with respect to the muscle's mechanical line of action. This property is known as muscle architecture. In this review, we describe the use of diffusion tensor (DT)-MRI muscle fiber tracking for the study of muscle architecture. In the first section, the importance of skeletal muscle architecture to function is discussed. In addition, traditional and complementary methods for the assessment of muscle architecture (brightness-mode ultrasound imaging and cadaver analysis) are presented. Next, DT-MRI is introduced and the structural basis for the reduced and anisotropic diffusion of water in muscle is discussed. The third section discusses issues related to the acquisition of skeletal muscle DT-MRI data and presents recommendations for optimal strategies. The fourth section discusses methods for the pre-processing of DT-MRI data, the available approaches for the calculation of the diffusion tensor and the seeding and propagating of fiber tracts, and the analysis of the tracking results to measure structural properties pertinent to muscle biomechanics. Lastly, examples are presented of how DT-MRI fiber tracking has been used to provide new insights into how muscles function, and important future research directions are highlighted. Copyright © 2016 John Wiley & Sons, Ltd.
Collapse
Affiliation(s)
- Bruce M. Damon
- Institute of Imaging Science, Vanderbilt University, Nashville TN USA
- Department of Radiology and Radiological Sciences, Vanderbilt University, Nashville TN USA
- Department of Biomedical Engineering, Vanderbilt University, Nashville TN USA
- Department of Molecular Physiology and Biophysics, Vanderbilt University, Nashville TN USA
| | - Martijn Froeling
- Department of Radiology, University Medical Center, Utrecht, the Netherlands
| | - Amanda K. W. Buck
- Institute of Imaging Science, Vanderbilt University, Nashville TN USA
- Department of Radiology and Radiological Sciences, Vanderbilt University, Nashville TN USA
- Department of Biomedical Engineering, Vanderbilt University, Nashville TN USA
| | - Jos Oudeman
- Department of Radiology, Academic Medical Center, Amsterdam, the Netherlands
| | - Zhaohua Ding
- Institute of Imaging Science, Vanderbilt University, Nashville TN USA
- Department of Radiology and Radiological Sciences, Vanderbilt University, Nashville TN USA
- Department of Biomedical Engineering, Vanderbilt University, Nashville TN USA
- Department of Electrical Engineering and Computer Engineering, Vanderbilt University, Nashville TN USA
| | - Aart J. Nederveen
- Department of Radiology, Academic Medical Center, Amsterdam, the Netherlands
| | - Emily C. Bush
- Institute of Imaging Science, Vanderbilt University, Nashville TN USA
| | - Gustav J. Strijkers
- Department of Biomedical Engineering and Physics, Academic Medical Center, Amsterdam, the Netherlands
| |
Collapse
|