Visualizing and interpreting cancer genomics data via the Xena platform

Quantitative acetylome analysis reveals histone modifications that may predict prognosis in hepatitis B-related hepatocellular carcinoma

Lysine acetylation (Kac) as an important posttranslational modification of histones is essential for the regulation of gene expression in hepatocellular carcinoma (HCC). However, the atlas of whole acetylated proteins in HCC tissues and the difference in protein acetylation between normal human tissues and HCC tissues are unknown. In this report, we characterized the proteome and acetyl proteome (acetylome) profile of normal, paracancerous, and HCC liver tissues in human clinical samples by quantitative proteomics techniques. We identified 6781 acetylation sites of 2582 proteins and quantified 2492 acetylation sites of 1190 proteins in normal, paracancerous, and HCC liver tissues. Among them, 15 proteins were multiacetylated with more than 10 lysine residues. The histone acetyltransferases p300 and CBP were found to be hyperacetylated in hepatitis B virus pathway. Moreover, we found that 250 Kac sites of 214 proteins were upregulated and 662 Kac sites of 451 proteins were downregulated in HCC compared with normal liver tissues. Additionally, the acetylation levels of lysine 120 in histone H2B (H2BK120ac), lysine 18 in histone H3.3 (H3.3K18ac), and lysine 77 in histone H4 (H4K77ac) were increased in HCC. Interestingly, the higher levels of H2BK120ac, H3.3K18ac, and H4K77ac were significantly associated with worse prognosis, such as poorer survival and higher recurrence in an independent clinical cohort of HCC patients. Overall, this study lays a foundation for understanding the functions of acetylation in HCC and provides potential prognostic factors for the diagnosis and therapy of HCC.

A Systems Biology Approach to Investigating the Interaction between Serotonin Synthesis by Tryptophan Hydroxylase and the Metabolic Homeostasis

Obesity has become a global public health and economic problem. Obesity is a major risk factor for a number of complications, such as type 2 diabetes, cardiovascular disease, fatty liver disease, and cancer. Serotonin (5-hydroxytryptamine [5-HT]) is a biogenic monoamine that plays various roles in metabolic homeostasis. It is well known that central 5-HT regulates appetite and mood. Several 5-HT receptor agonists and selective serotonin receptor uptake inhibitors (SSRIs) have shown beneficial effects on appetite and mood control in clinics. Although several genetic polymorphisms related to 5-HT synthesis and its receptors are strongly associated with obesity, there is little evidence of the role of peripheral 5-HT in human metabolism. In this study, we performed a systemic analysis of transcriptome data from the Genotype-Tissue Expression (GTEX) database. We investigated the expression of 5-HT and tryptophan hydroxylase (TPH), the rate-limiting enzyme of 5-HT biosynthesis, in the human brain and peripheral tissues. We also performed differential gene expression analysis and predicted changes in metabolites by comparing gene expressions of tissues with high TPH expression to the gene expressions of tissues with low TPH expression. Our analyses provide strong evidence that serotonin plays an important role in the regulation of metabolic homeostasis in humans.

Elevation in viral entry genes and innate immunity compromise underlying increased infectivity and severity of COVID-19 in cancer patients

Multiple studies have reported a doubling in risk of Coronavirus Disease-2019 (COVID-19) among cancer patients. Here, we examine the potential biological rationale behind this recurrent epidemiological observation. By leveraging large-scale genome-wide transcriptional data of normal and malignant tissues from adults and children, we found evidence of increased expression of SARS-CoV-2 viral entry genes in the cancer state, particularly in respiratory, gastrointestinal, and genitourinary tract tissues, with decreased expression in pediatric vs. adult samples. Additionally, by interrogating the temporal effects of radiotherapy on human peripheral blood mononuclear and mucosal cells, we observed important treatment-related alterations in host innate immunity, specifically type I interferon responses. Overall, cancers enhance expression of critical viral entry genes, and innate viral defenses can be dysregulated transiently during radiation treatments. These factors may contribute to the observed increased susceptibility to SARS-CoV-2 entry and severity of COVID-19 in cancer patients.

Upregulation of hsa-miR-196a-5p is associated with MIR196A2 methylation and affects the malignant biological behaviors of glioma

Hsa-miR-196a-5p is involved in tumorigenesis and progression. However, the driving factors for hsa-miR-196a-5p overexpression and its correlation with the clinicopathological features and prognosis of patients remain unclear in glioma. Thus, this study aimed to investigate the prognostic value of hsa-miR-196a-5p and its correlation with MIR196A2 methylation in glioma. We observed that hsa-miR-196a-5p expression was upregulated in glioma. Next, 112 patients were divided into high (n = 56) and low (n = 56) hsa-miR-196a-5p expression groups. The chi-square test showed that hsa-miR-196a-5p expression was significantly related to age, WHO grade, histopathology, IDH mutation status, and 1p/19q codeletion. Univariate and multivariate Cox regression analyses showed that hsa-miR-196a-5p expression was an independent prognostic factor. GO and KEGG enrichment analyses showed that hsa-miR-196a-5p may be involved in the MAPK signaling, focal adhesion and cancer-related pathways. Compared with the normal astrocyte cell line, glioma cell lines had an unregulated MIR196A2 methylation level, which was confirmed by TCGA data. The hypermethylated CpG sites of MIR196A2 were mainly concentrated in the gene body region, which was significantly associated with hsa-miR-196a-5p overexpression. Kaplan-Meier curves revealed that MIR196A2 hypermethylation was a poor prognostic factor. These findings suggest that hsa-miR-196a-5p overexpression may be involved in malignant biological behaviors, and MIR196A2 hypermethylation of the gene body was significantly associated with hsa-miR-196a-5p overexpression, which was a poor prognostic factor of glioma. Therefore, MIR196A2 hypermethylation may act as an early marker of prognosis of patients with glioma.

Novel Small Molecule Hsp90/Cdc37 Interface Inhibitors Indirectly Target K-Ras-Signaling

The ATP-competitive inhibitors of Hsp90 have been tested predominantly in kinase addicted cancers; however, they have had limited success. A mechanistic connection between Hsp90 and oncogenic K-Ras is not known. Here, we show that K-Ras selectivity is enabled by the loss of the K-Ras membrane nanocluster modulator galectin-3 downstream of the Hsp90 client HIF-1α. This mechanism suggests a higher drug sensitivity in the context of KRAS mutant, HIF-1α-high and/or Gal3-high cancer cells, such as those found, in particular, in pancreatic adenocarcinoma. The low toxicity of conglobatin further indicates a beneficial on-target toxicity profile for Hsp90/Cdc37 interface inhibitors. We therefore computationally screened >7 M compounds, and identified four novel small molecules with activities of 4 μM-44 μM in vitro. All of the compounds were K-Ras selective, and potently decreased the Hsp90 client protein levels without inducing the heat shock response. Moreover, they all inhibited the 2D proliferation of breast, pancreatic, and lung cancer cell lines. The most active compounds from each scaffold, furthermore, significantly blocked 3D spheroids and the growth of K-Ras-dependent microtumors. We foresee new opportunities for improved Hsp90/Cdc37 interface inhibitors in cancer and other aging-associated diseases.

Molecular Features of Cancer-associated Fibroblast Subtypes and their Implication on Cancer Pathogenesis, Prognosis, and Immunotherapy Resistance

PURPOSE

Cancer-associated fibroblasts (CAFs) are an important component of the tumor microenvironment, but a systematic investigation of their molecular characteristics and clinical relevance are lacking. Here, we sought to compare CAFs across multiple cancer types to identify critical molecular pathways activated in CAF subtypes, which may contribute to clinical outcome, disease progression, and immunotherapy resistance.

EXPERIMENTAL DESIGN

We performed integrated analysis of CAFs from melanoma, head and neck squamous cell carcinoma, and lung cancer, and identified the molecular characteristics that are distinctly active in each CAF subtype. Gene signatures for individual CAF subtypes were identified and used to study the association of subtype abundance with clinical outcome and immunotherapy resistance.

RESULTS

We identified six CAF subtypes (pan-CAF) shared across cancer types and uncovered the molecular characteristics and genetic pathways distinguishing them. Interestingly, these CAF subtypes express distinct immunosuppressive factors, such as CXCL12 and CXLC14, and stem cell-promoting factor IL6. In addition, we identified novel transcriptional drivers (MEF2C, TWIST1, NR1H3, RELB, and FOXM1) key to CAF heterogeneity. Furthermore, we showed that CAF subtypes were associated with different clinical outcomes and uncovered key molecular pathways that could activate or suppress cancer progression or were involved in resistance to anti-PD1 or anti-PD-L1 immunotherapy.

CONCLUSIONS

Our study identifies the molecular characteristics of CAF subtypes shared across several cancer types, implicates cancer types that may benefit from CAF subtype targeted therapies, and identifies specific CAF subtypes associated with immunotherapy resistance.

Functional relationship between CFTR and RAC3 expression for maintaining cancer cell stemness in human colorectal cancer

PURPOSE

CFTR mutations not only cause cystic fibrosis, but also increase the risk of colorectal cancer. A putative role of CFTR in colorectal cancer patients without cystic fibrosis has so far, however, not been investigated. RAC3 is a nuclear receptor coactivator that has been found to be overexpressed in several human tumors, and to be required for maintaining cancer stemness. Here, we investigated the functional relationship between CFTR and RAC3 for maintaining cancer stemness in human colorectal cancer.

METHODS

