Expression of FOXA1 and GATA-3 in breast cancer: the prognostic significance in hormone receptor-negative tumours

FOXA1 forms biomolecular condensates that unpack condensed chromatin to function as a pioneer factor

Eukaryotic DNA is packaged into chromatin in the nucleus, restricting the binding of transcription factors (TFs) to their target DNA sites. FOXA1 functions as a pioneer TF to bind condensed chromatin and initiate the opening of local chromatin for gene expression. However, the principles of FOXA1 recruitment and how it subsequently unpacks the condensed chromatin remain elusive. Here, we revealed that FOXA1 intrinsically forms submicron-sized condensates through its N- and C-terminal intrinsically disordered regions (IDRs). Notably, both IDRs enable FOXA1 to dissolve the condensed chromatin. In addition, the DNA-binding capacity of FOXA1 contributes to its ability to both form condensates and dissolve condensed chromatin. Further genome-wide investigation showed that IDRs enable FOXA1 to bind and unpack the condensed chromatin to regulate the proliferation and migration of breast cancer cells. This work provides a principle of how pioneer TFs function to initiate competent chromatin states using their IDRs.

Twist alters the breast tumor microenvironment via choline kinase to facilitate an aggressive phenotype

Twist (TWIST1) is a gene required for cell fate specification in embryos and its expression in mammary epithelium can initiate tumorigenesis through the epithelial-mesenchymal transition. To identify downstream target genes of Twist in breast cancer, we performed microarray analysis on the transgenic breast cancer cell line, MCF-7/Twist. One of the targets identified was choline kinase whose upregulation resulted in increased cellular phosphocholine and total choline containing compounds-a characteristic observed in highly aggressive metastatic cancers. To study the interactions between Twist, choline kinase, and their effect on the microenvironment, we used 1H magnetic resonance spectroscopy and found significantly higher phosphocholine and total choline, as well as increased phosphocholine/glycerophosphocholine ratio in MCF-7/Twist cells. We also observed significant increases in extracellular glucose, lactate, and [H +] ion concentrations in the MCF-7/Twist cells. Magnetic resonance imaging of MCF-7/Twist orthotopic breast tumors showed a significant increase in vascular volume and permeability surface area product compared to control tumors. In addition, by reverse transcription-quantitative polymerase chain reaction, we discovered that Twist upregulated choline kinase expression in estrogen receptor negative breast cancer cell lines through FOXA1 downregulation. Moreover, using The Cancer Genome Atlas database, we observed a significant inverse relationship between FOXA1 and choline kinase expression and propose that it could act as a modulator of the Twist/choline kinase axis. The data presented indicate that Twist is a driver of choline kinase expression in breast cancer cells via FOXA1 resulting in the generation of an aggressive breast cancer phenotype.

Recovering false negatives in CRISPR fitness screens with JLOE

It is widely accepted that pooled library CRISPR knockout screens offer greater sensitivity and specificity than prior technologies in detecting genes whose disruption leads to fitness defects, a critical step in identifying candidate cancer targets. However, the assumption that CRISPR screens are saturating has been largely untested. Through integrated analysis of screen data in cancer cell lines generated by the Cancer Dependency Map, we show that a typical CRISPR screen has a ∼20% false negative rate, in addition to library-specific false negatives. Replicability falls sharply as gene expression decreases, while cancer subtype-specific genes within a tissue show distinct profiles compared to false negatives. Cumulative analyses across tissues improves our understanding of core essential genes and suggest only a small number of lineage-specific essential genes, enriched for transcription factors that define pathways of tissue differentiation. To recover false negatives, we introduce a method, Joint Log Odds of Essentiality (JLOE), which builds on our prior work with BAGEL to selectively rescue the false negatives without an increased false discovery rate.

The role of the NDRG1 in the pathogenesis and treatment of breast cancer

Breast cancer (BC) is the leading cause of cancer death in women. This disease is heterogeneous, with clinical subtypes being estrogen receptor-α (ER-α) positive, having human epidermal growth factor receptor 2 (HER2) overexpression, or being triple-negative for ER-α, progesterone receptor, and HER2 (TNBC). The ER-α positive and HER2 overexpressing tumors can be treated with agents targeting these proteins, including tamoxifen and pertuzumab, respectively. Despite these treatments, resistance and metastasis are problematic, while TNBC is challenging to treat due to the lack of suitable targets. Many studies examining BC and other tumors indicate a role for N-myc downstream-regulated gene-1 (NDRG1) as a metastasis suppressor. The ability of NDRG1 to inhibit metastasis is due, in part, to the inhibition of the initial step in metastasis, namely the epithelial-to-mesenchymal transition. Paradoxically, there are also reports of NDRG1 playing a pro-oncogenic role in BC pathogenesis. The oncogenic effects of NDRG1 in BC have been reported to relate to lipid metabolism or the mTOR signaling pathway. The molecular mechanism(s) of how NDRG1 regulates the activity of multiple signaling pathways remains unclear. Therapeutic strategies that up-regulate NDRG1 have been developed and include agents of the di-2-pyridylketone thiosemicarbazone class. These compounds target oncogenic drivers in BC cells, suppressing the expression of multiple key hormone receptors including ER-α, progesterone receptor, androgen receptor, and prolactin receptor, and can also overcome tamoxifen resistance. Considering the varying role of NDRG1 in BC pathogenesis, further studies are required to examine what subset of BC patients would benefit from pharmacopeia that up-regulate NDRG1.

GATA3 Expression in Human Tumors: A Tissue Microarray Study on 16,557 Tumors

INTRODUCTION

GATA3 is a transcription factor involved in epithelial cell differentiation. GATA3 immunostaining is used as a diagnostic marker for breast and urothelial cancer but can also occur in other neoplasms.

METHODS

To evaluate GATA3 in normal and tumor tissues, a tissue microarray containing 16,557 samples from 131 different tumor types and subtypes and 608 samples of 76 different normal tissue types was analyzed by immunohistochemistry.

RESULTS

GATA3 positivity was found in 69 different tumor types including 23 types (18%) with at least one strongly positive tumor. Highest positivity rates occurred in noninvasive papillary urothelial carcinoma (92-99%), lobular carcinoma (98%), carcinoma of no special type of the breast (92%), basal cell carcinoma of the skin (97%), invasive urothelial carcinoma (73%), T-cell lymphoma (23%), adenocarcinoma of the salivary gland (16%), squamous cell carcinoma of the skin (16%), and colorectal neuroendocrine carcinoma (12%). In breast cancer, low GATA3 staining was linked to high pT stage (p = 0.03), high BRE grade (p < 0.0001), HER2 overexpression (p = 0.0085), estrogen and progesterone receptor negativity (p < 0.0001 each), and reduced survival (p = 0.03).

CONCLUSION

Our data demonstrate that GATA3 positivity can occur in various tumor entities. Low levels of GATA3 reflect cancer progression and poor patient prognosis in breast cancer.

A framework for clinical cancer subtyping from nucleosome profiling of cell-free DNA

Cell-free DNA (cfDNA) has the potential to inform tumor subtype classification and help guide clinical precision oncology. Here we develop Griffin, a framework for profiling nucleosome protection and accessibility from cfDNA to study the phenotype of tumors using as low as 0.1x coverage whole genome sequencing data. Griffin employs a GC correction procedure tailored to variable cfDNA fragment sizes, which generates a better representation of chromatin accessibility and improves the accuracy of cancer detection and tumor subtype classification. We demonstrate estrogen receptor subtyping from cfDNA in metastatic breast cancer. We predict estrogen receptor subtype in 139 patients with at least 5% detectable circulating tumor DNA with an area under the receive operator characteristic curve (AUC) of 0.89 and validate performance in independent cohorts (AUC = 0.96). In summary, Griffin is a framework for accurate tumor subtyping and can be generalizable to other cancer types for precision oncology applications.

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.

