Skip to main content

Main menu

  • Home
  • Content
    • Current
    • Ahead of print
    • Past Issues
    • JNM Supplement
    • SNMMI Annual Meeting Abstracts
    • Continuing Education
    • JNM Podcasts
  • Subscriptions
    • Subscribers
    • Institutional and Non-member
    • Rates
    • Journal Claims
    • Corporate & Special Sales
  • Authors
    • Submit to JNM
    • Information for Authors
    • Assignment of Copyright
    • AQARA requirements
  • Info
    • Reviewers
    • Permissions
    • Advertisers
  • About
    • About Us
    • Editorial Board
    • Contact Information
  • More
    • Alerts
    • Feedback
    • Help
    • SNMMI Journals
  • SNMMI
    • JNM
    • JNMT
    • SNMMI Journals
    • SNMMI

User menu

  • Subscribe
  • My alerts
  • Log in
  • My Cart

Search

  • Advanced search
Journal of Nuclear Medicine
  • SNMMI
    • JNM
    • JNMT
    • SNMMI Journals
    • SNMMI
  • Subscribe
  • My alerts
  • Log in
  • My Cart
Journal of Nuclear Medicine

Advanced Search

  • Home
  • Content
    • Current
    • Ahead of print
    • Past Issues
    • JNM Supplement
    • SNMMI Annual Meeting Abstracts
    • Continuing Education
    • JNM Podcasts
  • Subscriptions
    • Subscribers
    • Institutional and Non-member
    • Rates
    • Journal Claims
    • Corporate & Special Sales
  • Authors
    • Submit to JNM
    • Information for Authors
    • Assignment of Copyright
    • AQARA requirements
  • Info
    • Reviewers
    • Permissions
    • Advertisers
  • About
    • About Us
    • Editorial Board
    • Contact Information
  • More
    • Alerts
    • Feedback
    • Help
    • SNMMI Journals
  • View or Listen to JNM Podcast
  • Visit JNM on Facebook
  • Join JNM on LinkedIn
  • Follow JNM on Twitter
  • Subscribe to our RSS feeds
Research ArticleOncology

Associations Between Somatic Mutations and Metabolic Imaging Phenotypes in Non–Small Cell Lung Cancer

Stephen S.F. Yip, John Kim, Thibaud P. Coroller, Chintan Parmar, Emmanuel Rios Velazquez, Elizabeth Huynh, Raymond H. Mak and Hugo J.W.L. Aerts
Journal of Nuclear Medicine April 2017, 58 (4) 569-576; DOI: https://doi.org/10.2967/jnumed.116.181826
Stephen S.F. Yip
1Department of Radiation Oncology, Dana-Farber Cancer Institute, Brigham and Women’s Hospital, and Harvard Medical School, Boston, Massachusetts
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
John Kim
2Department of Radiology, University of Michigan Health System, Ann Arbor, Michigan; and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Thibaud P. Coroller
1Department of Radiation Oncology, Dana-Farber Cancer Institute, Brigham and Women’s Hospital, and Harvard Medical School, Boston, Massachusetts
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Chintan Parmar
1Department of Radiation Oncology, Dana-Farber Cancer Institute, Brigham and Women’s Hospital, and Harvard Medical School, Boston, Massachusetts
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Emmanuel Rios Velazquez
1Department of Radiation Oncology, Dana-Farber Cancer Institute, Brigham and Women’s Hospital, and Harvard Medical School, Boston, Massachusetts
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Elizabeth Huynh
1Department of Radiation Oncology, Dana-Farber Cancer Institute, Brigham and Women’s Hospital, and Harvard Medical School, Boston, Massachusetts
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Raymond H. Mak
1Department of Radiation Oncology, Dana-Farber Cancer Institute, Brigham and Women’s Hospital, and Harvard Medical School, Boston, Massachusetts
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Hugo J.W.L. Aerts
1Department of Radiation Oncology, Dana-Farber Cancer Institute, Brigham and Women’s Hospital, and Harvard Medical School, Boston, Massachusetts
3Department of Radiology, Dana-Farber Cancer Institute, Brigham and Women’s Hospital, and Harvard Medical School, Boston, Massachusetts
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Supplemental
  • Info & Metrics
  • PDF
Loading

Abstract

PET-based radiomics have been used to noninvasively quantify the metabolic tumor phenotypes; however, little is known about the relationship between these phenotypes and underlying somatic mutations. This study assessed the association and predictive power of 18F-FDG PET–based radiomic features for somatic mutations in non–small cell lung cancer patients. Methods: Three hundred forty-eight non–small cell lung cancer patients underwent diagnostic 18F-FDG PET scans and were tested for genetic mutations. Thirteen percent (44/348) and 28% (96/348) of patients were found to harbor epidermal growth factor receptor (EGFR) or Kristen rat sarcoma viral (KRAS) mutations, respectively. We evaluated 21 imaging features: 19 independent radiomic features quantifying phenotypic traits and 2 conventional features (metabolic tumor volume and maximum SUV). The association between imaging features and mutation status (e.g., EGFR-positive [EGFR+] vs. EGFR-negative) was assessed using the Wilcoxon rank-sum test. The ability of each imaging feature to predict mutation status was evaluated by the area under the receiver operating curve (AUC) and its significance was compared with a random guess (AUC = 0.5) using the Noether test. All P values were corrected for multiple hypothesis testing by controlling the false-discovery rate (FDRWilcoxon, FDRNoether) with a significance threshold of 10%. Results: Eight radiomic features and both conventional features were significantly associated with EGFR mutation status (FDRWilcoxon = 0.01–0.10). One radiomic feature (normalized inverse difference moment) outperformed all other features in predicting EGFR mutation status (EGFR+ vs. EGFR-negative, AUC = 0.67, FDRNoether = 0.0032), as well as differentiating between KRAS-positive and EGFR+ (AUC = 0.65, FDRNoether = 0.05). None of the features was associated with or predictive of KRAS mutation status (KRAS-positive vs. KRAS-negative, AUC = 0.50–0.54). Conclusion: Our results indicate that EGFR mutations may drive different metabolic tumor phenotypes that are captured in PET images, whereas KRAS-mutated tumors do not. This proof-of-concept study sheds light on genotype–phenotype interactions, using radiomics to capture and describe the phenotype, and may have potential for developing noninvasive imaging biomarkers for somatic mutations.

  • radiomics
  • somatic mutation
  • PET imaging
  • phenotype

