Cancer-specific mutation of GATA3 disrupts the transcriptional regulatory network governed by Estrogen Receptor alpha, FOXA1 and GATA3

Unravelling the therapeutic potential of forkhead box proteins in breast cancer: An update (Review)

Breast cancer, a prominent cause of mortality among women, develops from abnormal growth of breast tissue, thereby rendering it one of the most commonly detected cancers in the female population. Although numerous treatment strategies are available for breast cancer, discordance in terms of effective treatment and response still exists. Recently, the potential of signaling pathways and transcription factors has gained substantial attention in the cancer community; therefore, understanding their role will assist researchers in comprehending the onset and advancement of breast cancer. Forkhead box (FOX) proteins, which are important transcription factors, are considered crucial regulators of various cellular activities, including cell division and proliferation. The present study explored several subclasses of FOX proteins and their possible role in breast carcinogenesis, followed by the interaction between microRNA (miRNA) and FOX proteins. This interaction is implicated in promoting cell infiltration into surrounding tissues, ultimately leading to metastasis. The various roles that FOX proteins play in breast cancer development, their intricate relationships with miRNA, and their involvement in therapeutic resistance highlight the complexity of breast cancer dynamics. Therefore, recognizing the progress and challenges in current treatments is crucial because, despite advancements, persistent disparities in treatment effectiveness underscore the need for ongoing research, with future studies emphasizing the necessity for targeted strategies that account for the multifaceted aspects of breast cancer.

Progestogen-driven B7-H4 contributes to onco-fetal immune tolerance

Immune tolerance mechanisms are shared in cancer and pregnancy. Through cross-analyzing single-cell RNA-sequencing data from multiple human cancer types and the maternal-fetal interface, we found B7-H4 (VTCN1) is an onco-fetal immune tolerance checkpoint. We showed that genetic deficiency of B7-H4 resulted in immune activation and fetal resorption in allogeneic pregnancy models. Analogously, B7-H4 contributed to MPA/DMBA-induced breast cancer progression, accompanied by CD8+ T cell exhaustion. Female hormone screening revealed that progesterone stimulated B7-H4 expression in placental and breast cancer cells. Mechanistically, progesterone receptor (PR) bound to a newly identified -58 kb enhancer, thereby mediating B7-H4 transcription via the PR-P300-BRD4 axis. PR antagonist or BRD4 degrader potentiated immunotherapy in a murine B7-H4+ breast cancer model. Thus, our work unravels a mechanistic and biological connection of a female sex hormone (progesterone) to onco-fetal immune tolerance via B7-H4 and suggests that the PR-P300-BRD4 axis is targetable for treating B7-H4+ cancer.

Genomic transcription factor binding site selection is edited by the chromatin remodeling factor CHD4

Biologically precise enhancer licensing by lineage-determining transcription factors enables activation of transcripts appropriate to biological demand and prevents deleterious gene activation. This essential process is challenged by the millions of matches to most transcription factor binding motifs present in many eukaryotic genomes, leading to questions about how transcription factors achieve the exquisite specificity required. The importance of chromatin remodeling factors to enhancer activation is highlighted by their frequent mutation in developmental disorders and in cancer. Here, we determine the roles of CHD4 in enhancer licensing and maintenance in breast cancer cells and during cellular reprogramming. In unchallenged basal breast cancer cells, CHD4 modulates chromatin accessibility. Its depletion leads to redistribution of transcription factors to previously unoccupied sites. During cellular reprogramming induced by the pioneer factor GATA3, CHD4 activity is necessary to prevent inappropriate chromatin opening. Mechanistically, CHD4 promotes nucleosome positioning over GATA3 binding motifs to compete with transcription factor-DNA interaction. We propose that CHD4 acts as a chromatin proof-reading enzyme that prevents unnecessary gene expression by editing chromatin binding activities of transcription factors.

Redirecting the pioneering function of FOXA1 with covalent small molecules

Pioneer transcription factors (TFs) exhibit a specialized ability to bind to and open closed chromatin, facilitating engagement by other regulatory factors involved in gene activation or repression. Chemical probes are lacking for pioneer TFs, which has hindered their mechanistic investigation in cells. Here, we report the chemical proteomic discovery of electrophilic small molecules that stereoselectively and site-specifically bind the pioneer TF, FOXA1, at a cysteine (C258) within the forkhead DNA-binding domain. We show that these covalent ligands react with FOXA1 in a DNA-dependent manner and rapidly remodel its pioneer activity in prostate cancer cells reflected in redistribution of FOXA1 binding across the genome and directionally correlated changes in chromatin accessibility. Motif analysis supports a mechanism where the covalent ligands relax the canonical DNA binding preference of FOXA1 by strengthening interactions with suboptimal ancillary sequences in predicted proximity to C258. Our findings reveal a striking plasticity underpinning the pioneering function of FOXA1 that can be controlled by small molecules.

Pioneer factors: roles and their regulation in development

Pioneer factors are a subclass of transcription factors that can bind and initiate opening of silent chromatin regions. Pioneer factors subsequently regulate lineage-specific genes and enhancers and, thus, activate the zygotic genome after fertilization, guide cell fate transitions during development, and promote various forms of human cancers. As such, pioneer factors are useful in directed cell reprogramming. In this review, we define the structural and functional characteristics of pioneer factors, how they bind and initiate opening of closed chromatin regions, and the consequences for chromatin dynamics and gene expression during cell differentiation. We also discuss emerging mechanisms that modulate pioneer factors during development.

