Integrated Genomic Characterization of the Human Immunome in Cancer

Quantified pathway mutations associate epithelial-mesenchymal transition and immune escape with poor prognosis and immunotherapy resistance of head and neck squamous cell carcinoma

BACKGROUND

Pathway mutations have been calculated to predict the poor prognosis and immunotherapy resistance in head and neck squamous cell carcinoma (HNSCC). To uncover the unique markers predicting prognosis and immune therapy response, the accurate quantification of pathway mutations are required to evaluate epithelial-mesenchymal transition (EMT) and immune escape. Yet, there is a lack of score to accurately quantify pathway mutations.

MATERIAL AND METHODS

Firstly, we proposed Individualized Weighted Hallmark Gene Set Mutation Burden (IWHMB, https://github.com/YuHongHuang-lab/IWHMB ) which integrated pathway structure information and eliminated the interference of global Tumor Mutation Burden to accurately quantify pathway mutations. Subsequently, to further elucidate the association of IWHMB with EMT and immune escape, support vector machine regression model was used to identify IWHMB-related transcriptomic features (IRG), while Adversarially Regularized Graph Autoencoder (ARVGA) was used to further resolve IRG network features. Finally, Random walk with restart algorithm was used to identify biomarkers for predicting ICI response.

RESULTS

We quantified the HNSCC pathway mutation signatures and identified pathway mutation subtypes using IWHMB. The IWHMB-related transcriptomic features (IRG) identified by support vector machine regression were divided into 5 communities by ARVGA, among which the Community 1 enriching malignant mesenchymal components promoted EMT dynamically and regulated immune patterns associated with ICI responses. Bridge Hub Gene (BHG) identified by random walk with restart was key to IWHMB in EMT and immune escape, thus, more predictive for ICI response than other 70 public signatures.

CONCLUSION

In summary, the novel pathway mutation scoring-IWHMB suggested that the elevated malignancy mediated by pathway mutations is a major cause of poor prognosis and immunotherapy failure in HNSCC, and is capable of identifying novel biomarkers to predict immunotherapy response.

A Comprehensive Benchmark of Transcriptomic Biomarkers for Immune Checkpoint Blockades

Immune checkpoint blockades (ICBs) have revolutionized cancer therapy by inducing durable clinical responses, but only a small percentage of patients can benefit from ICB treatments. Many studies have established various biomarkers to predict ICB responses. However, different biomarkers were found with diverse performances in practice, and a timely and unbiased assessment has yet to be conducted due to the complexity of ICB-related studies and trials. In this study, we manually curated 29 published datasets with matched transcriptome and clinical data from more than 1400 patients, and uniformly preprocessed these datasets for further analyses. In addition, we collected 39 sets of transcriptomic biomarkers, and based on the nature of the corresponding computational methods, we categorized them into the gene-set-like group (with the self-contained design and the competitive design, respectively) and the deconvolution-like group. Next, we investigated the correlations and patterns of these biomarkers and utilized a standardized workflow to systematically evaluate their performance in predicting ICB responses and survival statuses across different datasets, cancer types, antibodies, biopsy times, and combinatory treatments. In our benchmark, most biomarkers showed poor performance in terms of stability and robustness across different datasets. Two scores (TIDE and CYT) had a competitive performance for ICB response prediction, and two others (PASS-ON and EIGS_ssGSEA) showed the best association with clinical outcome. Finally, we developed ICB-Portal to host the datasets, biomarkers, and benchmark results and to implement the computational methods for researchers to test their custom biomarkers. Our work provided valuable resources and a one-stop solution to facilitate ICB-related research.

Transcriptomic analyses of patient peripheral blood with hemoglobin depletion reveal glioblastoma biomarkers

Peripheral blood is gaining prominence as a noninvasive alternative to tissue biopsy to develop biomarkers for glioblastoma (GBM); however, widely utilized blood-based biomarkers in clinical settings have not yet been identified due to the lack of a robust detection approach. Here, we describe the application of globin reduction in RNA sequencing of whole blood (i.e., WBGR) and perform transcriptomic analysis to identify GBM-associated transcriptomic changes. By using WBGR, we improved the detection sensitivity of informatic reads and identified differential gene expression in GBM blood. By analyzing tumor tissues, we identified transcriptomic traits of GBM blood. Further functional enrichment analyses retained the most changed genes in GBM. Subsequent validation elicited a 10-gene panel covering mRNA, long noncoding RNA, and microRNA (i.e., GBM-Dx panel) that has translational potential to aid in the early detection or clinical management of GBM. Here, we report an integrated approach, WBGR, with comprehensive analytic capacity for blood-based marker identification.

Immune Pathways with Aging Characteristics Improve Immunotherapy Benefits and Drug Prediction in Human Cancer

(1) Background: Perturbation of immune-related pathways can make substantial contributions to cancer. However, whether and how the aging process affects immune-related pathways during tumorigenesis remains largely unexplored. (2) Methods: Here, we comprehensively investigated the immune-related genes and pathways among 25 cancer types using genomic and transcriptomic data. (3) Results: We identified several pathways that showed aging-related characteristics in various cancers, further validated by conventional aging-related gene sets. Genomic analysis revealed high mutation burdens in cytokines and cytokines receptors pathways, which were strongly correlated with aging in diverse cancers. Moreover, immune-related pathways were found to be favorable prognostic factors in melanoma. Furthermore, the expression level of these pathways had close associations with patient response to immune checkpoint blockade therapy in melanoma and non-small cell lung cancer. Applying a net-work-based method, we predicted immune- and aging-related genes in pan-cancer and utilized these genes for potential immunotherapy drug discovery. Mapping drug target data to our top-ranked genes identified potential drug targets, FYN, JUN, and SRC. (4) Conclusions: Taken together, our systematic study helped interpret the associations among immune-related pathways, aging, and cancer and could serve as a resource for promoting clinical treatment.

Analysis of Tumor Microenvironment Heterogeneity among Breast Cancer Subtypes to Identify Subtype-Specific Signatures

Breast cancer is one of the most frequent malignancies in women worldwide. According to 50-gene signature, Prediction Analysis of Microarray 50 (PAM50), breast cancer can be categorized into five molecular subtypes, and these subtypes are highly heterogeneous in different molecular characteristics. However, the landscape of their tumor microenvironment (TME) heterogeneity has not been fully researched. Using the multi-omics dataset of breast cancer from the METABRIC cohort (n = 1699), we conducted extensive analyses of TME-related features to investigate TME heterogeneity in each breast cancer subtype. We then developed a cell-based subtype set enrichment analysis to identify the subtype-specific TME cells, and further evaluate their prognostic effects. Our results illustrate that different breast cancer subtypes exhibit different TME patterns. Basal-like and HER2-enriched subtypes are associated with high immune scores, expression of most immune regulatory targets, and immune cell infiltration, suggesting that these subtypes could be defined as "immune hot" tumors and suitable for immune checkpoint blockade (ICB) therapy. In contrast, Luminal A and Luminal B subtypes are associated with low immune scores and immune cell infiltration, suggesting that these subtypes could be defined as "immune cold" tumors. Additionally, the Normal-like subtype has relatively high levels of both immune and stromal features, which indicates that the Normal-like subtype may be suitable for more diverse treatment strategies. Our study reveals the breast cancer tumor microenvironment heterogeneity across subtypes. The comprehensive analysis of breast cancer TME-related characteristics may help us to adopt a tailored treatment strategy for different subtypes of patients.

Combining Single-Cell and Transcriptomic Data Revealed the Prognostic Significance of Glycolysis in Pancreatic Cancer

Background: Pancreatic cancer (PC), the most common fatal solid malignancy, has a very dismal prognosis. Clinical computerized tomography (CT) and pathological TNM staging are no longer sufficient for determining a patient’s prognosis. Although numerous studies have suggested that glycolysis is important in the onset and progression of cancer, there are few publications on its impact on PC.

Methods: To begin, the single-sample gene set enrichment analysis (ssGSEA) approach was used to quantify the glycolysis pathway enrichment fraction in PC patients and establish its prognostic significance. The genes most related to the glycolytic pathway were then identified using weighted gene co-expression network analysis (WGCNA). The glycolysis-associated prognostic signature in PC patients was then constructed using univariate Cox regression and lasso regression methods, which were validated in numerous external validation cohorts. Furthermore, we investigated the activation of the glycolysis pathway in PC cell subtypes at the single-cell level, performed a quasi-time series analysis on the activated cell subtypes and then detected gene changes in the signature during cell development. Finally, we constructed a decision tree and a nomogram that could divide the patients into different risk subtypes, according to the signature score and their different clinical characteristics and assessed the prognosis of PC patients.

Results: Glycolysis plays a risky role in PC patients. Our glycolysis-related signature could effectively discriminate the high-risk and low-risk patients in both the trained cohort and the independent externally validated cohort. The survival analysis and multivariate Cox analysis indicated this gene signature to be an independent prognostic factor in PC. The prognostic ROC curve analysis suggested a high accuracy of this gene signature in predicting the patient prognosis in PC. The single-cell analysis suggested that the glycolytic pathway may be more activated in epithelial cells and that the genes in the signature were also mainly expressed in epithelial cells. The decision tree analysis could effectively identify patients in different risk subgroups, and the nomograms clearly show the prognostic assessment of PC patients.

Conclusion: Our study developed a glycolysis-related signature, which contributes to the risk subtype assessment of patients with PC and to the individualized management of patients in the clinical setting.

Pediatric Pan-Central Nervous System Tumor Methylome Analyses Reveal Immune-Related LncRNAs

Pediatric central nervous system (CNS) tumors are the second most common cancer diagnosis among children. Long noncoding RNAs (lncRNAs) emerge as critical regulators of gene expression, and they play fundamental roles in immune regulation. However, knowledge on epigenetic changes in lncRNAs in diverse types of pediatric CNS tumors is lacking. Here, we integrated the DNA methylation profiles of 2,257 pediatric CNS tumors across 61 subtypes with lncRNA annotations and presented the epigenetically regulated landscape of lncRNAs. We revealed the prevalent lncRNA methylation heterogeneity across pediatric pan-CNS tumors. Based on lncRNA methylation profiles, we refined 14 lncRNA methylation clusters with distinct immune microenvironment patterns. Moreover, we found that lncRNA methylations were significantly correlated with immune cell infiltrations in diverse tumor subtypes. Immune-related lncRNAs were further identified by investigating their correlation with immune cell infiltrations and potentially regulated target genes. LncRNA with methylation perturbations potentially regulate the genes in immune-related pathways. We finally identified several candidate immune-related lncRNA biomarkers (i.e., SSTR5-AS1, CNTN4-AS1, and OSTM1-AS1) in pediatric cancer for further functional validation. In summary, our study represents a comprehensive repertoire of epigenetically regulated immune-related lncRNAs in pediatric pan-CNS tumors, and will facilitate the development of immunotherapeutic targets.

Identification of immune infiltration-related genes as prognostic indicators for hepatocellular carcinoma

Hepatocellular carcinoma (HCC) has a high degree of malignancy and a poor prognosis. Immune infiltration-related genes have shown good predictive value in the prognosis of many solid tumours. In this study, we established and verified prognostic biomarkers consisting of immune infiltration-related genes in HCC. Gene expression data and clinical data were downloaded from The Cancer Genome Atlas (TCGA) database. Differential gene expression analysis, univariate Cox regression analysis and the least absolute shrinkage and selection operator (LASSO) regression algorithm were used to screen prognostic immune infiltration-related genes and to construct a risk scoring model. Kaplan-Meier (KM) survival plots and receiver operating characteristic (ROC) curve analysis were used to evaluate the prognostic performance of the risk scoring model in the TCGA-HCC cohort. In addition, a nomogram model with a risk score was established, and its predictive performance was verified by ROC analysis and calibration plot analysis in the TCGA-HCC cohort. Gene set enrichment analysis (GSEA) identified pathways and biological processes that may be enriched in the high-risk group. Finally, immune infiltration analysis was used to explore the characteristics of the tumour microenvironment related to the risk score. We identified 17 immune infiltration-related genes with prognostic value and constructed a risk scoring model. ROC analysis showed that the risk scoring model can accurately predict the 1-year, 3-year, and 5-year overall survival (OS) of HCC patients in the TCGA-HCC cohort. KM analysis showed that the OS of the high-risk group was significantly lower than that of the low-risk group (P < 0.001). The nomogram model effectively predicted the OS of HCC patients in the TCGA-HCC cohort. GSEA indicated that the immune infiltration-related genes may be involved in biological processes such as amino acid and lipid metabolism, matrisome and small molecule transportation, immune system regulation, and hepatitis virus infection. Immune infiltration analysis showed that the level of immune cell infiltration in the high-risk group was low, and the risk score was negatively correlated with infiltrating immune cells. Our prognostic model based on immune infiltration-related genes in HCC could help the prognostic assessment of HCC patients and provide potential targets for HCC inhibition.

Pan-cancer analysis of the DNA methylation patterns of long non-coding RNA

Long non-coding RNA (lncRNA) regulated by abnormal DNA methylation (ADM-lncRNA) emerges as a biomarker for cancer diagnosis and treatment. This study comprehensively described the methylation patterns of lncRNA in pan-cancer using the cancer data set in The Cancer Genome Atlas (TCGA). Based on the cancer heterogeneity of ADM-lncRNA in pan-cancer, we constructed a co-expression network of pan-cancer ADM-lncRNA (pADM-lncRNA) in 10 cancers, highlighting the combined action mode of abnormal DNA methylation, and indicating the internal connection among different cancers. Functional analysis revealed the pan-carcinogenic pathway of pADM-lncRNA and suggested potential factors for cancer heterogeneity and tumor immune microenvironment changes. Survival analysis showed the potential of pADM-lncRNA-mRNA co-expression pair as cancer biomarkers. Revealing the action mode of lncRNA and DNA methylation in cancer may help understand the key molecular mechanisms of cell carcinogenesis.

Integrative analysis of CRISPR screening data uncovers new opportunities for optimizing cancer immunotherapy

Background

In recent years, the application of functional genetic immuno-oncology screens has showcased the striking ability to identify potential regulators engaged in tumor-immune interactions. Although these screens have yielded substantial data, few studies have attempted to systematically aggregate and analyze them.

Methods

In this study, a comprehensive data collection of tumor immunity-associated functional screens was performed. Large-scale genomic data sets were exploited to conduct integrative analyses.

Results

We identified 105 regulator genes that could mediate resistance or sensitivity to immune cell-induced tumor elimination. Further analysis identified MON2 as a novel immune-oncology target with considerable therapeutic potential. In addition, based on the 105 genes, a signature named CTIS (CRISPR screening-based tumor-intrinsic immune score) for predicting response to immune checkpoint blockade (ICB) and several immunomodulatory agents with the potential to augment the efficacy of ICB were also determined.

Conclusion

Overall, our findings provide insights into immune oncology and open up novel opportunities for improving the efficacy of current immunotherapy agents.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12943-021-01462-z.

Significance of TP53 mutation in treatment and prognosis in head and neck squamous cell carcinoma

Background: TP53 is ranked as the most common mutated gene in head and neck squamous cell carcinoma (HNSCC). Results: The status of TP53 mutation was investigated on International Cancer Genome Consortium and The Cancer Genome Atlas database and TP53-related differentially expressed genes were screened out from transcriptome data and mutation information. A TP53-related prognostic gene signature (TIMP4, ONECUT2, CGNL1, DMRTA2 and NKX2.3) was constructed based on Cox regression analysis and LASSO algorithm. Univariate and multivariate analyses were carried out to identify promising prognosticators for HNSCC. Conclusion: Our findings provide a well-rounded landscape of TP53 mutation in HNSCC and pave the groundwork for developing innovative and effective cancer treatment methods for HNSCC.