Deciphering the tumor microenvironment through radiomics in non-small cell lung cancer: Correlation with immune profiles

Abstract

Growing evidence suggests that the efficacy of immunotherapy in non-small cell lung cancers (NSCLCs) is associated with the immune microenvironment within the tumor. We aimed to explore radiologic phenotyping using a radiomics approach to assess the immune microenvironment in NSCLC. Two independent NSCLC cohorts (training and test sets) were included. Single-sample gene set enrichment analysis was used to determine the tumor microenvironment, where type 1 helper T (Th1) cells, type 2 helper T (Th2) cells, and cytotoxic T cells were the targets for prediction with computed tomographic (CT) radiomic features. Multiple algorithms were in the modeling followed by final model selection. The training dataset comprised 89 NSCLCs and the test set included 60 cases of lung squamous cell carcinoma and adenocarcinoma. A total of 239 CT radiomic features were used. A linear discriminant analysis model was selected for the final model of Th2 cell group prediction. The area under the curve value of the final model on the test set was 0.684. Predictors of the linear discriminant analysis model were skewness (total and outer pixels), kurtosis, variance (subsampled from delta [subtraction inner pixels from outer pixels]), and informational measure of correlation. The performances of radiomics on test set of Th1 and cytotoxic T cell were not accurate enough to be predictable. A radiomics approach can be used to interrogate an entire tumor in a noninvasive manner and provide added diagnostic value to identify the immune microenvironment of NSCLC, in particular, Th2 cell signatures.

About the Authors:

Hyun Jung Yoon

Contributed equally to this work with: Hyun Jung Yoon, Jun Kang, Hyunjin Park

Roles Data curation, Methodology, Writing – original draft, Writing – review & editing

Affiliations Department of Radiology and Center for Imaging Science, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul, South Korea, Department of Radiology, Veterans Health Service Medical Center, Seoul, South Korea

Jun Kang

Contributed equally to this work with: Hyun Jung Yoon, Jun Kang, Hyunjin Park

Roles Conceptualization, Data curation, Formal analysis, Methodology, Supervision, Validation, Writing – original draft

Affiliation: Department of Hospital Pathology, Seoul St. Mary’s Hospital, College of Medicine, The Catholic University of Korea, Seoul, South Korea

ORCID logo http://orcid.org/0000-0002-7967-0917

Hyunjin Park

Contributed equally to this work with: Hyun Jung Yoon, Jun Kang, Hyunjin Park

Roles Conceptualization, Data curation, Formal analysis, Methodology, Software, Visualization, Writing – original draft

Affiliations School of Electronic and Electrical Engineering, Sungkyunkwan University, Suwon, South Korea, Center for Neuroscience Imaging Research, Institute for Basic Science, Suwon, South Korea

ORCID logo http://orcid.org/0000-0001-5681-8918

Insuk Sohn

Roles Data curation, Formal analysis

Affiliation: Statistics and Data Center, Samsung Medical Center, Seoul, South Korea

Seung-Hak Lee

Roles Data curation, Formal analysis

Affiliation: Department of Electrical and Computer Engineering, Sungkyunkwan University, Suwon, South Korea

Ho Yun Lee

Roles Conceptualization, Funding acquisition, Methodology, Supervision, Validation, Visualization, Writing – review & editing

* E-mail: [email protected]

Affiliations Department of Radiology and Center for Imaging Science, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul, South Korea, Department of Health Sciences and Technology, SAIHST, Sungkyunkwan University, Seoul, South Korea

ORCID logo http://orcid.org/0000-0001-9960-5648

Introduction