Cancer stemness was investigated by analysing the expression of stem cell markers, clonogenic growth and selective retention of fluorochrome, using stable transfection of shCFTR or shRAC3 in HCT116 colorectal cancer cells. In addition, we performed pathway enrichment and network analyses in both primary human colorectal cancer samples (TCGA, Xena platform) and Caco-2 colorectal cancer cells including (1) CD133+ or CD133- side populations and (2) CFTRwt or CFTRmut cells (ConsensusPathDB, STRING, Cytoscape, GeneMANIA).

RESULTS

We found that the CD133+ side population expresses higher levels of RAC3 and CFTR than the CD133- side population. RAC3 overexpression increased CFTR expression, whereas CFTR downregulation inhibited the cancer stem phenotype. CFTR mRNA levels were found to be increased in colorectal cancer samples from patients without cystic fibrosis compared to those with CFTR mutations, and this correlated with an increased expression of RAC3. The expression pattern of a gene set involved in inflammatory response and nuclear receptor modulation in CD133+ Caco-2 cells was found to be shared with that in CFTRwt Caco-2 cells. These genes may contribute to colorectal cancer development.

CONCLUSIONS

CFTR may play a non-tumor suppressor role in colorectal cancer development and maintenance involving enhancement of the expression of a set of genes related to cancer stemness and development in patients without CFTR mutations.

Integrated Analysis Revealed the MicroRNA-Based Prognostic Predicting Signature for Papillary Renal Cell Carcinoma

Renal cell carcinoma (RCC) is one of the most frequently occurring tumors worldwide. Herein, we established a microRNA (miRNA) predicting signature to assess the prognosis of papillary-type RCC (PRCC) patients. miR-1293, miR-34a, miR-551b, miR-937, miR-299, and miR-3199-2 were used in building the overall survival (OS)-related signature, whereas miR-7156, miR-211, and miR-301b were used to construct the formula of recurrence-free survival (RFS) with the help of LASSO Cox regression analysis. The Kaplan-Meier and receiver operating characteristic curves indicated good discrimination and efficiency of the two signatures. Functional annotation for the downstream genes of the OS/RFS-related miRNAs exposed the potential mechanisms of PRCC. Notably, the multivariate analyses suggested that the two signatures were independent risk factors for PRCC patients and had better prognostic capacity than any other classifier. In addition, the nomogram indicated synthesis effects and showed better predictive performance than clinicopathologic features and our signatures. We validated the OS and RFS prediction formulas in clinical samples and met our expectations. Finally, we established two novel miRNA-based OS and RFS predicting signatures for PRCC, which are reliable tools for assessing the prognosis of PRCC patients.

Immune classification of clear cell renal cell carcinoma

Since the outcome of treatments, particularly immunotherapeutic interventions, depends on the tumor immune micro-environment (TIM), several experimental and computational tools such as flow cytometry, immunohistochemistry, and digital cytometry have been developed and utilized to classify TIM variations. In this project, we identify immune pattern of clear cell renal cell carcinomas (ccRCC) by estimating the percentage of each immune cell type in 526 renal tumors using the new powerful technique of digital cytometry. The results, which are in agreement with the results of a large-scale mass cytometry analysis, show that the most frequent immune cell types in ccRCC tumors are CD8+ T-cells, macrophages, and CD4+ T-cells. Saliently, unsupervised clustering of ccRCC primary tumors based on their relative number of immune cells indicates the existence of four distinct groups of ccRCC tumors. Tumors in the first group consist of approximately the same numbers of macrophages and CD8+ T-cells and and a slightly smaller number of CD4+ T cells than CD8+ T cells, while tumors in the second group have a significantly high number of macrophages compared to any other immune cell type (P-value \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$<0.01$$\end{document}<0.01). The third group of ccRCC tumors have a significantly higher number of CD8+ T-cells than any other immune cell type (P-value \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$<0.01$$\end{document}<0.01), while tumors in the group 4 have approximately the same numbers of macrophages and CD4+ T-cells and a significantly smaller number of CD8+ T-cells than CD4+ T-cells (P-value \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$<0.01$$\end{document}<0.01). Moreover, there is a high positive correlation between the expression levels of IFNG and PDCD1 and the percentage of CD8+ T-cells, and higher stage and grade of tumors have a substantially higher percentage of CD8+ T-cells. Furthermore, the primary tumors of patients, who are tumor free at the last time of follow up, have a significantly higher percentage of mast cells (P-value \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$<0.01$$\end{document}<0.01) compared to the patients with tumors for all groups of tumors except group 3.

Immune classification of clear cell renal cell carcinoma

Since the outcome of treatments, particularly immunotherapeutic interventions, depends on the tumor immune micro-environment (TIM), several experimental and computational tools such as flow cytometry, immunohistochemistry, and digital cytometry have been developed and utilized to classify TIM variations. In this project, we identify immune pattern of clear cell renal cell carcinomas (ccRCC) by estimating the percentage of each immune cell type in 526 renal tumors using the new powerful technique of digital cytometry. The results, which are in agreement with the results of a large-scale mass cytometry analysis, show that the most frequent immune cell types in ccRCC tumors are CD8+ T-cells, macrophages, and CD4+ T-cells. Saliently, unsupervised clustering of ccRCC primary tumors based on their relative number of immune cells indicates the existence of four distinct groups of ccRCC tumors. Tumors in the first group consist of approximately the same numbers of macrophages and CD8+ T-cells and and a slightly smaller number of CD4+ T cells than CD8+ T cells, while tumors in the second group have a significantly high number of macrophages compared to any other immune cell type (P-value [Formula: see text]). The third group of ccRCC tumors have a significantly higher number of CD8+ T-cells than any other immune cell type (P-value [Formula: see text]), while tumors in the group 4 have approximately the same numbers of macrophages and CD4+ T-cells and a significantly smaller number of CD8+ T-cells than CD4+ T-cells (P-value [Formula: see text]). Moreover, there is a high positive correlation between the expression levels of IFNG and PDCD1 and the percentage of CD8+ T-cells, and higher stage and grade of tumors have a substantially higher percentage of CD8+ T-cells. Furthermore, the primary tumors of patients, who are tumor free at the last time of follow up, have a significantly higher percentage of mast cells (P-value [Formula: see text]) compared to the patients with tumors for all groups of tumors except group 3.

Gastric Carcinomas with Stromal B7-H3 Expression Have Lower Intratumoural CD8+ T Cell Density

CD8+ T cells are the main effector cells of anti-cancer immune response that can be regulated by various costimulatory and coinhibitory molecules, including members of the B7 family. B7 homolog 3 (B7-H3) appears as a promising marker for immunotherapy; however, its significance in gastric cancer (GC) is unclear yet. We evaluated the spatial distribution of CD8+ T cells in relation to the expression of B7-H3 by double immunohistochemical staining. The level of B7-H3 intensity was scored manually (0–3) and dichotomized into B7-H3-low and B7-H3-high groups. The distribution and density of CD8+ T cells was analysed using whole slide digital imaging. B7-H3 was expressed mainly in the stromal compartment of GC (n = 73, 76% of all cases). Tumours with high expression of B7-H3 showed larger spatial differences of CD8+ T cells (86.4/mm² in tumour centre vs. 414.9/mm² in invasive front) when compared to B7-H3-low group (157.7/mm² vs. 218.7/mm², respectively) (p < 0.001). This study provides insight into the expression pattern of B7-H3 in GC of Western origin. In GCs with higher level of B7-H3 expression, CD8+ T cells were spatially suppressed in the tumour centre suggesting that B7-H3 might be involved in tumour escape mechanisms from the immune response.

Elevated LINC01550 induces the apoptosis and cell cycle arrest of melanoma

Melanoma is a high-grade malignant subtype of human skin cancer with the highest mortality rate. Here we perform a bioinformatics analysis concerning human melanoma tissues by the Gene Expression Omnibus (GEO) database and Gene Expression Profiling Interactive Analysis (GEPIA) platform. We found that lncRNA LINC01550 was significantly down-regulated in the melanoma tissues as compared to the normal tissues. The low expression of LINC01550 was tightly associated with shorter overall survival and disease-free survival of patients with melanoma. LINC01550 expression is negatively associated with tumor cell proliferation and invasion abilities in melanoma as evidenced by the single-cell RNA sequencing (scRNA-seq) databases. LINC01550-overexpressing vectors were transferred into melanoma cells (WM35 and WM451). Up-regulation of LINC01550 significantly inhibited proliferation and invasion abilities, as well as induced cell apoptosis and G1 and S phase arrest of the melanoma cells. In conclusion, overexpression of LINC01550 may serve as a potential therapeutic target for melanoma.

CD96 Correlates With Immune Infiltration and Impacts Patient Prognosis: A Pan-Cancer Analysis

Background

Immunotherapy has significantly improved patient outcomes, but encountered obstacles recently. CD96, a novel immune checkpoint expressed on T cells and natural killer (NK) cells, is essential for regulating immune functions. However, how CD96 correlating with immune infiltration and patient prognosis in pan-cancer remains unclear.

Methods

HPA, TCGA, GEO, GTEx, Oncomine, TIMER2.0, PrognoScan, Linkedomics, Metascape, and GEPIA2 databases were used to analyze CD96 in cancers. Visualization of data was mostly achieved by R language, version 4.0.2.

Results

