Pan-cancer analyses of the nuclear receptor superfamily

NR3C2 affects the proliferation and invasiveness of colon cancer cells through the Wnt/β-Catenin signaling pathway

PURPOSE

The aim of this study was to explore the potential correlation between the nuclear receptor subfamily 3 group C member 2 (NR3C2) and outcomes of colon cancer, along with the mechanisms underlying this association.

METHOD

mRNA (messenger RNA) data and clinical records pertaining to colon cancer were retrieved from The Cancer Genome Atlas (TCGA) database. The analysis of NR3C2 expression discrepancies between normal colon and tumor tissues was conducted using R software. In addition, we also studied the relationship between NR3C2 expression and prognosis, pathological parameters. The relative role of NR3C2 were further predicted through bioinformatics methods and receiver operating characteristic (ROC) curve was used to evaluate the diagnostic value of NR3C2 in colon cancer. Single-cell data from colon cancer samples in the GEO (Gene Expression Omnibus) database further investigated the mechanism of the lower survival associated with NR3C2 dysregulation. NR3C2 expression in three fresh colon cancer samples and their respective paracancer samples was determined. Furthermore, colon cancer cell models overexpressing NR3C2 and with knockdown NR3C2 were constructed by lentiviral vector transfection. Cell Counting Kit-8 assay, transplantation of tumors in nude mice and transwell assays were used to examine the proliferation, migration and invasion of colon cancer cells. The effect on the Wnt/β-catenin pathway, activities of cellular autophagy and cell apoptosis were examined by assessing the expression levels of several key proteins, including Bcl-2, Bax, and LC3.

RESULTS

We found that NR3C2 was found a significantly lower level in colon cancer tissues than in adjacent tissues, which was associated with distant and lymphatic metastases, clinical stage, and poor clinical outcome, and it was an independent prognostic factor and potential marker of colon cancer. Single-cell transcriptome data identified the subset of circulating T and B cells with high expression of NR3C2, which is involved in TNF signaling pathway. Functional experiments show that downregulation of NR3C2 resultsed in the activation of the Wnt/β-catenin signaling pathway, and promotesd the proliferation and invasion of colon cancer cells while suppressing cell autophagy and apoptosis.

CONCLUSION

NR3C2 may regulate Wnt/β-catenin to affect the proliferation, invasion apoptosis and autophagy of colon cancer, and this axis is a potential target for the treatment of colon cancer.

MIF/NR3C2 axis regulates glucose metabolism reprogramming in pancreatic cancer through MAPK-ERK and AP-1 pathways

Inflammation and aberrant cellular metabolism are widely recognized as hallmarks of cancer. In pancreatic ductal adenocarcinoma (PDAC), inflammatory signaling and metabolic reprogramming are tightly interwoven, playing pivotal roles in the pathogenesis and progression of the disease. However, the regulatory functions of inflammatory mediators in metabolic reprogramming in pancreatic cancer have not been fully explored. Earlier, we demonstrated that pro-inflammatory mediator macrophage migration inhibitory factor (MIF) enhances disease progression by inhibiting its downstream transcriptional factor nuclear receptor subfamily 3 group C member 2 (NR3C2). Here, we provide evidence that MIF and NR3C2 interactively regulate metabolic reprogramming, resulting in MIF-induced cancer growth and progression in PDAC. MIF positively correlates with the HK1 (hexokinase 1), HK2 (hexokinase 2) and LDHA (lactate dehydrogenase) expression and increased pyruvate and lactate production in PDAC patients. Additionally, MIF augments glucose uptake and lactate efflux by upregulating HK1, HK2 and LDHA expression in pancreatic cancer cells in vitro and in mouse models of PDAC. Conversely, a reduction in HK1, HK2 and LDHA expression is observed in tumors with high NR3C2 expression in PDAC patients. NR3C2 suppresses HK1, HK2 and LDHA expression, thereby inhibiting glucose uptake and lactate efflux in pancreatic cancer. Mechanistically, MIF-mediated regulation of glycolytic metabolism involves the activation of the mitogen-activated protein kinase-ERK signaling pathway, whereas NR3C2 interacts with the activator protein 1 to regulate glycolysis. Our findings reveal an interactive role of the MIF/NR3C2 axis in regulating glucose metabolism supporting tumor growth and progression and may be a potential target for designing novel approaches for improving disease outcome.

Delineating the role of nuclear receptors in colorectal cancer, a focused review

Colorectal cancer (CRC) stands as one of the most prevalent form of cancer globally, causing a significant number of deaths, surpassing 0.9 million in the year 2020. According to GLOBOCAN 2020, CRC ranks third in incidence and second in mortality in both males and females. Despite extensive studies over the years, there is still a need to establish novel therapeutic targets to enhance the patients' survival rate in CRC. Nuclear receptors (NRs) are ligand-activated transcription factors (TFs) that regulate numerous essential biological processes such as differentiation, development, physiology, reproduction, and cellular metabolism. Dysregulation and anomalous expression of different NRs has led to multiple alterations, such as impaired signaling cascades, mutations, and epigenetic changes, leading to various diseases, including cancer. It has been observed that differential expression of various NRs might lead to the initiation and progression of CRC, and are correlated with poor survival outcomes in CRC patients. Despite numerous studies on the mechanism and role of NRs in this cancer, it remains of significant scientific interest primarily due to the diverse functions that various NRs exhibit in regulating key hallmarks of this cancer. Thus, modulating the expression of NRs with their agonists and antagonists, based on their expression levels, holds an immense prospect in the diagnosis, prognosis, and therapeutical modalities of CRC. In this review, we primarily focus on the role and mechanism of NRs in the pathogenesis of CRC and emphasized the significance of targeting these NRs using a variety of agents, which may represent a novel and effective strategy for the prevention and treatment of this cancer.

Mineralocorticoid receptor signaling inhibits bladder cancer progression

