Integrative genomic analysis in K562 chronic myelogenous leukemia cells reveals that proximal NCOR1 binding positively regulates genes that govern erythroid differentiation and Imatinib sensitivity

GPS2 promotes erythroid differentiation in K562 erythroleukemia cells primarily via NCOR1

G protein pathway suppressor 2 (GPS2) has been shown to play a pivotal role in human and mouse definitive erythropoiesis in an EKLF-dependent manner. However, whether GPS2 affects human primitive erythropoiesis is still unknown. This study demonstrated that GPS2 positively regulates erythroid differentiation in K562 cells, which have a primitive erythroid phenotype. Overexpression of GPS2 promoted hemin-induced hemoglobin synthesis in K562 cells as assessed by the increased percentage of benzidine-positive cells and the deeper red coloration of the cell pellets. In contrast, knockdown of GPS2 inhibited hemin-induced erythroid differentiation of K562 cells. GPS2 overexpression also enhanced erythroid differentiation of K562 cells induced by cytosine arabinoside (Ara-C). GPS2 induced hemoglobin synthesis by increasing the expression of globin and ALAS2 genes, either under steady state or upon hemin treatment. Promotion of erythroid differentiation of K562 cells by GPS2 mainly relies on NCOR1, as knockdown of NCOR1 or lack of the NCOR1-binding domain of GPS2 potently diminished the promotive effect. Thus, our study revealed a previously unknown role of GPS2 in regulating human primitive erythropoiesis in K562 cells.

Vitamin D and Systems Biology

The biological actions of the vitamin D receptor (VDR) have been investigated intensively for over 100 years and has led to the identification of significant insights into the repertoire of its biological actions. These were initially established to be centered on the regulation of calcium transport in the colon and deposition in bone. Beyond these well-known calcemic roles, other roles have emerged in the regulation of cell differentiation processes and have an impact on metabolism. The purpose of the current review is to consider where applying systems biology (SB) approaches may begin to generate a more precise understanding of where the VDR is, and is not, biologically impactful. Two SB approaches have been developed and begun to reveal insight into VDR biological functions. In a top-down SB approach genome-wide scale data are statistically analyzed, and from which a role for the VDR emerges in terms of being a hub in a biological network. Such approaches have confirmed significant roles, for example, in myeloid differentiation and the control of inflammation and innate immunity. In a bottom-up SB approach, current biological understanding is built into a kinetic model which is then applied to existing biological data to explain the function and identify unknown behavior. To date, this has not been applied to the VDR, but has to the related ERα and identified previously unknown mechanisms of control. One arena where applying top-down and bottom-up SB approaches may be informative is in the setting of prostate cancer health disparities.

Nuclear Receptor Coregulators in Hormone-Dependent Cancers

Nuclear receptors (NRs) function collectively as a transcriptional signaling network that mediates gene regulatory actions to either maintain cellular homeostasis in response to hormonal, dietary and other environmental factors, or act as orphan receptors with no known ligand. NR complexes are large and interact with multiple protein partners, collectively termed coregulators. Coregulators are essential for regulating NR activity and can dictate whether a target gene is activated or repressed by a variety of mechanisms including the regulation of chromatin accessibility. Altered expression of coregulators contributes to a variety of hormone-dependent cancers including breast and prostate cancers. Therefore, understanding the mechanisms by which coregulators interact with and modulate the activity of NRs provides opportunities to develop better prognostic and diagnostic approaches, as well as novel therapeutic targets. This review aims to gather and summarize recent studies, techniques and bioinformatics methods used to identify distorted NR coregulator interactions that contribute as cancer drivers in hormone-dependent cancers.

Reduced NCOR2 expression accelerates androgen deprivation therapy failure in prostate cancer

