FOXA1 represses the molecular phenotype of basal breast cancer cells

Increased prevalence of hybrid epithelial/mesenchymal state and enhanced phenotypic heterogeneity in basal breast cancer

Intra-tumoral phenotypic heterogeneity promotes tumor relapse and therapeutic resistance and remains an unsolved clinical challenge. Decoding the interconnections among different biological axes of plasticity is crucial to understand the molecular origins of phenotypic heterogeneity. Here, we use multi-modal transcriptomic data-bulk, single-cell, and spatial transcriptomics-from breast cancer cell lines and primary tumor samples, to identify associations between epithelial-mesenchymal transition (EMT) and luminal-basal plasticity-two key processes that enable heterogeneity. We show that luminal breast cancer strongly associates with an epithelial cell state, but basal breast cancer is associated with hybrid epithelial/mesenchymal phenotype(s) and higher phenotypic heterogeneity. Mathematical modeling of core underlying gene regulatory networks representative of the crosstalk between the luminal-basal and epithelial-mesenchymal axes elucidate mechanistic underpinnings of the observed associations from transcriptomic data. Our systems-based approach integrating multi-modal data analysis with mechanism-based modeling offers a predictive framework to characterize intra-tumor heterogeneity and identify interventions to restrict it.

Identification of lineage-specific epigenetic regulators FOXA1 and GRHL2 through chromatin accessibility profiling in breast cancer cell lines

Breast cancer is a heterogeneous disease, and breast cancer cell lines are invaluable for studying this heterogeneity. However, the epigenetic diversity across these cell lines remains poorly understood. In this study, we performed genome-wide chromatin accessibility analysis on 23 breast cancer cell lines, including 2 estrogen receptor (ER)-positive/human epidermal growth factor receptor 2 (HER2)-negative (ER+/HER2-), 3 ER+/HER2+, 3 HER2+, and 15 triple-negative breast cancer (TNBC) lines. These cell lines were classified into three groups based on their chromatin accessibility: the receptor-positive group (Group-P), TNBC basal group (Group-B), and TNBC mesenchymal group (Group-M). Motif enrichment analysis revealed that only Group-P exhibited coenrichment of forkhead box A1 (FOXA1) and grainyhead-like 2 (GRHL2) motifs, whereas Group-B was characterized by the presence of the GRHL2 motif without FOXA1. Notably, Group-M did not show enrichment of either FOXA1 or GRHL2 motifs. Furthermore, gene ontology analysis suggested that group-specific accessible regions were associated with their unique lineage characteristics. To investigate the epigenetic landscape regulatory roles of FOXA1 and GRHL2, we performed knockdown experiments targeting FOXA1 and GRHL2, followed by assay for transposase-accessible chromatin sequencing analysis. The findings revealed that FOXA1 maintains Group-P-specific regions while suppressing Group-B-specific regions in Group-P cells. In contrast, GRHL2 preserves commonly accessible regions shared between Group-P and Group-B in Group-B cells, suggesting that FOXA1 and GRHL2 play a pivotal role in preserving distinct chromatin accessibility patterns for each group. Specifically, FOXA1 distinguishes between receptor-positive and TNBC cell lines, whereas GRHL2 distinguishes between basal-like and mesenchymal subtypes in TNBC lines.

The role of transcription factor FOXA1/C2/M1/O3/P1/Q1 in breast cancer

Breast cancer is a common malignancy with the highest mortality rate among women worldwide. Its incidence is on the rise year after year, accounting for more than one-tenth of new cancers worldwide. Increasing evidence suggests that forkhead box (FOX) transcription factors play an important role in the occurrence and development of breast cancer. However, little is known about the relationship between the expression, prognostic value, function, and immune infiltration of FOX transcription factors in tumor microenvironment. We used bioinformatics to investigate expression and function of FOX factor in breast cancer. Our results revealed the expression levels of FOXA1 and FOXM1 were significantly higher in breast cancer tissues than in normal tissues. The high expression of mRNA in FOXA1 (P < .05), FOXM1 (P < .01), and FOXP1 (P < .05) groups was related to tumor stage. Survival analysis results showed that increased FOXP1 mRNA levels were significantly associated with overall survival (OS), recurrence-free survival (RFS), and distant metastasis-free survival (DMFS) in all patients with breast cancer (P < .05). Patients with the FOXA1 high-expression group had better RFS and DMFS than the low-expression group (P < .05), while patients with FOXM1 high-expression group had worse RFS, OS, and DMFS than the low-expression group (P < .05). Meanwhile, mutation analysis showed that genetic alterations in FOX transcription factors were significantly associated with shorter OS and progression-free survival (P < .05), but not with disease-free survival (P = .710) in patients with breast cancer. FOXP1, FOXA1, and FOXM1 may be used as potential biomarkers to predict the prognosis of patients with breast cancer. Functional enrichment indicated that FOX was mainly involved in cell division, cell senescence, cell cycle, and prolactin signaling pathway. In patients with breast cancer, FOXC2 expression was negatively correlated with the infiltration of B cells and positively correlated with the infiltration of neutrophils and dendritic cells. However, FOXM1 was negatively correlated with the infiltration of CD8 + T cells and macrophages and positively correlated with the infiltration of neutrophils and dendritic cells. These findings provided novel insights into the screening of prognostic biomarkers of the FOX family in breast cancer and laid a foundation for further research on the immune infiltration of the FOX transcription factor family members in tumors.

FOXA1 and FOXA2: the regulatory mechanisms and therapeutic implications in cancer

FOXA1 (Forkhead Box A1) and FOXA2 (Forkhead Box A2) serve as pioneering transcription factors that build gene expression capacity and play a central role in biological processes, including organogenesis and differentiation, glycolipid metabolism, proliferation, migration and invasion, and drug resistance. Notably, FOXA1 and FOXA2 may exert antagonistic, synergistic, or complementary effects in the aforementioned biological processes. This article focuses on the molecular mechanisms and clinical relevance of FOXA1 and FOXA2 in steroid hormone-induced malignancies and highlights potential strategies for targeting FOXA1 and FOXA2 for cancer therapy. Furthermore, the article describes the prospect of targeting upstream regulators of FOXA1/FOXA2 to regulate its expression for cancer therapy because of the drug untargetability of FOXA1/FOXA2.

Assessing Heterogeneity in the N-Telopeptides of Type I Collagen by Mass Spectrometry

Collagen cross-links created by the lysyl oxidase and lysyl hydroxylase families of enzymes are a significant contributing factor to the biomechanical strength and rigidity of tissues, which in turn influence cell signaling and ultimately cell phenotype. In the clinic, the proteolytically liberated N-terminal cross-linked peptide of collagen I (NTX) is used as a biomarker of bone and connective tissue turnover, which is altered in several disease processes. Despite the clinical utility of these collagen breakdown products, the majority of the cross-linked peptide species have not been identified in proteomic datasets. Here we evaluate several parameters for the preparation and identification of these peptides from the collagen I-rich Achilles tendon. Our refined approach involving chemical digestion for protein solubilization coupled with mass spectrometry allows for the identification of the NTX cross-links in a range of modification states. Based on the specificity of the enzymatic cross-linking reaction we utilized follow-up variable modification searches to facilitate identification with a wider range of analytical workflows. We then applied a spectral library approach to identify differences in collagen cross-links in bovine pulmonary hypertension. The presented method offers unique opportunities to understand extracellular matrix remodeling events in development, aging, wound healing, and fibrotic disease that modulate collagen architecture through lysyl-hydroxylase and lysyl-oxidase enzymes.

Subtype Transdifferentiation in Human Cancer: The Power of Tissue Plasticity in Tumor Progression

The classification of tumors into subtypes, characterized by phenotypes determined by specific differentiation pathways, aids diagnosis and directs therapy towards targeted approaches. However, with the advent and explosion of next-generation sequencing, cancer phenotypes are turning out to be far more heterogenous than initially thought, and the classification is continually being updated to include more subtypes. Tumors are indeed highly dynamic, and they can evolve and undergo various changes in their characteristics during disease progression. The picture becomes even more complex when the tumor responds to a therapy. In all these cases, cancer cells acquire the ability to transdifferentiate, changing subtype, and adapt to changing microenvironments. These modifications affect the tumor's growth rate, invasiveness, response to treatment, and overall clinical behavior. Studying tumor subtype transitions is crucial for understanding tumor evolution, predicting disease outcomes, and developing personalized treatment strategies. We discuss this emerging hallmark of cancer and the molecular mechanisms involved at the crossroads between tumor cells and their microenvironment, focusing on four different human cancers in which tissue plasticity causes a subtype switch: breast cancer, prostate cancer, glioblastoma, and pancreatic adenocarcinoma.