The biological or clinical significance of mineralocorticoid receptor (MR) in urothelial cancer remains largely unknown. The present study aimed to determine the functional role of MR in bladder cancer progression. In two of the human bladder cancer lines expressing MR, treatment with a natural MR ligand, aldosterone, significantly reduced cell proliferation and migration, which was restored by three MR antagonists clinically used, spironolactone (except colony formation of androgen receptor-positive cells cultured in the presence of androgens), eplerenone, and esaxerenone. Similarly, MR knockdown via shRNA virus infection resulted in significant increases in cell viability/migration, as well as colony formation, compared with control sublines. In addition, MR knockdown augmented the expression of β-catenin, c-fos, and N-cadherin, and lowered that of E-cadherin and p53, indicating the induction of the cadherin switching. Immunohistochemistry in surgical specimens detected MR signals in 58 (92.1%; 36.5% weakly-positive/1+, 44.4% moderately-positive/2+, and 11.1% strongly-positive/3+) of 63 muscle-invasive bladder cancers, which was significantly lower than in adjacent non-neoplastic urothelial tissues (100%; 15.7% 1+, 37.3% 2+, and 47.1% 3+). Moreover, patients with MR-high (3+) tumor had a significantly lower risk of cancer-specific mortality (P=0.039). Multivariable analysis further showed that strong MR expression was an independent predictor of cancer-specific survival in patients with muscle-invasive bladder cancer (hazard ratio 0.117, P=0.039). These findings suggest that MR signaling functions as a tumor suppressor in urothelial carcinoma and prevents tumor growth. Accordingly, there is a possibility that the concurrent use of anti-mineralocorticoids, particularly eplerenone and esaxerenone, in patients with bladder cancer rather contributes to the promotion of disease progression.

Epigenetic disruption of the RARγ complex impairs its function to bookmark AR enhancer interactions required for enzalutamide sensitivity in prostate cancer

The current study in prostate cancer (PCa) focused on the genomic mechanisms at the cross-roads of pro-differentiation signals and the emergence of lineage plasticity. We explored an understudied cistromic mechanism involving RARγ's ability to govern AR cistrome-transcriptome relationships, including those associated with more aggressive PCa features. The RARγ complex in PCa cell models was enriched for canonical cofactors, as well as proteins involved in RNA processing and bookmarking. Identifying the repertoire of miR-96 bound and regulated gene targets, including those recognition elements marked by m6A, revealed their significant enrichment in the RARγ complex. RARγ significantly enhanced the AR cistrome, particularly in active enhancers and super-enhancers, and overlapped with the binding of bookmarking factors. Furthermore, RARγ expression led to nucleosome-free chromatin enriched with H3K27ac, and significantly enhanced the AR cistrome in G2/M cells. RARγ functions also antagonized the transcriptional actions of the lineage master regulator ONECUT2. Similarly, gene programs regulated by either miR-96 or antagonized by RARγ were enriched in alternative lineages and more aggressive PCa phenotypes. Together these findings reveal an under-investigated role for RARγ, modulated by miR-96, to bookmark enhancer sites during mitosis. These sites are required by the AR to promote transcriptional competence, and emphasize luminal differentiation, while antagonizing ONECUT2.

Delineating the glioblastoma stemness by genes involved in cytoskeletal rearrangements and metabolic alterations

Literature data on glioblastoma ongoingly underline the link between metabolism and cancer stemness, the latter is one responsible for potentiating the resistance to treatment, inter alia due to increased invasiveness. In recent years, glioblastoma stemness research has bashfully introduced a key aspect of cytoskeletal rearrangements, whereas the impact of the cytoskeleton on invasiveness is well known. Although non-stem glioblastoma cells are less invasive than glioblastoma stem cells (GSCs), these cells also acquire stemness with greater ease if characterized as invasive cells and not tumor core cells. This suggests that glioblastoma stemness should be further investigated for any phenomena related to the cytoskeleton and metabolism, as they may provide new invasion-related insights. Previously, we proved that interplay between metabolism and cytoskeleton existed in glioblastoma. Despite searching for cytoskeleton-related processes in which the investigated genes might have been involved, not only did we stumble across the relation to metabolism but also reported genes that were found to be implicated in stemness. Thus, dedicated research on these genes in GSCs seems justifiable and might reveal novel directions and/or biomarkers that could be utilized in the future. Herein, we review the previously identified cytoskeleton/metabolism-related genes through the prism of glioblastoma stemness.

Protective role of mineralocorticoid receptor signaling in urothelial tumorigenesis

The expression status of mineralocorticoid receptor (MR) and its biological significance in human urothelial carcinoma remain unknown. The present study aimed to determine the functional role of MR in the development of urothelial cancer. In human normal urothelial SVHUC cells with exposure to a chemical carcinogen 3-methylcholanthrene (MCA), we assessed the effects of a natural MR ligand, aldosterone, and 3 MR antagonists, including spironolactone, eplerenone, and esaxerenone, as well as knockdown of MR via shRNA virus infection, on their neoplastic/malignant transformation. The in vitro system with carcinogen challenge showed that aldosterone and anti-mineralocorticoids significantly prevented and promoted, respectively, the neoplastic transformation of SVHUC cells. Similarly, MR knockdown in SVHUC cells considerably induced MCA-mediated neoplastic transformation, compared with a control subline. In addition, MR knockdown or antagonist treatment resulted in increases in the expression of β-catenin, c-Fos, and N-cadherin, and a decrease in that of E-cadherin. Meanwhile, spironolactone, which is known to possess anti-androgenic activity, rather suppressed the neoplastic transformation of a SVHUC subline stably expressing wild-type androgen receptor, indicating its dominant effect via the androgen receptor pathway. Immunohistochemistry in surgical specimens detected MR signals in 77 (98.7%; 23.1% weak/1+, 42.3% moderate/2+, and 33.3% strong/3+) of 78 non-invasive bladder tumors, which was significantly (P<0.001) lower than in adjacent non-neoplastic urothelial tissues (100%; 20.5% 2+ and 79.5% 3+). Moreover, the risks for disease recurrence after transurethral surgery were marginally lower in female patients with MR-high (2+/3+) tumor (P=0.068) and significantly lower in all patients with MR-high/glucocorticoid receptor-high tumor (P=0.025), compared with respective controls. These findings suggest that MR signaling functions as a suppressor for urothelial tumorigenesis.