In general, CD96 was differentially expressed between most cancer and adjacent normal tissues. CD96 significantly impacted the prognosis of diverse cancers. Especially, high CD96 expression was associated with poorer overall survival (OS) and disease-specific survival (DSS) in the TCGA lower grade glioma (LGG) cohort (OS, HR = 2.18, 95% CI = 1.79–2.66, P < 0.001). The opposite association was significantly observed in skin cutaneous melanoma (SKCM) cohort (OS, HR = 0.96, 95% CI = 0.94–0.98, P < 0.001). Notably, SKCM samples demonstrated the highest CD96 mutation frequency among all cancer types. Furthermore, in most cancers, CD96 expression level was significantly correlated with expression levels of recognized immune checkpoints and abundance of multiple immune infiltrates including CD8+ T cells, dendric cells (DCs), macrophages, monocytes, NK cells, neutrophils, regulatory T cells (Tregs), and follicular helper T cells (Tfh). CD96 was identified as a risk factor, protective factor, and irrelevant variable in LGG, SKCM and adrenocortical carcinoma (ACC), respectively. CD96 related genes were involved in negative regulation of leukocyte in LGG, however, involved in multiple positive immune processes in SKCM. Furthermore, CD96 was significantly associated with particular immune marker subsets. Importantly, it strongly correlated with markers of type 1 helper T cell (Th1) in SKCM, but not in LGG or ACC either.

Conclusions

CD96 participates in diverse immune responses, governs immune cell infiltration, and impacts malignant properties of various cancer types, thus standing as a potential biomarker for determining patient prognosis and immune infiltration in multiple cancers, especially in glioma and melanoma.

Colony-stimulating factor 3 signaling in colon and rectal cancers: Immune response and CMS classification in TCGA data

Colorectal cancer is the 2nd leading cause of cancer-related deaths in the world. The mechanisms underlying CRC development, progression, and resistance to treatment are complex and not fully understood. The immune response in the tumor microenvironment has been shown to play a significant role in many cancers, including colorectal cancer. Colony-stimulating factor 3 (CSF3) has been associated with changes to the immune environment in colorectal cancer animal models. We hypothesized that CSF3 signaling would correlate with pro-tumor tumor microenvironment changes associated with immune infiltrate and response. We utilized publicly available datasets to guide future mechanistic studies of the role CSF3 and its receptor (CSF3R) play in colorectal cancer development and progression. Here, we use bioinformatics data and mRNA from patients with colon (n = 242) or rectal (n = 92) cancers, obtained from The Cancer Genome Atlas Firehose Legacy dataset. We examined correlations of CSF3 and CSF3R expression with patient demographics, tumor stage and consensus molecular subtype classification. Gene expression correlations, cell type enrichment, Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data scores and Gene Ontology were used to analyze expression of receptor and ligand, tumor microenvironment infiltration of immune cells, and alterations in biological pathways. We found that CSF3 and CSF3R expression is highest in consensus molecular subtype 1 and consensus molecular subtype 4. Ligand and receptor expression are also correlated with changes in T cell and macrophage signatures. CSF3R significantly correlates with a large number of genes that are associated with poor colorectal cancer prognosis.

BUB1B Overexpression Is an Independent Prognostic Marker and Associated with CD44, p53, and PD-L1 in Renal Cell Carcinoma

INTRODUCTION

BUB1 mitotic checkpoint serine/threonine kinase B encoded by BUB1B gene is a member of the spindle assembly checkpoint family. Several reports have demonstrated that overexpression of BUB1B is associated with cancer progression and prognosis.

OBJECTIVE

This study aims to clarify the expression and function of BUB1B in renal cell carcinoma (RCC).

METHODS

The expression of BUB1B was determined using immunohistochemistry and bioinformatics analysis in RCC. The effects of BUB1B knockdown on cell growth and invasion were evaluated. We analyzed the interaction between BUB1B, cancer stem cell markers, p53, and PD-L1 in RCC.

RESULTS

In 121 cases of RCC, immunohistochemistry showed that 30 (25%) of the RCC cases were positive for BUB1B. High BUB1B expression was significantly correlated with high nuclear grade, T stage, and M stage. A Kaplan-Meier analysis showed that the high expression of BUB1B was associated with poor overall survival after nephrectomy. High BUB1B expression was associated with CD44, p53, and PD-L1 in RCC. Knockdown of BUB1B suppressed cell growth and invasion in RCC cell lines. Knockdown of BUB1B also suppressed the expression of CD44 and increased the expression of phospho-p53 (Ser15). In silico analysis showed that BUB1B was associated with inflamed CD8+, exhausted T-cell signature, IFN-γ signature, and the response to nivolumab.

CONCLUSION

These results suggest that BUB1B plays an oncogenic role and may be a promising predictive biomarker for survival in RCC.

Selective glutamine metabolism inhibition in tumor cells improves antitumor T lymphocyte activity in triple-negative breast cancer

Rapidly proliferating tumor and immune cells need metabolic programs that support energy and biomass production. The amino acid glutamine is consumed by effector T cells and glutamine-addicted triple-negative breast cancer (TNBC) cells, suggesting that a metabolic competition for glutamine may exist within the tumor microenvironment, potentially serving as a therapeutic intervention strategy. Here, we report that there is an inverse correlation between glutamine metabolic genes and markers of T cell-mediated cytotoxicity in human basal-like breast cancer (BLBC) patient data sets, with increased glutamine metabolism and decreased T cell cytotoxicity associated with poor survival. We found that tumor cell-specific loss of glutaminase (GLS), a key enzyme for glutamine metabolism, improved antitumor T cell activation in both a spontaneous mouse TNBC model and orthotopic grafts. The glutamine transporter inhibitor V-9302 selectively blocked glutamine uptake by TNBC cells but not CD8+ T cells, driving synthesis of glutathione, a major cellular antioxidant, to improve CD8+ T cell effector function. We propose a "glutamine steal" scenario, in which cancer cells deprive tumor-infiltrating lymphocytes of needed glutamine, thus impairing antitumor immune responses. Therefore, tumor-selective targeting of glutamine metabolism may be a promising therapeutic strategy in TNBC.

Molecular profiling of nucleocytoplasmic transport factor genes in breast cancer

Transport of functional molecules across the nuclear membrane of a eukaryotic cell is regulated by a dedicated set of transporter proteins that carry molecules into the nucleus or out of the nucleus to the cytoplasm for homeostasis of the cell. One of the categories of cargo molecules these transporters carry are the molecules for cell cycle regulation. Therefore, their role is critical in terms of cancer development. Any misregulation of the transport factors would means aberrant abundance of cell cycle regulators and might have consequences in cell cycle progression. While earlier studies have focussed on individual transport related molecules, a collective overview of how these molecules may be dysregulated in breast cancer is lacking.

Using genomic and transcriptomic datasets from TCGA (The Cancer Genome Atlas) and microarray platforms, we carried out bioinformatic analysis and provide a genetic and molecular profile of all the molecules directly related to nucleocytoplasmic shuttling of proteins and RNAs. Interestingly, we identified that many of these molecules are either mutated or have dysregulated expression in breast cancer. Strikingly, some of the molecules, namely, KPNA2, KPNA3, KPNA5, IPO8, TNPO1, XPOT, XPO7 and CSE1L were correlated with poor patient survival. This study provides a comprehensive genetic and molecular landscape of nucleocytoplasmic factors in breast cancer and points to the important roles of various nucleocytoplasmic factors in cancer progression. This data might have implications in prognosis and therapeutic targeting in breast cancer.

Low Expression of the Polycomb Protein RING1 Predicts Poor Prognosis in Human Breast Cancer

Background

To date, breast cancer remains the most common malignant tumor in women. In recent years, a growing number of studies on polycomb proteins have been conducted. The Ring finger protein1 (RING1), an essential component of the polycomb family of proteins, plays vital roles in the tumorigenesis of various cancer types. However, further research is required in determining RING1 expression and prognostic value in breast cancer.

Method

RING1 expression level in multiple cancer types was evaluated using the XENA and UALCAN databases. Real-time quantitative PCR (real-time qPCR) and immunohistochemistry (IHC) were used to confirm this expression. The prognostic value was analyzed using our follow-up data and the Kaplan-Meier plotter website. RING1 co-expressed genes and its promoter methylation level were calculated using the cBioPortal and UALCAN online tools. The gene ontology (GO) and the Kyoto encyclopedia of Genes and Genomes (KEGG) pathway enrichment were analyzed using the DAVID online analysis tool.

Result

RING1 expression was upregulated in CHOL (Bile Duct Cancer), ESCA (Esophageal Cancer), LIHC (Liver Cancer), and PCPG (Pheochromocytoma & Paraganglioma). However, its expression level was decreased in COAD (Colon Cancer), KICH (Kidney Chromophobe), KIRP (kidney papillary cell carcinoma), THCA (Thyroid Cancer), and BRCA (Breast carcinoma). RING1 low expression is an unfavorable prognostic factor in many cancer patients, especially in breast cancer patients. For breast cancer, the IHC result showed that RING1 protein expression significantly and negatively correlates with tumor size (P = 0.029), LNM (P = 0.017), TNM stage (P = 0.016), ER (P = 0.005), Ki67 (P = 0.015), and p53 status (P = 0.034). Moreover, the multivariate Cox regression model indicated that RING1 (P = 0.038) and ER (P = 0.029) expressions were independent prognostic markers for breast cancer. RING1 co-expressed genes were selected and included HDAC10, PIN1, CDK3, BAX, and BAD. GO analysis and KEGG pathway analyses revealed that RING1 related genes, were mainly enriched in "regulation of transcription", "apoptotic process", "protein transport", "protein binding", "Notch signaling pathway", and "Homologous recombination".

Conclusion

RING1 expression was downregulated in breast cancer, and its low expression was associated with worse disease outcomes. RING1 may act as a new prognostic biomarker for breast cancer.

Nongenetic Evolution Drives Lung Adenocarcinoma Spatial Heterogeneity and Progression