Staufen1 Represses the FOXA1-Regulated Transcriptome by Destabilizing FOXA1 mRNA in Colorectal Cancer Cells

Transcription factors play key roles in development and disease by controlling gene expression. Forkhead box A1 (FOXA1), is a pioneer transcription factor essential for mouse development and functions as an oncogene in prostate and breast cancer. In colorectal cancer (CRC), FOXA1 is significantly downregulated and high FOXA1 expression is associated with better prognosis, suggesting potential tumor suppressive functions. We therefore investigated the regulation of FOXA1 expression in CRC, focusing on well-differentiated CRC cells, where FOXA1 is robustly expressed. Genome-wide RNA stability assays identified FOXA1 as an unstable mRNA in CRC cells. We validated FOXA1 mRNA instability in multiple CRC cell lines and in patient-derived CRC organoids, and found that the FOXA1 3'UTR confers instability to the FOXA1 transcript. RNA pulldowns and mass spectrometry identified Staufen1 (STAU1) as a potential regulator of FOXA1 mRNA. Indeed, STAU1 knockdown resulted in increased FOXA1 mRNA and protein expression due to increased FOXA1 mRNA stability. Consistent with these data, RNA-seq following STAU1 knockdown in CRC cells revealed that FOXA1 targets were upregulated upon STAU1 knockdown. Collectively, this study uncovers a molecular mechanism by which FOXA1 is regulated in CRC cells and provides insights into our understanding of the complex mechanisms of gene regulation in cancer.

KBTBD2 promotes proliferation and migration of gastric cancer via activating EGFR signaling pathway

BACKGROUND

To explore the role of Kelch repeat and BTB (POZ) domain containing 2 (KBTBD2) in Gastric cancer(GC) via studying the level of KBTBD2 and its impact on GC cells and mice model.

METHODS

Expression of KBTBD2 in GC was analyzed by analysis of TCGA data, Western blotting and Real-time quantitative polymerasechain reaction (RT-qPCR). The role of KBTBD2 on GC cells proliferation, viability, invasion, migration and apoptosis in vitro were assessed by using western blotting，RT-qPCR，CCK-8, EDU, Colony Formation Assay, Wound healing assay, Transwell, JC-1 mitochondrial membrane potential and flow cytometry assay, respectively. And levels of Bcl-2, BAX, PARP, E-cadherin, Vimentin, N-cadherin, EGFR, SOS1, NROS, BRAF,ERK1/2 and GAPDH were tested by western blotting. Relation of KBTBD2 and epidermal growth factor receptor (EGFR) was predicted by KEGG analysis. KBTBD2 gene GSEA enrichment was analyzed by using R language. Moreover, CCK-8, western blotting, and wound healing assays were used to verify the correlation of KBTBD2 and EGFR pathway. Finally, tumor growth in mice was also investigated. Cells proliferation, migration and apoptosis were detected by Ki67 staining, Tunnel staining and mouse lung metastasis model.

RESULTS

KBTBD2 was highly expressed in GC, and was related to poor prognosis. Moreover, silencing KBTBD2 suppressed GC cell proliferation, migration and invasion, while also inhibited the EMT, but promoted apoptosis. At the same time, KBTBD2 overexpression showed opposite results. In addition, KBTBD2 regulated the EGFR pathway. Further, silencing KBTBD2 inhibited tumor growth, cell proliferation and migration but promoted apoptosis in vivo, and KBTBD2 overexpression showed opposite results.

CONCLUSIONS

KBTBD2 was highly expressed in GC. KBTBD2 promotes the progress of GC by activating EGFR signal pathway. KBTBD2 may thus be a novel target for treating GC.

Identification of aberrant luminal progenitors and mTORC1 as a potential breast cancer prevention target in BRCA2 mutation carriers

Inheritance of a BRCA2 pathogenic variant conveys a substantial life-time risk of breast cancer. Identification of the cell(s)-of-origin of BRCA2-mutant breast cancer and targetable perturbations that contribute to transformation remains an unmet need for these individuals who frequently undergo prophylactic mastectomy. Using preneoplastic specimens from age-matched, premenopausal females, here we show broad dysregulation across the luminal compartment in BRCA2mut/+ tissue, including expansion of aberrant ERBB3lo luminal progenitor and mature cells, and the presence of atypical oestrogen receptor (ER)-positive lesions. Transcriptional profiling and functional assays revealed perturbed proteostasis and translation in ERBB3lo progenitors in BRCA2mut/+ breast tissue, independent of ageing. Similar molecular perturbations marked tumours bearing BRCA2-truncating mutations. ERBB3lo progenitors could generate both ER+ and ER- cells, potentially serving as cells-of-origin for ER-positive or triple-negative cancers. Short-term treatment with an mTORC1 inhibitor substantially curtailed tumorigenesis in a preclinical model of BRCA2-deficient breast cancer, thus uncovering a potential prevention strategy for BRCA2 mutation carriers.

Epigenome editing in cancer: Advances and challenges for potential therapeutic options

Cancers are diseases caused by genetic and non-genetic environmental factors. Epigenetic alterations, some attributed to non-genetic factors, can lead to cancer development. Epigenetic changes can occur in tumor suppressors or oncogenes, or they may contribute to global cell state changes, making cells abnormal. Recent advances in gene editing technology show potential for cancer treatment. Herein, we will discuss our current knowledge of epigenetic alterations occurring in cancer and epigenetic editing technologies that can be applied to developing therapeutic options.

INO80 function is required for mouse mammary gland development, but mutation alone may be insufficient for breast cancer

The aberrant function of ATP-dependent chromatin remodeler INO80 has been implicated in multiple types of cancers by altering chromatin architecture and gene expression; however, the underlying mechanism of the functional involvement of INO80 mutation in cancer etiology, especially in breast cancer, remains unclear. In the present study, we have performed a weighted gene co-expression network analysis (WCGNA) to investigate links between INO80 expression and breast cancer sub-classification and progression. Our analysis revealed that INO80 repression is associated with differential responsiveness of estrogen receptors (ERs) depending upon breast cancer subtype, ER networks, and increased risk of breast carcinogenesis. To determine whether INO80 loss induces breast tumors, a conditional INO80-knockout (INO80 cKO) mouse model was generated using the Cre-loxP system. Phenotypic characterization revealed that INO80 cKO led to reduced branching and length of the mammary ducts at all stages. However, the INO80 cKO mouse model had unaltered lumen morphology and failed to spontaneously induce tumorigenesis in mammary gland tissue. Therefore, our study suggests that the aberrant function of INO80 is potentially associated with breast cancer by modulating gene expression. INO80 mutation alone is insufficient for breast tumorigenesis.

Multi-modal transcriptomic analysis unravels enrichment of hybrid epithelial/mesenchymal state and enhanced phenotypic heterogeneity in basal breast cancer

Intra-tumoral phenotypic heterogeneity promotes tumor relapse and therapeutic resistance and remains an unsolved clinical challenge. It manifests along multiple phenotypic axes and decoding the interconnections among these different axes is crucial to understand its molecular origins and to develop novel therapeutic strategies to control it. Here, we use multi-modal transcriptomic data analysis - bulk, single-cell and spatial transcriptomics - from breast cancer cell lines and primary tumor samples, to identify associations between epithelial-mesenchymal transition (EMT) and luminal-basal plasticity - two key processes that enable heterogeneity. We show that luminal breast cancer strongly associates with an epithelial cell state, but basal breast cancer is associated with hybrid epithelial/mesenchymal phenotype(s) and higher phenotypic heterogeneity. These patterns were inherent in methylation profiles, suggesting an epigenetic crosstalk between EMT and lineage plasticity in breast cancer. Mathematical modelling of core underlying gene regulatory networks representative of the crosstalk between the luminal-basal and epithelial-mesenchymal axes recapitulate and thus elucidate mechanistic underpinnings of the observed associations from transcriptomic data. Our systems-based approach integrating multi-modal data analysis with mechanism-based modeling offers a predictive framework to characterize intra-tumor heterogeneity and to identify possible interventions to restrict it.

A contrast set mining based approach for cancer subtype analysis

