Abstract

The aim of this study was to assess whether sequential multiparametric 18[F]fluoro-desoxy-glucose (18[F]FDG)/[18F]fluoromisonidazole ([18F]FMISO) PET-MRI in breast cancer patients is possible, facilitates information on tumor heterogeneity, and correlates with prognostic indicators. In this pilot study, IRB-approved, prospective study, nine patients with ten suspicious breast lesions (BIRADS 5) and subsequent breast cancer diagnosis underwent sequential combined [18F]FDG/[18F]FMISO PET-MRI. [18F]FDG was used to assess increased glycolysis, while [18F]FMISO was used to detect tumor hypoxia. MRI protocol included dynamic breast contrast-enhanced MRI (DCE-MRI) and diffusion-weighted imaging (DWI). Qualitative and quantitative multiparametric imaging findings were compared with pathological features (grading, proliferation, and receptor status) and clinical endpoints (recurrence/metastases and disease-specific death) using multiple correlation analysis. Histopathology was the standard of reference. There were several intermediate to strong correlations identified between quantitative bioimaging markers, histopathologic tumor characteristics, and clinical endpoints. Based on correlation analysis, multiparametric criteria provided independent information. The prognostic indicators proliferation rate, death, and presence/development of recurrence/metastasis correlated positively, whereas the prognostic indicator estrogen receptor status correlated negatively with PET parameters. The strongest correlations were found between disease-specific death and [18F]FDGmean (, ) and between the presence/development of metastasis and [18F]FDGmax (, ), respectively. This pilot study indicates that multiparametric [18F]FDG/[18F]FMISO PET-MRI might provide complementary quantitative prognostic information on breast tumors including clinical endpoints and thus might be used to tailor treatment for precision medicine in breast cancer.

1. Introduction

Breast cancer is a complex disease with remarkable intratumoral heterogeneity resulting in different clinical and phenotypic presentations, treatment responses, and outcomes [14]. Magnetic resonance imaging (MRI) and positron emission tomography (PET) play an important role in staging, therapy monitoring, and follow-up of breast cancer [5]. Recently, hybrid PET/MRI scanners have become clinically available. An inherent advantage of PET-MRI is that it provides noninvasive qualitative and quantitative in vivo information as well as spatiolongitudinal monitoring of the tumor microenvironment and its interaction, while invasive tissue sampling provides just “snapshots” of specific tumor regions.

The routinely used radiotracer [18F]FDG depicts increased tissue glycolysis indicative of malignancy and has high sensitivity for the detection of breast cancer. Initial studies have already demonstrated the potential of 2-deoxy-2-[18F]fluoro-D-glucose ([18F]FDG) PET/MRI for breast cancer diagnosis, characterization, and prognosis [6, 7]. However, for PET imaging, several other radiotracers that specifically target cancer hallmark processes relevant for cancer development, progression, and treatment resistance have also been developed, such as [18F]fluoromisonidazole ([18F]FMISO) for imaging hypoxia [8], [18F]fluorothymidine ([18F]FLT) for imaging proliferation [9], and [18F]fluoroestradiol ([18F]FES) for imaging hormone receptor status. The combined use of different radiotracers may add diagnostic, predictive, and prognostic information, and therefore, the full potential of PET-MRI in breast cancer is yet to be realized.

In breast cancer, tumor hypoxia has been recognized as an important feature and a key driver of intratumoral heterogeneity which in turn leads to the development of cell clones with an aggressive and treatment-resistant phenotype characterized by rapid progression and a poor prognosis [1012]. Hypoxia has been implicated as a confounder of the efficacy of cancer therapies as well as a prognostic factor for disease progression, metastases, and survival [1012]. PET imaging using [18F]FMISO has been shown to identify hypoxic tumor subvolumes and track spatiotemporal dynamics. Therefore, it might provide additional information to [18F]FDG particularly on breast cancer heterogeneity and have additional prognostic value. Therefore, the aim of this pilot study was twofold:(i)To assess whether combined sequential, dual tracer multiparametric PET-MRI at 3T in primary breast cancer patients is possible and facilitates information on tumor heterogeneity(ii)To investigate associations with recurrence and disease-specific patient survival

