Although immunotherapy has revolutionized cancer management, most patients do not derive benefits from it. Aiming to explore an appropriate strategy for immunotherapy efficacy prediction, we collected 6251 patients\' transcriptome data from multicohort population and analyzed the data using a machine learning algorithm. In this study, we found that patients from three immune gene clusters had different overall survival when treated with immunotherapy (P < 0.001), and that these clusters had differential states of hypoxia scores and metabolism functions. The immune gene score showed good immunotherapy efficacy prediction (AUC was 0.737 at 20 months), which was well validated. The immune gene score, tumor mutation burden, and long non-coding RNA score were further combined to build a tumor immune microenvironment signature, which correlated more strongly with overall survival (AUC, 0.814 at 20 months) than when using a single variable. Thus, we recommend using the characterization of the tumor immune microenvironment associated with immunotherapy efficacy via a multi-omics analysis of cancer.

Alternative splicing (AS) events modulate certain pathways and phenotypic plasticity in cancer. Although previous studies have computationally analyzed splicing events, it is still a challenge to uncover biological functions induced by reliable AS events from tremendous candidates. To provide essential splicing event signatures to assess pathway regulation, we developed a database by collecting two datasets: (i) reported literature and (ii) cancer transcriptome profile. The former includes knowledge-based splicing signatures collected from 63,229 PubMed abstracts using natural language processing, extracted for 202 pathways. The latter is the machine learning-based splicing signatures identified from pan-cancer transcriptome for 16 cancer types and 42 pathways. We established six different learning models to classify pathway activities from splicing profiles as a learning dataset. Top-ranked AS events by learning model feature importance became the signature for each pathway. To validate our learning results, we performed evaluations by (i) performance metrics, (ii) differential AS sets acquired from external datasets, and (iii) our knowledge-based signatures. The area under the receiver operating characteristic values of the learning models did not exhibit any drastic difference. However, random-forest distinctly presented the best performance to compare with the AS sets identified from external datasets and our knowledge-based signatures. Therefore, we used the signatures obtained from the random-forest model. Our database provided the clinical characteristics of the AS signatures, including survival test, molecular subtype, and tumor microenvironment. The regulation by splicing factors was additionally investigated. Our database for developed signatures supported retrieval and visualization system.

Cancer-associated fibroblasts (CAFs) can exert their immunosuppressive effects by secreting various effectors that are involved in the regulation of tumor-infiltrating immune cells as well as other immune components in the tumor immune microenvironment (TIME), thereby promoting tumorigenesis, progression, metastasis, and drug resistance. Although a large number of studies suggest that CAFs play a key regulatory role in the development of head and neck squamous cell carcinoma (HNSCC), there are limited studies on the relevance of CAFs to the prognosis of HNSCC. In this study, we identified a prognostic signature containing eight CAF-related genes for HNSCC by univariate Cox analysis, lasso regression, stepwise regression, and multivariate Cox analysis. Our validation in primary cultures of CAFs from human HNSCC and four human HNSCC cell lines confirmed that these eight genes are indeed characteristic markers of CAFs. Immune cell infiltration differences analysis between high-risk and low-risk groups according to the eight CAF-related genes signature hinted at CAFs regulatory roles in the TIME, further revealing its potential role on prognosis. The signature of the eight CAF-related genes was validated in different independent validation cohorts and all showed that it was a valid marker for prognosis. The significantly higher overall survival (OS) in the low-risk group compared to the high-risk group was confirmed by Kaplan-Meier (K-M) analysis, suggesting that the signature of CAF-related genes can be used as a non-invasive predictive tool for HNSCC prognosis. The low-risk group had significantly higher levels of tumor-killing immune cell infiltration, as confirmed by CIBERSORT analysis, such as CD8+ T cells, follicular helper T cells, and Dendritic cells (DCs) in the low-risk group. In contrast, the level of infiltration of pro-tumor cells such as M0 macrophages and activated Mast cells (MCs) was lower. It is crucial to delve into the complex mechanisms between CAFs and immune cells to find potential regulatory targets and may provide new evidence for subsequently targeted immunotherapy. These results suggest that the signature of the eight CAF-related genes is a powerful indicator for the assessment of the TIME of HNSCC. It may provide a new and reliable potential indicator for clinicians to predict the prognosis of HNSCC, which may be used to guide treatment and clinical decision-making in HNSCC patients. Meanwhile, CAF-related genes are expected to become tumor biomarkers and effective targets for HNSCC.