The task of detecting common and unique characteristics among different cancer subtypes is an important focus of research that aims to improve personalized therapies. Unlike current approaches mainly based on predictive techniques, our study aims to improve the knowledge about the molecular mechanisms that descriptively led to cancer, thus not requiring previous knowledge to be validated. Here, we propose an approach based on contrast set mining to capture high-order relationships in cancer transcriptomic data. In this way, we were able to extract valuable insights from several cancer subtypes in the form of highly specific genetic relationships related to functional pathways affected by the disease. To this end, we have divided several cancer gene expression databases by the subtype associated with each sample to detect which gene groups are related to each cancer subtype. To demonstrate the potential and usefulness of the proposed approach we have extensively analysed RNA-Seq gene expression data from breast, kidney, and colon cancer subtypes. The possible role of the obtained genetic relationships was further evaluated through extensive literature research, while its prognosis was assessed via survival analysis, finding gene expression patterns related to survival in various cancer subtypes. Some gene associations were described in the literature as potential cancer biomarkers while other results have been not described yet and could be a starting point for future research.

TLE3 Sustains Luminal Breast Cancer Lineage Fidelity to Suppress Metastasis

Breast cancer subtypes and their phenotypes parallel different stages of the mammary epithelial cell developmental hierarchy. Discovering mechanisms that control lineage identity could provide novel avenues for mitigating disease progression. Here we report that the transcriptional corepressor TLE3 is a guardian of luminal cell fate in breast cancer and operates independently of the estrogen receptor. In luminal breast cancer, TLE3 actively repressed the gene-expression signature associated with highly aggressive basal-like breast cancers (BLBC). Moreover, maintenance of the luminal lineage depended on the appropriate localization of TLE3 to its transcriptional targets, a process mediated by interactions with FOXA1. By repressing genes that drive BLBC phenotypes, including SOX9 and TGFβ2, TLE3 prevented the acquisition of a hybrid epithelial-mesenchymal state and reduced metastatic capacity and aggressive cellular behaviors. These results establish TLE3 as an essential transcriptional repressor that sustains the more differentiated and less metastatic nature of luminal breast cancers. Approaches to induce TLE3 expression could promote the acquisition of less aggressive, more treatable disease states to extend patient survival.

SIGNIFICANCE

Transcriptional corepressor TLE3 actively suppresses SOX9 and TGFβ transcriptional programs to sustain the luminal lineage identity of breast cancer cells and to inhibit metastatic progression.

Epigenomic mapping identifies an enhancer repertoire that regulates cell identity in bladder cancer through distinct transcription factor networks

Muscle-invasive bladder cancer (BLCA) is an aggressive disease. Consensus BLCA transcriptomic subtypes have been proposed, with two major Luminal and Basal subgroups, presenting distinct molecular and clinical characteristics. However, how these distinct subtypes are regulated remains unclear. We hypothesized that epigenetic activation of distinct super-enhancers could drive the transcriptional programs of BLCA subtypes. Through integrated RNA-sequencing and epigenomic profiling of histone marks in primary tumours, cancer cell lines, and normal human urothelia, we established the first integrated epigenetic map of BLCA and demonstrated the link between subtype and epigenetic control. We identified the repertoire of activated super-enhancers and highlighted Basal, Luminal and Normal-associated SEs. We revealed super-enhancer-regulated networks of candidate master transcription factors for Luminal and Basal subgroups including FOXA1 and ZBED2, respectively. FOXA1 CRISPR-Cas9 mutation triggered a shift from Luminal to Basal phenotype, confirming its role in Luminal identity regulation and induced ZBED2 overexpression. In parallel, we showed that both FOXA1 and ZBED2 play concordant roles in preventing inflammatory response in cancer cells through STAT2 inhibition. Our study furthers the understanding of epigenetic regulation of muscle-invasive BLCA and identifies a co-regulated network of super-enhancers and associated transcription factors providing potential targets for the treatment of this aggressive disease.

MicroRNAs within the Basal-like signature of Quadruple Negative Breast Cancer impact overall survival in African Americans

We previously found that QNBC tumors are more frequent in African Americans compared to TNBC tumors. To characterize this subtype further, we sought to determine the miRNA-mRNA profile in QNBC patients based on race. Both miRNA and mRNA expression data were analyzed from TCGA and validated using datasets from the METABRIC, TCGA proteomic, and survival analysis by KMPLOT. miRNA-mRNAs which include FOXA1 and MYC (mir-17/20a targets); GATA3 and CCNG2 (mir-135b targets); CDKN2A, CDK6, and B7-H3 (mir-29c targets); and RUNX3, KLF5, IL1-β, and CTNNB1 (mir-375 targets) were correlated with basal-like and immune subtypes in QNBC patients and associated with a worse survival. Thus, QNBC tumors have an altered gene signature implicated in racial disparity and poor survival.

Breast cancer plasticity is restricted by a LATS1-NCOR1 repressive axis

Breast cancer, the most frequent cancer in women, is generally classified into several distinct histological and molecular subtypes. However, single-cell technologies have revealed remarkable cellular and functional heterogeneity across subtypes and even within individual breast tumors. Much of this heterogeneity is attributable to dynamic alterations in the epigenetic landscape of the cancer cells, which promote phenotypic plasticity. Such plasticity, including transition from luminal to basal-like cell identity, can promote disease aggressiveness. We now report that the tumor suppressor LATS1, whose expression is often downregulated in human breast cancer, helps maintain luminal breast cancer cell identity by reducing the chromatin accessibility of genes that are characteristic of a "basal-like" state, preventing their spurious activation. This is achieved via interaction of LATS1 with the NCOR1 nuclear corepressor and recruitment of HDAC1, driving histone H3K27 deacetylation near NCOR1-repressed "basal-like" genes. Consequently, decreased expression of LATS1 elevates the expression of such genes and facilitates slippage towards a more basal-like phenotypic identity. We propose that by enforcing rigorous silencing of repressed genes, the LATS1-NCOR1 axis maintains luminal cell identity and restricts breast cancer progression.

The POLR3G Subunit of Human RNA Polymerase III Regulates Tumorigenesis and Metastasis in Triple-Negative Breast Cancer

RNA polymerase (Pol) III transcribes short untranslated RNAs that contribute to the regulation of gene expression. Two isoforms of human Pol III have been described that differ by the presence of the POLR3G/RPC32α or POLR3GL/RPC32β subunits. POLR3G was found to be expressed in embryonic stem cells and at least a subset of transformed cells, whereas POLR3GL shows a ubiquitous expression pattern. Here, we demonstrate that POLR3G is specifically overexpressed in clinical samples of triple-negative breast cancer (TNBC) but not in other molecular subtypes of breast cancer. POLR3G KO in the MDA-MB231 TNBC cell line dramatically reduces anchorage-independent growth and invasive capabilities in vitro. In addition, the POLR3G KO impairs tumor growth and metastasis formation of orthotopic xenografts in mice. Moreover, KO of POLR3G induces expression of the pioneer transcription factor FOXA1 and androgen receptor. In contrast, the POLR3G KO neither alters proliferation nor the expression of epithelial-mesenchymal transition marker genes. These data demonstrate that POLR3G expression is required for TNBC tumor growth, invasiveness and dissemination and that its deletion affects triple-negative breast cancer-specific gene expression.

Enhancer RNA (eRNA) in Human Diseases

Enhancer RNAs (eRNAs), a class of non-coding RNAs (ncRNAs) transcribed from enhancer regions, serve as a type of critical regulatory element in gene expression. There is increasing evidence demonstrating that the aberrant expression of eRNAs can be broadly detected in various human diseases. Some studies also revealed the potential clinical utility of eRNAs in these diseases. In this review, we summarized the recent studies regarding the pathological mechanisms of eRNAs as well as their potential utility across human diseases, including cancers, neurodegenerative disorders, cardiovascular diseases and metabolic diseases. It could help us to understand how eRNAs are engaged in the processes of diseases and to obtain better insight of eRNAs in diagnosis, prognosis or therapy. The studies we reviewed here indicate the enormous therapeutic potency of eRNAs across human diseases.

Loss of SNAI1 induces cellular plasticity in invasive triple-negative breast cancer cells