To reach this goal, we used multiple MRI parameters (high-resolution T2-weighted, dynamic contrast-enhanced (DCE)-MRI, and diffusion-weighted imaging (DWI) with apparent diffusion coefficient (ADC) mapping) and combined these parameters with that from PET using the radiotracers [18F]FDG for the assessment of glycolytic metabolic activity and [18F]FMISO for the detection of tumor hypoxia.

2. Materials and Methods

The Institutional Review Board approved this prospective, single-institution study, and all patients gave written, informed consent.

2.1. Patients

Nine patients (mean age of 53.2 years) with 10 breast cancer lesions who fulfilled the following inclusion criteria were included: 18 years or older; not pregnant; not breastfeeding; imaging finding at mammography or breast ultrasonography highly suggestive of malignancy (BI-RADS 5); no previous treatment; and no contraindications for MRI or MRI contrast agents (1). For all patients, the following information was recorded: age, histologic type, tumor grade, receptor status, tumor proliferation rate (Ki67), nodal status, presence of regional or distant metastases, date of progression (local recurrence, distant metastases) to determine duration (months) of recurrence-free survival (RFS), and date and cause of death or date of last follow-up to determine duration (months) of disease-specific death (DSS).

2.2. Imaging

All patients underwent sequential combined multiparametric dual tracer [18F]FDG/[18F]FMISO PET-MRI at 3T with examinations performed no more than 7 days apart.

2.3. MR Imaging

The MRI examinations were performed with the patient in the prone position using a 3T MRI (Tim Trio, Siemens, Erlangen, Germany) and a four-channel breast coil (InVivo, Orlando, FL, USA). The MRI protocol included the following sequences:(1)T2-weighted with fat suppression: axial turbo inversion recovery magnitude (TIRM) sequence: repetition time/echo time (TR/TE) 4800/59 msec; field of view (FOV) 340 mm; 44 slices at SI 4 mm; flip angle (FA) 120°; matrix 384 × 512; acquisition time (AT): 2 min 35 sec(2)T2-weighted with fat suppression: turbo spin echo sequence: TR/TE = 4630/194 msec; FOV 340 mm; 65 slices at slice thickness (SI) 2.5 mm; FA 128°; matrix 384 × 640; AT: 2 min 48 sec(3)DWI: axial three-acquisition trace diffusion-weighted, double-refocused, single-shot echo-planar imaging sequence with inversion recovery fat suppression: TR/TE/time of inversion (TI) 8000/59/210 msec; FOV 360 × 202 mm; 24 slices at 5 mm; matrix 172 × 96 (50% oversampling); b-values 50 and 850 s/mm, AT: 2 min 56 sec(4)DCE-MRI: transversal T1-weighted time-resolved angiography with stochastic trajectories (TWIST) with water excitation fat saturation: TR/TE 6.23/2.95 ms; FA 15°, FOV 196 × 330 mm 144 slices; spatial resolution 0.9 × 0.9 × 1 mm; temporal interpolation factor 2; temporal resolution 14 s; matrix 384 × 384; one average; center k-space region with a resampling rate of 23%; reacquisition density of peripheral k-space 20%; AT: 6 min 49 sec

A standard dose (0.1 mmol/kg of bodyweight) of Gd-DOTA Dotarem was administered intravenously as a bolus with a power injector (Spectris Solaris EP; Medrad) at 4 mL/s followed by a saline flush.

2.4. PET Imaging

All PET studies were acquired using a hybrid PET/CT (computed tomography) (Biograph 64 TruePoint PET/CT system, Siemens, Erlangen/Germany). For [18F]FDG PET/CT, patients fasted for five hours. Patients were injected with 3 MBq [18F]FDG and [18F]FMISO per kilogram bodyweight on different days. Scanning was started after an uptake time of 45 min for [18F]FDG and 210–240 min for [18F]FMISO. For both radiotracers, a prone PET dataset with low-dose unenhanced CT scans was recorded for attenuation correction. The same imaging and postreconstruction parameters were used for both PET/CT studies. Images were reconstructed using the iterative TrueX algorithm, which incorporates a specific correction for the point-spread function in addition to commonly used correction factors [13, 14]. Four iterations per 21 subsets were used, with a matrix size of 168 × 168, a transaxial FOV of 605 mm (pixel size 3.6 mm), and a section thickness of 5 mm. Further technical details are provided by the manufacturer [15]. For both examinations, patient positioning and the scan time took on average ∼40 minutes.

2.5. Image Fusion

