Breast Multiparametric MRI for Prediction of Neoadjuvant Chemotherapy Response in Breast Cancer: The BMMR2 Challenge.

Purpose To describe the design, conduct, and results of the Breast Multiparametric MRI for prediction of neoadjuvant chemotherapy Response (BMMR2) challenge. Materials and Methods The BMMR2 computational challenge opened on May 28, 2021, and closed on December 21, 2021. The goal of the challenge was to identify image-based markers derived from multiparametric breast MRI, including diffusion-weighted imaging (DWI) and dynamic contrast-enhanced (DCE) MRI, along with clinical data for predicting pathologic complete response (pCR) following neoadjuvant treatment. Data included 573 breast MRI studies from 191 women (mean age [±SD], 48.9 years ± 10.56) in the I-SPY 2/American College of Radiology Imaging Network (ACRIN) 6698 trial (ClinicalTrials.gov: NCT01042379). The challenge cohort was split into training (60%) and test (40%) sets, with teams blinded to test set pCR outcomes. Prediction performance was evaluated by area under the receiver operating characteristic curve (AUC) and compared with the benchmark established from the ACRIN 6698 primary analysis. Results Eight teams submitted final predictions. Entries from three teams had point estimators of AUC that were higher than the benchmark performance (AUC, 0.782 [95% CI: 0.670, 0.893], with AUCs of 0.803 [95% CI: 0.702, 0.904], 0.838 [95% CI: 0.748, 0.928], and 0.840 [95% CI: 0.748, 0.932]). A variety of approaches were used, ranging from extraction of individual features to deep learning and artificial intelligence methods, incorporating DCE and DWI alone or in combination. Conclusion The BMMR2 challenge identified several models with high predictive performance, which may further expand the value of multiparametric breast MRI as an early marker of treatment response. Clinical trial registration no. NCT01042379 Keywords: MRI, Breast, Tumor Response Supplemental material is available for this article. © RSNA, 2024.

The ISMRM Open Science Initiative for Perfusion Imaging (OSIPI): Results from the OSIPI-Dynamic Contrast-Enhanced challenge.

PURPOSE: K trans $$ {K}^{\mathrm{trans}} $$ has often been proposed as a quantitative imaging biomarker for diagnosis, prognosis, and treatment response assessment for various tumors. None of the many software tools for K trans $$ {K}^{\mathrm{trans}} $$ quantification are standardized. The ISMRM Open Science Initiative for Perfusion Imaging-Dynamic Contrast-Enhanced (OSIPI-DCE) challenge was designed to benchmark methods to better help the efforts to standardize K trans $$ {K}^{\mathrm{trans}} $$ measurement. METHODS: A framework was created to evaluate K trans $$ {K}^{\mathrm{trans}} $$ values produced by DCE-MRI analysis pipelines to enable benchmarking. The perfusion MRI community was invited to apply their pipelines for K trans $$ {K}^{\mathrm{trans}} $$ quantification in glioblastoma from clinical and synthetic patients. Submissions were required to include the entrants' K trans $$ {K}^{\mathrm{trans}} $$ values, the applied software, and a standard operating procedure. These were evaluated using the proposed OSIP I gold $$ \mathrm{OSIP}{\mathrm{I}}_{\mathrm{gold}} $$ score defined with accuracy, repeatability, and reproducibility components. RESULTS: Across the 10 received submissions, the OSIP I gold $$ \mathrm{OSIP}{\mathrm{I}}_{\mathrm{gold}} $$ score ranged from 28% to 78% with a 59% median. The accuracy, repeatability, and reproducibility scores ranged from 0.54 to 0.92, 0.64 to 0.86, and 0.65 to 1.00, respectively (0-1 = lowest-highest). Manual arterial input function selection markedly affected the reproducibility and showed greater variability in K trans $$ {K}^{\mathrm{trans}} $$ analysis than automated methods. Furthermore, provision of a detailed standard operating procedure was critical for higher reproducibility. CONCLUSIONS: This study reports results from the OSIPI-DCE challenge and highlights the high inter-software variability within K trans $$ {K}^{\mathrm{trans}} $$ estimation, providing a framework for ongoing benchmarking against the scores presented. Through this challenge, the participating teams were ranked based on the performance of their software tools in the particular setting of this challenge. In a real-world clinical setting, many of these tools may perform differently with different benchmarking methodology.