Cancer evolution determines molecular and morphologic intratumor heterogeneity and challenges the design of effective treatments. In lung adenocarcinoma, disease progression and prognosis are associated with the appearance of morphologically diverse tumor regions, termed histologic patterns. However, the link between molecular and histologic features remains elusive. Here, we generated multiomics and spatially resolved molecular profiles of histologic patterns from primary lung adenocarcinoma, which we integrated with molecular data from >2,000 patients. The transition from indolent to aggressive patterns was not driven by genetic alterations but by epigenetic and transcriptional reprogramming reshaping cancer cell identity. A signature quantifying this transition was an independent predictor of patient prognosis in multiple human cohorts. Within individual tumors, highly multiplexed protein spatial profiling revealed coexistence of immune desert, inflamed, and excluded regions, which matched histologic pattern composition. Our results provide a detailed molecular map of lung adenocarcinoma intratumor spatial heterogeneity, tracing nongenetic routes of cancer evolution. SIGNIFICANCE: Lung adenocarcinomas are classified based on histologic pattern prevalence. However, individual tumors exhibit multiple patterns with unknown molecular features. We characterized nongenetic mechanisms underlying intratumor patterns and molecular markers predicting patient prognosis. Intratumor patterns determined diverse immune microenvironments, warranting their study in the context of current immunotherapies.This article is highlighted in the In This Issue feature, p. 1307.

Identification of an Immune Gene-Associated Prognostic Signature and Its Association With a Poor Prognosis in Gastric Cancer Patients

The immune response plays a critical role in gastric cancer (GC) development, metastasis, and treatment. A better understanding of the tumor-immune system interactions in gastric cancer may provide promising diagnostic, prognostic, and therapeutic biomarkers for patients with this disease. In the present study, we aimed to identify a prognostic signature of GC through a comprehensive bioinformatics analysis on the tumor-immune interactions as well as the molecular characteristics. We firstly identified two immunophenotypes and immunological characteristics by employing multiple algorithms, such as the single sample Gene Sets Enrichment Analysis and Cell type Identification By Estimating Relative Subsets of RNA Transcripts. Next, we developed a six-immune-gene signature as a promising independent prognostic biomarker for GC using Lasso Cox regression and verified it via the external validation set and systematically correlated the immune signature with GC clinicopathologic features and genomic characteristics. Finally, a nomogram was successfully constructed based on the immune signature and clinical characteristics and showed a high potential for GC prognosis prediction. This study may shed light on the treatment strategies for GC patients from the perspective of immunology.

CTHRC1 Is a Prognostic Biomarker and Correlated With Immune Infiltrates in Kidney Renal Papillary Cell Carcinoma and Kidney Renal Clear Cell Carcinoma

Kidney renal clear cell carcinoma (KIRC) and kidney renal papillary cell carcinoma (KIRP) are the most common RCC types. RCC has high immune infiltration levels, and immunotherapy is currently one of the most promising treatments for RCC. Collagen triple helix repeat containing 1 (CTHRC1) is an extracellular matrix protein that regulates tumor invasion and modulates the tumor microenvironment. However, the association of CTHRC1 with the prognosis and tumor-infiltrating lymphocytes of KIRP and KIRC has not been reported. We examined the CTHRC1 expression differences in multiple tumor tissues and normal tissues via exploring TIMER, Oncomine, and UALCAN databases. Then, we searched the Kaplan-Meier plotter database to evaluate the correlation of CTHRC1 mRNA level with clinical outcomes. Subsequently, the TIMER platform and TISIDB website were chosen to assess the correlation of CTHRC1 with tumor immune cell infiltration level. We further explored the causes of aberrant CTHRC1 expression in tumorigenesis. We found that CTHRC1 level was significantly elevated in KIRP and KIRC tissues relative to normal tissues. CTHRC1 expression associates with tumor stage, histology, lymph node metastasis, and poor clinical prognosis in KIRP. The CTHRC1 level correlates to tumor grade, stage, nodal metastasis, and worse survival prognosis. Additionally, CTHRC1 is positively related to different tumor-infiltrating immune cells in KIRP and KIRC. Moreover, CTHRC1 was closely correlated with the gene markers of diverse immune cells. Also, high CTHRC1 expression predicted a worse prognosis in KIRP and KIRC based on immune cells. Copy number variations (CNV) and DNA methylation might contribute to the abnormal upregulation of CTHRC1 in KIRP and KIRC. In conclusion, CTHRC1 can serve as a biomarker to predict the prognosis and immune infiltration in KIRP and KIRC.

MiR-21 Is Required for the Epithelial-Mesenchymal Transition in MDA-MB-231 Breast Cancer Cells

Breast cancer (BCa) is one of the leading health problems among women. Although significant achievements have led to advanced therapeutic success with targeted therapy options, more efforts are required for different subtypes of tumors and according to genomic, transcriptomic, and proteomic alterations. This study underlines the role of microRNA-21 (miR-21) in metastatic MDA-MB-231 breast cancer cells. Following the knockout of miR-21 from MDA-MB-231 cells, which have the highest miR-21 expression levels compared to MCF-7 and SK-BR-3 BCa cells, a decrease in epithelial-mesenchymal transition (EMT) via downregulation of mesenchymal markers was observed. Wnt-11 was a critical target for miR-21, and the Wnt-11 related signaling axis was altered in the stable miR-21 knockout cells. miR-21 expression was associated with a significant increase in mesenchymal markers in MDA-MB-231 BCa cells. Furthermore, the release of extracellular vesicles (EVs) was significantly reduced in the miR-21 KO cells, alongside a significant reduction in relative miR-21 export in EV cargo, compared with control cells. We conclude that miR-21 is a leading factor involved in mesenchymal transition in MDA-MB-231 BCa. Future therapeutic strategies could focus on its role in the treatment of metastatic breast cancer.

A Risk Prediction Model for Breast Cancer Based on Immune Genes Related to Early Growth Response Proteins Family

Early growth response proteins (EGRs), a transcriptional regulatory family comprised of EGR1, EGR2, EGR3, and EGR 4, are reportedly involved in a vast array of functions. However, EGRs, as a whole, are rarely studied in breast cancer cases. This research was performed based on public datasets. The results demonstrated that, except EGR4, the other EGRs were differentially expressed genes in breast cancer. Subsequently, this study determined the prognosis significance of the EGR family, higher expression levels of EGRs indicating better overall survival (OS) and disease-free survival (DFS), except EGR4. So we attempted to explore the potential mechanism behind the prognostic value of EGRs. At the DNA level, however, neither DNA methylation status nor genetic alterations of EGRs contributed to the prognosis significance. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that EGRs were involved in several immune-related functions. Afterward, we assessed the correlation between EGRs and the immune system before establishing a risk prediction model with a 14-gene immune signature associated with EGRs, a prognostic nomogram predicting individuals’ 1-, 3-, and 5-year survival probabilities. The risk score was an independent prognosis predictor in the breast cancer cohorts. This study evidenced EGRs’ significance for tumor immunity, demonstrating that the EGR family may be a potential immunotherapeutic target for breast cancer. The 14-gene immune signature is a promising prognostic biomarker in breast cancer.

EPHA7 mutation as a predictive biomarker for immune checkpoint inhibitors in multiple cancers

BACKGROUND

A critical and challenging process in immunotherapy is to identify cancer patients who could benefit from immune checkpoint inhibitors (ICIs). Exploration of predictive biomarkers could help to maximize the clinical benefits. Eph receptors have been shown to play essential roles in tumor immunity. However, the association between EPH gene mutation and ICI response is lacking.

METHODS

Clinical data and whole-exome sequencing (WES) data from published studies were collected and consolidated as a discovery cohort to analyze the association between EPH gene mutation and efficacy of ICI therapy. Another independent cohort from Memorial Sloan Kettering Cancer Center (MSKCC) was adopted to validate our findings. The Cancer Genome Atlas (TCGA) cohort was used to perform anti-tumor immunity and pathway enrichment analysis.

RESULTS

Among fourteen EPH genes, EPHA7-mutant (EPHA7-MUT) was enriched in patients responding to ICI therapy (FDR adjusted P < 0.05). In the discovery cohort (n = 386), significant differences were detected between EPHA7-MUT and EPHA7-wildtype (EPHA7-WT) patients regarding objective response rate (ORR, 52.6% vs 29.1%, FDR adjusted P = 0.0357) and durable clinical benefit (DCB, 70.3% vs 42.7%, FDR adjusted P = 0.0200). In the validation cohort (n = 1144), significant overall survival advantage was observed in EPHA7-MUT patients (HR = 0.62 [95% confidence interval, 0.39 to 0.97], multivariable adjusted P = 0.0367), which was independent of tumor mutational burden (TMB) and copy number alteration (CNA). Notably, EPHA7-MUT patients without ICI therapy had significantly worse overall survival in TCGA cohort (HR = 1.33 [95% confidence interval, 1.06 to 1.67], multivariable adjusted P = 0.0139). Further gene set enrichment analysis revealed enhanced anti-tumor immunity in EPHA7-MUT tumor.

CONCLUSIONS

EPHA7-MUT successfully predicted better clinical outcomes in ICI-treated patients across multiple cancer types, indicating that EPHA7-MUT could serve as a potential predictive biomarker for immune checkpoint inhibitors.

GPRC5A: An emerging prognostic biomarker for predicting malignancy of Pancreatic Cancer based on bioinformatics analysis

