Integrated Analysis of the ETS Family in Melanoma Reveals a Regulatory Role of ETV7 in the Immune Microenvironment

Single-Nuclei Transcriptome Profiling Reveals Intra-Tumoral Heterogeneity and Characterizes Tumor Microenvironment Architecture in a Murine Melanoma Model

Malignant melanoma is an aggressive cancer, with a high risk of metastasis and mortality rates, characterized by cancer cell heterogeneity and complex tumor microenvironment (TME). Single cell biology is an ideal and powerful tool to address these features at a molecular level. However, this approach requires enzymatic cell dissociation that can influence cellular coverage. By contrast, single nucleus RNA sequencing (snRNA-seq) has substantial advantages including compatibility with frozen samples and the elimination of a dissociation-induced, transcriptional stress response. To better profile and understand the functional diversity of different cellular components in melanoma progression, we performed snRNA-seq of 16,839 nuclei obtained from tumor samples along the growth of murine syngeneic melanoma model carrying a BRAFV600E mutation and collected 9 days or 23 days after subcutaneous cell injection. We defined 11 different subtypes of functional cell clusters among malignant cells and 5 different subsets of myeloid cells that display distinct global transcriptional program and different enrichment in early or advanced stage of tumor growth, confirming that this approach was useful to accurately identify intratumor heterogeneity and dynamics during tumor evolution. The current study offers a deep insight into the biology of melanoma highlighting TME reprogramming through tumor initiation and progression, underlying further discovery of new TME biomarkers which may be potentially druggable.

Functional interaction between receptor tyrosine kinase MET and ETS transcription factors promotes prostate cancer progression

Prostate cancer, the most common malignancy in men, has a relatively favourable prognosis. However, when it spreads to the bone, the survival rate drops dramatically. The development of bone metastases leaves patients with aggressive prostate cancer, the leading cause of death in men. Moreover, bone metastases are incurable and very painful. Hepatocyte growth factor receptor (MET) and fusion of genes encoding E26 transformation-specific (ETS) transcription factors are both involved in the progression of the disease. ETS gene fusions, in particular, have the ability to induce the migratory and invasive properties of prostate cancer cells, whereas MET receptor, through its signalling cascades, is able to activate transcription factor expression. MET signalling and ETS gene fusions are intimately linked to high-grade prostate cancer. However, the collaboration of these factors in prostate cancer progression has not yet been investigated. Here, we show, using cell models of advanced prostate cancer, that ETS translocation variant 1 (ETV1) and transcriptional regulator ERG (ERG) transcription factors (members of the ETS family) promote tumour properties, and that activation of MET signalling enhances these effects. By using a specific MET tyrosine kinase inhibitor in a humanised hepatocyte growth factor (HGF) mouse model, we also establish that MET activity is required for ETV1/ERG-mediated tumour growth. Finally, by performing a comparative transcriptomic analysis, we identify target genes that could play a relevant role in these cellular processes. Thus, our results demonstrate for the first time in prostate cancer models a functional interaction between ETS transcription factors (ETV1 and ERG) and MET signalling that confers more aggressive properties and highlight a molecular signature characteristic of this combined action.

Identifying cell type-specific transcription factor-mediated activity immune modules reveal implications for immunotherapy and molecular classification of pan-cancer

