FOXA1 defines cancer cell specificity

A novel nuclear receptor NR1D1 suppresses HSD17B12 transcription to regulate granulosa cell apoptosis and autophagy via the AMPK pathway in sheep

Dominant follicular development and atresia are governed by the proliferation of granulosa cells (GCs), a process influenced by the delicate balance between apoptosis and autophagy. Oxidative stress, a pivotal catalyst of GCs apoptosis, modulates gene expression through epigenetic mechanisms, including chromatin remodeling. Nevertheless, the regulatory mechanisms underpinning GCs functionality in relation to prolificacy remain inadequately elucidated. In this study, we discovered that the chromatin accessibility of nuclear receptor subfamily 1 group D member 1 (NR1D1) was markedly enhanced in dominant follicular GCs from low-prolificacy sheep, as evidenced by Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq), which correlated with elevated NR1D1 transcript levels. Remarkably, NR1D1 emerged as a novel regulator of follicular development, exhibiting heightened expression in dominant follicles. The overexpression of NR1D1 induced cell cycle arrest, autophagy activation, and mitochondrial dysfunction via the AMPK pathway, while its knockdown fostered GCs survival and functionality. Furthermore, NR1D1 inhibits the transcription of HSD17B12, thereby contributing to oxidative stress (ROS)-induced apoptosis, as demonstrated by CUT&Tag-qPCR and dual luciferase assays. The downregulation of HSD17B12 partially alleviated the effects of NR1D1 knockdown on GCs functionality. These findings indicate that NR1D1 orchestrates GCs proliferation and apoptosis through the suppression of HSD17B12 and the activation of the AMPK pathway, establishing NR1D1 as a novel transcription factor implicated in follicular development and ovarian function, with significant implications for prolificacy.

Amogel: a multi-omics classification framework using associative graph neural networks with prior knowledge for biomarker identification

The advent of high-throughput sequencing technologies, such as DNA microarray and DNA sequencing, has enabled effective analysis of cancer subtypes and targeted treatment. Furthermore, numerous studies have highlighted the capability of graph neural networks (GNN) to model complex biological systems and capture non-linear interactions in high-throughput data. GNN has proven to be useful in leveraging multiple types of omics data, including prior biological knowledge from various sources, such as transcriptomics, genomics, proteomics, and metabolomics, to improve cancer classification. However, current works do not fully utilize the non-linear learning potential of GNN and lack of the integration ability to analyse high-throughput multi-omics data simultaneously with prior biological knowledge. Nevertheless, relying on limited prior knowledge in generating gene graphs might lead to less accurate classification due to undiscovered significant gene-gene interactions, which may require expert intervention and can be time-consuming. Hence, this study proposes a graph classification model called associative multi-omics graph embedding learning (AMOGEL) to effectively integrate multi-omics datasets and prior knowledge through GNN coupled with association rule mining (ARM). AMOGEL employs an early fusion technique using ARM to mine intra-omics and inter-omics relationships, forming a multi-omics synthetic information graph before the model training. Moreover, AMOGEL introduces multi-dimensional edges, with multi-omics gene associations or edges as the main contributors and prior knowledge edges as auxiliary contributors. Additionally, it uses a gene ranking technique based on attention scores, considering the relationships between neighbouring genes. Several experiments were performed on BRCA and KIPAN cancer subtypes to demonstrate the integration of multi-omics datasets (miRNA, mRNA, and DNA methylation) with prior biological knowledge of protein-protein interactions, KEGG pathways and Gene Ontology. The experimental results showed that the AMOGEL outperformed the current state-of-the-art models in terms of classification accuracy, F1 score and AUC score. The findings of this study represent a crucial step forward in advancing the effective integration of multi-omics data and prior knowledge to improve cancer subtype classification.

Pan-cancer genomic analysis reveals FOXA1 amplification is associated with adverse outcomes in non-small cell lung, prostate, and breast cancers

BACKGROUND

Alterations in forkhead box A1 (FOXA1), a pioneer transcription factor, are associated with poor prognosis in breast cancer and prostate cancer. We characterized FOXA1 genomic alterations and their clinical impacts in a large pan-cancer cohort from the American Association for Cancer Research Genomics, Evidence, Neoplasia, Information, Exchange database.

METHODS

FOXA1 alterations were characterized across more than 87 000 samples from more than 30 cancer types for primary and metastatic tumors alongside patient characteristics and clinical outcomes. FOXA1 alterations were queried in the Memorial Sloan Kettering - Metastatic Events and Tropisms (MSK-MET) cohort (a GENIE subset), allowing definition of hazard ratios (HRs) and survival estimates based on Cox proportional hazard models.

RESULTS

FOXA1 was altered in 1869 (2.1%) samples, with distinct patterns across different cancers: prostate cancer enriched with indel-inframe alterations, breast cancer with missense mutations, and lung cancers with copy number amplifications. Of 74 715 samples with FOXA1 copy number profiles, amplification was detected in 834 (1.1%). Amplification was most common in non-small cell lung cancer (NSCLC; 3% in primary; 6% in metastatic) and small cell lung cancer (4.1% primary; 3.5% metastatic), followed by breast cancer (2% primary; 1.6% metastatic) and prostate cancer (2.2% primary; 1.6% metastatic). Copy number amplifications were associated with decreased overall survival in NSCLC (HR = 1.45, 95% confidence interval [CI] = 1.06 to 1.99; P = .02), breast cancer (HR = 3.04, 95% CI = 1.89 to 4.89; P = 4e-6), and prostate cancer (HR = 1.94, 95% CI = 1.03 to 3.68; P = .04). Amplifications were associated with widespread metastases in NSCLC, breast cancer, and prostate cancer.

CONCLUSIONS

FOXA1 demonstrates distinct alteration profiles across cancer sites. Our findings suggest an association between FOXA1 amplification and enhanced metastatic potential and decreased survival, highlighting prognostic and therapeutic potential in breast cancer, prostate cancer, and NSCLC.

The EstroGene2.0 database for endocrine therapy response and resistance in breast cancer

