Publications

Also see our publication highlights for a featured list.

Publications

  • Circumventing the curse of dimensionality in magnetic resonance fingerprinting through a deep learning approach. NMR in biomedicine Barbieri, M., Lee, P. K., Brizi, L., Giampieri, E., Solera, F., Castellani, G., Hargreaves, B. A., Testa, C., Lodi, R., Remondini, D. 1800: e4670

    Abstract

    Magnetic resonance fingerprinting (MRF) is a rapidly developing approach for fast quantitative MRI. A typical drawback of dictionary-based MRF is an explosion of the dictionary size as a function of the number of reconstructed parameters, according to the "curse of dimensionality", which determines an explosion of resource requirements. Neural networks (NNs) have been proposed as a feasible alternative, but this approach is still in its infancy. In this work, we design a deep learning approach to MRF using a fully connected network (FCN). In the first part we investigate, by means of simulations, how the NN performance scales with the number of parameters to be retrieved in comparison with the standard dictionary approach. Four MRF sequences were considered: IR-FISP, bSSFP, IR-FISP-B1 , and IR-bSSFP-B1 , the latter two designed to be more specific for B 1 + parameter encoding. Estimation accuracy, memory usage, and computational time required to perform the estimation task were considered to compare the scalability capabilities of the dictionary-based and the NN approaches. In the second part we study optimal training procedures by including different data augmentation and preprocessing strategies during training to achieve better accuracy and robustness to noise and undersampling artifacts. The study is conducted using the IR-FISP MRF sequence exploiting both simulations and in vivo acquisitions. Results demonstrate that the NN approach outperforms the dictionary-based approach in terms of scalability capabilities. Results also allow us to heuristically determine the optimal training strategy to make an FCN able to predict T1 , T2 , and M0 maps that are in good agreement with those obtained with the original dictionary approach. k-SVD denoising is proposed and found to be critical as a preprocessing step to handle undersampled data.

    View details for DOI 10.1002/nbm.4670

    View details for PubMedID 35088466

  • Volumetric and multispectral DWI near metallic implants using a non-linear phase Carr-Purcell-Meiboom-Gill diffusion preparation. Magnetic resonance in medicine Lee, P. K., Yoon, D., Sandberg, J. K., Vasanawala, S. S., Hargreaves, B. A. 1800

    Abstract

    PURPOSE: DWI near metal implants has not been widely explored due to substantial challenges associated with through-slice and in-plane distortions, the increased encoding requirement of different spectral bins, and limited SNR. There is no widely adopted clinical protocol for DWI near metal since the commonly used EPI trajectory fails completely due to distortion from extreme off-resonance ranging from 2 to 20 kHz. We present a sequence that achieves DWI near metal with moderate b-values (400-500 s/mm2 ) and volumetric coverage in clinically feasible scan times.THEORY AND METHODS: Multispectral excitation with Cartesian sampling, view angle tilting, and kz phase encoding reduce in-plane and through-plane off-resonance artifacts, and Carr-Purcell-Meiboom-Gill (CPMG) spin-echo refocusing trains counteract T2* effects. The effect of random phase on the refocusing train is eliminated using a stimulated echo diffusion preparation. Root-flipped Shinnar-Le Roux refocusing pulses permits preparation of a high spectral bandwidth, which improves imaging times by reducing the number of excitations required to cover the desired spectral range. B1 sensitivity is reduced by using an excitation that satisfies the CPMG condition in the preparation. A method for ADC quantification insensitive to background gradients is presented.RESULTS: Non-linear phase refocusing pulses reduces the peak B1 by 46% which allows RF bandwidth to be doubled. Simulations and phantom experiments show that a non-linear phase CPMG pulse pair reduces B1 sensitivity. Application in vivo demonstrates complementary contrast to conventional multispectral acquisitions and improved visualization compared to DW-EPI.CONCLUSION: Volumetric and multispectral DW imaging near metal can be achieved with a 3D encoded sequence.

    View details for DOI 10.1002/mrm.29153

    View details for PubMedID 35014729

  • Biopsy Marker Localization with Thermo-Acoustic Ultrasound for Lumpectomy Guidance. Medical physics Dixit, N., Daniel, B. L., Hargreaves, B. A., Pauly, J. M., Scott, G. C. 2021

    Abstract

    PURPOSE: Almost one in four lumpectomies fails to fully remove cancerous tissue from the breast, requiring reoperation. This high failure rate suggests that existing lumpectomy guidance methods are inadequate for allowing surgeons to consistently identify the proper volume of tissue for excision. Current guidance techniques either provide little information about the tumor position or require surgeons to frequently switch between making incisions and manually probing for a marker placed at the lesion site. This article explores the feasibility of thermo-acoustic ultrasound (TAUS) to enable hands-free localization of metallic biopsy markers throughout surgery, which would allow for continuous visualization of the lesion site in the breast without the interruption of surgery. In a TAUS-based localization system, microwave excitations would be transmitted into the breast, and the amplification in microwave absorption around the metallic markers would generate acoustic signals from the marker sites through the thermo-acoustic effect. Detection and ranging of these signals by multiple acoustic receivers on the breast could then enable marker localization through acoustic multilateration.METHODS: Physics simulations were used to characterize the TAUS signals generated from different markers by microwave excitations. First, electromagnetic simulations determined the spatial pattern of the amplification in microwave absorption around the markers. Then, acoustic simulations characterized the acoustic fields generated from these markers at various acoustic frequencies. TAUS-based one-dimensional (1D) ranging of two metallic markers - including a biopsy marker that is FDA-approved for clinical use - immersed in saline was also performed using a bench-top setup. To perform TAUS acquisitions, a microwave applicator was driven by 2.66 GHz microwave signals that were amplitude-modulated by chirps at the desired acoustic excitation frequencies, and the resulting TAUS signal from the markers was detected by an ultrasonic transducer.RESULTS: The simulation results show that the geometry of the marker strongly impacts the quantity and spatial pattern of both the microwave absorption around the marker and the resulting TAUS signal generated from the marker. The simulated TAUS signal maps and acoustic frequency responses also make clear that the marker geometry plays an important role in determining the overall system response. Using the bench-top setup, TAUS detection and 1D localization of the markers was successfully demonstrated for multiple different combinations of microwave applicator and metallic marker. These initial results indicate that TAUS-based localization of biopsy markers is feasible.CONCLUSIONS: Through microwave excitations and acoustic detection, TAUS can be used to localize metallic biopsy markers. With further development, TAUS opens new avenues to enable a more intuitive lumpectomy guidance system that could help to achieve better lumpectomy outcomes.

    View details for DOI 10.1002/mp.15115

    View details for PubMedID 34287972

  • Clinical utility of accelerated MAVRIC-SL with robust-PCA compared to conventional MAVRIC-SL in evaluation of total hip arthroplasties. Skeletal radiology Doyle, Z., Yoon, D., Lee, P. K., Rosenberg, J., Hargreaves, B. A., Beaulieu, C. F., Stevens, K. J. 2021

    Abstract

    OBJECTIVE: To compare the diagnostic performance of a conventional metal artifact suppression sequence MAVRIC-SL (multi-acquisition variable-resonance image combination selective) and a novel 2.6-fold faster sequence employing robust principal component analysis (RPCA), in the MR evaluation of hip implants at 3T.MATERIALS AND METHODS: Thirty-six total hip implants in 25 patients were scanned at 3T using a conventional MAVRIC-SL proton density-weighted sequence and an RPCA MAVRIC-SL proton density-weighted sequence. Comparison was made of image quality, geometric distortion, visualization around acetabular and femoral components, and conspicuity of abnormal imaging findings using the Wilcoxon signed-rank test and a non-inferiority test. Abnormal findings were correlated with subsequent clinical management and intraoperative findings if the patient underwent subsequent surgery.RESULTS: Mean scores for conventional MAVRIC-SL were better than RPCA MAVRIC-SL for all qualitative parameters (p<0.05), although the probability of RPCA MAVRIC-SL being clinically useful was non-inferior to conventional MAVRIC-SL (within our accepted 10% difference, p<0.05), except for visualization around the acetabular component. Abnormal imaging findings were seen in 25 hips, and either equally visible or visible but less conspicuous on RPCA MAVRIC-SL in 21 out of 25 cases. In 4 cases, a small joint effusion was queried on MAVRIC-SL but not RPCA MAVRIC-SL, but the presence or absence of a small effusion did not affect subsequent clinical management and patient outcome.CONCLUSION: While the overall image quality is reduced, RPCA MAVRIC-SL allows for significantly reduced scan time and maintains almost equal diagnostic performance.

    View details for DOI 10.1007/s00256-021-03848-y

    View details for PubMedID 34223946

  • Zero echo time pediatric musculoskeletal magnetic resonance imaging: initial experience. Pediatric radiology Sandberg, J. K., Young, V. A., Yuan, J., Hargreaves, B. A., Wishah, F., Vasanawala, S. S. 2021

    Abstract

    BACKGROUND: Projection radiography (XR) is often supplemented by both CT (to evaluate osseous structures with ionizing radiation) and MRI (for marrow and soft-tissue assessment). Zero echo time (ZTE) MR imaging produces a "CT-like" osseous contrast that might obviate CT.OBJECTIVE: This study investigated our institution's initial experience in implementing an isotropic ZTE MR imaging sequence for pediatric musculoskeletal examinations.MATERIALS AND METHODS: Pediatric patients referred for extremity MRI at 3tesla (T) underwent ZTE MR imaging to yield images with contrast similar to that of CT. A radiograph-like image was also created with ray-sum image processing. We assessed ZTE-CT/XR anatomical image quality (Sanat) from 0 (nondiagnostic) to 5 (outstanding). Further, we made image comparisons on a 5-point scale (Scomp) (range of -2 = conventional CT/XR greater anatomical delineation to +2 = ZTE-CT/XR greater anatomical delineation; 0=same) for three cohorts: (1) ZTE-XR to conventional radiography, (2) ZTE-CT to conventional CT and (3) pathological lesion assessment on ZTE-XR to conventional radiography. We measured cortical thickness of ZTE-XR and ZTE-CT and compared these with conventional imaging. We calculated confidence interval of proportions, Wilcoxon rank sum test and intraclass correlation coefficients for inter-reader agreement.RESULTS: Cohorts 1, 2 and 3 consisted of 40, 20 and 35 cases, respectively (age range 0.6-23.0years). ZTE-CT versus CT and ZTE-XR versus radiography of cortical thicknesses were not significantly different (P=0.55 and P=0.31, respectively). Cortical delineation was rated diagnostic or better (score of 3, 4 or 5) in all cases (confidence interval of proportions = 100%) for ZTE-CT/XR. Similarly, intramedullary cavity delineation was rated diagnostic or better in all cases for ZTE-CT, and ZTE-XR was at least diagnostic in 58-63% of cases. For cohort 2, cortex and intramedullary cavity Scomp for ZTE-CT was comparable to those of conventional CT, with confidence interval of proportion (sum of score of -1 to +2) of 93-100% and 95%, respectively. Pathology visualized on ZTE-CT/XR was comparable; Scomp confidence interval of proportions was 95%/97-100%, with improved delineation of non-displaced fractures on ZTE-XR. Readers had moderate to near-perfect intraclass correlation coefficient (range=0.60-0.93).CONCLUSION: Implementation of a diagnostic-quality ZTE MRI sequence in the pediatric population is feasible and can be performed as a complementary pulse sequence to enhance musculoskeletal MRI studies. Compared to conventional CT, ZTE has comparable cortical delineation, intramedullary cavity and pathology visualization. While not intended as a replacement for conventional radiography, ZTE-XR provides similar visualization of pathology.

    View details for DOI 10.1007/s00247-021-05125-5

    View details for PubMedID 34156504

  • Non-contrast MRI of synovitis in the knee using quantitative DESS. European radiology Thoenen, J., Stevens, K. J., Turmezei, T. D., Chaudhari, A., Watkins, L. E., McWalter, E. J., Hargreaves, B. A., Gold, G. E., MacKay, J. W., Kogan, F. 2021

    Abstract

    OBJECTIVES: To determine whether synovitis graded by radiologists using hybrid quantitative double-echo in steady-state (qDESS) images can be utilized as a non-contrast approach to assess synovitis in the knee, compared against the reference standard of contrast-enhanced MRI (CE-MRI).METHODS: Twenty-two knees (11 subjects) with moderate to severe osteoarthritis (OA) were scanned using CE-MRI, qDESS with a high diffusion weighting (qDESSHigh), and qDESS with a low diffusion weighting (qDESSLow). Four radiologists graded the overall impression of synovitis, their diagnostic confidence, and regional grading of synovitis severity at four sites (suprapatellar pouch, intercondylar notch, and medial and lateral peripatellar recesses) in the knee using a 4-point scale. Agreement between CE-MRI and qDESS, inter-rater agreement, and intra-rater agreement were assessed using a linearly weighted Gwet's AC2.RESULTS: Good agreement was seen between CE-MRI and both qDESSLow (AC2=0.74) and qDESSHigh (AC2=0.66) for the overall impression of synovitis, but both qDESS sequences tended to underestimate the severity of synovitis compared to CE-MRI. Good inter-rater agreement was seen for both qDESS sequences (AC2=0.74 for qDESSLow, AC2=0.64 for qDESSHigh), and good intra-rater agreement was seen for both sequences as well (qDESSLow AC2=0.78, qDESSHigh AC2=0.80). Diagnostic confidence was moderate to high for qDESSLow (mean=2.36) and slightly less than moderate for qDESSHigh (mean=1.86), compared to mostly high confidence for CE-MRI (mean=2.73).CONCLUSIONS: qDESS shows potential as an alternative MRI technique for assessing the severity of synovitis without the use of a gadolinium-based contrast agent.KEY POINTS: The use of the quantitative double-echo in steady-state (qDESS) sequence for synovitis assessment does not require the use of a gadolinium-based contrast agent. Preliminary results found that low diffusion-weighted qDESS (qDESSLow) shows good agreement to contrast-enhanced MRI for characterization of the severity of synovitis, with a relative bias towards underestimation of severity. Preliminary results also found that qDESSLow shows good inter- and intra-rater agreement for the depiction of synovitis, particularly for readers experienced with the sequence.

    View details for DOI 10.1007/s00330-021-08025-2

    View details for PubMedID 33993332

  • The International Workshop on Osteoarthritis Imaging Knee MRI Segmentation Challenge: A Multi-Institute Evaluation and Analysis Framework on a Standardized Dataset. Radiology. Artificial intelligence Desai, A. D., Caliva, F., Iriondo, C., Mortazi, A., Jambawalikar, S., Bagci, U., Perslev, M., Igel, C., Dam, E. B., Gaj, S., Yang, M., Li, X., Deniz, C. M., Juras, V., Regatte, R., Gold, G. E., Hargreaves, B. A., Pedoia, V., Chaudhari, A. S., IWOAI Segmentation Challenge Writing Group 2021; 3 (3): e200078

    Abstract

    Purpose: To organize a multi-institute knee MRI segmentation challenge for characterizing the semantic and clinical efficacy of automatic segmentation methods relevant for monitoring osteoarthritis progression.Materials and Methods: A dataset partition consisting of three-dimensional knee MRI from 88 retrospective patients at two time points (baseline and 1-year follow-up) with ground truth articular (femoral, tibial, and patellar) cartilage and meniscus segmentations was standardized. Challenge submissions and a majority-vote ensemble were evaluated against ground truth segmentations using Dice score, average symmetric surface distance, volumetric overlap error, and coefficient of variation on a holdout test set. Similarities in automated segmentations were measured using pairwise Dice coefficient correlations. Articular cartilage thickness was computed longitudinally and with scans. Correlation between thickness error and segmentation metrics was measured using the Pearson correlation coefficient. Two empirical upper bounds for ensemble performance were computed using combinations of model outputs that consolidated true positives and true negatives.Results: Six teams (T 1-T 6) submitted entries for the challenge. No differences were observed across any segmentation metrics for any tissues (P = .99) among the four top-performing networks (T 2, T 3, T 4, T 6). Dice coefficient correlations between network pairs were high (> 0.85). Per-scan thickness errors were negligible among networks T 1-T 4 (P = .99), and longitudinal changes showed minimal bias (< 0.03 mm). Low correlations (rho < 0.41) were observed between segmentation metrics and thickness error. The majority-vote ensemble was comparable to top-performing networks (P = .99). Empirical upper-bound performances were similar for both combinations (P = .99).Conclusion: Diverse networks learned to segment the knee similarly, where high segmentation accuracy did not correlate with cartilage thickness accuracy and voting ensembles did not exceed individual network performance.See also the commentary by Elhalawani and Mak in this issue.Keywords: Cartilage, Knee, MR-Imaging, Segmentation © RSNA, 2020Supplemental material is available for this article.

    View details for DOI 10.1148/ryai.2021200078

    View details for PubMedID 34235438

  • Detection of knee synovitis using non-contrast-enhanced qDESS compared with contrast-enhanced MRI. Arthritis research & therapy de Vries, B. A., Breda, S. J., Sveinsson, B., McWalter, E. J., Meuffels, D. E., Krestin, G. P., Hargreaves, B. A., Gold, G. E., Oei, E. H. 2021; 23 (1): 55

    Abstract

    BACKGROUND: To assess diagnostic accuracy of quantitative double-echo in steady-state (qDESS) MRI for detecting synovitis in knee osteoarthritis (OA).METHODS: Patients with different degrees of radiographic knee OA were included prospectively. All underwent MRI with both qDESS and contrast-enhanced T1-weighted magnetic resonance imaging (CE-MRI). A linear combination of the two qDESS images can be used to create an image that displays contrast between synovium and the synovial fluid. Synovitis on both qDESS and CE-MRI was assessed semi-quantitatively, using a whole-knee synovitis sum score, indicating no/equivocal, mild, moderate, and severe synovitis. The correlation between sum scores of qDESS and CE-MRI (reference standard) was determined using Spearman's rank correlation coefficient and intraclass correlation coefficient for absolute agreement. Receiver operating characteristic analysis was performed to assess the diagnostic performance of qDESS for detecting different degrees of synovitis, with CE-MRI as reference standard.RESULTS: In the 31 patients included, very strong correlation was found between synovitis sum scores on qDESS and CE-MRI (rho=0.96, p<0.001), with high absolute agreement (0.84 (95%CI 0.14-0.95)). Mean sum score (SD) values on qDESS 5.16 (3.75) were lower than on CE-MRI 7.13 (4.66), indicating systematically underestimated synovitis severity on qDESS. For detecting mild synovitis or higher, high sensitivity and specificity were found for qDESS (1.00 (95%CI 0.80-1.00) and 0.909 (0.571-1.00), respectively). For detecting moderate synovitis or higher, sensitivity and specificity were good (0.727 (95%CI 0.393-0.927) and 1.00 (0.800-1.00), respectively).CONCLUSION: qDESS MRI is able to, however with an underestimation, detect synovitis in patients with knee OA.

    View details for DOI 10.1186/s13075-021-02436-8

    View details for PubMedID 33581741

  • Characterizing the transient response of knee cartilage to running: Decreases in cartilage T2 of female recreational runners. Journal of orthopaedic research : official publication of the Orthopaedic Research Society Crowder, H. A., Mazzoli, V. n., Black, M. S., Watkins, L. E., Kogan, F. n., Hargreaves, B. A., Levenston, M. E., Gold, G. E. 2021

    Abstract

    Cartilage transmits and redistributes biomechanical loads in the knee joint during exercise. Exercise-induced loading alters cartilage hydration and is detectable using quantitative MRI, where T2 relaxation time (T2 ) is influenced by cartilage collagen composition, fiber orientation, and changes in extracellular matrix. This study characterized short-term transient responses of healthy knee cartilage to running-induced loading using bilateral scans and image registration. Eleven healthy female recreational runners (33.73±4.22 years) and four healthy female controls (27.25±1.38 years) were scanned on a 3T GE MRI scanner with qDESS before running over-ground (runner group) or resting (control group) for 40 minutes. Subjects were scanned immediately post-activity at five-minute intervals for 60 minutes. T2 times were calculated for femoral, tibial, and patellar cartilage at each time point and analyzed using a mixed-effects model and Bonferroni post-hoc. There were immediate decreases in T2 (mean±SEM) post-run in superficial femoral cartilage of at least 3.3±0.3% (P=0.002) between baseline and Time 0 that remained for 25 minutes, a decrease in superficial tibial cartilage T2 of 2.9±0.4% (P=0.041) between baseline and Time 0, and a decrease in superficial patellar cartilage T2 of 3.6±0.3% (P=0.020) 15 minutes post-run. There were decreases in the medial posterior region of superficial femoral cartilage T2 of at least 5.3±0.2% (P=0.022) within five minutes post-run that remained at 60 minutes post-run. Clinical Significance: These results increase understanding of transient responses of healthy cartilage to repetitive, exercise-induced loading and establish preliminary recommendations for future definitive studies of cartilage response to running. This article is protected by copyright. All rights reserved.

    View details for DOI 10.1002/jor.24994

    View details for PubMedID 33483997

  • Diffusion Tensor Imaging of Skeletal Muscle Contraction Using Oscillating Gradient Spin Echo. Frontiers in neurology Mazzoli, V. n., Moulin, K. n., Kogan, F. n., Hargreaves, B. A., Gold, G. E. 2021; 12: 608549

    Abstract

    Diffusion tensor imaging (DTI) measures water diffusion in skeletal muscle tissue and allows for muscle assessment in a broad range of neuromuscular diseases. However, current DTI measurements, typically performed using pulsed gradient spin echo (PGSE) diffusion encoding, are limited to the assessment of non-contracted musculature, therefore providing limited insight into muscle contraction mechanisms and contraction abnormalities. In this study, we propose the use of an oscillating gradient spin echo (OGSE) diffusion encoding strategy for DTI measurements to mitigate the effect of signal voids in contracted muscle and to obtain reliable diffusivity values. Two OGSE sequences with encoding frequencies of 25 and 50 Hz were tested in the lower leg of five healthy volunteers with relaxed musculature and during active dorsiflexion and plantarflexion, and compared with a conventional PGSE approach. A significant reduction of areas of signal voids using OGSE compared with PGSE was observed in the tibialis anterior for the scans obtained in active dorsiflexion and in the soleus during active plantarflexion. The use of PGSE sequences led to unrealistically elevated axial diffusivity values in the tibialis anterior during dorsiflexion and in the soleus during plantarflexion, while the corresponding values obtained using the OGSE sequences were significantly reduced. Similar findings were seen for radial diffusivity, with significantly higher diffusivity measured in plantarflexion in the soleus muscle using the PGSE sequence. Our preliminary results indicate that DTI with OGSE diffusion encoding is feasible in human musculature and allows to quantitatively assess diffusion properties in actively contracting skeletal muscle. OGSE holds great potential to assess microstructural changes occurring in the skeletal muscle during contraction, and for non-invasive assessment of contraction abnormalities in patients with muscle diseases.

    View details for DOI 10.3389/fneur.2021.608549

    View details for PubMedID 33658976

    View details for PubMedCentralID PMC7917051