Immune checkpoint blockade therapy, an anti-cancer treatment that potentiates the ability of the immune system to recognize and destroy cancer cells, has become a standard in the care of patients with non-small cell lung cancers (NSCLCs) [1–5]. The expression of the programmed death-ligand 1 (PD-L1) gene in the tumor has been associated with a response to checkpoint blockade therapy and prognosis, to some extent [6, 7], but the understanding of the complex interactions between tumors and their microenvironment remains insufficient to distinguish between patients who will respond to therapy and those who should be offered alternative treatment. Growing evidence suggests that the efficacy of immunotherapy in cancers is associated with the immune microenvironment of tumors, thereby type 1 and 2 immune responses are believed to be important predictive biomarkers [8–11]. Generally, type 1 immune responses are thought to be the most relevant to antitumor immunity, and many cancer immunotherapies are designed to augment type 1 immune responses that involve cytotoxic T cells and type 1 helper T (Th1) cells to eliminate tumors [9], whereas a type 2 immune response characterized by type 2 helper T (Th2) cells is often associated with a tumor-permissive environment [8]. Recent work revealed that the inhibition of melanoma tumor formation occurred in an activated type 2 immune response that was distinct from the classical tumor microenvironment [10].

Immunohistochemistry after tumor resection is the only technique that yields quantitative information on the immune microenvironment and immune response type assessment. Therefore, efficient and reliable noninvasive biomarkers are urgently needed. Radiomics has recently been identified as a way to use accurate quantitative imaging descriptors in line with advances in image processing techniques [12, 13]. High-dimensional imaging data allows an in-depth characterization of tumor phenotypes, with the underlying hypothesis that imaging reflects not only macroscopic but also the cellular and molecular properties of tissues. Radiomic features are complementary to biopsies and have the advantage of being noninvasive, which allows the evaluation of a tumor and its microenvironment, characterization of spatial heterogeneity, and longitudinal assessment of disease evolution [13, 14]. Since radiomic features are composed of hundreds of variables, a predictive modeling algorithm can be implemented for regularization, feature selection, and tree-based algorithms.

Accordingly, we have conducted a study to demonstrate the feasibility of radiomic prediction of the immune microenvironment by using computed tomographic (CT) imaging in NSCLC. Our ultimate goal was to identify useful predictive radiomic characteristics of the immune microenvironment and to further develop treatment strategies.

Materials and methods

Datasets

Two independent datasets were used for the training and test sets. The training dataset labeled “Lung3”, contains data from patients treated at the MAASTRO Clinic, the Netherlands [13]. The test set was a dataset from The Cancer Genome Atlas (TCGA), which included both squamous cell carcinoma and adenocarcinoma cases [15, 16]. Both training and test sets were publicly available CT image and gene expression data. Thus, no Institutional Review Board approval specific to this study was required.

Radiomic feature extraction for prediction variables