All images were coregistered and fused to allow delineation of regions of interest (ROIs). An automatic, region-based, rigid, fine image registration protocol was used in the Mirada RTx software (Mirada Medical Ltd.) and followed with manual adjustments if necessary (5–10 min per pair of images). The MRI-PET registration was performed in two steps: first, the T2-weighted MRI and the corresponding CT were registered, and second, that transformation was applied to the PET dataset (Figure S1). Within the MRI series, the initial frame of reference was used as the starting point for rigid registration, correcting for possible patient movement. At the end of the process, all images were aligned with the T2-weighted MRI (Figure 1). For the purpose of voxel-by-voxel analysis, all datasets were also resampled to the resolution of the PET: 4.1 × 4.1 × 3.0 mm.

2.6. Data Analysis

For each patient, the tumor was delineated by an experienced breast radiologist and nuclear medicine physician, and the following parameters were derived from the delineated regions of interest (ROI): volume (cm3); mean signal intensity on T2-weighted images with and without fat saturation as well as on native, early, and delayed DCE-MR images (DCEnative, DCEearly, and DCEdelayed, respectively); initial enhancement (IE) ratio and washout (WO) ratio (IE ratio defined as (DCEearly−DCEnative)/DCEnative and WO ratio defined as (DCEdelayed−DCEearly)/DCEnative)); mean ADC value on DWI; mean and maximum standard uptake value (SUV) for each tracer; and tumor to background (TBR) ratio (calculated as the ROI SUV normalized to the SUV measured in the patients’ aorta) for each tracer. For all PET-derived parameters, a ratio between [18F]FDG and respective [18F]FMISO parameter was calculated.

2.7. Statistical Analysis

Imaging findings and calculated ratios were compared with clinical data (i.e., histopathologic tumor characteristics including hormone receptor status (estrogen, progesterone and human epidermal growth factor receptor 2 receptor) and proliferation rate (Ki67)) and patient outcomes (i.e., disease-specific death (DSS) and presence/development of metastases) using multiple correlation analysis.

Person correlation coefficients (R) were also calculated for each parameter pair. Statistically significant correlations were color-coded (blue for negative correlations and green for positive correlations). values were marked in green for significance at the ≦0.05 level (2-tailed) and red for that at the ≦0.01 level. For image-derived parameters, Person correlation coefficients were additionally calculated on a voxel-by-voxel level. For PET-derived TBR, a value in each voxel was considered in relation to the SUVaorta; hence, TBRmean and TBRmax were not used [1618]. Voxel-by-voxel analysis was additionally performed for [18F]FMISOTBR > 1.4 and [18F]FDGTBR > 2.0 parameters, where only the voxels in which TBR was larger than 1.4 or 2.0, respectively, were considered, and only moderate and strong correlations were color-coded (Figure S2).

3. Results

Relevant clinical data, pathological tumor characteristics, and clinical endpoints for each patient are summarized in Table 1. The mean follow-up time for all patients was 31.6 months (range, 3–60 months). One patient presented with one breast cancer in each breast which were of different tumor biology; we therefore considered this patient as two separate cases (ID 6: IDC ER/PR negative, Her2 positive, and Ki67 90%; ID 7: IDC ER/PR positive, Her2 negative, and Ki67 20%). There were six (60%) cases of distant metastases present at diagnosis, two (20%) cases of distant metastases present at six and 17 months, respectively, and two (20%) cases of no recurrence. In five cases (50%), the patient died of breast cancer during the follow-up period at a median interval of 10 months (range, 3–19 months).

Several intermediate to strong correlations were identified between quantitative imaging markers, histopathologic tumor characteristics, and outcome measures as summarized in Figure 1. Signal intensity on T2-weighted images and ADCmean values correlated negatively with [18F]FMISO and [18F]FDG parameters; the strongest correlations were found between signal intensity on T2-weighted images and [18F]FMISOmean (, ) and between ADCmean and [18F]FDGTBRmax (, ). Among the quantitative MR imaging parameters, the strongest correlation was found between signal intensity on T2-weighted fat-saturated images and IE ratio (, ). Among the PET parameters, a noticeable trend of moderate correlation between the PET markers was observed in the middle cluster of the graph in Figure 2. A moderate-to-strong positive correlation was found between [18F]FMISOmax and lesion volume (, ) and between [18F]FDGTBRmax and [18F]FMISOTBRmax (, ).