17-β Estradiol (E2) distinctly regulates the expression of IL-4 and IL-13 in Th2 cells via modulating the interplay between GATA3 and PU.1

17-β Estradiol (E2) has long standing known functions in regulating human physiology as well as immune system. E2 is known to elicit a protective role against experimental autoimmune encephalomyelitis (EAE) and has been used as a drug for treatment against multiple sclerosis. Moreover, E2 regulates the adaptive immune system by directly affecting the T helper cell subsets differentiation and antibody secretion mediated by B cells. Reports have shown that E2 promotes Th1 and Treg cell differentiation; whereas it attenuated the Th17 and Tfh cell differentiation. Albeit multiple and contrasting studies, the mechanisms of behind E2 action on Th2 cells remained understudied. Hence, we sought to dissect the impact of E2 in Th2 cell differentiation. In this study, we elucidated the molecular mechanisms behind E2-mediated regulation of the differentiation of Th2 cells. We observed that E2 significantly attenuated the IL-4-secreting Th2 population in an ERα-dependent manner. We validated these findings using ICI 182, 780, an antagonist to ERα, not ERβ and ectopically overexpressing ERα in Th2 cells. We further determined that ERα alters the recruitment of GATA3 and PU.1 to Il4 gene by directly interacting with them. This altered recruitment was observed to be stronger at Il4 than Il13 locus. Interestingly, we detected a distinct recruitment of GATA3 and PU.1 at Il13 gene; however, there was no E2-mediated broad alteration in the recruitment of histone-modifiers at Il13 locus. These findings suggest that E2 regulates Il4 in a distinctly separate mechanism as opposed to Il13 locus in Th2 cells.

Interpretable predictive models of genome-wide aryl hydrocarbon receptor-DNA binding reveal tissue-specific binding determinants

The aryl hydrocarbon receptor (AhR) is an inducible transcription factor whose ligands include the potent environmental contaminant 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Ligand-activated AhR binds to DNA at dioxin response elements (DREs) containing the core motif 5'-GCGTG-3'. However, AhR binding is highly tissue specific. Most DREs in accessible chromatin are not bound by TCDD-activated AhR, and DREs accessible in multiple tissues can be bound in some and unbound in others. As such, AhR functions similarly to many nuclear receptors. Given that AhR possesses a strong core motif, it is suited for a motif-centered analysis of its binding. We developed interpretable machine learning models predicting the AhR binding status of DREs in MCF-7, GM17212, and HepG2 cells, as well as primary human hepatocytes. Cross-tissue models predicting transcription factor (TF)-DNA binding generally perform poorly. However, reasons for the low performance remain unexplored. By interpreting the results of individual within-tissue models and by examining the features leading to low cross-tissue performance, we identified sequence and chromatin context patterns correlated with AhR binding. We conclude that AhR binding is driven by a complex interplay of tissue-agnostic DRE flanking DNA sequence and tissue-specific local chromatin context. Additionally, we demonstrate that interpretable machine learning models can provide novel and experimentally testable mechanistic insights into DNA binding by inducible TFs.

A comprehensive analysis of GATA3 expression in carcinomas of various origins with emphasis on lung carcinomas

GATA3 is a transcription factor involved in embryogenesis of multiple human tissues and in maintaining cell differentiation and tissue homeostasis in the adult organism. GATA3 is also involved in carcinogenesis and regarded as a sensitive marker for urothelial and breast carcinomas, albeit expression in carcinomas of non-breast/urothelial origin has been frequently reported. We sought to examine the extent and intensity of GATA3 expression in various carcinomas, mainly lung, urothelial, and breast and various other primary sites. Patients with breast carcinoma (N=40), carcinoma of the urinary bladder/renal pelvis (N=40), lung carcinoma (N=110) and various other origins (N=45) were included in the study. One hundred and sixty-five patients had a primary tumor diagnosis, and 70 cases had a metastatic tumor diagnosis. Our results showed that GATA3 expression was significantly more common in carcinomas of the breast, urinary bladder and renal pelvis compared to all other origins. All primary and 93% of metastatic urinary bladder carcinomas and 94% of the primary and 80% of metastatic breast carcinomas expressed GATA3. Expression was lower in non-urothelial histology of urinary primaries and in triple negative breast carcinomas. Focal staining, mostly faint, was seen in 5.6% of the primary lung adenocarcinomas and 35% of the primary lung squamous cell carcinomas. More extensive and intense staining was seen in 3.7% of the primary lung adenocarcinomas and 12% of the primary lung squamous cell carcinomas. Expression, mostly focal was also seen in 30% of the metastatic lung carcinomas. Finally, high expression was seen in 12.5% of the other tumors (one metastatic pancreatic carcinoma, one metastatic salivary gland adenocarcinoma NOS, one metastatic squamous cell carcinoma of the skin, one primary uterine cervix serous carcinoma, and one squamous cell carcinoma of the head and neck) and focal expression was present in another 22% of them. No ideal cut-off for positivity for GATA3 staining could be identified. In conclusion our study shows that GATA3 staining has two caveats in its use: the first is that non classical histologies of urothelial carcinomas and TNBC may be negative for the marker and secondly carcinomas of various origins may show (although rarely) intense positivity.