Value of Proline, Glutamic Acid, and Leucine-Rich Protein 1 and GATA Binding Protein 3 Expression in Breast Cancer: An Immunohistochemical study

GATA binding protein 3 was more sensitive than traditional markers such as gross cystic disease fluid protein 15 and mammaglobin for identifying primary and metastatic breast carcinomas, but its significance decreased in triple-negative breast cancer. Recent studies showed a high expression rate of proline glutamic acid and leucine-rich protein in breast cancer and their superiority over GATA3 in triple-negative breast cancer. Our study provided new insights into the diagnostic and prognostic roles of PELP1 and GATA3 in primary and metastatic breast cancer. An immunohistochemical assay was carried out using PELP1 and GATA3 in 60 cases of primary breast cancer and 15 metastatic. Invasive carcinoma of no special type was the predominant type (80%). The majority of cases were grade 3 (68.3%). GATA3 expression was 83.3% positive in primary breast carcinomas and 73.5% positive in metastatic breast carcinomas. In comparison, PELP1 had a 96.7% positive expression rate in primary breast carcinomas and an 86.7% positive expression rate in metastasis. There was a statistically significant agreement between GATA3 and PELP1 in the diagnosis of the cases. PELP1 is a significantly higher proportion of both primary and metastatic breast carcinomas than GATA3. In breast cancer, there was a strong association between favorable prognostic factors and GATA3 expression, with evidence of an inverse association with Ki-67 overexpression.

Decidua Parietalis Mesenchymal Stem/Stromal Cells and Their Secretome Diminish the Oncogenic Properties of MDA231 Cells In Vitro

Mesenchymal stem cells (MSCs) have been shown to suppress tumor growth, inhibit angiogenesis, regulate cellular signaling, and induce apoptosis in cancer cells. We have earlier reported that placenta-derived decidua parietalis mesenchymal stem/stromal cells (DPMSCs) not only retained their functional characteristics in the cancer microenvironment but also exhibited increased expression of anti-apoptotic genes, demonstrating their anti-tumor properties in the tumor setting. In this study, we have further evaluated the effects of DPMSCs on the functional outcome of human breast cancer cell line MDA231. MDA231 cells were exposed to DPMSCs, and their biological functions, including adhesion, proliferation, migration, and invasion, were evaluated. In addition, genomic and proteomic modifications of the MDA231 cell line, in response to the DPMSCs, were also evaluated. MDA231 cells exhibited a significant reduction in proliferation, migration, and invasion potential after their treatment with DPMSCs. Furthermore, DPMSC treatment diminished the angiogenic potential of MDA231 cells. DPMSC treatment modulated the expression of various pro-apoptotic as well as oncogenes in MDA231 cells. The properties of DPMSCs to inhibit the invasive characteristics of MDA231 cells demonstrate that they may be a useful candidate in a stem-cell-based therapy against cancer.

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.

ZEB2 regulates endocrine therapy sensitivity and metastasis in luminal a breast cancer cells through a non-canonical mechanism

PURPOSE

The transcription factors ZEB1 and ZEB2 mediate epithelial-to-mesenchymal transition (EMT) and metastatic progression in numerous malignancies including breast cancer. ZEB1 and ZEB2 drive EMT through transcriptional repression of cell-cell junction proteins and members of the tumor suppressive miR200 family. However, in estrogen receptor positive (ER +) breast cancer, the role of ZEB2 as an independent driver of metastasis has not been fully investigated.

METHODS

In the current study, we induced exogenous expression of ZEB2 in ER + MCF-7 and ZR-75-1 breast cancer cell lines and examined EMT gene expression and metastasis using dose-response qRT-PCR, transwell migration assays, proliferation assays with immunofluorescence of Ki-67 staining. We used RNA sequencing to identify pathways and genes affected by ZEB2 overexpression. Finally, we treated ZEB2-overexpressing cells with 17β-estradiol (E2) or ICI 182,780 to evaluate how ZEB2 affects estrogen response.

RESULTS

Contrary to expectation, we found that ZEB2 did not increase canonical epithelial nor decrease mesenchymal gene expressions. Furthermore, ZEB2 overexpression did not promote a mesenchymal cell morphology. However, ZEB1 and ZEB2 protein expression induced significant migration of MCF-7 and ZR-75-1 breast cancer cells in vitro and MCF-7 xenograft metastasis in vivo. Transcriptomic (RNA sequencing) pathway analysis revealed alterations in estrogen signaling regulators and pathways, suggesting a role for ZEB2 in endocrine sensitivity in luminal A breast cancer. Expression of ZEB2 was negatively correlated with estrogen receptor complex genes in luminal A patient tumors. Furthermore, treatment with 17β-estradiol (E2) or the estrogen receptor antagonist ICI 182,780 had no effect on growth of ZEB2-overexpressing cells.

CONCLUSION

ZEB2 is a multi-functional regulator of drug sensitivity, cell migration, and metastasis in ER + breast cancer and functions through non-canonical mechanisms.

Analysis of differential expression of hypoxia-inducible microRNA-210 gene targets in mild and severe preeclamptic patients

Preeclampsia (PE) is a multi-system disorder that is specific to human pregnancy. Inadequate oxygenation of uterus and placenta is considered as one of the leading causes for the disease. MicroRNA-210(miR-210) is one of the prime molecules that has emerged in response to hypoxia. The objective of this study was to determine miR-210 expression patterns in plasma from severe PE and mild PE patients, and how that affects the expression of miR-210 target genes. The expression levels of miR-210 were validated using reverse transcription-quantitative PCR in plasma of severe PE (15) and mild PE (15) patients in comparison to controls subjects (15) with normal pregnancy. Then, the association between miR-210 and its downstream genes was validated by using human miR-210 targets RT2 profiler PCR Array. Both the categories (mild and severe) showed significantly high miR-210 expression levels. Also out of the 84 hypoxia miR-210 associated genes screened using mRNA, 18 genes were found to be differentially expressed in severe PE whereas 16 genes in mild PE cases with varying magnitude. All the genes in both the PE groups were found downregulated in comparison to controls. These downregulated genes expressed in both the cases were shown to be participating in immunosuppression, apoptosis, cell growth, signaling, angiogenesis, DNA repair. This study provides novel data on the genes that work downstream of miR-210 and how dysregulated expression of miR-210 can affect their expression and in turn functioning which can be associated with PE risk and severity. This study is the very first to determine the effect of miR-210 expression levels on associated genes in plasma samples.

GATA3 as an Adjunct Prognostic Factor in Breast Cancer Patients with Less Aggressive Disease: A Study with a Review of the Literature

BACKGROUND

GATA binding protein 3 (GATA3) expression is positively correlated with estrogen receptor (ER) expression, but its prognostic value as an independent factor remains unclear. Thus, we undertook the current study to evaluate the expression of GATA3 and its prognostic value in a large series of breast carcinomas (BCs) with long-term follow-up.

METHODS

A total of 702 consecutive primary invasive BCs resected between 1989 and 1993 in our institution were arranged in tissue microarrays, immunostained for ER, progesterone receptor (PR), ki-67, HER2, p53, and GATA3, and scored. Clinico-pathological data were retrospectively collected.

RESULTS

GATA3 was evaluable in 608 (87%) of the 702 cases; it was positive in 413 (68%) cases and negative in 195 (32%) cases. GATA3 positivity was significantly associated with lower grade (p < 0.0001), size (p = 0.0463), stage (p = 0.0049), ER+ (p < 0.0001), PR+ (p < 0.0001), HER2- (p = 0.0175), and p53 wild-type pattern (p < 0.0001). The median follow-up was 183 months, GATA3 positivity was associated with better overall survival (HR 0.70, p = 0.001), and its prognostic value was retained in a multivariate analysis. The association with better overall survival was stronger in patients with grade 1-2, pT1-2, pN0, stage I-II, ER+, PR+, ki-67 < 20%, HER2-, a wild-type p53 immunohistochemical pattern, and in luminal B BC.