Estrogen-Related Receptors Gene Expression and Copy Number Alteration Association With the Clinicopathologic Characteristics of Breast Cancer

Purpose:

It has been suggested that dysregulation of transcription factors expression or activity plays significant roles in breast cancer (BC) severity and poor prognosis. Therefore, our study aims to thoroughly evaluate the estrogen-related receptor isoforms (ESRRs) expression and copy number alteration (CNA) status and their association with clinicopathologic characteristics in BC.

Methods:

A METABRIC dataset consist of 2509 BC patients’ samples was obtained from the cBioPortal public domain. The gene expression, putative CNA, and relevant tumor information of ESRRs were retrieved. ESRRs messenger RNA (mRNA) expression in BC cell lines was obtained from the Cancer Cell Line Encyclopedia (CCLE). Association and correlation analysis of ESRRs expression with BC clinicopathologic characteristics and molecular subtype were performed. Kaplan–Meier survival analysis was conducted to evaluate the prognostic value of ESRRs expression on patient survival.

Results:

ESRRα expression correlated negatively with patients’ age and overall survival, whereas positively correlated with tumor size, the number of positive lymph nodes, and Nottingham prognostic index (NPI). Conversely, ESRRγ expression was positively correlated with patients’ age and negatively correlated with NPI. ESRRα and ESRRγ expression were significantly associated with tumor grade, expression of hormone receptors, human epidermal growth factor receptor 2 (HER2), and molecular subtype, whereas ESRRβ was only associated with tumor stage. A significant and distinct association of each of ESRRs CNA with various clinicopathologic and prognostic factors was also observed. Kaplan–Meier survival analysis demonstrated no significant difference for survival curves among BC patients with high or low expression of ESRRα, β, or γ. On stratification, high ESRRα expression significantly reduced survival among premenopausal patients, patients with grade I/II, and early-stage disease. In BC cell lines, only ESRRα expression was significantly higher in HER2-positive cells. No significant association was observed between ESRRβ expression and any of the clinicopathologic characteristics examined.

Conclusions:

In this clinical dataset, ESRRα and ESRRγ mRNA expression and CNA show a significant correlation and association with distinct clinicopathologic and prognostic parameters known to influence treatment outcomes; however, ESRRβ failed to show a robust role in BC pathogenesis. ESRRα and ESRRγ can be employed as therapeutic targets in BC-targeted therapy. However, the role of ESRRβ in BC pathogenesis remains unclear.

Molecular Signatures of Tumour and Its Microenvironment for Precise Quantitative Diagnosis of Oral Squamous Cell Carcinoma: An International Multi-Cohort Diagnostic Validation Study

BACKGROUND

Heterogeneity in oral potentially malignant disorder (OPMD) poses a problem for accurate prognosis that impacts on treatment strategy and patient outcome. A holistic assessment based on gene expression signatures from both the tumour cells and their microenvironment is necessary to provide a more precise prognostic assessment than just tumour cell signatures alone.

METHODS

We reformulated our previously established multigene qPCR test, quantitative Malignancy Index Diagnostic System (qMIDS) with new genes involved in matrix/stroma and immune modulation of the tumour microenvironment. An algorithm calculates and converts a panel of 16 gene mRNA expression levels into a qMIDS index to quantify risk of malignancy for each sample.

RESULTS

The new qMIDS^V2 assay was validated in a UK oral squamous cell carcinoma (OSCC) cohort (n = 282) of margin and tumour core samples demonstrating significantly better diagnostic performance (AUC = 0.945) compared to previous qMIDS^V1 (AUC = 0.759). Performance of qMIDS^V2 were independently validated in Chinese (n = 35; AUC = 0.928) and Indian (n = 95; AUC = 0.932) OSCC cohorts. Further, 5-year retrospective analysis on an Indian dysplastic lesion cohort (n = 30) showed that qMIDS^V2 was able to significantly differentiate between lesions without transformation and those with malignant transformation.

CONCLUSIONS

This study validated a novel multi-gene qPCR test on a total of 535 tissue specimens from UK, China and India, demonstrating a rapid minimally invasive method that has a potential application for dysplasia risk stratification. Further study is required to establish if qMIDS^V2 could be used to improve OPMD patient management, guide treatment strategy and reduce oral cancer burden.

A Comprehensive Bioinformatic Analysis of NOTCH Pathway Involvement in Stomach Adenocarcinoma

Background

Activation of NOTCH signaling pathways, which are key regulators of multiple cellular functions, has been frequently implicated in cancer pathogenesis, and NOTCH inhibitors have received much recent focus in the context of cancer therapeutics. However, the role and possible involvement of NOTCH pathways in stomach adenocarcinoma (STAD) are unclear. Here, putative regulatory mechanisms and functions of NOTCH pathways in STAD were investigated.

Methods

Publicly available data from the TCGA-STAD database were utilized to explore the involvement of canonical NOTCH pathways in STAD by analyzing RNA expression levels of NOTCH receptors, ligands, and downstream genes. Statistical analysis of the data pertaining to cancer and noncancerous samples was performed using R software packages and public databases/webservers.

Results

Significant differential gene expression between control and STAD samples was noted for all NOTCH receptors (NOTCH1, 2, 3, and 4), the delta-like NOTCH ligands (DLL-3 and 4), and typical downstream genes (HES1 and HEY1). Four genes (NOTCH1, NOTCH2, NOTCH3, and HEY1) presented prognostic values for the STAD outcome in terms of overall survival. Functional enrichment analysis indicated that NOTCH family genes-strongly correlated genes were mainly enriched in several KEGG signaling pathways such as the PI3K-Akt signaling pathway, human papillomavirus infection, focal adhesion, Rap1 signaling pathway, and ECM-receptor interaction. Gene set enrichment analysis (GSEA) results showed that NOTCH family genes-significantly correlated genes were mainly enriched in four signaling pathways, ECM (extracellular matrix), tumor angiogenesis, inflammatory response, and immune regulation.

Conclusions