Different chromatin-scanning modes lead to targeting of compacted chromatin by pioneer factors FOXA1 and SOX2

Pioneer transcription factors interact with nucleosomes to scan silent, compact chromatin, enabling cooperative events that modulate gene activity. While at a subset of sites pioneer factors access chromatin by assisted loading with other transcription factors, the nucleosome-binding properties of pioneer factors enable them to initiate zygotic genome activation, embryonic development, and cellular reprogramming. To better understand nucleosome targeting in vivo, we assess whether pioneer factors FoxA1 and Sox2 target stable or unstable nucleosomes and find that they target DNase-resistant, stable nucleosomes, whereas HNF4A, a non-nucleosome binding factor, targets open, DNase-sensitive chromatin. Despite FOXA1 and SOX2 targeting similar proportions of DNase-resistant chromatin, using single-molecule tracking, we find that FOXA1 uses lower nucleoplasmic diffusion and longer residence times while SOX2 uses higher nucleoplasmic diffusion and shorter residence times to scan compact chromatin, while HNF4 scans compact chromatin much less efficiently. Thus, pioneer factors target compact chromatin through distinct processes.

The mutational landscape of a US Midwestern breast cancer cohort reveals subtype-specific cancer drivers and prognostic markers

BACKGROUND

Female breast cancer remains the second leading cause of cancer-related death in the USA. The heterogeneity in the tumor morphology across the cohort and within patients can lead to unpredictable therapy resistance, metastasis, and clinical outcome. Hence, supplementing classic pathological markers with intrinsic tumor molecular markers can help identify novel molecular subtypes and the discovery of actionable biomarkers.

METHODS

We conducted a large multi-institutional genomic analysis of paired normal and tumor samples from breast cancer patients to profile the complex genomic architecture of breast tumors. Long-term patient follow-up, therapeutic regimens, and treatment response for this cohort are documented using the Breast Cancer Collaborative Registry. The majority of the patients in this study were at tumor stage 1 (51.4%) and stage 2 (36.3%) at the time of diagnosis. Whole-exome sequencing data from 554 patients were used for mutational profiling and identifying cancer drivers.

RESULTS

We identified 54 tumors having at least 1000 mutations and 185 tumors with less than 100 mutations. Tumor mutational burden varied across the classified subtypes, and the top ten mutated genes include MUC4, MUC16, PIK3CA, TTN, TP53, NBPF10, NBPF1, CDC27, AHNAK2, and MUC2. Patients were classified based on seven biological and tumor-specific parameters, including grade, stage, hormone receptor status, histological subtype, Ki67 expression, lymph node status, race, and mutational profiles compared across different subtypes. Mutual exclusion of mutations in PIK3CA and TP53 was pronounced across different tumor grades. Cancer drivers specific to each subtype include TP53, PIK3CA, CDC27, CDH1, STK39, CBFB, MAP3K1, and GATA3, and mutations associated with patient survival were identified in our cohort.

CONCLUSIONS

This extensive study has revealed tumor burden, driver genes, co-occurrence, mutual exclusivity, and survival effects of mutations on a US Midwestern breast cancer cohort, paving the way for developing personalized therapeutic strategies.

GATA3 mediates doxorubicin resistance by inhibiting CYB5R2-catalyzed iron reduction in breast cancer cells

AIMS

Neoadjuvant chemotherapy (NAC) is the primary preoperative therapy for breast cancer. The luminal subtype of breast cancer shows less NAC response than the basal subtype, with an inefficient NAC treatment effect. Understanding of the molecular and cellular mechanisms responsible for this chemoresistance is an important issue when determining optimal treatment.

METHODS

Doxorubicin-induced apoptosis and ferroptosis was investigated using cytotoxicity, western blotting, and flow cytometry assays. The role of GATA3 in modulating doxorubicin-induced cell death was investigated both in vitro and in vivo. RNA-seq, qPCR, ChIP, and luciferase assay and association analyses were performed to investigate the regulation of CYB5R2 by GATA3. The function of GATA3 and CYB5R2 in regulating doxorubicin-induced ferroptosis was evaluated with iron, ROS, and lipid peroxidation detection assays. Immunohistochemistry was performed for results validation.

RESULTS

Doxorubicin-induced basal breast cancer cell death is dependent on iron-mediated ferroptosis. Overexpression of the luminal signature transcriptional factor GATA3 mediates doxorubicin resistance. GATA3 promotes cell viability by decreasing ferroptosis-related gene CYB5R2 expression and by maintaining iron homeostasis. Analyzing data from the public and our cohorts demonstrates that GATA3 and CYB5R2 are associated with NAC response.

CONCLUSIONS

GATA3 promotes doxorubicin resistance by inhibiting CYB5R2-mediated iron metabolism and ferroptosis. Therefore, patients with breast cancer who display high GATA3 expression do not benefit from doxorubicin-based NAC regimens.

P4HTM: A Novel Downstream Target of GATA3 in Breast Cancer