CONCLUSIONS

Our findings indicate that GATA3 is a positive prognostic marker in BC patients, especially in patients with biologically less aggressive BC. Incorporating GATA3 immunohistochemistry into routine practice could help further stratify BC patients for their risk.

GATA binding protein 3 (GATA3) as a marker for metaplastic spindle cell carcinoma of the breast

Spindle cell lesions of the breast comprise a diverse set of tumors; harboring significant histological and immunohistochemical (IHC) overlap. Accurate diagnosis and classification of spindle cell lesions in the breast remains challenging, especially in core biopsies. In the current study, we evaluated a spectrum of spindle cell lesion of the breast with a panel of IHC antibodies in an effort to differentiate metaplastic spindle cell carcinoma from its benign and malignant mimickers. Our study included 92 patients who underwent breast core biopsies or breast resections at Northwell Health who were diagnosed with benign and malignant tumor/tumor-like spindle cell lesions. Tumors subtypes in this the study included: angiosarcoma, nodular fasciitis, fibromatosis, myofibroblastoma, phyllodes tumors (benign, borderline and malignant), primary sarcomas and metaplastic spindle cell carcinoma. Our biomarker panel included high molecular weight keratin (HMWK), CAM5.2, AE1/AE3, p63, CD34 and GATA3. GATA3 expression was significantly higher in metaplastic carcinomas (88.9 % vs 4.1 %, p < 0.001), when compared to other spindle cell lesions. The sensitivity and specificity for detecting metaplastic carcinomas reached 84.2 % and 97.3 %, respectively. Regarding cytokeratin panels, none of the three individual markers were as sensitive or specific for metaplastic breast carcinoma. GATA3 is the most specific and sensitive marker forfor the identification of metaplastic spindle cell carcinoma of the breast.

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.

Time-varying effects of FOXA1 on breast cancer prognosis

PURPOSE

Results of previous studies on the associations between Forkhead box A1 (FOXA1) expression in breast cancer tissues and the prognosis varied depending on the follow-up durations. The present study would investigate whether there is a time-varying effect of FOXA1 in breast cancer tissues on the prognosis.

METHODS

FOXA1 expressions were evaluated in 1041 primary invasive breast tumors with tissue microarrays by immunohistochemistry. Cox models with restricted cubic splines and Kaplan-Meier survival analysis were used to examine the associations between FOXA1 and the prognosis. Flexible parametric models were applied to explore the time-varying effect of FOXA1.

RESULTS

Overall, the association between FOXA1 expression and the prognosis was not significant but varied on the time of follow-up. Compared to FOXA1 ≤ 270 of H-score, the hazard ratios (HRs) of death for those with 271-285 of FOXA1 expression increased from 0.35 (95% CI 0.14-0.86) at 6 months after diagnosis to 2.88 (95% CI 1.35-6.15) at 120 months with a crossover at around 36 months. Similar patterns were also observed for FOXA1 > 285 of H-score and for progression free survival (PFS). Moreover, when allowed both FOXA1 and estrogen receptor (ER) to change over time in the model (considering that ER had a similar time-varying effect), these time-varying effects remained for FOXA1 on both overall survival (OS) (P < 0.01) and PFS (P = 0.01) but were attenuated for ER (P = 0.13 for OS).

CONCLUSIONS

This study revealed an independent time-varying effect of FOXA1 on breast cancer prognosis, which would provide an insight into the roles of FOXA1 as a marker of breast cancer prognosis and may help optimize the medication strategies.

GATA3 somatic mutations are associated with clinicopathological features and expression profile in TCGA breast cancer patients

The effect of somatic mutations and the gene expression profiles on the prognosis is well documented in cancer research. This study was conducted to evaluate the association of GATA3 somatic mutations with tumor features, survival, and expression profiles in breast cancer. Clinicopathological information was compared between TCGA-BRCA patients with GATA3-mutant and non-mutant tumors in all patients as well as in ER-positive subgroup. Cox-regression method was used to evaluate the association of the GATA3 mutation status with overall survival time. Differential gene expression, functional annotation, and protein-protein interaction analyses were performed using edgeR, Metascape, DAVID, STRING and CytoNCA. GATA3-mutant and non-mutant samples had significantly different clinicopathological features (p < 0.05). While GATA3 mutation status was not associated with the overall survival in the entire cohort (p_adj = 0.52), the GATA3-wild type ER-positive cases had a better prognosis than mutant ones (p_adj = 0.04). GATA3 expression was higher in tumors than normal tissues. Several pathways were different between mutant and non-mutant groups (p < 0.05). Interleukin-6 was found as the highest scored gene in both comparisons (normal vs. mutant and normal vs. non-mutant groups) in the entire patient and in the ER-positive subgroup, suggesting the association of IL6 with breast tumorigenesis. These findings suggest that GATA3 mutations can be associated with several tumor characteristics and influence the pattern of gene expression. However, GATA3 mutation status seems to be a prognostic factor for the disease only in ER-positive patients.

Functionally analyzing the important roles of hepatocyte nuclear factor 3 (FoxA) in tumorigenesis

Transcriptional factors (TFs) play a central role in governing gene expression under physiological conditions including the processes of embryonic development, metabolic homeostasis and response to extracellular stimuli. Conceivably, the aberrant dysregulations of TFs would dominantly result in various human disorders including tumorigenesis, diabetes and neurodegenerative diseases. Serving as the most evolutionarily reserved TFs, Fox family TFs have been explored to exert distinct biological functions in neoplastic development, by manipulating diverse gene expression. Recently, among the Fox family members, the pilot roles of FoxAs attract more attention due to their functions as both pioneer factor and transcriptional factor in human tumorigenesis, particularly in the sex-dimorphism tumors. Therefore, the pathological roles of FoxAs in tumorigenesis have been well-explored in modulating inflammation, immune response and metabolic homeostasis. In this review, we comprehensively summarize the impressive progression of FoxA functional annotation, clinical relevance, upstream regulators and downstream effectors, as well as valuable animal models, and highlight the potential strategies to target FoxAs for cancer therapies.

Bioinformatics analyses and biological function of lncRNA ZFPM2-AS1 and ZFPM2 gene in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) remains one of the most lethal malignant tumors worldwide; however, the etiology of HCC still remains poorly understood. In the present study, cancer-omics databases, including The Cancer Genome Atlas, GTEx and Gene Expression Omnibus, were systematically analyzed in order to investigate the role of the long non-coding RNA (lncRNA) zinc finger protein, FOG family member 2-antisense 1 (ZFPM2-AS1) and the zinc finger protein, FOG family member 2 (ZFPM2) gene in the occurrence and progression of HCC. It was identified that the expression levels of lncRNA ZFPM2-AS1 were significantly increased in HCC tissues, whereas expression levels of the ZFPM2 gene were significantly decreased in HCC tissues compared with normal liver tissues. Higher expression levels of ZFPM2-AS1 were significantly associated with a less favorable prognosis of HCC, whereas higher expression levels of the ZFPM2 gene were associated with a more favorable prognosis of HCC. Genetic alterations in the ZFPM2 gene may contribute to a worse prognosis of HCC. Validation of the GSE14520 dataset also demon stared that ZFPM2 gene expression levels were significantly decreased in HCC tissues (P<0.001). The receiver operating characteristic (ROC) analysis of the ZFPM2 gene indicated high accuracy of this gene in distinguishing between HCC tissues and non-tumor tissues. The areas under the ROC curves were >0.8. Using integrated strategies, the present study demonstrated that lncRNA ZFPM2-AS1 and the ZFPM2 gene may contribute to the occurrence and prognosis of HCC. These findings may provide a novel understanding of the molecular mechanisms underlying the occurrence and prognosis of HCC.