NOTCH family genes may play an essential role in the progression of STAD by modulating immune cells and mediating ECM synthesis, angiogenesis, focal adhesion, and PI3K-Akt signaling. Multiple NOTCH family genes are valuable candidate biomarkers or therapeutic targets for the management of STAD.

Systematic pan-cancer characterization of nuclear receptors identifies potential cancer biomarkers and therapeutic targets

The nuclear receptor (NR) superfamily is one of the major druggable gene families, representing targets of approximately 13.5% of approved drugs. Certain NRs, such as estrogen receptor and androgen receptor, have been well demonstrated to be functionally involved in cancer and serve as informative biomarkers and therapeutic targets in oncology. However, the spectrum of NR dysregulation across cancers remains to be comprehensively characterized. Through computational integration of genetic, genomic, and pharmacologic profiles, we characterized the expression, recurrent genomic alterations, and cancer dependency of NRs at a large scale across primary tumor specimens and cancer cell lines. Expression levels of NRs were highly cancer-type specific and globally downregulated in tumors compared to corresponding normal tissue. Although the majority of NRs showed copy number losses in cancer, both recurrent focal gains and losses were identified in select NRs. Recurrent mutations and transcript fusions of NRs were observed in a small portion of cancers, serving as actionable genomic alterations. Analysis of large-scale CRISPR and RNAi screening datasets identified 10 NRs as strongly selective essential genes for cancer cell growth. In a subpopulation of tumor cells, growth dependencies correlated significantly with expression or genomic alterations. Overall, our comprehensive characterization of NRs across cancers may facilitate the identification and prioritization of potential biomarkers and therapeutic targets, as well as the selection of patients for precision cancer treatment.

A network pharmacology perspective for deciphering potential mechanisms of action of Solanum nigrum L. in bladder cancer

BACKGROUND

Solanum nigrum L. decoction has been used as a folklore medicine in China to prevent the postoperative recurrence of bladder cancer (BC). However, there are no previous pharmacological studies on the protective mechanisms of this activity of the plant. Thus, this study aimed to perform a systematic analysis and to predict the potential action mechanisms underlying S. nigrum activity in BC based on network pharmacology.

METHODS

Based on network pharmacology, the active ingredients of S. nigrum and the corresponding targets were identified using the Traditional Chinese Medicines for Systems Pharmacology Database and Analysis Platform database, and BC-related genes were screened using GeneCards and the Online Mendelian Inheritance in Man database. In addition, ingredient-target (I-T) and protein-protein interaction (PPI) networks were constructed using STRING and Cytoscape, Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were conducted, and then the pathways directly related to BC were integrated manually to reveal the pharmacological mechanism underlying S. nigrum-medicated therapeutic effects in BC.

RESULTS

Seven active herbal ingredients from 39 components of S. nigrum were identified, which shared 77 common target genes related to BC. I-T network analysis revealed that quercetin was associated with all targets and that NCOA2 was targeted by four ingredients. Besides, interleukin 6 had the highest degree value in the PPI network, indicating a hub role. A subsequent gene enrichment analysis yielded 86 significant GO terms and 89 significant pathways, implying that S. nigrum had therapeutic benefits in BC through multi-pathway effects, including the HIF-1, TNF, P53, MAPK, PI3K/Akt, apoptosis and bladder cancer pathway.

CONCLUSIONS

S. nigrum may mediate pharmacological effects in BC through multi-target and various signaling pathways. Further validation is required experimentally. Network pharmacology approach provides a predicative novel strategy to reveal the holistic mechanism of action of herbs.

Nuclear receptors: Lipid and hormone sensors with essential roles in the control of cancer development

Nuclear receptors (NRs) are a superfamily of ligand-activated transcription factors that act as biological sensors and use a combination of mechanisms to modulate positively and negatively gene expression in a spatial and temporal manner. The highly orchestrated biological actions of several NRs influence the proliferation, differentiation, and apoptosis of many different cell types. Synthetic ligands for several NRs have been the focus of extensive drug discovery efforts for cancer intervention. This review summarizes the roles in tumour growth and metastasis of several relevant NR family members, namely androgen receptor (AR), estrogen receptor (ER), glucocorticoid receptor (GR), thyroid hormone receptor (TR), retinoic acid receptors (RARs), retinoid X receptors (RXRs), peroxisome proliferator-activated receptors (PPARs), and liver X receptors (LXRs). These studies are key to develop improved therapeutic agents based on novel modes of action with reduced side effects and overcoming resistance.

Performance-weighted-voting model: An ensemble machine learning method for cancer type classification using whole-exome sequencing mutation

Background

With improvements in next-generation DNA sequencing technology, lower cost is needed to collect genetic data. More machine learning techniques can be used to help with cancer analysis and diagnosis.

Methods

We developed an ensemble machine learning system named performance-weighted-voting model for cancer type classification in 6,249 samples across 14 cancer types. Our ensemble system consists of five weak classifiers (logistic regression, SVM, random forest, XGBoost and neural networks). We first used cross-validation to get the predicted results for the five classifiers. The weights of the five weak classifiers can be obtained based on their predictive performance by solving linear regression functions. The final predicted probability of the performance-weighted-voting model for a cancer type can be determined by the summation of each classifier's weight multiplied by its predicted probability.

Results

Using the somatic mutation count of each gene as the input feature, the overall accuracy of the performance-weighted-voting model reached 71.46%, which was significantly higher than the five weak classifiers and two other ensemble models: the hard-voting model and the soft-voting model. In addition, by analyzing the predictive pattern of the performance-weighted-voting model, we found that in most cancer types, higher tumor mutational burden can improve overall accuracy.

Conclusion

This study has important clinical significance for identifying the origin of cancer, especially for those where the primary cannot be determined. In addition, our model presents a good strategy for using ensemble systems for cancer type classification.

Identification of transcription factor co-regulators that drive prostate cancer progression

