Deep segmentation networks predict survival of non-small cell lung cancer
Stephen Baek Yusen He Bryan G. Allen John M. Buatti Brian J. Smith Ling Tong Zhiyu Sun Jia Wu Maximilian Diehn Billy W. Loo Kristin A. Plichta Steven N. Seyedin Maggie Gannon Katherine R. Cabel Yusung Kim Xiaodong Wu


Towards a more complex description of chemical profiles in exoplanet retrievals: A 2-layer parameterisation
    [
    December 30, 2023
===========================================================================================================





 empty





§ INTRODUCTION


According to the World Health Organization (WHO), lung cancer remains the leading cause of cancer-related deaths worldwide, with 2.1 million new cases diagnosed and 1.8 million deaths in 2018<cit.>. NSCLC accounts for 80-85% of lung cancer diagnoses<cit.> and the five-year survival rate of NSCLC remains low (23%) compared to other leading cancer sites such as colorectal (64.5%), breast (89.6%), and prostate (98.2%)<cit.>. Historically, the tumor, nodes, and metastases (TNM) staging system has served as the major prognostic factor in predicting therapeutic outcomes, but it does not differentiate responders and non-responders within the same stage<cit.>. The maximum and the mean of standard uptake values (SUV_MAX and SUV_MEAN) have been reported for their correlation with survival<cit.>, but are of limited clinical value due to their unsatisfactory predictive power and lack of robustness<cit.>. Other prognostic markers have also been studied, including TLG, which incorporates metabolic tumor volume (MTV) and metabolic activity (TLG = MTV × SUV_MEAN). Reports<cit.> suggest that TLG may have better prognostic power than SUV_MAX or SUV_MEAN. These metrics, however, are not optimal and do not provide a comprehensive image-based analysis of tumors<cit.>. More recently, radiomics approaches, which employ semi-automated analysis based on a few hand-crafted imaging features describing intratumoral heterogeneity, demonstrated higher prognostic power<cit.>. However, these features still have limited predictive power ranging between 0.5 and 0.79 in terms of the area under the ROC (AUC)<cit.>. Recent literature on CNNs demonstrates their strong potential for cancer prognostication<cit.>, however the clinical implications of deep learning remain questioned due to the limited interpretability of CNNs.

Here, we propose an interpretable and highly accurate framework to solve this problem by capitalizing on the unprecedented success of deep convolutional neural networks (CNN). More specifically, we investigate U-Net<cit.>, a convolutional encoder-decoder network that has demonstrated exceptional performance in tumor detection and segmentation tasks. Illustrated in Fig, <ref>a, U-Net takes a three-dimensional (3D) volume image as an input, processes it through a “bottleneck layer” where the image features are compressed, and reconstructs the image into a binary segmentation map indicating a pixel-wise tumor classification result. Here, we focused on the information encoded at the bottleneck layer, which contains rich visual characteristics of the tumor and hypothesized that the encoded information at this layer might be relevant to the tumor malignancy, and thus cancer survival, which is the central hypothesis of this paper.





§ RESULTS

In prior studies<cit.>, we analyzed PET/CT images of 96 NSCLC patients that were obtained within 3 months prior to SBRT, whose summary statistics are illustrated in Extended Data Fig. <ref> and Extended Data Table <ref>. For each volume image, the region of interest (ROI) with a dimension of 96 mm × 96 mm × 48 mm was set around each SBRT treated tumor location and the image was cropped to the ROI volume. Two separate U-Net models were trained to perform tumor segmentation in PET and CT images, respectively. Each of the models was trained with 38 ROI images and the corresponding physician contours, but no other information such as survival time was provided. Segmentation performance was tested on 22 independent ROI images that were not included in training, and the average Sørensen-Dice similarity coefficients (DSC) were 0.861 ± 0.037 (mean ± std) and 0.828 ± 0.087 for CT and PET, respectively. After training, each U-Net model learned to encode 55,296 features at the bottleneck layer for each patient, resulting in a total of 110,592 features per patient.








These features are an intermediate throughput of U-Net, and are then decoded to generate an automated segmentation in the network. It is likely that these features summarize some rich structural and functional geometry of the intratumoral and peritumoral area, some of which might be relevant to cancer survival.  We test this proposition by conducting validation studies and examining the statistical prognostic power of those features. One challenge here is that the number of features (d=110,592) is substantially larger than the number of cases, making a statistical analysis prone to overfitting. To this end, we first reduce the number of features via an unsupervised feature selection process in which the survival information is hidden. The k-medoids clustering method <cit.> is employed to serve this purpose, as the method is known to be able to select representative features from a large pool of inter-correlated features in similar literature <cit.>. For the feature similarity measure in k-medoids clustering, we use the Pearson correlation distance defined as 1-R, where R is the Pearson correlation coefficient. The optimal number of clusters is determined by the Silhouette method <cit.>. Finally, given the optimal k clusters, the medoids of the clusters are selected to be the representative features of the clusters and form an unsupervised feature pool. We choose the final set of features via the least absolute shrinkage and selection operator (LASSO)<cit.> on the preliminary feature pool obtained from the k-medoids method. LASSO attempts to select a few features among the cluster centers that have a strong relationship with the survival outcomes using the fixed regularization parameter. The selected features then serve as regressors in the logistic regression model to predict the cancer survival outcome. Since the performance of k-medoids clustering is also dependent on a random initialization of the cluster centers, we tested 10 times for each k and computed the mean of the summation of the inner cluster distances to ensure consistency of the clustering outcome. We then computed the curve of the Silhouette value with respect to the number of clusters in total to determine the optimal number of clusters. The medoids of the feature clusters are selected to form a reduced feature set.










We then examined the selected features and their prognostic power via cross-validation on the NSCLC data set collected at the University of Iowa Hospitals and Clinics (UIHC). As summarized in Extended Data Fig <ref> and Extended Data Table <ref>, the UIHC data set is comprised of primary and follow-up PET/CT images of total 96 NSCLC cases with their survival status and other clinical meta data. Information on the cause of death is also available so that the deceased cases can be further broken down to overall and disease-specific deaths. On this data set, we aim to predict four survival categories, namely, 2-year overall survival (2OS), 5-year overall survival (5OS), 2-year disease-specific survival (2DS), and 5-year disease-specific survival (5DS). Total N=96 cases qualify for 2OS category, 74 for 5OS, 92 for 2DS, and 45 for 5DS, depending on the survival status and the cause of death. Other clinical meta data, such as sex, age, smoking history, and tumor, lymph nodes, and metastasis (TNM) staging, exist in the data set and may provide an improved prognostic power when added as regressors. However, we exclude all other parameters but U-Net encoded image features, in order to focus the analysis on the image features only. Only primary images are encoded via U-Net and used for prediction. Follow-up images are reserved for comparison and discussion later in this paper. Each cross validation experiment is comprised of unsupervised feature selection using k-medoids clustering, LASSO-based feature selection, and training of the logistic regression model. These tasks are performed independently from the other cross-validation experiments. In each experiment, we measure the accuracy, sensitivity, specificity, and the area under the receiver operating characteristic curve (AUC) per each survival category. We compute the average and the standard deviation of these performance metrics across the cross-validation experiments to derive the estimated performance metrics and their 95% confidence intervals (95% CI).

In such cross-validation experiments, estimated AUC of the proposed prediction model is 0.88 (95% CI: 0.80-0.96) for the prediction of 2OS. For other survival criteria, the estimated AUCs are similar, namely 0.89 (95% CI: 0.85-0.93) for 5OS, 0.86 (95% CI: 0.81-0.91) for 2DS, and 0.88 (95% CI: 0.81-0.95) for 5DS. Note that the estimated AUC values of the conventional TLG and other radiomics markers <cit.> range between 0.60 and 0.83 on the same data set. A more recent deep learning method reported in Hosny <cit.> produces AUCs between 0.70∼0.73 on the same data set. A graphical illustration of the result as well as the full set of performance metrics are available in Fig. <ref> and Extended Data Table <ref>, respectively.





Moreover, we further validated the result on an extramural dataset provided by the Stanford Cancer Institute. The Stanford data set is comprised of primary CT images of 26 NSCLC cases which received SBRT, of which 18 survived and 8 died according to 2 year OS and 1 survived and 25 died according to 5 year OS. Neither PET images nor disease-specific survival information were available in the Stanford data set. Training of the U-Net, feature selection (k-medoids and LASSO), and construction of the survival model are performed only on the UIHC data set, and the Stanford data set is used for validation only. In this setting, the estimated AUC is 0.87 (95% CI: 0.80-0.94) for 2 year OS, and 0.90 (95% CI: 0.82-0.98) for 5 year OS. More detailed results are illustrated in Fig. <ref> and Extended Data Table <ref>.






Meanwhile, we visualized the features learned by U-Net to develop an intuitive understanding of what those prognostic markers represent. The image features learned by U-Net are essentially artificial neurons in deep neural networks. In principle, we can visualize a neuron by showing multitudes of image patterns and observing which image pattern activates the neuron the most. Practically, we employ an optimization-based approach<cit.> where the objective is to maximize an individual neuron's activation value by manipulating the input image pattern:

    𝐗^*= max_𝐗 q(𝐗 | 𝐖, 𝐛),

where q( · | 𝐖, 𝐛) is the U-Net encoder with the trained model parameters 𝐖 and 𝐛, and 𝐗 is the input image pattern. Displayed in Fig. <ref> and Extended Data Fig. <ref> are visualizations of features captured by artificial neurons.




We also visualized which regions in the patient images predicted low survival probability. We employed a guided gradient backpropagation approach<cit.>. The main idea of the guided backpropagation algorithm is to compute ∂ P/∂ x_i,j,k where P is the probability of death and x_i,j,k is a voxel value at the position (i,j,k) in the patient image. The gradient ∂ P/∂ x_i,j,k can be interpreted as the change of the death probability when the voxel x_i,j,k changes to a different value. If the voxel was not so significant in predicting death, the gradient value would be small, where as if the voxel played an important role in predicting high probability of death, the gradient value would be greater. Displayed in Fig. <ref> and Extended Data Figs. <ref>-<ref> are heat maps representing the gradient. Heated regions (red) are the areas that lowered the probability of survival whereas the other areas (blue) are the ones that had negligible effect on the survival.

















§ DISCUSSION


As illustrated Fig. <ref> and Extended Data Table <ref>, the contrast in the predictive performance between the U-Net features and the conventional imaging features (SUV_MEAN, SUV_MAX and TLG) was evident, quantitatively proving the strong prognostic power of the U-Net features. There also was a noticeable enhancement of prognostic performance when the U-Net features were compared with a recent deep learning approach as in Hosny <cit.> The same trend could be observed in Fig. <ref> and Extended Data Table <ref> where we validated the U-Net features against an extramural data set, confirming the enhanced prognostic performance of the U-Net features. Here, it is worth reemphasizing that the U-Net was trained without any survival-related information, and, thus it is highly unlikely that the U-Net-learned features were overfitted to the survival data or biased towards them. Nonetheless, while these U-Net features were identified independently from the survival data, the U-Net features demonstrated strong evidence of correlation and thus prognostic power for NSCLC survival as discussed above.

Here, a justification should be necessary on the rationale behind taking a “detour” by training a segmentation network first, then extracting prognostic image features, and training a survival model, as opposed to directly training a CNN model to the survival data as in other literature.<cit.> To this end, we make the following arguments. First, a CNN model trained on a segmentation task is more robust to overfitting and is more generalizable, as the image features themselves are developed from an unsupervised training. A direct prediction of survival from an image tends to be less robust and less generalizable, as the millions of parameters in a CNN can easily misconstrue the trend, and this process is difficult to control. Secondly, segmentation-trained features contain more structural, geometric features that are human-interpretable, whereas direct-trained features tend to be biased towards less-intuitive and fuzzy texture patterns. It is, in fact, well-understood now that CNNs tend to be biased towards texture rather than shape, while the human cognitive process works in the opposite.<cit.> For example, Geirhos et al.<cit.> conducted an experiment where the texture of a cat in an image was replaced by the texture of an elephant—CNNs started to discern the image as an elephant while humans still categorized the image as a cat. From this, we can infer that the direct training of a CNN for survival prediction would result in a development of less sensical features lacking geometric information. In contrast, segmentation-trained features are incentivized to focus on geometric structures and shapes rather than textures, and thus are more sensical.

The above arguments are indeed supported by the visualization results in Fig. <ref> and Extended Data Fig. <ref>, where the benchmark CNN method appears to be capturing some noisy texture patterns. In contrast, many of the U-Net features appear to be capturing more sensical and interpretable geometric shapes from the image, such as tumor-like blobs (e.g. C00048, C25988, P39051, P47258) or heterogeneity of the tumor (e.g. C08680). Interestingly, some of the U-Net features, for example C01777 and C37399, were looking for tube-like structures nearby the tumor-like blobs, which might be capturing blood vessels and lymphatics in the peritumoral area. This is consistent with the widely accepted clinical knowledge that tumors can show enhanced growth in the presence of nearby vessels and lymphatics as they carry nutrition to supply the tumoral growth.


Meanwhile, we further visualized the proposed prediction model by creating a heat map highlighting regions that predicted low survival probability. The map generated by the guided back-propagation algorithm confirmed that the prediction model was looking at clinically sound regions for making predictions by responding more to tumor and surrounding tumor regions than other regions. In addition to the results presented above, this is yet more evident that the prediction model has produced generalizable features and rules for making prognostic prediction, and is thus not overfitted. Moreover, through comparison with post-SBRT CT images and clinical records of the patients, we observe that the heat map has the potential to identify regions of progressions or recurrence. For the case illustrated in Fig. <ref> as an example, a small nodule had been found initially at the progression region, but was not treated by SBRT as the region was not originally identified as tumor by PET. However, this region was identified as a high risk region on heat map by the U-Net as well as the primary tumor region of SBRT. Tumor recurrence was found on post-SBRT CT images and matched the highlighted risk region on pre-SBRT U-Net generated heat map. From such observations, the heat map visualization has the potential to identify regions at high risk for tumor progression or recurrence that could be utilized for the purpose of assisting patient-tailored treatment planning in the future. For this reason, we believe this indicates more rigorous risk map developments and requires quantitative follow-up validations, which will be our subsequent project.

In summary, we discovered that the U-Net segmentation algorithm trained for automated tumor segmentation on PET/CT, codifies rich structural and functional geometry at the bottleneck layer. These codified features, in turn, could be used for survival prediction in cancer patients even though the U-Net was trained without any survival-related information. The survival model based on such U-Net features demonstrated significantly higher predictive power compared to conventional PET-based, metabolic burden metrics such as TLG or relatively recent hand-crafted radiomics approaches. The validity of this discovery was further confirmed by the validation on an extramural data set provided by the Stanford Cancer Institute. Furthermore, we visualized the survival-related U-Net features and observed that they were indeed depicting intratumoral and/or peritumoral structures that had been previously acknowledged as potentially relevant to cancer survival. Our approach awaits a further validation against a larger number of observations and in a larger variety of cancer types. However, there was not enough clinical evidence to conclude that the visualization of the U-Net features may identify potential regions of recurrence and metastasis and, thus, a follow-up study is suggested. Our findings may suggest a new starting point for quantitative image-based cancer prognosis with a great deal of important new knowledge to be discovered.







§ METHODS


Subjects     PET-CT images of NSCLC patients who received SBRT at the University of Iowa Hospitals and Clinics were investigated in this study. The images were obtained using a dual PET-CT scanner (Biograph 40 PET/CT, Siemens Medical Solutions USA, Inc., Hoffman Estates, IL). CT image used for SBRT planning were co-registered with CT images of PET-CT datasets using a deformable image registration (DIR). Using the deformation map (vectors), PET images were resampled with the primary CT images of SBRT plan. Post-SBRT follow-up images were used only for qualitative validation of the visualization results. The gross tumor volume (GTV) for each of the images was delineated by radiation oncologists on both CT and PET images with the guidance of the corresponding images in the other modality. All tumor contouring was completed using VelocityAI (Varian Medical System, Inc., Palo Alto, CA).

A total of 96 cases (Male=44; Feamle=52) were investigated in this study. Patients' group stages vary in sub-categories including 41 in stage I, 10 in stage II, 16 in stage III, 29 in stage IV. Meanwhile, the histology was confirmed by a thoracic pathologist based off visual interrogation of the biopsied or surgically resected specimen. Histologies among the 96 patients include 48 adenocarcinoma, 41 squamous cell, 1 adenosquamous, 3 metastases from previous NSCLC, and 1 without biopsy. For overall survival from SBRT, the 2-year survival rate is 54% and the 5-year survival rate is 6%. For overall survival from diagnosis, the 2-year survival rate is 66% and the 5-year survival rate is 28%. For disease-specific survival from SBRT, the 2-year survival rate is 61% and the 5-year survival rate is 23%. For disease-specific survival from diagnosis, the 2-year disease specific survival rate is 73% and the 5-year disease specific survival rate is 51%. Among a total 96 patients, the qualified portion utilized in this research in each survival sub-category has been visualized and discussed in Extended Data Fig. <ref> respectively.

In this study, the total of 96 utilized NSCLC patients were retrospectively analyzed after approval from the University of Iowa Institutional Review Board (IRB: 200503706; Name: The utility of imaging in cancer: staging, assessment of treatment response, and surveillance). All data collection, experimental procedures, and methods applied are in accordance with the relevant guidelines and regulations. All patients consented for the use of their clinical information and medical images. All participants enrolled in this study signed an informed consent developed and approved by the University of Iowa Institutional Review Board. All scans were in digital imaging and communications in medicine (DICOM) format and de-codified.

Data Processing     Each pair of the co-registered PET-CT images were resampled with an isotropic spacing in all three dimensions. Each image was then cropped into a fixed size of 96 × 96 × 48 voxels where each voxel represents 1 cubic millimeter (mm^3) after resampling. Within the resampled voxels, intensity values were clipped to the range of [-500, 200] for all CTs and [0.01, 20] for all PETs, to remove outliers. More details are described in our previous work<cit.>. 

U-Net Features     In our previous work<cit.>, two independent 3D U-Net networks were constructed and trained for automated tumor segmentation in PET and CT images, respectively. The U-Net networks comprise two major components: the encoder network and the decoder network (Fig. <ref>). The encoder network takes a 96 × 96 × 48 volume image as an input. The first convolution layer produces 32 features attached to each voxel, representing low-level visual cues such as lines, edges, and blobs. These features are then down-sampled by half in all three dimensions and the down-sampled volume is fed into the second convolution layer, which then produces 64 features per each voxel. This is repeated three more times increasing the number of features to 128, 256, and finally 512, while the volume size is reduced by half in all three dimensions each time. The final 6 × 6 × 3 × 512 (=55,296) features that the encoder network produces are an abstract, high-level summary of the input image, which is then decoded by the symmetric decoder to produce the binary segmentation map (1: tumor, 0: none). The convolutional kernels are of a size 3 × 3 × 3 across all layers and the max-pooling layers of a 2 × 2 × 2 window size with a stride of 2 were used for down-sampling. 

 To train the PET-CT segmentation U-Net networks, we utilized co-registered PET-CT scan pairs from 60 patients with primary NSCLC. For each PET-CT image, the slice image size is 512 × 512 and the number of slices varies from 112 to 293. The tumor contour on each of the PET and CT scans were labeled by physicians as groundtruth. In data preprocessing, all PET-CT images are resampled with an isotropic spacing of 1 × 1 × 1 in voxels and then cropped at a fixed size of 3D volumes (96 × 96 × 48) centered on the mass gravity of each tumor. 

All 60 PET-CT scan pairs were split into a training data set with 38 patients and testing data set with 22 patients. Data augmentation was performed using simple translation, rotation and flip operations and the augmented training set has over 3000 3D PET-CT scan pairs respectively.

The 3D-UNet was trained using open source TensorFlow package and ran on NVIDIA GeForce GTX 1080
Ti GPU with 11GB of memory. The Adam optimization method was utilized with a mini-batch size of 1 and for 20 epochs. To prevent overfitting, the weight decay and early-stop techniques were adopted to obtain the best performance on the test set where the DSC value was computed.



After the U-Nets had been trained, PET and CT images of each patient are fed into the U-Nets to produce 55,296 features per imaging modality. These 55,296 features extracted from the U-Net encoder were used for the analyses throughout the paper.  The schematic diagram with respect to the methodologies applied in this research for the feature analysis has been illustrated in Fig. <ref>.

k-medoids feature clustering     In this research, several k-medoids clustering experiments were conducted on the training dataset for the cross-validation experiments. Pearson's correlation distance was employed as the distance metric for clustering the features as expressed in:

    𝒟(𝐗,𝐘) =  1- Cov(𝐗,𝐘)/√(Var(𝐗)Var(𝐘))

where X and Y denote the two distinct features; Cov() denotes the covariance of the two features and Var() is the variance of the underlying feature. The sum of inner cluster distances were computed by setting various k values, and the optimal number of clusters were determined by the Silhouette method <cit.>. The medoids of all clusters were selected as candidate features to construct the survival prediction models.


Feature Selection     The LASSO regression algorithm was employed to narrow down the scope of analysis to survival-related features from the medoids of clusters obtained from clustering results. The LASSO algorithm used in this study is expressed as:

    min_β y - β X ^2_2 + λβ_1

where y denotes the survival outcome (1: alive, 0: dead), X is a vector containing all latent variables extracted from the U-Net network, β denotes the coefficient of regression, and λ is the penalty coefficient. The L1-norm in the second term penalizes the selection of redundant variables. The parameter λ was determined via cross validation on the training dataset. Latent variables that survived the L1-penalty with the best λ were selected for the logistic regression model to predict the survival outcome. Using the LASSO-selected variables, we applied logistic regression to estimate the coefficients and predict the survival outcome.

Logistic Regression     We formulate survival as a dummy variable. The task of predicting survival outcome can then be formulated as a binary class probability prediction problem and we select the linear logistic regression model as our statistical model:

    y = { 1 + e^-(β_0 + ∑_i=1^pβ_i x_i)}^-1

where y denotes the predicted survival probability, x_i is the i-th LASSO-selected U-Net feature, and β_i is the regression coefficient. The performance of a prediction model was measured via 6-fold cross validation. In experiment, the samples were split into training and test sets. The models were trained with the training set and the test sets were left out for validation. The proportion of survival and death cases were controlled to be equal in the test sets. Reported performance metrics in this paper are based on the statistics of the test set validations across 6-fold cross validation.

Visualization     An activation maximization scheme<cit.> was employed to visualize the LASSO-selected U-Net features. For a trained U-Net encoder q( · | 𝐖, 𝐛), neurons at the bottleneck layer corresponding to the LASSO-selected features were denoted as q_i. Then, Eq. <ref> was solved for each individual neuron via gradient ascent:

    𝐗^(k+1) = 𝐗^(k) + γ^(k)∇ q_i (𝐗^(k)),

where 𝐗^(k) is the current solution at k-th iteration and γ^(k) is a step length. We set γ^(k) as 1/σ^(k) where σ^(k) denotes the standard deviation of the gradients. The gradient ∇ q_i was computed using the standard backpropagation algorithm. The initial image 𝐗^(0) was initialized with random voxel values following the Gaussian distribution 𝒩(128,1). Displayed in Fig. <ref>a-b are the final solution 𝐗^* after 20 iterations.

Moreover, we also visualized a risk map by evaluating each voxel's contribution to the prediction of survival. We employed a guided backpropagation approach similar to Selvaraju <cit.>. For each voxel in the input image, with marginal change of the survival probability with respect to the voxel's intensity, defined as ∂ P/∂ x_i,j,k, where P is the probability of death and x_i,j,k is a voxel value at the position (i,j,k). In the guided backpropagation process, we rectified the gradient by dropping the negative gradient values to focus on the “risk”. This was achieved by applying rectified linear unit (ReLU) activation when the values were backpropagated from node to node:

    α^(m) =  max( ∂(P)/∂A^(m)_i,j,k,0
        )

where A^(m) denotes the activation map corresponding to the m-th convolutional kernel at the bottleneck encoding. Note that only the LASSO-selected features were involved in the survival model P such that ∂(P)/∂A^(m)_i,j,k is zero most of the time. Finally, the risk map ℛ was defined as a linear combination of all activation maps at the bottleneck layer with the coefficients α^(m) obtained from the above:

    ℛ(𝐗) = ∑_mα^(m)A^(m)(𝐗).






§ ACKNOWLEDGEMENTS

The authors want to express their gratitude to Dr. Ruijiang Li at Stanford University for his assistance in coordinating the collaboration with Stanford University and for providing the extramural data set. They also thank Dr. Sanjay Aneja at Yale University for providing insight and expertise that greatly assisted the research.
Research reported in this publication was supported by the National Cancer Institute (NCI) of the National Institutes of Health (NIH) under award number 1R21CA209874 and partially by U01CA140206 and P30CA086862.




§ AUTHOR CONTRIBUTIONS STATEMENT


Y.K. and X.W. conceived and managed the experiments, S.B., Y.H. and L.T. conducted and analyzed the experiments, Y.H. and Z.S. performed visualization experiments, S.B. and B.J.S. advised the statistical analysis, Y.H., M.G., and K.R.C. processed the raw data, B.G.A., J.M.B., and K.A.P. provided the contours for training of the U-net, B.G.A., J.M.B., and S.N.S. advised clinical discussions. R.L., J.W., M.D., and B.L. provided and processed the extramural data set. All authors reviewed the manuscript.









§ ADDITIONAL INFORMATION




 §.§ Competing interests


Dr. Buatti's work has been funded in part by the National Institutes of Health/National Cancer Institute grants P01 CA217797001A1, UL1 TR002537,
U01 CA140206, and 1R21CA209874.
Dr. Kim's work has been funded in part by the National Institutes of Health/National Cancer Institute grant 1R21CA209874.
Dr. Wu's work has been funded in part by the National Institutes of Health grants 1R21CA209874 and R01EB020665.
Drs. Baek, Wu, Kim, Smith, Allen, Buatti and Mr. He are the co-inventors of the provisional patent (U.S. App. No. 62/811,326. Systems And Methods For Image Segmentation And Survival Prediction Using Convolutional Neural Networks) based upon the work of this manuscript.
Dr. Plichta, Dr. Seyedin, Ms. Gannon, and Ms. Cabel declare no potential conflict of interest.



§ EXTENDED DATA

[table]name=Extended Data Table
[figure]name=Extended Data Figure





[ht]


Predictive power of linear model on 2-yr./5-yr. survival outcome


    2c
    2*Features     2c(a) TLG    2c(b) Radiomic<cit.>    2c(c) Benchmark<cit.>    2c(c) CT U-Net    2c(d) PET U-Net    2c(e) All U-Net

            2-yr.     5-yr.    2-yr.     5-yr.    2-yr.     5-yr.    2-yr.     5-yr.     2-yr.     5-yr.     2-yr.     5-yr.
    
4*Accuracy 
 (95% CI)     2*OS     0.53     0.55     0.62     0.65     0.69     0.75     0.84     0.86     0.83     0.87     0.86     0.89

           (0.46-0.60)     (0.47-0.63)     (0.54-0.70)     (0.55-0.75)     (0.65-0.73)     (0.70-0.80)     (0.81-0.87)     (0.80-0.92)     (0.79-0.87)     (0.80-0.94)     (0.79-0.93)     (0.85-0.94)

         2*DS     0.57     0.57     0.64     0.80     0.73     0.76     0.86     0.89     0.87     0.88     0.88     0.88

           (0.54-0.60)     (0.52-0.62)     (0.59-0.69)     (0.75-0.85)     (0.67-0.80)     (0.69-0.83)     (0.80-0.92)     (0.84-0.94)     (0.81-0.93)     (0.83-0.93)     (0.81-0.95)     (0.82-0.94)

    4*Sensitivity 
 (95% CI)     2*OS       0.40       0.32      0.44      0.60      0.63      0.69      0.72      0.81      0.75      0.82      0.74      0.82

           (0.27-0.53)     (0.20-0.44)     (0.34-0.54)     (0.53-0.67)     (0.57-0.70)     (0.62-0.76)     (0.61-0.83)     (0.76-0.88)     (0.68-0.82)     (0.74-0.90)     (0.67-0.81)     (0.76-0.88)

         2*DS      0.17      0.33      0.40      0.91      0.68      0.70      0.77      0.84      0.79      0.85      0.79      0.86

           (0.11-0.23)     (0.23-0.43)     (0.31-0.49)     (0.86-0.96)     (0.60-0.76)     (0.65-0.75)    
    (0.67-0.87)     (0.78-0.90)     (0.70-0.88)     (0.78-0.92)     (0.73-0.85)     (0.80-0.92)

    4*Specificity 
 (95% CI)     2*OS      0.66      0.78      0.80      0.68      0.77      0.80      0.91      0.90      0.90      0.93      0.91      0.92

           (0.55-0.77)     (0.71-0.85)     (0.68-0.92)     (0.60-0.76)     (0.70-0.84)     (0.75-0.85)     (0.78-0.99)     (0.83-0.97)     (0.80-0.99)     (0.84-0.99)     (0.84-0.98)     (0.84-0.99)

         2*DS      0.97      0.81      0.88      0.70      0.78      0.81      0.92      0.92      0.90      0.91      0.91      0.91

           (0.95-0.99)     (0.71-0.91)     (0.82-0.94)     (0.62-0.78)     (0.72-0.84)     (0.76-0.86)     (0.81-0.99)     (0.86-0.98)     (0.85-0.95)     (0.84-0.98)     (0.83-0.99)     (0.86-0.96)

    4*AUC 
 (95% CI)     2*OS     0.60     0.63     0.61     0.69      0.70      0.72      0.83      0.85      0.82      0.87      0.88      0.89

           (0.46-0.74)     (0.56-0.70)     (0.48-0.74)     (0.60-0.78)     (0.65-0.75)     (0.66-0.78)     (0.75-0.91)     (0.75-0.95)     (0.73-0.89)     (0.79-0.95)     (0.80-0.96)     (0.85-0.93)

         2*DS     0.68     0.68     0.64     0.83      0.72      0.73      0.85      0.85      0.86      0.88      0.86      0.88

           (0.61-0.75)     (0.60-0.76)     (0.57-0.71)     (0.79-0.87)     (0.67-0.77)     (0.68-0.78)     (0.75-0.95)     (0.77-0.92)     (0.79-0.93)     (0.80-0.96)     (0.81-0.91)     (0.81-0.95)

    




