Learning common and specific patterns from data of multiple interrelated biological scenarios with matrix factorization

High-throughput biological technologies (e.g. ChIP-seq, RNA-seq and single-cell RNA-seq) rapidly accelerate the accumulation of genome-wide omics data in diverse interrelated biological scenarios (e.g. cells, tissues and conditions). Integration and differential analysis are two common paradigms for exploring and analyzing such data. However, current integrative methods usually ignore the differential part, and typical differential analysis methods either fail to identify combinatorial patterns of difference or require matched dimensions of the data. Here, we propose a flexible framework CSMF to combine them into one paradigm to simultaneously reveal Common and Specific patterns via Matrix Factorization from data generated under interrelated biological scenarios. We demonstrate the effectiveness of CSMF with four representative applications including pairwise ChIP-seq data describing the chromatin modification map between K562 and Huvec cell lines; pairwise RNA-seq data representing the expression profiles of two different cancers; RNA-seq data of three breast cancer subtypes; and single-cell RNA-seq data of human embryonic stem cell differentiation at six time points. Extensive analysis yields novel insights into hidden combinatorial patterns in these multi-modal data. Results demonstrate that CSMF is a powerful tool to uncover common and specific patterns with significant biological implications from data of interrelated biological scenarios.

Avoidance and Period-Shortening of Neoadjuvant Chemotherapy Against Triple-Negative Breast Cancer in Stages I and II: Importance of Ki-67 Labeling Index and the Recognition of Apocrine-Type Lesions

BACKGROUND

Triple-negative breast cancer encompasses heterogeneous subtypes. Neoadjuvant chemotherapy is ineffective against some triple-negative breast cancers, while others show a favorable prognosis despite chemoresistance.

METHODS

A total of 51 cases with stages I and II triple-negative breast cancer were analyzed; 34 triple-negative breast cancers treated with neoadjuvant chemotherapy were divided into "good responders" (n = 22), showing therapeutic effect G2b or G3 in surgical specimens, and "poor responders" with therapeutic effect G0, G1a, G1b, and G2a (n = 12). Neoadjuvant chemotherapy was spared in 17 cases (non-neoadjuvant chemotherapy group). Apocrine-type triple-negative breast cancer was defined as triple-negative breast cancer immunoreactive for both androgen receptor and forkhead-box protein A1. Triple-negative breast cancer other than apocrine-type (n = 16) and special types (myoepithelial, medullary, adenoid cystic, and spindle cell carcinomas, n = 6) was categorized as basal-like subtype (n = 29). Prognosis was evaluated in each category.

RESULTS

Neoadjuvant chemotherapy provoked significant effects against basal-like triple-negative breast cancer with high Ki-67 labeling (≧50%), and tumor-infiltrating lymphocytes predicted high chemosensitivity. Neoadjuvant chemotherapy was avoidable in triple-negative breast cancer of apocrine- and special types showing low (<50%) Ki-67 labeling. Ten (59%) lesions in the non-neoadjuvant chemotherapy group belonged to the apocrine-type. When clinical complete remission shown by contrast-enhanced magnetic resonance imaging was reached in the course of neoadjuvant chemotherapy against basal-like triple-negative breast cancer, the neoadjuvant chemotherapy period was shortened in 14 (64%) of 22 good responders. Disease-free and overall survival rates were excellent in all groups.

CONCLUSIONS

The following 2 hypothetical proposals should be proven by large-scale clinical trials. Immunohistochemical recognition of apocrine-type triple-negative breast cancer with low Ki-67 labeling is important for avoiding ineffective/unnecessary neoadjuvant chemotherapy. By employing appropriate clinical imaging, period-shortening is achievable in basal-like triple-negative breast cancer with high Ki-67 labeling.

Androgen receptor and FOXA1 coexpression define a "luminal-AR" subtype of feline mammary carcinomas, spontaneous models of breast cancer

Background

Invasive mammary carcinomas that spontaneously develop in female cats are associated with high mortality, and resemble the most aggressive human breast cancers, especially triple-negative breast cancer (TNBC). Transcriptome studies showed that TNBCs are a heterogeneous group that includes a potentially hormone-dependent subtype named luminal-AR. Some authors proposed an immunohistochemical definition of the luminal-AR subtype, which is not only positive for Androgen Receptor (AR), but also either positive for the transcription factor Forkhead box A1 (FOXA1), or negative for basal markers. The objectives of this study were to describe AR and FOXA1 expressions in feline mammary carcinomas (FMCs), their prognostic value, and if their coexpression could define a “luminal-AR” subtype of triple-negative mammary carcinomas in cats.

Methods

In a previously described retrospective cohort of 180 female cats with FMCs, with a 2-year follow-up post-mastectomy, we assessed AR, FOXA1, ER, PR, Ki-67, HER2, and CK14 expressions by automated immunohistochemistry.

Results

Of the 180 FMCs, 57 (32%) were luminal; i.e., ER and/or PR positive, and 123 (68%) were triple-negative (ER–, PR– and HER2–) FMCs. AR overexpression (found in 33 cases/180, 18%) and FOXA1 index ≥1% (64/180, 36%) were associated with a longer disease-free interval, overall survival, and cancer-specific survival in cats with FMC. Analysis of AR, FOXA1 and CK14 coexpression in triple-negative FMCs showed that AR+ triple-negative FMCs were heterogeneous: there existed an AR+ FOXA1+ CK14– subgroup (n = 7) associated with a better cancer-specific survival by multivariate survival analysis (HR = 0.26, 95% CI: 0.07–0.89, p = 0.03) compared to AR+ FOXA1–CK14+ triple-negative FMCs (n = 46) (HR = 1.00), independently of the pathologic tumor size and pathologic nodal stage. The non-basal-like subtype of triple-negative FMCs that coexpresses AR and FOXA1 (the AR+ FOXA1+ CK14– subgroup) could represent the equivalent of the luminal-AR subgroup of human triple-negative breast cancer.

Conclusions

We identified an AR+ FOXA1+ CK14– subgroup of triple-negative FMCs that might correspond to the luminal-AR subgroup of human triple-negative breast cancers. Cats with FMC may be interesting spontaneous animal models to investigate new strategies targeting the androgen receptor, especially in the aggressive subtype of AR+ basal-like triple-negative mammary carcinomas with loss of FOXA1 expression (the AR+ FOXA1–CK14+ subgroup).

Comprehensive analysis of the GATA transcription factor gene family in breast carcinoma using gene microarrays, online databases and integrated bioinformatics

Integrated studies of accumulated data can be performed to obtain more reliable information and more feasible measures for investigating the potential diagnostic and prognostic biomarkers of breast cancer and exploring related molecular mechanisms. Our study aimed to explore the GATA family members involved in breast cancer by integrating data from The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO) and other online databases. We performed an integrated analysis of published studies from GEO and analyzed clinical data from TCGA and GTEx to evaluate the clinical significance and prognosis values of the GATA family in breast cancer. GATA3 was found to be upregulated and exhibited a favorable value in the diagnosis and prognosis of breast cancer. Through this study, we identified possible GATA3-correlated genes and core pathways that play an important role, which requires further investigation in breast cancer.

Germline Genetic Variants in GATA3 and Breast Cancer Treatment Outcomes in SWOG S8897 Trial and the Pathways Study

INTRODUCTION

GATA3 is a critical transcription factor in maintaining the differentiated state of luminal mammary epithelial cells. We sought to determine the prognostic and predictive roles of GATA3 genotypes for breast cancer.

PATIENTS AND METHODS

Twelve single nucleotide polymorphisms (SNPs) were genotyped in 2 breast cancer cohorts, including the SWOG S8897 trial where patients were treated with adjuvant chemotherapy (CAF [cyclophosphamide, doxorubicin, 5-fluorouracil] vs. CMF [cyclophosphamide, methotrexate, 5-fluorouracil]) or untreated, and the observational Pathways Study.