Non–small cell lung cancer (NSCLC) accounts for more than 80% of lung cancers and affects more than 220,000 people in the United States. It is the leading cause of cancer-related death, with a 5-y survival rate of only approximately 15% (1). Advances in cancer genomics have demonstrated that NSCLC is driven by somatic mutations in key oncogenes, such as epidermal growth factor receptor (EGFR) and Kristen rat sarcoma viral (KRAS) (2,3). These discoveries have led to the use and regulatory approval of EGFR-specific tyrosine kinase inhibitors for treatment of patients who harbor an EGFR mutation (4–6). Furthermore, KRAS mutants have been found to respond poorly to both tyrosine kinase inhibitors and conventional cisplatin-based chemotherapy, resulting in poor treatment outcomes (7,8). Therefore, given the driving force of mutation status in the efficacy of NSCLC therapies, identification of mutation status is crucial for selecting the most effective treatments, and ultimately, dictating patient outcomes. The current standard of care uses molecular testing to identify mutation status based on biopsies of tumor tissue or surgical resection; however, molecular testing can be limited by invasive procedures, long processing times, sampling error, and tissue samples that are not always readily available (9–11).

Noninvasive 18F-FDG PET/CT is increasingly used for imaging of glucose metabolism and is part of the standard initial workup for NSCLC patients in the United States (12,13). Furthermore, the metabolic tumor phenotype captured in PET images may be an indication of the underlying biology associated with somatic mutations. For example, a mutation in EGFR activates the Akt signaling pathway, which promotes glucose consumption in tumor cells for their continuous growth and survival (14,15). In addition, increased glucose transport and glycolysis have been observed in KRAS-mutated colorectal and pancreatic cancer cell lines (16,17). Thus, previous studies have investigated the associations between these mutations and SUV measures from PET images; however, there have been conflicting findings (18–22). These inconsistent reports may be due to the fact that simple SUV measures fail to capture the spatial relationships between image voxels, which may be more informative of the biology of these mutations and describe the degree of tumor heterogeneity (23,24). Highly heterogeneous tumors are often associated with poor prognosis because resistance to therapy can develop in subpopulations of tumor cells (25,26). Therefore, accurate quantification of tumor heterogeneity from PET images may provide important information for the identification of mutation status and precision medicine. Heterogeneity in the tumor phenotype can be quantitatively described through radiomic features (23,27,28), which use advanced mathematic models to quantify the spatial relationship between image voxels (29).

In this proof-of-concept study, we assessed the associations between PET radiomic features with somatic mutations and their ability to predict mutation status in 348 NSCLC patients. To our knowledge, this is the first study to investigate the relationship between somatic mutations and the metabolic phenotype, which may provide valuable information for developing noninvasive imaging biomarkers to complement molecular testing in determining mutation status.

MATERIALS AND METHODS

Patient Imaging

This retrospective study was conducted under a Dana-Farber/Harvard Cancer Center Institutional Review Board–approved protocol. This study included 348 patients with NSCLC who underwent a diagnostic 18F-FDG PET/CT scan before treatment between September 2003 and December 2013. Patient characteristics are shown in Table 1.

View this table:
  • View inline
  • View popup
TABLE 1

Patient Characteristics

Patients were injected with 344–921 MBq (9.3–24.9 mCi) of 18F-FDG and scanned approximately 65 min after injection on Discovery (GE Healthcare), Biograph (Siemens AG), or GEMINI TF (Philips) PET/CT scanners (Table 1). Attenuation correction was performed on the PET images using the corresponding CT images. The acquisition time was 3–5 min per bed position for a whole-body scan.

Mutation Status Analysis and Patient Classification

Tissue samples of primary tumors were acquired through biopsy or surgical resection. Somatic mutations were tested using a polymerase chain reaction–based method or PROFILE Oncomap (30). The nucleotide sequence encoding the kinase domain (exons 18–24) of EGFR and exons 2–3 of KRAS were analyzed by a polymerase chain reaction–based method and capillary gel electrophoresis. PROFILE Oncomap is a mass spectrometry genotyping technique that analyzes more than 470 unique mutations in 41 oncogenes. Fifty-one percent (178/348) of patients had their mutation status identified by polymerase chain reaction and 49% (170/348) using PROFILE Oncomap.

Of the 348 patients who were tested for EGFR mutations, 44 patients tested positively for an EGFR mutation (EGFR+), whereas 304 patients were EGFR-negative (EGFR–). Among the 44 EGFR+ patients, 19 (43%) and 2 (5%) had a L858R or L861Q substitution mutation in exon 21, respectively. Two (5%) patients had a G719C and 2 (5%) had a G719S substitution mutation in exon 18. A deletion mutation in exon 19 was found in 19 (43%) EGFR+ patients. Three hundred seventeen patients were tested for KRAS mutations. Although 96 patients harbored a KRAS mutation (KRAS-positive [KRAS+]), 221 patients were KRAS-negative (KRAS–). Of the 96 KRAS+ patients, 47 (49%), 48 (50%), and 1 (1%) had a substitution mutation in exon 1, 2, or 3, respectively.

One hundred eighty-five patients tested negatively for EGFR and KRAS mutations (EGFR– and KRAS–). Thirty-one patients were tested only for an EGFR mutation, but not KRAS, because of the lack of medical necessity. In this study, we performed 3 comparisons: EGFR+ versus EGFR–, EGFR+ versus KRAS+, and KRAS+ versus KRAS–. A subset analysis was also performed on patients with adenocarcinoma histology (251 patients). Of these 251 patients, 36 and 84 patients harbored EGFR and KRAS mutations, respectively.