Breast cancer continues to be a major cause of death among women. The GATA3 gene is often overexpressed in breast cancer and is widely used to support a diagnosis. However, lower expression of GATA3 has been linked to poorer prognosis along with frequent gene mutations. Therefore, the role of GATA3 in breast cancer appears to be context specific. This study aims to identify a new downstream target of GATA3 to better understand its regulatory network. Clinical data analysis identified the prolyl 4-hydroxylase transmembrane protein (P4HTM) as one of the most highly co-expressed genes with GATA3. Immunohistochemical staining of breast tumors confirms co-expression between GATA3 and P4HTM at the protein level. Similar to GATA3, P4HTM expression levels are linked to patient prognosis, with lower levels indicating poorer survival. Genomics data found that GATA3 binds to the P4HTM locus, and that ectopic expression of GATA3 in basal breast cancer cells increases the P4HTM transcript level. These results collectively suggest that P4HTM is a novel downstream target of GATA3 in breast cancer and is involved in tumor progression.

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.

Chromatin Variants Reveal the Genetic Determinants of Oncogenesis in Breast Cancer

Breast cancer presents as multiple distinct disease entities. Each tumor harbors diverse cell populations defining a phenotypic heterogeneity that impinges on our ability to treat patients. To date, efforts mainly focused on genetic variants to find drivers of inter- and intratumor phenotypic heterogeneity. However, these efforts have failed to fully capture the genetic basis of breast cancer. Through recent technological and analytical approaches, the genetic basis of phenotypes can now be decoded by characterizing chromatin variants. These variants correspond to polymorphisms in chromatin states at DNA sequences that serve a distinct role across cell populations. Here, we review the function and causes of chromatin variants as they relate to breast cancer inter- and intratumor heterogeneity and how they can guide the development of treatment alternatives to fulfill the goal of precision cancer medicine.

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.

Comparison of the mutational profiles of neuroendocrine breast tumours, invasive ductal carcinomas and pancreatic neuroendocrine carcinomas

The pathophysiology and the optimal treatment of breast neuroendocrine tumours (NETs) are unknown. We compared the mutational profiles of breast NETs (n = 53) with those of 724 publicly available invasive ductal carcinoma (IDC) and 98 pancreatic NET (PNET) cases. The only significantly different pathogenetic or unknown variant rate between breast NETs and IDCs was detected in the TP53 (11.3% in breast NETs and 41% in IDCs, adjusted p value 0.027) and ADCK2 (9.4% in breast NETs vs. 0.28% in IDCs, adjusted p value 0.045) genes. Between breast NETs and PNETs, different pathogenetic or unknown variant frequencies were detected in 30 genes. For example, MEN1 was mutated in only 6% of breast NETs and 37% in PNETs (adjusted p value 0.00050), and GATA3 pathogenetic or unknown variants were only found in 17.0% of breast NETs and 0% in PNETs (adjusted p value 0.0010). The most commonly affected oncogenic pathways in the breast NET cases were PI3K/Akt/mTOR, NOTCH and RTK-RAS pathways. Breast NETs had typically clock-like mutational signatures and signatures associated with defective DNA mismatch repair in their mutational landscape. Our results suggest that the breast NET mutational profile more closely resembles that of IDCs than that of PNETs. These results also revealed several potentially druggable targets, such as MMRd, in breast NETs. In conclusion, breast NETs are indeed a separate breast cancer entity, but their optimal treatment remains to be elucidated.

DNA- and RNA-Binding Proteins Linked Transcriptional Control and Alternative Splicing Together in a Two-Layer Regulatory Network System of Chronic Myeloid Leukemia

DNA- and RNA-binding proteins (DRBPs) typically possess multiple functions to bind both DNA and RNA and regulate gene expression from more than one level. They are controllers for post-transcriptional processes, such as splicing, polyadenylation, transportation, translation, and degradation of RNA transcripts in eukaryotic organisms, as well as regulators on the transcriptional level. Although DRBPs are reported to play critical roles in various developmental processes and diseases, it is still unclear how they work with DNAs and RNAs simultaneously and regulate genes at the transcriptional and post-transcriptional levels. To investigate the functional mechanism of DRBPs, we collected data from a variety of databases and literature and identified 118 DRBPs, which function as both transcription factors (TFs) and splicing factors (SFs), thus called DRBP-SF. Extensive investigations were conducted on four DRBP-SFs that were highly expressed in chronic myeloid leukemia (CML), heterogeneous nuclear ribonucleoprotein K (HNRNPK), heterogeneous nuclear ribonucleoprotein L (HNRNPL), non-POU domain–containing octamer–binding protein (NONO), and TAR DNA-binding protein 43 (TARDBP). By integrating and analyzing ChIP-seq, CLIP-seq, RNA-seq, and shRNA-seq data in K562 using binding and expression target analysis and Statistical Utility for RBP Functions, we discovered a two-layer regulatory network system centered on these four DRBP-SFs and proposed three possible regulatory models where DRBP-SFs can connect transcriptional and alternative splicing regulatory networks cooperatively in CML. The exploration of the identified DRBP-SFs provides new ideas for studying DRBP and regulatory networks, holding promise for further mechanistic discoveries of the two-layer gene regulatory system that may play critical roles in the occurrence and development of CML.