RESULTS

In the S8897 trial, rs3802604 and rs568727 were associated with disease-free survival and overall survival in the treated group, regardless of chemotherapy regimen. The GG genotype of rs3802604 conferred poorer overall survival (adjusted hazard ratio, 2.45; 95% confidence interval, 1.48-4.05) and disease-free survival (adjusted hazard ratio, 1.95; 95% confidence interval, 1.27-2.99) compared with the AA genotype. Similar associations were found for rs568727. In contrast, no association with either SNP was found in the untreated group. Subgroup analyses indicated that these 2 SNPs more strongly influenced outcomes in the patients who also received tamoxifen. However, the associations in the subgroup with tamoxifen treatment were not replicated in the Pathways Study, possibly owing to substantial differences between the 2 patient cohorts, such as chemotherapy regimen and length of follow-up. Results from joint analyses across these 2 cohorts were marginally significant, driven by the results in S8897. Bioinformatic analyses support potential functional disruption of the GATA3 SNPs in breast tissue.

CONCLUSIONS

The present study provides some evidence for the predictive value of GATA3 genotypes for breast cancer adjuvant therapies. Future replication studies in appropriate patient populations are warranted.

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.

A DNA Methylation-Based Test for Breast Cancer Detection in Circulating Cell-Free DNA

BACKGROUND

Breast cancer (BrC) is the most frequent neoplasm in women. New biomarkers, including aberrant DNA methylation, may improve BrC management. Herein, we evaluated the detection and prognostic performance of seven genes' promoter methylation (APC, BRCA1, CCND2, FOXA1, PSAT1, RASSF1A and SCGB3A1).

METHODS