An untrained deep learning method for reconstructing dynamic MR images from accelerated model-based data.

PURPOSE: To implement physics-based regularization as a stopping condition in tuning an untrained deep neural network for reconstructing MR images from accelerated data. METHODS: The ConvDecoder (CD) neural network was trained with a physics-based regularization term incorporating the spoiled gradient echo equation that describes variable-flip angle data. Fully-sampled variable-flip angle k-space data were retrospectively accelerated by factors of R = {8, 12, 18, 36} and reconstructed with CD, CD with the proposed regularization (CD + r), locally low-rank (LR) reconstruction, and compressed sensing with L1-wavelet regularization (L1). Final images from CD + r training were evaluated at the "argmin" of the regularization loss; whereas the CD, LR, and L1 reconstructions were chosen optimally based on ground truth data. The performance measures used were the normalized RMS error, the concordance correlation coefficient, and the structural similarity index. RESULTS: The CD + r reconstructions, chosen using the stopping condition, yielded structural similarity indexs that were similar to the CD (p = 0.47) and LR structural similarity indexs (p = 0.95) across R and that were significantly higher than the L1 structural similarity indexs (p = 0.04). The concordance correlation coefficient values for the CD + r T1 maps across all R and subjects were greater than those corresponding to the L1 (p = 0.15) and LR (p = 0.13) T1 maps, respectively. For R ≥ 12 (≤4.2 min scan time), L1 and LR T1 maps exhibit a loss of spatially refined details compared to CD + r. CONCLUSION: The use of an untrained neural network together with a physics-based regularization loss shows promise as a measure for determining the optimal stopping point in training without relying on fully-sampled ground truth data.

Analysis of simplicial complexes to determine when to sample for quantitative DCE MRI of the breast.

PURPOSE: A method is presented to select the optimal time points at which to measure DCE-MRI signal intensities, leaving time in the MR exam for high-spatial resolution image acquisition. THEORY: Simplicial complexes are generated from the Kety-Tofts model pharmacokinetic parameters Ktrans and ve . A geometric search selects optimal time points for accurate estimation of perfusion parameters. METHODS: The DCE-MRI data acquired in women with invasive breast cancer (N = 27) were used to retrospectively compare parameter maps fit to full and subsampled time courses. Simplicial complexes were generated for a fixed range of Kety-Tofts model parameters and for the parameter ranges weighted by estimates from the fully sampled data. The largest-area manifolds determined the optimal three time points for each case. Simulations were performed along with retrospectively subsampled data fits. The agreement was computed between the model parameters fit to three points and those fit to all points. RESULTS: The optimal three-point sample times were from the data-informed simplicial complex analysis and determined to be 65, 204, and 393 s after arrival of the contrast agent to breast tissue. In the patient data, tumor-median parameter values fit using all points and the three selected time points agreed with concordance correlation coefficients of 0.97 for Ktrans and 0.67 for ve . CONCLUSION: It is possible to accurately estimate pharmacokinetic parameters from three properly selected time points inserted into a clinical DCE-MRI breast exam. This technique can provide guidance on when to capture images for quantitative data between high-spatial-resolution DCE-MRI images.

Integrating mechanism-based modeling with biomedical imaging to build practical digital twins for clinical oncology.