Nuclear-localized, iron-bound superoxide dismutase-2 antagonizes epithelial lineage programs to promote stemness of breast cancer cells via a histone demethylase activity

The dichotomous behavior of superoxide dismutase-2 (SOD2) in cancer biology has long been acknowledged and more recently linked to different posttranslational forms of the enzyme. However, a distinctive activity underlying its tumor-promoting function is yet to be described. Here, we report that acetylation, one of such posttranslational modifications (PTMs), increases SOD2 affinity for iron, effectively changing the biochemical function of this enzyme from that of an antioxidant to a demethylase. Acetylated, iron-bound SOD2 localizes to the nucleus, promoting stem cell gene expression via removal of suppressive epigenetic marks such as H3K9me3 and H3K927me3. Particularly, H3K9me3 was specifically removed from regulatory regions upstream of Nanog and Oct-4, two pluripotency factors involved in cancer stem cell reprogramming. Phenotypically, cells expressing nucleus-targeted SOD2 (NLS-SOD2) have increased clonogenicity and metastatic potential. FeSOD2 operating as H3 demethylase requires H2O2 as substrate, which unlike cofactors of canonical demethylases (i.e., oxygen and 2-oxoglutarate), is more abundant in tumor cells than in normal tissue. Therefore, our results indicate that FeSOD2 is a demethylase with unique activities and functions in the promotion of cancer evolution toward metastatic phenotypes.

Upregulated GATA3/miR205-5p Axis Inhibits MFNG Transcription and Reduces the Malignancy of Triple-Negative Breast Cancer

Simple Summary

Triple-negative cancer (TNBC) is a deadly disease that presents a potential health threat to women worldwide. It is the most aggressive and presents a poor prognosis among all breast cancer subgroups. We previously demonstrated that the elevated expression of manic fringe (MFNG) plays a pivotal role in breast cancer. However, the mechanism through which MFNG is regulated remains obscure. The study presented here set out to determine the mechanism by which MFNG expression is regulated in TNBC. Our findings revealed that GATA3 and miR-205-p cooperatively block the transcription of MFNG leading to the inhibition of cell migration and tumor growth in vitro and in vivo. Our study uncovers a novel GATA3/miR-205-p/MFNG feed-forward loop and miR205-5p could be adopted as a potential therapeutic strategy of TNBC.

Abstract

Triple-negative breast cancer (TNBC) accounts for approximately 20% of all breast carcinomas and has the worst prognosis of all breast cancer subtypes due to the lack of an effective target. Therefore, understanding the molecular mechanism underpinning TNBC progression could explore a new target for therapy. While the Notch pathway is critical in the development process, its dysregulation leads to TNBC initiation. Previously, we found that manic fringe (MFNG) activates the Notch signaling and induces breast cancer progression. However, the underlying molecular mechanism of MFNG upstream remains unknown. In this study, we explore the regulatory mechanisms of MFNG in TNBC. We show that the increased expression of MFNG in TNBC is associated with poor clinical prognosis and significantly promotes cell growth and migration, as well as Notch signaling activation. The mechanistic studies reveal that MFNG is a direct target of GATA3 and miR205-5p and demonstrate that GATA3 and miR205-5p overexpression attenuate MFNG oncogenic effects, while GATA3 knockdown mimics MFNG phenotype to promote TNBC progression. Moreover, we illustrate that GATA3 is required for miR205-5p activation to inhibit MFNG transcription by binding to the 3′ UTR region of its mRNA, which forms the GATA3/miR205-5p/MFNG feed-forward loop. Additionally, our in vivo data show that the miR205-5p mimic combined with polyetherimide-black phosphorus (PEI-BP) nanoparticle remarkably inhibits the growth of TNBC-derived tumors which lack GATA3 expression. Collectively, our study uncovers a novel GATA3/miR205-5p/MFNG feed-forward loop as a pathway that could be a potential therapeutic target for TNBC.

Breast Cancer Prognosis Prediction and Immune Pathway Molecular Analysis Based on Mitochondria-Related Genes

Background

Mitochondria play an important role in breast cancer (BRCA). We aimed to build a prognostic model based on mitochondria-related genes.

Method

Univariate Cox regression analysis, random forest, and the LASSO method were performed in sequence on pretreated TCGA BRCA datasets to screen out genes from a Gene Set Enrichment Analysis, Gene Ontology: biological process gene set to build a prognosis risk score model. Survival analyses and ROC curves were performed to verify the model by using the GSE103091 dataset. The BRCA datasets were equally divided into high- and low-risk score groups. Comparisons between clinical features and immune infiltration related to different risk scores and gene mutation analysis and drug sensitivity prediction were performed for different groups.

Result

Four genes, MRPL36, FEZ1, BMF, and AFG1L, were screened to construct our risk score model in which the higher the risk score, the poorer the prognosis. Univariate and multivariate analyses showed that the risk score was significantly associated with age, M stage, and N stage. The gene mutation probability in the high-risk score group was significantly higher than that in the low-risk score group. Patients with higher risk scores were more likely to die. Drug sensitivity prediction in different groups indicated that PF-562271 and AS601245 might be new inhibitors of BRCA.

Conclusion

We developed a new workable risk score model based on mitochondria-related genes for BRCA prognosis and identified new targets and drugs for BRCA research.