This study addresses the roles of nuclear receptor corepressor 2 (NCOR2) in prostate cancer (PC) progression in response to androgen deprivation therapy (ADT). Reduced NCOR2 expression significantly associates with shorter disease-free survival in patients with PC receiving adjuvant ADT. Utilizing the CWR22 xenograft model, we demonstrate that stably reduced NCOR2 expression accelerates disease recurrence following ADT, associates with gene expression patterns that include neuroendocrine features, and induces DNA hypermethylation. Stably reduced NCOR2 expression in isogenic LNCaP (androgen-sensitive) and LNCaP-C4–2 (androgen-independent) cells revealed that NCOR2 reduction phenocopies the impact of androgen treatment and induces global DNA hypermethylation patterns. NCOR2 genomic binding is greatest in LNCaP-C4–2 cells and most clearly associates with forkhead box (FOX) transcription factor FOXA1 binding. NCOR2 binding significantly associates with transcriptional regulation most when in active enhancer regions. These studies reveal robust roles for NCOR2 in regulating the PC transcriptome and epigenome and underscore recent mutational studies linking NCOR2 loss of function to PC disease progression.

Long et al. show that reduced levels of NCOR2 lead to accelerated prostate cancer recurrence during androgen withdrawal in a patient-derived xenograft model. NCOR2 reduction is characterized by incomplete response to androgen withdrawal, and recurrent tumors show increased neuroendocrine traits. These phenotypic changes are associated with hypermethylated enhancers.

GPS2 promotes erythroid differentiation by control of the stability of EKLF protein

Erythropoiesis is a complex multistage process that involves differentiation of early erythroid progenitors to enucleated mature red blood cells, in which lineage-specific transcription factors play essential roles. Erythroid Krüppel-like factor (EKLF/KLF1) is a pleiotropic erythroid transcription factor that is required for the proper maturation of the erythroid cells, whose expression and activation are tightly controlled in a temporal and differentiation stage-specific manner. Here, we uncover a novel role of G-protein pathway suppressor 2 (GPS2), a subunit of the nuclear receptor corepressor/silencing mediator of retinoic acid and thyroid hormone receptor corepressor complex, in erythrocyte differentiation. Our study demonstrates that knockdown of GPS2 significantly suppresses erythroid differentiation of human CD34+ cells cultured in vitro and xenotransplanted in nonobese diabetic/severe combined immunodeficiency/interleukin-2 receptor γ-chain null mice. Moreover, global deletion of GPS2 in mice causes impaired erythropoiesis in the fetal liver and leads to severe anemia. Flow cytometric analysis and Wright-Giemsa staining show a defective differentiation at late stages of erythropoiesis in Gps2-/- embryos. Mechanistically, GPS2 interacts with EKLF and prevents proteasome-mediated degradation of EKLF, thereby increasing EKLF stability and transcriptional activity. Moreover, we identify the amino acids 191-230 region in EKLF protein, responsible for GPS2 binding, that is highly conserved in mammals and essential for EKLF protein stability. Collectively, our study uncovers a previously unknown role of GPS2 as a posttranslational regulator that enhances the stability of EKLF protein and thereby promotes erythroid differentiation.

The nuclear receptor corepressor NCoR1 regulates hematopoiesis and leukemogenesis in vivo

Enhanced understanding of normal and malignant hematopoiesis pathways should facilitate the development of effective clinical treatment strategies for hematopoietic malignancies. Nuclear receptor corepressor 1 (NCoR1) has been implicated in transcriptional repression and embryonic organ development, but its role in hematopoiesis is yet to be fully elucidated. Here, we showed that hematopoietic-specific loss of NCoR1 leads to expansion of the hematopoietic stem cell (HSC) pool due to aberrant cell cycle entry of long-term HSCs under steady-state conditions. Moreover, NCoR1-deficient HSCs exhibited normal self-renewal capacity but severely impaired lymphoid-differentiation potential in competitive hematopoietic-reconstitution assays. Transcriptome analysis further revealed that several hematopoiesis-associated genes are regulated by NCoR1. In addition, NCoR1 deficiency in hematopoietic cells delayed the course of leukemia and promoted leukemia cell differentiation in an MLL-AF9-induced mouse model. NCoR1 and its partner, histone deacetylase 3, can modulate histone acetylation and gene transcription through binding the promoter regions of myeloid-differentiation genes. Our collective results support the critical involvement of NCoR1 in normal and malignant hematopoiesis in vivo.