Digital twins employ mathematical and computational models to virtually represent a physical object (e.g., planes and human organs), predict the behavior of the object, and enable decision-making to optimize the future behavior of the object. While digital twins have been widely used in engineering for decades, their applications to oncology are only just emerging. Due to advances in experimental techniques quantitatively characterizing cancer, as well as advances in the mathematical and computational sciences, the notion of building and applying digital twins to understand tumor dynamics and personalize the care of cancer patients has been increasingly appreciated. In this review, we present the opportunities and challenges of applying digital twins in clinical oncology, with a particular focus on integrating medical imaging with mechanism-based, tissue-scale mathematical modeling. Specifically, we first introduce the general digital twin framework and then illustrate existing applications of image-guided digital twins in healthcare. Next, we detail both the imaging and modeling techniques that provide practical opportunities to build patient-specific digital twins for oncology. We then describe the current challenges and limitations in developing image-guided, mechanism-based digital twins for oncology along with potential solutions. We conclude by outlining five fundamental questions that can serve as a roadmap when designing and building a practical digital twin for oncology and attempt to provide answers for a specific application to brain cancer. We hope that this contribution provides motivation for the imaging science, oncology, and computational communities to develop practical digital twin technologies to improve the care of patients battling cancer.

Quantitative multiparametric MRI predicts response to neoadjuvant therapy in the community setting.

BACKGROUND: The purpose of this study was to determine whether advanced quantitative magnetic resonance imaging (MRI) can be deployed outside of large, research-oriented academic hospitals and into community care settings to predict eventual pathological complete response (pCR) to neoadjuvant therapy (NAT) in patients with locally advanced breast cancer. METHODS: Patients with stage II/III breast cancer (N = 28) were enrolled in a multicenter study performed in community radiology settings. Dynamic contrast-enhanced (DCE) and diffusion-weighted (DW)-MRI data were acquired at four time points during the course of NAT. Estimates of the vascular perfusion and permeability, as assessed by the volume transfer rate (Ktrans) using the Patlak model, were generated from the DCE-MRI data while estimates of cell density, as assessed by the apparent diffusion coefficient (ADC), were calculated from DW-MRI data. Tumor volume was calculated using semi-automatic segmentation and combined with Ktrans and ADC to yield bulk tumor blood flow and cellularity, respectively. The percent change in quantitative parameters at each MRI scan was calculated and compared to pathological response at the time of surgery. The predictive accuracy of each MRI parameter at different time points was quantified using receiver operating characteristic curves. RESULTS: Tumor size and quantitative MRI parameters were similar at baseline between groups that achieved pCR (n = 8) and those that did not (n = 20). Patients achieving a pCR had a larger decline in volume and cellularity than those who did not achieve pCR after one cycle of NAT (p < 0.05). At the third and fourth MRI, changes in tumor volume, Ktrans, ADC, cellularity, and bulk tumor flow from baseline (pre-treatment) were all significantly greater (p < 0.05) in the cohort who achieved pCR compared to those patients with non-pCR. CONCLUSIONS: Quantitative analysis of DCE-MRI and DW-MRI can be implemented in the community care setting to accurately predict the response of breast cancer to NAT. Dissemination of quantitative MRI into the community setting allows for the incorporation of these parameters into the standard of care and increases the number of clinical community sites able to participate in novel drug trials that require quantitative MRI.

Quantitative magnetic resonance imaging and tumor forecasting of breast cancer patients in the community setting.

This protocol describes a complete data acquisition, analysis and computational forecasting pipeline for employing quantitative MRI data to predict the response of locally advanced breast cancer to neoadjuvant therapy in a community-based care setting. The methodology has previously been successfully applied to a heterogeneous patient population. The protocol details how to acquire the necessary images followed by registration, segmentation, quantitative perfusion and diffusion analysis, model calibration, and prediction. The data collection portion of the protocol requires ~25 min of scanning, postprocessing requires 2-3 h, and the model calibration and prediction components require ~10 h per patient depending on tumor size. The response of individual breast cancer patients to neoadjuvant therapy is forecast by application of a biophysical, reaction-diffusion mathematical model to these data. Successful application of the protocol results in coregistered MRI data from at least two scan visits that quantifies an individual tumor's size, cellularity and vascular properties. This enables a spatially resolved prediction of how a particular patient's tumor will respond to therapy. Expertise in image acquisition and analysis, as well as the numerical solution of partial differential equations, is required to carry out this protocol.

