A Predictive Network-Based Immune Checkpoint Blockade Immunotherapeutic Signature Optimizing Patient Selection and Treatment Strategies.

Immune checkpoint blockade (ICB) therapy has brought significant advancements to the field of oncology. However, the diverse responses among patients highlight the need for more accurate predictive tools. In this study, insights are drawn from tumor-immunology pathways, and a novel network-based ICB immunotherapeutic signature, termed ICBnetIS, is constructed. The signature is derived from advanced biological network-based computational strategies involving co-expression networks and molecular interactions networks. The efficacy of ICBnetIS is established through its association with enhanced patient survival and a robust immune response characterized by diverse immune cell infiltration and active anti-tumor immune pathways. The validation process positions ICBnetIS as an effective tool in predicting responses to ICB therapy, analyzing ICB data from a broad collection of over 700 samples from multiple cancer types of more than 15 datasets. It achieves an aggregated prediction AUC of 0.784, which outperforms the other nine renowned immunotherapeutic signatures, indicating the superior predictive capability of ICBnetIS. To sum up, the findings suggest ICBnetIS as a potent tool in predicting ICB therapy responses, offering significant implications for patient selection and treatment optimization in oncology. The study highlights the role of ICBnetIS in advancing personalized treatment strategies, potentially transforming the clinical landscape of ICB therapy.

A method for predicting drugs that can boost the efficacy of immune checkpoint blockade.

Combination therapy is a promising therapeutic strategy to enhance the efficacy of immune checkpoint blockade (ICB); however, predicting drugs for effective combination is challenging. Here we developed a general data-driven method called CM-Drug for screening compounds that can boost ICB treatment efficacy based on core and minor gene sets identified between responsive and nonresponsive samples in ICB therapy. The CM-Drug method was validated using melanoma and lung cancer mouse models, with combined therapeutic efficacy demonstrated in eight of nine predicted compounds. Among these compounds, taltirelin had the strongest synergistic effect. Mechanistic analysis and experimental verification demonstrated that taltirelin can stimulate CD8+ T cells and is mediated by the induction of thyroid-stimulating hormone. This study provides an effective and general method for predicting and evaluating drugs for combination therapy and identifies candidate compounds for future ICB combination therapy.

Optimization of FFPE preparation and identification of gene attributes associated with RNA degradation.

Formalin-fixed paraffin-embedded (FFPE) tissues are widely available specimens for clinical studies. However, RNA degradation in FFPE tissues often restricts their utility. In this study, we determined optimal FFPE preparation conditions, including tissue ischemia at 4°C (<48 h) or 25°C for a short time (0.5 h), 48-h fixation at 25°C and sampling from FFPE scrolls instead of sections. Notably, we observed an increase in intronic reads and a significant change in gene rank based on expression level in the FFPE as opposed to fresh-frozen (FF) samples. Additionally, we found that more reads were mapped to genes associated with chemical stimulus in FFPE samples. Furthermore, we demonstrated that more degraded genes in FFPE samples were enriched in genes with short transcripts and high free energy. Besides, we found 40 housekeeping genes exhibited stable expression in FF and FFPE samples across various tissues. Moreover, our study showed that FFPE samples yielded comparable results to FF samples in dimensionality reduction and pathway analyses between case and control samples. Our study established the optimal conditions for FFPE preparation and identified gene attributes associated with degradation, which would provide useful clues for the utility of FFPE tissues in clinical practice and research.

Identification of P21 (CDKN1A) Activated Kinase 4 as a Susceptibility Gene for Familial Non-Medullary Thyroid Carcinoma.