In prostate cancer (PCa), and many other hormone-dependent cancers, there is clear evidence for distorted transcriptional control as disease driver mechanisms. Defining which transcription factor (TF) and coregulators are altered and combine to become oncogenic drivers remains a challenge, in part because of the multitude of TFs and coregulators and the diverse genomic space on which they function. The current study was undertaken to identify which TFs and coregulators are commonly altered in PCa. We generated unique lists of TFs (n = 2662), coactivators (COA; n = 766); corepressors (COR; n = 599); mixed function coregulators (MIXED; n = 511), and to address the challenge of defining how these genes are altered we tested how expression, copy number alterations and mutation status varied across seven prostate cancer (PCa) cohorts (three of localized and four advanced disease). Testing of significant changes was undertaken by bootstrapping approaches and the most significant changes were identified. For one commonly and significantly altered gene were stably knocked-down expression and undertook cell biology experiments and RNA-Seq to identify differentially altered gene networks and their association with PCa progression risks. COAS, CORS, MIXED and TFs all displayed significant down-regulated expression (q.value < 0.1) and correlated with protein expression (r 0.4-0.55). In localized PCa, stringent expression filtering identified commonly altered TFs and coregulator genes, including well-established (e.g. ERG) and underexplored (e.g. PPARGC1A, encodes PGC1α). Reduced PPARGC1A expression significantly associated with worse disease-free survival in two cohorts of localized PCa. Stable PGC1α knockdown in LNCaP cells increased growth rates and invasiveness and RNA-Seq revealed a profound basal impact on gene expression (~ 2300 genes; FDR < 0.05, logFC > 1.5), but only modestly impacted PPARγ responses. GSEA analyses of the PGC1α transcriptome revealed that it significantly altered the AR-dependent transcriptome, and was enriched for epigenetic modifiers. PGC1α-dependent genes were overlapped with PGC1α-ChIP-Seq genes and significantly associated in TCGA with higher grade tumors and worse disease-free survival. These methods and data demonstrate an approach to identify cancer-driver coregulators in cancer, and that PGC1α expression is clinically significant yet underexplored coregulator in aggressive early stage PCa.

The Modulation of Phase II Drug-Metabolizing Enzymes in Proliferating and Differentiated CaCo-2 Cells by Hop-Derived Prenylflavonoids

Prenylflavonoids in the human organism exhibit various health-beneficial activities, although they may interfere with drugs via the modulation of the expression and/or activity of drug-metabolizing enzymes. As intestinal cells are exposed to the highest concentrations of prenylflavonoids, we decided to study the cytotoxicity and modulatory effects of the four main hop-derived prenylflavonoids on the activities and mRNA expression of the main drug-conjugating enzymes in human CaCo-2 cells. Proliferating CaCo-2 cells were used for these purposes as a model of colorectal cancer cells, and differentiated CaCo-2 cells were used as an enterocyte-like model. All the tested prenylflavonoids inhibited the CaCo-2 cells proliferation, with xanthohumol proving the most effective (IC₅₀ 8.5 µM). The prenylflavonoids modulated the activities and expressions of the studied enzymes to a greater extent in the differentiated, as opposed to the proliferating, CaCo-2 cells. In the differentiated cells, all the prenylflavonoids caused a marked increase in glutathione S-transferase and catechol-O-methyltransferase activities, while the activity of sulfotransferase was significantly inhibited. Moreover, the prenylflavonoids upregulated the mRNA expression of uridine diphosphate (UDP)-glucuronosyl transferase 1A6 and downregulated that of glutathione S-transferase 1A1/2.

LINC01128 resisted acute myeloid leukemia through regulating miR-4260/NR3C2

Acute myeloid leukemia (AML) is a prevalent class of blood disease with a high occurrence rate and relapse rate. The role of dysregulated microRNAs (miRNAs) in AML is emerging. MiR-4260 was identified to be a carcinogenic miRNA in colorectal cancer, but never has it been reported in AML. We aimed to study the function and mechanism of miR-4260 in AML. The miR-4260 level was higher in AML cell lines than the normal cell lines. Inhibition of miR-4260 hindered proliferation and increased apoptosis of AML cells. Mechanistically, long intergenic non-protein coding RNA 1128 (LINC01128) competed with nuclear receptor subfamily 3 group C member 2 (NR3C2) for miR-4260 so as to upregulate NR3C2. We identified the reduced levels of LINC01128 and NR3C2 in AML and it was suggested through rescue assays that LINC01128 repressed AML progression through regulating miR-4260/NR3C2 axis. In conclusion, we firstly uncovered that LINC01128 resisted acute myeloid leukemia through regulating miR-4260/NR3C2, providing novel clues for the treatment improvement of AML.

Pan-cancer analyses of human nuclear receptors reveal transcriptome diversity and prognostic value across cancer types

The human nuclear receptor (NR) superfamily comprises 48 ligand-dependent transcription factors that play regulatory roles in physiology and pathophysiology. In cancer, NRs have long served as predictors of disease stratification, treatment response, and clinical outcome. The Cancer Genome Atlas (TCGA) Pan-Cancer project provides a wealth of genetic data for a large number of human cancer types. Here, we examined NR transcriptional activity in 8,526 patient samples from 33 TCGA ‘Pan-Cancer’ diseases and 11 ‘Pan-Cancer’ organ systems using RNA sequencing data. The web-based Kaplan-Meier (KM) plotter tool was then used to evaluate the prognostic potential of NR gene expression in 21/33 cancer types. Although, most NRs were significantly underexpressed in cancer, NR expression (moderate to high expression levels) was predominantly restricted (46%) to specific tissues, particularly cancers representing gynecologic, urologic, and gastrointestinal ‘Pan-Cancer’ organ systems. Intriguingly, a relationship emerged between recurrent positive pairwise correlation of Class IV NRs in most cancers. NR expression was also revealed to play a profound effect on patient overall survival rates, with ≥5 prognostic NRs identified per cancer type. Taken together, these findings highlighted the complexity of NR transcriptional networks in cancer and identified novel therapeutic targets for specific cancer types.

The orphan nuclear receptor estrogen-related receptor beta (ERRβ) in triple-negative breast cancer

PURPOSE

Triple-negative breast cancer (TNBC)/basal-like breast cancer (BLBC) is a highly aggressive form of breast cancer. We previously reported that a small molecule agonist ligand for the orphan nuclear receptor estrogen-related receptor beta (ERRβ or ESRRB) has growth inhibitory and anti-mitotic activity in TNBC cell lines. In this study, we evaluate the association of ESRRB mRNA, copy number levels, and protein expression with demographic, clinicopathological, and gene expression features in breast tumor clinical specimens.