Background: Pancreatic cancer (PaCa) is a highly lethal malignancy. The treatment options for PaCa lack efficacy. The study aimed to explore the molecular biomarkers for predicting survival of PaCa and identify the potential carcinogenic mechanisms of the selected gene. Methods: Based on public databases of PaCa, differentially expressed genes (DEGs) were identified using Networkanalyst. Survival analyses were exerted on GEPIA. Oncomine and The Human Protein Atlas were used for verifying the expression on mRNA and protein levels. Enrichment analyses were generated on Metascape and gene set enrichment analysis (GSEA). Univariate analyses were performed to determine the clinical factors associated with the expression of GPRC5A. Results: GPRC5A was identified as the most valuable gene in predicting survival of PaCa patients. Patients with high expression of GPRC5A showed larger tumor size, higher TNM stages, higher tumor grade, and more positive resection margin. In mutant KRAS, TP53, CDKN2A and SMAD4 group, the expression of GPRC5A was higher than non-mutant group. Mechanistically, GPRC5A may promote metastasis of PaCa mainly via regulating epithelial-mesenchymal transition (EMT) and neuroactive ligand-receptor interaction. Conclusion: GPRC5A may act as an oncogene in the progression of PaCa and could be a prognostic biomarker in predicting survival of PaCa.

The expression patterns and prognostic significance of pleckstrin homology-like domain family A (PHLDA) in lung cancer and malignant mesothelioma

Background

Pleckstrin homology domain family A (PHLDA) genes play important roles in cancer cellular processes, including inhibiting Akt activation, repressing growth factor signaling, inhibiting the negative feedback of EGFR/ErbB2 signaling cells, and inducing apoptosis. However, the prognostic significance of PHLDA in non-small cell lung cancer (NSCLC) and malignant pleural mesothelioma (MM) remains unclear. The present study investigates the associations between PHLDA expression patterns and their prognostic value in lung adenocarcinoma (LUAD) and MM.

Methods

We analyzed PHLDA family members at the genomic level in silico to explore their mRNA expression pattern and predictive significance in LUAD and MM. We then created a PHLDA–drug interaction network and a protein-protein interaction (PPI) network using different databases. Finally, we immunohistochemically assessed the protein expression of each PHLDA family member on tissue microarrays (TMAs) in both LUAD and MM cohorts with long-term follow-up.

Results

While PHLDA1 mRNA expression in both LUAD and MM was lower than that of normal tissue, PHLDA2 mRNA was significantly overexpressed in LUAD, and PHLDA3 mRNA was overexpressed in MM. In NSCLC, both low PHLDA1 mRNA expression and high PHLDA3 mRNA expression correlated with worse overall survival (OS) (P<0.01), whereas high PHLDA2 mRNA expression was associated with better OS (P<0.01). In MM, patients presenting high PHLDA1 and PHLDA2 mRNA expression had poor OS (P=0.01 and P<0.01, respectively). In addition, the PHLDA-drug interaction network indicated that several common drugs could potentially modulate PHLDA expression, and the PPI network suggested that PHLDA1 interacts with Notch family members, whereas PHLDA3 interacts with TP53. Our results also showed that the expression of PHLDA2 and PHLDA3 was significantly higher in LUAD and MM than that of PHLDA1 (P<0.05) and was associated with the risk of death. While patients with PHLDA2 >85.09 cells/mm² had a low risk of death (P=0.01) and a median survival time of 48 months, those with PHLDA3 <70.38 cells/mm² had a high risk of death (P=0.03) and a median survival time of 34 months.

Conclusions

We shed light on the role of the PHLDA family as promising predictive biomarkers and potential therapeutic targets in LUAD and MM.

Exploring the Differential Expression and Prognostic Significance of the COL11A1 Gene in Human Colorectal Carcinoma: An Integrated Bioinformatics Approach

Colorectal cancer is one of the most common cancers of humans and the second highest in cancer-related death. Genes used as prognostic biomarkers play an imperative role in cancer detection and may direct the development of appropriate therapeutic strategies. Collagen type XI alpha 1 (COL11A1) is a minor fibrillary collagen that has an essential role in the regulation of cell division, differentiation, proliferation, migration, growth, and apoptosis of intestinal and colon cells. The present study seeks to evaluate the significance of the COL11A1 gene in the progression of colorectal cancer in humans across the various parameters using advanced bioinformatics approaches. The application of various databases and servers like ONCOMINE, UALCAN, and GEPIA were accessed for analyzing the differential expression of the COLL11A1 gene and its relative influence over the survival of the transformed subjects. In addition, oncogenomics of COL11A1 gene, mutations associated with this gene and interacting partners of the gene in the context of oncogenesis were studied using COSMIC, cBioPortal, GeneMANIA, and NetworkAnalyst. Our experimental data indicate that the COL11A1 gene is overexpressed in the transformed tissues across the various clinicopathological parameters reduces the probability of survival in both overall and disease-specific survival cases. Mutational studies imply that it can induce perturbations in various signaling pathways viz. RTK-RAS-PI3K, Wnt, TGF-β, and TP53 pathways influencing cancer development. Also, a positive association and correlation amongst the THBS2, COL10A1, COL5A2, and COL1A2 genes were observed, which most likely to contribute to the upregulation of carcinogenesis. Conclusively, this comprehensive study indicates the COL11A1 gene to be a significant contributor in the etiology of colorectal cancer, henceforth this gene can be considered as a prognostic biomarker for the conception of diagnostic and therapeutic strategies against colorectal cancer in the near future.

RNA m6A Methylation Regulators Subclassify Luminal Subtype in Breast Cancer

RNA N⁶-methyladenosine (m⁶A) methylation is the most prevalent epitranscriptomic modification in mammals, with a complex and fine-tuning regulatory system. Recent studies have illuminated the potential of m⁶A regulators in clinical applications including diagnosis, therapeutics, and prognosis. Based on six datasets of breast cancer in The Cancer Genome Atlas (TCGA) database and two additional proteomic datasets, we provide a comprehensive view of all the known m⁶A regulators in their gene expression, copy number variations (CNVs), DNA methylation status, and protein levels in breast tumors and their association with prognosis. Among four breast cancer subtypes, basal-like subtype exhibits distinct expression and genomic alteration in m⁶A regulators from other subtypes. Accordingly, four representative regulators (IGF2BP2, IGF2BP3, YTHDC2, and RBM15) are identified as basal-like subtype-featured genes. Notably, luminal A/B samples are subclassified into two clusters based on the methylation status of those four genes. In line with its similarity to basal-like subtype, cluster1 shows upregulation in immune-related genes and cell adhesion molecules, as well as an increased number of tumor-infiltrating lymphocytes. Besides, cluster1 has worse disease-free and progression-free survival, especially among patients diagnosed with stage II and luminal B subtype. Together, this study highlights the potential functions of m⁶A regulators in the occurrence and malignancy progression of breast cancer. Given the heterogeneity within luminal subtype and high risk of recurrence and metastasis in a portion of patients, the prognostic stratification of luminal A/B subtypes utilizing basal-featured m⁶A regulators may help to improve the accuracy of diagnosis and therapeutics of breast cancer.

MethCORR infers gene expression from DNA methylation and allows molecular analysis of ten common cancer types using fresh-frozen and formalin-fixed paraffin-embedded tumor samples

Background

Transcriptional analysis is widely used to study the molecular biology of cancer and hold great biomarker potential for clinical patient stratification. Yet, accurate transcriptional profiling requires RNA of a high quality, which often cannot be retrieved from formalin-fixed, paraffin-embedded (FFPE) tumor tissue that is routinely collected and archived in clinical departments. To overcome this roadblock to clinical testing, we previously developed MethCORR, a method that infers gene expression from DNA methylation data, which is robustly retrieved from FFPE tissue. MethCORR was originally developed for colorectal cancer and with this study, we aim to: (1) extend the MethCORR method to 10 additional cancer types and (2) to illustrate that the inferred gene expression is accurate and clinically informative.

Results

Regression models to infer gene expression information from DNA methylation were developed for ten common cancer types using matched RNA sequencing and DNA methylation profiles (HumanMethylation450 BeadChip) from The Cancer Genome Atlas Project. Robust and accurate gene expression profiles were inferred for all cancer types: on average, the expression of 11,000 genes was modeled with good accuracy and an intra-sample correlation of R² = 0.90 between inferred and measured gene expression was observed. Molecular pathway analysis and transcriptional subtyping were performed for breast, prostate, and lung cancer samples to illustrate the general usability of the inferred gene expression profiles: overall, a high correlation of r = 0.96 (Pearson) in pathway enrichment scores and a 76% correspondence in molecular subtype calls were observed when using measured and inferred gene expression as input. Finally, inferred expression from FFPE tissue correlated better with RNA sequencing data from matched fresh-frozen tissue than did RNA sequencing data from FFPE tissue (P < 0.0001; Wilcoxon rank-sum test).

Conclusions

In all cancers investigated, MethCORR enabled DNA methylation-based transcriptional analysis, thus enabling future analysis of cancer in situations where high-quality DNA, but not RNA, is available. Here, we provide the framework and resources for MethCORR modeling of ten common cancer types, thereby widely expanding the possibilities for transcriptional studies of archival FFPE material.

Five Circular RNAs in Metabolism Pathways Related to Prostate Cancer

Prostate cancer (PCa) is the most common malignant tumor in men, and its incidence increases with age. Serum prostate-specific antigen and tissue biopsy remain the standard for diagnosis of suspected PCa. However, these clinical indicators may lead to aggressive overtreatment in patients who have been treated sufficiently with active surveillance. Circular RNAs (circRNAs) have been recently recognized as a new type of regulatory RNA that is not easily degraded by RNases and other exonucleases because of their covalent closed cyclic structure. Thus, we utilized high-throughput sequencing data and bioinformatics analysis to identify specifically expressed circRNAs in PCa and filtered out five specific circRNAs for further analysis—hsa_circ_0006410, hsa_circ_0003970, hsa_circ_0006754, hsa_circ_0005848, and a novel circRNA, hsa_circ_AKAP7. We constructed a circRNA-miRNA regulatory network and used miRNA and differentially expressed mRNA interactions to predict the function of the selected circRNAs. Furthermore, survival analysis of their cognate genes and PCR verification of these five circRNAs revealed that they are closely related to well-known PCa pathways such as the MAPK signaling pathway, P53 pathway, androgen receptor signaling pathway, cell cycle, hormone-mediated signaling pathway, and cellular lipid metabolic process. By understanding the related metabolism of circRNAs, these circRNAs could act as metabolic biomarkers, and monitoring their levels could help diagnose PCa. Meanwhile, the exact regulatory mechanism for AR-related regulation in PCa is still unclear. The circRNAs we found can provide new solutions for research in this field.