The transcription factor SNAI1 mediates epithelial-mesenchymal transition, fibroblast activation and controls inter-tissue migration. High SNAI1 expression characterizes metastatic triple-negative breast carcinomas, and its knockout by CRISPR/Cas9 uncovered an epithelio-mesenchymal phenotype accompanied by reduced signaling by the cytokine TGFβ. The SNAI1 knockout cells exhibited plasticity in differentiation, drifting towards the luminal phenotype, gained stemness potential and could differentiate into acinar mammospheres in 3D culture. Loss of SNAI1 de-repressed the transcription factor FOXA1, a pioneering factor of mammary luminal progenitors. FOXA1 induced a specific gene program, including the androgen receptor (AR). Inhibiting AR via a specific antagonist regenerated the basal phenotype and blocked acinar differentiation. Thus, loss of SNAI1 in the context of triple-negative breast carcinoma cells promotes an intermediary luminal progenitor phenotype that gains differentiation plasticity based on the dual transcriptional action of FOXA1 and AR. This function of SNAI1 provides means to separate cell invasiveness from progenitor cell de-differentiation as independent cellular programs.

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.

Increased FOXA1 levels induce apoptosis and inhibit proliferation in FOXA1-low expressing basal breast cancer cells

The transcription factor FOXA1, which is a member of the forkhead class of DNA-binding proteins, interacts with Estrogen Receptor (ER) to mediate breast cancer progression. However, its role in basal breast cancer cells remains unclear. Although the overall levels of FOXA1 are decreased in the basal subtype of clinical TCGA breast cancer samples, the high levels of FOXA1 improve the survival of the patients from this subtype. This clinical phenomenon is consistent with that of FOXA1 stimulating apoptosis in FOXA1-low expressing basal breast cancer cells, such as MDA-MB-231, and MDA-MB-468 cells. In this study, we have constructed an inducible expression system of FOXA1 and demonstrated the induced expression of FOXA1 resulting in apoptosis and cell cycle arrest in MDA-MB-231 cells, as confirmed by transcriptomic analysis and in vivo tumor-grafted models. Furthermore, the low levels of Estrogen Receptor-1 (ESR1) are critical for FOXA1 in terms of its repressive roles in the cells, as evidenced by clinical data analysis indicating that the high levels of FOXA1 improve the survival of ESR1^Low patients, but worsen the survival of ESR1^High patients of breast cancer. When introduced into MDA-MB-231 cells, ESR1 counteracts the tumor suppressor roles of FOXA1 by altering the FOXA1-regulated gene transcription and the two proteins together maintain the tumor progression in vivo. Our cumulative results suggest that FOXA1 suppresses the basal breast cancer cells with FOXA1-low expressing status independent of ESR1 by inducing apoptosis and inhibiting cell proliferation, thereby implicating its potential therapeutic role in this group of breast cancer.

Gene expression signatures of individual ductal carcinoma in situ lesions identify processes and biomarkers associated with progression towards invasive ductal carcinoma

Ductal carcinoma in situ (DCIS) is considered a non-invasive precursor to breast cancer, and although associated with an increased risk of developing invasive disease, many women with DCIS will never progress beyond their in situ diagnosis. The path from normal duct to invasive ductal carcinoma (IDC) is not well understood, and efforts to do so are hampered by the substantial heterogeneity that exists between patients, and even within patients. Here we show gene expression analysis from > 2,000 individually micro-dissected ductal lesions representing 145 patients. Combining all samples into one continuous trajectory we show there is a progressive loss in basal layer integrity heading towards IDC, coupled with two epithelial to mesenchymal transitions, one early and a second coinciding with the convergence of DCIS and IDC expression profiles. We identify early processes and potential biomarkers, including CAMK2N1, MNX1, ADCY5, HOXC11 and ANKRD22, whose reduced expression is associated with the progression of DCIS to invasive breast cancer.

Progression from ductal carcinoma in situ (DCIS) to invasive ductal carcinoma (IDC) remains poorly understood. Here, the authors analyse over 2700 micro-dissected samples using transcriptomics to identify genes that characterise different stages of DCIS to IDC progression, and identify IDC-associated markers within early-stage lesions.

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.

A Transcriptional Link between HER2, JAM-A and FOXA1 in Breast Cancer

Overexpression of the human epidermal growth factor receptor-2 (HER2) is associated with aggressive disease in breast and certain other cancers. At a cellular level, the adhesion protein Junctional Adhesion Molecule-A (JAM-A) has been reported to regulate the expression of HER3 via a transcriptional pathway involving FOXA1. Since FOXA1 is also a suggested transcription factor for HER2, this study set out to determine if JAM-A regulates HER2 expression via a similar mechanism. An integrated tripartite approach was taken, involving cellular expression studies after targeted disruption of individual players in the putative pathway, in silico identification of relevant HER2 promoter regions and, finally, interrogation of cancer patient survival databases to deconstruct functionally important links between HER2, JAM-A and FOXA1 gene expression. The outcome of these investigations revealed a unidirectional pathway in which JAM-A expression transcriptionally regulates that of HER2 by influencing the binding of FOXA1 to a specific site in the HER2 gene promoter. Moreover, a correlation between JAM-A and HER2 gene expression was identified in 75% of a sample of 40 cancer types from The Cancer Genome Atlas, and coincident high mean mRNA expression of JAM-A, HER2 and FOXA1 was associated with poorer survival outcomes in HER2-positive (but not HER2-negative) patients with either breast or gastric tumors. These investigations provide the first evidence of a transcriptional pathway linking JAM-A, HER2 and FOXA1 in cancer settings, and support potential future pharmacological targeting of JAM-A as an upstream regulator of HER2.

Subtype-Independent ANP32E Reduction During Breast Cancer Progression in Accordance with Chromatin Relaxation

BACKGROUND

Chromatin state provides a clear decipherable blueprint for maintenance of transcriptional patterns, exemplifying a mitotically stable form of cellular programming in dividing cells. In this regard, genomic studies of chromatin states within cancerous tissues have the potential to uncover novel aspects of tumor biology and unique mechanisms associated with disease phenotypes and outcomes. The degree to which chromatin state differences occur in accordance with breast cancer features has not been established.

METHODS

We applied a series of unsupervised computational methods to identify chromatin and molecular differences associated with discrete physiologies across human breast cancer tumors.

RESULTS

Chromatin patterns alone are capable of stratifying tumors in association with cancer subtype and disease progression. Major differences occur at DNA motifs for the transcription factor FOXA1, in hormone receptor-positive tumors, and motifs for SOX9 in Basal-like tumors. We find that one potential driver of this effect, the histone chaperone ANP32E, is inversely correlated with tumor progression and relaxation of chromatin at FOXA1 binding sites. Tumors with high levels of ANP32E exhibit an immune response and proliferative gene expression signature, whereas tumors with low ANP32E levels appear programmed for differentiation.

CONCLUSIONS

Our results indicate that ANP32E may function through chromatin state regulation to control breast cancer differentiation and tumor plasticity. This study sets a precedent for future computational studies of chromatin changes in carcinogenesis.

Positive Regulation of Estrogen Receptor Alpha in Breast Tumorigenesis

Estrogen receptor alpha (ERα, NR3A1) contributes through its expression in different tissues to a spectrum of physiological processes, including reproductive system development and physiology, bone mass maintenance, as well as cardiovascular and central nervous system functions. It is also one of the main drivers of tumorigenesis in breast and uterine cancer and can be targeted by several types of hormonal therapies. ERα is expressed in a subset of luminal cells corresponding to less than 10% of normal mammary epithelial cells and in over 70% of breast tumors (ER+ tumors), but the basis for its selective expression in normal or cancer tissues remains incompletely understood. The mapping of alternative promoters and regulatory elements has delineated the complex genomic structure of the ESR1 gene and shed light on the mechanistic basis for the tissue-specific regulation of ESR1 expression. However, much remains to be uncovered to better understand how ESR1 expression is regulated in breast cancer. This review recapitulates the current body of knowledge on the structure of the ESR1 gene and the complex mechanisms controlling its expression in breast tumors. In particular, we discuss the impact of genetic alterations, chromatin modifications, and enhanced expression of other luminal transcription regulators on ESR1 expression in tumor cells.

Analysis of multiple gene co-expression networks to discover interactions favoring CFTR biogenesis and ΔF508-CFTR rescue

BACKGROUND

We previously reported that expression of a miR-138 mimic or knockdown of SIN3A in primary cultures of cystic fibrosis (CF) airway epithelia increased ΔF508-CFTR mRNA and protein levels, and partially restored CFTR-dependent chloride transport. Global mRNA transcript profiling in ΔF508-CFBE cells treated with miR-138 mimic or SIN3A siRNA identified two genes, SYVN1 and NEDD8, whose inhibition significantly increased ΔF508-CFTR trafficking, maturation, and function. Little is known regarding the dynamic changes in the CFTR gene network during such rescue events. We hypothesized that analysis of condition-specific gene networks from transcriptomic data characterizing ΔF508-CFTR rescue could help identify dynamic gene modules associated with CFTR biogenesis.