METHODS

ESRRB mRNA-level expression and clinical associations were analyzed using RNAseq data. Array-based comparative genomic hybridization determined ESRRB copy number in African-American and Caucasian women. Transcription factor activity was measured using promoter-reporter luciferase assays in TNBC cell lines. Semi-automatic quantification of immunohistochemistry measured ERRβ protein expression on a 150-patient tissue microarray series.

RESULTS

ESRRB mRNA expression is significantly lower in TNBC/BLBC versus other breast cancer subtypes. There is no evidence of ESRRB copy number loss. ESRRB mRNA expression is correlated with the expression of genes associated with neuroactive ligand-receptor interaction, metabolic pathways, and deafness. These genes contain G/C-rich transcription factor binding motifs. The ESRRB message is alternatively spliced into three isoforms, which we show have different transcription factor activity in basal-like versus other TNBC cell lines. We further show that the ERRβ2 and ERRβsf isoforms are broadly expressed in breast tumors at the protein level.

CONCLUSIONS

Decreased ESRRB mRNA expression and distinct patterns of ERRβ isoform subcellular localization and transcription factor activity are key features in TNBC/BLBC.

Vitamin D-VDR Signaling Inhibits Wnt/β-Catenin-Mediated Melanoma Progression and Promotes Antitumor Immunity

1α,25-Dihydroxyvitamin D3 signals via the vitamin D receptor (VDR). Higher serum vitamin D is associated with thinner primary melanoma and better outcome, although a causal mechanism has not been established. As patients with melanoma commonly avoid sun exposure, and consequent vitamin D deficiency might worsen outcomes, we interrogated 703 primary melanoma transcriptomes to understand the role of vitamin D-VDR signaling and replicated the findings in The Cancer Genome Atlas metastases. VDR expression was independently protective for melanoma-related death in both primary and metastatic disease. High tumor VDR expression was associated with upregulation of pathways mediating antitumor immunity and corresponding with higher imputed immune cell scores and histologically detected tumor-infiltrating lymphocytes. High VDR-expressing tumors had downregulation of proliferative pathways, notably Wnt/β-catenin signaling. Deleterious low VDR levels resulted from promoter methylation and gene deletion in metastases. Vitamin D deficiency (<25 nmol/L ∼ 10 ng/mL) shortened survival in primary melanoma in a VDR-dependent manner. In vitro functional validation studies showed that elevated vitamin D-VDR signaling inhibited Wnt/β-catenin signaling genes. Murine melanoma cells overexpressing VDR produced fewer pulmonary metastases than controls in tail-vein metastasis assays. In summary, vitamin D-VDR signaling contributes to controlling pro-proliferative/immunosuppressive Wnt/β-catenin signaling in melanoma and this is associated with less metastatic disease and stronger host immune responses. This is evidence of a causal relationship between vitamin D-VDR signaling and melanoma survival, which should be explored as a therapeutic target in primary resistance to checkpoint blockade. SIGNIFICANCE: VDR expression could potentially be used as a biomarker to stratify patients with melanoma that may respond better to immunotherapy.

Retinoid X receptor gamma (RXRG) is an independent prognostic biomarker in ER-positive invasive breast cancer

BACKGROUND

Retinoid X Receptor Gamma (RXRG) is a member of the nuclear receptor superfamily and plays a role in tumour suppression. This study aims to explore the prognostic significance of RXRG in breast cancer.

METHODS

Primary breast cancer tissue microarrays (n = 923) were immuno-stained for RXRG protein and correlated with clinicopathological features, and patient outcome.

RESULTS

Nuclear RXRG expression was significantly associated with smaller tumour size (p = 0.036), lower grade (p < 0.001), lobular histology (p = 0.016), lower Nottingham Prognostic Index (p = 0.04) and longer breast cancer-specific survival (p < 0.001), and longer time to distant metastasis (p = 0.002). RXRG expression showed positive association with oestrogen receptor (ER)-related biomarkers: GATA3, FOXA1, STAT3 and MED7 (all p < 0.001) and a negative correlation with the Ki67 proliferation marker. Multivariate analysis demonstrated RXRG protein as an independent predictor of longer breast cancer-specific survival and distant metastasis-free survival. In the external validation cohorts, RXRG expression was associated with improved patients' outcome (p = 0.025). In ER-positive tumours, high expression of RXRG was associated with better patient outcome regardless of adjuvant systemic therapy. ER signalling pathway was the top predicted master regulator of RXRG protein expression (p = 0.005).

CONCLUSION

This study provides evidence for the prognostic value of RXRG in breast cancer particularly the ER-positive tumours.

Role of glucocorticoid signaling in urothelial tumorigenesis: Inhibition by prednisone presumably through inducing glucocorticoid receptor transrepression

Glucocorticoids, including dexamethasone (DEX) and prednisone (PRED), have been prescribed in patients with neoplastic disease as cytotoxic agents or comedications. Nonetheless, it remains uncertain whether they have an impact on the development of bladder cancer. We, therefore, assessed the functional role of the glucocorticoid-mediated glucocorticoid receptor (GR) signaling in urothelial tumorigenesis. Tumor formation was significantly delayed in xenograft-bearing mice with implantation of control bladder cancer UMUC3 cells or nonneoplastic urothelial SVHUC cells undergoing malignant transformation induced by a chemical carcinogen 3-methylcholanthrene (MCA), compared with respective GR knockdown xenografts. Using the in vitro system with MCA-SVHUC cells, we screened 11 GR ligands, including DEX, and found significant inhibitory effects of PRED on their neoplastic transformation. The effects of PRED were restored by a GR antagonist RU486 in GR-positive MCA-SVHUC cells, while PRED failed to inhibit the neoplastic transformation of GR knockdown cells. Significant decreases in the expression levels of oncogenes (c-Fos/c-Jun) and significant increases in those of a tumor suppressor UGT1A were seen in MCA-SVHUC-control cells (vs GR-short hairpin RNA) or PRED-treated MCA-SVHUC-control cells (vs mock). In addition, N-butyl-N-(4-hydroxybutyl) nitrosamine induced bladder cancer in all of eight mock-treated mice vs seven (87.5%) of DEX-treated (P = .302) or four (50%) of PRED-treated (P = .021) animals. Finally, DEX was found to considerably induce both transactivation (activation of glucocorticoid-response element mediated transcription and expression of its targets) and transrepression (suppression of nuclear factor-kappa B transactivation and expression of its regulated genes) of GR in SVHUC cells, while PRED more selectively induced GR transrepression. These findings suggest that PRED could prevent urothelial tumorigenesis presumably via inducing GR transrepression.

