Background: In addition to directly lysing tumors, oncolytic viruses also induce antitumor immunity by recruiting and activating immune cells in the local tumor microenvironment. However, the activation of the immune cells induced by oncolytic viruses is always accompanied by high-level expression of immune checkpoints in these cells, which may reduce the efficacy of the oncolytic viruses. The aim of this study is to arm the oncolytic vaccinia virus (VV) with immune checkpoint blockade to enhance its antitumor efficacy. Methods: Through homologous recombination with the parental VV, an engineered VV-scFv-TIGIT was produced, which encodes a single-chain variable fragment (scFv) targeting T-cell immunoglobulin and ITIM domain (TIGIT). The antitumor efficacy of the VV-scFv-TIGIT was explored in several subcutaneous and ascites tumor models. The antitumor efficacy of VV-scFv-TIGIT combined with programmed cell death 1 (PD-1) or lymphocyte-activation gene 3 (LAG-3) blockade was also investigated. Results: The VV-scFv-TIGIT effectively replicated in tumor cells and lysed them, and prompt the infected tumor cells to secret the functional scFv-TIGIT. Compared with control VV, intratumoral injection of VV-scFv-TIGIT in several mouse subcutaneous tumor models showed superior antitumor efficacy, accompanied by more T cell infiltration and a higher degree of CD8⁺ T cells activation. Intraperitoneal injection of VV-scFv-TIGIT in a mouse model of malignant ascites also significantly improved T cell infiltration and CD8⁺ T cell activation, resulting in more than 90% of the tumor-bearing mice being cured. Furthermore, the antitumor immune response induced by VV-scFv-TIGIT was dependent on CD8⁺ T cells which mediated a long-term immunological memory and a systemic antitumor immunity against the same tumor. Finally, the additional combination of PD-1 or LAG-3 blockade further enhanced the antitumor efficacy of VV-scFv-TIGIT, increasing the complete response rate of tumor-bearing mice. Conclusions: Oncolytic virotherapy using engineered VV-scFv-TIGIT was an effective strategy for cancer immunotherapy. Administration of VV-scFv-TIGIT caused a profound reshaping of the suppressive tumor microenvironment from ‘cold’ to ‘hot’ status. VV-scFv-TIGIT also synergized with PD-1 or LAG-3 blockade to achieve a complete response to tumors with poor response to VV or immune checkpoint blockade monotherapy.

Background: The current staging system is imprecise for prognostic prediction of early-stage non–small cell lung cancer (NSCLC). This study aimed to develop a robust prognostic signature for early-stage NSCLC, allowing classification of patients with a high risk of poor outcome and specific treatment decision. Method: In the present study, a comprehensive genome-wide profiling analysis was conducted using a retrospective pool of early-stage NSCLC patient data from the previous datasets of Gene Expression Omnibus (GEO) including GSE31210, GSE37745, and GSE50081 and The Cancer Genome Atlas (TCGA). Cox proportional hazards models were implemented to determine the association between gene expression levels and overall patient survival in each dataset. The common genes among all datasets were selected as candidate prognostic genes. A risk score model was developed and validated using four independent datasets and the entire cohort. The Kaplan-Meier with log-rank test was used to assess survival difference. Results: A univariate Cox proportional hazards regression analysis for each dataset showed that a total of 2280 genes in GSE31210, 762 genes in GSE37745, 871 genes in GSE50081, and 666 genes in TCGA were identified as candidate protective genes, while overall 2131 genes in GSE31210, 913 in GSE37745, 1107 in GSE50081, and 997 in TCGA were identified as candidate risky genes. There were 8 common genes associated with overall survival, including 7 mRNA and 1 lncRNA. By using the Step-wise multivariate Cox analysis, an 8-gene prognostic signature (CDCP1, HMMR, TPX2, CIRBP, HLF, KBTBD7, SEC24B-AS1, and SH2B1) for early-stage NSCLC was developed. Patients in the high-risk group had shorter overall survival than those in the low-risk group. Multivariate regression and stratified analysis suggested that the prognostic power of the 8-gene signature was independent of other clinical factors. Furthermore, the 8-gene signature achieved AUC values of 0.726, 0.701, 0.725 and 0.650 in GSE31210, GSE37745, GSE50081 and TCGA, respectively. Moreover, the combination of the 8-gene signature and the stage resulted to a better patient classification for survival prediction and treatment decision. Conclusion: This study developed a robust gene signature with great value for prognostic prediction in early-stage NSCLC, which may contribute to patient classification and personalized treatment decisions.