Characterizing Errors in Pharmacokinetic Parameters from Analyzing Quantitative Abbreviated DCE-MRI Data in Breast Cancer.

This study characterizes the error that results when performing quantitative analysis of abbreviated dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) data of the breast with the Standard Kety-Tofts (SKT) model and its Patlak variant. More specifically, we used simulations and patient data to determine the accuracy with which abbreviated time course data could reproduce the pharmacokinetic parameters, Ktrans (volume transfer constant) and ve (extravascular/extracellular volume fraction), when compared to the full time course data. SKT analysis of simulated abbreviated time courses (ATCs) based on the imaging parameters from two available datasets (collected with a 3T MRI scanner) at a temporal resolution of 15 s (N = 15) and 7.23 s (N = 15) found a concordance correlation coefficient (CCC) greater than 0.80 for ATCs of length 3.0 and 2.5 min, respectively, for the Ktrans parameter. Analysis of the experimental data found that at least 90% of patients met this CCC cut-off of 0.80 for the ATCs of the aforementioned lengths. Patlak analysis of experimental data found that 80% of patients from the 15 s resolution dataset and 90% of patients from the 7.27 s resolution dataset met the 0.80 CCC cut-off for ATC lengths of 1.25 and 1.09 min, respectively. This study provides evidence for both the feasibility and potential utility of performing a quantitative analysis of abbreviated breast DCE-MRI in conjunction with acquisition of current standard-of-care high resolution scans without significant loss of information in the community setting.

The rate of breast fibroglandular enhancement during dynamic contrast-enhanced MRI reflects response to neoadjuvant therapy.

PURPOSE: This study assesses the rate of enhancement of breast fibroglandular tissue after administration of a magnetic resonance imaging (MRI) gadolinium-based contrast agent and determines its relationship with response to neoadjuvant therapy (NAT) in women with breast cancer. METHOD: Women with locally advanced breast cancer (N = 19) were imaged four times over the course of NAT. Dynamic contrast-enhanced (DCE) MRI was acquired after administration of a gadolinium-based contrast agent with a temporal resolution of 7.27 s. The tumor, fibroglandular tissue, and adipose tissue were semi-automatically segmented using a manually drawn region of interest encompassing the tumor followed by fuzzy c-means clustering. The rate and relative intensity of signal enhancement were calculated for each voxel within the tumor and fibroglandular tissue. RESULTS: The rate of fibroglandular tissue enhancement after contrast agent injection declined by an average of 29 % over the course of NAT. This decline was present in 16 of the 19 patients in the study. The rate of enhancement is significantly higher in women who achieve pathological complete response (pCR) after both 1 cycle (68 % higher, p < 0.05) and after 3-5 cycles of NAT (58 % higher; p < 0.05). The relative intensity of fibroglandular enhancement correlates with the rate of enhancement (R2 = 0.64, p < 0.001) and is higher in women who achieve pCR after both 1 cycle and after 3-5 cycles of NAT (p < 0.05, both timepoints). CONCLUSION: The rate of fibroglandular tissue enhancement declines over the course of therapy, provides novel information not reflected by tumoral measures, and may predict pathological response early in the course of therapy, with smaller declines in enhancement in women who achieve favorable response.

Evaluating patient-specific neoadjuvant regimens for breast cancer via a mathematical model constrained by quantitative magnetic resonance imaging data.

