Due to the complex mechanisms affecting anti-tumor immune response, a single biomarker is insufficient to identify patients who will benefit from immune checkpoint inhibitors (ICIs) treatment. Therefore, a comprehensive predictive model is urgently required to predict the response to ICIs. A total of 162 non-small-cell lung cancer (NSCLC) patients undergoing ICIs treatment from three independent cohorts were enrolled and used as training and test cohorts (training cohort = 69, test cohort1 = 72, test cohort2 = 21). Eight genomic markers were extracted or calculated for each patient. Ten machine learning classifiers, such as the gaussian process classifier, random forest, and support vector machine (SVM), were evaluated. Three genomic biomarkers, namely tumor mutation burden, intratumoral heterogeneity, and loss of heterozygosity in human leukocyte antigen were screened out, and the SVM_poly method was adopted to construct a durable clinical benefit (DCB) prediction model. Compared with a single biomarker, the DCB multi-feature model exhibits better predictive value with the area under the curve values equal to 0.77 and 0.78 for test cohort1 and cohort2, respectively. The patients predicted to have DCB showed improved median progression-free survival (mPFS) and median overall survival (mOS) than those predicted to have non-durable clinical benefit.Copyright © 2022 Wang, Zhang, Pan, Yi, Cui, Li, Wang, Gao, Wu, Chen, Jiang and Chu.