Mammaglobin, GATA-binding protein 3 (GATA3), and epithelial growth factor receptor (EGFR) expression in different breast cancer subtypes and their clinical significance

Increasing evidence has shown that mammaglobin, GATA-binding protein 3 (GATA3), and epithelial growth factor receptor (EGFR) have unique clinical implications for breast cancer subtyping and classification, as well as for breast cancer targeted therapy. It is particularly important to clarify the correlation between their expression and different molecular breast carcinoma subtypes to better understand the molecular basis of the subtypes and to identify effective therapeutic targets for the disease. This study aimed to evaluate mammaglobin, GATA3, and EGFR expression in different breast cancer subtypes, as well as their clinical significance. Subjects of the study included 228 patients with breast cancer at The First Affiliated Hospital of University of Science and Technology of China. They were divided into triple negative (TN), Luminal A, Luminal B, and HER-2 positive (HER-2.P) breast cancer groups based on molecular classification. Immunohistochemical methods were used to detect mammaglobin, GATA3, and EGFR expression in cases of different molecular subtypes before determining the correlation between protein expression and subtype. Mammaglobin and GATA3 expression levels were found to significantly vary with respect to histopathological grade, lymph node status, and molecular subtype; EGFR expression was significantly correlated with breast cancer histopathological grade and molecular subtype. For breast cancer, the expression levels of mammaglobin and GATA3, as well as mammaglobin and EGFR, were significantly correlated. In addition, there was a significantly negative correlation between the expression levels of GATA3 and EGFR in breast cancer tissue samples, especially in HER-2.P samples. These findings provide a theoretical basis for assessing breast cancer clinical prognosis based on the cancer subtype, and hence, have significant practical value.

The hub ten gene-based risk score system using RNA m6A methylation regulator features and tumor immune microenvironment in breast cancer

BACKGROUND

RNA N⁶-methyladenosine (m⁶A) modification is primarily regulated by m⁶A regulators, which play significant epigenetic regulatory roles in tumorigenesis, tumor development, and tumor immune microenvironment. However, the correlation between m⁶A regulators and immune cell infiltration in breast cancer remains unclear.

METHODS

In this study, m⁶A modification patterns were evaluated based on 31 m⁶A modification regulators. m⁶A clusters were determined by consensus clustering. Immune landscape and immune cell infiltration subgroups were characterized by m⁶A clusters. Key module and hub genes related to m⁶A regulators and immune infiltration cells were identified by WGCNA. LASSO algorithm was applied to select prognostic signatures. Multivariate Cox regression analysis was applied to assess the prognostic value of gene signatures.

RESULTS

Two distinct m⁶A clusters were determined based on the expression of 31 m⁶A modification regulators and characterized by two tumor immune microenvironment (TIME) immune cell infiltration subgroups. Further, a total of 1971 differentially expressed genes between breast cancer patients and healthy controls were screened, nine modules associated with clinical characteristics of breast cancer patients were identified. Later, one key module and 13 hub genes correlated with m⁶A regulators and immune infiltration cells were identified. LASSO Cox regression analysis selected and constructed a ten-gene prognostic model to build a risk score system for individual breast cancer patient prognosis. The performance of the ten-gene-based risk score system was further validated in an independent dataset with an AUC of 0.659.

CONCLUSIONS

This study revealed that m⁶A modification regulators played a significant role in the TIME regulation of breast cancer. The hub ten gene-based risk score system is valuable in predicting the prognosis of breast cancer patients, which may provide potential significance for breast cancer diagnosis, prognosis, and immunotherapy in the future.

GATA3 Truncation Mutants Alter EMT Related Gene Expression via Partial Motif Recognition in Luminal Breast Cancer Cells

GATA3 is known to be one of the most frequently mutated genes in breast cancer. More than 10% of breast tumors carry mutations in this gene. However, the functional consequence of GATA3 mutations is still largely unknown. Clinical data suggest that different types of GATA3 mutations may have distinct roles in breast cancer characterization. In this study, we have established three luminal breast cancer cell lines that stably express different truncation mutants (X308 splice site deletion, C321 frameshift, and A333 frameshift mutants) found in breast cancer patients. Transcriptome analysis identified common and distinct gene expression patterns in these GATA3 mutant cell lines. In particular, the impacts on epithelial-to-mesenchymal transition (EMT) related genes are similar across these mutant cell lines. Chromatin localization of the mutants is highly overlapped and exhibits non-canonical motif enrichment. Interestingly, the A333 frameshift mutant expressed cells displayed the most significant impact on the GATA3 binding compared to X308 splice site deletion and C321fs mutants expressed cells. Our results suggest the common and different roles of GATA3 truncation mutations during luminal breast cancer development.

Mutations in Noncoding Cis-Regulatory Elements Reveal Cancer Driver Cistromes in Luminal Breast Cancer