Background: Familial non-medullary thyroid carcinoma (FNMTC) is a genetically predisposed disease with unclear genetic mechanisms. This makes research on susceptibility genes important for the diagnosis and treatment options. Methods: This study included a five-member family affected by papillary thyroid carcinoma. The candidate genes were identified through whole-exome sequencing and Sanger sequencing in family members, other FNMTC patients, and sporadic non-medullary thyroid carcinoma patients. The pathogenicity of the mutation was predicted using in silico tools. Cell phenotype experiments in vitro and models of lung distant metastasis in vivo were conducted to confirm the characteristics of the mutation. Transcriptome sequencing and mechanistic validation were employed to compare the disparities between PAK4 wild-type (WT) and PAK4 mutant (MUT) cell lines. Results: This mutation alters the protein structure, potentially increasing instability by affecting hydrophobicity, intra-molecular hydrogen bonding, and phosphorylation sites. It specifically promotes phosphorylated PAK4 nuclear translocation and expression in thyroid tissue and cell lines. Compared with the WT cells line, PAK4 I417T demonstrates enhanced proliferation, invasiveness, accelerated cell division, and inhibition of cell apoptosis in vitro. In addition, it exhibits a significant propensity for metastasis in vivo. It activates tumor necrosis factor signaling through increased phosphorylation of PAK4, JNK, NFκB, and c-Jun, unlike the WT that activates it via the PAK4-NFκ-MMP9 axis. In addition, PAK4 MUT protein interacts with matrix metalloproteinase (MMP)3 and regulates MMP3 promoter activity, which is not observed in the WT. Conclusions: Our study identified PAK4: c.T1250C: p.I417T as a potential susceptibility gene for FNMTC. The study concludes that the mutant form of PAK4 exhibits oncogenic function, suggesting its potential as a novel diagnostic molecular marker for FNMTC.

ncRNADrug: a database for validated and predicted ncRNAs associated with drug resistance and targeted by drugs.

Drug resistance is a major barrier in cancer treatment and anticancer drug development. Growing evidence indicates that non-coding RNAs (ncRNAs), especially microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs), play pivotal roles in cancer progression, therapy, and drug resistance. Furthermore, ncRNAs have been proven to be promising novel therapeutic targets for cancer treatment. Reversing dysregulated ncRNAs by drugs holds significant potential as an effective therapeutic strategy for overcoming drug resistance. Therefore, we developed ncRNADrug, an integrated and comprehensive resource that records manually curated and computationally predicted ncRNAs associated with drug resistance, ncRNAs targeted by drugs, as well as potential drug combinations for the treatment of resistant cancer. Currently, ncRNADrug collects 29 551 experimentally validated entries involving 9195 ncRNAs (2248 miRNAs, 4145 lncRNAs and 2802 circRNAs) associated with the drug resistance of 266 drugs, and 32 969 entries involving 10 480 ncRNAs (4338 miRNAs, 6087 lncRNAs and 55 circRNAs) targeted by 965 drugs. In addition, ncRNADrug also contains associations between ncRNAs and drugs predicted from ncRNA expression profiles by differential expression analysis. Altogether, ncRNADrug surpasses the existing related databases in both data volume and functionality. It will be a useful resource for drug development and cancer treatment. ncRNADrug is available at http://www.jianglab.cn/ncRNADrug.

ICBcomb: a comprehensive expression database for immune checkpoint blockade combination therapy.