Both [18F]FMISOTBRmean (, ) and [18F]FDGmean (, ) correlated strongly with Ki67. There was a negative correlation of PET parameters with ER; the strongest correlation was found between [18F]FDGmean and ER (, ) (Table 2).

Clinical endpoints such as death and presence/development of metastasis correlated positively with PET parameters. The strongest correlation was found between DSS and [18F]FDGmean (, ) and between metastases and [18F]FDGmax (, ). There were also moderate correlations for DSS with both [18F]FMISOTBRmean and [18F]FMISOTBRmax. The calculated ratio of [18F]FDGTBRmax/[18F]FMISOTBRmax had a strong correlation with the presence/development of metastases (, ).

In voxel-by-voxel analysis, no significant correlations were identified between MRI parameters or MRI and PET parameters. Only a weak correlation between two PET tracers without restricted TBR was found (Figure S2).

4. Discussion

Results of high-throughput molecular profiling studies have revealed that breast cancer is a disease with a remarkable tumor heterogeneity resulting in varying genetic, phenotypic, and behavioral characteristics; clinical presentations; and treatment responses [1, 49]. This recognized tumor heterogeneity, and the lack of understanding thereof significantly contributes to treatment failures and breast cancer deaths. Breast cancer heterogeneity is driven by tumor genomic instability and selective pressures from the tumor microenvironment with hypoxia being one of the most significant factors. Tumor hypoxia leads to the development of cell clones with an aggressive and treatment-resistant phenotype characterized by rapid progression and a poor prognosis [1012]. Hypoxia has been implicated as a confounder of the efficacy of cancer therapies as well as a prognostic factor for disease progression, metastases, and survival [1012]. To date, current invasive tools and conventional imaging technologies cannot provide a comprehensive assessment of breast cancer heterogeneity of the whole tumor. In this context, there is a unique opportunity for noninvasive functional imaging with PET-MRI in breast cancer using multiple MRI parameters, yet to date its potential has not been explored in detail. To our knowledge, this is the first study that investigates whether combined sequential, dual tracer multiparametric PET-MRI at 3T in primary breast cancer patients is possible, facilitates information on tumor heterogeneity, and correlates with survival outcomes. In this pilot study performed in a limited number of patients, we show that a comprehensive noninvasive assessment of breast cancer biology, heterogeneity, hypoxic tumor microenvironment, and the degree to which breast cancer is successful in adapting to sustain growth is feasible with in vivo functional imaging.

We demonstrated the feasibility and potential of multiparametric [18F]FDG/[18F]FMISO PET-MRI in breast cancer patients. Our pilot study showed several intermediate to strong correlations between quantitative imaging biomarkers, histopathologic tumor characteristics, and clinical endpoints. Based on correlation analysis, multiparametric criteria provided independent information on tumor heterogeneity that is known to significantly contribute to treatment failures and breast cancer deaths; in the future, this information can be used to personalize cancer treatment.

Multiparametric MRI of the breast is already the state-of-the-art for the staging of breast cancer patients, and hybrid PET/MR scanners are now clinically available; therefore, this information will readily be available. In this pilot study, there was no correlation of MRI parameters and clinical endpoints, which might be due to the small sample size as prior studies have already demonstrated both a prognostic and predictive potential [1926]. For instance, Durando et al. [25] showed that ADC can be a prognostic factor for malignancy aggressiveness; Dietzel et al. [20] and Alduk et al. [19] showed that MRI-derived parameters associate with nodal status, predict positive axillary lymph nodes, or larger tumor sizes; and Bae et al. [23] found significant association between pretreatment T2 MRI features with recurrence-free survival. On the contrary, [18F]FDG and [18F]FMISO PET parameters showed a moderate-to-strong positive correlation with tumor proliferation rate (Ki67) (a poor prognostic indicator [27]) and DSS. [18F]FDG demonstrated strong positive correlation with the presence/development of metastasis as well as a strong negative correlation with ER positivity (a good prognostic indicator [28]). The results of the current study support findings by others that have demonstrated that [18F]FDG and MRI parameters may aid in the assessment of tumor aggressiveness and metastatic potential, which further highlights the value of noninvasive functional imaging for improved risk assessment [7].