Lung adenocarcinoma (LUAD) is the most prevalent lung cancer and one of the leading causes of death. Previous research found a link between LUAD and Aldehyde Dehydrogenase 2 (ALDH2), a member of aldehyde dehydrogenase gene (ALDH) superfamily. In this study, we identified additional useful prognostic markers for early LUAD identification and targeting LUAD therapy by analyzing the expression level, epigenetic mechanism, and signaling activities of ALDH2 in LUAD patients. The obtained results demonstrated that ALDH2 gene and protein expression significantly downregulated in LUAD patient samples. Furthermore, The American Joint Committee on Cancer (AJCC) reported that diminished ALDH2 expression was closely linked to worse overall survival (OS) in different stages of LUAD. Considerably, ALDH2 showed aberrant DNA methylation status in LUAD cancer. ALDH2 was found to be downregulated in the proteomic expression profile of several cell biology signaling pathways, particularly stem cell-related pathways. Finally, the relationship of ALDH2 activity with stem cell-related factors and immune system were reported. In conclusion, the downregulation of ALDH2, abnormal DNA methylation, and the consequent deficit of stemness signaling pathways are relevant prognostic and therapeutic markers in LUAD.

UNASSIGNED: Dysregulation of long noncoding RNAs (lncRNAs) has been reported to be associated with multiple tumors where they act as tumor suppressors or accelerators. The lncRNA CYTOR was identified as an oncogene involved in many cancers, such as gastric cancer, colorectal cancer, hepatocellular carcinoma, and renal cell carcinoma. However, the role of CYTOR in bladder cancer (BCa) has rarely been reported.
UNASSIGNED: Using cancer datasets from The Cancer Genome Atlas (TCGA) program, we analyzed the association between CYTOR expression and prognostic value, oncogenic pathways, antitumor immunity and immunotherapy response in BCa. The influence of CYTOR on the immune infiltration pattern in the urothelial carcinoma microenvironment was further verified in our dataset. Single-cell analysis revealed the role of CYTOR in the tumor microenvironment (TME) of BCa. Finally, we evaluated the expression of CYTOR in BCa in the Peking University First Hospital (PKU-BCa) dataset and its correlation with the malignant phenotype of BCa in vitro and in vivo.
UNASSIGNED: The results indicated that CYTOR was highly expressed in multiple cancer samples, including BCa, and increased CYTOR expression contributed to poor overall survival (OS). Additionally, elevated CYTOR expression was significantly correlated with clinicopathological features of BCa, such as female sex, advanced TNM stage, high histological grade and non-papillary subtype. Functional characterization revealed that CYTOR may be involved in immune-related pathways and the epithelial mesenchymal transformation (EMT) process. Moreover, CYTOR had a significant association with infiltrating immune cells, including M2 macrophages and regulatory T cells (Tregs). CYTOR facilitates the crosstalk between cancer-associated fibroblasts (CAFs) and macrophages, and mediates M2 polarization of macrophages. Correlation analysis revealed a positive correlation between CYTOR expression and programmed cell death-1 (PD-1)/programmed death ligand 1 (PD-L1)/expression and other targets for specific immunotherapy in BCa, which are recognized to predict the efficacy of immunotherapy.
UNASSIGNED: These results suggest that CYTOR serves as a potential biomarker for predicting survival outcome, TME cell infiltration characteristics and immunotherapy response in BCa.

Transcriptome analysis of head and neck squamous cell carcinoma (HNSCC) has been pivotal to comprehending the convoluted biology of HNSCC tumors. MAPKAPK2 or MK2 is a critical modulator of the mRNA turnover of crucial genes involved in HNSCC progression. However, MK2-centric transcriptome profiles of tumors are not well known. This study delves into HNSCC progression with MK2 at the nexus to delineate the biological relevance and intricate crosstalk of MK2 in the tumor milieu. We performed next-generation sequencing-based transcriptome profiling of HNSCC cells and xenograft tumors to ascertain mRNA expression profiles in MK2-wild type and MK2-knockdown conditions. The findings were validated using gene expression assays, immunohistochemistry, and transcript turnover studies. Here, we identified a pool of crucial MK2-regulated candidate genes by annotation and differential gene expression analyses. Regulatory network and pathway enrichment revealed their significance and involvement in the HNSCC pathogenesis. Additionally, 3\'-UTR-based filtering recognized important MK2-regulated downstream target genes and validated them by nCounter gene expression assays. Finally, immunohistochemistry and transcript stability studies revealed the putative role of MK2 in regulating the transcript turnover of IGFBP2, MUC4, and PRKAR2B in HNSCC. Conclusively, MK2-regulated candidate genes were identified in this study, and their plausible involvement in HNSCC pathogenesis was elucidated. These genes possess investigative values as targets for diagnosis and therapeutic interventions for HNSCC.

