Patients
Consecutive patients that were suspected to have lung cancer who received both CT and 18F-FDG PET/MRI before surgery from November 2015 to May 2019 at National Cancer Center Hospital were enrolled in this study. The inclusion criteria were as follows: 1) pathologically proven non-small cell lung cancer, 2) both CT and PET/MRI were performed before thoracic surgery, and 3) follow-up after the operation was performed at the same institute. The exclusion criteria were as follows: 1) preoperative distant metastatic disease, 2) advanced primary cancer other than NSCLC, 3) treatment with chemotherapy before surgery, 4) multiple primary tumors, 5) invasive mucinous adenocarcinoma, and 6) small nodular lesion less than 5 mm on the long axis by thin-section CT images and/or PET/MRI. Since this was a retrospective single-center study, the informed consent requirement was waived by the institutional review board (research proposal number 2018-048).
18F-FDG PET/MRI acquisition and data analysis
Before the injection of 18F-FDG, patients fasted for at least 4 h. If plasma glucose levels were less than 200 mg/dL, patients were injected intravenously with 3-4 MBq/kg of the radiotracer, depending on patient height and body weight. After injection, patients rested for 50–70 min before image acquisition. 18F-FDG PET/MRI imaging was performed using a 3T PET/MRI (SIGNA PET/MRI, GE Healthcare, Milwaukee, WI, USA). Patients were positioned in the supine headfirst position. The acquisition was started in the pelvic region and moved toward the head. A localizer MRI scan was performed to define the bed positions. The axial field of view contained the body volume from the head to the thigh.
The adaptive PET/MRI protocol for the lung after whole-body PET/MRI was performed using a dedicated coil. The standard pulse sequences consisted of multiplanar 2D T2-weighted fast spin-echo (transaxial and coronal) with a 5-mm slice thickness. To avoid respiratory motion-induced artifacts, all lung surveys were performed with breath-holding. Furthermore, the protocol included a transaxial 3D T1-weighted spoiled gradient-echo sequence, with and without contrast enhancement, with 5-mm slice thickness (axial and coronal) covering the primary site of the lung cancer. Intravenous contrast agents (meglumine gadoterate; Guerbet, Tokyo, Japan) were administered.
Both the whole-body MRI and the adaptive PET/MRI images were reconstructed according to the ordered-subsets expectation maximization using 2 iterations, 16 subsets, with a matrix size of 192 × 192, and 5.0 mm3 post-filter at a dedicated workstation with MRI-based attenuation correction, scatter, and decay corrections. The voxel size for PET was 3.125 × 3.125 × 2.780 mm3 and that of MR images was 2.0 × 1.6 × 2.0 mm3 for whole-body imaging and 0.58 × 0.58 × 2.5 mm3 for adaptive chest imaging. After correct positioning of the spatial acquisition windows had been ensured, the whole-body PET/MRI acquisition was initiated at a 2-min acquisition time per bed position. For attenuation correction of the PET data from the PET/MRI scanner, the reconstruction software provided by the manufacturer used attenuation maps generated based on the 2-point Dixon MRI sequences obtained for every bed position. This approach has recently been demonstrated to provide results comparable to those of conventional attenuation correction by low dose CT [13, 14]. The procedure was implemented in the post-processing software of the scanner and operated automatically. The Dixon fat- and water-weighted images were used to create an attenuation map with 4 distinct tissue-classes: background, lungs, fat, and soft tissue. The lungs were identified by connected component analysis of the air in the inner part of the body. By application of a morphologic closing filter, virtual air artifacts induced by the absence of an MRI signal in cortical bone, heart, and aorta (due to blood flow) were corrected. Attenuation of the PET signal caused by instrumentation such as the patient bed and the fixed MRI coils was automatically integrated into the attenuation maps [15] without contrast enhancement, at 5-mm slice thickness (axial and coronal).
CT imaging
All CT scans of the study were performed on a multi-detector CT scanner (Aquilion PRIME, Aquilion Precision, Aquilion 64; Canon Medical Systems Corporation). Scans were performed supine after full inspiration with caudocranial scan direction including the entire ribcage and upper abdomen with Tube-current modulation (range, 85-545 mA), and 120 kV. An intravenous contrast agent (iopamidol; FujiPharma, Toyama, Japan) was used in most patients. Contiguous 5-mm section images were obtained and reconstruction of 1-mm thin-slice section images of the lung around the tumor was added.
All CT images were viewed using standard lung windows (window level, -600 HU; window width, 1900 HU). All image data were stored in DICOM format on PACS.
Image analysis
TNM staging was determined according to the UICC TNM Classification of Malignant Tumors, Eighth Edition. Location of primary tumors, tumor size in the greatest dimension, density of tumors, atelectasis, separate tumor nodules, and tumor invasion into the other structures were registered for clinical T-factor. Considering the quality of PET/MRI images, the size criteria of T1 and T2 in the TNM classification were only defined as ≤3 cm and 3 cm to ≤5 cm, respectively, and no subdivisions were used. Density criteria were assessed by only thin-section CT alone and not by PET/MRI. Location of primary tumors, atelectasis, separate tumor nodules, and tumor invasion into other structures were assessed based on visual assessments of each modality.
One board-certified radiologist (HW) specialized in chest imaging with 25 years of experience reviewed all chest CT images. Tumor size was assessed based on linear measurement of the largest diameter in axial thin-section CT-lung window. The reader used the monitor (Eizo RadiForce RX440) to measure the tumor size. The diagnosis of T3 was based on the thin-section CT-mediastinal window. In cases where the tumor diameter was less than 5 cm, infiltration was considered positive if the contact between the tumor and the parietal pleura was more than 3 cm, if there was thickening of the pleura in contact with the tumor, or if the tumor clearly formed a mass on the chest wall. In addition, nodular lesions were considered as intrapulmonary metastases if they were found in the same lobe by the CT-lung field, except for nodules with no change over time, or in cases with obvious calcification. For the diagnosis of N-factor, lymph node metastasis was considered positive if the short diameter was 10 mm or more in thin-section CT-mediastinal window.
One nuclear medicine physician (KI), board-certified in diagnostic radiology with 16 years of experience, interpreted all PET/MRI images. An FDG–avid lesion was defined as focal, abnormally increased 18F-FDG uptake versus background, with or without a corresponding anatomic lesion on the MR scan and suggestive of lung lesion and metastasis. PET/MRI images were analyzed using PET VCAR software by visually examining all the images on a computer display at the workstation (Advantage Workstation; GE Healthcare). The maximum standardized uptake value (SUVmax) of each lesion was measured visually using a circular region-of-interest. Tumor size was assessed by measuring the maximum diameter on the axial section of the MRI T2-weighted image. For the evaluation of the N factor, lymph nodes with a short diameter of 10 mm or more on MRI T2-weighted images or high signal on diffusion-weighted images and visually higher abnormal accumulation than surrounding tissues coincident with the lymph nodes on PET/MRI fusion images were considered positive for lymph node metastasis.
Pathological evaluation
A board-certified pathologist reviewed hematoxylin-eosin stained, formalin-fixed paraffin-embedded specimens, histologically confirmed diagnosis, and evaluated histopathological findings of each case of lung cancer.
Statistical analysis
Statistical analyses were performed using R [16]. Data were presented as mean ± standard deviation [SD], and a P value of 0.05 or less was considered significant. κ statistics were performed to analyze the concordance of PET/MRI, CT, and pathological staging. Bland-Altman plots were used to examine the concordance of tumor size measurements performed based on PET/MRI, CT, and pathological specimens. For each type of measurement, the percentage of the relative difference between the tumor size measurements was plotted by using the average of the two measurements. The limits of agreement were then calculated by taking the mean of the percentage of relative differences between the two measurements and two standard deviations of these differences. The cutoff value of SUVmax was determined by the receiver operating characteristic curve.
Disease-free survival (DFS) was defined as the time from surgery until an occurrence of recurrence or metastasis or last follow-up visit. Overall survival (OS) was defined as the time from surgery to death from any cause or last follow-up. Surviving patients were censored at the last follow-up period. The last follow-up date for DFS and OS calculation was April 30, 2021. For DFS and OS analysis, the data were dichotomized by TNM staging (less than Stage I and higher than Stage II). DFS analysis of N-factor (N0 or higher than N1) was also performed. The log-rank (Mantel-Cox) test was used to evaluate the difference between Kaplan–Meier curves. To calculate the risk ratios and 95% confidence intervals (CIs), univariate analysis was used to identify factors associated with DFS and OS. Then, factors found to be significant by univariate analysis (P < 0.05) were entered into a Cox multivariate regression analysis model. For the univariate analysis, we used dummy variables of 1 for the following factors: age ≥ 75 y, male, ex-smoker or current smoker, and pathology other than squamous cell carcinoma. Then, forward stepwise multivariate regression analysis was performed to identify factors correlated with DFS and OS based on calculating hazard ratios (HRs) and 95% CIs.