Methylation levels were assessed in primary BrC tissues by quantitative methylation-specific polymerase chain reaction (QMSP) and in circulating cell-free DNA (ccfDNA) by multiplex QMSP from two independent cohorts of patients (Cohort #1, n = 137; and Cohort #2, n = 44). Receiver operating characteristic (ROC) curves were constructed, and log-rank test and Cox regression were performed to assess the prognostic value of genes' methylation levels.

RESULTS

The gene-panel APC, FOXA1, RASSF1A, SCGB3A1 discriminated normal from cancerous tissue with high accuracy (95.55%). In multivariable analysis, high PSAT1-methylation levels [>percentile 75 (P75)] associated with longer disease-free survival, whereas higher FOXA1-methylation levels (>P75) associated with shorter disease-specific survival. The best performing panel in ccfDNA (APC, FOXA1 and RASSF1A) disclosed a sensitivity, specificity and accuracy over 70%.

CONCLUSIONS

This approach enables BrC accurate diagnosis and prognostic stratification in tissue samples, and allows for early detection in liquid biopsies, thus suggesting a putative value for patient management.

Regulation of the oncogenic phenotype by the nuclear body protein ZC3H8

BACKGROUND

The Zc3h8 gene encodes a protein with three zinc finger motifs in the C-terminal region. The protein has been identified as a component of the Little Elongation Complex, involved in transcription of small nuclear RNAs. ZC3H8 is overexpressed in a number of human and mouse breast cancer cell lines, and elevated mRNA levels are associated with a poorer prognosis for women with breast cancer.

METHODS

We used RNA silencing to decrease levels of expression in mouse mammary tumor cells and overexpression of ZC3H8 in cells derived from the normal mouse mammary gland. We measured characteristics of cell behavior in vitro, including proliferation, migration, invasion, growth in soft agar, and spheroid growth. We assessed the ability of these cells to form tumors in syngeneic BALB/c mice. ZC3H8 protein was visualized in cells using confocal microscopy.

RESULTS

Tumor cells with lower ZC3H8 expression exhibited decreased proliferation rates, slower migration, reduced ability to invade through a basement membrane, and decreased anchorage independent growth in vitro. Cells with lower ZC3H8 levels formed fewer and smaller tumors in animals. Overexpression of ZC3H8 in non-tumorigenic COMMA-D cells led to an opposite effect. ZC3H8 protein localized to both PML bodies and Cajal bodies within the nucleus. ZC3H8 has a casein kinase 2 (CK2) phosphorylation site near the N-terminus, and a CK2 inhibitor caused the numerous PML bodies and ZC3H8 to coalesce to a few larger bodies. Removal of the inhibitor restored PML bodies to their original state. A mutant ZC3H8 lacking the predicted CK2 phosphorylation site showed localization and numbers of ZC3H8/PML bodies similar to wild type. In contrast, a mutant constructed with a glutamic acid in place of the phosphorylatable threonine showed dramatically increased numbers of smaller nuclear foci.

CONCLUSIONS

These experiments demonstrate that Zc3h8 expression contributes to aggressive tumor cell behavior in vitro and in vivo. Our studies show that ZC3H8 integrity is key to maintenance of PML bodies. The work provides a link between the Little Elongation Complex, PML bodies, and the cancer cell phenotype.

The impact of the effectiveness of GATA3 as a prognostic factor in breast cancer

The transcription factor GATA3 plays a significant role in mammary gland development and differentiation. We analyzed expression of GATA3 in breast cancer (BC) cell lines and clinical specimens from BC patients in Taiwan. Semiquantitative reverse-transcription polymerase chain reaction (RT-PCR), quantitative real-time PCR were carried out to determine the mRNA level of GATA3 from 241 pairs of matched tumor and adjacent normal tissues from anonymous female donors. GATA3 immunohistochemistry (IHC) staining and H-score were performed (n = 25). Inducing and silencing of GATA3 were done by exposure MCF-7 cell line to nicotine or curcumin, respectively. GATA3 expression was detected in most of the estrogen receptor-positive (ER+) tumor specimens (176/241, 73%) compared with paired normal tissues (65/241, 27%) (P < .001). The GATA3 level was highest in Luminal A, and independent t-tests revealed higher GATA3 was associated with ER+ (P = .018) and BC stages (stage II, and stage IV). Nuclear protein expression of GATA3 was detected in tumor tissues (P < .001) with higher H-score in Luminal A patients (P = .012). Kaplan-Meier survival analyses showed that ER+/progesterone receptor (PgR)+ and lower grade BC patients with relatively high GATA3 had better clinical overall survival (OS). GATA3 regulates ERα and BCL-2 as BC luminal subtype markers. Cox univariate and multivariate analyses demonstrated that the expression of GATA3 was an effective predictor of the risk of death. We demonstrated a correlation between GATA3 expression and only ER+ and suggest that a higher GATA3 expression is a good prognostic factor for OS for ER+ BC patients.

Immunohistochemical Analysis of the Expression of Breast Markers in Basal-like Breast Carcinomas Defined as Triple Negative Cancers Expressing Keratin 5

Estrogen and progesterone receptors are possible markers for suggesting a mammary origin of metastatic carcinoma, but are useless in cases of triple negative breast cancers (TNBC). Five other potential markers of breast origin were investigated on tissue microarrays in a series of TNBCs showing keratin 5 expression, consistent with a basal-like phenotype. GATA-3 staining was observed in 82 of 115 triple negative cases (71.3%) including 23 cases with >5% staining. Mammaglobin staining was detected in 30 cases (26.0%) including 12 with >5% staining. GCDFP-15 was seen in 23 cases (20.0%) including 9 with >5% staining. NY-BR-1 positivity was present in 7 cases (6.0%) including 3 patients with >5% staining. BCA-225 staining was observed in 74 cases (64.3%); however this latter marker lacks also specificity owing to the reported widespread staining in other malignancies. GATA-3, mammaglobin and GCDFP-15 coexpression was seen in one case (0.9%), whereas GATA-3 and mammaglobin or mammaglobin and GCDFP-15 coexpression was present in 2 and 2 cases (1.7%), respectively. Using at least 5% staining as cut-off, the expression of any of the last 4 markers was 34.7%. The expression of GATA-3, mammaglobin, GCDFP-15 and NY-BR-1 is lower in TNBC-s than in breast carcinomas in general, and this may be even lower in basal-like carcinomas. Although these markers are not fully specific, by using them, a subset of basal-like TNBC-s can be identified as of mammary origin. However, a substantial proportion will not show any staining with any of these markers.

A functional variant rs4442975 modulating FOXA1-binding affinity does not influence the risk or progression of breast cancer in Chinese Han population

The DNA-binding protein FOXA1 has been shown to regulate nearly all estrogen receptor-chromatin interactions, thereby influencing target gene expression levels in breast cancer (BC) cells. Recently, the rs4442975 T-allele, which disrupts the recruitment of FOXA1 and interacts with the IGFBP5 promoter, was associated to BC susceptibility in a European population. We conducted a hospital-based case-control study that included 1227 cases and 1285 controls to explore the potential association between rs4442975 and BC risk in Chinese Han population, and the effect of this SNP on BC progression was also observed in cases. No significant associations between rs4442975 and BC risk were observed under any genetic models, with an odds ratio of 0.96 (95% confidence interval = 0.81-1.15) under the additive model. When stratified based on estrogen or progesterone receptor expression, smoking or drinking habits, or menopausal status, similar negative associations were observed for all subgroups. No significant associations were observed between rs4442975 and traditional progression factors such as tumor size, nodal status, distant metastasis, or TNM staging. These results reveal that rs4442975 may not confer a risk of BC occurrence or progression in the Chinese Han population, which indicates a distinct association related to genetic heterogeneity across ethnic populations.

Role of transcription factor FOXA1 in non‑small cell lung cancer

In our previous study, stable subpopulations of the A549 lung cancer cell line with high/low invasive potential (H/L‑INV) were obtained. In the present study, microarray analysis of the H/L‑INV A549 subpopulations was performed to evaluate genes associated with high invasiveness. Forkhead box protein A1 (FOXA1) was selected for further investigation. The expression levels of FOXA1 in the primary lesion and metastatic lymph nodes were assessed using reverse transcription‑quantitative polymerase chain reaction (RT‑qPCR) analysis. In addition, the mRNA and protein expression levels of FOXA1 were examined in H‑INV A549 cells transfected with a specific FOXA1 small interfering RNA (siRNA), and the role of FOXA1 in the proliferation, invasion and metastasis of non‑small cell lung cancer (NSCLC) cells was evaluated. FOXA1 was overexpressed in metastatic lymph nodes, compared with its expression in NSCLC primary tumours. The results of western blot and RT‑qPCR analyses confirmed that FOXA1 siRNA transfection led to a decrease in the expression of FOXA1 in H‑INV A549 cells. FOXA1 siRNA transfection caused G0/G1 phase cell cycle arrest, and also reduced the invasion, migration and proliferation abilities of the H‑INV A549 cells. In conclusion, the results of the present study suggested that FOXA1 is a potential oncogene in NSCLC; therefore, specific interference of the expression of FOXA1 may represent a novel approach for the treatment of NSCLC.

Discovery and validation of information theory-based transcription factor and cofactor binding site motifs

Data from ChIP-seq experiments can derive the genome-wide binding specificities of transcription factors (TFs) and other regulatory proteins. We analyzed 765 ENCODE ChIP-seq peak datasets of 207 human TFs with a novel motif discovery pipeline based on recursive, thresholded entropy minimization. This approach, while obviating the need to compensate for skewed nucleotide composition, distinguishes true binding motifs from noise, quantifies the strengths of individual binding sites based on computed affinity and detects adjacent cofactor binding sites that coordinate with the targets of primary, immunoprecipitated TFs. We obtained contiguous and bipartite information theory-based position weight matrices (iPWMs) for 93 sequence-specific TFs, discovered 23 cofactor motifs for 127 TFs and revealed six high-confidence novel motifs. The reliability and accuracy of these iPWMs were determined via four independent validation methods, including the detection of experimentally proven binding sites, explanation of effects of characterized SNPs, comparison with previously published motifs and statistical analyses. We also predict previously unreported TF coregulatory interactions (e.g. TF complexes). These iPWMs constitute a powerful tool for predicting the effects of sequence variants in known binding sites, performing mutation analysis on regulatory SNPs and predicting previously unrecognized binding sites and target genes.

A novel FOXA1/ESR1 interacting pathway: A study of Oncomine™ breast cancer microarrays

Forkhead box protein A1 (FOXA1) is essential for the growth and differentiation of breast epithelium, and has a favorable outcome in breast cancer (BC). Elevated FOXA1 expression in BC also facilitates hormone responsiveness in estrogen receptor (ESR)-positive BC. However, the interaction between these two pathways is not fully understood. FOXA1 and GATA binding protein 3 (GATA3) along with ESR1 expression are responsible for maintaining a luminal phenotype, thus suggesting the existence of a strong association between them. The present study utilized the Oncomine™ microarray database to identify FOXA1:ESR1 and FOXA1:ESR1:GATA3 co-expression co-regulated genes. Oncomine™ analysis revealed 115 and 79 overlapping genes clusters in FOXA1:ESR1 and FOXA1:ESR1:GATA3 microarrays, respectively. Five ESR1 direct target genes [trefoil factor 1 (TFF1/PS2), B-cell lymphoma 2 (BCL2), seven in absentia homolog 2 (SIAH2), cellular myeloblastosis viral oncogene homolog (CMYB) and progesterone receptor (PGR)] were detected in the co-expression clusters. To further investigate the role of FOXA1 in ESR1-positive cells, MCF7 cells were transfected with a FOXA1 expression plasmid, and it was observed that the direct target genes of ESR1 (PS2, BCL2, SIAH2 and PGR) were significantly regulated upon transfection. Analysis of one of these target genes, PS2, revealed the presence of two FOXA1 binding sites in the vicinity of the estrogen response element (ERE), which was confirmed by binding assays. Under estrogen stimulation, FOXA1 protein was recruited to the FOXA1 site and could also bind to the ERE site (although in minimal amounts) in the PS2 promoter. Co-transfection of FOXA1/ESR1 expression plasmids demonstrated a significantly regulation of the target genes identified in the FOXA1/ESR1 multi-arrays compared with only FOXA1 transfection, which was suggestive of a synergistic effect of ESR1 and FOXA1 on the target genes. In summary, the present study identified novel FOXA1, ESR1 and GATA3 co-expressed genes that may be involved in breast tumorigenesis.

Characterisation of GATA3 expression in invasive breast cancer: differences in histological subtypes and immunohistochemically defined molecular subtypes

AIMS

GATA-binding protein 3 (GATA3) is a sensitive and relatively specific marker in breast and urothelial carcinomas. Its diagnostic utility in primary and metastatic breast cancers has been explored and confirmed. However, the relationship between GATA3 expression and different breast carcinoma intrinsic subtypes has not been specifically defined in the literature despite a few reports with a small number of cases. The aim of the current investigation is to clarify GATA3 expression among different histological subtypes and surrogate molecular breast carcinoma subtypes in a large series of cases.

METHODS

Immunohistochemical staining of GATA3, GCDFP15 and mammaglobin was performed in a cohort of 1637 cases of primary invasive breast carcinoma. The association of GATA3 expression with different histological and surrogate intrinsic subtypes was assessed and compared with the expression of GCDFP15 and mammaglobin.

RESULTS

The overall positivity of GATA3 across the various immunohistochemistry-based surrogate intrinsic subtypes was 99.51% for luminal A-like, 97.70% for luminal B-like, 68.50% for HER2 overexpression and 20.16% for triple negative breast cancers. GATA3 expression was positively correlated with estrogen receptor (ER)-positive (luminal subtypes) breast carcinomas. For luminal-like and HER2 overexpression subtypes, GATA3 was much more sensitive than GCDFP15 and mammaglobin. For triple negative tumours, GATA3 was less sensitive than GCDFP15.

CONCLUSIONS

GATA3 exhibits a relatively high sensitivity for breast carcinomas. It is more sensitive than GCDFP15 and mammaglobin in luminal-like and HER2 overexpression subtypes. GATA3 expression is associated with breast carcinomas of luminal subtype and low histological grade.

Prognostic and clinicopathological value of GATA binding protein 3 in breast cancer: A systematic review and meta-analysis

The potential prognostic value of GATA binding protein 3 (GATA3) in breast cancer has recently increased, although the evidence is inconclusive. This meta-analysis of 10 articles involving 5,080 breast cancer patients explored the prognostic and clinicopathological value of GATA3 in breast cancer. Time to tumor progression (TTP) and overall survival (OS) were primary endpoints. Pooled hazard ratio (HR), pooled risk ratio (RR), and 95% confidence interval (CI) were calculated to evaluate the association between GATA3, prognosis, and clinicopathological parameters. High GATA3 expression predicts breast cancer, with a HR (HR = 0.671; 95% CI = 0.475–0.947; P = 0.023) of TTP, but is not associated with OS (HR = 0.889; 95% CI = 0.789–1.001; P = 0.052). GATA3 overexpression is associated with positive ER (RR = 3.155; 95% CI = 1.680–5.923; P = 0.000), positive PR (RR = 3.949; 95% CI = 1.567–9.954, P = 0.004), lower nuclear grade (RR = 0.435; 95% CI = 0.369–0.514; P = 0.000), and smaller tumor size (RR = 0.816; 95% CI = 0.709–0.940; P = 0.005). High GATA3 expression may predict TTP in breast cancer, and such patients may show better clinicopathological features.

Plk1 Regulates the Repressor Function of FoxM1b by inhibiting its Interaction with the Retinoblastoma Protein

FoxM1b is a cell cycle-regulated transcription factor, whose over-expression is a marker for poor outcome in cancers. Its transcriptional activation function requires phosphorylation by Cdk1 or Cdk2 that primes FoxM1b for phosphorylation by Plk1, which triggers association with the co-activator CBP. FoxM1b also possesses transcriptional repression function. It represses the mammary differentiation gene GATA3 involving DNMT3b and Rb. We investigated what determines the two distinct functions of FoxM1b: activation and repression. We show that Rb binds to the C-terminal activation domain of FoxM1b. Analyses with phospho-defective and phospho-mimetic mutants of FoxM1b identified a critical role of the Plk1 phosphorylation sites in regulating the binding of FoxM1b to Rb and DNMT3b. That is opposite of what was seen for the interaction of FoxM1b with CBP. We show that, in addition to GATA3, FoxM1b also represses the mammary luminal differentiation marker FoxA1 by promoter-methylation, and that is regulated by the Plk1 phosphorylation sites in FoxM1b. Our results show that the Plk1 phosphorylation sites in FoxM1b serve as a regulator for its repressor function, and they provide insights into how FoxM1b inhibits differentiation genes and activates proliferation genes during cancer progression.

Characterisation of male breast cancer: a descriptive biomarker study from a large patient series

Male breast cancer (MBC) is rare. We assembled 446 MBCs on tissue microarrays and assessed clinicopathological information, together with data from 15 published studies, totalling 1984 cases. By immunohistochemistry we investigated 14 biomarkers (ERα, ERβ1, ERβ2, ERβ5, PR, AR, Bcl-2, HER2, p53, E-cadherin, Ki67, survivin, prolactin, FOXA1) for survival impact. The main histological subtype in our cohort and combined analyses was ductal (81%, 83%), grade 2; (40%, 44%), respectively. Cases were predominantly ERα (84%, 82%) and PR positive (74%, 71%), respectively, with HER2 expression being infrequent (2%, 10%), respectively. In our cohort, advanced age (>67) was the strongest predictor of overall (OS) and disease free survival (DFS) (p = 0.00001; p = 0.01, respectively). Node positivity negatively impacted DFS (p = 0.04). FOXA1 p = 0.005) and AR p = 0.009) were both positively prognostic for DFS, remaining upon multivariate analysis. Network analysis showed ERα, AR and FOXA1 significantly correlated. In summary, the principle phenotype of MBC was luminal A, ductal, grade 2. In ERα+ MBC, only AR had prognostic significance, suggesting AR blockade could be employed therapeutically.

