- 1Department of Ultrasound, The Affiliated Drum Tower Hospital of Nanjing University Medical School, Nanjing, China
- 2Department of General Surgery, Nanfang Hospital, Southern Medical University, Guangdong Provincial Engineering Technology Research Center of Minimally Invasive Surgery, Guangzhou, China
- 3School of Biomedical Engineering, Southern Medical University, Guangzhou, China
- 4Department of Radiology, Southern Medical University, Guangzhou, China
Objective: The aim of this study is to identify prognostic imaging biomarkers and create a radiogenomics nomogram to predict overall survival (OS) in gastric cancer (GC).
Material: RNA sequencing data from 407 patients with GC and contrast-enhanced computed tomography (CECT) imaging data from 46 patients obtained from The Cancer Genome Atlas (TCGA) and The Cancer Imaging Archive (TCIA) were utilized to identify radiogenomics biomarkers. A total of 392 patients with CECT images from the Nanfang Hospital database were obtained to create and validate a radiogenomics nomogram based on the biomarkers.
Methods: The prognostic imaging features that correlated with the prognostic gene modules (selected by weighted gene coexpression network analysis) were identified as imaging biomarkers. A nomogram that integrated the radiomics score and clinicopathological factors was created and validated in the Nanfang Hospital database. Nomogram discrimination, calibration, and clinical usefulness were evaluated.
Results: Three prognostic imaging biomarkers were identified and had a strong correlation with four prognostic gene modules (P < 0.05, FDR < 0.05). The radiogenomics nomogram (AUC = 0.838) resulted in better performance of the survival prediction than that of the TNM staging system (AUC = 0.765, P = 0.011; Delong et al.). In addition, the radiogenomics nomogram exhibited good discrimination, calibration, and clinical usefulness in both the training and validation cohorts.
Conclusions: The novel prognostic radiogenomics nomogram that was constructed achieved excellent correlation with prognosis in both the training and validation cohort of Nanfang Hospital patients with GC. It is anticipated that this work may assist in clinical preferential treatment decisions and promote the process of precision theranostics in the future.
Introduction
Although the incidence of gastric cancer (GC) has decreased over the last 3 decades, it remains the second leading cause of cancer-related death worldwide and the most prevalent cancer in eastern Asia, whose burden is still substantial (1). Patients experience a high rate of tumor recurrence (20%–40%) although following treatment with traditional standard therapies—surgical resection plus adjuvant chemotherapy or radiochemotherapy (2, 3). The American Joint Committee on Cancer (AJCC) tumor-node-metastasis (TNM) staging system is currently the most commonly used tumor staging system worldwide and is considered to be the most valuable method for assessing the prognosis of malignant tumors (4). However, large variations in clinical outcomes have been shown in patients with the same stage and similar treatment regimens (5, 6). These findings suggest that the present GC staging system provides inadequate prognostic information and does not reflect the biological heterogeneity of GC (7). Therefore, identifying additional effective prognostic model that owns incremental value to the TNM staging system is necessary to achieve individualized medical treatments.
A radiogenomic analysis aimed at investigating molecular biomarkers can identify imaging traits that correspond to different molecular phenotypes with clinical and biologic relevance. It represents the evolution of the radiology–pathology correlation from the tissue level to the subcellular level and can compensate for the deficiency that the biological interpretations of imaging traits are lacking. Radiogenomics analysis can hold the promise to meet the clinical demand noninvasively and comprehensively, compared to biopsy with limited specimens, because the extensive imaging features for each patient can comprehensively describe the tumor phenotype characteristics. This non-invasive method can be performed repeatedly and is therefore eminently suitable for treatment follow-up.
Recent advances in the radiogenomics of other cancers, such as hepatocellular carcinoma (8), non–small cell lung cancer (9), glioblastoma multiforme (10), breast cancer (11), and liver cancer (12), have confirmed the potential synergy of integrating imaging and genomic data. However, thus far, few studies have investigated the radiogenomics of GC.
A nomogram is based on multivariate regression analysis and includes important influencing factors related to tumor prognosis. The nomogram has become the focus of interest in cancer research in recent years and is considered a useful tool for quantifying risk (13, 14). Therefore, this study aimed to develop a predictive model of the overall survival (OS) nomogram based on the radiogenomic features combined with the clinicopathological characteristics to predict prognosis and validate its incremental prognostic value to the traditional TNM staging system.
Materials and Methods
Data Source
We studied 407 patients with GC whose data (including RNA sequencing data and clinical information) was obtained from The Cancer Genome Atlas (TCGA) (15)—a landmark cancer genomics program that contains a large amount of genomic and clinical datasets. RNA sequencing data were normalized using the fragments per kilobase transcriptome per million reads method. Among the 407 patients with GC, 46 patients with contrast-enhanced computed tomography (CECT) imaging data were obtained from The Cancer Imaging Archive (TCIA) (16)—a service that hosts a large archive of medical images of cancer accessible for public download. Most patients underwent imaging on a multi-detector row CT scanners: 4-, 16-, or 64-slice CT scanners with slice thickness of 2.5~5 mm at 120 kVp and 200~500 mA.
In addition, because the number of the patients in TCGA database was limited (there were only 46 patients who owned imaging data), we obtained external 424 patients with GC with imaging data in Nanfang Hospital to develop a prognostic nomogram model. Thirty-two patients were excluded for the invisible lesion or insufficient image quality, and then 392 patients with GC were included in further research. All of these patients underwent partial or total radical gastrectomy and did not receive other therapy preoperatively. All of the contrast-enhanced abdominal CT images were acquired within 30 days before surgical resection; the detailed inclusion and exclusion criteria were described in our previous study (17). Patient characteristics are detailed in Table 1.
Imaging Feature Extraction
To avoid artificial and subjective deviation, two professional abdominal radiologists cooperated to delineate each tumor using the Medical Imaging Interaction Toolkit (MITK) (version 3.6.0 Apr. 1, 2017) without referring to any clinical information. Four samples are shown in Supplementary Figure S1. Forty-four imaging features were extracted, including four first-order statistics and 40 texture features, from the delineated tumor outlines for each person using MATLAB automatically. Four different matrices, including the gray-level cooccurrence matrix (GLCM) (18), the neighborhood gray-tone difference matrix (NGTDM) (19), the gray-level size zone matrix (GLSZM) (20), and the gray-level run-length matrix (GLRLM) (21), were utilized to evaluate the 40 texture features of the regions of interest (ROIs). A detailed description of the four matrices and the features is shown in Supplementary Table S1.
Intra- and interclass correlation coefficients (ICCs) were used to assess the intra- and interobserver reproducibility of radiomics feature extraction. Thirty images were chosen randomly for ROI segmentation and feature extraction. The ROI segmentation was performed by two radiologists with, respectively, 12 (reader 1) and 8 years (reader 2) of experience in gastric CT interpretation. Reader 1 then repeated the same procedure 1 week later. An ICC greater than 0.75 indicates good agreement of the feature extraction. To ensure the independence of imaging features, we removed the features with high correlation. We calculated the mean absolute correlation using the “findCorrelation” function in the caret package (version 6.0-84) in R software. When the pairwise correlation was greater than 0.8, the feature was identified as the one with redundancy and will be eliminated.
Gene Module Clustering
Compared with traditional analysis methods, which focus on a single gene or a few genes, we used weighted gene coexpression network analysis (WGCNA) (22), a systematic biological method to describe the pattern of gene association between different samples and obtain highly synergistic gene sets. The distinct advantage is that WGCNA converts an enormous amount of highly correlated genes into modules defined as clusters of densely interconnected genes to eliminate the problem of multiple hypothesis testing corrections.
WGCNA methodology can be used to construct a gene coexpression network, identify modules using hierarchical clustering and the dynamic tree cut method, study module relationships using eigengene networks, and obtain hub genes—the key regulators of each module. Because of the early advantages, the application of the WGCNA method has therefore been extended to variable types of high-throughput datasets, such as proteomic and metabolomic datasets, in recent years (23–25). To ensure the reliability of the prognostic gene modules, all the 407 patients with gene expression data were included.
Survival Analysis and Functional Enrichment of Gene Modules
The survival analysis of a time-to-event outcome is a widely applied method to predict survival from a set of patient predictors or covariates. Univariate Cox proportional hazards regression (26) was used to identify the gene modules with significant prognostic ability in OS, defined as the time from pathology diagnosis to the date of death or the last follow-up. For the module eigengene of gene modules that were continuous, the contributions to OS were investigated by Kaplan–Meier survival analysis (log-rank test).
We annotated the biological functions of the survival-related modules of GC using Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, allowing the identification of common biologic pathways from a public database. GO and KEGG analyses were performed using the clusterProfiler package (27, 28) in R. In addition, we calculated the false discovery rate (FDR) using the Benjamini–Hochberg method (29) to correct the multiple hypothesis testing of the survival analysis.
Development and Ability Evaluation of Imaging Features
We utilized Spearman rank correlations to build a radiogenomic map combining selected imaging features and prognostic gene modules. The FDR was used again to correct for multiple testing. The imaging features that showed high correlation with the prognostic gene modules (P < 0.05, FDR < 0.05) were further analyzed by the univariate Cox proportional hazards model to assess the prognostic significance. Afterward, a bootstrap approach (30) was used to assess whether the prognostic signature had significant power compared with random chance (C-index = 0.5). For 100 times, we calculated the C-index of the imaging signature based on 35 randomly selected samples with correct outcome data and on 35 randomly chosen samples with random outcome data. The Wilcoxon test was applied to assess the differences between the two distributions (imaging signature and random chances).The prognostic signature with significant power (P < 0.05) was identified as the imaging biomarker.
Generation of the Radiomics Score
On the basis of the prognostic imaging biomarkers identified from the TCGA database, the radiomics score (RIS) was constructed using a multivariate Cox proportional hazards regression model. The equation of RIS is detailed in the Supplementary Material. Because the median RIS is convenient for clinical application, it was used as a cutoff value to divide the patients into low‐risk and high‐risk groups. The log-rank test (Mantel–Cox test) was used to compare the difference in the survival status between two groups. We calculated RIS both in the TCGA database and Nanfang Hospital database, and log-rank test was used to identify their prognostic value. To discuss the prognostic significance of the imaging biomarkers and RIS, we tried to explain it from the perspective of immunity. We utilized the CIBERSORT algorithm (an algorithm that calculates the ratios of infiltrated immune cells) to investigate the immune infiltration level in patients with GC.
Construction and Validation of the Nomogram
To compensate for the limitation of the number of patients in TCGA database, we utilized the database from Nanfang Hospital to construct a prognostic nomogram model for its large samples. We divided patients from Nanfang Hospital into training and validation cohorts in a ratio of 1:1 using the stratified randomization method. The training cohort was used to construct the nomogram that integrated both the RIS and significant clinicopathological characteristics, and the validation cohort was used to identify the predictive accuracy of the nomogram. The nomogram was constructed using “rms” package (31) of R software. With the application of the bootstrap method (1,000 replicates), a calibration curve was used to visualize the deviation of predicted probabilities from what actually happened. Moreover, decision curve analysis (32) was employed to quantify the clinical utility and compare the nomogram that we created with the TNM staging system. Finally, the receiver operating characteristic curve (ROC) analysis (33) was performed to measure the predictive performance of the nomogram score and the TNM stage. The Delong test was used to assess the differences between the two models.
Statistical Analysis
Only two-sided P-value < 0.05 was considered to have necessary to statistic. The statistical analysis was implemented on R software (version 4.0.4). Several packages including ggplot2, rms, survival, survivalROC, survminer, nomogramFormula, and rmda were used to perform other statistical calculations and graphical work.
Results
Imaging Feature Extraction
We extracted 44 imaging features from each person, and 16 imaging features were removed because their pairwise correlation was greater than 0.8. The correlation matrix is shown in Supplementary Figure S2. The intraobserver ICCs ranged from 0.812 to 0.978, and the interobserver ICCs ranged from 0.786 to 0.912, indicating a favorable intra- and interobserver feature extraction reproducibility.
Gene Module Clustering and Function Annotation
We grouped transcripts with correlated expression levels into gene coexpression modules using the WGCNA approach. The cluster dendrogram and gene counts for each module are shown in Supplementary Figure S3A. Analysis of network topology for various soft-thresholding powers is shown in Figures S3B, C. The soft thresholding power β was set at 9 because the scale independence reached 0.8 (Figure S3B) and had a relatively high-average connectivity (Figure S3C). At last, we identified 33 gene modules among 58,428 genes. Across the 33 modules, four prognostic modules (MEgreen, MEwhite, MEdarkturquoise, and MEroyalblue) were identified using univariate Cox proportional hazards regression (Table 2), and the Kaplan–Meier plots of the four prognostic modules are shown in Figure 1. Afterward, four modules were tested for the enrichment of specific biological functions in GO terms and KEGG pathways. The results are shown in Supplementary Table S2. For example, the gene module MEdarkturquoise, enriched in protein kinase activator activity, was strongly correlated with a poor prognosis (P < 0.05). The gene module MEwhite, enriched in NADH dehydrogenase activity, was strongly correlated with positive survival.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/882786/fonc-12-882786-HTML-r1/image_m/fonc-12-882786-g001.jpg)
Figure 1 Kaplan–Meier plot with the univariate survival analysis of four gene modules. (A) The Kaplan–Meier plot shows that patients with a higher “MEdarkturquoise” value have shorter OS (blue lines), and patients with a lower “MEdarkturquoise” value have longer OS (red lines). (B) The Kaplan–Meier plot shows that patients with a higher “MEgreen” value have shorter OS (blue lines), and patients with a lower “MEgreen” value have longer OS (red lines). (C) The Kaplan–Meier plot shows that patients with a higher “MEroyalblue” value have shorter OS (blue lines), and patients with a lower “MEroyalblue” value have longer OS (red lines). (D) The Kaplan–Meier plot shows that patients with a higher “MEwhite” value have longer OS (blue lines), and patients with a lower “MEwhite” value have shorter OS (red lines).
Radiogenomic Correlation
We created the radiogenomic map that combined the remaining imaging features and four prognostic gene modules using Spearman rank correlations (Figure 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/882786/fonc-12-882786-HTML-r1/image_m/fonc-12-882786-g002.jpg)
Figure 2 Radiogenomic map revealing 19 statistically significant associations between 15 CT semantic features and four prognostic gene modules in GC. Significant pairwise correlations (P < 0.05) are indicated with star symbols. The corresponding correlation coefficients are displayed in a heatmap; the red color indicates positive correlations, and the blue color indicates negative correlations.
One first-order statistic (volume.eccentricity) and 14 texture features (GLCM.Energy, GLCM.Contrast, GLCM.Entropy, GLCM.Variance, GLRLM.GLN, GLRLM.RP, GLRLM.SRLGE, GLRLM.LRLGE, GLSZM.GLN, GLSZM.ZSN, GLSZM.ZP, GLSZM.HGZE, GLSZM.LZLGE, and NGTDM.Complexity) were identified to have a strong correlation with the gene modules (P < 0.05, FDR < 0.05). There were 19 pairs of significant correlations among the features and gene modules. For example, the gene module MEdarkturquoise was negatively correlated with GLCM.Energy (P = 0.0379), GLRLM.GLN (P = 0.0161), GLSZM.ZSN (P = 0.0181), and GLSZM.ZP (P=0.0224) and positively correlated with GLCM.Contrast (P = 0.0079) and GLRLM.SRLGE (P = 0.0115).
Identification of Imaging Biomarkers and Radiomics Score
Univariate Cox proportional hazards regression was employed to assess the survival relationship of the 15 imaging features with survival outcomes. Three texture features, namely, NGTDM.Complexity, GLRLM.LRLGE, and GLSZM.ZP (shown in Table 3), were significantly correlated with OS. By using the bootstrap method, the Wilcoxon test showed significant differences between the imaging features and random chances. Thus, the three imaging features were identified as the imaging biomarkers of GC. In the TCGA database, the RIS was constructed using a multivariate Cox proportional hazards regression model based on three imaging features and their corresponding coefficients (Table S3). Then, patients were divided into a low-RIS group (n = 22) and a high-RIS group (n = 23) according to the median RIS. Compared with the low-RIS group, high-RIS group had a significantly lower OS rate [Hazard Ratio (HR) = 3.05, 95% confidence interval (CI): 1.35–6.92] (Figure 3A). To verify the ability of the RIS in predicting the OS of patients with GC, we further validated our findings in the database from Nanfang hospital, which yielded the similar results as above (HR = 3.41,95% CI: 2.41–4.82) (Figure 3B). Furthermore, multivariate Cox regression analysis was performed to determine whether the RIS was an independent prognostic factor for patients’ OS in both databases. In TCGA database, the P-value of RIS was 0.003 (HR = 1.99, 95% CI: 1.27–3.12), and in the database from Nanfang hospital, the P-value of RIS was less than 0.001(HR = 2.04, 95% CI: 1.52–2.75). Thus, RIS could be identified as an independent prognostic factor. Correlation between the imaging biomarkers/RIS and immune infiltration level is shown in the Supplementary Material.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/882786/fonc-12-882786-HTML-r1/image_m/fonc-12-882786-g003.jpg)
Figure 3 Kaplan–Meier estimates of the patients’ survival using the radiomics score. The Kaplan–Meier plots were used to visualize the patients’ survival probabilities for the low-RIS versus high-RIS group of patients based on the median radiomics score. (A) Kaplan–Meier curves for TCGA database patients (N = 45). (B) Kaplan–Meier curves for patients from Nanfang Hospital database (N = 392). The differences between the two curves were determined by the two-side log-rank test.
Building a Predictive Radiomics Nomogram for OS Prediction in GC
In the training cohort of Nanfang Hospital database, the multivariate Cox analysis was employed to identify the independent prognostic factor in the RIS and clinicopathological characteristics. RIS, tumor metastasis, pathological lymph node positive detection rate, primary tumor, and age were significantly associated with OS (Table 4), and they were integrated to construct a nomogram (Figure 4A). After adding the sum of the points of the five variables, we can obtain the patients’ OS probabilities. Patients were divided into high– and low–nomogram point groups based on the median point. The Kaplan–Meier curve showed that the OS of the high–nomogram point group was significantly poorer [ hazarad ratio (HR) = 8.58, 95% CI: 5.17–14.23, P < 0.001] (Figure 4B). The calibration curve of our nomogram is shown in Figures 4C, D, which demonstrated the accurate predictive ability of the nomogram. Moreover, the nomogram model showed a better net benefit and broader threshold probability than the AJCC TNM staging system in the decision curves (Figure 4E), and if the threshold probability is between 0 and 0.53, then using the radiogenomics nomogram to predict OS adds more benefit than AJCC stage. ROC analysis revealed that the nomogram model (AUC = 0.838) had a better prognostic value compared to the TNM stage (AUC = 0.765) (Figure 4F). Moreover, the AUC value of 1-, 3-, and 5-year survival were 0.801 (95% CI: 0.707–0.895), 0.829 (95% CI: 0.755–0.904), and 0.809 (95% CI: 0.688–0.931), respectively (Figure 4G). It indicated that the nomogram model was a powerful tool for predicting patients’ survival.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/882786/fonc-12-882786-HTML-r1/image_m/fonc-12-882786-t004.jpg)
Table 4 Multivariable Cox regression analysis of clinical pathological parameters in training cohort and the whole cohort.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/882786/fonc-12-882786-HTML-r1/image_m/fonc-12-882786-g004.jpg)
Figure 4 Development and validation of the raidiomics nomogram. (A) Nomogram constructed in conjunction with the radiomics score and clinical characterization that predict the 3- and 5-year overall survival of patients with gastric cancer. (B) Kaplan–Meier curves for patients with high and low nomogram score in the training cohort. (C, D) Plots depict the calibration of radiomics nomograms in terms of agreement between predicted and observed 3-year (C) and 5-year (D) outcomes. (E) Decision curves of the nomogram model and TNM stage for the survival predictions of patients with GC. (F) The ROC comparation between the nomogram model and the TNM stage. (G) Time‐dependent ROC analysis of the nomogram model for OS prediction in the training cohort. The area under the ROC curve was 0.803, 0.838, and 0.811 for the nomogram score at 1, 3, and 5 years, respectively.
Validation of the Radiomics Nomogram in the Validation Cohort
To further identify the predictive ability of the nomogram, we performed similar analysis in the validation cohort. The nomogram score for each patient was calculated in the same way (detailed in Table S4). Consisted with the findings mentioned above, Kaplan–Meier analysis showed that the patients in high–nomogram point group had a poor OS (HR = 4.93, 95% CI: 3.05–7.95, P < 0.001) (Figure S4A). Moreover, ROC analysis demonstrated that the nomogram (AUC = 0.791) model improved prognostic value in the validation compared with the TNM stage (AUC = 0.686) (Figure S4B). In addition, the time-dependent ROC showed that the nomgram model had a good accuracy with 0.843 (95% CI: 0.772–0.913) in 1 year, 0.772 (95% CI: 0.686–0.857) in 3 years and 0.767 (95% CI: 0.664–0.869) in 5 years (Figure S4C).
Discussion
In this study, we integrated quantitative CT imaging features with RNA sequencing data to create a radiogenomic map and identified radiogenomic biomarkers of GC. The results revealed three prognostic imaging biomarkers and 19 significant pairwise associations between the imaging features and prognostic gene modules. Then, a RIS based on the three biomarkers was developed and proved to be independently associated with OS. A nomogram integrating both the RIS and clinicopathological characteristics performed better than the traditional TNM staging system, which demonstrated well the incremental value of the radiogenomic biomarker for individualized OS estimation.
The advantages of our radiogenomics analysis are as follows. First, radiogenomics analysis can hold the promise to meet the clinical demand noninvasively and comprehensively. Although the molecular biomarkers gained from preoperative biopsy have clinical benefit for prognosis of GC, the invasive approach may increase the potential risk to the patients. Moreover, compared with biopsy with limited specimens, we extracted extensive imaging features for each patient to comprehensively describe the tumor phenotype characteristics. In addition, there are barriers of high cost and technical problems when implementing invasive biopsy (34). Second, as for traditional radiomics research studies, the imaging features lack biological interpretations. Radiogenomics, combining the imaging features and molecular signatures, can compensate for this deficiency. In our study, prognostic value for clinical imaging with relevant molecular biology information is provided. Third, in contrast to traditional gene analysis methods, which focus on a single gene or a few genes and fail to reflect the association between genes, we employed WGCNA to divide 58,428 genes of each GC sample into 33 modules according to the correlation between genes. Therefore, complicated combinations and multiple testing problems can be avoided. In addition, we studied molecular prognostic significance by measuring OS predictors from public cohorts, which showed the prognostic associations of each gene module in GC. The three imaging biomarkers and correlations between the prognostic gene modules and imaging features were discussed in the Supplementary Material.
In the nomogram, RIS, tumor metastasis, lymph node positive detection rate, primary tumor, and age were kept after multivariate Cox analysis. The AUC of this radiogenomic nomogram is 0.838 in the training cohort and 0.791 in the validation cohort, which is higher than the previously published genomic nomogram (AUC = 0.744) by Chen et al. (35) and radiomic nomogram (AUC = 0.771) by Wang et al. (36). In addition, combining radiogenomic, clinical, and pathological information together produced this nomogram with better performance than using TNM staging information alone (Figure 4F). Thus, this comprehensive and individualized risk score calculation method could be used as stratification criteria in randomized studies and clinical trials.
Our study still had some limitations. First, the sample size of the CT images corresponding to gene expression profiles in TCGA database was not large enough, which may limit the reliability (8). Therefore, further studies with larger sample sizes are needed. However, to some extent, external validation from Nanfang Hospital was used to improve the dependability in this study. Second, the image datasets downloaded from the TCIA were extremely heterogeneous in terms of scanner modalities, manufacturers, and acquisition protocols, which may have resulted in adverse effects (37). Third, because of the limited experimental conditions and lack of enough funding, the significant correlations between imaging features and gene data could not be validated in molecular biology experiment and clinical trials. However, to our knowledge, this is the first study to predict OS for patients with GC by radiogenomics, which might provide valuable references in this field.
Conclusion
In conclusion, in this study, we presented three imaging biomarkers with significant prognostic value and a nomogram prognostic model for patients with GC. The imaging biomarkers correlated with RNA expression profiles can provide possible biological interpretations and may facilitate a better understanding of the molecular characterization of GC. The proposed nomogram showed better prognostic performance than the TNM staging system. Researchers, clinicians, and patients can easily predict the survival probability using this nomogram.
Data Availability Statement
The data that support the findings of this study are available from the corresponding author, TC, upon reasonable request.
Author Contributions
Conceptualization: TC, HL, YyW, YqL, and WtK; Formal analysis: HL, YyW and LC; Funding acquisition: TC and WtK; Methodology, TC, HL, YyW, JyY, and ZC; Software: HL, YyW, DyL, YyZ, and YL; Supervision: TC, HL, JgY, and WtK; Writing—original draft: HL, YyW, WtK, and TC; Writing—review and editing: HL, TC, and WtK. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by Guangdong Provincial Key Laboratory of Precision Medicine for Gastrointestinal Cancer (2020B121201004), the Guangdong Natural Science Foundation General Project (2019A1515011520), Guangdong Natural Science Foundation Outstanding Youth Project (2021B1515020055), and Guangdong Provincial Major Talents Project (2019JC05Y361).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We are grateful to all of the patients and cooperators for participating in this retrospective study, especially for the data from TCGA and TCIA.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.882786/full#supplementary-material
References
1. Bray F, Ren JS, Masuyer E, Ferlay J. Global Estimates of Cancer Prevalence for 27 Sites in the Adult Population in 2008. Int J Cancer (2013) 132:1133–45. doi: 10.1002/ijc.27711
2. Cunningham D, Allum WH, Stenning SP, Thompson JN, Van deVelde CJH, Nicolson M, et al. Perioperative Chemotherapy Versus Surgery Alone for Resectable Gastroesophageal Cancer. N Engl J Med (2006) 355:11–20. doi: 10.1056/NEJMoa055531
3. Macdonald JS, Smalley SR, Benedetti J, Hundahl SA, Estes NC, Stemmermann GN, et al. Chemoradiotherapy After Surgery Compared With Surgery Alone for Adenocarcinoma of the Stomach or Gastroesophageal Junction. N Engl J Med (2001) 345:725–30. doi: 10.1056/NEJMoa010187
4. Yang B, Zhong J, Zhong J, Ma L, Li A, Ji H, et al. Development and Validation of a Radiomics Nomogram Based on (18)F-Fluorodeoxyglucose Positron Emission Tomography/Computed Tomography and Clinicopathological Factors to Predict the Survival Outcomes of Patients With Non-Small Cell Lung Cancer. Front Oncol (2020) 10:1042. doi: 10.3389/fonc.2020.01042
5. Bang YJ, Kim Y-W, Yang H-K, Chung HC, Park Y-K, Lee KH, et al. Adjuvant Capecitabine and Oxaliplatin for Gastric Cancer After D2 Gastrectomy (CLASSIC): A Phase 3 Open-Label, Randomised Controlled Trial. Lancet (2012) 379:315–21. doi: 10.1016/S0140-6736(11)61873-4
6. Jackson A, O'Connor JP, Parker GJ, Jayson GC. Imaging Tumor Vascular Heterogeneity and Angiogenesis Using Dynamic Contrast-Enhanced Magnetic Resonance Imaging. Clin Cancer Res (2007) 13:3449–59. doi: 10.1158/1078-0432.CCR-07-0238
7. Noh SH, Park SR, Yang H-K, Chung HC, Chung I-J, Kim S-W, et al. Adjuvant Capecitabine Plus Oxaliplatin for Gastric Cancer After D2 Gastrectomy (CLASSIC): 5-Year Follow-Up of an Open-Label, Randomised Phase 3 Trial. Lancet Oncol (2014) 15:1389–96. doi: 10.1016/S1470-2045(14)70473-5
8. Xia W, Chen Y, Zhang R, Yan Z, Zhou X, Zhang B, et al. Radiogenomics of Hepatocellular Carcinoma: Multiregion Analysis-Based Identification of Prognostic Imaging Biomarkers by Integrating Gene Data-a Preliminary Study. Phys Med Biol (2018) 63:35044. doi: 10.1088/1361-6560/aaa609
9. Zhou M, Leung A, Echegaray S, Gentles A, Shrager JB, Jensen KC, et al. Non-Small Cell Lung Cancer Radiogenomics Map Identifies Relationships Between Molecular and Imaging Phenotypes With Prognostic Implications. Radiology (2018) 286:307–15. doi: 10.1148/radiol.2017161845
10. Diehn M, Nardini C, Wang DS, McGovern S, Jayaraman M, Liang Y, et al. Identification of Noninvasive Imaging Surrogates for Brain Tumor Gene-Expression Modules. Proc Natl Acad Sci USA (2008) 105:5213–8. doi: 10.1073/pnas.0801279105
11. Yamamoto S, Han W, Kim Y, Du L, Jamshidi N, Huang D, et al. Breast Cancer: Radiogenomic Biomarker Reveals Associations Among Dynamic Contrast-Enhanced MR Imaging, Long Noncoding RNA, and Metastasis. Radiology (2015) 275:384–92. doi: 10.1148/radiol.15142698
12. Segal E, Sirlin CB, Ooi C, Adler AS, Gollub J, Chen X, et al. Decoding Global Gene Expression Programs in Liver Cancer by Noninvasive Imaging. Nat Biotechnol (2007) 25:675–80. doi: 10.1038/nbt1306
13. Han DS, Suh Y-S, Kong S-H, Lee H-J, Choi Y, Aikou S, et al. Nomogram Predicting Long-Term Survival After D2 Gastrectomy for Gastric Cancer. J Clin Oncol (2012) 30:3834–40. doi: 10.1200/JCO.2012.41.8343
14. Tang XR, Li Y-Q, Liang S-B, Jiang W, Liu F, Ge W-X, et al. Development and Validation of a Gene Expression-Based Signature to Predict Distant Metastasis in Locoregionally Advanced Nasopharyngeal Carcinoma: A Retrospective, Multicentre, Cohort Study. Lancet Oncol (2018) 19:382–93. doi: 10.1016/S1470-2045(18)30080-9
15. Tomczak K, Czerwinska P, Wiznerowicz M. The Cancer Genome Atlas (TCGA): An Immeasurable Source of Knowledge. Contemp Oncol (Pozn) (2015) 19:A68–77. doi: 10.5114/wo.2014.47136
16. Prior FW, Clark K, Commean P, Freymann J, Jaffe C, Kirby J, et al. TCIA: An Information Resource to Enable Open Science. Annu Int Conf IEEE Eng Med Biol Soc (2013) 2013:1282–5. doi: 10.1109/EMBC.2013.6609742
17. Chen T, Li X, Mao Q, Wang Y, Li H, Wang C, et al. An Artificial Intelligence Method to Assess the Tumor Microenvironment With Treatment Outcomes for Gastric Cancer Patients After Gastrectomy. J Trans Med (2022) 20(1):1–12. doi: 10.1186/s12967-022-03298-7
18. Haralick RM, Shanmugam K, Dinstein I. Textural Features for Image Classification. IEEE Trans Systems Man Cybernetics (1973) SMC-3:610–21. doi: 10.1109/TSMC.1973.4309314
19. Ama Da Sun M, King R. Textural Features Corresponding to Textural Properties. IEEE Trans Syst Man Cybernetics (1989) 19:1264–74. doi: 10.1109/21.44046
20. Thibault G, Angulo J, Meyer F. Advanced Statistical Matrices for Texture Characterization: Application to Cell Classification. IEEE T. Bio-Med Eng (2014) 61:630–7. doi: 10.1109/TBME.2013.2284600
21. Galloway M. Texture Analysis Using Gray Level Run Lengths. Comput Graphics Image Processing (1975) 4:172–9. doi: 10.1016/S0146-664X(75)80008-6
22. Langfelder P, Horvath S. WGCNA: An R Package for Weighted Correlation Network Analysis. BMC Bioinf (2008) 9:559. doi: 10.1186/1471-2105-9-559
23. Su Y, Wang J, Shi M, Niu X, Yu X, Gao L, et al. Metabolomic and Network Analysis of Astaxanthin-Producing Haematococcus Pluvialis Under Various Stress Conditions. Bioresour Technol (2014) 170:522–9. doi: 10.1016/j.biortech.2014.08.018
24. Kunowska N, Rotival M, Yu L, Choudhary J, Dillon N. Identification of Protein Complexes That Bind to Histone H3 Combinatorial Modifications Using Super-SILAC and Weighted Correlation Network Analysis. Nucleic Acids Res (2015) 43:1418–32. doi: 10.1093/nar/gku1350
25. Pei G, Chen L, Wang J, Qiao J, Zhang W. Protein Network Signatures Associated With Exogenous Biofuels Treatments in Cyanobacterium Synechocystis Sp. PCC 6803. Front Bioeng Biotechnol (2014) 2:48. doi: 10.3389/fbioe.2014.00048
26. van Dijk PC, Jager KJ, Zwinderman AH, Zoccali C, Dekker FW. The Analysis of Survival Data in Nephrology: Basic Concepts and Methods of Cox Regression. Kidney Int (2008) 74:705–9. doi: 10.1038/ki.2008.294
27. Yu G, Li F, Qin Y, Bo X, Wu Y, Wang S, et al. GOSemSim: An R Package for Measuring Semantic Similarity Among GO Terms and Gene Products. Bioinformatics (2010) 26:976–8. doi: 10.1093/bioinformatics/btq064
28. Wu T, Hu E, Xu S, Chen M, Guo P, Dai Z, et al. Clusterprofiler 4.0: A Universal Enrichment Tool for Interpreting Omics Data. Innovation (2021) 2(3):11. doi: 10.1016/j.xinn.2021.100141
29. Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J R Stat Society Ser B: Methodol (1995) 57:289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
30. Aerts HJ, Velazquez ER, Leijenaar RTH, Parmar C, Grossmann P, Carvalho S, et al. Decoding Tumour Phenotype by Noninvasive Imaging Using a Quantitative Radiomics Approach. Nat Commun (2014) 5:4006. doi: 10.1038/ncomms5006
31. Iasonos A, Schrag D, Raj GV, Panageas KS. How to Build and Interpret a Nomogram for Cancer Prognosis. J Clin Oncol (2008) 26:1364–70. doi: 10.1200/JCO.2007.12.9791
32. Vickers AJ, Elkin EB. Decision Curve Analysis: A Novel Method for Evaluating Prediction Models. Med Decis Making (2006) 26:565–74. doi: 10.1177/0272989X06295361
33. Blanche P, Dartigues JF, Jacqmin-Gadda H. Estimating and Comparing Time-Dependent Areas Under Receiver Operating Characteristic Curves for Censored Event Times With Competing Risks. Stat Med (2013) 32:5381–97. doi: 10.1002/sim.5958
34. Pope WB. Genomics of Brain Tumor Imaging. Neuroimaging Clin N Am (2015) 25:105–19. doi: 10.1016/j.nic.2014.09.006
35. Chen Z, Liu B, Yi M, Qiu H, Yuan X. A Prognostic Nomogram Model Based on mRNA Expression of DNA Methylation-Driven Genes for Gastric Cancer. Front Oncol.584733 (2020) 10. doi: 10.3389/fonc.2020.584733
36. Wang S, Feng C, Dong D, Li H, Zhou J, Ye Y, et al. Preoperative Computed Tomography-Guided Disease-Free Survival Prediction in Gastric Cancer: A Multicenter Radiomics Study. Med Phys (2020) 47:4862–71. doi: 10.1002/mp.14350
Keywords: radiogenomic, gastric cancer, nomogram, radiomic, prognosis
Citation: Liu H, Wang Y, Liu Y, Lin D, Zhang C, Zhao Y, Chen L, Li Y, Yuan J, Chen Z, Yu J, Kong W and Chen T (2022) Contrast-Enhanced Computed Tomography–Based Radiogenomics Analysis for Predicting Prognosis in Gastric Cancer. Front. Oncol. 12:882786. doi: 10.3389/fonc.2022.882786
Received: 24 February 2022; Accepted: 16 May 2022;
Published: 22 June 2022.
Edited by:
Changqiang Wu, North Sichuan Medical College, ChinaReviewed by:
Dillip Muduly, All India Institute of Medical Sciences Bhubaneswar, IndiaZixu Yuan, The Sixth Affiliated Hospital of Sun Yat-sen University, China
Copyright © 2022 Liu, Wang, Liu, Lin, Zhang, Zhao, Chen, Li, Yuan, Chen, Yu, Kong and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tao Chen, drchentao@163.com; Wentao Kong, breezewen@163.com
†These authors have contributed equally to this work