The ability to accurately predict response and then rigorously optimize a therapeutic regimen on a patient-specific basis, would transform oncology. Toward this end, we have developed an experimental-mathematical framework that integrates quantitative magnetic resonance imaging (MRI) data into a biophysical model to predict patient-specific treatment response of locally advanced breast cancer to neoadjuvant therapy. Diffusion-weighted and dynamic contrast-enhanced MRI data is collected prior to therapy, after 1 cycle of therapy, and at the completion of the first therapeutic regimen. The model is initialized and calibrated with the first 2 patient-specific MRI data sets to predict response at the third, which is then compared to patient outcomes (N = 18). The model's predictions for total cellularity, total volume, and the longest axis at the completion of the regimen are significant within expected measurement precision (P< 0.05) and strongly correlated with measured response (P < 0.01). Further, we use the model to investigate, in silico, a range of (practical) alternative treatment plans to achieve the greatest possible tumor control for each individual in a subgroup of patients (N = 13). The model identifies alternative dosing strategies predicted to achieve greater tumor control compared to the standard of care for 12 of 13 patients (P < 0.01). In summary, a predictive, mechanism-based mathematical model has demonstrated the ability to identify alternative treatment regimens that are forecasted to outperform the therapeutic regimens the patients clinically. This has important implications for clinical trial design with the opportunity to alter oncology care in the future.

Sustained improvements in MRI outcomes with abatacept following the withdrawal of all treatments in patients with early, progressive rheumatoid arthritis.

OBJECTIVES: To assess structural damage progression with subcutaneous abatacept (ABA) in the Assessing Very Early Rheumatoid arthritis Treatment (AVERT) trial following abrupt withdrawal of all rheumatoid arthritis (RA) medication in patients achieving Disease Activity Score (DAS)-defined remission or low disease activity. METHODS: Patients with early, active RA were randomised to ABA plus methotrexate (ABA/MTX) 125 mg/week, ABA 125 mg/week or MTX for 12 months. All RA treatments were withdrawn after 12 months in patients with DAS28 (C reactive protein (CRP)) <3.2. Adjusted mean changes from baseline in MRI-based synovitis, osteitis and erosion were calculated for the intention-to-treat population. RESULTS: 351 patients were randomised and treated: ABA/MTX (n=119), ABA (n=116) or MTX (n=116). Synovitis and osteitis improved, and progression of erosion was statistically less with ABA/MTX versus MTX at month 12 (-2.35 vs -0.68, -2.58 vs -0.68, 0.19 vs 1.53, respectively; p<0.01 for each) and month 18 (-1.34 vs -0.49 -2.03 vs 0.34, 0.13 vs 2.0, respectively; p<0.01 for erosion); ABA benefits were numerically intermediate to those for ABA/MTX and MTX. CONCLUSIONS: Structural benefits with ABA/MTX or ABA may be maintained 6 months after withdrawal of all treatments in patients who have achieved remission or low disease activity. TRIAL REGISTRATION NUMBER: NCT01142726; Results.

Monitoring cartilage loss in the hands and wrists in rheumatoid arthritis with magnetic resonance imaging in a multi-center clinical trial: IMPRESS (NCT00425932).