Whole-genome sequencing of primary breast tumors enabled the identification of cancer driver genes and noncoding cancer driver plexuses from somatic mutations. However, differentiating driver from passenger events among noncoding genetic variants remains a challenge. Herein, we reveal cancer-driver cis-regulatory elements linked to transcription factors previously shown to be involved in development of luminal breast cancers by defining a tumor-enriched catalogue of approximately 100,000 unique cis-regulatory elements from 26 primary luminal estrogen receptor (ER)⁺ progesterone receptor (PR)⁺ breast tumors. Integrating this catalog with somatic mutations from 350 publicly available breast tumor whole genomes, we uncovered cancer driver cistromes, defined as the sum of binding sites for a transcription factor, for ten transcription factors in luminal breast cancer such as FOXA1 and ER, nine of which are essential for growth in breast cancer with four exclusive to the luminal subtype. Collectively, we present a strategy to find cancer driver cistromes relying on quantifying the enrichment of noncoding mutations over cis-regulatory elements concatenated into a functional unit. IMPLICATIONS: Mapping the accessible chromatin of luminal breast cancer led to discovery of an accumulation of mutations within cistromes of transcription factors essential to luminal breast cancer. This demonstrates coopting of regulatory networks to drive cancer and provides a framework to derive insight into the noncoding space of cancer.

Mutational signatures in GATA3 transcription factor and its DNA binding domain that stimulate breast cancer and HDR syndrome

Transcription factors (TFs) play important roles in many biochemical processes. Many human genetic disorders have been associated with mutations in the genes encoding these transcription factors, and so those mutations became targets for medications and drug design. In parallel, since many transcription factors act either as tumor suppressors or oncogenes, their mutations are mostly associated with cancer. In this perspective, we studied the GATA3 transcription factor when bound to DNA in a crystal structure and assessed the effect of different mutations encountered in patients with different diseases and phenotypes. We generated all missense mutants of GATA3 protein and DNA within the adjacent and the opposite GATA3:DNA complex models. We mutated every amino acid and studied the new binding of the complex after each mutation. Similarly, we did for every DNA base. We applied Poisson-Boltzmann electrostatic calculations feeding into free energy calculations. After analyzing our data, we identified amino acids and DNA bases keys for binding. Furthermore, we validated those findings against experimental genetic data. Our results are the first to propose in silico modeling for GATA:DNA bound complexes that could be used to score effects of missense mutations in other classes of transcription factors involved in common and genetic diseases.

SARS-CoV-2 infection mediates differential expression of human endogenous retroviruses and long interspersed nuclear elements

SARS-CoV-2 promotes an imbalanced host response that underlies the development and severity of COVID-19. Infections with viruses are known to modulate transposable elements (TEs), which can exert downstream effects by modulating host gene expression, innate immune sensing, or activities encoded by their protein products. We investigated the impact of SARS-CoV-2 infection on TE expression using RNA-Seq data from cell lines and from primary patient samples. Using a bioinformatics tool, Telescope, we showed that SARS-CoV-2 infection led to upregulation or downregulation of TE transcripts, a subset of which differed from cells infected with SARS, Middle East respiratory syndrome coronavirus (MERS-CoV or MERS), influenza A virus (IAV), respiratory syncytial virus (RSV), and human parainfluenza virus type 3 (HPIV3). Differential expression of key retroelements specifically identified distinct virus families, such as Coronaviridae, with unique retroelement expression subdividing viral species. Analysis of ChIP-Seq data showed that TEs differentially expressed in SARS-CoV-2 infection were enriched for binding sites for transcription factors involved in immune responses and for pioneer transcription factors. In samples from patients with COVID-19, there was significant TE overexpression in bronchoalveolar lavage fluid and downregulation in PBMCs. Thus, although the host gene transcriptome is altered by infection with SARS-CoV-2, the retrotranscriptome may contain the most distinctive features of the cellular response to SARS-CoV-2 infection.

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.

The dynamic broad epigenetic (H3K4me3, H3K27ac) domain as a mark of essential genes