CT images were gathered from The Cancer Imaging Archive (TCIA, http://cancerimagingarchive.net/) [17]. All images from TCIA were acquired using CT scans at following variable image resolutions; For “Lung3”, mean in-plane resolution, 0.913 ± 0.204 mm (range, 0.602–1.367 mm) and mean slice thickness, 3.989 ± 1.236 mm. For “TCGA”, mean in-plane resolution, 0.718 ± 0.136 mm (range, 0.547–0.977 mm) and mean slice thickness, 3.930 ± 1.429 mm. The values following the mean are standard deviation (SD) values. In "Lung 3" cohort, 62 (69.7%) cases were scanned by contrast-enhanced CT and 27 (30.3%) patients were scanned by non-contrast CT. In "TCGA", 25 (41.7%) cases were scanned by contrast-enhanced CT and 35 (58.3%) cases were scanned by non-contrast CT. For the radiomic analysis, regions of interest (ROIs) were delineated on axial CT images to generate a volume of interest that included the entire target lesion (Fig 1). The quantitative CT analysis was based on physical, histogram-based, shape, local, filter-based, and fractal model-based features from the manually derived ROI. All radiomic features were calculated using the open-source code (PyRadiomics) [18] and in-house code using MATLAB. Features unavailable in PyRadiomics were implemented using the in-house code and there are several published articles which have been used the same software [19, 20]. Details are described in S1 File.

[Figure omitted. See PDF.]

Fig 1. Development of the tumor microenvironment prediction model.

https://doi.org/10.1371/journal.pone.0231227.g001

Immune microenvironment profiling for target variables

The training set gene expression data were generated using the Rosetta/Merck Human RSTA Custom Affymetrix 2.0 microarray (GPL15048) and stored in the Gene Expression Omnibus (GEO accession: GSE58661). The test set gene expression data were generated using RNAseq and downloaded from the Broad GDAC Firehose. The immune cell signature was estimated from gene expression data using single-sample gene set enrichment analysis (ssGSEA) to determine the immune microenvironment of lung cancer. The GSVA R package was used for the ssGSEA [21]. A set of 28 immune cell gene lists was selected from a previous study conducted by Bindea et al. [22]. We focused on Th1, Th2, and cytotoxic T cell signatures. The immune cell signature of TCGA data was initially obtained from all 1016 cases. The test set included 60 TCGA cases that were available with CT images. Principal component analysis (PCA) of immune cell signatures was done with 1016 TCGA non-small cell carcinomas (Fig 2A). The immune cell signature that gave a positive principal component 1 (PC1) score included Th1 cells, CD8 T cells, macrophages, and dendritic cells. The immune cell signature that gave PC2 positive and PC1 negative scores included Th2 cells and cancer cells. The immune cell signature that gave PC2 negative and PC1 positive scores included Th17 cells, eosinophils, and NK CD56 bright cells. This result means that the Th1 cell and Th2 cell signatures are almost independent. The PCA results with “Lung3” were similar to TCGA (Fig 2B). Thus, signatures of Th1, Th2, and cytotoxic T cells were selected as target variables and grouped by their mean cut-off values (high vs. low). The cut-off value was set as the average of the value of immune cell signatures. The training set and the test set applied their own average cut-off value instead of the same cut-off value, because the distributions of Th1, Th2, and cytotoxic T cells signature of the test sets (TCGA cohort) shift further to the right than the training set (“Lung3” cohort) (Fig 3).

[Figure omitted. See PDF.]

Fig 2.

Principal component analysis of immune cell signatures with 1016 TCGA (A) and “Lung3” (B) non-small cell carcinomas.

https://doi.org/10.1371/journal.pone.0231227.g002

[Figure omitted. See PDF.]

Fig 3. Distributions of type 1 helper T-cell, type 2 helper T-cell, and cytotoxic T cell signatures of the test set (TCGA cohort) and training set (“Lung3” cohort).

https://doi.org/10.1371/journal.pone.0231227.g003

Modeling for the prediction of tumor immune microenvironment

The Th1, Th2, and cytotoxic T cell groups were predicted using machine learning algorithms, which included penalized logistic regression, penalized discriminant analysis, sparse discriminant analysis, linear discriminant analysis, naive Bayes, classification and regression tree (CART), bagged CART, and random forest. R (version 3.6.0) and the caret R package were used for the prediction modeling [23]. Data pre-processing methods included removing zero variance variables, NA imputation (KNN), resolving skewness (Yeo-Johnson), centering, and scaling (Z-score normalization). Data preprocessing was applied to all predictor variables with the recipe (v.0.1.9) rpackage. The detailed preprocessing parameters set with the default value. Resampling was done using the 0.632 bootstrap estimator.

Results

Demographic information and tumor characteristics for both the training and test sets are presented in Table 1. The training dataset comprised 89 NSCLCs and the test set included 60 cases from all 1016 lung squamous cell carcinoma and adenocarcinomas whose CT images are available. For radiomic feature extraction, we evaluated a total of 239 CT radiomic features, which were divided into seven groups as follows: three physical features, 59 histogram-based features, 10 shape features, 95 local features, 63 filter-based features (LoG filter), three fractal model-based features, and six sigmoid features. Details are provided in S1 Table.

[Figure omitted. See PDF.]

Table 1. Demographic information of the “Lung3” and TCGA dataset.

https://doi.org/10.1371/journal.pone.0231227.t001

Training set prediction

Classification trees and rule-based models and recursive feature elimination methods performed better than penalized linear classification models. Random forest, linear discriminant analysis, and penalized logistic regression models performed well in the training set of Th2 cell prediction (area under the curve, AUC = 0.795, 0.772, and 0.754, respectively) (Table 2 and Fig 4). Random forest, penalized discriminant analysis, and bagged CART models performed well in the training set of Th1 cells (AUC = 0.751, 0.711, and 0.741, respectively) and cytotoxic T cells (AUC = 0.681, 0.674, and 0.647, respectively) predictions (Tables 3 and 4).

[Figure omitted. See PDF.]

Fig 4. Receiver operating characteristic curve of top 3 models on the training set of type 2 helper T cells.

https://doi.org/10.1371/journal.pone.0231227.g004

[Figure omitted. See PDF.]

Table 2. Type 2 helper T-cell signature training set prediction.

https://doi.org/10.1371/journal.pone.0231227.t002

[Figure omitted. See PDF.]

Table 3. Type 1 helper T-cell signature training set prediction.

https://doi.org/10.1371/journal.pone.0231227.t003

[Figure omitted. See PDF.]

Table 4. Cytotoxic T cell signature training set prediction.

https://doi.org/10.1371/journal.pone.0231227.t004

Test set prediction

We selected a linear discriminant analysis model for Th2 cell prediction and penalized discriminant analysis model for Th1 cell prediction because both models are simpler than bagged CART or random forest. Although the tree based ensemble methods including random forest and bagged CART showed best performance, but the performance was not significantly different between the ensemble methods and linear models including linear discriminant analysis and penalized discriminant analysis models. The prediction performance on the test set of Th2 cells was good (AUC = 0.684) (Table 5 and Fig 5) but not on Th1 cells (AUC = 0.536–0.564) and cytotoxic T cells (AUC = 0.533–0.612) (S2 and S3 Tables). Predictors of linear discriminant analysis model were skewness (total and outer pixels), kurtosis, variance (subsampled from delta [subtraction inner pixels from outer pixels]), and informational measure of correlation (IMC) (Table 6).

[Figure omitted. See PDF.]

Fig 5. Receiver operating characteristic curve of top 3 models on the test set of type 2 helper T cells.

https://doi.org/10.1371/journal.pone.0231227.g005

[Figure omitted. See PDF.]

Table 5. Performance of prediction on test set of type 2 helper T-cells.

https://doi.org/10.1371/journal.pone.0231227.t005

[Figure omitted. See PDF.]

Table 6. Linear discriminant model predictor variables for type 2 helper T-cells.

https://doi.org/10.1371/journal.pone.0231227.t006

Discussion

Several previous reports suggest that malignant phenotypes and efficacy of immunotherapy in cancers are defined not only by the intrinsic activities of cancer cells but also by components in the tumor microenvironment, especially tumor-infiltrating immune cells [24–30]. Despite the volume of previous research in this setting, the complex interactions between tumors and their microenvironment remain poorly understood. In addition, these studies have evaluated tumor-infiltrating immune cells using immunohistochemistry-based analyses alone. Most histologic approaches only involve small biopsies or surgical specimens and are therefore limited due to the heterogeneity of tumors. Thus, there is no valid noninvasive method that describes the whole tumor microenvironment, raising a compelling need for the development of a “translator” to predict the tumor microenvironment. In this situation, radiomics could offer an alternative. However, there are very few studies that have related radiomic analysis to the tumor microenvironment [31, 32].

In this study, based on the deconvolution of gene expression data from 89 patients in the “Lung3” cohort using ssGSEA and immune microenvironment profiling with the use of 28 immune cell signatures from gene expression, we found associations with certain radiomic features for subsets of Th1, Th2, and cytotoxic T cells in NSCLCs. Through a prediction performance test using a TCGA test set, we also demonstrated that the radiomic prediction for Th2 cell signatures of NSCLCs was feasible (AUC = 0.684), even though the performances of radiomics on the test set of Th1 and cytotoxic T cells were not accurate enough to be predictable. Our results indicate that skewness (total and outer pixels), kurtosis, variance (subsampled from delta), and IMC correlated with Th2 cells. Our study is in line with but distinguishable from previous reports [31–33]. The other reports showed an association between radiomics and certain immune cell signatures such as CD8 cell, CD3 cell, or tumor PD-L1 expression, and demonstrated model creation and validation in independent cohorts. And Grossmann at al. [33] used the same data set (“Lung3”) for validation cohort, and found a relationship between radiomic features, immune response, inflammation, and survival, which was further validated by immunohistochemical staining. As in our study, they showed that imaging can provide a promising way to predict the immune phenotype of tumors and to infer clinical outcomes for patients with cancer who had been treated with immunotherapy. Compared with these studies, our study had distinguishable points in that we used more integral radiomic features including filter-based and fractal model-based features reflecting the characteristics of tumor margins, as well as first and second-order or textural features, and estimated the comprehensive immune microenvironment with 28 immune cell signatures. In addition, we homogenized our data by using only NSCLC cases whereas previous two studies dealt with heterogeneous cohorts of different tumor histopathology and location [31, 32]. Even though our results do not provide good generalization across tumor types and locations, they predict the immune microenvironment from radiomic information concentrated on NSCLC, in particular, and therefore we believe our results are more specific.

Most past studies dealing with antitumor immunity have focused on CD8+ cytotoxic T cells, as their cytotoxic activity was deemed to be critical for tumor eradication. Recent work indicates that specific subsets of CD4+ T cells (namely, Th1 cells), B cells, macrophages as well as dendritic cells provide an important contribution to antitumor immune responses [28, 29]. These cytotoxic T cells and Th1 cells are the main cell signatures of the type 1 immune response [9]. In addition, other studies have demonstrated that the levels of Th2 and regulatory T (Treg) cells, which are important elements of the type 2 immune response, are associated with negative clinical outcomes [30]. In contrary to this traditional scheme of pro- and antitumor microenvironment, a recent study demonstrated that the type 2 immune response characterized by Th2 cells may inhibit tumor formation in melanoma [10]. In line with this research, the Th1, Th2, and the cytotoxic T cell signature prediction test and the demonstration that the radiomic prediction of the Th2 cell signature is feasible in our study are clinically meaningful because these three cells are subsets of influential immune cells in tumor immune responses as mentioned above.

We found that skewness (total and outer pixels), kurtosis, variance (subsampled from delta), and IMC were correlated with Th2 cells. Among these features, outer pixel skewness refers to skewness of the peripheral portion of the tumor. More interestingly, given that this radiomic feature reflects the peritumoral radiologic phenotype, this finding may indicate that the immune microenvironment of NSCLCs is closely related with certain interactions between the tumors and surrounding normal lung tissues. In other words, these edges or peritumoral radiologic findings offer information that can help predict tumor immune activity in the development of immunotherapy strategies. There has been no previous research regarding the relationship between the five radiomic features, including skewness of outer pixels and subsets of immune cells, therefore further verification of this finding in larger cohorts is needed.

Despite the advantages of utilizing a different cohort for external validation, this analysis has several limitations. First, the training data are retrospective and limited to 89 patients from a single institution. This deficiency could be addressed in future work by using a larger patient cohort. Second, CT scans we obtained for radiomic feature selection were partly contrast-enhanced and partly non-contrast and there are differences in image resolution and even may be in acquisition protocol between training and test data (Scan protocol including tube voltage and tube current was not provided). Contrast agent may obscure the radiomic textural features and variable acquisition protocol might lead to be unstable in features because textural- and intensity-based features may be affected by intensity and scanner variability. However, it is rather difficult to have the same acquisition protocol between multi-center data in practice. Testing using independent data with some differences in acquisition reflects the real clinical practice and the performance measured under such circumstance is a good measure of generalization of the proposed method. The data are obtained from the public database and many studies using the data still used the data as a whole despite the heterogeneity [13, 34, 35]. Nevertheless, further verification for reproducibility of radiomic features according to acquisition protocol is needed. Lastly, this study includes multiple histologic types in both the “Lung3” and TCGA cohorts. The histologic type could serve as a confounding factor when predicting tumor immune microenvironment status using radiomic features. However, we believe our findings and the comprehensive radiomics approach described herein are meaningful in terms of building baseline research data for the next relevant study. In addition, this study shows the usefulness of radiomics in the tumor microenvironment domain, especially for NSCLCs. Above all, we attempted to perform external validation and our results were favorable for Th2 cell signature prediction, which is clinically important.

Conclusions

In conclusion, radiomic prediction for the immune microenvironment of NSCLC was feasible and found a relationship between radiomic features and Th2 cell signatures. This study demonstrated the potential of radiomic features as noninvasive biomarkers to capture the tumor microenvironment properties and to predict the tumor-suppressive or tumor-permissive status of NSCLC. The results of this study may help define categories of tumor immune activity for patients with NSCLC and develop immunotherapy strategies.

Supporting information

[Figure omitted. See PDF.]

S1 File. Details about the different groups of radiomic features assessed in our study.

https://doi.org/10.1371/journal.pone.0231227.s001

(DOCX)

S1 Table. Radiomic features used to predict the tumor microenvironment.

https://doi.org/10.1371/journal.pone.0231227.s002

(DOCX)

S2 Table. Performance of prediction on the test set of type 1 helper T cells.

https://doi.org/10.1371/journal.pone.0231227.s003

(DOCX)

S3 Table. Performance of prediction on the test set of cytotoxic T cells.

https://doi.org/10.1371/journal.pone.0231227.s004

(DOCX)

Citation: Yoon HJ, Kang J, Park H, Sohn I, Lee S-H, Lee HY (2020) Deciphering the tumor microenvironment through radiomics in non-small cell lung cancer: Correlation with immune profiles. PLoS ONE 15(4): e0231227. https://doi.org/10.1371/journal.pone.0231227

References

1. Borghaei H, Paz-Ares L, Horn L, Spigel DR, Steins M, Ready NE, et al. Nivolumab versus docetaxel in advanced nonsquamous non-small-cell lung cancer. The New England journal of medicine. 2015;373(17):1627–39. Epub 2015/09/29. pmid:26412456; PubMed Central PMCID: PMC5705936.

2. Brahmer J, Reckamp KL, Baas P, Crino L, Eberhardt WE, Poddubskaya E, et al. Nivolumab versus docetaxel in advanced squamous-cell non-small-cell lung cancer. The New England journal of medicine. 2015;373(2):123–35. Epub 2015/06/02. pmid:26028407; PubMed Central PMCID: PMC4681400.

3. Herbst RS, Baas P, Kim DW, Felip E, Perez-Gracia JL, Han JY, et al. Pembrolizumab versus docetaxel for previously treated, PD-L1-positive, advanced non-small-cell lung cancer (KEYNOTE-010): a randomised controlled trial. Lancet (London, England). 2016;387(10027):1540–50. Epub 2015/12/30. pmid:26712084.

4. Rizvi NA, Hellmann MD, Snyder A, Kvistborg P, Makarov V, Havel JJ, et al. Cancer immunology. Mutational landscape determines sensitivity to PD-1 blockade in non-small cell lung cancer. Science (New York, NY). 2015;348(6230):124–8. Epub 2015/03/15. pmid:25765070; PubMed Central PMCID: PMC4993154.

5. Gubens MA, Davies M. NCCN Guidelines Updates: New Immunotherapy Strategies for Improving Outcomes in Non-Small Cell Lung Cancer. Journal of the National Comprehensive Cancer Network: JNCCN. 2019;17(5.5):574–8. Epub 2019/05/23. pmid:31117034.

6. Herbst RS, Soria JC, Kowanetz M, Fine GD, Hamid O, Gordon MS, et al. Predictive correlates of response to the anti-PD-L1 antibody MPDL3280A in cancer patients. Nature. 2014;515(7528):563–7. Epub 2014/11/28. pmid:25428504; PubMed Central PMCID: PMC4836193.

7. Shukuya T, Carbone DP. Predictive markers for the efficacy of anti-PD-1/PD-L1 antibodies in lung cancer. Journal of thoracic oncology: official publication of the International Association for the Study of Lung Cancer. 2016;11(7):976–88. Epub 2016/03/06. pmid:26944305.

8. DeNardo DG, Barreto JB, Andreu P, Vasquez L, Tawfik D, Kolhatkar N, et al. CD4(+) T cells regulate pulmonary metastasis of mammary carcinomas by enhancing protumor properties of macrophages. Cancer cell. 2009;16(2):91–102. Epub 2009/08/04. pmid:19647220; PubMed Central PMCID: PMC2778576.

9. Mellman I, Coukos G, Dranoff G. Cancer immunotherapy comes of age. Nature. 2011;480(7378):480–9. Epub 2011/12/24. pmid:22193102; PubMed Central PMCID: PMC3967235.

10. Wolf MT, Ganguly S, Wang TL, Anderson CW, Sadtler K, Narain R, et al. A biologic scaffold-associated type 2 immune microenvironment inhibits tumor formation and synergizes with checkpoint immunotherapy. Science translational medicine. 2019;11(477). Epub 2019/02/01. pmid:30700576.

11. Sadtler K, Estrellas K, Allen BW, Wolf MT, Fan H, Tam AJ, et al. Developing a pro-regenerative biomaterial scaffold microenvironment requires T helper 2 cells. Science (New York, NY). 2016;352(6283):366–70. Epub 2016/04/16. pmid:27081073; PubMed Central PMCID: PMC4866509.

12. Lee G, Lee HY, Ko ES, Jeong WK. Radiomics and imaging genomics in precision medicine. Precis Future Med. 2017;1(1):10–31.

13. Aerts HJWL, Velazquez ER, Leijenaar RTH, Parmar C, Grossmann P, Carvalho S, et al. Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach. Nature communications. 2014;5:4006. pmid:24892406

14. Limkin EJ, Sun R, Dercle L, Zacharaki EI, Robert C, Reuze S, et al. Promises and challenges for the implementation of computational medical imaging (radiomics) in oncology. Annals of oncology: official journal of the European Society for Medical Oncology. 2017;28(6):1191–206. Epub 2017/02/09. pmid:28168275.

15. Network TCGAR. Comprehensive genomic characterization of squamous cell lung cancers. Nature. 2012;489(7417):519–25. pmid:22960745

16. Network TCGAR. Comprehensive molecular profiling of lung adenocarcinoma. Nature. 2014;511(7511):543–50. pmid:25079552

17. Prior FW, Clark K, Commean P, Freymann J, Jaffe C, Kirby J, et al. TCIA: An information resource to enable open science. Conference proceedings: Annual International Conference of the IEEE Engineering in Medicine and Biology Society IEEE Engineering in Medicine and Biology Society Annual Conference. 2013;2013:1282–5. Epub 2013/10/11. pmid:24109929; PubMed Central PMCID: PMC4257783.

18. van Griethuysen JJM, Fedorov A, Parmar C, Hosny A, Aucoin N, Narayan V, et al. Computational Radiomics System to Decode the Radiographic Phenotype. Cancer research. 2017;77(21):e104–e7. Epub 2017/11/03. pmid:29092951; PubMed Central PMCID: PMC5672828.

19. Cho HH, Lee G, Lee HY, Park H. Marginal radiomics features as imaging biomarkers for pathological invasion in lung adenocarcinoma. European radiology. 2020. Epub 2020/01/23. pmid:31965255.

20. Song SH, Park H, Lee G, Lee HY, Sohn I, Kim HS, et al. Imaging Phenotyping Using Radiomics to Predict Micropapillary Pattern within Lung Adenocarcinoma. Journal of thoracic oncology: official publication of the International Association for the Study of Lung Cancer. 2017;12(4):624–32. Epub 2016/12/08. pmid:27923715.

21. Hanzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics. 2013;14:7. Epub 2013/01/18. pmid:23323831; PubMed Central PMCID: PMC3618321.

22. Bindea G, Mlecnik B, Tosolini M, Kirilovsky A, Waldner M, Obenauf AC, et al. Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. Immunity. 2013;39(4):782–95. Epub 2013/10/22. pmid:24138885.

23. Kuhn M. Building Predictive Models in R Using the caret Package. Journal of Statistical Software. 2008;28(1):1–26.

24. Candido J, Hagemann T. Cancer-related inflammation. Journal of clinical immunology. 2013;33 Suppl 1:S79–84. Epub 2012/12/12. pmid:23225204.

25. Mantovani A, Allavena P, Sica A, Balkwill F. Cancer-related inflammation. Nature. 2008;454(7203):436–44. Epub 2008/07/25. pmid:18650914.

26. DeNardo DG, Andreu P, Coussens LM. Interactions between lymphocytes and myeloid cells regulate pro- versus anti-tumor immunity. Cancer metastasis reviews. 2010;29(2):309–16. Epub 2010/04/21. pmid:20405169 PubMed Central PMCID: PMC2865635.

27. Jochems C, Schlom J. Tumor-infiltrating immune cells and prognosis: the potential link between conventional cancer therapy and immunity. Experimental biology and medicine (Maywood, NJ). 2011;236(5):567–79. Epub 2011/04/14. pmid:21486861; PubMed Central PMCID: PMC3229261.

28. Gentles AJ, Newman AM, Liu CL, Bratman SV, Feng W, Kim D, et al. The prognostic landscape of genes and infiltrating immune cells across human cancers. Nature medicine. 2015;21(8):938–45. Epub 2015/07/21. pmid:26193342; PubMed Central PMCID: PMC4852857.

29. Gu-Trantien C, Willard-Gallo K. Tumor-infiltrating follicular helper T cells: The new kids on the block. Oncoimmunology. 2013;2(10):e26066. Epub 2013/11/19. pmid:24244900; PubMed Central PMCID: PMC3825729.

30. Şenbabaoğlu Y, Winer AG, Gejman RS, Liu M, Luna A, Ostrovnaya I, et al. The landscape of T cell infiltration in human cancer and its association with antigen presenting gene expression. bioRxiv. 2015:025908.

31. Sun R, Limkin EJ, Vakalopoulou M, Dercle L, Champiat S, Han SR, et al. A radiomics approach to assess tumour-infiltrating CD8 cells and response to anti-PD-1 or anti-PD-L1 immunotherapy: an imaging biomarker, retrospective multicohort study. The Lancet Oncology. 2018;19(9):1180–91. Epub 2018/08/19. pmid:30120041.

32. Tang C, Hobbs B, Amer A, Li X, Behrens C, Canales JR, et al. Development of an immune-pathology informed radiomics model for non-small cell lung cancer. Scientific reports. 2018;8(1):1922. Epub 2018/02/02. pmid:29386574; PubMed Central PMCID: PMC5792427.

33. Grossmann P, Stringfield O, El-Hachem N, Bui MM, Rios Velazquez E, Parmar C, et al. Defining the biological basis of radiomic phenotypes in lung cancer. eLife. 2017;6. Epub 2017/07/22. pmid:28731408; PubMed Central PMCID: PMC5590809.

34. Hosny A, Parmar C, Coroller TP, Grossmann P, Zeleznik R, Kumar A, et al. Deep learning for lung cancer prognostication: A retrospective multi-cohort radiomics study. PLoS medicine. 2018;15(11):e1002711. Epub 2018/12/01. pmid:30500819; PubMed Central PMCID: PMC6269088 following competing interests: RHM reports personal fees (Scientific Advisory Board) from AstraZeneca, personal fees from New RT (lecture honorarium), outside the submitted work. RJG is a consultant, shareholder and grantee of HealthMyne, Inc., and a consultant, shareholder and grantee of Helix Biopharma. Neither Helix nor HealthMyne have a financial stake in the results of this study. All other authors declare no conflict of interest.

35. Chaddad A, Desrosiers C, Toews M, Abdulkarim B. Predicting survival time of lung cancer patients using radiomic analysis. Oncotarget. 2017;8(61):104393–407. Epub 2017/12/22. pmid:29262648; PubMed Central PMCID: PMC5732814.