INTRODUCTION: Magnetic resonance imaging (MRI) is increasingly being used in clinical trials of rheumatoid arthritis (RA) because of its superiority over x-ray radiography (XR) in detecting and monitoring change in bone erosion, osteitis and synovitis. However, in contrast to XR, the MRI scoring method that was used in most clinical trials did not include cartilage loss. This limitation has been an obstacle to accepting MRI as a potential alternative to XR in clinical trials. Cross-sectional studies have shown MRI to be sensitive for cartilage loss in the hands and wrist; although, longitudinal sensitivity to change has not yet been confirmed. In this study we examined the ability of MRI to monitor change in cartilage loss in patients with RA in a multi-site clinical trial setting. METHODS: Thirty-one active RA patients from a clinical trial (IMPRESS) who were randomized equally into treatment with either rituximab + methotrexate or placebo + methotrexate had MRI of the dominant hand/wrist at baseline, 12 weeks and 24 weeks at 3 clinical sites in the US. Twenty-seven of these patients also had XR of both hands/wrists and both feet at baseline and 24 weeks. One radiologist scored all XR images using the van der Heijde-modified Sharp method blinded to visit order. The same radiologist scored MR images for cartilage loss using a previously validated 9-point scale, and bone erosion using the Outcome Measures in Rheumatology Clinical Trials (OMERACT) RA MRI Score (RAMRIS) blinded to visit order and XR scores. Data from the two treatment arms were pooled for this analysis. RESULTS: Mean MRI cartilage score increased at 12 and 24 weeks, and reached statistical significance at 24 weeks. XR total Sharp score, XR erosion score and XR joint-space narrowing (JSN) score all increased at 24 weeks, but only XR total Sharp score increased significantly. CONCLUSIONS: To our knowledge, this is the first publication of a study demonstrating MRI's ability to monitor cartilage loss in a multi-site clinical trial. Combined with MRI's established performance in monitoring bone erosions in RA, these findings suggest that MRI may offer a superior alternative to XR in multi-site clinical trials of RA.

Evaluating joint-space narrowing and cartilage loss in rheumatoid arthritis by using MRI.

INTRODUCTION: Magnetic resonance imaging (MRI) has been shown to be superior to radiography (XR) for assessing synovitis, osteitis, and bone erosion in rheumatoid arthritis (RA), particularly in clinical trials. However, relatively little has been reported on the ability of MRI to evaluate articular cartilage loss, or joint-space narrowing (JSN), in the hands and wrists. In a previous study, we adapted the nine-point Genant-modified Sharp XR-JSN score for use with MRI (MRI-JSN). In this study, we compare MRI-JSN with XR-JSN by using images from two multicenter clinical trials. METHODS: Baseline XR and 1.5-Tesla MR images of one hand and wrist from each of 47 subjects with RA enrolled in one of two multicenter clinical trials were evaluated by using the XR-JSN and MRI-JSN methods by a single radiologist experienced in the two methods. Radiographs and MR images were read independently on different occasions. RESULTS: In total, 575 of 611 joints were compared (one metacarpophalangeal joint of the thumb and 35 proximal interphalangeal joints were outside the MRI field of view and could not be assessed). The 22 (47%) subjects showed JSN with both XR and MRI, and 25 (53%) subjects showed no JSN with either method. No subject showed JSN with only one or the other method. MRI showed high agreement with XR (intraclass correlation coefficient = 0.83). Sensitivity of MRI for JSN, by using XR as the gold standard, was 0.94; specificity was 0.91; accuracy was 0.91; positive predictive value was 0.64; and negative predictive value was 0.99. CONCLUSIONS: This validation exercise suggests that MRI JSN scoring may offer a viable alternative to XR JSN scoring in multicenter clinical trials of RA. However, the relative longitudinal sensitivity of MRI to change and the ability to discriminate therapeutic effect on JSN were not evaluated in this study.

Variable-density one-shot Fourier velocity encoding.

In areas of highly pulsatile and turbulent flow, real-time imaging with high temporal, spatial, and velocity resolution is essential. The use of 1D Fourier velocity encoding (FVE) was previously demonstrated for velocity measurement in real time, with fewer effects resulting from off-resonance. The application of variable-density sampling is proposed to improve velocity measurement without a significant increase in readout time or the addition of aliasing artifacts. Two sequence comparisons are presented to improve velocity resolution or increase the velocity field of view (FOV) to unambiguously measure velocities up to 5 m/s without aliasing. The results from a tube flow phantom, a stenosis phantom, and healthy volunteers are presented, along with a comparison of measurements using Doppler ultrasound (US). The studies confirm that variable-density acquisition of kz-kv space improves the velocity resolution and FOV of such data, with the greatest impact on the improvement of FOV to include velocities in stenotic ranges.