Transcriptionally active chromatin is marked by tri-methylation of histone H3 at lysine 4 (H3K4me3) located after first exons and around transcription start sites. This epigenetic mark is typically restricted to narrow regions at the 5`end of the gene body, though a small subset of genes have a broad H3K4me3 domain which extensively covers the coding region. Although most studies focus on the H3K4me3 mark, the broad H3K4me3 domain is associated with a plethora of histone modifications (e.g., H3 acetylated at K27) and is therein termed broad epigenetic domain. Genes marked with the broad epigenetic domain are involved in cell identity and essential cell functions and have clinical potential as biomarkers for patient stratification. Reducing expression of genes with the broad epigenetic domain may increase the metastatic potential of cancer cells. Enhancers and super-enhancers interact with the broad epigenetic domain marked genes forming a hub of interactions involving nucleosome-depleted regions. Together, the regulatory elements coalesce with transcription factors, chromatin modifying/remodeling enzymes, coactivators, and the Mediator and/or Integrator complex into a transcription factory which may be analogous to a liquid–liquid phase-separated condensate. The broad epigenetic domain has a dynamic chromatin structure which supports frequent transcription bursts. In this review, we present the current knowledge of broad epigenetic domains.

Cell-specific expression of ENACα gene by FOXA1 in the glucocorticoid receptor pathway

Introduction:

The glucocorticoid receptor (GR) is one of the most widely studied ligand-dependent nuclear receptors. The combination of transcriptional regulatory factors required for the expression of individual genes targeted by GR varies across cell types; however, the mechanisms underlying this cell type–specific regulation of gene expression are not yet clear.

Methods:

Here, we investigated genes regulated by GR in two different cell lines, A549 and ARPE-19, and examined how gene expression varied according to the effect of pioneer factors using RNA-seq and RT-qPCR.

Results:

Our RNA-seq results identified 19 and 63 genes regulated by GR that are ARPE-19-specific and A549-specific, respectively, suggesting that GR induces the expression of different sets of genes in a cell type–specific manner. RT-qPCR confirmed that the epithelial sodium channel (ENACα) gene is an ARPE-19 cell-specific GR target gene, whereas the FK506 binding protein 5 (FKBP5) gene was A549 cell-specific. There was a significant decrease in ENACα expression in FOXA1-deficient ARPE-19 cells, suggesting that FOXA1 might function as a pioneer factor enabling the selective expression of ENACα in ARPE-19 cells but not in A549 cells.

Conclusion:

These findings indicate that ENACα expression in ARPE-19 cells is regulated by FOXA1 and provide insights into the molecular mechanisms of cell type–specific expression of GR-regulated genes.

The estrogen receptor/GATA3/FOXA1 transcriptional network: lessons learned from breast cancer

Cellular identity and physiologic function in mammary epithelial cells and in many breast cancers flow from the action of a network of master transcriptional regulators including estrogen receptor alpha, GATA3, and FOXA1. The last decade has seen the completion of multiple large sequencing projects focusing on breast cancer. These massive compendia of sequence data have provided a wealth of new information linking mutation in these transcription factors to alterations in tumor biology and transcriptional program. The emerging details on mutation in cancer, and direct experimental exploration of hypotheses based on it, are now providing a wealth of new information on the roles played by estrogen receptor alpha, GATA3, and FOXA1 in regulating gene transcription and how their combined action contributes to a network shaping cell function in both physiologic and disease states.

The molecular basis of mammary gland development and epithelial differentiation

Our understanding of the molecular events underpinning the development of mammalian organ systems has been increasing rapidly in recent years. With the advent of new and improved next-generation sequencing methods, we are now able to dig deeper than ever before into the genomic and epigenomic events that play critical roles in determining the fates of stem and progenitor cells during the development of an embryo into an adult. In this review, we detail and discuss the genes and pathways that are involved in mammary gland development, from embryogenesis, through maturation into an adult gland, to the role of pregnancy signals in directing the terminal maturation of the mammary gland into a milk producing organ that can nurture the offspring. We also provide an overview of the latest research in the single-cell genomics of mammary gland development, which may help us to understand the lineage commitment of mammary stem cells (MaSCs) into luminal or basal epithelial cells that constitute the mammary gland. Finally, we summarize the use of 3D organoid cultures as a model system to study the molecular events during mammary gland development. Our increased investigation of the molecular requirements for normal mammary gland development will advance the discovery of targets to predict breast cancer risk and the development of new breast cancer therapies.

LncSEA: a platform for long non-coding RNA related sets and enrichment analysis

Long non-coding RNAs (lncRNAs) have been proven to play important roles in transcriptional processes and various biological functions. Establishing a comprehensive collection of human lncRNA sets is urgent work at present. Using reference lncRNA sets, enrichment analyses will be useful for analyzing lncRNA lists of interest submitted by users. Therefore, we developed a human lncRNA sets database, called LncSEA, which aimed to document a large number of available resources for human lncRNA sets and provide annotation and enrichment analyses for lncRNAs. LncSEA supports >40 000 lncRNA reference sets across 18 categories and 66 sub-categories, and covers over 50 000 lncRNAs. We not only collected lncRNA sets based on downstream regulatory data sources, but also identified a large number of lncRNA sets regulated by upstream transcription factors (TFs) and DNA regulatory elements by integrating TF ChIP-seq, DNase-seq, ATAC-seq and H3K27ac ChIP-seq data. Importantly, LncSEA provides annotation and enrichment analyses of lncRNA sets associated with upstream regulators and downstream targets. In summary, LncSEA is a powerful platform that provides a variety of types of lncRNA sets for users, and supports lncRNA annotations and enrichment analyses. The LncSEA database is freely accessible at http://bio.liclab.net/LncSEA/index.php.

Estrogen Receptor on the move: Cistromic plasticity and its implications in breast cancer

Estrogen Receptor (ERα) is a hormone-driven transcription factor, critically involved in driving tumor cell proliferation in the vast majority of breast cancers (BCas). ERα binds the genome at cis-regulatory elements, dictating the expression of a large spectrum of responsive genes in 3D genomic space. While initial reports described a rather static ERα chromatin binding repertoire, we now know that ERα DNA interactions are highly versatile, altered in breast tumor development and progression, and deviate between tumors from patients with differential outcome. Multiple cellular signaling cascades are known to impinge on ERα genomic function, changing its cistrome to retarget the receptor to other regions of the genome and reprogram its impact on breast cell biology. This review describes the current state-of-the-art on which factors manipulate the ERα cistrome and how this alters the response to both endogenous and exogenous hormonal stimuli, ultimately impacting BCa cell progression and response to commonly used therapeutic interventions. Novel insights in ERα cistrome dynamics may pave the way for better patient diagnostics and the development of novel therapeutic interventions, ultimately improving cancer care and patient outcome.