DIGGER: exploring the functional role of alternative splicing in protein interactions

Alternative splicing plays a major role in regulating the functional repertoire of the proteome. However, isoform-specific effects to protein-protein interactions (PPIs) are usually overlooked, making it impossible to judge the functional role of individual exons on a systems biology level. We overcome this barrier by integrating protein-protein interactions, domain-domain interactions and residue-level interactions information to lift exon expression analysis to a network level. Our user-friendly database DIGGER is available at https://exbio.wzw.tum.de/digger and allows users to seamlessly switch between isoform and exon-centric views of the interactome and to extract sub-networks of relevant isoforms, making it an essential resource for studying mechanistic consequences of alternative splicing.

High LARGE1 Expression May Predict Benefit from Adjuvant Chemotherapy in Resected Non-Small-Cell Lung Cancer

Background

LARGE1 plays a pivotal role in glycosylation of alpha-Dystroglycan (α-DG) and is aberrantly downregulated in cell lines originating from epithelium-derived cancers including lung cancer. However, the expression of LARGE1 and its clinical significance in NSCLC are not clear.

Materials and Methods

The data were collected from the TCGA database to investigate LARGE1 expression in stage I–III NSCLC and explore its associations with clinicopathological parameters and overall survival of patients. The prognostic role of LARGE1 was examined in subgroups according to clinical features and treatments. The results were validated in external cohorts from the NCBI GEO database. Gene Set Enrichment Analysis (GSEA) was performed to investigate the potential molecular mechanisms during LARGE1 alteration in NSCLC.

Results

LARGE1 was aberrantly downregulated in NSCLC compared with adjacent tissues and normal lung tissues and in tumors with advanced stage compared with early stage. There was only a trend of association between high LARGE1 with OS in multivariate analysis. Surprisingly, high LARGE1 was significantly associated with improved OS in a subgroup of the patients with adjuvant chemotherapy (ACT) and a significant interaction between LARGE1 expression and ACT was found. Improved OS after ACT was also found in patients with high LARGE1 compared to those with low LARGE1. When combining LARGE1 expression and ACT, compared with patients with non-ACT, HR of low LARGE1/ACT was 0.592 (95% CI=0.432–0.813, P=0.0012), and HR of high LARGE1/ACT was 0.124 (95% CI=0.031–0.505, P=0.0036). The results were verified in two external cohorts from the GEO database. GSEA indicated that LARGE1 might downregulate cell cycle pathway to improve ACT sensitivity and subsequently the prognosis in NSCLC.

Conclusion

High LARGE1 can be used to identify the patients with resected stage I–III NSCLC most likely to benefit from adjuvant chemotherapy.

Oncometabolite L-2-hydroxyglurate directly induces vasculogenic mimicry through PHLDB2 in renal cell carcinoma

Metabolism reprograming is a hallmark of cancer and plays an important role in tumor progression. The aberrant metabolism in renal cell carcinoma (RCC) leads to accumulation of the oncometabolite l‐2‐hydroxyglurate (L‐2HG). L‐2HG has been reported to inhibit the activity of some α‐ketoglutarate‐dependent dioxygenases such as TET enzymes, which mediate epigenetic alteration, including DNA and histone demethylation. However, the detailed functions of L‐2HG in renal cell carcinoma have not been investigated thoroughly. In our study, we found that L‐2HG was significantly elevated in tumor tissues compared to adjacent tissues. Furthermore, we demonstrated that L‐2HG promoted vasculogenic mimicry (VM) in renal cancer cell lines through reducing the expression of PHLDB2. A mechanism study revealed that activation of the ERK1/2 pathway was involved in L‐2HG‐induced VM formation. In conclusion, these findings highlighted the pathogenic link between L‐2HG and VM and suggested a novel therapeutic target for RCC.

What's new?

Metabolic reprograming, a hallmark of cancer, influences tumor progression. In the case of renal cell carcinoma (RCC) specifically, progression appears to be facilitated by the oncometabolite L‐2‐hydroxyglurate (L‐2HG), though underlying mechanisms remain enigmatic. Here, the authors investigated the ability of L‐2HG in RCC to promote vasculogenic mimicry (VM), in which aggressive cancer cells form vessel‐like networks that support tumor growth. Analyses of RCC patient tissues revealed elevated L‐2HG levels, wherein tumor cells with greater L‐2HG levels exhibited more VM structures. TCGA data and high‐throughput sequencing analyses further show that L‐2HG contributes to VM formation via reduction of PHLDB2 levels.

MuSA: a graphical user interface for multi-OMICs data integration in radiogenomic studies

Analysis of large-scale omics data along with biomedical images has gaining a huge interest in predicting phenotypic conditions towards personalized medicine. Multiple layers of investigations such as genomics, transcriptomics and proteomics, have led to high dimensionality and heterogeneity of data. Multi-omics data integration can provide meaningful contribution to early diagnosis and an accurate estimate of prognosis and treatment in cancer. Some multi-layer data structures have been developed to integrate multi-omics biological information, but none of these has been developed and evaluated to include radiomic data. We proposed to use MultiAssayExperiment (MAE) as an integrated data structure to combine multi-omics data facilitating the exploration of heterogeneous data. We improved the usability of the MAE, developing a Multi-omics Statistical Approaches (MuSA) tool that uses a Shiny graphical user interface, able to simplify the management and the analysis of radiogenomic datasets. The capabilities of MuSA were shown using public breast cancer datasets from TCGA-TCIA databases. MuSA architecture is modular and can be divided in Pre-processing and Downstream analysis. The pre-processing section allows data filtering and normalization. The downstream analysis section contains modules for data science such as correlation, clustering (i.e., heatmap) and feature selection methods. The results are dynamically shown in MuSA. MuSA tool provides an easy-to-use way to create, manage and analyze radiogenomic data. The application is specifically designed to guide no-programmer researchers through different computational steps. Integration analysis is implemented in a modular structure, making MuSA an easily expansible open-source software.

Synergism of Proneurogenic miRNAs Provides a More Effective Strategy to Target Glioma Stem Cells

Simple Summary

miRNAs function as critical regulators of gene expression and have been defined as contributors of cancer phenotypes by acting as oncogenes or tumor suppressors. Based on these findings, miRNA-based therapies have been explored in the treatment of many different malignancies. The use of single miRNAs has faced some challenges and showed limited success. miRNAs cooperate to regulate distinct biological processes and pathways and, therefore, combination of related miRNAs could amplify the repression of oncogenic factors and the effect on cancer relevant pathways. We established that the combination of tumor suppressor miRNAs miR-124, miR-128, and miR-137 is much more effective than single miRNAs in disrupting proliferation and survival of glioma stem cells and neuroblastoma lines and promoting differentiation and response to radiation. Subsequent genomic analyses showed that other combinations of tumor suppressor miRNAs could be equally effective, and its use could provide new routes to target in special cancer-initiating cell populations.

Abstract

Tumor suppressor microRNAs (miRNAs) have been explored as agents to target cancer stem cells. Most strategies use a single miRNA mimic and present many disadvantages, such as the amount of reagent required and the diluted effect on target genes. miRNAs work in a cooperative fashion to regulate distinct biological processes and pathways. Therefore, we propose that miRNA combinations could provide more efficient ways to target cancer stem cells. We have previously shown that miR-124, miR-128, and miR-137 function synergistically to regulate neurogenesis. We used a combination of these three miRNAs to treat glioma stem cells and showed that this treatment was much more effective than single miRNAs in disrupting cell proliferation and survival and promoting differentiation and response to radiation. Transcriptomic analyses indicated that transcription regulation, angiogenesis, metabolism, and neuronal differentiation are among the main biological processes affected by transfection of this miRNA combination. In conclusion, we demonstrated the value of using combinations of neurogenic miRNAs to disrupt cancer phenotypes and glioma stem cell growth. The synergistic effect of these three miRNA amplified the repression of oncogenic factors and the effect on cancer relevant pathways. Future therapeutic approaches would benefit from utilizing miRNA combinations, especially when targeting cancer-initiating cell populations.

Serine-70 phosphorylated Bcl-2 prevents oxidative stress-induced DNA damage by modulating the mitochondrial redox metabolism

Bcl-2 phosphorylation at serine-70 (S70pBcl2) confers resistance against drug-induced apoptosis. Nevertheless, its specific mechanism in driving drug-resistance remains unclear. We present evidence that S70pBcl2 promotes cancer cell survival by acting as a redox sensor and modulator to prevent oxidative stress-induced DNA damage and execution. Increased S70pBcl2 levels are inversely correlated with DNA damage in chronic lymphocytic leukemia (CLL) and lymphoma patient-derived primary cells as well as in reactive oxygen species (ROS)- or chemotherapeutic drug-treated cell lines. Bioinformatic analyses suggest that S70pBcl2 is associated with lower median overall survival in lymphoma patients. Empirically, sustained expression of the redox-sensitive S70pBcl2 prevents oxidative stress-induced DNA damage and cell death by suppressing mitochondrial ROS production. Using cell lines and lymphoma primary cells, we further demonstrate that S70pBcl2 reduces the interaction of Bcl-2 with the mitochondrial complex-IV subunit-5A, thereby reducing mitochondrial complex-IV activity, respiration and ROS production. Notably, targeting S70pBcl2 with the phosphatase activator, FTY720, is accompanied by an enhanced drug-induced DNA damage and cell death in CLL primary cells. Collectively, we provide a novel facet of the anti-apoptotic Bcl-2 by demonstrating that its phosphorylation at serine-70 functions as a redox sensor to prevent drug-induced oxidative stress-mediated DNA damage and execution with potential therapeutic implications.