METHODS

We applied a computational method, termed M-module, to analyze multiple gene networks, each of which exhibited differential activity compared to a baseline condition. In doing so, we identified both unique and shared gene pathways across multiple differential networks. To construct differential networks, gene expression data from CFBE cells were divided into three groups: (1) siRNA inhibition of NEDD8 and SYVN1; (2) miR-138 mimic and SIN3A siRNA; and (3) temperature (27 °C for 24 h, 40 °C for 24 h, and 27 °C for 24 h followed by 40 °C for 24 h).

RESULTS

Interrogation of individual networks (e.g., NEDD8/SYVN1 network), combinations of two networks (e.g., NEDD8/SYVN1 + temperature networks), and all three networks yielded sets of 1-modules, 2-modules, and 3-modules, respectively. Gene ontology analysis revealed significant enrichment of dynamic modules in pathways including translation, protein metabolic/catabolic processes, protein complex assembly, and endocytosis. Candidate CFTR effectors identified in the analysis included CHURC1, GZF1, and RPL15, and siRNA-mediated knockdown of these genes partially restored CFTR-dependent transepithelial chloride current to ΔF508-CFBE cells.

CONCLUSIONS

The ability of the M-module to identify dynamic modules involved in ΔF508 rescue provides a novel approach for studying CFTR biogenesis and identifying candidate suppressors of ΔF508.

FOXA1: A Pioneer of Nuclear Receptor Action in Breast Cancer

The pioneering function of FOXA1 establishes estrogen-responsive transcriptomes in luminal breast cancer. Dysregulated FOXA1 chromatin occupancy through focal amplification, mutation, or cofactor recruitment modulates estrogen receptor (ER) transcriptional programs and drives endocrine-resistant disease. However, ER is not the sole nuclear receptor (NR) expressed in breast cancers, nor is it the only NR for which FOXA1 serves as a licensing factor. Receptors for androgens, glucocorticoids, and progesterone are also found in the majority of breast cancers, and their functions are also impacted by FOXA1. These NRs interface with ER transcriptional programs and, depending on their activation level, can reprogram FOXA1-ER cistromes. Thus, NR interplay contributes to endocrine therapy response and resistance and may provide a vulnerability for future therapeutic benefit in patients. Herein, we review what is known regarding FOXA1 regulation of NR function in breast cancer in the context of cell identity, endocrine resistance, and NR crosstalk in breast cancer progression and treatment.

Gene expression analysis of invasive breast carcinoma yields differential patterns in luminal subtypes of breast cancer

Breast cancer is a heterogeneous disease, and new biomarkers are needed for more accurate classification and prediction of prognosis. The goal of this study is to assess the expression of breast cancer classification genes, to identify new molecular signatures in different intrinsic subtypes of breast cancer and to correlate their expression with different clinical variables. The study included 84 female patients newly diagnosed with non-metastatic breast cancer at the outpatient clinic at the National Cancer Institute, Cairo University, Egypt. Detection of 17 breast cancer classification genes was done using RT-PCR in tumor and normal tissues. Estrogen receptor (ER), progesterone receptor (PR), HER2, and Ki67 expression were assessed using IHC assay for intrinsic subtyping. Combined expression of FOXA1 and GATA3 was statistically higher in luminal subtypes in comparison to non-luminal subtypes. In Luminal A subtype; GRB7, EGFR, PTGS2, ID1, and KRT5 were significantly downregulated. FOXA1 and GATA3 were significantly upregulated in luminal B subtype, where EGFR and PTGS2 were significantly downregulated. While ESR1, EGFR, KRT5 and PTGS2 showed significantly low expression in tumor tissue in Her2 enriched subtype, TFF3 was significantly downregulated in triple negative subtype. GATA3 and FOXA1 expression exhibited significant correlation with tumor grade. Furthermore, GATA3, FOXA1, ESR1, and ID1 were also correlated significantly with staging of the tumor. Combined expression of ESR1, FOXA1 and GATA3 represents a molecular signature of luminal subtypes. Long term follow-up is needed to investigate the prognostic effect of breast cancer classification genes found in this study.

lncRNA NEAT1 promotes the Taxol resistance of breast cancer via sponging the miR-23a-3p-FOXA1 axis

Breast cancer is the most prevalent malignancy among women worldwide. Paclitaxel (Taxol) is a widely applied chemotherapeutic agent against breast cancer. Although Taxol therapy has achieved improvements recently, development of chemoresistance of breast cancer patients is a major obstacle, leading to therapeutic failure. Long non-coding RNAs (lncRNAs) play pivotal roles in tumorigenesis and progresses of breast cancer. However, the biological roles and molecular targets of lncRNA NEAT1 in Taxol-resistant breast cancer remain unclear. Here, we report that NEAT1 is significantly upregulated in breast tumors and cell lines. In addition, silencing NEAT1 effectively sensitizes breast cancer cells to Taxol. Bioinformatical analysis and luciferase assay demonstrated that miR-23a-3p could be sponged and downregulated by NEAT1. We demonstrated that miR-23a-3p was downregulated and functioned as a tumor suppressor in breast cancer. Furthermore, in the established Taxol-resistant MDA-MB-231 breast cancer cell line, we detected significantly increased NEAT1 expression and downregulated miR-23a-3p expression. Importantly, FOXA1 was identified and validated as a direct target of miR-23a-3p in breast cancer cells. Rescue experiments demonstrated that the restoration of miR-23a-3p in NEAT1-overexpressing Taxol-resistant breast cancer cells successfully overcame the NEAT1-promoted Taxol resistance. Taken together, our results revealed the clinical roles and molecular mechanisms for the NEAT1-mediated chemoresistance, providing new insights into the development of non-coding RNA-based therapeutic strategies for enhancing the anti-cancer effects of traditional chemotherapeutic drugs.

AGR2 and AGR3 play an important role in the clinical characterization and prognosis of basal like breast cancer

INTRODUCTION

Basal-like Breast Cancer (BLBC) represents an important molecular subtype of breast cancer characterized by an aggressive behavior, molecular pathology poorly understood and a limited treatment.

OBJECTIVE

We aim to search for molecular differences between non-BLBC and BLBC tumors in order to propose possible diagnostic and prognostic biomarkers using databases. Metodology: Microarray processed data were downloaded from GEO database considering non-BLBC and BLBC. Enrichment analysis was evaluated using GO consortium and Ingenuity, protein-protein interaction, gene Ontology and co-expression analysis using STRING. Gene expression data was extracted using TCGA, METABRIC and Breast Cancer Gene-Expression Miner v4.2 databases. The Survival was evaluated using The Kaplan-Meier plotter.

RESULTS

Were identified 58 upregulated and 58 downregulated genes enriched in signaling pathways like PDGF, Angiogenesis, Integrin and WNT. AGR2 and AGR3 expression were reduced in BLBC in relation to non-BLBC tumors, patients aged ≤51 years, and with negativity of ER, PR and HER-2 and nodal status. Low expression of AGR2 and AGR3 were associated with worse OS and RFS for all breast cancer cases. But according to the molecular stratification, low AGR2 conferred worst OS in luminal A, worst RFS in BLBC and good OS and RFS in luminal B. High AGR3 conferred worse OS and RFS in BLBC, but low AGR3 attributed worse OS in luminal A.

CONCLUSION

AGR2 and AGR3 expression were able to differentiate non-BLBC from BLBC. Downregulation of AGR2 and AGR3 was associated with BLBC clinical phenotype. Furthermore, both genes behave different when considering prognosis and molecular stratification.

Transcriptional regulation of human DUSP4 gene by cancer-related transcription factors

Dual specificity phosphatase 4 (DUSP4), a member of the dual specificity phosphatase family, is responsible for the dephosphorylation and inactivation of ERK, JNK and p38, which are mitogen-activated protein kinases involved in cell proliferation, differentiation and apoptosis, but also in inflammation processes. Given its importance for cellular signalling, DUSP4 is subjected to a tight regulation and there is growing evidence that its expression is dysregulated in several tumours. However, the mechanisms underlying DUSP4 transcriptional regulation remain poorly understood. Here, we analysed the regulation of the human DUSP4 promoters 1 and 2, located upstream of exons 1 and 2, respectively, by the cancer-related transcription factors (TFs) STAT3, FOXA1, CTCF and YY1. The presence of binding sites for these TFs was predicted in both promoters through the in silico analysis of DUSP4, and their functionality was assessed through luciferase activity assays. Regulatory activity of the TFs tested was found to be promoter-specific. While CTCF stimulated the activity of promoter 2 that controls the transcription of variants 2 and X1, STAT3 stimulated the activity of promoter 1 that controls the transcription of variant 1. YY1 positively regulated both promoters, although to different extents. Through site-directed mutagenesis, the functionality of YY1 binding sites present in promoter 2 was confirmed. This study provides novel insights into the transcriptional regulation of DUSP4, contributing to a better comprehension of the mechanisms of its dysregulation observed in several types of cancer.