The success of immune checkpoint blockade (ICB) promotes the immunotherapy to be a new pillar in cancer treatment. However, the low response rate of the ICB therapy limits its application. To increase the response rate and enhance efficacy, the ICB combination therapy has emerged and its clinical trials are increasing. Nevertheless, the gene expression profile and its pattern of ICB combination were not comprehensively studied, which limits the understanding of the ICB combination therapy and the identification of new drugs. Here, we constructed ICBcomb (http://bioinfo.life.hust.edu.cn/ICBcomb/), a comprehensive database, by analyzing the human and mouse expression data of the ICB combination therapy and comparing them between groups treated with ICB, other drugs or their combinations. ICBcomb contains 1399 samples across 29 cancer types involving 52 drugs. It provides a user-friendly web interface for demonstrating the results of the available comparisons in the ICB combination therapy datasets with five functional modules: [1, 2] the 'Dataset/Disease' modules for browsing the expression, enrichment and comparison results in each dataset or disease; [3] the 'Gene' module for inputting a gene symbol and displaying its expression and comparison results across datasets/diseases; [4] the 'Gene Set' module for GSVA/GSEA enrichment analysis on the built-in gene sets and the user-input gene sets in different comparisons; [5] the 'Immune Cell' module for immune cell infiltration comparison between different groups by immune cell abundance analysis. The ICBcomb database provides the first resource for gene expression profile and comparison in ICB combination therapy, which may provide clues for discovering the mechanism of effective combination strategies and new combinatory drugs.

Revealing the impact of CD70 expression on the manufacture and functions of CAR-70 T-cells based on single-cell transcriptomics.

BACKGROUND: Chimeric antigen receptor-modified T cells (CAR T-cells) have shown exhilarative clinical efficacy for hematological malignancies. However, a shared antigen pool between healthy and malignant T-cells remains a concept to be technically and clinically explored for CAR T-cell therapy in T-cell cancers. No guidelines for engineering CAR T-cells targeting self-expressed antigens are currently available. METHOD: Based on anti-CD70 CAR (CAR-70) T-cells, we constructed CD70 knock-out and wild-type CAR (CAR-70KO and CAR-70WT) T-cells and evaluated their manufacturing and anti-tumor capability. Single-cell RNA sequencing and TCR sequencing were performed to further reveal the underlying differences between the two groups of CAR T-cells. RESULTS: Our data showed that the disruption of target genes in T-cells before CAR transduction advantaged the expansion and cell viability of CAR T-cells during manufacturing periods, as well as the degranulation, anti-tumor efficacy, and proliferation potency in response to tumor cells. Meanwhile, more naïve and central memory phenotype CAR+ T-cells, with higher TCR clonal diversity, remained in the final products in KO samples. Gene expression profiles revealed a higher activation and exhaustion level of CAR-70WT T-cells, while signaling transduction pathway analysis identified a higher level of the phosphorylation-related pathway in CAR-70KO T-cells. CONCLUSION: This study evidenced that CD70 stimulation during manufacturing process induced early exhaustion of CAR-70 T-cells. Knocking-out CD70 in T-cells prevented the exhaustion and led to a better-quality CAR-70 T-cell product. Our research will contribute to good engineering CAR T-cells targeting self-expressed antigens.

scLiverDB: a Database of Human and Mouse Liver Transcriptome Landscapes at Single-Cell Resolution.

The liver is critical for the digestive and immune systems. Although the physiology and pathology of liver have been well studied and many scRNA-seq data are generated, a database and landscape for characterizing cell types and gene expression in different liver diseases or developmental stages at single-cell resolution are lacking. Hence, scLiverDB is developed, a specialized database for human and mouse liver transcriptomes to unravel the landscape of liver cell types, cell heterogeneity and gene expression at single-cell resolution across various liver diseases/cell types/developmental stages. To date, 62 datasets including 9,050 samples and 1,741,734 cells is curated. A uniform workflow is used, which included quality control, dimensional reduction, clustering, and cell-type annotation to analyze datasets on the same platform; integrated manual and automatic methods for accurate cell-type identification and provided a user-friendly web interface with multiscale functions. There are two case studies to show the usefulness of scLiverDB, which identified the LTB (lymphotoxin Beta) gene as a potential biomarker of lymphoid cells differentiation and showed the expression changes of Foxa3 (forkhead box A3) in liver chronic progressive diseases. This work provides a crucial resource to resolve molecular and cellular information in normal, diseased, and developing human and mouse livers.

Small-sample learning reveals propionylation in determining global protein homeostasis.

Proteostasis is fundamental for maintaining organismal health. However, the mechanisms underlying its dynamic regulation and how its disruptions lead to diseases are largely unclear. Here, we conduct in-depth propionylomic profiling in Drosophila, and develop a small-sample learning framework to prioritize the propionylation at lysine 17 of H2B (H2BK17pr) to be functionally important. Mutating H2BK17 which eliminates propionylation leads to elevated total protein level in vivo. Further analyses reveal that H2BK17pr modulates the expression of 14.7-16.3% of genes in the proteostasis network, and determines global protein level by regulating the expression of genes involved in the ubiquitin-proteasome system. In addition, H2BK17pr exhibits daily oscillation, mediating the influences of feeding/fasting cycles to drive rhythmic expression of proteasomal genes. Our study not only reveals a role of lysine propionylation in regulating proteostasis, but also implements a generally applicable method which can be extended to other issues with little prior knowledge.

Platelet RNA enables accurate detection of ovarian cancer: an intercontinental, biomarker identification study.

Platelets are reprogrammed by cancer via a process called education, which favors cancer development. The transcriptional profile of tumor-educated platelets (TEPs) is skewed and therefore practicable for cancer detection. This intercontinental, hospital-based, diagnostic study included 761 treatment-naïve inpatients with histologically confirmed adnexal masses and 167 healthy controls from nine medical centers (China, n = 3; Netherlands, n = 5; Poland, n = 1) between September 2016 and May 2019. The main outcomes were the performance of TEPs and their combination with CA125 in two Chinese (VC1 and VC2) and the European (VC3) validation cohorts collectively and independently. Exploratory outcome was the value of TEPs in public pan-cancer platelet transcriptome datasets. The AUCs for TEPs in the combined validation cohort, VC1, VC2, and VC3 were 0.918 (95% CI 0.889-0.948), 0.923 (0.855-0.990), 0.918 (0.872-0.963), and 0.887 (0.813-0.960), respectively. Combination of TEPs and CA125 demonstrated an AUC of 0.922 (0.889-0.955) in the combined validation cohort; 0.955 (0.912-0.997) in VC1; 0.939 (0.901-0.977) in VC2; 0.917 (0.824-1.000) in VC3. For subgroup analysis, TEPs exhibited an AUC of 0.858, 0.859, and 0.920 to detect early-stage, borderline, non-epithelial diseases and 0.899 to discriminate ovarian cancer from endometriosis. TEPs had robustness, compatibility, and universality for preoperative diagnosis of ovarian cancer since it withstood validations in populations of different ethnicities, heterogeneous histological subtypes, and early-stage ovarian cancer. However, these observations warrant prospective validations in a larger population before clinical utilities.

A comprehensive landscape of transcription profiles and data resources for human leukemia.

As a heterogeneous group of hematologic malignancies, leukemia has been widely studied at the transcriptome level. However, a comprehensive transcriptomic landscape and resources for different leukemia subtypes are lacking. Thus, in this study, we integrated the RNA sequencing data sets of >3000 samples from 14 leukemia subtypes and 53 related cell lines via a unified analysis pipeline. We depicted the corresponding transcriptomic landscape and developed a user-friendly data portal LeukemiaDB. LeukemiaDB was designed with 5 main modules: protein-coding gene, long noncoding RNA (lncRNA), circular RNA, alternative splicing, and fusion gene modules. In LeukemiaDB, users can search and browse the expression level, regulatory modules, and molecular information across leukemia subtypes or cell lines. In addition, a comprehensive analysis of data in LeukemiaDB demonstrates that (1) different leukemia subtypes or cell lines have similar expression distribution of the protein-coding gene and lncRNA; (2) some alternative splicing events are shared among nearly all leukemia subtypes, for example, MYL6 in A3SS, MYB in A5SS, HMBS in retained intron, GTPBP10 in mutually exclusive exons, and POLL in skipped exon; (3) some leukemia-specific protein-coding genes, for example, ABCA6, ARHGAP44, WNT3, and BLACE, and fusion genes, for example, BCR-ABL1 and KMT2A-AFF1 are involved in leukemogenesis; (4) some highly correlated regulatory modules were also identified in different leukemia subtypes, for example, the HOXA9 module in acute myeloid leukemia and the NOTCH1 module in T-cell acute lymphoblastic leukemia. In summary, the developed LeukemiaDB provides valuable insights into oncogenesis and progression of leukemia and, to the best of our knowledge, is the most comprehensive transcriptome resource of human leukemia available to the research community.

GSCA: an integrated platform for gene set cancer analysis at genomic, pharmacogenomic and immunogenomic levels.

Cancer initiation and progression are likely caused by the dysregulation of biological pathways. Gene set analysis (GSA) could improve the signal-to-noise ratio and identify potential biological insights on the gene set level. However, platforms exploring cancer multi-omics data using GSA methods are lacking. In this study, we upgraded our GSCALite to GSCA (gene set cancer analysis, http://bioinfo.life.hust.edu.cn/GSCA) for cancer GSA at genomic, pharmacogenomic and immunogenomic levels. In this improved GSCA, we integrated expression, mutation, drug sensitivity and clinical data from four public data sources for 33 cancer types. We introduced useful features to GSCA, including associations between immune infiltration with gene expression and genomic variations, and associations between gene set expression/mutation and clinical outcomes. GSCA has four main functional modules for cancer GSA to explore, analyze and visualize expression, genomic variations, tumor immune infiltration, drug sensitivity and their associations with clinical outcomes. We used case studies of three gene sets: (i) seven cell cycle genes, (ii) tumor suppressor genes of PI3K pathway and (iii) oncogenes of PI3K pathway to prove the advantage of GSCA over single gene analysis. We found novel associations of gene set expression and mutation with clinical outcomes in different cancer types on gene set level, while on single gene analysis level, they are not significant associations. In conclusion, GSCA is a user-friendly web server and a useful resource for conducting hypothesis tests by using GSA methods at genomic, pharmacogenomic and immunogenomic levels.

AnimalTFDB 4.0: a comprehensive animal transcription factor database updated with variation and expression annotations.

Transcription factors (TFs) are proteins that interact with specific DNA sequences to regulate gene expression and play crucial roles in all kinds of biological processes. To keep up with new data and provide a more comprehensive resource for TF research, we updated the Animal Transcription Factor Database (AnimalTFDB) to version 4.0 (http://bioinfo.life.hust.edu.cn/AnimalTFDB4/) with up-to-date data and functions. We refined the TF family rules and prediction pipeline to predict TFs in genome-wide protein sequences from Ensembl. As a result, we predicted 274 633 TF genes and 150 726 transcription cofactor genes in AnimalTFDB 4.0 in 183 animal genomes, which are 86 more species than AnimalTFDB 3.0. Besides double data volume, we also added the following new annotations and functions to the database: (i) variations (including mutations) on TF genes in various human cancers and other diseases; (ii) predicted post-translational modification sites (including phosphorylation, acetylation, methylation and ubiquitination sites) on TFs in 8 species; (iii) TF regulation in autophagy; (iv) comprehensive TF expression annotation for 38 species; (v) exact and batch search functions allow users to search AnimalTFDB flexibly. AnimalTFDB 4.0 is a useful resource for studying TF and transcription regulation, which contains comprehensive annotation and classification of TFs and transcription cofactors.

lncRNASNP v3: an updated database for functional variants in long non-coding RNAs.

Long non-coding RNAs (lncRNAs) act as versatile regulators of many biological processes and play vital roles in various diseases. lncRNASNP is dedicated to providing a comprehensive repository of single nucleotide polymorphisms (SNPs) and somatic mutations in lncRNAs and their impacts on lncRNA structure and function. Since the last release in 2018, there has been a huge increase in the number of variants and lncRNAs. Thus, we updated the lncRNASNP to version 3 by expanding the species to eight eukaryotic species (human, chimpanzee, pig, mouse, rat, chicken, zebrafish, and fruitfly), updating the data and adding several new features. SNPs in lncRNASNP have increased from 11 181 387 to 67 513 785. The human mutations have increased from 1 174 768 to 2 387 685, including 1 031 639 TCGA mutations and 1 356 046 CosmicNCVs. Compared with the last release, updated and new features in lncRNASNP v3 include (i) SNPs in lncRNAs and their impacts on lncRNAs for eight species, (ii) SNP effects on miRNA-lncRNA interactions for eight species, (iii) lncRNA expression profiles for six species, (iv) disease & GWAS-associated lncRNAs and variants, (v) experimental & predicted lncRNAs and drug target associations and (vi) SNP effects on lncRNA expression (eQTL) across tumor & normal tissues. The lncRNASNP v3 is freely available at http://gong_lab.hzau.edu.cn/lncRNASNP3/.

Regulation of IFN-γ-mediated PD-L1 expression by MYC in colorectal cancer with wild-type KRAS and TP53 and its clinical implications.

Introduction: In the tumor microenvironment, interferon gamma (IFN-γ) secreted by tumor infiltrating lymphocytes can upregulate programmed cell death 1 ligand 1 (PD-L1) expression in many cancers. The present study evaluated the expression of PD-L1 in selected colorectal cancer cell lines with IFN-γ treatment and explored the correlation between programmed cell death 1 ligand 1 expression and KRAS/TP53 mutation status. Methods: The selected colorectal cancer cell lines had known KRAS mutations or TP53 mutations. TCGA data analysis were used to investigate the correlation between overall survival of patient with anti-PD-1/PD-L1 immunotherapy and KRAS/TP53 mutation status. Besides, the correlation between PD-L1 expression and KRAS/TP53 mutation status were also investigated by using TCGA data analysis. In vitro experiments were used to explore the mechanism underlying KRAS- and TP53-related PD-L1 expression. Results: Firstly, TCGA data analysis for gene expression and overall survival and an in vitro study revealed that the wild-type KRAS/TP53 cell lines exhibited hyperresponsiveness to interferon gamma exposure and correlated with better survival in patients receiving anti-PD-1/PD-L1 treatment. Secondly, experimental data revealed that interferon gamma induced the upregulation of programmed cell death 1 ligand 1 mainly through regulating MYC in wild-type KRAS and TP53 colorectal cancers. Discussion: Our findings revealed that the response to anti-PD-1/PD-L1 cancer immunotherapy frequently happened in wild-type KRAS and TP53 colorectal cancers, which were also found to show higher programmed cell death 1 ligand 1 expression. Our results indicate that the wild-type KRAS/TP53 colorectal cancer cell lines may respond better to interferon gamma treatment, which causes increased programmed cell death 1 ligand 1 expression and may be a mechanism underlying the better responses to anti-PD-1/PD-L1 therapies in wild-type KRAS and wild-type TP53 colorectal cancer. Furthermore, the experimental results suggest that interferon gamma regulated programmed cell death 1 ligand 1 expression through the regulation of MYC, which may further affect the response to PD-1/PD-L1 cancer immunotherapy. These results suggest a novel potential treatment strategy for enhancing the efficacy of PD-1/PD-L1 blockade immunotherapy in most colorectal cancer patients.

Single-Cell Transcriptomics Reveals Longevity Immune Remodeling Features Shared by Centenarians and Their Offspring.

Centenarians, who show mild infections and low incidence of tumors, are the optimal model to investigate healthy aging. However, longevity related immune characteristics has not been fully revealed largely due to lack of appropriate controls. In this study, single-cell transcriptomic analysis of peripheral blood mononuclear cells (PBMCs) derived from seven centenarians (CEN), six centenarians' offspring (CO), and nine offspring spouses or neighbors (Control, age-matched to CO) are performed to investigate the shared immune features between CEN and CO. The results indicate that among all 12 T cell clusters, the cytotoxic-phenotype-clusters (CPC) and the naïve-phenotype-clusters (NPC) significantly change between CEN and ontrol. Compared to Control, both CEN and CO are characterized by depleted NPC and increased CPC, which is dominated by CD8+ T cells. Furthermore, CPC from CEN and CO share enhanced signaling pathways and transcriptional factors associated with immune response, and possesse similar T-cell-receptor features, such as high clonal expansion. Interestingly, rather than a significant increase in GZMK+ CD8 cells during aging, centenarians show accumulation of GZMB+ and CMC1+ CD8 T cells. Collectively, this study unveils an immune remodeling pattern reflected by both quantitative increase and functional reinforcement of cytotoxic T cells which are essential for healthy aging.

Three Major Gastrointestinal Cancers Could Be Distinguished through Subclass-Specific IgG Glycosylation.

Esophageal cancer (EC), gastric cancer (GC), and colorectal cancer (CRC) are three major digestive tract tumors with higher morbidity and mortality due to significant molecular heterogeneity. Altered IgG glycosylation has been observed in inflammatory activities and disease progression, and the IgG glycome profile could be used for disease stratification. However, IgG N-glycome profiles in these three cancers have not been systematically investigated. Herein, subclass-specific IgG glycosylation in CRC, GC, and EC was comprehensively characterized by liquid chromatography-tandem mass spectrometry. It was found that IgG1 sialylation was decreased in all three cancers, and the alterations in CRC and EC may be subclass-specific. IgG4 mono-galactosylation was increased in all three cancers, which was a subclass-specific change in all of them. Additionally, glycopeptides of IgG1-H5N5, IgG2-H4N3F1, and IgG4-H4N4F1 could distinguish all three cancer groups from controls with fair diagnostic performance. Furthermore, bioinformatics verified the differential expression of relevant glycosyltransferase genes in cancer progression. Significantly, those three gastrointestinal cancers could be distinguished from each other using subclass-specific IgG glycans. These findings demonstrated the spatial and temporal diversity of IgG N-glycome among digestive cancers, increasing our understanding of the molecular mechanisms of EC, GC, and CRC pathogenesis.

Multi-omics characterization of autophagy-related molecular features for therapeutic targeting of autophagy.

Autophagy is a major contributor to anti-cancer therapy resistance. Many efforts have been made to understand and overcome autophagy-mediated therapy resistance, but these efforts have been unsuccessful in clinical applications. In this study, we establish an autophagy signature to estimate tumor autophagy status. We then classify approximately 10,000 tumor samples across 33 cancer types from The Cancer Genome Atlas into autophagy score-high and autophagy score-low groups. We characterize the associations between multi-dimensional molecular features and tumor autophagy, and further analyse the effects of autophagy status on drug response. In contrast to the conventional view that the induction of autophagy serves as a key resistance mechanism during cancer therapy, our analysis reveals that autophagy induction may also sensitize cancer cells to anti-cancer drugs. We further experimentally validate this phenomenon for several anti-cancer drugs in vitro and in vivo, and reveal that autophagy inducers potentially sensitizes tumor cells to etoposide through downregulating the expression level of DDIT4. Our study provides a comprehensive landscape of molecular alterations associated with tumor autophagy and highlights an opportunity to leverage multi-omics analysis to utilize multiple drug sensitivity induced by autophagy.

ICBatlas: A Comprehensive Resource for Depicting Immune Checkpoint Blockade Therapy Characteristics from Transcriptome Profiles.

Immune checkpoint blockade (ICB) therapy provides remarkable clinical benefits for multiple cancer types. Much work is currently being conducted to investigate the mechanisms of ICB therapy at the transcriptional level. Integrating the data produced by these studies will help us give more insight into the transcriptomic features of ICB therapy. We collected the transcriptome and clinical data of ICB-treated patient samples from the Gene Expression Omnibus, ArrayExpress, The Cancer Genome Atlas, and dbGaP databases. On the basis of the clinical information, all samples are initially classified into response/nonresponse or pretreatment/on-treatment groups. Differential expression, pathway enrichment, and immune cell infiltration analyses are performed between the samples from different groups. We also introduce the Response Score (RS) calculated by integrating the variability degree and the frequency of the dysregulated genes in the responders to evaluate the impact of gene expression on the response. Finally, all the abovementioned contents are integrated into the ICBatlas database. ICBatlas provides the transcriptome features of ICB therapy through the analysis of 1,515 ICB-treated samples from 25 studies across nine cancer types. The data in ICBatlas include clinical outcomes, treatment-related genes, biological pathways, and immune cell infiltration. Users can investigate the abovementioned transcriptome features in the response (R vs. NR) or treatment (Pre vs. On) modules at the data set, cancer type, or immune checkpoint level and compare the degree of gene impact on the response in the RS module. ICBatlas is the first database to show the transcriptome features on ICB therapy in human cancers and freely available at http://bioinfo.life.hust.edu.cn/ICBatlas/.