Although, in this pilot study, patients with high [18F]FMISO uptake did not show the tendency to present with or develop metastasis, there was a positive correlation yet weaker correlation than [18F]FDGmean of the [18F]FDGTBRmax/[18F]FMISOTBRmax ratio with presence of metastasis. These interesting findings suggest that whenever within a metabolically active tumor there is significant hypoxia, which is a key driver for the emergence of treatment resistant and aggressive clones, these patients are at an increased risk of metastases and might benefit from intensified treatment or tighter follow-up. In addition, there was a correlation of increased [18F]FMISO uptake with death, which further suggests that these tumors are more treatment-resistant. This is in agreement with findings from Asano et al. who demonstrated that higher [18F] FMISOTBR and ER negativity were independent predictors of shorter disease-free survival, where higher [18F] FMISOTBR was associated with higher plasma levels of angiogenic hypoxic markers [29]. Cheng et al. also showed that [18F]FMISO PET/CT can be used to predict primary endocrine resistance in ER-positive breast cancer [8].

Voxel-by-voxel analysis showed only low correlations between two PET tracers, whereas descriptive statistics suggested higher correlations. This implies that tumors characterized by regions of increased metabolism have at the same time highly hypoxic subregions that are spatially disjointed. It has to be noted that the voxel-wise analysis is per se prone to spatial inaccuracies (i.e., dependent on image registration and temporal difference between the PET scans). Despite the relatively low resolution of the PET scans (resulting from the physics of the annihilation process) combined with the resampling of the image dataset, the voxel-by-voxel analysis error might be locally exceeding the dimension of one voxel. Therefore, our findings will have to be further validated, preferentially using a simultaneous hybrid PET/MRI.

The main limitation of this pilot study, by virtue of the novelty of the imaging protocol, is the small patient cohort, limiting statistical analysis, and therefore, no strong conclusions can be drawn. In addition, due to this fact, we did not perform modeling to assess which parameters from [18F]FDG/[18F]FMISO PET-MRI are most strongly associated with RFS and DSS. Further larger studies, preferably with simultaneous hybrid PET/MRI scanners, are necessary to confirm these findings and to clarify which parameters independently or jointly provide prognostic information on breast tumors, including clinical endpoints. In this study, we performed both qualitative and quantitative analysis of [18F]FDG/[18F]FMISO PET-MRI data. The T2w MRI signal intensity was not normalized, and as it can vary from scan to scan depending on the gain settings, it cannot be considered quantitative.

5. Conclusions

In conclusion, according to these preliminary results, multiparametric [18F]FDG/[18F]FMISO PET-MRI provides complementary quantitative prognostic information on breast tumors including clinical endpoints and thus might be used to tailor treatment for precision medicine in breast cancer.

Data Availability

The imaging data used to support the findings of this study are available from the corresponding author upon request.

Disclosure

THH, DG, MH, and WW are members of the Medical Imaging Cluster (MIC) of the Medical University Vienna.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

Funding was provided by the Austrian Nationalbank “Jubiläumsfond” (Project Nos. 16219 and 15082), the 2020–Research and Innovation Framework Programme (PHC-11-2015 No. 667211-2), and seed grants from Novomed, Austria, and Guerbet, France. The financial support by the Federal Ministry of Science, Research and Economy and the National Foundation for Research, Technology and Development is gratefully acknowledged. Katja Pinker was also funded in part through the NIH/NCI Cancer Center Support Grant P30 CA008748. The authors acknowledge the support in manuscript writing and editing from Joanne Chin.

Supplementary Materials

Figure S1: example of the registration of the CTFDG, CTFMISO, and T2w MRI: (A) CTFDG; (B) T2w MRI; (C) fusion of registered A and B; (D) CTFDG; (E) CTFMISO; (F) fusion of registered D and E. Figure S2: graph showing numeric and color-coded Person’s correlation coefficients of all imaging parameters performed voxel-by-voxel within each patient’s delineated ROIs. Abbreviations: T2w, T2-weighted magnetic resonance imaging (MRI); TIRM, turbo inversion recovery magnitude MRI; ADC, apparent diffusion coefficient; IE, initial enhancement ratio; WO washout ratio; FMISO/FDGTBR, tumor to background ratio based on standard uptake value (SUV) in a certain voxel normalized to mean SUV in the aorta; FMISOTBR > 1.4 and FDGTBR > 2.0, as defined previously but only taking into account voxels with TBR > 1.4 and 2.0, respectively. (Supplementary Materials)