LPCAT1 promotes brain metastasis of lung adenocarcinoma by up-regulating PI3K/AKT/MYC pathway

Background

Brain metastasis (BM) is associated with poor prognosis, recurrence, and death in patients with non-small cell lung cancer (NSCLC). Lysophosphatidylcholine acyltransferase 1 (LPCAT1) has been reported to be involved in the progression, metastasis and recurrence of malignancies. However, the potential role of LPCAT1 in NSCLC remains poorly understood. This study was aimed to identify genes involved in lung adenocarcinoma (LUAD) brain metastasis, and look into the role of LPCAT1 in LUAD progression.

Methods

We used integrative genomic analysis to identify genes involved in lung adenocarcinomas. LPCAT1 expression was evaluated in tumor tissues from LUAD patients and LUAD cell lines. The role of LPCAT1 was subsequently investigated both in vitro and in vivo. The mechanism underlying the involvement of LPCAT1 in LUAD progression was explored with the activator of PI3K/AKT pathway. RNA sequencing was performed to confirm the involvement of LPCAT1 and associated pathway in LUAD brain metastasis.

Results

LPCAT1 was up-regulated in LUAD tissues and cell lines. shRNA-mediated depletion of LPCAT1 not only abrogated cell proliferation, migration and invasion in vitro, but also arrested tumor growth and brain metastases in vivo. Notably, LPCAT1 at least partially influenced LUAD progression through PI3K/AKT signal pathway by targeting MYC transcription. Moreover, expression of LPCAT1 was higher in tissues of LUAD patients with BM than those without BM as revealed by IHC staining, RNA-Sequencing and qPCR analysis. Finally, elevated LPCAT1 expression in patients with lung adenocarcinomas was associated with a poor clinical outcome.

Conclusions

This study showed that LPCAT1 works as a regulator of cell metastasis and may serve as a novel therapeutic target for BM in lung adenocarcinoma.

Electronic supplementary material

The online version of this article (10.1186/s13046-019-1092-4) contains supplementary material, which is available to authorized users.

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.

Principled multi-omic analysis reveals gene regulatory mechanisms of phenotype variation

Recent studies have analyzed large-scale data sets of gene expression to identify genes associated with interindividual variation in phenotypes ranging from cancer subtypes to drug sensitivity, promising new avenues of research in personalized medicine. However, gene expression data alone is limited in its ability to reveal cis-regulatory mechanisms underlying phenotypic differences. In this study, we develop a new probabilistic model, called pGENMi, that integrates multi-omic data to investigate the transcriptional regulatory mechanisms underlying interindividual variation of a specific phenotype—that of cell line response to cytotoxic treatment. In particular, pGENMi simultaneously analyzes genotype, DNA methylation, gene expression, and transcription factor (TF)-DNA binding data, along with phenotypic measurements, to identify TFs regulating the phenotype. It does so by combining statistical information about expression quantitative trait loci (eQTLs) and expression-correlated methylation marks (eQTMs) located within TF binding sites, as well as observed correlations between gene expression and phenotype variation. Application of pGENMi to data from a panel of lymphoblastoid cell lines treated with 24 drugs, in conjunction with ENCODE TF ChIP data, yielded a number of known as well as novel (TF, Drug) associations. Experimental validations by TF knockdown confirmed 41% of the predicted and tested associations, compared to a 12% confirmation rate of tested nonassociations (controls). An extensive literature survey also corroborated 62% of the predicted associations above a stringent threshold. Moreover, associations predicted only when combining eQTL and eQTM data showed higher precision compared to an eQTL-only or eQTM-only analysis using pGENMi, further demonstrating the value of multi-omic integrative analysis.

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.