Endocrine therapies targeting the estrogen receptor (ER/ESR1) are the cornerstone to treat ER-positive breast cancers patients, but resistance often limits their effectiveness. Notable progress has been made although the fragmented way data is reported has reduced their potential impact. Here, we introduce EstroGene2.0, an expanded database of its precursor 1.0 version. EstroGene2.0 focusses on response and resistance to endocrine therapies in breast cancer models. Incorporating multi-omic profiling of 361 experiments from 212 studies across 28 cell lines, a user-friendly browser offers comprehensive data visualization and metadata mining capabilities ( https://estrogeneii.web.app/ ). Taking advantage of the harmonized data collection, our follow-up meta-analysis revealed transcriptomic landscape and substantial diversity in response to different classes of ER modulators. Endocrine-resistant models exhibit a spectrum of transcriptomic alterations including a contra-directional shift in ER and interferon signalings, which is recapitulated clinically. Dissecting multiple ESR1-mutant cell models revealed the different clinical relevance of cell model engineering and identified high-confidence mutant-ER targets, such as NPY1R. These examples demonstrate how EstroGene2.0 helps investigate breast cancer's response to endocrine therapies and explore resistance mechanisms.

EstroGene2.0: A multi-omic database of response to estrogens, ER-modulators, and resistance to endocrine therapies in breast cancer

Endocrine therapies targeting the estrogen receptor (ER/ESR1) are the cornerstone to treat ER-positive breast cancers patients, but resistance often limits their effectiveness. Understanding the molecular mechanisms is thus key to optimize the existing drugs and to develop new ER-modulators. Notable progress has been made although the fragmented way data is reported has reduced their potential impact. Here, we introduce EstroGene2.0, an expanded database of its precursor 1.0 version. EstroGene2.0 focusses on response and resistance to endocrine therapies in breast cancer models. Incorporating multi-omic profiling of 361 experiments from 212 studies across 28 cell lines, a user-friendly browser offers comprehensive data visualization and metadata mining capabilities (https://estrogeneii.web.app/). Taking advantage of the harmonized data collection, our follow-up meta-analysis revealed substantial diversity in response to different classes of ER-modulators including SERMs, SERDs, SERCA and LDD/PROTAC. Notably, endocrine resistant models exhibit a spectrum of transcriptomic alterations including a contra-directional shift in ER and interferon signaling, which is recapitulated clinically. Furthermore, dissecting multiple ESR1-mutant cell models revealed the different clinical relevance of genome-edited versus ectopic overexpression model engineering and identified high-confidence mutant-ER targets, such as NPY1R. These examples demonstrate how EstroGene2.0 helps investigate breast cancer's response to endocrine therapies and explore resistance mechanisms.

Cis-regulatory control of transcriptional timing and noise in response to estrogen

Cis-regulatory elements control transcription levels, temporal dynamics, and cell-cell variation or transcriptional noise. However, the combination of regulatory features that control these different attributes is not fully understood. Here, we used single-cell RNA-seq during an estrogen treatment time course and machine learning to identify predictors of expression timing and noise. We found that genes with multiple active enhancers exhibit faster temporal responses. We verified this finding by showing that manipulation of enhancer activity changes the temporal response of estrogen target genes. Analysis of transcriptional noise uncovered a relationship between promoter and enhancer activity, with active promoters associated with low noise and active enhancers linked to high noise. Finally, we observed that co-expression across single cells is an emergent property associated with chromatin looping, timing, and noise. Overall, our results indicate a fundamental tradeoff between a gene's ability to quickly respond to incoming signals and maintain low variation across cells.

Cis-regulatory control of transcriptional timing and noise in response to estrogen

Cis-regulatory elements control transcription levels, temporal dynamics, and cell-cell variation or transcriptional noise. However, the combination of regulatory features that control these different attributes is not fully understood. Here, we used single cell RNA-seq during an estrogen treatment time course and machine learning to identify predictors of expression timing and noise. We find that genes with multiple active enhancers exhibit faster temporal responses. We verified this finding by showing that manipulation of enhancer activity changes the temporal response of estrogen target genes. Analysis of transcriptional noise uncovered a relationship between promoter and enhancer activity, with active promoters associated with low noise and active enhancers linked to high noise. Finally, we observed that co-expression across single cells is an emergent property associated with chromatin looping, timing, and noise. Overall, our results indicate a fundamental tradeoff between a gene's ability to quickly respond to incoming signals and maintain low variation across cells.

Role of Fork-Head Box Genes in Breast Cancer: From Drug Resistance to Therapeutic Targets

Breast cancer has been acknowledged as one of the most notorious cancers, responsible for millions of deaths around the globe. Understanding the various factors, genetic mutations, comprehensive pathways, etc., that are involved in the development of breast cancer and how these affect the development of the disease is very important for improving and revitalizing the treatment of this global health issue. The forkhead-box gene family, comprising 19 subfamilies, is known to have a significant impact on the growth and progression of this cancer. The article looks into the various forkhead genes and how they play a role in different types of cancer. It also covers their impact on cancer drug resistance, interaction with microRNAs, explores their potential as targets for drug therapies, and their association with stem cells.

ERα prevents tumorigenesis of both liver and breast cancer cells through CCN5

Estrogen receptor alpha (ERα)-mediated estrogen signaling has also shown to prevent hepatic tumorigenesis in mice. Consistent with this, hormone replacement therapy with estrogen supplementation dramatically reduced the risk of hepatocellular carcinoma. Silencing of ERα is also a key event for the transformation of ERα-positive breast cancer cells into malignant triple-negative breast cancer cells. However, the mechanisms underlying ERα-mediated prevention of both hepatic and mammary tumorigenesis in humans are still unclear. Here, we present a functional genomics study of ERα targeting by comparing human liver cancer cells with human breast cancer cells using "loss or gain of function" genetic assays of ERα in vitro and in vivo. We discover that cellular communication network factor 5 (CCN5) is a direct downstream target of ERα; ERα suppresses growth and prevents tumorigenesis and malignant transformation of both liver and breast cancer cells through CCN5 in humans. The ERα-CCN5 regulatory axis functions as suppressors for both hepatic and mammary tumors, which is a common mechanism of preventing tumorigenesis for both liver cancer and breast cancer in humans.

FOXA1/2 depletion drives global reprogramming of differentiation state and metabolism in a human liver cell line and inhibits differentiation of human stem cell-derived hepatic progenitor cells

FOXA factors are critical members of the developmental gene regulatory network (GRN) composed of master transcription factors (TF) which regulate murine cell fate and metabolism in the gut and liver. How FOXA factors dictate human liver cell fate, differentiation, and simultaneously regulate metabolic pathways is poorly understood. Here, we aimed to determine the role of FOXA2 (and FOXA1 which is believed to compensate for FOXA2) in controlling hepatic differentiation and cell metabolism in a human hepatic cell line (HepG2). siRNA mediated knockdown of FOXA1/2 in HepG2 cells significantly downregulated albumin (p < .05) and GRN TF gene expression (HNF4α, HEX, HNF1ß, TBX3) (p < .05) and significantly upregulated endoderm/gut/hepatic endoderm markers (goosecoid [GSC], FOXA3, and GATA4), gut TF (CDX2), pluripotent TF (NANOG), and neuroectodermal TF (PAX6) (p < .05), all consistent with partial/transient reprograming. shFOXA1/2 targeting resulted in similar findings and demonstrated evidence of reversibility of phenotype. RNA-seq followed by bioinformatic analysis of shFOXA1/2 knockdown HepG2 cells demonstrated 235 significant downregulated genes and 448 upregulated genes, including upregulation of markers for alternate germ layers lineages (cardiac, endothelial, muscle) and neurectoderm (eye, neural). We found widespread downregulation of glycolysis, citric acid cycle, mitochondrial genes, and alterations in lipid metabolism, pentose phosphate pathway, and ketogenesis. Functional metabolic analysis agreed with these findings, demonstrating significantly diminished glycolysis and mitochondrial respiration, with concomitant accumulation of lipid droplets. We hypothesized that FOXA1/2 inhibit the initiation of human liver differentiation in vitro. During human pluripotent stem cells (hPSC)-hepatic differentiation, siRNA knockdown demonstrated de-differentiation and unexpectedly, activation of pluripotency factors and neuroectoderm. shRNA knockdown demonstrated similar results and activation of SOX9 (hepatobiliary). These results demonstrate that FOXA1/2 controls hepatic and developmental GRN, and their knockdown leads to reprogramming of both differentiation and metabolism, with applications in studies of cancer, differentiation, and organogenesis.

FOXA1 in Breast Cancer: A Luminal Marker with Promising Prognostic and Predictive Impact

The present review focuses on the function of the forkhead protein FOXA1 in breast cancer (BC) in relation to steroid hormone receptors. We explored the currently available analytic approaches for FOXA1 assessment both at gene and protein levels, comparing the differences between the available techniques used for its diagnostic assessment. In addition, we elaborated on data regarding the prognostic and predictive role of this marker in BC based on several studies that evaluated its expression in relation to the outcome and/or response to therapy. FOXA1, similar to the androgen receptor (AR), may have a dual role in BC according to hormonal status. In luminal cancers, its expression contributes to a better prognosis, while in triple-negative breast cancers (TNBC), it implies an adverse outcome. Consequently, we observed that FOXA1-positive expression in a neoadjuvant setting may predict a lack of response in luminal BC as opposed to TNBC, in which FOXA1 allegedly increases its chemosensitivity. In conclusion, considering its accessible and convenient identification by immunohistochemistry, its important impact on prognosis, and its suitability to identify patients with different responses to chemotherapy, we propose that FOXA1 could be tested in routine diagnostics as an additional prognostic and predictive marker in BC.

Elevated FOXA1 Expression Indicates Poor Prognosis in Liver Cancer due to Its Effects on Cell Proliferation and Metastasis

Purpose

Studying the pathogenesis of liver cancer is conducive to the exploration of effective diagnostic and prognostic biomarkers. In this study, we investigated the expression of FOXA1 and its oncogenic role in hepatocellular carcinoma (HCC).

Methods

Transcriptome data of HCC tissues were downloaded from The Cancer Genome Atlas (TCGA) and GEO databases and analyzed using R software. We also upregulated FOXA1 expression in HCC cells and investigated the role of FOXA1 in the proliferation and migration of HCC cells through proliferation, colony formation, wound healing, and Transwell assays.

Results

An analysis of the transcriptome data in TCGA database revealed found that FOXA1 is highly expressed in HCC tissues and that patients with low FOXA1 expression have a better prognosis. High FOXA1 expression was mainly associated with extracellular matrix organization, cancer, and mitosis. The results of an immunohistochemistry (IHC) assay showed that FOXA1 protein was highly expressed in HCC tissues, and patients with low FOXA1 expression showed longer disease-specific survival times and progression-free intervals. The results from quantitative reverse transcription–PCR (RT–qPCR) and Western blot experiments showed that the expression of FOXA1 in liver cancer cell lines was higher than that in immortalized human liver cell lines. Proliferation, wound healing, and Transwell experiments showed that FOXA1 enhanced the proliferation and migration abilities of liver cancer and immortalized human cell lines.

Conclusion

Our research suggests that FOXA1 plays an important role in promoting the recurrence and metastasis of HCC by increasing cell proliferation and metastasis.

CRISPR/Cas9 system in breast cancer therapy: advancement, limitations and future scope

Cancer is one of the major causes of mortality worldwide, therefore it is considered a major health concern. Breast cancer is the most frequent type of cancer which affects women on a global scale. Various current treatment strategies have been implicated for breast cancer therapy that includes surgical removal, radiation therapy, hormonal therapy, chemotherapy, and targeted biological therapy. However, constant effort is being made to introduce novel therapies with minimal toxicity. Gene therapy is one of the promising tools, to rectify defective genes and cure various cancers. In recent years, a novel genome engineering technology, namely the clustered regularly interspaced short palindromic repeat (CRISPR)-associated protein-9 (Cas9) has emerged as a gene-editing tool and transformed genome-editing techniques in a wide range of biological domains including human cancer research and gene therapy. This could be attributed to its versatile characteristics such as high specificity, precision, time-saving and cost-effective methodologies with minimal risk. In the present review, we highlight the role of CRISPR/Cas9 as a targeted therapy to tackle drug resistance, improve immunotherapy for breast cancer.

CREBBP/EP300 acetyltransferase inhibition disrupts FOXA1-bound enhancers to inhibit the proliferation of ER+ breast cancer cells

Therapeutic targeting of the estrogen receptor (ER) is a clinically validated approach for estrogen receptor positive breast cancer (ER+ BC), but sustained response is limited by acquired resistance. Targeting the transcriptional coactivators required for estrogen receptor activity represents an alternative approach that is not subject to the same limitations as targeting estrogen receptor itself. In this report we demonstrate that the acetyltransferase activity of coactivator paralogs CREBBP/EP300 represents a promising therapeutic target in ER+ BC. Using the potent and selective inhibitor CPI-1612, we show that CREBBP/EP300 acetyltransferase inhibition potently suppresses in vitro and in vivo growth of breast cancer cell line models and acts in a manner orthogonal to directly targeting ER. CREBBP/EP300 acetyltransferase inhibition suppresses ER-dependent transcription by targeting lineage-specific enhancers defined by the pioneer transcription factor FOXA1. These results validate CREBBP/EP300 acetyltransferase activity as a viable target for clinical development in ER+ breast cancer.

Compressive stress-mediated p38 activation required for ERα + phenotype in breast cancer

Breast cancer is now globally the most frequent cancer and leading cause of women's death. Two thirds of breast cancers express the luminal estrogen receptor-positive (ERα + ) phenotype that is initially responsive to antihormonal therapies, but drug resistance emerges. A major barrier to the understanding of the ERα-pathway biology and therapeutic discoveries is the restricted repertoire of luminal ERα + breast cancer models. The ERα + phenotype is not stable in cultured cells for reasons not fully understood. We examine 400 patient-derived breast epithelial and breast cancer explant cultures (PDECs) grown in various three-dimensional matrix scaffolds, finding that ERα is primarily regulated by the matrix stiffness. Matrix stiffness upregulates the ERα signaling via stress-mediated p38 activation and H3K27me3-mediated epigenetic regulation. The finding that the matrix stiffness is a central cue to the ERα phenotype reveals a mechanobiological component in breast tissue hormonal signaling and enables the development of novel therapeutic interventions. Subject terms: ER-positive (ER + ), breast cancer, ex vivo model, preclinical model, PDEC, stiffness, p38 SAPK.

KDM6A-ARHGDIB axis blocks metastasis of bladder cancer by inhibiting Rac1

Background

KDM6A, a histone demethylase, is frequently mutated in bladder cancer (BCa). However, the role and detailed molecular mechanism of KDM6A involved in bladder cancer progression remains unknown.

Methods

Tissue specimens were used to determine the expression levels and prognostic values of KDM6A and ARHGDIB. The MTT, colony formation, wound healing and Transwell migration and invasion assays were employed to detect the BCa cell proliferation, migration and invasion, respectively. Chemotaxis of macrophages was used to evaluate the ability of KDM6A to recruit macrophages. A subcutaneous tumour model and tail vein tumour injection in nude mice were used to assess the role of KDM6A in vivo. RNA sequencing, qPCR, Western blot, ChIP and phalloidin staining assay were performed to investigate the molecular functions of KDM6A. Dual-luciferase reporter assay was used to determine the effects of KDM6A and FOXA1 on the promoters of the ARHGDIB and KDM6A.

Results

We showed that the KDM6A inhibited the motility and invasiveness of the BCa cells. Mechanistically, KDM6A promotes the transcription of ARHGDIB by demethylating histone H3 lysine di/trimethylation (H3K27me2/3) and consequently leads to inhibition of Rac1. EZH2, which catalyses the methylation of H3K27, functions to silence ARHGDIB expression, and an EZH2 inhibitor can neutralize the metastatic effect caused by KDM6A deficiency. Furthermore, we demonstrated that FOXA1 directly binds to the KDM6A promoter and thus transactivates KDM6A, leading to diminished metastatic potential.

Conclusion

Our findings establish the critical role of the FOXA1-KDM6A-ARHGDIB axis in restraining the malignancy of BCa and identify KDM6A and EZH2 as potential therapeutic targets in the management of BCa.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12943-021-01369-9.

Knockdown of Forkhead box A1 suppresses the tumorigenesis and progression of human colon cancer cells through regulating the phosphatase and tensin homolog/Akt pathway

Background

This study aimed to evaluate the role and the underlying mechanisms of Forkhead box A1 (encoded by FOXA1) in colon cancer.

Methods

We analyzed FOXA1 mRNA and protein expression in colon cancer tissues and cell lines. We also silenced FOXA1 expression in HCT116 and SW480 cells to evaluate the effects on cell proliferation, cell cycle, migration, and invasion by using MTT, colony formation, flow cytometry, and the Transwell assay, respectively.

Results

FOXA1 immunostaining was higher in colon cancer tissues than adjacent healthy tissues. FOXA1 mRNA and protein expression was significantly increased in human colon cancer cells compared with a normal colonic cell line. FOXA1 expression was also significantly higher in colorectal cancer tissues from TCGA data sets and was associated with worse prognosis in the R2 database. FOXA1 expression was negatively correlated with the extent of its methylation, and its knockdown reduced proliferation, migration, and invasion, and induced G2/M phase arrest in HCT116 and SW480 cells by suppressing the phosphatase and tensin homolog/Akt signaling pathway and inhibiting epithelial–mesenchymal transition.

Conclusion

FOXA1 may act as an oncogene in colon cancer tumorigenesis and development.

Subtype-associated epigenomic landscape and 3D genome structure in bladder cancer

Muscle-invasive bladder cancers are characterized by their distinct expression of luminal and basal genes, which could be used to predict key clinical features such as disease progression and overall survival. Transcriptionally, FOXA1, GATA3, and PPARG are shown to be essential for luminal subtype-specific gene regulation and subtype switching, while TP63, STAT3, and TFAP2 family members are critical for regulation of basal subtype-specific genes. Despite these advances, the underlying epigenetic mechanisms and 3D chromatin architecture responsible for subtype-specific regulation in bladder cancer remain unknown. RESULT: We determine the genome-wide transcriptome, enhancer landscape, and transcription factor binding profiles of FOXA1 and GATA3 in luminal and basal subtypes of bladder cancer. Furthermore, we report the first-ever mapping of genome-wide chromatin interactions by Hi-C in both bladder cancer cell lines and primary patient tumors. We show that subtype-specific transcription is accompanied by specific open chromatin and epigenomic marks, at least partially driven by distinct transcription factor binding at distal enhancers of luminal and basal bladder cancers. Finally, we identify a novel clinically relevant transcription factor, Neuronal PAS Domain Protein 2 (NPAS2), in luminal bladder cancers that regulates other subtype-specific genes and influences cancer cell proliferation and migration. CONCLUSION: In summary, our work identifies unique epigenomic signatures and 3D genome structures in luminal and basal urinary bladder cancers and suggests a novel link between the circadian transcription factor NPAS2 and a clinical bladder cancer subtype.

Genome-Wide Analysis of the FOXA1 Transcriptional Network Identifies Novel Protein-Coding and Long Noncoding RNA Targets in Colorectal Cancer Cells

Differentiation status of tumors is correlated with metastatic potential and malignancy. FOXA1 (forkhead box A1) is a transcription factor known to regulate differentiation in certain tissues. Here, we investigate FOXA1 function in human colorectal cancer (CRC). We found that FOXA1 is robustly expressed in the normal human colon but significantly downregulated in colon adenocarcinoma. Applying FOXA1 chromatin immunoprecipitation coupled with deep sequencing and transcriptome analysis upon FOXA1 knockdown in well-differentiated CRC cells and FOXA1 overexpression in poorly differentiated CRC cells, we identified novel protein-coding and lncRNA genes regulated by FOXA1. Among the numerous novel FOXA1 targets we identified, we focused on CEACAM5, a tumor marker and facilitator of cell adhesion. We show that FOXA1 binds to a distal enhancer downstream of CEACAM5 and strongly activates its expression. Consistent with these data, we show that FOXA1 inhibits anoikis in CRC cells. Collectively, our results uncover novel protein-coding and noncoding targets of FOXA1 and suggest a vital role of FOXA1 in enhancing CEACAM5 expression and anoikis resistance in CRC cells.

TADs enriched in histone H1.2 strongly overlap with the B compartment, inaccessible chromatin, and AT-rich Giemsa bands

Giemsa staining of metaphase chromosomes results in a characteristic banding useful for identification of chromosomes and its alterations. We have investigated in silico whether Giemsa bands (G bands) correlate with epigenetic and topological features of the interphase genome. Staining of G-positive bands decreases with GC content; nonetheless, G-negative bands are GC heterogeneous. High GC bands are enriched in active histone marks, RNA polymerase II, and SINEs and associate with gene richness, gene expression, and early replication. Low GC bands are enriched in repressive marks, lamina-associated domains, and LINEs. Histone H1 variants distribute heterogeneously among G bands: H1X is enriched at high GC bands and H1.2 is abundant at low GC, compacted bands. According to epigenetic features and H1 content, G bands can be organized in clusters useful to compartmentalize the genome. Indeed, we have obtained Hi-C chromosome interaction maps and compared topologically associating domains (TADs) and A/B compartments to G banding. TADs with high H1.2/H1X ratio strongly overlap with B compartment, late replicating, and inaccessible chromatin and low GC bands. We propose that GC content is a strong driver of chromatin compaction and 3D genome organization, that Giemsa staining recapitulates this organization denoted by high-throughput techniques, and that H1 variants distribute at distinct chromatin domains. DATABASES: Hi-C data on T47D breast cancer cells have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE147627.

The ubiquitinase ZFP91 promotes tumor cell survival and confers chemoresistance through FOXA1 destabilization

The forkhead box A1 (FOXA1), one of the forkhead class of DNA-binding proteins, functions as a transcription factor and plays a vital role in cellular control of embryonic development and cancer progression. Downregulation of FOXA1 has reported in several types of cancer, which contributes to cancer cell survival and chemoresistance. However, the mechanism for FOXA1 downregulation in cancer remains unclear. Here, we report that the ubiquitination enzyme zinc finger protein 91 (ZFP91) ubiquitinates and destabilizes FOXA1, which promotes cancer cell growth. High level of ZFP91 expression correlates with low level of FOXA1 protein in human gastric cancer (GC) cell lines and patient samples. Furthermore, ZFP91 knockdown reduces FOXA1 polyubiquitination, which decreases FOXA1 turnover and enhances cellular sensitivity to chemotherapy. Taken together, our findings reveal ZFP91-FOXA1 axis plays an important role in promoting GC progression and provides us a potential therapeutic intervention in the treatment of GC.

Shaping Chromatin States in Prostate Cancer by Pioneer Transcription Factors

The androgen receptor (AR) is a critical therapeutic target in prostate cancer that responds to antagonists in primary disease, but inevitably becomes reactivated, signaling onset of the lethal castration-resistant prostate cancer (CRPC) stage. Epigenomic investigation of the chromatin environment and interacting partners required for AR transcriptional activity has uncovered three pioneer factors that open up chromatin and facilitate AR-driven transcriptional programs. FOXA1, HOXB13, and GATA2 are required for normal AR transcription in prostate epithelial development and for oncogenic AR transcription during prostate carcinogenesis. AR signaling is dependent upon these three pioneer factors both before and after the clinical transition from treatable androgen-dependent disease to untreatable CRPC. Agents targeting their respective DNA binding or downstream chromatin-remodeling events have shown promise in preclinical studies of CRPC. AR-independent functions of FOXA1, HOXB13, and GATA2 are emerging as well. While all three pioneer factors exert effects that promote carcinogenesis, some of their functions may inhibit certain stages of prostate cancer progression. In all, these pioneer factors represent some of the most promising potential therapeutic targets to emerge thus far from the study of the prostate cancer epigenome.

Noncoding mutations target cis-regulatory elements of the FOXA1 plexus in prostate cancer

Prostate cancer is the second most commonly diagnosed malignancy among men worldwide. Recurrently mutated in primary and metastatic prostate tumors, FOXA1 encodes a pioneer transcription factor involved in disease onset and progression through both androgen receptor-dependent and androgen receptor-independent mechanisms. Despite its oncogenic properties however, the regulation of FOXA1 expression remains unknown. Here, we identify a set of six cis-regulatory elements in the FOXA1 regulatory plexus harboring somatic single-nucleotide variants in primary prostate tumors. We find that deletion and repression of these cis-regulatory elements significantly decreases FOXA1 expression and prostate cancer cell growth. Six of the ten single-nucleotide variants mapping to FOXA1 regulatory plexus significantly alter the transactivation potential of cis-regulatory elements by modulating the binding of transcription factors. Collectively, our results identify cis-regulatory elements within the FOXA1 plexus mutated in primary prostate tumors as potential targets for therapeutic intervention.

FOXA1 pioneer transcription factor is recurrently mutated in primary and metastatic prostate tumors. Here, authors identify a set of six cis-regulatory elements in the FOXA1 regulatory plexus harboring somatic SNVs in primary prostate tumors and characterize their role in regulating FOXA1 expression and prostate cancer cell growth.

Chromatin interactome mapping at 139 independent breast cancer risk signals

BACKGROUND

Genome-wide association studies have identified 196 high confidence independent signals associated with breast cancer susceptibility. Variants within these signals frequently fall in distal regulatory DNA elements that control gene expression.

RESULTS

We designed a Capture Hi-C array to enrich for chromatin interactions between the credible causal variants and target genes in six human mammary epithelial and breast cancer cell lines. We show that interacting regions are enriched for open chromatin, histone marks for active enhancers, and transcription factors relevant to breast biology. We exploit this comprehensive resource to identify candidate target genes at 139 independent breast cancer risk signals and explore the functional mechanism underlying altered risk at the 12q24 risk region.

CONCLUSIONS

Our results demonstrate the power of combining genetics, computational genomics, and molecular studies to rationalize the identification of key variants and candidate target genes at breast cancer GWAS signals.

The Many Faces of Gene Regulation in Cancer: A Computational Oncogenomics Outlook

Cancer is a complex disease at many different levels. The molecular phenomenology of cancer is also quite rich. The mutational and genomic origins of cancer and their downstream effects on processes such as the reprogramming of the gene regulatory control and the molecular pathways depending on such control have been recognized as central to the characterization of the disease. More important though is the understanding of their causes, prognosis, and therapeutics. There is a multitude of factors associated with anomalous control of gene expression in cancer. Many of these factors are now amenable to be studied comprehensively by means of experiments based on diverse omic technologies. However, characterizing each dimension of the phenomenon individually has proven to fall short in presenting a clear picture of expression regulation as a whole. In this review article, we discuss some of the more relevant factors affecting gene expression control both, under normal conditions and in tumor settings. We describe the different omic approaches that we can use as well as the computational genomic analysis needed to track down these factors. Then we present theoretical and computational frameworks developed to integrate the amount of diverse information provided by such single-omic analyses. We contextualize this within a systems biology-based multi-omic regulation setting, aimed at better understanding the complex interplay of gene expression deregulation in cancer.

Cerumenogram: a new frontier in cancer diagnosis in humans

Cancer is the deadliest human disease and the development of new diagnosis methods is important to increase the chances of a cure. In this work it was developed a new method, named here for the first time as cerumenogram, using cerumen (earwax) as a new biomatrix for diagnosis. Earwax samples collected from cancer patients (cancer group) and cancer-free patients (control group) were analyzed by Headspace/Gas Chromatography-Mass Spectrometry (HS/GC-MS), following with multivariate analysis steps to process the raw data generated. In total, 158 volatile organic metabolites (VOMs) were identified in the cerumen samples. The 27 selected as potential VOMs biomarkers for cancer provided 100% discrimination between the cancer and control groups. This new test can thus be routinely employed for cancer diagnoses that is non-invasive, fast, cheap, and highly accurate.

Forkhead‑box A1 regulates tumor cell growth and predicts prognosis in colorectal cancer

Forkhead box A1 (FOXA1) functions as a tumor suppressor gene or an oncogene in various types of cancer; however, the distinct function of FOXA1 in colorectal cancer is unclear. The present study aimed to evaluate whether FOXA1 affects the oncogenic behavior of colorectal cancer cells, and to investigate its prognostic value in colorectal cancer. The impact of FOXA1 on tumor cell behavior was investigated using small interfering RNA and the pcDNA6‑myc vector in human colorectal cancer cell lines. To investigate the role of FOXA1 in the progression of human colorectal cancer, an immunohistochemical technique was used to localize FOXA1 protein in paraffin‑embedded tissue blocks obtained from 403 patients with colorectal cancer. Tumor cell apoptosis and proliferation were evaluated using a terminal deoxynucleotidyl transferase‑mediated dUTP nick‑end labeling assay and Ki‑67 immunohistochemical staining, respectively. FOXA1 knockdown inhibited tumor cell invasion in colorectal cancer cells, and induced apoptosis and cell cycle arrest. FOXA1 knockdown activated cleaved caspase‑poly (ADP‑ribose) polymerase, upregulated the expression of p53 upregulated modulator of apoptosis, and downregulated BH3 interacting domain death agonist and myeloid cell leukemia‑1, leading to the induction of apoptosis. FOXA1 knockdown increased the phosphorylation level of signal transducer and activator of tran-scription‑3. By contrast, these results were reversed following the overexpression of FOXA1. The overexpression of FOXA1 was associated with differentiation, lymphovascular invasion, advanced tumor stage, depth of invasion, lymph node metastasis and poor survival rate. The mean Ki‑67 labeling index value of FOXA1‑positive tumors was significantly higher than that of FOXA1‑negative tumors. However, no significant association was observed between the expression of FOXA1 and the mean apoptotic index value. These results indicate that FOXA1 is associated with tumor progression via the modulation of tumor cell survival in human colorectal cancer.

The prognostic relevance of FOXA1 and Nestin expression in breast cancer metastases: a retrospective study of 164 cases during a 10-year period (2004-2014)

Background

Current prognostic markers cannot adequately predict the clinical outcome of breast cancer patients. Therefore, additional biomarkers need to be included in routine immune panels. FOXA1 was a significant predictor of favorable outcome in primary breast cancer, while Nestin expression is preferentially found in triple-negative tumors with increased rate of nodal metastases, and reduced survival. No studies have investigated the prognostic value of FOXA1 and Nestin expression in breast cancer metastases.

Methods

Breast cancer metastases (n = 164) from various anatomical sites were retrospectively analyzed by immunohistochemistry for FOXA1, Nestin and GATA3 expression. Cox regression analysis assessed the prognostic value of FOXA1 and Nestin expression.

Results

In breast cancer metastases, FOXA1 expression was associated with Nestin-negativity, GATA3-positivity, ER-positivity, HER2-positivity and non-triple-negative status (P < 0.05). In contrast, Nestin expression was associated with FOXA1-negative, GATA3-negative, ER-negative, and triple-negative metastases (P < 0.05). Univariate Cox regression analysis showed FOXA1 expression was predictive of overall survival (OS, P = 0.00048) and metastasis-free survival (DMFS, P = 0.0011), as well as, distant metastasis-free survival in ER-positive patients (P = 0.036) and overall survival in ER-negative patients (P = 0.024). Multivariate analysis confirmed the significance of FOXA1 for both survival endpoints in metastatic breast cancer patients (OS, P = 0.0033; DMFS, P = 0.015).

Conclusions

In our study, FOXA1 was expressed mostly in ER-positive breast cancer metastases. Expression of Nestin was related to triple-negative metastases, where brain was the most frequent metastatic site. These findings highlight the clinical utility of FOXA1 and Nestin expression and warrant their inclusion in routine immunohistochemical panels for breast carcinoma.

The Dominant Role of Forkhead Box Proteins in Cancer

Forkhead box (FOX) proteins are multifaceted transcription factors that are significantly implicated in cancer, with various critical roles in biological processes. Herein, we provide an overview of several key members of the FOXA, FOXC, FOXM1, FOXO and FOXP subfamilies. Important pathophysiological processes of FOX transcription factors at multiple levels in a context-dependent manner are discussed. We also specifically summarize some major aspects of FOX transcription factors in association with cancer research such as drug resistance, tumor growth, genomic alterations or drivers of initiation. Finally, we suggest that targeting FOX proteins may be a potential therapeutic strategy to combat cancer.

Conservation of epigenetic regulation by the MLL3/4 tumour suppressor in planarian pluripotent stem cells

Currently, little is known about the evolution of epigenetic regulation in animal stem cells. Here we demonstrate, using the planarian stem cell system to investigate the role of the COMPASS family of MLL3/4 histone methyltransferases that their function as tumor suppressors in mammalian stem cells is conserved over a long evolutionary distance. To investigate the potential conservation of a genome-wide epigenetic regulatory program in animal stem cells, we assess the effects of Mll3/4 loss of function by performing RNA-seq and ChIP-seq on the G2/M planarian stem cell population, part of which contributes to the formation of outgrowths. We find many oncogenes and tumor suppressors among the affected genes that are likely candidates for mediating MLL3/4 tumor suppression function. Our work demonstrates conservation of an important epigenetic regulatory program in animals and highlights the utility of the planarian model system for studying epigenetic regulation.

Aging Mouse Models Reveal Complex Tumor-Microenvironment Interactions in Cancer Progression

Mouse models and genetically engineered mouse models (GEMM) are essential experimental tools for the understanding molecular mechanisms within complex biological systems. GEMM are especially useful for inferencing phenocopy information to genetic human diseases such as breast cancer. Human breast cancer modeling in mice most commonly employs mammary epithelial-specific promoters to investigate gene function(s) and, in particular, putative oncogenes. Models are specifically useful in the mammary epithelial cell in the context of the complete mammary gland environment. Gene targeted knockout mice including conditional targeting to specific mammary cells can reveal developmental defects in mammary organogenesis and demonstrate the importance of putative tumor suppressor genes. Some of these models demonstrate a non-traditional type of tumor suppression which involves interplay between the tumor susceptible cell and its host/environment. These GEMM help to reveal the processes of cancer progression beyond those intrinsic to cancer cells. Furthermore, the, analysis of mouse models requires appropriate consideration of mouse strain, background, and environmental factors. In this review, we compare aging-related factors in mouse models for breast cancer. We introduce databases of GEMM attributes and colony functional variations.

Break Breast Cancer Addiction by CRISPR/Cas9 Genome Editing

Breast cancer is the leading diagnosed cancer for women globally. Evolution of breast cancer in tumorigenesis, metastasis and treatment resistance appears to be driven by the aberrant gene expression and protein degradation encoded by the cancer genomes. The uncontrolled cancer growth relies on these cellular events, thus constituting the cancerous programs and rendering the addiction towards them. These programs are likely the potential anticancer biomarkers for Personalized Medicine of breast cancer. This review intends to delineate the impact of the CRSPR/Cas-mediated genome editing in identification and validation of these anticancer biomarkers. It reviews the progress in three aspects of CRISPR/Cas9-mediated editing of the breast cancer genomes: Somatic genome editing, transcription and protein degradation addictions.

SNAIL1-mediated downregulation of FOXA proteins facilitates the inactivation of transcriptional enhancer elements at key epithelial genes in colorectal cancer cells

Phenotypic conversion of tumor cells through epithelial-mesenchymal transition (EMT) requires massive gene expression changes. How these are brought about is not clear. Here we examined the impact of the EMT master regulator SNAIL1 on the FOXA family of transcription factors which are distinguished by their particular competence to induce chromatin reorganization for the activation of transcriptional enhancer elements. We show that the expression of SNAIL1 and FOXA genes is anticorrelated in transcriptomes of colorectal tumors and cell lines. In cellular EMT models, ectopically expressed Snail1 directly represses FOXA1 and triggers downregulation of all FOXA family members, suggesting that loss of FOXA expression promotes EMT. Indeed, cells with CRISPR/Cas9-induced FOXA-deficiency acquire mesenchymal characteristics. Furthermore, ChIP-seq data analysis of FOXA chromosomal distribution in relation to chromatin structural features which characterize distinct states of transcriptional activity, revealed preferential localization of FOXA factors to transcriptional enhancers at signature genes that distinguish epithelial from mesenchymal colon tumors. To validate the significance of this association, we investigated the impact of FOXA factors on structure and function of enhancers at the CDH1, CDX2 and EPHB3 genes. FOXA-deficiency and expression of dominant negative FOXA2 led to chromatin condensation at these enhancer elements. Site-directed mutagenesis of FOXA binding sites in reporter gene constructs and by genome-editing in situ impaired enhancer activity and completely abolished the active chromatin state of the EPHB3 enhancer. Conversely, expression of FOXA factors in cells with inactive CDX2 and EPHB3 enhancers led to chromatin opening and de novo deposition of the H3K4me1 and H3K27ac marks. These findings establish the pioneer function of FOXA factors at enhancer regions of epithelial genes and demonstrate their essential role in maintaining enhancer structure and function. Thus, by repressing FOXA family members, SNAIL1 targets transcription factors at strategically important positions in gene-regulatory hierarchies, which may facilitate transcriptional reprogramming during EMT.

Cancer patient mortality is overwhelmingly due to distant organ metastases. Epithelial-mesenchymal transition is a process thought to facilitate local invasion and dissemination of cancer cells, thereby promoting metastasis. The conversion of epithelial cells into mesenchymal, fibroblast-like cells requires profound gene expression changes. A few transcription factors like SNAIL1 can initiate these changes, but are unlikely to be solely responsible for all of them. In our study we asked, whether destabilization of epithelial gene expression programs could involve FOXA transcription factors. FOXA factors represent a special subgroup of regulatory proteins, so-called pioneer factors, with unique roles in the activation of transcriptional enhancers which are key regulatory DNA elements that orchestrate spatio-temporal gene expression. In a model of colorectal cancer we found that SNAIL1 represses FOXA factors, and demonstrate that FOXA factors are associated with enhancer elements at epithelial signature genes. Indeed, FOXA factors are sufficient to initiate enhancer activation and necessary to maintain their activity. Our findings indicate that SNAIL1 induces pervasive repression of epithelial genes through a hierarchical scheme of alterations in transcription factor expression which may be applicable to other instances of cell fate changes and transcriptional reprogramming.

Pathobiology of the 129:Stat1 -/- mouse model of human age-related ER-positive breast cancer with an immune infiltrate-excluded phenotype

Background

Stat1 gene-targeted knockout mice (129S6/SvEvTac-Stat1^tm1Rds) develop estrogen receptor-positive (ER⁺), luminal-type mammary carcinomas at an advanced age. There is evidence for both host environment as well as tumor cell-intrinsic mechanisms to initiate tumorigenesis in this model. In this report, we summarize details of the systemic and mammary pathology at preneoplastic and tumor-bearing time points. In addition, we investigate tumor progression in the 129:Stat1^−/− host compared with wild-type 129/SvEv, and we describe the immune cell reaction to the tumors.

Methods

Mice housed and treated according to National Institutes of Health guidelines and Institutional Animal Care and Use Committee-approved methods were evaluated by histopathology, and their tissues were subjected to immunohistochemistry with computer-assisted quantitative image analysis. Tumor cell culture and conditioned media from cell culture were used to perform macrophage (RAW264.7) cell migration assays, including the 129:Stat1^−/−-derived SSM2 cells as well as control Met1 and NDL tumor cells and EpH4 normal cells.

Results

Tumorigenesis in 129:Stat1^−/− originates from a population of FoxA1⁺ large oval pale cells that initially appear and accumulate along the mammary ducts in segments or regions of the gland prior to giving rise to mammary intraepithelial neoplasias. Progression to invasive carcinoma is accompanied by a marked local stromal and immune cell response composed predominantly of T cells and macrophages. In conditioned media experiments, cells derived from 129:Stat1^−/− tumors secrete both chemoattractant and chemoinhibitory factors, with greater attraction in the extracellular vesicular fraction and inhibition in the soluble fraction. The result appears to be recruitment of the immune reaction to the periphery of the tumor, with exclusion of immune cell infiltration into the tumor.

Conclusions

129:Stat1^−/− is a unique model for studying the critical origins and risk reduction strategies in age-related ER⁺ breast cancer. In addition, it can be used in preclinical trials of hormonal and targeted therapies as well as immunotherapies.

Electronic supplementary material

The online version of this article (doi:10.1186/s13058-017-0892-8) contains supplementary material, which is available to authorized users.

The Emerging Roles of Forkhead Box (FOX) Proteins in Osteosarcoma

Osteosarcoma is the most common bone cancer primarily occurring in children and young adults. Over the past few years, the deregulation of a superfamily transcription factors, known as forkhead box (FOX) proteins, has been demonstrated to contribute to the pathogenesis of osteosarcoma. Molecular mechanism studies have demonstrated that FOX family proteins participate in a variety of signaling pathways and that their expression can be regulated by multiple factors. The dysfunction of FOX genes can alter osteosarcoma cell differentiation, metastasis and progression. In this review, we summarized the evidence that FOX genes play direct or indirect roles in the development and progression of osteosarcoma, and evaluated the emerging role of FOX proteins as targets for therapeutic intervention.

Foxa1 gene and protein in developing rat eccrine sweat glands

To investigate the development of eccrine sweat glands and the expression of Foxa1 genes and proteins in the course of development, the footpads from E15.5 to E21.5, P1-P12, P14, P21, P28 and P56 rats were subjected to immunofluorescence staining of FoxA1 and double immunofluorescence staining of K14/α-SMA, FoxA1/K7 and FoxA1/α-SMA, and were processed for Foxa1 gene detection by RT-qPCR. The results showed that rat eccrine sweat gland germs was first observed emerging from the basal layer of epidermis at E19.5, and then elongated downward into the dermis, forming straight ducts by E21.5. Early development of the secretory segments appeared at P1. The Foxa1 gene was not expressed in rat footpads until P2, but from P2 to P5, its expression up-regulated sharply, and thereafter maintained at a high level until adulthood. FoxA1 protein was first observed at P6 in eccrine sweat glands, four days after initial detection of Foxa1 gene transcripts. In skin, FoxA1-positive cells were present exclusively in secretory coils, with 95% being K7-positive secretory cells and 5% being α-SMA-positive myoepithelial cells. We conclude that Foxa1 can be used as a marker of eccrine sweat glands in skin and also as a marker of secretory coils, and Foxa1 is related to the development of secretory coils.