PET Feature Extraction and Selection

The metabolic tumor volume (MTV) was delineated on PET images with a fixed SUV threshold of 2.5. A nuclear medicine board-certified radiologist was masked from the mutation data and reviewed the delineations of all 348 MTVs. All radiomic features were implemented in MATLAB (The MathWorks Inc.) and computed within the MTV. Before computation of the radiomic features, the image voxel intensities were resampled into equally spaced bin widths of 0.1 (31).

In total, 68 radiomic features were extracted from the PET images, including 7 shape features (29), 13 histogram-based features (29), and 48 texture features. The texture features included 22 gray level cooccurrence matrix (32), 11 run length matrix (33), 10 size zone matrix (34), and 5 neighborhood gray-tone difference matrix (35) features (Supplemental Table 1; supplemental materials are available at http://jnm.snmjournals.org). In addition, 5 conventional features (MTV, SUVmax, SUVpeak, SUVmean, and SUVtotal) were computed for comparison with the radiomic features.

Radiomics Feature Selection

The Spearman correlation coefficient (R) was used to assess the correlation between all radiomic and conventional features in R software (version 3.2). Feature pairs with a |R| greater than 0.95 were considered to be highly correlated and likely to provide redundant rather than complementary information about the mutation status. Features with the highest average |R| were excluded. As a result, 21 features (19 radiomic [1 shape feature, 3 histogram-based features, and 15 texture features] and 2 conventional) were included in the analysis. The selected features are shown in Table 2. All other 51 features were correlated to at least 1 of these 21 features with an R greater than 0.95 (Supplemental Figs. 1 and 2). A Pearson correlation coefficient was used to assess the correlation between the selected 21 features, which ranged from 0.05 to 0.95 as shown in Supplemental Figure 2.

View this table:
  • View inline
  • View popup
TABLE 2

Description of Selected PET-Based Features That Were Included in Analysis

Data Analysis

The association of PET features with mutation status was evaluated by comparing the distribution of each feature within the following groups: EGFR+ versus EGFR–, KRAS+ versus KRAS–, and EGFR+ versus KRAS+. The Wilcoxon rank-sum test was used to determine whether there was a significant difference in the feature values between each mutation status. All P values were corrected for multiple hypothesis testing using the Benjamini–Hochberg method (false-discovery rate [FDR]) (36) (FDRWilcoxon), with a significance threshold of 10%.

Univariate analysis was performed with R software (version 3.2) using the pROC and survcomp packages from Bioconductor (37). We evaluated the performance of each feature in classifying patients according to their mutation status (i.e., EGFR+ from EGFR–, KRAS+ from KRAS–, and EGFR+ from KRAS+). The performance was quantified using the AUC. The AUC is interpreted as the probability of correctly classifying the patients into different mutation categories and ranges from 0.5 to 1.0, where a value of 0.5 indicates a random guess and a value of 1.0 indicates a perfect classification. The Noether test was used to determine whether the AUCs were significantly greater than 0.5. P values were FDR-corrected (FDRNoether), with a significance threshold of 10%.

To assess whether the radiomic features significantly outperformed the conventional features in predicting mutation status, the receiver-operating-characteristic (ROC) curves of the radiomic features with the highest AUCs were compared with the ROC curves of MTV and SUVmax using a bootstrap test (38). P values of less than 0.05 were considered statistically significant.

RESULTS

This study assessed the association of PET-based radiomic features with NSCLC mutation status and evaluated their ability to predict mutation status in a large patient cohort. Of 348 patients who were tested for somatic mutations, 13% (44/348) of patients were EGFR+, whereas 28% (96/348) of patients were KRAS+ (Table 1). Most patients in the total cohort were women (61%), Caucasian (91%), and current/former smokers (82%) and had adenocarcinoma (72%). Approximately 40% and 60% of patients in the total cohort were overall stage I/II and stage III/IV, respectively. We observed that patients who harbored an EGFR mutation (EGFR+) generally had smaller and more compact tumors with lower 18F-FDG uptake than EGFR– (Figs. 1 and 2).

FIGURE 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 1.

From left to right are patients with EGFR mutation, KRAS mutation, and EGFR– and KRAS– tumors. Stage I and III tumors are shown in top and bottom rows, respectively. Arrows indicate locations of lung tumors.

FIGURE 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 2.

Comparison of PET features between mutation statuses. Wilcoxon rank-sum test was used to determine whether there was a significant difference in the PET feature between the mutation statuses. FDRWilcoxon ≤ 0.10 is indicated by *. Values of all PET features were normalized using z-transformation. Entries in columns of EGFR+ vs. EGFR–, EGFR+ vs. KRAS+, and KRAS+ vs. KRAS− represent differences of medians of transformed measures. For example, in EGFR+ vs. EGFR– column, entry values < 0 indicate that median value of PET feature for EGFR+ is lower than EGFR–. LGSZE = low gray small zone emphasis; LGLZE = low gray large zone emphasis; LRE = long run emphasis; LRLGE = long run low gray emphasis; LZE = large zone emphasis; RunPct = run percentage; SZE = small zone emphasis; SZV = size zone variability.

The total set of PET features extracted from the tumor volume was reduced to a small subset of descriptive features using an unsupervised selection method that reduced redundancy among the features. This method resulted in 21 features (19 radiomic and 2 conventional) that were analyzed for their associations with and predictive power for mutation status.

Comparison of PET Features Between Mutations

To assess the associations between PET features and mutation status, we compared the distribution of values of each feature between the following groups: EGFR+ versus EGFR–, EGFR+ versus KRAS+, and KRAS+ versus KRAS– (Fig. 2). A significant difference between the medians of the feature values for each mutation status indicated an association between the feature and mutation. Significant associations with 8 radiomic features (compactness, SUVmin, MaxProb, InvDiffmomnor, long run low gray emphasis, size zone variability, low gray size zone emphasis, and busyness), as well as both conventional features (MTV, SUVmax) and EGFR mutation status was observed (EGFR+ vs. EGFR–, FDRWilcoxon = 0.01–0.10). A significant association was found in the radiomic feature normalized inverse difference moment (InvDiffmomnor) between EGFR+ and KRAS+ (FDRWilcoxon = 0.08). None of the features was found to be significantly different between KRAS+ and KRAS– (FDRWilcoxon ≥ 0.92). Compactness was significantly associated with the EGFR mutation (EGFR+ vs. wild-type [WT], FDRWilcoxon = 0.08) (Supplemental Fig. 3). KRAS mutants demonstrated no distinctive imaging features. In particular, none of the PET features was significantly associated with the KRAS mutation (FDRWilcoxon ≥ 0.67) (Supplemental Fig. 3).

Prediction of Mutation Status

The 21 PET features were evaluated for their ability to predict mutation status and were assessed by the AUC. Eight radiomic features (compactness, SUVmin, MaxProb, InvDiffmomnor, long run low gray emphasis, size zone variability, low gray size zone emphasis, and busyness) and both conventional features (MTV, SUVmax) were predictive of EGFR mutation status, with AUCs ranging from 0.59 to 0.67 (FDRNoether = 0.0032–0.09) (Fig. 3). The radiomic feature InvDiffmomnor (AUC = 0.67, FDRNoether = 0.0032) had the highest AUC for discriminating EGFR+ from EGFR–, and long run low gray emphasis (FDRNoether = 0.06) and busyness (FDRNoether = 0.09) had AUCs of 0.59. Five radiomic features, including compactness, MaxProb, InvDiffmomnor, size-zone-variability, and low gray size zone emphasis, were able to significantly discriminate EGFR+ from the EGFR-KRAS WT (AUC = 0.61–0.67, FDRNoether = 0.003–0.04) (Supplemental Fig. 4). In addition, conventional features MTV and SUVmax were significantly predictive of EGFR+ and WT differentiation (AUC = 0.62–0.64, FDRNoether = 0.01–0.02). InvDiffmomnor outperformed all other measures in EGFR mutation prediction.

FIGURE 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 3.

AUC. * indicates that AUC is significantly > 0.50 (random guessing) assessed with Noether test (FDRNoether ≤ 0.10). Many of the features significantly predict EGFR+ tumors; however, they are not able to predict KRAS+ tumors. LGSZE = low gray small zone emphasis; LGLZE = low gray large zone emphasis; LRE = long run emphasis; LRLGE = long run low gray emphasis; LZE = large zone emphasis; RunPct = run percentage; SZE = small zone emphasis; SZV = size zone variability.

Only 1 radiomic feature, InvDiffmomnor (AUC = 0.65, FDRNoether = 0.05), was able to significantly differentiate EGFR+ from KRAS+ (Fig. 3). In particular, none of the conventional PET features was able to differentiate EGFR+ from KRAS+ (AUC = 0.56–0.60, FDRNoether ≥ 0.29). All other features had AUCs ranging from 0.50 to 0.61 (FDRNoether ≥ 0.29) for EGFR+ and KRAS+ discrimination. None of the features was significantly predictive of KRAS mutation status (KRAS+ vs. KRAS–, AUC < 0.55, FDRNoether ≥ 0.92) (Fig. 3) or discriminating KRAS+ from WT (AUC = 0.50–0.57, FDRNoether ≥ 0.65) (Supplemental Fig. 4).

The predictive performances of the significant radiomic features and conventional features for mutation status were compared directly to assess the value of radiomics over conventional features. The radiomic feature InvDiffmomnor significantly outperformed the conventional feature MTV in differentiating EGFR+ from EGFR– (P = 0.02) but did not significantly outperform the conventional feature SUVmax (P = 0.14).

In patients with only adenocarcinoma histology, InvDiffmomnor (AUC = 0.66, FDRNoether = 0.01) was observed to outperform all other features (AUC ≤ 0.64, FDRNoether ≥ 0.03) in differentiating EGFR+ from EGFR– (Supplemental Fig. 5), as observed with the full cohort. Adenocarcinoma EGFR mutants had a significantly higher SUVmax than WTs (FDRWilcoxon = 0.06).

DISCUSSION

Identifying mutation status in NSCLC patients is an important component of selecting an optimal treatment plan for the patient. Clinically, a patient’s mutation status is identified using molecular testing based on pathologic examination of tumor tissues; however, these data may not always be readily available for all patients. We hypothesized that the genotype of tumors may be reflected in the phenotype that is captured in medical images. PET imaging is commonly used in the clinic for cancer management and can provide a noninvasive depiction of the tumor metabolic phenotype. We investigated the association between somatic mutations and 19 novel PET radiomic features and 2 conventional features (MTV and SUVmax) in 348 NSCLC patients. We observed that radiomic features are strongly associated with EGFR mutation status and could significantly discriminate EGFR+ from EGFR– and EGFR+ from KRAS+.

We found that EGFR mutants demonstrated distinctive imaging features compared with tumors without an EGFR mutation. For example, the metabolic region of EGFR+ tumors was smaller and had a lower SUVmax than EGFR– tumors. Because of their low SUV uptake and smaller MTV, our results are consistent with the notion that EGFR mutants tend to be more indolent than other lung cancer types (39,40). In addition, 8 radiomic features were significantly associated with and predictive of EGFR mutation status, where InvDiffmomnor was the most predictive for EGFR mutation status and significantly outperformed the conventional measure MTV (P ≤ 0.02). Furthermore, InvDiffmomnor was also found to outperform all other PET-based features in predicting EGFR mutations in patients with adenocarcinoma histology and was also predictive for differentiating EGFR+ from KRAS+ (Supplemental Fig. 5). However, only 5 features were shown to significantly discriminate EGFR+ from WT (Supplemental Figs. 3 and 4). This may be due to the fact that only 31 patients were tested only for EGFR mutations, but not KRAS mutations, because of a lack of medical necessity. Seventy-four percent (23/31) of these patients were EGFR–, whereas their KRAS mutation status was unknown. As it was unclear if those 23 patients were also KRAS–, they were excluded from the EGFR+ and WT differentiation analysis, which may explain why SUVmin, busyness, and long run low gray emphasis could significantly differentiate EGFR+ from EGFR–, but not EGFR+ from WT.

Our results showed that a lower SUVmax was associated with EGFR+ tumors. Previous studies have also shown a substantial association between SUVmax and EGFR mutation status in NSCLC patients, where lower SUVmax was associated with EGFR+ (18,20,41); however, there have been conflicting conclusions with other studies. Conversely, studies by Huang et al. (21) and Lee et al. (42) identified that higher SUV measures, such as SUVmax, were a strong predictor of EGFR+ tumors. In other studies, Lee et al. (43) and Chung et al. (44) found that none of their PET features (e.g., SUVmax, MTV, and total glycolysis lesion) was significantly associated with or predicted EGFR mutation status.

The differences observed in these aforementioned studies may be attributed to the patient demographics for each study. Our study contained the largest patient cohort (348 patients) compared with the other studies (206 patients in Lee et al. (43), 106 patients in Chung et al. (44), 77 patients in Huang et al. (21), and 214 patients in Lee et al. (42)). Furthermore, the previous studies were from Asian institutions whereas our study was from a North American cohort, which resulted in more than 90% of our patients being Caucasian (Table 1). Another study that had a predominantly Caucasian cohort was reported by Mak et al. (18), who reported that a lower SUVmax was associated with EGFR mutants, which is consistent with our current study.

The frequency of EGFR mutations has been found to be substantially higher in Asian countries than Western countries (45). Therefore, the effect of ethnicity is reflected in the frequency of EGFR+ mutations in each cohort, where the previous studies had a higher percentage of EGFR+ mutants than in our study (23% in the cohort of Lee et al. (43), 40% in Chung et al. (44), 64% in Huang et al. (21), and 24% in Lee et al. (42) vs. only 13% in our study). However, despite the differences between each study, it is clear that the conclusions on using conventional PET-based features for prediction of mutation status remain debatable. The limitation of conventional PET features may be that they do not adequately quantify the tumor phenotype that reflects mutation status. Using a radiomics approach with greater descriptive power to quantify the tumor phenotype, we found that not only were many radiomic features predictive of EGFR mutation status, they also could be better predictors than conventional PET features. Although this was an initial exploratory study to demonstrate that radiomics could be applied for mutation status prediction using PET images, these findings need to be further investigated in multiple large datasets acquired from different countries.

None of the conventional or radiomic features was associated with or predictive of KRAS mutation status (KRAS+ vs. KRAS−). We found that KRAS+ tumors tended to be bigger and had greater values in SUV measures than EGFR+ tumors, suggesting that KRAS+ tumors may be more aggressive (Figs. 1 and 2). Only 1 radiomic feature (InvDiffmomnor) could significantly discriminate KRAS+ from EGFR+, whereas other features, including conventional features, were unable to.

Other groups have also investigated the ability of conventional PET features to identify KRAS mutation status. Caicedo et al. (19) investigated whether SUV measures, including maximum, peak, and average SUV, could predict EGFR and KRAS mutation status in 102 patients with NSCLC. Unlike our results, they found that KRAS+ tumors had significantly higher values in SUV measures than both EGFR mutants and KRAS-EGFR WTs. One reason for this difference between our study and their study may be because a proportion of the KRAS– tumors in our study were also EGFR+ (36/221 [16%] KRAS– tumors were EGFR mutants). In our analysis of KRAS+ and WT, we found that all the features performed poorly in prediction for the KRAS mutation. Moreover, we also found that the feature distributions between KRAS+ and KRAS– were not significantly different for any of the PET features. However, these differences may also be attributed to the different cohort size because we investigated the association between PET imaging features and mutation status in 348 tumors, whereas the study by Caicedo et al. (19) investigated 102 tumors.

There are several limitations of our study. First, because of the retrospective nature of the study, the imaging protocols were not standardized for all the patients, resulting in different acquisition and reconstruction parameters (Table 1). Although these differences may lead to variability in quantification of SUV (46–48), the distributions of scanners used for PET/CT image acquisition were similar for all patient groups (Table 1). Second, partial-volume effects resulting from limited PET spatial resolution may lead to an underestimation of the metabolic tumor region and SUV measures (49). Furthermore, because of insufficient data acquisition and limited reconstruction techniques, respiratory motion can induce image blurring in the static (3-dimensional PET) image acquisition, subsequently affecting the quantification of the PET features (50). Despite these limitations, many radiomic features were still found to be significantly associated with and predictive of EGFR mutation status. Whether PET standardization, correction of partial-volume effects, and the use of respiratory-gated PET/CT imaging can improve the value of PET features in mutation identification needs to be further investigated. Although this study focused on the use of radiomic features extracted from the untreated primary tumors for mutation prediction, in the future, it would be interesting to also investigate the predictive value of features extracted from metastatic lesions.

CONCLUSION

Tumor metabolic phenotypes that are driven by EGFR mutations could be quantified by radiomic features. Several radiomic features were strongly associated with the EGFR mutation status; in particular, InvDiffmomnor outperformed all other PET features (radiomic and conventional) in EGFR mutation prediction (EGFR+ vs. EGFR− and EGFR+ vs. KRAS+). All PET features were poorly associated with KRAS mutations, potentially indicating that tumors of KRAS mutants demonstrated weak distinctive imaging features.

Radiomic features extracted from diagnostic PET images have potential for identifying EGFR mutations in NSCLC. Our study may help develop an imaging biomarker to noninvasively identify EGFR mutation status using PET imaging to complement, but not to replace, molecular testing. Prospective studies with masked mutation status and independent datasets will be needed to further validate the predictive power of the radiomic features. Furthermore, future studies will need to investigate how the EGFR mutation gives rise to certain phenotypic traits that are quantified by radiomic features.

DISCLOSURE

We acknowledge support from the National Institutes of Health (award nos. U01CA190234 and U24CA194354) and research seed funding grant from the American Association of Physicists in Medicine. No other potential conflict of interest relevant to this article was reported.

Acknowledgments

We thank the PROFILE team for their help with somatic mutation testing.

Footnotes

  • ↵* Contributed equally to this work.

  • Published online Sep. 29, 2016.

  • © 2017 by the Society of Nuclear Medicine and Molecular Imaging.

REFERENCES

  1. 1.↵
    1. Siegel RL,
    2. Miller KD,
    3. Jemal A
    . Cancer statistics, 2015. CA Cancer J Clin. 2015;65:5–29.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Pao W,
    2. Girard N
    . New driver mutations in non-small-cell lung cancer. Lancet Oncol. 2011;12:175–180.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Pao W,
    2. Hutchinson KE
    . Chipping away at the lung cancer genome. Nat Med. 2012;18:349–351.
    OpenUrlCrossRefPubMed
  4. 4.↵
    1. Paez JG,
    2. Jänne PA,
    3. Lee JC,
    4. et al
    . EGFR mutations in lung cancer: correlation with clinical response to gefitinib therapy. Science. 2004;304:1497–1500.
    OpenUrlAbstract/FREE Full Text
  5. 5.
    1. Maemondo M,
    2. Inoue A,
    3. Kobayashi K,
    4. et al
    . Gefitinib or chemotherapy for non–small-cell lung cancer with mutated EGFR. N Engl J Med. 2010;362:2380–2388.
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Zhou C,
    2. Wu Y-L,
    3. Chen G,
    4. et al
    . Erlotinib versus chemotherapy as first-line treatment for patients with advanced EGFR mutation-positive non-small-cell lung cancer (OPTIMAL, CTONG-0802): a multicentre, open-label, randomised, phase 3 study. Lancet Oncol. 2011;12:735–742.
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. Loriot Y,
    2. Mordant P,
    3. Deutsch E,
    4. Olaussen KA,
    5. Soria J-C
    . Are RAS mutations predictive markers of resistance to standard chemotherapy? Nat Rev Clin Oncol. 2009;6:528–534.
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. Pao W,
    2. Wang TY,
    3. Riely GJ,
    4. et al
    . KRAS mutations and primary resistance of lung adenocarcinomas to gefitinib or erlotinib. PLoS Med. 2005;2:e17.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Lindeman NI,
    2. Cagle PT,
    3. Beasley MB,
    4. et al
    . Molecular testing guideline for selection of lung cancer patients for EGFR and ALK tyrosine kinase inhibitors: guideline from the college of American pathologists, international association for the study of lung cancer, and association for molecular pathology. J Mol Diagn. 2013;15:415–453.
    OpenUrlCrossRefPubMed
  10. 10.
    1. Pao W,
    2. Ladanyi M
    . Epidermal growth factor receptor mutation testing in lung cancer: searching for the ideal method. Clin Cancer Res. 2007;13:4954–4955.
    OpenUrlFREE Full Text
  11. 11.↵
    1. Gerlinger M,
    2. Rowan AJ,
    3. Horswell S,
    4. et al
    . Intratumor heterogeneity and branched evolution revealed by multiregion sequencing. N Engl J Med. 2012;366:883–892.
    OpenUrlCrossRefPubMed
  12. 12.↵
    1. Cuaron J,
    2. Dunphy M,
    3. Rimner A
    . Role of FDG-PET scans in staging, response assessment, and follow-up care for non-small cell lung cancer. Front Oncol. 2013;2:208.
    OpenUrlPubMed
  13. 13.↵
    1. Weber WA,
    2. Petersen V,
    3. Schmidt B,
    4. et al
    . Positron emission tomography in non–small-cell lung cancer: prediction of response to chemotherapy by quantitative assessment of glucose use. J Clin Oncol. 2003;21:2651–2657.
    OpenUrlAbstract/FREE Full Text
  14. 14.↵
    1. Gan Y,
    2. Shi C,
    3. Inge L,
    4. Hibner M,
    5. Balducci J,
    6. Huang Y
    . Differential roles of ERK and Akt pathways in regulation of EGFR-mediated signaling and motility in prostate cancer cells. Oncogene. 2010;29:4947–4958.
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Elstrom RL,
    2. Bauer DE,
    3. Buzzai M,
    4. et al
    . Akt stimulates aerobic glycolysis in cancer cells. Cancer Res. 2004;64:3892–3899.
    OpenUrlAbstract/FREE Full Text
  16. 16.↵
    1. Yun J,
    2. Rago C,
    3. Cheong I,
    4. et al
    . Glucose deprivation contributes to the development of KRAS pathway mutations in tumor cells. Science. 2009;325:1555–1559.
    OpenUrlAbstract/FREE Full Text
  17. 17.↵
    1. Ying H,
    2. Kimmelman Alec C,
    3. Lyssiotis Costas A,
    4. et al
    . Oncogenic Kras maintains pancreatic tumors through regulation of anabolic glucose metabolism. Cell. 2012;149:656–670.
    OpenUrlCrossRefPubMed
  18. 18.↵
    1. Mak RH,
    2. Digumarthy SR,
    3. Muzikansky A,
    4. et al
    . Role of 18F-dluorodeoxyglucose positron emission tomography in predicting epidermal growth factor receptor mutations in non-small cell lung cancer. Oncologist. 2011;16:319–326.
    OpenUrlAbstract/FREE Full Text
  19. 19.↵
    1. Caicedo C,
    2. Garcia-Velloso M,
    3. Lozano M,
    4. et al
    . Role of [18F]FDG PET in prediction of KRAS and EGFR mutation status in patients with advanced non-small-cell lung cancer. Eur J Nucl Med Mol Imaging. 2014;41:2058–2065.
    OpenUrlCrossRefPubMed
  20. 20.↵
    1. Na II,
    2. Byun BH,
    3. Kim KM,
    4. et al
    . 18F-FDG uptake and EGFR mutations in patients with non-small cell lung cancer: a single-institution retrospective analysis. Lung Cancer. 2010;67:76–80.
    OpenUrlCrossRefPubMed
  21. 21.↵
    1. Huang C-T,
    2. Yen R-F,
    3. Cheng M-F,
    4. et al
    . Correlation of F-18 fluorodeoxyglucose-positron emission tomography maximal standardized uptake value and EGFR mutations in advanced lung adenocarcinoma. Med Oncol. 2010;27:9–15.
    OpenUrlCrossRefPubMed
  22. 22.↵
    1. Kawada K,
    2. Nakamoto Y,
    3. Kawada M,
    4. et al
    . Relationship between 18F-dluorodeoxyglucose accumulation and KRAS/BRAF mutations in colorectal cancer. Clin Cancer Res. 2012;18:1696–1703.
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    1. Cook GR,
    2. Siddique M,
    3. Taylor B,
    4. Yip C,
    5. Chicklore S,
    6. Goh V
    . Radiomics in PET: principles and applications. Clin Transl Imaging. 2014;2:269–276.
    OpenUrl
  24. 24.↵
    1. Rahim MK,
    2. Kim S,
    3. So H,
    4. et al
    . Recent trends in PET image interpretations using volumetric and texture-based quantification methods in nuclear oncology. Nucl Med Mol Imaging. 2014;48:1–15.
    OpenUrlCrossRef
  25. 25.↵
    1. Burrell RA,
    2. McGranahan N,
    3. Bartek J,
    4. Swanton C
    . The causes and consequences of genetic heterogeneity in cancer evolution. Nature. 2013;501:338–345.
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. Marusyk A,
    2. Almendro V,
    3. Polyak K
    . Intra-tumour heterogeneity: a looking glass for cancer? Nat Rev Cancer. 2012;12:323–334.
    OpenUrlCrossRefPubMed
  27. 27.↵
    1. Henriksson E,
    2. Kjellen E,
    3. Wahlberg P,
    4. Ohlsson T,
    5. Wennerberg J,
    6. Brun E
    . 2-deoxy-2-[18F]fluoro-d-glucose uptake and correlation to intratumoral heterogeneity. Anticancer Res. 2007;27:2155–2159.
    OpenUrlAbstract/FREE Full Text
  28. 28.↵
    1. Yip SS,
    2. Aerts JWLH
    . Applications and limitations of radiomics. Phys Med Biol. 2016;61:R150–R166.
    OpenUrlPubMed
  29. 29.↵
    1. Aerts HJ,
    2. Velazquez ER,
    3. Leijenaar RT,
    4. et al
    . Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach. Nat Commun. 2014;5:4006.
    OpenUrlPubMed
  30. 30.↵
    1. MacConaill LE,
    2. Garcia E,
    3. Shivdasani P,
    4. et al
    . Prospective enterprise-level molecular genotyping of a cohort of cancer patients. J Mol Diagn. 2014;16:660–672.
    OpenUrlCrossRefPubMed
  31. 31.↵
    1. Leijenaar RT,
    2. Nalbantov G,
    3. Carvalho S,
    4. et al
    . The effect of SUV discretization in quantitative FDG-PET radiomics: the need for standardized methodology in tumor texture analysis. Sci Rep. 2015;5:11075.
    OpenUrlCrossRefPubMed
  32. 32.↵
    1. Haralick RM,
    2. Shanmugam K,
    3. Dinstein IH
    . Textural features for image classification. IEEE Trans Syst Man Cybern Syst. 1973;SMC-3:610–621.
    OpenUrlCrossRef
  33. 33.↵
    1. Galloway MM
    . Texture analysis using gray level run lengths. Comput Vision Graph. 1975;4:172–179.
    OpenUrl
  34. 34.↵
    1. Thibault G,
    2. Fertil B,
    3. Navarro C,
    4. et al
    . Shape and texture indexes application to cell nuclei classification. Int J Pattern Recogn. 2013;27:1357002.
    OpenUrl
  35. 35.↵
    1. Amadasun M,
    2. King R
    . Textural features corresponding to textural properties. IEEE Trans Syst Man Cybern Syst. 1989;19:1264–1274.
    OpenUrl
  36. 36.↵
    1. Benjamini Y,
    2. Hochberg Y
    . Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol. 1995;57:289–300.
    OpenUrl
  37. 37.↵
    1. Gentleman RC,
    2. Carey VJ,
    3. Bates DM,
    4. et al
    . Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5:R80.
    OpenUrlCrossRefPubMed
  38. 38.↵
    1. Hanley JA,
    2. McNeil BJ
    . A method of comparing the areas under receiver operating characteristic curves derived from the same cases. Radiology. 1983;148:839–843.
    OpenUrlCrossRefPubMed
  39. 39.↵
    1. Chmielecki J,
    2. Foo J,
    3. Oxnard GR,
    4. et al
    . Optimization of dosing for EGFR-mutant non–small cell lung cancer with evolutionary cancer modeling. Sci Transl Med. 2011;3:90ra59.
    OpenUrlAbstract/FREE Full Text
  40. 40.↵
    1. Rekhtman N,
    2. Ang DC,
    3. Riely GJ,
    4. Ladanyi M,
    5. Moreira AL
    . KRAS mutations are associated with solid growth pattern and tumor-infiltrating leukocytes in lung adenocarcinoma. Mod Pathol. 2013;26:1307–1319.
    OpenUrlCrossRefPubMed
  41. 41.↵
    1. Choi Y-J,
    2. Cho B,
    3. Jeong Y,
    4. et al
    . Correlation between 18F-fluorodeoxyglucose uptake and epidermal growth factor receptor mutations in advanced lung cancer. Nucl Med Mol Imaging. 2012;46:169–175.
    OpenUrl
  42. 42.↵
    1. Lee Y,
    2. Lee H-J,
    3. Kim YT,
    4. et al
    . Imaging characteristics of stage I non-small cell lung cancer on CT and FDG-PET: relationship with epidermal growth factor receptor protein expression status and survival. Korean J Radiol. 2013;14:375–383.
    OpenUrl
  43. 43.↵
    1. Lee SM,
    2. Bae SK,
    3. Jung SJ,
    4. Kim CK
    . FDG uptake in non–small cell lung cancer is not an independent predictor of EGFR or KRAS mutation status: a retrospective analysis of 206 patients. Clin Nucl Med. 2015;40:950–958.
    OpenUrl
  44. 44.↵
    1. Chung HW,
    2. Lee KY,
    3. Kim HJ,
    4. Kim WS,
    5. So Y
    . FDG PET/CT metabolic tumor volume and total lesion glycolysis predict prognosis in patients with advanced lung adenocarcinoma. J Cancer Res Clin Oncol. 2014;140:89–98.
    OpenUrlCrossRefPubMed
  45. 45.↵
    1. Dearden S,
    2. Stevens J,
    3. Wu Y-L,
    4. Blowers D
    . Mutation incidence and coincidence in non small-cell lung cancer: meta-analyses by ethnicity and histology (mutMap). Ann Oncol. 2013;24:2371–2376.
    OpenUrlAbstract/FREE Full Text
  46. 46.↵
    1. Boellaard R
    . Standards for PET image acquisition and quantitative data analysis. J Nucl Med. 2009;50:11S–20S.
    OpenUrlAbstract/FREE Full Text
  47. 47.
    1. Galavis PE,
    2. Hollensen C,
    3. Jallow N,
    4. Paliwal B,
    5. Jeraj R
    . Variability of textural features in FDG PET images due to different acquisition modes and reconstruction parameters. Acta Oncol. 2010;49:1012–1016.
    OpenUrlCrossRefPubMed
  48. 48.↵
    1. Leijenaar RTH,
    2. Carvalho S,
    3. Velazquez ER,
    4. et al
    . Stability of FDG-PET radiomics features: an integrated analysis of test-retest and inter-observer variability. Acta Oncol. 2013;52:1391–1397.
    OpenUrlCrossRefPubMed
  49. 49.↵
    1. Soret M,
    2. Bacharach SL,
    3. Buvat I
    . Partial-volume effect in PET tumor imaging. J Nucl Med. 2007;48:932–945.
    OpenUrlAbstract/FREE Full Text
  50. 50.↵
    1. Yip S,
    2. McCall K,
    3. Aristophanous M,
    4. Chen AB,
    5. Aerts HJWL,
    6. Berbeco R
    . Comparison of texture features derived from static and respiratory-gated PET images in non-small cell lung cancer. PLoS One. 2014;9:e115510.
    OpenUrlCrossRefPubMed
  • Received for publication July 27, 2016.
  • Accepted for publication September 3, 2016.
PreviousNext
Back to top

In this issue

Journal of Nuclear Medicine: 58 (4)
Journal of Nuclear Medicine
Vol. 58, Issue 4
April 1, 2017
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
Print
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for your interest in spreading the word on Journal of Nuclear Medicine.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Associations Between Somatic Mutations and Metabolic Imaging Phenotypes in Non–Small Cell Lung Cancer
(Your Name) has sent you a message from Journal of Nuclear Medicine
(Your Name) thought you would like to see the Journal of Nuclear Medicine web site.
Citation Tools
Associations Between Somatic Mutations and Metabolic Imaging Phenotypes in Non–Small Cell Lung Cancer
Stephen S.F. Yip, John Kim, Thibaud P. Coroller, Chintan Parmar, Emmanuel Rios Velazquez, Elizabeth Huynh, Raymond H. Mak, Hugo J.W.L. Aerts
Journal of Nuclear Medicine Apr 2017, 58 (4) 569-576; DOI: 10.2967/jnumed.116.181826

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
Associations Between Somatic Mutations and Metabolic Imaging Phenotypes in Non–Small Cell Lung Cancer
Stephen S.F. Yip, John Kim, Thibaud P. Coroller, Chintan Parmar, Emmanuel Rios Velazquez, Elizabeth Huynh, Raymond H. Mak, Hugo J.W.L. Aerts
Journal of Nuclear Medicine Apr 2017, 58 (4) 569-576; DOI: 10.2967/jnumed.116.181826
Twitter logo Facebook logo LinkedIn logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One
Bookmark this article

Jump to section

  • Article
    • Abstract
    • MATERIALS AND METHODS
    • RESULTS
    • DISCUSSION
    • CONCLUSION
    • DISCLOSURE
    • Acknowledgments
    • Footnotes
    • REFERENCES
  • Figures & Data
  • Supplemental
  • Info & Metrics
  • PDF

Related Articles

  • This Month in JNM
  • PubMed
  • Google Scholar

Cited By...

  • Performance and Robustness of Machine Learning-based Radiomic COVID-19 Severity Prediction
  • Evolutionary dynamics at the tumor edge reveal metabolic imaging biomarkers
  • Radiomics Analysis of Clinical Myocardial Perfusion Stress SPECT Images to Identify Coronary Artery Calcification
  • Introduction to Radiomics
  • Metabolic Biomarker-Based BRAFV600 Mutation Association and Prediction in Melanoma
  • Identifying relationships between imaging phenotypes and lung cancer-related mutation status: EGFR and KRAS
  • A Postreconstruction Harmonization Method for Multicenter Radiomic Studies in PET
  • Somatic Mutations Drive Distinct Imaging Phenotypes in Lung Cancer
  • Google Scholar

More in this TOC Section

Oncology

  • The Evolving Role of FDG-PET/CT in Assessing Primary Cutaneous Lymphoma
  • Molecular imaging of Renal clear cell carcinoma subtypes as metabolic diseases
  • Role of PET/CT in the management of multiple myeloma
Show more Oncology

Clinical

  • The Evolving Role of FDG-PET/CT in Assessing Primary Cutaneous Lymphoma
  • Molecular imaging of Renal clear cell carcinoma subtypes as metabolic diseases
  • Role of PET/CT in the management of multiple myeloma
Show more Clinical

Similar Articles

Keywords

  • radiomics
  • somatic mutation
  • PET imaging
  • Phenotype
SNMMI

© 2025 SNMMI

Powered by HighWire