UNASSIGNED: To develop models based on radiomics and genomics for predicting the histopathologic nuclear grade with localized clear cell renal cell carcinoma (ccRCC) and to assess whether macro-radiomics models can predict the microscopic pathological changes.
UNASSIGNED: In this multi-institutional retrospective study, a computerized tomography (CT) radiomic model for nuclear grade prediction was developed. Utilizing a genomics analysis cohort, nuclear grade-associated gene modules were identified, and a gene model was constructed based on top 30 hub mRNA to predict the nuclear grade. Using a radiogenomic development cohort, biological pathways were enriched by hub genes and a radiogenomic map was created.
UNASSIGNED: The four-features-based SVM model predicted nuclear grade with an area under the curve (AUC) score of 0.94 in validation sets, while a five-gene-based model predicted nuclear grade with an AUC of 0.73 in the genomics analysis cohort. A total of five gene modules were identified to be associated with the nuclear grade. Radiomic features were only associated with 271 out of 603 genes in five gene modules and eight top 30 hub genes. Differences existed in the enrichment pathway between associated and un-associated with radiomic features, which were associated with two genes of five-gene signatures in the mRNA model.
UNASSIGNED: The CT radiomics models exhibited higher predictive performance than mRNA models. The association between radiomic features and mRNA related to nuclear grade is not universal.

E3 ubiquitin ligases (E3s) and deubiquitinating enzymes (DUBs) play key roles in protein degradation. However, a large number of E3 substrate interactions (ESIs) and DUB substrate interactions (DSIs) remain elusive. Here, we present DeepUSI, a deep learning-based framework to identify ESIs and DSIs using the rich information present in protein sequences. Utilizing the collected golden standard dataset, key hyperparameters in the process of model training, including the ones relevant to data sampling and number of epochs, have been systematically assessed. The performance of DeepUSI was thoroughly evaluated by multiple metrics, based on internal and external validation. Application of DeepUSI to cancer-associated E3 and DUB genes identified a list of druggable substrates with functional implications, warranting further investigation. Together, DeepUSI presents a new framework for predicting substrates of E3 ubiquitin ligases and deubiquitinates.

Head and neck squamous cell carcinoma (HNSC) is one of most common malignancies with high mortality worldwide. Importantly, the molecular heterogeneity of HNSC complicates the clinical diagnosis and treatment, leading to poor overall survival outcomes. To dissect the complex heterogeneity, recent studies have reported multiple molecular subtyping systems. For instance, HNSC can be subdivided to four distinct molecular subtypes: atypical, basal, classical, and mesenchymal, of which the mesenchymal subtype is characterized by upregulated epithelial-mesenchymal transition (EMT) and associated with poorer survival outcomes. Despite a wealth of studies into the complex molecular heterogeneity, the regulatory mechanism specific to this aggressive subtype remain largely unclear. Herein, we developed a network-based bioinformatics framework that integrates lncRNA and mRNA expression profiles to elucidate the subtype-specific regulatory mechanisms. Applying the framework to HNSC, we identified a clinically relevant lncRNA LNCOG as a key master regulator mediating EMT underlying the mesenchymal subtype. Five genes with strong prognostic values, namely ANXA5, ITGA5, CCBE1, P4HA2, and EPHX3, were predicted to be the putative targets of LNCOG and subsequently validated in other independent datasets. By integrative analysis of the miRNA expression profiles, we found that LNCOG may act as a ceRNA to sponge miR-148a-3p thereby upregulating ITGA5 to promote HNSC progression. Furthermore, our drug sensitivity analysis demonstrated that the five putative targets of LNCOG were also predictive of the sensitivities of multiple FDA-approved drugs. In summary, our bioinformatics framework facilitates the dissection of cancer subtype-specific lncRNA regulatory mechanisms, providing potential novel biomarkers for more optimized treatment of HNSC.

Brain glioma is one of the cancer types with worst prognosis, and LMO2 has been reported to play oncogenic functions in brain gliomas. Herein, analysis of datasets from The Cancer Genome Atlas (TCGA) indicated that higher LMO2 level in patient samples indicated worse prognosis in lower grade gliomas (LGG) but not glioblastoma multiforme (GBM). Further, in tumor tissues consisting of a variety of cell types, LMO2 level indicated intratumoral endothelium and pattern recognition receptor (PRR) response in both LGGs and GBMs, and additionally indicated cytotoxic T-lymphocyte, M2 macrophage infiltration and fibroblast specifically in LGGs. Moreover, only in LGGs these aspects were significantly associated with patient survival, in either risky or protective manner, and these dissected associations can give a better prediction on patient prognosis than LMO2 alone. This study not only provided more detailed understandings of LMO2 functional representatives in brain gliomas but also demonstrated that dealing with certain gene (LMO2 in this study) in transcriptome data with the Weighted Gene Co-Expression Network Analysis (WGCNA) method was a robust strategy for dissecting exact and reasonable gene functions/associations in a complicated tumor environment.