Deletion of Foxa1 in the mouse mammary gland results in abnormal accumulation of luminal progenitor cells: a link between reproductive factors and ER-/TNBC breast cancer?

In humans, parity without breastfeeding increases risk of estrogen receptor-negative (ER-) breast cancer and is associated with hypermethylation of FOXA1, a pioneer factor regulating lineage commitment of mammary gland luminal progenitor cells. We postulate that pregnancy-associated repression of FOXA1 results in the accumulation of aberrant, differentiation-arrested luminal progenitor cells which, following additional genetic and epigenetic insults, may give rise to ER- tumors. Consistent with this hypothesis, we show that deletion of Foxa1 in the mouse mammary gland results in a two-fold increase in the proportion of luminal progenitor cells and a reduction in mammary gland epithelial cells that stain positive for ER. These results provide compelling support for the notion that reduced Foxa1 expression is sufficient to alter mammary gland luminal cell fate determination in vivo, which could be a mechanism linking parity with ER- breast cancer.

FOXA1 is a determinant of drug resistance in breast cancer cells

PURPOSE

Breast cancer is one of the most commonly diagnosed cancers in women. Five subtypes of breast cancer differ in their genetic expression profiles and carry different prognostic values, with no treatments available for some types, such as triple-negative, due to the absence of genetic signatures that could otherwise be targeted by molecular therapies. Although endocrine treatments are largely successful for estrogen receptor (ER)-positive cancers, a significant proportion of patients with metastatic tumors fail to respond and acquire resistance to therapy. FOXA1 overexpression mediates endocrine therapy resistance in ER-positive breast cancer, although the regulation of chemotherapy response by FOXA1 has not been addressed previously. FOXA1, together with EP300 and RUNX1, regulates the expression of E-cadherin, and is expressed in luminal, but absent in triple-negative and basal-like breast cancers. We have previously determined that EP300 regulates drug resistance and tumor initiation capabilities in breast cancer cells.

METHODS

Here we describe the generation of breast cancer cell models in which FOXA1 expression has been modulated either by expression of hairpins targeting FOXA1 mRNA or overexpression plasmids.

RESULTS

Upon FOXA1 knockdown in luminal MCF-7 and T47D cells, we found an increase in doxorubicin and paclitaxel sensitivity as well as a decrease in anchorage independence. Conversely, upregulation of FOXA1 in basal-like MDA-MB-231 cells led to an increase in drug resistance and anchorage independence.

CONCLUSION

Together, these data suggest that FOXA1 plays a role in making tumors more aggressive.

FOXQ1 is Differentially Expressed Across Breast Cancer Subtypes with Low Expression Associated with Poor Overall Survival

Purpose

Forkhead box Q1 (FOXQ1) has been shown to contribute to the development and progression of cancers, including ovarian and breast cancer (BC). However, research exploring FOXQ1 expression, copy number variation (CNV), and prognostic value across different BC subtypes is limited. Our purpose was to evaluate FOXQ1 mRNA expression, CNV, and prognostic value across BC subtypes.

Materials and Methods

We determined FOXQ1 expression and CNV in BC patient tumors using RT-qPCR and qPCR, respectively. We also analyzed FOXQ1 expression and CNV in BC cell lines in the CCLE database using K-means clustering. The prognostic value of FOXQ1 expression in the TCGA-BRCA database was assessed using univariate and multivariate Cox's regression analysis as well as using the online tools OncoLnc, GEPIA, and UALCAN.

Results

Our analyses reveal that FOXQ1 mRNA is differentially expressed between different subtypes of BC and is significantly decreased in luminal BC and HER2 patients when compared to normal breast tissue samples. Furthermore, analysis of BC cell lines showed that FOXQ1 mRNA expression was independent of CNV. Moreover, patients with low FOXQ1 mRNA expression had significantly poorer overall survival compared to those with high FOXQ1 mRNA expression. Finally, low FOXQ1 expression had a critical impact on the prognostic values of BC patients and was an independent predictor of overall survival when it was adjusted for BC subtypes and to two other FOX genes, FOXF2 and FOXM1.

Conclusion

Our study reveals for the first time that FOXQ1 is differentially expressed across BC subtypes and that low expression of FOXQ1 is indicative of poor prognosis in patients with BC.

Human Epidermal Growth Factor Receptor-3 Expression Is Regulated at Transcriptional Level in Breast Cancer Settings by Junctional Adhesion Molecule-A via a Pathway Involving Beta-Catenin and FOXA1

Simple Summary

Signaling from the human epidermal growth factor receptor (HER) family of proteins increases in many cancers, including breast. HER2-high breast cancers are successfully treated with anti-HER2 therapies, but these drugs are limited by the fact that patients frequently develop resistance to them. One common mechanism by which resistance develops is when tumors acquire high levels of a family member called HER3. We had previously shown that a protein called JAM-A regulates the level of HER2 in breast cancer cells, and is associated with the development of resistance to HER2-targeted therapies. In this study we show for the first time that JAM-A levels also regulate those of HER3. Using breast cancer cell and tissue models and culminating in patient tissue material, we provide evidence that JAM-A regulates HER3 expression via a pathway involving the transcription factors β-catenin and FOXA1. We suggest that JAM-A merits future investigation as a novel drug target for its potential to reduce HER3 tumorigenic signaling and to offset the development of resistance to HER2-targeted therapies.

Abstract

The success of breast cancer therapies targeting the human epidermal growth factor receptor-2 (HER2) is limited by the development of drug resistance by mechanisms including upregulation of HER3. Having reported that HER2 expression and resistance to HER2-targeted therapies can be regulated by Junctional Adhesion Molecule-A (JAM-A), this study investigated if JAM-A regulates HER3 expression. Expressional alteration of JAM-A in breast cancer cells was used to test expressional effects on HER3 and its effectors, alongside associated functional behaviors, in vitro and semi-in vivo. HER3 transcription factors were identified and tested for regulation by JAM-A. Finally a patient tissue microarray was used to interrogate connections between putative pathway components connecting JAM-A and HER3. This study reveals for the first time that HER3 and its effectors are regulated at gene/protein expression level by JAM-A in breast cancer cell lines; with functional consequences in in vitro and semi-in vivo models. In bioinformatic, cellular and patient tissue models, this was associated with regulation of the HER3 transcription factor FOXA1 by JAM-A via a pathway involving β-catenin. Our data suggest a novel model whereby JAM-A expression regulates β-catenin localization, in turn regulating FOXA1 expression, which could drive HER3 gene transcription. JAM-A merits investigation as a novel target to prevent upregulation of HER3 during the development of resistance to HER2-targeted therapies, or to reduce HER3-dependent tumorigenic signaling.

Triple‑negative breast cancer therapy: Current and future perspectives (Review)

Triple-negative breast cancer (TNBC) accounts for 10-15% of all breast cancer cases. TNBCs lack estrogen and progesterone receptors and express low levels of HER2, and therefore do not respond to hormonal or anti-HER2 therapies. TNBC is a particularly aggressive form of breast cancer that generally displays poorer prognosis compared to other breast cancer subtypes. TNBC is chemotherapy sensitive, and this treatment remains the standard of care despite its limited benefit. Recent advances with novel agents have been made for specific subgroups with PD-L1⁺ tumors or germline Brca-mutated tumors. However, only a fraction of these patients responds to immune checkpoint or PARP inhibitors and even those who do respond often develop resistance and relapse. Various new agents and combination strategies have been explored to further understand molecular and immunological aspects of TNBC. In this review, we discuss clinical trials in the management of TNBC as well as perspectives for potential future treatments.

Deficit of human ornithine aminotransferase in gyrate atrophy: Molecular, cellular, and clinical aspects