Data Mining for Identification of Molecular Targets in Ovarian Cancer

Ovarian cancer is possibly the sixth most common malignancy worldwide, in Mexico representing the fourth leading cause of gynecological cancer death more than 70% being diagnosed at an advanced stage and the survival being very poor. Ovarian tumors are classified according to histological characteristics, epithelial ovarian cancer as the most common (~80%). We here used high-density microarrays and a systems biology approach to identify tissue-associated deregulated genes. Non-malignant ovarian tumors showed a gene expression profile associated with immune mediated inflammatory responses (28 genes), whereas malignant tumors had a gene expression profile related to cell cycle regulation (1,329 genes) and ovarian cell lines to cell cycling and metabolism (1,664 genes).

GATA3 expression and its relationship with clinicopathological parameters in invasive breast carcinomas

GATA3, as a transcription factor, is associated with estrogen receptor (ER) expression and necessary for luminal cell differentiation in mammary glands. Association of GATA3 expression with clinicopathological parameters, molecular subtypes of tumors, disease free survival (DFS) and overall survival (OS) for breast carcinoma patients were evaluated in this study. We immunohistochemically stained GATA3, CK5/6, EGFR, CK14 and vimentin on tissue microarray blocks of 457 invasive breast carcinomas. Tumors are sub-classified as luminal A, luminal B, HER2 expressing, basal-like and null type according to their hormonal status with cerbB2, CK 5/6 and EGFR expressions. Follow-up data for 254 cases were obtained. 215/457 (47%) tumors were GATA3 positive. GATA3 expression was inversely correlated with mitotic count (p<0.0001), nuclear grade (p=0.001), histological grade (p=0.001), tumor necrosis (p=0.001), stromal lymphocytic response (p<0.01), nipple invasion (p=0.01), metastasis (p=0.03), vimentin (p=0.0003), EGFR (p=0.015) and CK14 (p=0.001) expressions; and directly associated with ER (p<0.0001) and progesterone receptor (PR) (p<0.0001) expressions. Luminal A carcinomas had the highest frequency for GATA-3 (140/245), however basal-like carcinomas had the lowest (1/42) (p<0.0001). None of the medullary and metaplastic carcinomas expressed GATA3. GATA3 was associated with good DFS and OS (p=0.001 and p=0.0009) and was an independent prognostic factor for OS. GATA3 expression, regardless of the subtype, may have a prognostic significance for breast carcinomas through its ability to promote the differentiation of luminal progenitor cells.

Identifying biological mechanisms for favorable cancer prognosis using non-hypothesis-driven iterative survival analysis

Survival analyses based on the Kaplan-Meier estimate have been pervasively used to support or validate the relevance of biological mechanisms in cancer research. Recently, with the appearance of gene expression high-throughput technologies, this kind of analysis has been applied to tumor transcriptomics data. In a 'bottom-up' approach, gene-expression profiles that are associated with a deregulated pathway hypothetically involved in cancer progression are first identified and then subsequently correlated with a survival effect, which statistically supports or requires the rejection of such a hypothesis. In this work, we propose a 'top-down' approach, in which the clinical outcome (survival) is the starting point that guides the identification of deregulated biological mechanisms in cancer by a non-hypothesis-driven iterative survival analysis. We show that the application of our novel method to a population of ~2,000 breast cancer patients of the METABRIC consortium allows the identification of several well-known cancer mechanisms, such as ERBB4, HNF3A and TGFB pathways, and the investigation of their paradoxical dual effect. In addition, several novel biological mechanisms are proposed as potentially involved in cancer progression. The proposed exploratory methodology can be considered both alternative and complementary to classical 'bottom-up' approaches for validation of biological hypotheses. We propose that our method may be used to better characterize cancer, and may therefore impact the future design of therapies that are truly molecularly tailored to individual patients. The method, named SURCOMED, was implemented as a web-based tool, which is publicly available at http://surcomed.vital-it.ch. R scripts are also available at http://surcomed.sourceforge.net).