Systematic investigation of tumor-infiltrating immune (TII) cells is important to the development of immunotherapies, and the clinical response prediction in cancers. There exists complex transcriptional regulation within TII cells, and different immune cell types display specific regulation patterns. To dissect transcriptional regulation in TII cells, we first integrated the gene expression profiles from single-cell datasets, and proposed a computational pipeline to identify TII cell type-specific transcription factor (TF) mediated activity immune modules (TF-AIMs). Our analysis revealed key TFs, such as BACH2 and NFKB1 play important roles in B and NK cells, respectively. We also found some of these TF-AIMs may contribute to tumor pathogenesis. Based on TII cell type-specific TF-AIMs, we identified eight CD8+ T cell subtypes. In particular, we found the PD1 + CD8+ T cell subset and its specific TF-AIMs associated with immunotherapy response. Furthermore, the TII cell type-specific TF-AIMs displayed the potential to be used as predictive markers for immunotherapy response of cancer patients. At the pan-cancer level, we also identified and characterized six molecular subtypes across 9680 samples based on the activation status of TII cell type-specific TF-AIMs. Finally, we constructed a user-friendly web interface CellTF-AIMs (http://bio-bigdata.hrbmu.edu.cn/CellTF-AIMs/) for exploring transcriptional regulatory pattern in various TII cell types. Our study provides valuable implications and a rich resource for understanding the mechanisms involved in cancer microenvironment and immunotherapy.

ETV7 promotes colorectal cancer progression through upregulation of IFIT3

Members of the E26 transformation-specific (ETS) variant transcription factor family act as either tumor suppressors or oncogenic factors in numerous types of cancer. ETS variant transcription factor 7 (ETV7) participates in the development of malignant tumors, whereas its involvement in colorectal cancer (CRC) is less clear. In this study, The Cancer Genome Atlas (TCGA) and immunochemistry staining were applied to check the clinical relevance of ETV7 and interferon-induced protein with tetratricopeptide repeats 3 (IFIT3) in CRC patients. Overexpression and knockdown of ETV7 and IFIT3 were conducted by transfecting the cells with pCDNA3.1 plasmids and siRNAs, respectively. Western blotting was used to detect the protein expression of ETV7 in CRC cells. Cell Counting Kit-8, cell colony formation, and Transwell assays, as well as flow cytometry, were used to evaluate the proliferation, migration, cell cycle, and apoptosis of CRC cells. Furthermore, western blotting, RT-qPCR, and luciferase assay were used to explore the regulation of ETV7 on IFIT3. Rescue assay was used to investigate the significance of ETV7/IFIT3 axis on CRC progression. We found that ETV7 was upregulated in CRC tissues and cells. Overexpression of ETV7 stimulated the proliferation, migration, and cell cycle amplification, and reduced the apoptosis of CRC cells. Downregulation of ETV7 exerted the opposite effect on CRC cell progression. Moreover, we demonstrated that ETV7 stimulated the transcription activity, the mRNA and protein expression of IFIT3 in CRC cells. There was a positive correlation between ETV7 and IFIT3 in CRC patients. IFIT3 knockdown reversed the promotive effect exerted by overexpression of ETV7 on the amplification and migration of CRC cells. By contrast, overexpression of IFIT3 blocked the inhibitory effect of ETV7-targeting siRNA. In summary, ETV7 induces progression of CRC by activating the transcriptional expression of IFIT3. The EVT7/IFIT3 axis may be a novel target for CRC therapy.

Tannins in Phyllanthus emblica L. improves cisplatin efficacy in lung cancer cells by boosting endoplasmic reticulum stress to trigger immunogenic cell death

BACKGROUND

Lung cancer is one of the deadliest cancers world-wide and immunotherapy has been considered as a promising therapeutic strategy. Previously, our study found that tannins in Phyllanthus emblica L. (PTF) could inhibit the growth of tumor by activating the immune response in liver cancer, and also exhibited a cytotoxicity on human lung cancer cells A549, H460, H1703 in vitro.

OBJECTIVE

To explore whether PTF inhibited the growth of lung cancer through its immune-regulating function and to clarify underlying mechanisms.

METHODS

The induction of immunogenic cell death (ICD) were characterized by calreticulin exposure, extracellular ATP secretion, and High Mobility Group Box 1(HMGB1) release both in vivo using LLC-derived xenograft tumor model and in vitro using both mouse LLC and human A549 cancer cells.

RESULTS

PTF inhibited lung cancer cells growth and tumorigenesis in vivo/vitro and promoted anti-tumor immune responses. We further found that PTF could induce ICD, which then activated Type I interferon responses and CXCL9/10-mediated chemotaxis. Mechanistically, PTF induced the formation of intracellular protein aggregates and following activation of PERK/ATF4/CHOP-dependent endoplasmic reticulum stress-related ICD. Moreover, PTF improved the antitumor efficacy of cisplatin by inducing ICD both in vitro and in vivo. Finally, we screened out 5 components from PTF, including gallocatechin, gallic acid, methyl gallate, ethyl gallate and ellagic acid, which could induce ICD in vitro and might be considered as the potential antitumor pharmacodynamic substances.

CONCLUSION

In conclusion, PTF inhibits the growth of lung cancer by triggering ICD and remodeling the tumor microenvironment, suggesting that PTF may have promising prospects as an adjacent immunotherapy for cancers.

Functional analysis of recurrent CDC20 promoter variants in human melanoma

Small nucleotide variants in non-coding regions of the genome can alter transcriptional regulation, leading to changes in gene expression which can activate oncogenic gene regulatory networks. Melanoma is heavily burdened by non-coding variants, representing over 99% of total genetic variation, including the well-characterized TERT promoter mutation. However, the compendium of regulatory non-coding variants is likely still functionally under-characterized. We developed a pipeline to identify hotspots, i.e. recurrently mutated regions, in melanoma containing putatively functional non-coding somatic variants that are located within predicted melanoma-specific regulatory regions. We identified hundreds of statistically significant hotspots, including the hotspot containing the TERT promoter variants, and focused on a hotspot in the promoter of CDC20. We found that variants in the promoter of CDC20, which putatively disrupt an ETS motif, lead to lower transcriptional activity in reporter assays. Using CRISPR/Cas9, we generated an indel in the CDC20 promoter in human A375 melanoma cell lines and observed decreased expression of CDC20, changes in migration capabilities, increased growth of xenografts, and an altered transcriptional state previously associated with a more proliferative and less migratory state. Overall, our analysis prioritized several recurrent functional non-coding variants that, through downregulation of CDC20, led to perturbation of key melanoma phenotypes.

A novel CD8+ T cell-related gene signature for predicting the prognosis and immunotherapy efficacy in bladder cancer

OBJECTIVE

To identify CD8+ T cell-related molecular clusters and establish a novel gene signature for predicting the prognosis and efficacy of immunotherapy in bladder cancer (BCa).

METHODS

Transcriptome and clinical data of BCa samples were obtained from the Cancer Genome Atlas (TCGA) and GEO databases. The CD8+ T cell-related genes were screened through the CIBERSORT algorithm and correlation analysis. Consensus clustering analysis was utilized to identified CD8+ T cell-related molecular clusters. A novel CD8+ T cell-related prognostic model was developed using univariate Cox regression analysis and Lasso regression analysis. Internal and external validations were performed and the validity of the model was validated in a real-world cohort. Finally, preliminary experimental verifications were carried out to verify the biological functions of SH2D2A in bladder cancer.

RESULTS

A total of 52 CD8+ T cell-related prognostic genes were screened and two molecular clusters with notably diverse immune cell infiltration, prognosis and clinical features were developed. Then, a novel CD8+ T cell-related prognostic model was constructed. The patients with high-risk scores exhibited a significantly worse overall survival in training, test, whole TCGA and validating cohort. The AUC was 0.766, 0.725, 0.739 and 0.658 in the four cohorts sequentially. Subgroup analysis suggested that the novel prognostic model has a robust clinical application for selecting high-risk patients. Finally, we confirmed that patients in the low-risk group might benefit more from immunotherapy or chemotherapy, and validated the prognostic model in a real-world immunotherapy cohort. Preliminary experiment showed that SH2D2A was capable of attenuating proliferation, migration and invasion of BCa cells.

CONCLUSIONS

CD8+ T cell-related molecular clusters were successfully identified. Besides, a novel CD8+ T cell-related prognostic model with an excellent predictive performance in predicting survival rates and immunotherapy efficacy of BCa was developed.

RNAseq profiling of blood from patients with coronary artery disease: Signature of a T cell imbalance

Background

Cardiovascular disease had a global prevalence of 523 million cases and 18.6 million deaths in 2019. The current standard for diagnosing coronary artery disease (CAD) is coronary angiography either by invasive catheterization (ICA) or computed tomography (CTA). Prior studies employed single-molecule, amplification-independent RNA sequencing of whole blood to identify an RNA signature in patients with angiographically confirmed CAD. The present studies employed Illumina RNAseq and network co-expression analysis to identify systematic changes underlying CAD.

Methods

Whole blood RNA was depleted of ribosomal RNA (rRNA) and analyzed by Illumina total RNA sequencing (RNAseq) to identify transcripts associated with CAD in 177 patients presenting for elective invasive coronary catheterization. The resulting transcript counts were compared between groups to identify differentially expressed genes (DEGs) and to identify patterns of changes through whole genome co-expression network analysis (WGCNA).

Results

The correlation between Illumina amplified RNAseq and the prior SeqLL unamplified RNAseq was quite strong (r = 0.87), but there was only 9 % overlap in the DEGs identified. Consistent with the prior RNAseq, the majority (93 %) of DEGs were down-regulated ~1.7-fold in patients with moderate to severe CAD (>20 % stenosis). DEGs were predominantly related to T cells, consistent with known reductions in Tregs in CAD. Network analysis did not identify pre-existing modules with a strong association with CAD, but patterns of T cell dysregulation were evident. DEGs were enriched for transcripts associated with ciliary and synaptic transcripts, consistent with changes in the immune synapse of developing T cells.

Conclusions

These studies confirm and extend a novel mRNA signature of a Treg-like defect in CAD. The pattern of changes is consistent with stress-related changes in the maturation of T and Treg cells, possibly due to changes in the immune synapse.

ETV7 reduces inflammatory responses in breast cancer cells by repressing the TNFR1/NF-κB axis

The transcription factor ETV7 is an oncoprotein that is up-regulated in all breast cancer (BC) types. We have recently demonstrated that ETV7 promoted breast cancer progression by increasing cancer cell proliferation and stemness and was also involved in the development of chemo- and radio-resistance. However, the roles of ETV7 in breast cancer inflammation have yet to be studied. Gene ontology analysis previously performed on BC cells stably over-expressing ETV7 demonstrated that ETV7 was involved in the suppression of innate immune and inflammatory responses. To better decipher the involvement of ETV7 in these signaling pathways, in this study, we identified TNFRSF1A, encoding for the main receptor of TNF-α, TNFR1, as one of the genes down-regulated by ETV7. We demonstrated that ETV7 directly binds to the intron I of this gene, and we showed that the ETV7-mediated down-regulation of TNFRSF1A reduced the activation of NF-κB signaling. Furthermore, in this study, we unveiled a potential crosstalk between ETV7 and STAT3, another master regulator of inflammation. While it is known that STAT3 directly up-regulates the expression of TNFRSF1A, here we demonstrated that ETV7 reduces the ability of STAT3 to bind to the TNFRSF1A gene via a competitive mechanism, recruiting repressive chromatin remodelers, which results in the repression of its transcription. The inverse correlation between ETV7 and TNFRSF1A was confirmed also in different cohorts of BC patients. These results suggest that ETV7 can reduce the inflammatory responses in breast cancer through the down-regulation of TNFRSF1A.

ETS transcription factors: Multifaceted players from cancer progression to tumor immunity

The E26 transformation specific (ETS) family comprises 28 transcription factors, the majority of which are involved in tumor initiation and development. Serving as a group of functionally heterogeneous gene regulators, ETS factors possess a structurally conserved DNA-binding domain. As one of the most prominent families of transcription factors that control diverse cellular functions, ETS activation is modulated by multiple intracellular signaling pathways and post-translational modifications. Disturbances in ETS activity often lead to abnormal changes in oncogenicity, including cancer cell survival, growth, proliferation, metastasis, genetic instability, cell metabolism, and tumor immunity. This review systematically addresses the basics and advances in studying ETS factors, from their tumor relevance to clinical translational utility, with a particular focus on elucidating the role of ETS family in tumor immunity, aiming to decipher the vital role and clinical potential of regulation of ETS factors in the cancer field.

Exploration of immune response mechanisms in cadmium and copper co-exposed juvenile golden cuttlefish (Sepia esculenta) based on transcriptome profiling

Sepia esculenta is a popular economic cephalopod with high yield, delicious meat, and rich nutrition. With the rapid development of heavy industry and medical industry, a large amount of waste has been released into the ocean recklessly in recent years, inducing a significant increase in the content of heavy metals, especially cadmium (Cd) and copper (Cu), in the ocean. This phenomenon significantly affects the growth and development of S. esculenta, causing a serious blow to its artificial breeding. In this study, transcriptome analysis is used to initially explore immune response mechanisms of Cd and Cu co-exposed juvenile S. esculenta. The results show that 1,088 differentially expressed genes (DEGs) are identified. And DEGs functional enrichment analysis results suggests that co-exposure may promote inflammatory and innate immune responses in juvenile S. esculenta. Fifteen key genes that might regulate the immunity of S. esculenta are identified using protein-protein interaction (PPI) network and KEGG enrichment analyses, of which the three genes with the highest number of interactions or involve in more KEGG pathways are identified as hub genes that might significantly affect the immune response processes. Comprehensive analysis of PPI network and KEGG signaling pathway is used for the first time to explore co-exposed S. esculenta juvenile immune response processes. Our results preliminarily reveal immune response mechanisms of cephalopods exposed to heavy metals and provide a valuable resource for further understanding of mollusk immunity.

Core immune cell infiltration signatures identify molecular subtypes and promote precise checkpoint immunotherapy in cutaneous melanoma

Yutao Wang, China Medical University, ChinaThe tumor microenvironment (TME) has been shown to impact the prognosis of tumors in patients including cutaneous melanoma (CM); however, not all components of TME are important. Given the aforementioned situation, the functional immune cell contents correlated with CM patient prognosis are needed to optimize present predictive models and reflect the overall situation of TME. We developed a novel risk score named core tumor-infiltrating immune cell score (cTICscore), which showed certain advantages over existing biomarkers or TME-related signatures in predicting the prognosis of CM patients. Furthermore, we explored a new gene signature named cTILscore−related module gene score (cTMGs), based on four identified TME-associated genes (GCH1, GZMA, PSMB8, and PLAAT4) showing a close correlation with the cTICscore, which was generated by weighted gene co-expression network analysis and least absolute shrinkage and selection operator analysis to facilitate clinical application. Patients with low cTMGs had significantly better overall survival (OS, P = 0.002,< 0.001, = 0.002, and = 0.03, respectively) in the training and validating CM datasets. In addition, the area under the curve values used to predict the immune response in four CM cohorts were 0.723, 0.723, 0.754, and 0.792, respectively, and that in one gastric cohort was 0.764. Therefore, the four-gene signature, based on cTICscore, might improve prognostic information, serving as a predictive tool for CM patients receiving immunotherapy.cutaneous melanoma, tumor microenvironment, prognosis, immunotherapy, cTICscore

Bioinformatics-based analysis of SUMOylation-related genes in hepatocellular carcinoma reveals a role of upregulated SAE1 in promoting cell proliferation

The function of small ubiquitin-like modifier (SUMO)-related genes in hepatocellular carcinoma (HCC) remains unclear. This study aimed to analyze the expression profile and prognostic relevance of SUMO-related genes using publicly available data. A set of bioinformatics tools and experiments were integrated to explore the mechanism of the genes of interest. The least absolute shrinkage and selection operator Cox regression analysis was used to construct a prognostic model. SUMO-2 and SUMO-activating enzyme subunit 1 (SAE1) were upregulated in HCC. The enrichment analysis indicated that SUMO-2 and SAE1 might regulate the cell cycle. The downregulation of SAE1 inhibited the proliferation of HCC cells, whereas the upregulation of the gene promoted cell proliferation. IGF2BP3 contributed to the upregulation of SAE1 in an N6-methyladenosine (m6A)-dependent way. Eventually, an SAE1-related risk score (SRRS) was developed and validated in HCC. SRRS could serve as an independent prognostic factor and predict the efficiency of transarterial chemoembolization in patients with HCC.

Cuproptosis-Related Risk Score Predicts Prognosis and Characterizes the Tumor Microenvironment in Hepatocellular Carcinoma

Aims

Cuproptosis is a recently identified form of programmed cell death; however, its role in hepatocellular carcinoma (HCC) remains unclear.

Methods

A set of bioinformatic tools was integrated to analyze the expression and prognostic significance of ferredoxin 1 (FDX1), the key regulator of cuproptosis. A cuproptosis-related risk score (CRRS) was developed via correlation analyses, least absolute shrinkage and selection operator (LASSO) Cox regression, and multivariate Cox regression. The metabolic features, mutation signatures, and immune profile of CRRS-classified HCC patients were investigated, and the role of CRRS in therapy guidance was analyzed.

Results

FDX1 was significantly downregulated in HCC, and its high expression was associated with longer survival time. HCC patients in the high-CRRS group showed a significantly lower overall survival (OS) and enriched in cancer-related pathways. Mutation analyses revealed that the high-CRRS HCC patients had a high mutational frequency of some tumor suppressors such as tumor protein P53 (TP53) and Breast-cancer susceptibility gene 1 (BRCA1)-associated protein 1 (BAP1) and a low frequency of catenin beta 1 (CTNNB1). Besides, HCC patients with high CRRS showed an increase of protumor immune infiltrates and a high expression of immune checkpoints. Moreover, the area under the curve (AUC) values of CRRS in predicting the efficiency of sorafenib and the non-responsiveness to transcatheter arterial chemoembolization (TACE) in HCC patients reached 0.877 and 0.764, respectively.

Significance

The cuproptosis-related signature is helpful in prognostic prediction and in guiding treatment for HCC patients.

Comprehensive Analysis Identified ETV7 as a Potential Prognostic Biomarker in Bladder Cancer

Background

The tumor microenvironment (TME) plays a crucial role in the initiation and progression of cancer. Bladder cancer (BLCA) is a malignant tumor of the genitourinary system. Its heterogeneity results in significant differences in the prognosis of patients. To date, this is still a huge challenge for clinical treatment. In recent years, more and more evidence showed that dysregulation of transcription factors (TFs) plays an important role in tumor progression, invasion, and metastasis. Unfortunately, the role of TFs on the tumor microenvironment in bladder cancer is unclear.

Methods

The original data of BLCA and corresponding adjacent tissues were obtained from The Cancer Genome Atlas (TCGA) database. TFs were downloaded from the Animal Transcription Factor DataBase (Animal TFDB). Intersection analysis was used to obtain TFs that were differentially expressed between tumor and adjacent tissues. Gene Set Cancer Analysis (GSCALite) and CIBERSORT software were used to reveal the key differentially expressed TFs (DE-TFs). Subsequently, UALCAN and Human Protein Atlas (HPA) databases were used to disclose the expression of key DE-TFs in BLCA. The K-M curve divulged the relationship between the key DE-TFs and the patient's overall survival (OS), and the univariate and multivariate Cox regression analyses were conducted to explore independent prognostic factors. The cluster profiler package and Gene Set Enrichment Analysis (GSEA) were used for functional enrichment of genes related to the key DE-TFs. Finally, CIBERSORT software analyzed the immune landscape of BLCA.

Results

We obtained a total of 117 BLCA-related DE-TFs. Among them, ETV7 was identified as the key DE-TFs due to its association with the autophagy activation pathway and various immune cells in cancer. Online databases of UALCAN and HPA indicated that ETV7 was overexpressed in tumors and negatively correlated with tumor severity. The K-M curve showed that the OS of patients with high expression of ETV7 was poor, which indicated that it was an independent prognostic factor. Functional enrichment of 87 DEGs between ETV7-high and -low expression groups indicated that it was closely related to the immune response and the functions of a variety of immune cells. Finally, CIBERSORT results proved that the high and low expression of ETV7 also caused significant differences in the tumor immune microenvironment of patients.

Conclusion

Overall, we proved that the transcription factor ETV7 was a novel prognostic factor, which may improve the individualized outcome prediction in BLCA by regulating the tumor immune microenvironment.

Tumor purity as a prognosis and immunotherapy relevant feature in cervical cancer

Background: Tumor purity plays a vital role in the biological process of solid tumors, but its function in gynecologic cancers remains unclear. This study explored the correlation between tumor purity and immune function of gynecological cancers and its reliability as a prognostic indicator of immunotherapy.

Methods: Gynecological cancer-related datasets were downloaded from The Cancer Genome Atlas (TCGA). Tumor purity was calculated by the ESTIMATE algorithm. A LASSO Cox regression analysis was performed to construct the risk score model. A Kaplan–Meier Plotter was used to explore the relationships between tumor purity and cancer prognosis. We performed the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Set Enrichment Analysis (GSEA) to explore the pathways in the subgroups. A nomogram was used to quantitatively assess the cancer prognosis.

Results: Tumor purity was negatively correlated with B cell infiltration in cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC). Approximately 420 genes were positively associated with B cell infiltration and CESC prognosis and were enriched in immune-related signaling pathways. There were 11 key genes used to construct a risk score model. The low-risk group had a higher immune score and better prognosis than the high-risk group. A nomogram based on risk score, T stage, and clinical-stage had good predictive value in quantitatively evaluating CESC prognosis.