Effect of the Large and Small T-Antigens of Human Polyomaviruses on Signaling Pathways

Viruses are intracellular parasites that require a permissive host cell to express the viral genome and to produce new progeny virus particles. However, not all viral infections are productive and some viruses can induce carcinogenesis. Irrespective of the type of infection (productive or neoplastic), viruses hijack the host cell machinery to permit optimal viral replication or to transform the infected cell into a tumor cell. One mechanism viruses employ to reprogram the host cell is through interference with signaling pathways. Polyomaviruses are naked, double-stranded DNA viruses whose genome encodes the regulatory proteins large T-antigen and small t-antigen, and structural proteins that form the capsid. The large T-antigens and small t-antigens can interfere with several host signaling pathways. In this case, we review the interplay between the large T-antigens and small t-antigens with host signaling pathways and the biological consequences of these interactions.

Exploiting vulnerabilities of cancer by targeting nuclear receptors of stromal cells in tumor microenvironment

The tumor microenvironment is a complex and dynamic cellular community comprising the tumor epithelium and various tumor-supporting cells such as immune cells, fibroblasts, immunosuppressive cells, adipose cells, endothelial cells, and pericytes. The interplay between the tumor microenvironment and tumor cells represents a key contributor to immune evasiveness, physiological hardiness and the local and systemic invasiveness of malignant cells. Nuclear receptors are master regulators of physiological processes and are known to play pro-/anti-oncogenic activities in tumor cells. However, the actions of nuclear receptors in tumor-supporting cells have not been widely studied. Given the excellent druggability and extensive regulatory effects of nuclear receptors, understanding their biological functionality in the tumor microenvironment is of utmost importance. Therefore, the present review aims to summarize recent evidence about the roles of nuclear receptors in tumor-supporting cells and their implications for malignant processes such as tumor proliferation, evasion of immune surveillance, angiogenesis, chemotherapeutic resistance, and metastasis. Based on findings derived mostly from cell culture studies and a few in vivo animal cancer models, the functions of VDR, PPARs, AR, ER and GR in tumor-supporting cells are relatively well-characterized. Evidence for other receptors, such as RARβ, RORγ, and FXR, is limited yet promising. Hence, the nuclear receptor signature in the tumor microenvironment may harbor prognostic value. The clinical prospects of a tumor microenvironment-oriented cancer therapy exploiting the nuclear receptors in different tumor-supporting cells are also encouraging. The major challenge, however, lies in the ability to develop a highly specific drug delivery system to facilitate precision medicine in cancer therapy.

The miR-96 and RARγ signaling axis governs androgen signaling and prostate cancer progression

Expression levels of retinoic acid receptor gamma (NR1B3/RARG, encodes RARγ) are commonly reduced in prostate cancer (PCa). Therefore, we sought to establish the cellular and gene regulatory consequences of reduced RARγ expression, and determine RARγ regulatory mechanisms. RARG shRNA approaches in non-malignant (RWPE-1 and HPr1-AR) and malignant (LNCaP) prostate models revealed that reducing RARγ levels, rather than adding exogenous retinoid ligand, had the greatest impact on prostate cell viability and gene expression. ChIP-Seq defined the RARγ cistrome, which was significantly enriched at active enhancers associated with AR binding sites. Reflecting a significant genomic role for RARγ to regulate androgen signaling, RARγ knockdown in HPr1-AR cells significantly regulated the magnitude of the AR transcriptome. RARγ downregulation was explained by increased miR-96 in PCa cell and mouse models, and TCGA PCa cohorts. Biochemical approaches confirmed that miR-96 directly regulated RARγ expression and function. Capture of the miR-96 targetome by biotin-miR-96 identified that RARγ and a number of RARγ interacting co-factors including TACC1 were all targeted by miR-96, and expression of these genes were prominently altered, positively and negatively, in the TCGA-PRAD cohort. Differential gene expression analyses between tumors in the TCGA-PRAD cohort with lower quartile expression levels of RARG and TACC1 and upper quartile miR-96, compared to the reverse, identified a gene network including several RARγ target genes (e.g., SOX15) that significantly associated with worse disease-free survival (hazard ratio 2.23, 95% CI 1.58 to 2.88, p = 0.015). In summary, miR-96 targets a RARγ network to govern AR signaling, PCa progression and disease outcome.

Pan-cancer analysis of transcriptional metabolic dysregulation using The Cancer Genome Atlas

Understanding metabolic dysregulation in different disease settings is vital for the safe and effective incorporation of metabolism-targeted therapeutics in the clinic. Here, using transcriptomic data for 10,704 tumor and normal samples from The Cancer Genome Atlas, across 26 disease sites, we present a novel bioinformatics pipeline that distinguishes tumor from normal tissues, based on differential gene expression for 114 metabolic pathways. We confirm pathway dysregulation in separate patient populations, demonstrating the robustness of our approach. Bootstrapping simulations were then applied to assess the biological significance of these alterations. We provide distinct examples of the types of analysis that can be accomplished with this tool to understand cancer specific metabolic dysregulation, highlighting novel pathways of interest, and patterns of metabolic flux, in both common and rare disease sites. Further, we show that Master Metabolic Transcriptional Regulators explain why metabolic differences exist, can segregate patient populations, and predict responders to different metabolism-targeted therapeutics.

Vitamin D Receptor Signaling and Cancer