Global gene expression profiles of hematopoietic stem and progenitor cells from patients with chronic myeloid leukemia: the effect of in vitro culture with or without imatinib

In this study, we determined the gene expression profiles of bone marrow‐derived cell fractions, obtained from normal subjects and Chronic Myeloid Leukemia (CML) patients, that were highly enriched for hematopoietic stem (HSCs) and progenitor (HPCs) cells. Our results indicate that the profiles of CML HSCs and HPCs were closer to that of normal progenitors, whereas normal HSCs showed the most different expression profile of all. We found that the expression profiles of HSCs and HPCs from CML marrow were closer to each other than those of HSCs and HPCs from normal marrow. The major biologic processes dysregulated in CML cells included DNA repair, cell cycle, chromosome condensation, cell adhesion, and the immune response. We also determined the genomic changes in both normal and CML progenitor cells under culture conditions, and found that several genes involved in cell cycle, steroid biosynthesis, and chromosome segregation were upregulated, whereas genes involved in transcription regulation and apoptosis were downregulated. Interestingly, these changes were the same, regardless of the addition of Imatinib (IM) to the culture. Finally, we identified three genes—PIEZO2, RXFP1, and MAMDC2‐ that are preferentially expressed by CML primitive cells and that encode for cell membrane proteins; thus, they could be used as biomarkers for CML stem cells.

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.

Touch and go: nuclear proteolysis in the regulation of metabolic genes and cancer

The recruitment of transcription factors to promoters and enhancers is a critical step in gene regulation. Many of these proteins are quickly removed from DNA after they completed their function. Metabolic genes in particular are dynamically regulated and continuously adjusted to cellular requirements. Transcription factors controlling metabolism are therefore under constant surveillance by the ubiquitin–proteasome system, which can degrade DNA‐bound proteins in a site‐specific manner. Several of these metabolic transcription factors are critical to cancer cells, as they promote uncontrolled growth and proliferation. This review highlights recent findings in the emerging field of nuclear proteolysis and outlines novel paradigms for cancer treatment, with an emphasis on multiple myeloma.

Pan-cancer analyses of the nuclear receptor superfamily

Nuclear receptors (NR) act as an integrated conduit for environmental and hormonal signals to govern genomic responses, which relate to cell fate decisions. We review how their integrated actions with each other, shared co-factors and other transcription factors are disrupted in cancer. Steroid hormone nuclear receptors are oncogenic drivers in breast and prostate cancer and blockade of signaling is a major therapeutic goal. By contrast to blockade of receptors, in other cancers enhanced receptor function is attractive, as illustrated initially with targeting of retinoic acid receptors in leukemia. In the post-genomic era large consortia, such as The Cancer Genome Atlas, have developed a remarkable volume of genomic data with which to examine multiple aspects of nuclear receptor status in a pan-cancer manner. Therefore to extend the review of NR function we have also undertaken bioinformatics analyses of NR expression in over 3000 tumors, spread across six different tumor types (bladder, breast, colon, head and neck, liver and prostate). Specifically, to ask how the NR expression was distorted (altered expression, mutation and CNV) we have applied bootstrapping approaches to simulate data for comparison, and also compared these NR findings to 12 other transcription factor families. Nuclear receptors were uniquely and uniformly downregulated across all six tumor types, more than predicted by chance. These approaches also revealed that each tumor type had a specific NR expression profile but these were most similar between breast and prostate cancer. Some NRs were down-regulated in at least five tumor types (e.g. NR3C2/MR and NR5A2/LRH-1)) whereas others were uniquely down-regulated in one tumor (e.g. NR1B3/RARG). The downregulation was not driven by copy number variation or mutation and epigenetic mechanisms maybe responsible for the altered nuclear receptor expression.