Exploration of Coding and Non-coding Variants in Cancer Using GenomePaint

GenomePaint (https://genomepaint.stjude.cloud/) is an interactive visualization platform for whole-genome, whole-exome, transcriptome, and epigenomic data of tumor samples. Its design captures the inter-relatedness between DNA variations and RNA expression, supporting in-depth exploration of both individual cancer genomes and full cohorts. Regulatory non-coding variants can be inspected and analyzed along with coding variants, and their functional impact further explored by examining 3D genome data from cancer cell lines. Further, GenomePaint correlates mutation and expression patterns with patient outcomes, and supports custom data upload. We used GenomePaint to unveil aberrant splicing that disrupts the RING domain of CREBBP, discover cis activation of the MYC oncogene by duplication of the NOTCH1-MYC enhancer in B-lineage acute lymphoblastic leukemia, and explore the inter- and intra-tumor heterogeneity at EGFR in adult glioblastomas. These examples demonstrate that deep multi-omics exploration of individual cancer genomes enabled by GenomePaint can lead to biological insights for follow-up validation.

GRNdb: decoding the gene regulatory networks in diverse human and mouse conditions

Gene regulatory networks (GRNs) formed by transcription factors (TFs) and their downstream target genes play essential roles in gene expression regulation. Moreover, GRNs can be dynamic changing across different conditions, which are crucial for understanding the underlying mechanisms of disease pathogenesis. However, no existing database provides comprehensive GRN information for various human and mouse normal tissues and diseases at the single-cell level. Based on the known TF-target relationships and the large-scale single-cell RNA-seq data collected from public databases as well as the bulk data of The Cancer Genome Atlas and the Genotype-Tissue Expression project, we systematically predicted the GRNs of 184 different physiological and pathological conditions of human and mouse involving >633 000 cells and >27 700 bulk samples. We further developed GRNdb, a freely accessible and user-friendly database (http://www.grndb.com/) for searching, comparing, browsing, visualizing, and downloading the predicted information of 77 746 GRNs, 19 687 841 TF-target pairs, and related binding motifs at single-cell/bulk resolution. GRNdb also allows users to explore the gene expression profile, correlations, and the associations between expression levels and the patient survival of diverse cancers. Overall, GRNdb provides a valuable and timely resource to the scientific community to elucidate the functions and mechanisms of gene expression regulation in various conditions.

Interferon Regulatory Factor 9 Promotes Lung Cancer Progression via Regulation of Versican

Simple Summary

Lung cancer is the leading cause of cancer-related deaths worldwide, accounting for more than 1.6 million deaths per year. The tumor microenvironment (TME) has been shown to play a crucial role in tumor progression and metastasis, and transcription factors link TME signaling to oncogenesis. Type I interferons (IFNs) are strong immune modulators that possess antiproliferative and proapoptotic properties. In this study, we investigated the role of the transcription factor interferon regulatory factor 9 (IRF9) in the IFN pathway in lung cancer. We performed in vitro and in vivo experiments to reveal the oncogenic properties of IRF9, which was highly upregulated in lung adenocarcinoma. For the first time, we showed that IRF9 binds to the promoter of the known oncogene versican, regulates its expression, and thereby promotes oncogenic activity.

Abstract

Transcription factors can serve as links between tumor microenvironment signaling and oncogenesis. Interferon regulatory factor 9 (IRF9) is recruited and expressed upon interferon stimulation and is dependent on cofactors that exert in tumor-suppressing or oncogenic functions via the JAK-STAT pathway. IRF9 is frequently overexpressed in human lung cancer and is associated with decreased patient survival; however, the underlying mechanisms remain to be elucidated. Here, we used stably transduced lung adenocarcinoma cell lines (A549 and A427) to overexpress or knockdown IRF9. Overexpression led to increased oncogenic behavior in vitro, including enhanced proliferation and migration, whereas knockdown reduced these effects. These findings were confirmed in vivo using lung tumor xenografts in nude mice, and effects on both tumor growth and tumor mass were observed. Using RNA sequencing, we identified versican (VCAN) as a novel downstream target of IRF9. Indeed, IRF9 and VCAN expression levels were found to be correlated. We showed for the first time that IRF9 binds at a newly identified response element in the promoter region of VCAN to regulate its transcription. Using an siRNA approach, VCAN was found to enable the oncogenic properties (proliferation and migration) of IRF9 transduced cells, perhaps with CDKN1A involvement. The targeted inhibition of IRF9 in lung cancer could therefore be used as a new treatment option without multimodal interference in microenvironment JAK-STAT signaling.

CD73, Tumor Plasticity and Immune Evasion in Solid Cancers

Simple Summary

Tumors are ecosystems composed of cancer cells and non-tumor stroma together in a hypoxic environment often described as wounds that do not heal. Accumulating data suggest that solid tumors hijack cellular plasticity possibly to evade detection by the immune system. CD73-mediated generation of the purine nucleoside adenosine, is an important biochemical constituent of the immunosuppressive tumor microenvironment. In this review, the association between CD73 expression and features associated with cellular plasticity involving stemness, epithelial-to-mesenchymal transition and metastasis together with immune infiltration is summarized for a wide range of solid tumor types. Our analyses demonstrate that CD73 correlates with signatures associated with cellular plasticity in solid tumors. In addition, there are strong associations between CD73 expression and type of infiltrating lymphocytes. Collectively, the observations suggest a biomarker-based stratification to identify CD73-adenosinergic rich tumors may help identify patients with solid cancers who will respond to a combinatorial strategy that includes targeting CD73.

Abstract

Regulatory networks controlling cellular plasticity, important during early development, can re-emerge after tissue injury and premalignant transformation. One such regulatory molecule is the cell surface ectoenzyme ecto-5′-nucleotidase that hydrolyzes the conversion of extracellular adenosine monophosphate to adenosine (eADO). Ecto-5′-nucleotidase (NT5E) or cluster of differentiation 73 (CD73), is an enzyme that is encoded by NT5E in humans. In normal tissue, CD73-mediated generation of eADO has important pleiotropic functions ranging from the promotion of cell growth and survival, to potent immunosuppression mediated through purinergic G protein-coupled adenosine receptors. Importantly, tumors also utilize several mechanisms mediated by CD73 to resist therapeutics and in particular, evade the host immune system, leading to undesired resistance to targeted therapy and immunotherapy. Tumor cell CD73 upregulation is associated with worse clinical outcomes in a variety of cancers. Emerging evidence indicates a link between tumor cell stemness with a limited host anti-tumor immune response. In this review, we provide an overview of a growing body of evidence supporting the pro-tumorigenic role of CD73 and adenosine signaling. We also discuss data that support a link between CD73 expression and tumor plasticity, contributing to dissemination as well as treatment resistance. Collectively, targeting CD73 may represent a novel treatment approach for solid cancers.

SPONGEdb: a pan-cancer resource for competing endogenous RNA interactions

microRNAs (miRNAs) are post-transcriptional regulators involved in many biological processes and human diseases, including cancer. The majority of transcripts compete over a limited pool of miRNAs, giving rise to a complex network of competing endogenous RNA (ceRNA) interactions. Currently, gene-regulatory networks focus mostly on transcription factor-mediated regulation, and dedicated efforts for charting ceRNA regulatory networks are scarce. Recently, it became possible to infer ceRNA interactions genome-wide from matched gene and miRNA expression data. Here, we inferred ceRNA regulatory networks for 22 cancer types and a pan-cancer ceRNA network based on data from The Cancer Genome Atlas. To make these networks accessible to the biomedical community, we present SPONGEdb, a database offering a user-friendly web interface to browse and visualize ceRNA interactions and an application programming interface accessible by accompanying R and Python packages. SPONGEdb allows researchers to identify potent ceRNA regulators via network centrality measures and to assess their potential as cancer biomarkers through survival, cancer hallmark and gene set enrichment analysis. In summary, SPONGEdb is a feature-rich web resource supporting the community in studying ceRNA regulation within and across cancer types.

Sex differences in expression of immune elements emerge in children, young adults and mice with osteosarcoma

BACKGROUND

Males < 40 years old are more likely to be diagnosed with and die from osteosarcoma (OS). The underlying mechanisms may depend on sex differences in immune response.

METHODS

We used SEER data to estimate survival differences between males and females aged < 40 years at OS diagnosis. In NCI TARGET-OS cases, we determined sex differences in gene expression, conducted Gene Set Enrichment Analysis (GSEA), and applied the LM22 signature to identify biologic sex differences. We compared sex differences in gene expression profiles in TARGET-OS to those observed in Sleeping Beauty (SB) transposon mutagenesis accelerated Trp53R270H-mutant mouse-OS and healthy adult osteoblasts.

RESULTS

Males had worse 17-year overall survival than females (SEER p < 0.0001). From 87 TARGET-OS cases, we observed 1018 genes and 69 pathways that differed significantly by sex (adjusted p < 0.05). Pathway and gene lists overlapped with those from mice (p = 0.03) and healthy osteoblasts (p = 0.017), respectively. Pathways that differed significantly by sex were largely immune-based and included the PD-1/PD-L1 immunotherapy pathway. We observed sex differences in M2 macrophages (LM22; p = 0.056) and M1-M2 macrophage transition (GSEA; p = 0.037) in TARGET-OS. LM22 trends were similar in mice. Twenty-four genes differentially expressed by sex in TARGET-OS had existing cancer therapies.

CONCLUSIONS

Sex differences in OS gene expression were similar across species and centered on immune pathways. Identified sex-specific therapeutic targets may improve outcomes in young individuals with OS.

Upregulation of TMEM45A Promoted the Progression of Clear Cell Renal Cell Carcinoma in vitro

BACKGROUND

Clear cell renal cell carcinoma (ccRCC) is the most common and aggressive type of primary kidney cancer worldwide. Transmembrane protein 45A (TMEM45A) has been reported to be closely associated with the progression of several cancers. However, the role of TMEM45A in ccRCC remains unclear. Our study intended to explore the potential role of TMEM45A in ccRCC.

METHODS

Data on the expression of TMEM45A were obtained from multiple databases, including UCSC, GEPIA2, Oncomine and TIMER. Real-world samples of ccRCC and paired normal renal tissues were used to confirm the information obtained from the databases. In addition, the prognostic value of TMEM45A was evaluated. Loss-of-function assays were performed using TMEM45A-targeting lentivirus to evaluate the biological role of TMEM45A in renal cancer cells. Gene set enrichment analysis (GSEA) was performed to investigate the potential molecular mechanisms.

RESULTS

TMEM45A was significantly overexpressed in patients with ccRCC and correlated with poor overall survival and disease-free survival. In addition, the expression of TMEM45A was closely associated with various clinicopathological parameters such as histological grade and TNM stage. Knockdown of TMEM45A inhibited the proliferation and migration and promoted the apoptosis of ccRCC cells in vitro. The results of the GSEA suggested that TMEM45A was potentially involved in the promotion of epithelial-mesenchymal transition (EMT) and inflammatory response in ccRCC.

CONCLUSION

TMEM45A was overexpressed and associated with poor survival and acted as a tumour promoter in ccRCC; therefore, might be a potential prognostic marker and therapeutic target.

PTGS2 Over-Expression: A Colorectal Carcinoma Initiator not an Invasive Factor

BACKGROUND

Cyclooxygenase-2 (COX-2) main product is Prostaglandin E2 (PGE2) which cause mitogenesis and inflammation. COX-2 is the product of prostaglandin-endoperoxide synthase 2 (PTGS2) gene expression. COX-2 dysregulation can cause angiogenesis, differentiation, and promotion of cancer and its suppression related to control of the tumor's size, number, and cell shape. This study focused on the association of COX-2 expression with colorectal carcinoma (CRC) among Iranian patients on mRNA level and in the Cancer Genome Atlas Program (TCGA) colon and rectum RNAseq dataset, and its relation with pathological features.

METHODS

PTGS2 expression was assayed by quantitative-PCR method from 90 tissue samples collected from 45 participants. The control samples come from the non-tumor area of the same patients. The data analyzed based on ΔΔCq. The PTGS2-RNAseq data extracted and analyzed by UCSC Xena browser, and its association assessed the occurrence of CRC and invasive-features.

RESULTS

PTGS2 showed very significant over-expression in tumor tissues (p< 0.0001) with an N-fold expression of 2.25. But, there was not any significant association between PTGS2 and CRC invasive-pathological features such as Lymphatic, vascular and perineural invasion, the Grades of cancer, and Pathologic-M in both parts of this study.

CONCLUSION

The increase in PTGS2 is related to the occurrence of CRC among patient samples. But in both part of this study, PTGS2 is not an invasive factor, and it does not affect the cell differentiation of tumors and metastasis. Based on the high N-fold for patient samples, it can be a strong candidate as a CRC initiator biomarker.

LINC00467 Is Upregulated by DNA Copy Number Amplification and Hypomethylation and Shows ceRNA Potential in Lung Adenocarcinoma

Lung adenocarcinoma (LUAD) is the most common histological lung cancer, and it is the leading cause of cancer-related deaths worldwide. Long noncoding RNAs (lncRNAs) have been implicated in tumorigenesis. LINC00467 is a novel lncRNA that is abnormally expressed in several cancer types including LUAD. However, its function and regulatory mechanism in LUAD progression remain unclear. In this study, based on The Cancer Genome Atlas data mining, we demonstrated that DNA copy number amplification and hypomethylation was positively correlated with LINC00467 expression in LUAD. In addition, DNA copy number amplification was significantly associated with distant metastasis, immune infiltration and poor survival. Microarray analysis demonstrated that LINC00467 knockdown in the LUAD A549 cell line led to a distinct microRNA expression profile that impacted various target genes involved in multiple biological processes. This finding suggests that LINC00467 may regulate LUAD progression by functioning as a competing endogenous RNA (ceRNA). Finally, we constructed a ceRNA network that included two microRNAs (hsa-miR-1225-5p, hsa-miR-575) and five mRNAs (BARX2, BCL9, KCNK1, KIAA1324, TMEM182) specific to LINC00467 in LUAD. Subsequent Kaplan-Meier survival analysis in both The Cancer Genome Atlas and Gene Expression Omnibus databases revealed that two genes, BARX2 and BCL9, were potential prognostic biomarkers for LUAD patients. In conclusion, our data provide possible mechanisms underlying the abnormal upregulation of LINC00467 as well as a comprehensive view of the LINC00467-mediated ceRNA network in LUAD, thereby highlighting its potential role in diagnosis and therapy.

Immunological Role and Prognostic Value of APBB1IP in Pan-Cancer Analysis

Objective: APBB1IP is a Rap1-binding protein that mainly acts as a regulator of leukocyte recruitment and pathogen clearance through complement-mediated phagocytosis. However, the role of APBB1IP in tumor immunity remains unclear. This study was carried out to evaluate the prognostic landscape of APBB1IP in pan-cancer analysis and investigate the relationship between APBB1IP expression and immune infiltration. Methods: We explored the expression pattern and prognostic value of APBB1IP in pan-cancer analysis through Kaplan-Meier Plotter and multiple databases, including TCGA, Oncomine. We then assessed the correlation between APBB1IP expression and immune cell infiltration using the TIMER database. Furthermore, we identified the proteins that interact with APBB1IP and performed epigenetic and transcriptional analyses. Multivariate Cox regression analyses were applied to construct a prognostic model, which consisted of APBB1IP and its interacting proteins, based on the lung cancer cohorts from the Gene Expression Omnibus (GEO) database. Results: The expression of APBB1IP was correlated with the prognosis of several types of cancer. APBB1IP upregulation was found to be associated with increased immune cell infiltration, especially for CD8+ T cells, natural killer (NK) cells, and immune regulators. A link was found between APBB1IP and immune-related proteins including RAP1A/B, TLN1/2 and VCL in the interaction network. Conclusion: APBB1IP can serve as a prognostic biomarker in pan-cancer analysis. APBB1IP upregulation was correlated with increased immune-cell infiltration, and the expression APBB1IP in different tumors might be related to the tumor immune microenvironment.

Super enhancers define regulatory subtypes and cell identity in neuroblastoma

Half of the children diagnosed with neuroblastoma (NB) have high-risk disease, disproportionately contributing to overall childhood cancer-related deaths. In addition to recurrent gene mutations, there is increasing evidence supporting the role of epigenetic deregulation in disease pathogenesis. Yet, comprehensive cis-regulatory network descriptions from NB are lacking. Here, using genome-wide H3K27ac profiles across 60 NBs, covering the different clinical and molecular subtypes, we identified four major super-enhancer-driven epigenetic subtypes and their underlying master regulatory networks. Three of these subtypes recapitulated known clinical groups; namely, MYCN-amplified, MYCN non-amplified high-risk and MYCN non-amplified low-risk NBs. The fourth subtype, exhibiting mesenchymal characteristics, shared cellular identity with multipotent Schwann cell precursors, was induced by RAS activation and was enriched in relapsed disease. Notably, CCND1, an essential gene in NB, was regulated by both mesenchymal and adrenergic regulatory networks converging on distinct super-enhancer modules. Overall, this study reveals subtype-specific super-enhancer regulation in NBs.

LINC00467, Driven by Copy Number Amplification and DNA Demethylation, Is Associated with Oxidative Lipid Metabolism and Immune Infiltration in Breast Cancer

Breast cancer (BRCA) is a malignant tumor with a high incidence and poor prognosis in females. However, its pathogenesis remains unclear. In this study, based on bioinformatic analysis, we found that LINC00467 was highly expressed in BRCA and was associated with tumor metastasis and poor prognosis. The genomic and epigenetic analysis showed that LINC00467 may also be regulated by copy number amplification (CNA), chromatin openness, and DNA methylation. In vitro experiments showed that it could promote the proliferation, migration, and invasion of BRCA cells. Competitive endogenous RNA (ceRNA) regulatory network analysis and weighted gene coexpression network analysis (WGCNA) suggested that LINC00467 may play a role in signaling pathways of peroxisomal lipid metabolism, immunity, and others through microRNAs (miRNAs) targeting transforming growth factor beta 2 (TGFB2). In addition, copy number amplification and high expression of LINC00467 were associated with the low infiltration of CD8+ and CD4+ T cells. In conclusion, we found that LINC00467, driven by copy number amplification and DNA demethylation, may be a potential biomarker for the diagnosis and prognosis of BRCA and a tumor promoter acting as a potential therapeutic target for BRCA as well.