The vitamin D receptor (VDR) binds the secosteroid hormone 1,25(OH)₂D₃ with high affinity and regulates gene programs that control a serum calcium levels, as well as cell proliferation and differentiation. A significant focus has been to exploit the VDR in cancer settings. Although preclinical studies have been strongly encouraging, to date clinical trials have delivered equivocal findings that have paused the clinical translation of these compounds. However, it is entirely possible that mining of genomic data will help to refine precisely what are the key anticancer actions of vitamin D compounds and where these can be used most effectively.

Integrative genomic approaches to dissect clinically-significant relationships between the VDR cistrome and gene expression in primary colon cancer

Recently, we undertook a pan-cancer analyses of the nuclear hormone receptor (NR) superfamily in The Cancer Genome Atlas (TCGA), and revealed that the vitamin D receptor (NR1I1/VDR) was commonly and significantly down-regulated specifically in colon adenocarcinoma cohort (COAD). To examine the consequence of down-regulated VDR expression we re-analyzed VDR chromatin immunoprecipitation sequencing (ChIP-Seq) data from LS180 colon cancer cells (GSE31939). This analysis identified 1809 loci that displayed significant (p.adj<0.01) differential binding of the VDR in response 1,25(OH)₂D₃ treatment; 947 peaks annotated to 672 genes. We examined expression patterns in the COAD cohort of 286 tumors compared to 41 normal samples and revealed that VDR bound genes were significantly positively correlated to VDR expression compared to the background transcriptome, suggesting direct regulation by VDR. Gene set enrichment analyses revealed significant enrichment for genes known to be regulated by a number of other transcription factors including SMADs and JUN. Filtering VDR associated genes for those that were commonly and significantly altered in COAD revealed a cohort of 27 differentially expressed genes. The expression patterns of these genes clustered tumors and significantly associated with disease free survival. For instance, males with low expression of Lectin, Galactoside Binding Soluble 4 (LGALS4, encodes the colon tumor suppressor, Galactin 4) had significantly shorted disease free survival. These analyses suggest that reduced expression of VDR in colon cancer (but neither loss nor mutation) changes the actions of the VDR by both dampening the expression of tumor suppressors (e.g. LGALS4) whilst either stabilizing or not down-regulating expression of oncogenes (e.g. Carbonic Anhydrase 9 (CA9)). These integrative genomic approaches are relatively generic and applicable to the study of any transcription factor.

Bioinformatic approaches to interrogating vitamin D receptor signaling

Bioinformatics applies unbiased approaches to develop statistically-robust insight into health and disease. At the global, or "20,000 foot" view bioinformatic analyses of vitamin D receptor (NR1I1/VDR) signaling can measure where the VDR gene or protein exerts a genome-wide significant impact on biology; VDR is significantly implicated in bone biology and immune systems, but not in cancer. With a more VDR-centric, or "2000 foot" view, bioinformatic approaches can interrogate events downstream of VDR activity. Integrative approaches can combine VDR ChIP-Seq in cell systems where significant volumes of publically available data are available. For example, VDR ChIP-Seq studies can be combined with genome-wide association studies to reveal significant associations to immune phenotypes. Similarly, VDR ChIP-Seq can be combined with data from Cancer Genome Atlas (TCGA) to infer the impact of VDR target genes in cancer progression. Therefore, bioinformatic approaches can reveal what aspects of VDR downstream networks are significantly related to disease or phenotype.

Gene expression-based biomarkers for discriminating early and late stage of clear cell renal cancer

In this study, an attempt has been made to identify expression-based gene biomarkers that can discriminate early and late stage of clear cell renal cell carcinoma (ccRCC) patients. We have analyzed the gene expression of 523 samples to identify genes that are differentially expressed in the early and late stage of ccRCC. First, a threshold-based method has been developed, which attained a maximum accuracy of 71.12% with ROC 0.67 using single gene NR3C2. To improve the performance of threshold-based method, we combined two or more genes and achieved maximum accuracy of 70.19% with ROC of 0.74 using eight genes on the validation dataset. These eight genes include four underexpressed (NR3C2, ENAM, DNASE1L3, FRMPD2) and four overexpressed (PLEKHA9, MAP6D1, SMPD4, C11orf73) genes in the late stage of ccRCC. Second, models were developed using state-of-art techniques and achieved maximum accuracy of 72.64% and 0.81 ROC using 64 genes on validation dataset. Similar accuracy was obtained on 38 genes selected from subset of genes, involved in cancer hallmark biological processes. Our analysis further implied a need to develop gender-specific models for stage classification. A web server, CancerCSP, has been developed to predict stage of ccRCC using gene expression data derived from RNAseq experiments.

The Genomic Impact of DNA CpG Methylation on Gene Expression; Relationships in Prostate Cancer

The process of DNA CpG methylation has been extensively investigated for over 50 years and revealed associations between changing methylation status of CpG islands and gene expression. As a result, DNA CpG methylation is implicated in the control of gene expression in developmental and homeostasis processes, as well as being a cancer-driver mechanism. The development of genome-wide technologies and sophisticated statistical analytical approaches has ushered in an era of widespread analyses, for example in the cancer arena, of the relationships between altered DNA CpG methylation, gene expression, and tumor status. The remarkable increase in the volume of such genomic data, for example, through investigators from the Cancer Genome Atlas (TCGA), has allowed dissection of the relationships between DNA CpG methylation density and distribution, gene expression, and tumor outcome. In this manner, it is now possible to test that the genome-wide correlations are measurable between changes in DNA CpG methylation and gene expression. Perhaps surprisingly is that these associations can only be detected for hundreds, but not thousands, of genes, and the direction of the correlations are both positive and negative. This, perhaps, suggests that CpG methylation events in cancer systems can act as disease drivers but the effects are possibly more restricted than suspected. Additionally, the positive and negative correlations suggest direct and indirect events and an incomplete understanding. Within the prostate cancer TCGA cohort, we examined the relationships between expression of genes that control DNA methylation, known targets of DNA methylation and tumor status. This revealed that genes that control the synthesis of S-adenosyl-l-methionine (SAM) associate with altered expression of DNA methylation targets in a subset of aggressive tumors.