An accurate and robust gene signature is of the utmost importance in assisting oncologists to make a more accurate evaluation in clinical practice. In our study, we extracted key mRNAs significantly related to colorectal cancer (CRC) prognosis and we constructed an expression-based gene signature to predict CRC patients’ survival. mRNA expression profiles and clinicopathological data of colon adenocarcinoma (COAD) cases and rectum adenocarcinoma (READ) were collected from The Cancer Genome Atlas database to investigate gene expression alteration associated to the prognosis of CRC. Differentially expressed mRNAs (DEMs) were detected between COAD/READ and normal tissue samples. Relying on a univariate and multivariate Cox regression analyses, a mRNA panel signature was established and used for predicting the overall survival (OS) in CRC patients. Receiver operating characteristic curve was used to evaluate the prognosis performance of our model through calculating the AUC values corresponding to the 3-year and 5-year survival. To assess the performance of gene signature in the given cancer subgroups (CRC entire cohort, COAD cohort, and READ cohort), a stratified analysis was carried out according to clinical factors. A total of 5341 and 5594 DEMs were collected from COAD vs. normal tissue samples, and READ vs. normal samples respectively. A univariate regression analysis for the common DEMs between COAD and READ cohorts resulted in 14 common mRNAs related to OS. The multivariate Cox regression analysis revealed that 6 of these mRNAs (EPHA6, TIMP1, IRX6, ART5, HIST3H2BB, and FOXD1) had significant prognostic value allowing the discrimination between high- and low-risk patients, implying poor and good outcomes, respectively. The stratified analysis identified 6-gene signature as an independent prognostic signature in predicting CRC patients’ survival. The 6-gene signature could act as an independent biomarker for survival prediction of CRC patients.

Glioblastoma (GBM) is an aggressive tumor of the central nervous system that has poor prognosis despite extensive therapy. Therefore, it is essential to identify a gene expression-based signature for predicting GBM prognosis. The RNA sequencing data of GBM patients from the Chinese Glioma Genome Atlas (CGGA) and The Cancer Genome Atlas (TCGA) databases were employed in our study. The univariate and multivariate regression models were utilized to assess the relative contribution of each gene to survival prediction in both cohorts, and the common genes in two cohorts were identified as a final prognostic model. A prognostic risk score was calculated based on the prognostic gene signature. This prognostic signature stratified the patients into the low- and high-risk groups. Multivariate regression and stratification analyses were implemented to determine whether the gene signature was an independent prognostic factor. We identified a 6-gene signature through univariate and multivariate regression models. This prognostic signature stratified the patients into the low- and high-risk groups, implying improved and poor outcomes respectively. Multivariate regression and stratification analyses demonstrated that the predictive value of the 6-gene signature was independent of other clinical factors. This study highlights the significant implications of having a gene signature as a prognostic predictor in GBM, and its potential application in personalized therapy.

The high mortality of patients with non-small cell lung cancer (NSCLC) emphasizes the necessity of identifying a robust and reliable prognostic signature for NSCLC patients. This study aimed to identify and validate a prognostic signature for the prediction of both disease-free survival (DFS) and overall survival (OS) of NSCLC patients by integrating multiple datasets. We firstly downloaded three independent datasets under the accessing number of GSE31210, GSE37745 and GSE50081, and then performed an univariate regression analysis to identify the candidate prognostic genes from each dataset, and identified the gene signature by overlapping the candidates. Then, we built a prognostic model to predict DFS and OS using a risk score method. Kaplan–Meier curve with log-rank test was used to determine the prognostic significance. Univariate and multivariate Cox proportional hazard regression models were implemented to evaluate the influences of various variables on DFS and OS. The robustness of the prognostic gene signature was evaluated by re-sampling tests based on the combined GEO dataset (GSE31210, GSE37745 and GSE50081). Furthermore, a The Cancer Genome Atlas (TCGA)-NSCLC cohort was utilized to validate the prediction power of the gene signature. Finally, the correlation of the risk score of the gene signature and the Gene set variation analysis (GSVA) score of cancer hallmark gene sets was investigated. We identified and validated a six-gene prognostic signature in this study. This prognostic signature stratified NSCLC patients into the low-risk and high-risk groups. Multivariate regression and stratification analyses demonstrated that the six-gene signature was an independent predictive factor for both DFS and OS when adjusting for other clinical factors. Re-sampling analysis implicated that this six-gene signature for predicting prognosis of NSCLC patients is robust. Moreover, the risk score of the gene signature is correlated with the GSVA score of 7 cancer hallmark gene sets. This study provided a robust and reliable gene signature that had significant implications in the prediction of both DFS and OS of NSCLC patients, and may provide more effective treatment strategies and personalized therapies.