Conclusions: This study is the first to reveal the correlation between tumor purity and immunity in CESC and suggests that low-risk patients may be more sensitive to immunotherapy. This provides a theoretical basis for the clinical treatment of CESC.

ImmReg: the regulon atlas of immune-related pathways across cancer types

Immune system gene regulation perturbation has been found to be a major cause of the development of various types of cancer. Numbers of mechanisms contribute to gene expression regulation, thus, systematically identification of potential regulons of immune-related pathways is critical to cancer immunotherapy. Here, we comprehensively chart the landscape of transcription factors, microRNAs, RNA binding proteins and long noncoding RNAs regulation in 17 immune-related pathways across 33 cancers. The potential immunology regulons are likely to exhibit higher expressions in immune cells, show expression perturbations in cancer, and are significantly correlated with immune cell infiltrations. We also identify a panel of clinically relevant immunology regulons across cancers. Moreover, the regulon atlas of immune-related pathways helps prioritizing cancer-related genes (i.e. ETV7, miR-146a-5p, ZFP36 and HCP5). We further identified two molecular subtypes of glioma (cold and hot tumour phenotypes), which were characterized by differences in immune cell infiltrations, expression of checkpoints, and prognosis. Finally, we developed a user-friendly resource, ImmReg (http://bio-bigdata.hrbmu.edu.cn/ImmReg/), with multiple modules to visualize, browse, and download immunology regulation. Our study provides a comprehensive landscape of immunology regulons, which will shed light on future development of RNA-based cancer immunotherapies.

Multi-Omics Analysis of Anlotinib in Pancreatic Cancer and Development of an Anlotinib-Related Prognostic Signature

Aberrant regulation of angiogenesis involves in the growth and metastasis of tumors, but angiogenesis inhibitors fail to improve overall survival of pancreatic cancer patients in previous phase III clinical trials. A comprehensive knowledge of the mechanism of angiogenesis inhibitors against pancreatic cancer is helpful for clinical purpose and for the selection of patients who might benefit from the inhibitors. In this work, multi-omics analyses (transcriptomics, proteomics, and phosphoproteomics profiling) were carried out to delineate the mechanism of anlotinib, a novel angiogenesis inhibitor, against pancreatic cancer cells. The results showed that anlotinib exerted noteworthy cytotoxicity on pancreatic cancer cells. Multi-omics analyses revealed that anlotinib had a profound inhibitory effect on ribosome, and regulated cell cycle, RNA metabolism and lysosome. Based on the multi-omics results and available data deposited in public databases, an anlotinib-related gene signature was further constructed to identify a subgroup of pancreatic cancer patients who had a dismal prognosis and might be responsive to anlotinib.