Integrated computational approach to the analysis of RNA-seq data reveals new transcriptional regulators of psoriasis

Psoriasis is a common inflammatory skin disease with complex etiology and chronic progression. To provide novel insights into the regulatory molecular mechanisms of the disease, we performed RNA sequencing analysis of 14 pairs of skin samples collected from patients with psoriasis. Subsequent pathway analysis and extraction of the transcriptional regulators governing psoriasis-associated pathways was executed using a combination of the MetaCore Interactome enrichment tool and the cisExpress algorithm, followed by comparison to a set of previously described psoriasis response elements. A comparative approach allowed us to identify 42 core transcriptional regulators of the disease associated with inflammation (NFκB, IRF9, JUN, FOS, SRF), the activity of T cells in psoriatic lesions (STAT6, FOXP3, NFATC2, GATA3, TCF7, RUNX1), the hyperproliferation and migration of keratinocytes (JUN, FOS, NFIB, TFAP2A, TFAP2C) and lipid metabolism (TFAP2, RARA, VDR). In addition to the core regulators, we identified 38 transcription factors previously not associated with the disease that can clarify the pathogenesis of psoriasis. To illustrate these findings, we analyzed the regulatory role of one of the identified transcription factors (TFs), FOXA1. Using ChIP-seq and RNA-seq data, we concluded that the atypical expression of the FOXA1 TF is an important player in the disease as it inhibits the maturation of naive T cells into the (CD4+FOXA1+CD47+CD69+PD-L1(hi)FOXP3-) regulatory T cell subpopulation, therefore contributing to the development of psoriatic skin lesions.

GATA3 expression in clinically useful groups of breast carcinoma: a comparison with GCDFP15 and mammaglobin for identifying paired primary and metastatic tumors

GATA3 has been recognized as the novel marker for identifying primary and metastatic breast carcinomas, consistently showing that GATA3 was significantly more sensitive than traditional markers gross cystic disease fluid protein 15 (GCDFP15) and mammaglobin (MGB). However, clinically useful groups of breast carcinomas status were not identified, which were determining appropriate treatment strategy, affecting the prognosis. In this study, we undertook a comparative study of the marker GATA3 and GCDFP15 and MGB in clinically useful groups of paired primary and metastatic breast cancer. We retrieved 64 cases of matched primary and metastatic breast cancer from the surgical pathology archive at our institution. According to the emerging 2015 St. Gallen Consensus, the clinically useful groups were divided into ER and/or PR (+), HER2 (-), abbreviated as A; ER and/or PR (+), HER2 (+), abbreviated as B; ER and PR (-), HER2 (+), abbreviated as C; ER, PR and HER2 (-), abbreviated as D; each group contained 16 cases (n=16). Tissue microarrays were created, with three 1-mm punch specimens from each case. The tissue microarrays were cut at 4-μm thickness and stained with monoclonal antibodies to GATA3, GCDFP15, and MGB. Staining intensity (0-3+) and extent (0%-100%) were scored with an H-score calculated (range, 0-300). Sensitivities by varying H-score cutoffs (any; ≥50; ≥150) for a positive result in the clinically useful groups of matched primary or metastatic breast cancer among GATA3, GCDFP15, and MGB. GATA3 was significantly more sensitive than GCDFP15 and MGB A and B groups (P<.05) rather than C and D groups (P>.05). However, GATA3 in conjunction with GCDFP15 and MGB detection could improve the sensitivity of C group (P<.05) rather than D group (P>.05). Significantly, good coincidence was observed between primary and metastatic tumor GATA3 expression (κ value = 0.826 >0.75) as compared with the coincidence of GCDFP15 (κ value =0.492 <0.75) and MGB (κ value =0.593 <0.75) (both P<.05). In conclusion, GATA3 expression did not show the same sensitivity for the clinically useful groups of breast cancer. GATA3 expression is positively correlated with ER-positive, PR-positive, and HER2-positive carcinomas. In addition, the matched primary and metastatic tumor expression of GATA3 shows good coincidence. We propose the careful selection of GATA3 for identifying hormone receptor negativity of breast cancer, especially in the case of triple-negative breast cancer.

Genomic Changes in Normal Breast Tissue in Women at Normal Risk or at High Risk for Breast Cancer

Sporadic breast cancer develops through the accumulation of molecular abnormalities in normal breast tissue, resulting from exposure to estrogens and other carcinogens beginning at adolescence and continuing throughout life. These molecular changes may take a variety of forms, including numerical and structural chromosomal abnormalities, epigenetic changes, and gene expression alterations. To characterize these abnormalities, a review of the literature has been conducted to define the molecular changes in each of the above major genomic categories in normal breast tissue considered to be either at normal risk or at high risk for sporadic breast cancer. This review indicates that normal risk breast tissues (such as reduction mammoplasty) contain evidence of early breast carcinogenesis including loss of heterozygosity, DNA methylation of tumor suppressor and other genes, and telomere shortening. In normal tissues at high risk for breast cancer (such as normal breast tissue adjacent to breast cancer or the contralateral breast), these changes persist, and are increased and accompanied by aneuploidy, increased genomic instability, a wide range of gene expression differences, development of large cancerized fields, and increased proliferation. These changes are consistent with early and long-standing exposure to carcinogens, especially estrogens. A model for the breast carcinogenic pathway in normal risk and high-risk breast tissues is proposed. These findings should clarify our understanding of breast carcinogenesis in normal breast tissue and promote development of improved methods for risk assessment and breast cancer prevention in women.

Coexpression of androgen receptor and FOXA1 in nonmetastatic triple-negative breast cancer: ancillary study from PACS08 trial

AIM

Microarray studies identified a subgroup of molecular apocrine tumors (estrogen receptor [ER] negative/androgen receptor [AR] positive) that express luminal genes including FOXA1. FOXA1 may direct AR to sites normally occupied by ER in luminal tumors, inducing an estrogen-like gene program that stimulated proliferation.

MATERIALS & METHODS

Expression of AR and FOXA1 was evaluated by immunohistochemistry in 592 patients with nonmetastatic triple-negative breast cancer (TNBC).

RESULTS

Coexpression of AR and FOXA1 was found in 15.2% of patients. These tumors were more frequently lobular, found in older patients and exhibited a lower nuclear grade and a greater degree of node involvement. They less often exhibited lymphocytic infiltrate, pushing margins, syncytial architecture, central fibrosis or necrosis.

CONCLUSION

TNBC with coexpression of AR and FOXA1 seems to behave like luminal tumors with a morphological profile distinct from other TNBC. These biomarkers could be useful to identify a subgroup of TNBC and could have future therapeutic implications.

Clinicopathological significance of forkhead box protein A1 in breast cancer: A meta-analysis

The aim of the present study was to investigate the associations between the expression of forkhead box protein A1 (FOXA1) and differential clinicopathological characteristics in breast cancer using a meta-analysis method. Eligible studies that investigated the correlation between FOXA1 expression and the clinical characteristics of breast cancer were collected through searching numerous databases, including PubMed, EMBASE, the Chinese National Knowledge Infrastructure and the VIP database. In total, eight studies were included in the meta-analysis. Following a systematic analysis, the expression of FOXA1 was found to be significantly associated with the estrogen receptor α status, the progesterone receptor status, lymph node metastasis and the histological grade in breast cancer. However, no statistically significant association was observed between FOXA1 expression and the human epidermal growth factor receptor-2 status in breast cancer patients.