Gyrate Atrophy (GA) of the choroid and retina (MIM# 258870) is an autosomal recessive disorder due to mutations of the OAT gene encoding ornithine-delta-aminotransferase (OAT), associated with progressive retinal deterioration and blindness. The disease has a theoretical global incidence of approximately 1:1,500,000. OAT is mainly involved in ornithine catabolism in adults, thus explaining the hyperornithinemia as hallmark of the disease. Patients are treated with an arginine-restricted diet, to limit ornithine load, or the administration of Vitamin B6, a precursor of the OAT coenzyme pyridoxal phosphate. Although the clinical and genetic aspects of GA are known for many years, the enzymatic phenotype of pathogenic variants and their response to Vitamin B6, as well as the molecular mechanisms explaining retinal damage, are poorly clarified. Herein, we provide an overview of the current knowledge on the biochemical properties of human OAT and on the molecular, cellular, and clinical aspects of GA.

LncRNA DSCAM-AS1 interacts with YBX1 to promote cancer progression by forming a positive feedback loop that activates FOXA1 transcription network

Rationale: The forkhead box A1 (FOXA1) is a crucial transcription factor in initiation and development of breast, lung and prostate cancer. Previous studies about the FOXA1 transcriptional network were mainly focused on protein-coding genes. Its regulatory network of long non-coding RNAs (lncRNAs) and their role in FOXA1 oncogenic activity remains unknown.

Methods: The Cancer Genome Atlas (TCGA) data, RNA-seq and ChIP-seq data were used to analyze FOXA1 regulated lncRNAs. RT-qPCR was used to detect the expression of DSCAM-AS1, RT-qPCR and Western blotting were used to determine the expression of FOXA1, estrogen receptor α (ERα) and Y box binding protein 1 (YBX1). RNA pull-down and RIP-qPCR were employed to investigate the interaction between DSCAM-AS1 and YBX1. The effect of DSCAM-AS1 on malignant phenotypes was examined through in vitro and in vivo assays.

Results: In this study, we conducted a global analysis of FOXA1 regulated lncRNAs. For detailed analysis, we chose lncRNA DSCAM-AS1, which is specifically expressed in lung adenocarcinoma, breast and prostate cancer. The expression level of DSCAM-AS1 is regulated by two super-enhancers (SEs) driven by FOXA1. High expression levels of DSCAM-AS1 was associated with poor prognosis. Knockout experiments showed DSCAM-AS1 was essential for the growth of xenograft tumors. Moreover, we demonstrated DSCAM-AS1 can regulate the expression of the master transcriptional factor FOXA1. In breast cancer, DSCAM-AS1 was also found to regulate ERα. Mechanistically, DSCAM-AS1 interacts with YBX1 and influences the recruitment of YBX1 in the promoter regions of FOXA1 and ERα.

Conclusion: Our study demonstrated that lncRNA DSCAM-AS1 was transcriptionally activated by super-enhancers driven by FOXA1 and exhibited lineage-specific expression pattern. DSCAM-AS1 can promote cancer progression by interacting with YBX1 and regulating expression of FOXA1 and ERα.

Triggering a switch from basal- to luminal-like breast cancer subtype by the small-molecule diptoindonesin G via induction of GABARAPL1

Breast cancer is a heterogeneous disease that includes different molecular subtypes. The basal-like subtype has a poor prognosis and a high recurrence rate, whereas the luminal-like subtype confers a more favorable patient prognosis partially due to anti-hormone therapy responsiveness. Here, we demonstrate that diptoindonesin G (Dip G), a natural product, exhibits robust differentiation-inducing activity in basal-like breast cancer cell lines and animal models. Specifically, Dip G treatment caused a partial transcriptome shift from basal to luminal gene expression signatures and prompted sensitization of basal-like breast tumors to tamoxifen therapy. Dip G upregulated the expression of both GABARAPL1 (GABA_A receptor-associated protein-like 1) and ERβ. We revealed a previously unappreciated role of GABARAPL1 as a regulator in the specification of breast cancer subtypes that is dependent on ERβ levels. Our findings shed light on new therapeutic opportunities for basal-like breast cancer via a phenotype switch and indicate that Dip G may serve as a leading compound for the therapy of basal-like breast cancer.

Relationships between Breast Feeding and Breast Cancer Subtypes: Lessons Learned from Studies in Humans and in Mice

There are differential risk relationships between parity and breast cancer according to estrogen receptor (ER) status, with an increased risk of ER^- disease reduced by breastfeeding. This may be particularly relevant for understanding the higher incidence of ER^- tumors in Black women, who are more likely to be parous and less likely to breastfeed than other U.S. groups. Potential mechanisms for these relationships may include effects of disordered breast involution on inflammatory milieu in the breast as well as epigenetic reprogramming in the mammary gland, which can affect cell fate decisions in progenitor cell pools. In normal breast tissue, parity has been associated with hypermethylation of FOXA1, a pioneer transcription factor that promotes the luminal phenotype in luminal progenitors, while repressing the basal phenotype. In breast tumors, relationships between FOXA1 methylation and parity were strongest among women who did not breastfeed. Here, we summarize the epidemiologic literature regarding parity, breastfeeding, and breast cancer subtypes, and review potential mechanisms whereby these factors may influence breast carcinogenesis, with a focus on effects on progenitor cell pools in the mammary gland.

Hybrid analysis of gene dynamics predicts context-specific expression and offers regulatory insights

Motivation

To understand the regulatory pathways underlying diseases, studies often investigate the differential gene expression between genetically or chemically differing cell populations. Differential expression analysis identifies global changes in transcription and enables the inference of functional roles of applied perturbations. This approach has transformed the discovery of genetic drivers of disease and possible therapies. However, differential expression analysis does not provide quantitative predictions of gene expression in untested conditions. We present a hybrid approach, termed Differential Expression in Python (DiffExPy), that uniquely combines discrete, differential expression analysis with in silico differential equation simulations to yield accurate, quantitative predictions of gene expression from time-series data.

Results

To demonstrate the distinct insight provided by DiffExpy, we applied it to published, in vitro, time-series RNA-seq data from several genetic PI3K/PTEN variants of MCF10a cells stimulated with epidermal growth factor. DiffExPy proposed ensembles of several minimal differential equation systems for each differentially expressed gene. These systems provide quantitative models of expression for several previously uncharacterized genes and uncover new regulation by the PI3K/PTEN pathways. We validated model predictions on expression data from conditions that were not used for model training. Our discrete, differential expression analysis also identified SUZ12 and FOXA1 as possible regulators of specific groups of genes that exhibit late changes in expression. Our work reveals how DiffExPy generates quantitatively predictive models with testable, biological hypotheses from time-series expression data.

Availability and implementation

DiffExPy is available on GitHub (https://github.com/bagherilab/diffexpy).

Supplementary information

Supplementary data are available at Bioinformatics online.

Gene expression profiles for low-dose exposure to diethyl phthalate in rodents and humans: a translational study with implications for breast carcinogenesis

Phthalates are commonly included as ingredients in personal care products such as cosmetics, shampoos and perfumes. Diethyl phthalate (DEP) has been found to be anti-androgenic and linked with adverse reproductive effects on males, but effects on females are poorly understood. We designed an integrative and translational study to experimentally examine the effects of DEP exposure at a human-equivalent dose on the mammary transcriptome in rats and to subsequently examine the DEP gene signature in breast tissues (both pre-malignant and tumor) from a population study. In Sprague-Dawley rats treated orally with DEP from birth to adulthood, we identified a signature panel of 107 genes predominantly down-regulated by DEP exposure. Univariate analysis of this 107 DEP gene signature in pre-malignant breast tissues revealed that six genes (P4HA1, MPZL3, TMC4, PLEKHA6, CA8, AREG) were inversely associated with monoethyl phthalate (MEP; the urinary metabolite of DEP) concentration (p < 0.05) among postmenopausal women; all six genes loaded on to one of seven factors identified by factor analysis. Transcription factor enrichment analysis revealed that genes in this factor were enriched for androgen receptor binding sites. These six genes were also significantly down-regulated in pre-malignant adjacent tissues compared to the corresponding tumor tissues in pair-wise analyses (p < 0.05). Results from our translational study indicate that low level exposure to diethyl phthalate results in measurable genomic changes in breast tissue with implications in breast carcinogenesis.

Nuclear accumulation of MKL1 in luminal breast cancer cells impairs genomic activity of ERα and is associated with endocrine resistance

Estrogen receptor (ERα) is central in driving the development of hormone-dependent breast cancers. A major challenge in treating these cancers is to understand and overcome endocrine resistance. The Megakaryoblastic Leukemia 1 (MKL1, MRTFA) protein is a master regulator of actin dynamic and cellular motile functions, whose nuclear translocation favors epithelial-mesenchymal transition. We previously demonstrated that nuclear accumulation of MKL1 in estrogen-responsive breast cancer cell lines promotes hormonal escape. In the present study, we confirm through tissue microarray analysis that nuclear immunostaining of MKL1 is associated with endocrine resistance in a cohort of breast cancers and we decipher the underlining mechanisms using cell line models. We show through gene expression microarray analysis that the nuclear accumulation of MKL1 induces dedifferentiation leading to a mixed luminal/basal phenotype and suppresses estrogen-mediated control of gene expression. Chromatin immunoprecipitation of DNA coupled to high-throughput sequencing (ChIP-Seq) shows a profound reprogramming in ERα cistrome associated with a massive loss of ERα binding sites (ERBSs) generally associated with lower ERα-binding levels. Novel ERBSs appear to be associated with EGF and RAS signaling pathways. Collectively, these results highlight a major role of MKL1 in the loss of ERα transcriptional activity observed in certain cases of endocrine resistances, thereby contributing to breast tumor cells malignancy.

Cooperative driver pathway discovery via fusion of multi-relational data of genes, miRNAs and pathways

Discovering driver pathways is an essential step to uncover the molecular mechanism underlying cancer and to explore precise treatments for cancer patients. However, due to the difficulties of mapping genes to pathways and the limited knowledge about pathway interactions, most previous work focus on identifying individual pathways. In practice, two (or even more) pathways interplay and often cooperatively trigger cancer. In this study, we proposed a new approach called CDPathway to discover cooperative driver pathways. First, CDPathway introduces a driver impact quantification function to quantify the driver weight of each gene. CDPathway assumes that genes with larger weights contribute more to the occurrence of the target disease and identifies them as candidate driver genes. Next, it constructs a heterogeneous network composed of genes, miRNAs and pathways nodes based on the known intra(inter)-relations between them and assigns the quantified driver weights to gene-pathway and gene-miRNA relational edges. To transfer driver impacts of genes to pathway interaction pairs, CDPathway collaboratively factorizes the weighted adjacency matrices of the heterogeneous network to explore the latent relations between genes, miRNAs and pathways. After this, it reconstructs the pathway interaction network and identifies the pathway pairs with maximal interactive and driver weights as cooperative driver pathways. Experimental results on the breast, uterine corpus endometrial carcinoma and ovarian cancer data from The Cancer Genome Atlas show that CDPathway can effectively identify candidate driver genes [area under the receiver operating characteristic curve (AUROC) of $\geq $0.9] and reconstruct the pathway interaction network (AUROC of>0.9), and it uncovers much more known (potential) driver genes than other competitive methods. In addition, CDPathway identifies 150% more driver pathways and 60% more potential cooperative driver pathways than the competing methods. The code of CDPathway is available at http://mlda.swu.edu.cn/codes.php?name=CDPathway.

Racial Differences in the Association Between Luminal Master Regulator Gene Expression Levels and Breast Cancer Survival

PURPOSE

Compared with their European American (EA) counterparts, African American (AA) women are more likely to die from breast cancer in the United States. This disparity is greatest in hormone receptor-positive subtypes. Here we uncover biological factors underlying this disparity by comparing functional expression and prognostic significance of master transcriptional regulators of luminal differentiation.

EXPERIMENTAL DESIGN

Data and biospecimens from 262 AA and 293 EA patients diagnosed with breast cancer from 2001 to 2010 at a major medical center were analyzed by IHC for functional biomarkers of luminal differentiation, including estrogen receptor (ESR1) and its pioneer factors, FOXA1 and GATA3. Integrated comparison of protein levels with network-level gene expression analysis uncovered predictive correlations with race and survival.

RESULTS

Univariate or multivariate HRs for overall survival, estimated from digital IHC scoring of nuclear antigen, show distinct differences in the magnitude and significance of these biomarkers to predict survival based on race: ESR1 [EA HR = 0.47; 95% confidence interval (CI), 0.31-0.72 and AA HR = 0.77; 95% CI, 0.48-1.18]; FOXA1 (EA HR = 0.38; 95% CI, 0.23-0.63 and AA HR = 0.53; 95% CI, 0.31-0.88), and GATA3 (EA HR = 0.36; 95% CI, 0.23-0.56; AA HR = 0.57; CI, 0.56-1.4). In addition, we identify genes in the downstream regulons of these biomarkers highly correlated with race and survival.

CONCLUSIONS

Even within clinically homogeneous tumor groups, regulatory networks that drive mammary luminal differentiation reveal race-specific differences in their association with clinical outcome. Understanding these biomarkers and their downstream regulons will elucidate the intrinsic mechanisms that drive racial disparities in breast cancer survival.

FOXA1 Protein Expression in ER+ and ER- Breast Cancer in Relation to Parity and Breastfeeding in Black and White Women

BACKGROUND

Forkhead box protein A1 (FOXA1) promotes luminal differentiation, and hypermethylation of the gene can be a mechanism of developing estrogen receptor-negative (ER^-) breast cancer. We examined FOXA1 in breast tumor and adjacent normal tissue in relation to reproductive factors, particularly higher parity and no breastfeeding, that are associated with ER^- tumors.

METHODS

We performed IHC for FOXA1 in breast tumors (n = 1,329) and adjacent normal tissues (n = 298) in the Women's Circle of Health Study (949 Blacks and 380 Whites). Protein expression levels were summarized by histology (H) scores. Generalized linear models were used to assess FOXA1 protein expression in relation to reproductive factors by ER status.

RESULTS

ER-positive (ER⁺) versus ER^- tumors had higher FOXA1 protein expression (P < 0.001). FOXA1 expression was higher in tumor versus paired adjacent normal tissue in women with ER⁺ or non-triple-negative cancer (both P < 0.001), but not in those with ER^- or triple-negative cancer. Higher number of births (1, 2, and 3+) was associated with lower FOXA1 protein expression in ER⁺ tumors [differences in H score, or β = -8.5; 95% confidence interval (CI), -15.1 to -2.0], particularly among parous women who never breastfed (β = -10.4; 95% CI, -19.7 to -1.0), but not among those who breastfed (β = -7.5; 95% CI, -16.9 to 1.8). The associations for ER^- tumors were similar, although they were not statistically significant.

CONCLUSIONS

In this tumor-based study, higher parity was associated with lower FOXA1 expression in ER⁺ tumors, and breastfeeding may ameliorate the influence.

IMPACT

These findings contribute to our understanding of FOXA1 methylation and breast cancer etiology.

Novel Breast-Specific Long Non-coding RNA LINC00993 Acts as a Tumor Suppressor in Triple-Negative Breast Cancer

Background: Triple-negative breast cancer (TNBC) was characterized by breast cancers that do not express estrogen receptor (ER), progesterone receptor (PR), or human epidermal growth factor receptor (HER)-2 genes. TNBC patients are associated with a shorter median time to relapse and death for the lack of available treatment targets. Long non-coding RNAs (LncRNAs) have been reported to play an important role in the development of TNBC. We identified a novel breast-specific long non-coding RNA LINC00993, but less was known about its expression pattern and functional role in TNBC.

Methods: LINC00993 RNA expression was detected across different types of clinical breast cancer samples by using qRT-PCR. Bioinformatic methods “guilt by association” and gene set enrichment analysis (GSEA) were used to predict LINC00993 functions. Subcellular localization of LINC00993 in cells was detected by RNA fluorescence in situ hybridization (FISH). Effect of LINC00993 on cell growth was measured by plate colony formation assays, typical growth curve, and an in vivo tumor model. Cell cycle analysis was done by flow cytometry analysis. Key cell cycle regulators were detected by Western blot.

Results: LINC00993 was largely downregulated in TNBC, and higher expression indicated better outcome. LINC00993 located mainly in the nucleus. LINC00993 suppressed TNBC growth both in vitro and in vivo. LINC00993 was predicted to be involved in cell cycle pathways by using “guilt by association” and GSEA methods. Key cell cycle regulators like p16^INK4A, p14^ARF, p53, and p21 were affected by LINC00993 overexpression.

Conclusions: A new breast-specific lincRNA LINC00993 was identified with a tumor-suppressive feature and with prognostic value. This is the first research on LINC00993 function. Our results suggest that controlling LINC00993 level may be beneficial for breast cancer treatment.