Whole-genome cartography of estrogen receptor alpha binding sites

Relaxin Modulates the Genomic Actions and Biological Effects of Estrogen in the Myometrium

Estradiol (E2) and relaxin (Rln) are steroid and polypeptide hormones, respectively, with important roles in the female reproductive tract, including myometrium. Some actions of Rln, which are mediated by its membrane receptor RXFP1, require or are augmented by E2 signaling through its cognate nuclear steroid receptor, estrogen receptor alpha (ERα). In contrast, other actions of Rln act in opposition to the effects of E2. Here we explored the molecular and genomic mechanisms that underlie the functional interplay between E2 and Rln in the myometrium. We used both ovariectomized female mice and immortalized human myometrial cells expressing wild-type or mutant ERα (hTERT-HM-ERα cells). Our results indicate that Rln modulates the genomic actions and biological effects of estrogen in the myometrium and myometrial cells by reducing phosphorylation of ERα on serine 118 (S118), as well as by reducing the E2-dependent binding of ERα across the genome. These effects were associated with changes in the hormone-regulated transcriptome, including a decrease in the E2-dependent expression of some genes and enhanced expression of others. The inhibitory effects of Rln cotreatment on the E2-dependent phosphorylation of ERα required the nuclear dual-specificity phosphatases DUSP1 and DUSP5. Moreover, the inhibitory effects of Rln were reflected in a concomitant inhibition of the E2-dependent contraction of myometrial cells. Collectively, our results identify a pathway that integrates Rln/RXFP1 and E2/ERα signaling, resulting in a convergence of membrane and nuclear signaling pathways to control genomic and biological outcomes.

Genes Co-Expressed with ESR2 Influence Clinical Outcomes in Cancer Patients: TCGA Data Analysis

ERβ has been assigned a tumor suppressor role in many cancer types. However, as conflicting findings emerge, ERβ's tissue-specific expression and functional role have remained elusive. There remains a notable gap in compact and comprehensive analyses of ESR2 mRNA expression levels across diverse tumor types coupled with an exploration of its potential gene network. In this study, we aim to address these gaps by presenting a comprehensive analysis of ESR2 transcriptomic data. We distinguished cancer types with significant changes in ESR2 expression levels compared to corresponding healthy tissue and concluded that ESR2 influences patient survival. Gene Set Enrichment Analysis (GSEA) distinguished molecular pathways affected by ESR2, including oxidative phosphorylation and epithelial-mesenchymal transition. Finally, we investigated genes displaying similar expression patterns as ESR2 in tumor tissues, identifying potential co-expressed genes that may exert a synergistic effect on clinical outcomes, with significant results, including the expression of ACIN1, SYNE2, TNFRSF13C, and MDM4. Collectively, our results highlight the significant influence of ESR2 mRNA expression on the transcriptomic landscape and the overall metabolism of cancerous cells across various tumor types.

REThinking the role of the RET oncogene in breast cancer

The REarranged during Transfection (RET) receptor tyrosine kinase plays a crucial role in the development of various anatomical structures during embryogenesis and it is involved in many physiological cellular processes. This protein is also associated with the initiation of various cancer types, such as thyroid cancer, non-small cell lung cancer, and multiple endocrine neoplasms. In breast cancer, and especially in the estrogen receptor-positive (ER+) subtype, the activity of RET is of notable importance. Indeed, RET seems to be involved in tumor progression, resistance to therapies, and cellular proliferation. Nevertheless, the ways RET alterations could impact the prognosis of breast cancer and its response to treatment remain only partially elucidated. Several inhibitors of RET kinase have been developed thus far, with various degrees of selectivity toward RET inhibition. These molecules showed notable efficacy in the treatment of RET-driven tumors, including some breast cancer cases. Despite these encouraging results, further investigation is needed to fully understand the potential role RET inhibition in breast cancer. This review aims to recapitulate the existing evidence about the role of RET oncogene in breast cancer, from its pathogenic and potentially prognostic role, to the clinical applications of RET inhibitors.

PSMD14 stabilizes estrogen signaling and facilitates breast cancer progression via deubiquitinating ERα

The over-activation of ERα signaling is regarded as the major driver for luminal breast cancers, which could be effective controlled via selective estrogen receptor modulators (SERM), such as tamoxifen. The endocrine resistance is still a challenge for breast cancer treatment, while recently studies implicate the post-translational modification on ERα play important roles in endocrine resistance. The stability of ERα protein and ERα transcriptome are subject to a balance between E3 ubiquitin ligases and deubiquitinases. Through deubiquitinases siRNA library screening, we discover PSMD14 as a critical deubiquitinase for ERα signaling and breast cancer progression. PSMD14 could facilitate breast cancer progression through ERα signaling in vitro and in vivo, while pharmaceutical inhibition of PSMD14 via Thiolutin could block the tumorigenesis in breast cancer. In endocrine resistant models, PSMD14 inhibition could de-stabilize the resistant form of ERα (Y537S) and restore tamoxifen sensitivity. Molecular studies reveal that PSMD14 could inhibition K48-linked poly-ubiquitination on ERα, facilitate ERα transcriptome. Interestingly, ChIP assay shows that ERα could bind to the promoter region of PSMD14 and facilitate its gene transcription, which indicates PSMD14 is both the upstream modulator and downstream target for ERα signaling in breast cancer. In general, we identified a novel positive feedback loop between PSMD14 and ERα signaling in breast cancer progression, while blockade of PSMD14 could be a plausible strategy for luminal breast cancer.

Gene transcription regulation by ER at the single cell and allele level

In this short review we discuss the current view of how the estrogen receptor (ER), a pivotal member of the nuclear receptor superfamily of transcription factors, regulates gene transcription at the single cell and allele level, focusing on in vitro cell line models. We discuss central topics and new trends in molecular biology including phenotypic heterogeneity, single cell sequencing, nuclear phase separated condensates, single cell imaging, and image analysis methods, with particular focus on the methodologies and results that have been reported in the last few years using microscopy-based techniques. These observations augment the results from biochemical assays that lead to a much more complex and dynamic view of how ER, and arguably most transcription factors, act to regulate gene transcription.

Regulation of mRNA translation by estrogen receptor in breast cancer

Breast cancer is one of the leading causes of cancer-related fatalities and the most often diagnosed malignancy in women globally. Dysregulation of sex hormone signaling pathways mediated by the estrogen receptor (ER) in breast cancer is well characterized. Although ER is known to promote cell growth and survival by altering gene transcription, recent research suggests that its effects in cancers are also mediated through dysregulation of protein synthesis. This implies that ER can coordinately affect gene expression through both translational and transcriptional pathways, leading to the development of malignancy. In this review, we will cover the current understanding of how the ER controls mRNA translation in breast cancer and discuss any potential clinical implications of this phenomenon.

Ovarian ERβ cistrome and transcriptome reveal chromatin interaction with LRH-1

BACKGROUND

Estrogen receptor beta (ERβ, Esr2) plays a pivotal role in folliculogenesis and ovulation, yet its exact mechanism of action is mainly uncharacterized.

RESULTS

We here performed ERβ ChIP-sequencing of mouse ovaries followed by complementary RNA-sequencing of wild-type and ERβ knockout ovaries. By integrating the ERβ cistrome and transcriptome, we identified its direct target genes and enriched biological functions in the ovary. This demonstrated its strong impact on genes regulating organism development, cell migration, lipid metabolism, response to hypoxia, and response to estrogen. Cell-type deconvolution analysis of the bulk RNA-seq data revealed a decrease in luteal cells and an increased proportion of theca cells and a specific type of cumulus cells upon ERβ loss. Moreover, we identified a significant overlap with the gene regulatory network of liver receptor homolog 1 (LRH-1, Nr5a2) and showed that ERβ and LRH-1 extensively bound to the same chromatin locations in granulosa cells. Using ChIP-reChIP, we corroborated simultaneous ERβ and LRH-1 co-binding at the ERβ-repressed gene Greb1 but not at the ERβ-upregulated genes Cyp11a1 and Fkbp5. Transactivation assay experimentation further showed that ERβ and LRH-1 can inhibit their respective transcriptional activity at classical response elements.

CONCLUSIONS

By characterizing the genome-wide endogenous ERβ chromatin binding, gene regulations, and extensive crosstalk between ERβ and LRH-1, along with experimental corroborations, our data offer genome-wide mechanistic underpinnings of ovarian physiology and fertility.

International Union of Basic and Clinical Pharmacology CXIII: Nuclear Receptor Superfamily-Update 2023

The NR superfamily comprises 48 transcription factors in humans that control a plethora of gene network programs involved in a wide range of physiologic processes. This review will summarize and discuss recent progress in NR biology and drug development derived from integrating various approaches, including biophysical techniques, structural studies, and translational investigation. We also highlight how defective NR signaling results in various diseases and disorders and how NRs can be targeted for therapeutic intervention via modulation via binding to synthetic lipophilic ligands. Furthermore, we also review recent studies that improved our understanding of NR structure and signaling. SIGNIFICANCE STATEMENT: Nuclear receptors (NRs) are ligand-regulated transcription factors that are critical regulators of myriad physiological processes. NRs serve as receptors for an array of drugs, and in this review, we provide an update on recent research into the roles of these drug targets.

A novel approach for predicting upstream regulators (PURE) that affect gene expression

External factors such as exposure to a chemical, drug, or toxicant (CDT), or conversely, the lack of certain chemicals can cause many diseases. The ability to identify such causal CDTs based on changes in the gene expression profile is extremely important in many studies. Furthermore, the ability to correctly infer CDTs that can revert the gene expression changes induced by a given disease phenotype is a crucial step in drug repurposing. We present an approach for Predicting Upstream REgulators (PURE) designed to tackle this challenge. PURE can correctly infer a CDT from the measured expression changes in a given phenotype, as well as correctly identify drugs that could revert disease-induced gene expression changes. We compared the proposed approach with four classical approaches as well as with the causal analysis used in Ingenuity Pathway Analysis (IPA) on 16 data sets (1 rat, 5 mouse, and 10 human data sets), involving 8 chemicals or drugs. We assessed the results based on the ability to correctly identify the CDT as indicated by its rank. We also considered the number of false positives, i.e. CDTs other than the correct CDT that were reported to be significant by each method. The proposed approach performed best in 11 out of the 16 experiments, reporting the correct CDT at the very top 7 times. IPA was the second best, reporting the correct CDT at the top 5 times, but was unable to identify the correct CDT at all in 5 out of the 16 experiments. The validation results showed that our approach, PURE, outperformed some of the most popular methods in the field. PURE could effectively infer the true CDTs responsible for the observed gene expression changes and could also be useful in drug repurposing applications.

Dax1 modulates ERα-dependent hypothalamic estrogen sensing in female mice

Coupling the release of pituitary hormones to the developmental stage of the oocyte is essential for female fertility. It requires estrogen to restrain kisspeptin (KISS1)-neuron pulsatility in the arcuate hypothalamic nucleus, while also exerting a surge-like effect on KISS1-neuron activity in the AVPV hypothalamic nucleus. However, a mechanistic basis for this region-specific effect has remained elusive. Our genomic analysis in female mice demonstrate that some processes, such as restraint of KISS1-neuron activity in the arcuate nucleus, may be explained by region-specific estrogen receptor alpha (ERα) DNA binding at gene regulatory regions. Furthermore, we find that the Kiss1-locus is uniquely regulated in these hypothalamic nuclei, and that the nuclear receptor co-repressor NR0B1 (DAX1) restrains its transcription specifically in the arcuate nucleus. These studies provide mechanistic insight into how ERα may control the KISS1-neuron, and Kiss1 gene expression, to couple gonadotropin release to the developmental stage of the oocyte.

Predicting cross-tissue hormone-gene relations using balanced word embeddings

MOTIVATION

Inter-organ/inter-tissue communication is central to multi-cellular organisms including humans, and mapping inter-tissue interactions can advance system-level whole-body modeling efforts. Large volumes of biomedical literature have fostered studies that map within-tissue or tissue-agnostic interactions, but literature-mining studies that infer inter-tissue relations, such as between hormones and genes are solely missing.

RESULTS

We present a first study to predict from biomedical literature the hormone-gene associations mediating inter-tissue signaling in the human body. Our BioEmbedS* models use neural network-based Biomedical word Embeddings with a Support Vector Machine classifier to predict if a hormone-gene pair is associated or not, and whether an associated gene is involved in the hormone's production or response. Model training relies on our unified dataset Hormone-Gene version 1 of ground-truth associations between genes and endocrine hormones, which we compiled and carefully balanced in the embedded space to handle data disparities, such as between poorly- versus well-studied hormones. Our BioEmbedS model recapitulates known gene mediators of tissue-tissue signaling with 70.4% accuracy; predicts novel inter-tissue communication genes in humans, which are enriched for hormone-related disorders; and generalizes well to mouse, thereby holding promise for its extension to other multi-cellular organisms as well.

AVAILABILITY AND IMPLEMENTATION

Freely available at https://cross-tissue-signaling.herokuapp.com are our model predictions & datasets; https://github.com/BIRDSgroup/BioEmbedS has all relevant code.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

SOPHIE: Generative Neural Networks Separate Common and Specific Transcriptional Responses

Genome-wide transcriptome profiling identifies genes that are prone to differential expression (DE) across contexts, as well as genes with changes specific to the experimental manipulation. Distinguishing genes that are specifically changed in a context of interest from common differentially expressed genes (DEGs) allows more efficient prediction of which genes are specific to a given biological process under scrutiny. Currently, common DEGs or pathways can only be identified through the laborious manual curation of experiments, an inordinately time-consuming endeavor. Here we pioneer an approach, Specific cOntext Pattern Highlighting In Expression data (SOPHIE), for distinguishing between common and specific transcriptional patterns using a generative neural network to create a background set of experiments from which a null distribution of gene and pathway changes can be generated. We apply SOPHIE to diverse datasets including those from human, human cancer, and bacterial pathogen Pseudomonas aeruginosa. SOPHIE identifies common DEGs in concordance with previously described, manually and systematically determined common DEGs. Further molecular validation indicates that SOPHIE detects highly specific but low-magnitude biologically relevant transcriptional changes. SOPHIE's measure of specificity can complement log2 fold change values generated from traditional DE analyses. For example, by filtering the set of DEGs, one can identify genes that are specifically relevant to the experimental condition of interest. Consequently, these results can inform future research directions. All scripts used in these analyses are available at https://github.com/greenelab/generic-expression-patterns. Users can access https://github.com/greenelab/sophie to run SOPHIE on their own data.

An approach to uncover the relationship between 17b-estradiol and ESR1/ESR2 ratio in the regulation of canine corpus luteum

The canine corpus luteum (CL) is able to synthetise, activate and deactivate 17b-estradiol (E2) and also expresses nuclear estrogen receptors in a time-dependent manner during diestrus. Nevertheless, we are still missing a better comprehension of E2 functions in the canine CL, especially regarding the specific roles of estrogen receptor alpha (ERa) and ERb, encoded by ESR1 and 2, respectively. For that purpose, we analyzed transcriptomic data of canine non-pregnant CL collected on days 10, 20, 30, 40, 50 and 60 of diestrus and searched for differentially expressed genes (DEG) containing predicted transcription factor binding sites (TFBS) for ESR1 or ESR2. Based on biological functions of DEG presenting TFBS, expression of select transcripts and corresponding proteins was assessed. Additionally, luteal cells were collected across specific time points during diestrus and specificity of E2 responses was tested using ERa and/or ERb inhibitors. Bioinformatic analyses revealed 517 DEGs containing TFBS, from which 67 for both receptors. In general, abundance of predicted ESR1 targets was greater in the beginning, while abundance of ESR2 targets was greater in the end of diestrus. ESR1/ESR2 ratio shifted from an increasing to a decreasing pattern from day 30 to 40 post ovulation. Specific receptor inhibition suggested an ERa-mediated positive regulation of CL function at the beginning of diestrus and an ERb-mediated effect contributing to luteal regression. In conclusion, our data points toward a broad spectrum of action of E2 and its nuclear receptors, which can also act as transcription factors for other genes regulating canine CL function.

FL-DTD: an integrated pipeline to predict the drug interacting targets by feedback loop-based network analysis

Drug target discovery is an essential step to reveal the mechanism of action (MoA) underlying drug therapeutic effects and/or side effects. Most of the approaches are usually labor-intensive while unable to identify the tissue-specific interacting targets, especially the targets with weaker drug binding affinity. In this work, we proposed an integrated pipeline, FL-DTD, to predict the drug interacting targets of novel compounds in a tissue-specific manner. This method was built based on a hypothesis that cells under a status of homeostasis would take responses to drug perturbation by activating feedback loops. Therefore, the drug interacting targets can be predicted by analyzing the network responses after drug perturbation. We evaluated this method using the expression data of estrogen stimulation, gene manipulation and drug perturbation and validated its good performance to identify the annotated drug targets. Using STAT3 as a target protein, we applied this method to drug perturbation data of 500 natural compounds and predicted five compounds with STAT3 interacting activities. Experimental assay validated the STAT3-interacting activities of four compounds. Overall, our evaluation suggests that FL-DTD predicts the drug interacting targets with good accuracy and can be used for drug target discovery.

Allosteric interactions prime androgen receptor dimerization and activation

The androgen receptor (AR) is a nuclear receptor that governs gene expression programs required for prostate development and male phenotype maintenance. Advanced prostate cancers display AR hyperactivation and transcriptome expansion, in part, through AR amplification and interaction with oncoprotein cofactors. Despite its biological importance, how AR domains and cofactors cooperate to bind DNA has remained elusive. Using single-particle cryo-electron microscopy, we isolated three conformations of AR bound to DNA, showing that AR forms a non-obligate dimer, with the buried dimer interface utilized by ancestral steroid receptors repurposed to facilitate cooperative DNA binding. We identify novel allosteric surfaces which are compromised in androgen insensitivity syndrome and reinforced by AR's oncoprotein cofactor, ERG, and by DNA-binding motifs. Finally, we present evidence that this plastic dimer interface may have been adopted for transactivation at the expense of DNA binding. Our work highlights how fine-tuning AR's cooperative interactions translate to consequences in development and disease.

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

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

Estrogens and phytoestrogens in body functions

Estrogens are the hormones of reproduction in women as well as of many other important functions in the male and female body. They undergo significant changes in the different phases of life, e.g. during puberty, pregnancy or at menopause/andropause. Phytoestrogens are natural non-steroidal phenolic plant compounds that can mimic the activity of estrogens and their beneficial effects in women and in men. This narrative review summarizes the literature on the physiological role of estrogens and the several potential health benefits of phytoestrogens, with particular attention given to the possible role of phytoestrogens in aging.

Oestrogen engages brain MC4R signalling to drive physical activity in female mice

Oestrogen depletion in rodents and humans leads to inactivity, fat accumulation and diabetes^1,2, underscoring the conserved metabolic benefits of oestrogen that inevitably decrease with age. In rodents, the preovulatory surge in 17β-oestradiol (E2) temporarily increases energy expenditure to coordinate increased physical activity with peak sexual receptivity. Here we report that a subset of oestrogen-sensitive neurons in the ventrolateral ventromedial hypothalamic nucleus (VMHvl)^3-7 projects to arousal centres in the hippocampus and hindbrain, and enables oestrogen to rebalance energy allocation in female mice. Surges in E2 increase melanocortin-4 receptor (MC4R) signalling in these VMHvl neurons by directly recruiting oestrogen receptor-α (ERα) to the Mc4r gene. Sedentary behaviour and obesity in oestrogen-depleted female mice were reversed after chemogenetic stimulation of VMHvl neurons expressing both MC4R and ERα. Similarly, a long-term increase in physical activity is observed after CRISPR-mediated activation of this node. These data extend the effect of MC4R signalling - the most common cause of monogenic human obesity⁸ - beyond the regulation of food intake and rationalize reported sex differences in melanocortin signalling, including greater disease severity of MC4R insufficiency in women⁹. This hormone-dependent node illuminates the power of oestrogen during the reproductive cycle in motivating behaviour and maintaining an active lifestyle in women.

FOXA1: A Pioneer of Nuclear Receptor Action in Breast Cancer

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

Sex differences in skeletal muscle revealed through fiber type, capillarity, and transcriptomics profiling in mice

Skeletal muscle anatomy and physiology are sexually dimorphic but molecular underpinnings and muscle-specificity are not well-established. Variances in metabolic health, fitness level, sedentary behavior, genetics, and age make it difficult to discern inherent sex effects in humans. Therefore, mice under well-controlled conditions were used to determine female and male (n = 19/sex) skeletal muscle fiber type/size and capillarity in superficial and deep gastrocnemius (GA-s, GA-d), soleus (SOL), extensor digitorum longus (EDL), and plantaris (PLT), and transcriptome patterns were also determined (GA, SOL). Summed muscle weight strongly correlated with lean body mass (r2  = 0.67, p < 0.0001, both sexes). Other phenotypes were muscle-specific: e.g., capillarity (higher density, male GA-s), myofiber size (higher, male EDL), and fiber type (higher, lower type I and type II prevalences, respectively, in female SOL). There were broad differences in transcriptomics, with >6000 (GA) and >4000 (SOL) mRNAs differentially-expressed by sex; only a minority of these were shared across GA and SOL. Pathway analyses revealed differences in ribosome biology, transcription, and RNA processing. Curation of sexually dimorphic muscle transcripts shared in GA and SOL, and literature datasets from mice and humans, identified 11 genes that we propose are canonical to innate sex differences in muscle: Xist, Kdm6a, Grb10, Oas2, Rps4x (higher, females) and Ddx3y, Kdm5d, Irx3, Wwp1, Aldh1a1, Cd24a (higher, males). These genes and those with the highest "sex-biased" expression in our study do not contain estrogen-response elements (exception, Greb1), but a subset are proposed to be regulated through androgen response elements. We hypothesize that innate muscle sexual dimorphism in mice and humans is triggered and then maintained by classic X inactivation (Xist, females) and Y activation (Ddx3y, males), with coincident engagement of X encoded (Kdm6a) and Y encoded (Kdm5d) demethylase epigenetic regulators that are complemented by modulation at some regions of the genome that respond to androgen.

Estrogen receptor-α signaling in post-natal mammary development and breast cancers

17β-estradiol controls post-natal mammary gland development and exerts its effects through Estrogen Receptor ERα, a member of the nuclear receptor family. ERα is also critical for breast cancer progression and remains a central therapeutic target for hormone-dependent breast cancers. In this review, we summarize the current understanding of the complex ERα signaling pathways that involve either classical nuclear “genomic” or membrane “non-genomic” actions and regulate in concert with other hormones the different stages of mammary development. We describe the cellular and molecular features of the luminal cell lineage expressing ERα and provide an overview of the transgenic mouse models impacting ERα signaling, highlighting the pivotal role of ERα in mammary gland morphogenesis and function and its implication in the tumorigenic processes. Finally, we describe the main features of the ERα-positive luminal breast cancers and their modeling in mice.

Global investigation of estrogen-responsive genes regulating lipid metabolism in the liver of laying hens

BACKGROUND

Estrogen plays an essential role in female development and reproductive function. In chickens, estrogen is critical for lipid metabolism in the liver. The regulatory molecular network of estrogen in chicken liver is poorly understood. To identify estrogen-responsive genes and estrogen functional sites on a genome-wide scale, we determined expression profiles of mRNAs, lncRNAs, and miRNAs in estrogen-treated ((17β-estradiol)) and control chicken livers using RNA-Sequencing (RNA-Seq) and studied the estrogen receptor α binding sites by ChIP-Sequencing (ChIP-Seq).

RESULTS

We identified a total of 990 estrogen-responsive genes, including 962 protein-coding genes, 11 miRNAs, and 17 lncRNAs. Functional enrichment analyses showed that the estrogen-responsive genes were highly enriched in lipid metabolism and biological processes. Integrated analysis of the data of RNA-Seq and ChIP-Seq, identified 191 genes directly targeted by estrogen, including 185 protein-coding genes, 4 miRNAs, and 2 lncRNAs. In vivo and in vitro experiments showed that estrogen decreased the mRNA expression of PPARGC1B, which had been reported to be linked with lipid metabolism, by directly increasing the expression of miR-144-3p.

CONCLUSIONS

These results increase our understanding of the functional network of estrogen in chicken liver and also reveal aspects of the molecular mechanism of estrogen-related lipid metabolism.

Serum C-X-C motif chemokine ligand 14 levels are associated with serum C-peptide and fatty liver index in type 2 diabetes mellitus patients

AIMS/INTRODUCTION

Recent studies have suggested C-X-C motif chemokine ligand 14 (CXCL14), secreted from adipose tissue, to play an important role in the pathogenesis of metabolic syndrome. However, the clinical significance of CXCL14 in humans has not been elucidated. This study aimed to assess correlations between serum CXCL14 levels and clinical parameters in patients with type 2 diabetes mellitus.

MATERIALS AND METHODS

In total, 176 individuals with type 2 diabetes mellitus were recruited. Serum CXCL14 concentrations were determined by enzyme-linked immunosorbent assay. We examined the associations of serum CXCL14 levels with laboratory values, abdominal computed tomography image information, surrogate markers used for evaluating the pathological states of diabetes, obesity and atherosclerosis.

RESULTS

Serum CXCL14 levels correlated positively with body mass index, waist circumference, subcutaneous and visceral fat areas, and serum alanine transaminase, uric acid, total cholesterol, low-density lipoprotein cholesterol, triglycerides and C-peptide (CPR) levels. In contrast, CXCL14 levels correlated inversely with age, pulse wave velocity and serum adiponectin levels. Multiple linear regression analysis showed serum levels of CPR (β = 0.227, P = 0.038) and the fatty liver index (β = 0.205, P = 0.049) to be the only parameters showing independent statistically significant associations with serum CXCL14 levels.

CONCLUSIONS

Serum CXCL14 levels were independently associated with serum CPR and fatty liver index in patients with type 2 diabetes mellitus. In these patients, a high serum CPR concentration might reflect insulin resistance rather than β-cell function, because CXCL14 showed simple correlations with obesity-related parameters. Collectively, these data suggested that serum CXCL14 levels in type 2 diabetes patients might be useful predictors of elevated serum CPR and hepatic steatosis.

PredTAD: A machine learning framework that models 3D chromatin organization alterations leading to oncogene dysregulation in breast cancer cell lines

Topologically associating domains, or TADs, play important roles in genome organization and gene regulation; however, they are often altered in diseases. High-throughput chromatin conformation capturing assays, such as Hi-C, can capture domains of increased interactions, and TADs and boundaries can be identified using well-established analytical tools. However, generating Hi-C data is expensive. In our study, we addressed the relationship between multi-omics data and higher-order chromatin structures using a newly developed machine-learning model called PredTAD. Our tool uses already-available and cost-effective datatypes such as transcription factor and histone modification ChIPseq data. Specifically, PredTAD utilizes both epigenetic and genetic features as well as neighboring information to classify the entire human genome as boundary or non-boundary regions. Our tool can predict boundary changes between normal and breast cancer genomes. Among the most important features for predicting boundary alterations were CTCF, subunits of cohesin (RAD21 and SMC3), and chromosome number, suggesting their roles in conserved and dynamic boundaries formation. Upon further analysis, we observed that genes near altered TAD boundaries were found to be involved in several important breast cancer signaling pathways such as Ras, Jak-STAT, and estrogen signaling pathways. We also discovered a TAD boundary alteration that contributes to RET oncogene overexpression. PredTAD can also successfully predict TAD boundary changes in other conditions and diseases. In conclusion, our newly developed machine learning tool allowed for a more complete understanding of the dynamic 3D chromatin structures involved in signaling pathway activation, altered gene expression, and disease state in breast cancer cells.

Genome-wide estrogen receptor β chromatin binding in human colon cancer cells reveals its tumor suppressor activity

Colorectal cancer (CRC) is the third leading cause of cancer death in the western world. In women, menopausal hormone therapy has been shown to reduce CRC incidence by 20%. Studies demonstrate that estrogen activating estrogen receptor beta (ERβ) protects against CRC. ERβ is a nuclear receptor that regulates gene expression through interactions with the chromatin. This molecular mechanism is, however, not well characterized in colon. Here, we present for the first time, the cistrome of ERβ in different colon cancer cell lines. We use cell lines engineered to express ERβ, optimize and validate an ERβ antibody for chromatin-immunoprecipitation (ChIP), and perform ChIP-Seq. We identify key binding motifs, including ERE, AP-1, and TCF sites, and we determine enrichment of binding to cis-regulatory chromatin sites of genes involved in tumor development, cell migration, cell adhesion, apoptosis, and Wnt signaling pathways. We compare the corresponding cistromes of colon and breast cancer and find that they are conserved for about a third of genes, including GREB1, but that ERβ tethering to TCF and KLF family motifs is characteristic for colon. We exemplify upregulation of putative CRC tumor suppressor gene CST5 where ERβ in colon cells binds to cis-regulatory regions nearby (-351 bp) the transcriptional start site. Our work provides a foundation for understanding the mechanism of action of ERβ in CRC prevention.

TET2 is a component of the estrogen receptor complex and controls 5mC to 5hmC conversion at estrogen receptor cis-regulatory regions

Estrogen receptor-α (ER) drives tumor development in ER-positive (ER+) breast cancer. The transcription factor GATA3 has been closely linked to ER function, but its precise role in this setting remains unclear. Quantitative proteomics was used to assess changes to the ER complex in response to GATA3 depletion. Unexpectedly, few proteins were lost from the ER complex in the absence of GATA3, with the only major change being depletion of the dioxygenase TET2. TET2 binding constituted a near-total subset of ER binding in multiple breast cancer models, with loss of TET2 associated with reduced activation of proliferative pathways. TET2 knockdown did not appear to change global methylated cytosine (5mC) levels; however, oxidation of 5mC to 5-hydroxymethylcytosine (5hmC) was significantly reduced, and these events occurred at ER enhancers. These findings implicate TET2 in the maintenance of 5hmC at ER sites, providing a potential mechanism for TET2-mediated regulation of ER target genes.

ERα in the Control of Mitochondrial Function and Metabolic Health

Decrements in metabolic health elevate disease risk, including type 2 diabetes, heart disease, and certain cancers. Thus, treatment strategies to combat metabolic dysfunction are needed. Reduced ESR1 (estrogen receptor, ERα) expression is observed in muscle from women, men, and animals presenting clinical features of the metabolic syndrome. Human studies of natural expression of ESR1 in metabolic tissues show that muscle expression of ESR1 is positively correlated with markers of metabolic health, including insulin sensitivity. Herein, we highlight the important impact of ERα on mitochondrial form and function and present how these actions of the receptor govern metabolic homeostasis. Studies identifying ERα-regulated pathways for disease prevention will lay the foundation for the design of novel therapeutics to improve the health of women while limiting secondary complications that have plagued traditional hormone replacement interventions.

Estrogen-induced hypomethylation and overexpression of YAP1 facilitate breast cancer cell growth and survival

Increased expression of Yes-associated protein-1 (YAP1) was shown to correlate with reduced survival in breast cancer (BC) patients. However, the exact mechanism of YAP1 regulation in BC cells remains ambiguous. Genomic sequence search showed that the promoter region of the YAP1 gene contains CpG Islands, hence the likelihood of epigenetic regulation by DNA methylation. To address this possibility, the effect of estrogen (17β estradiol; E2) on YAP1 gene expression and YAP1 promoter methylation status was evaluated in BC cells. The functional consequences of E2 treatment in control and YAP1-silenced BC cells were also investigated. Our data showed that E2 modulates YAP1 expression by hypomethylation of its promoter region via downregulation of DNA methyltransferase 3B (DNMT3B); an effect that seems to facilitate tumor progression in BC cells. Although the effect of E2 on YAP1 expression was estrogen receptor (ER) dependent, E2 treatment also upregulated YAP1 expression in MDA-MB231 and SKBR3 cells, which are known ER-negative BC cell lines but expresses ERα. Functionally, E2 treatment resulted in increased cell proliferation, decreased apoptosis, cell cycle arrest, and autophagic flux in MCF7 cells. The knockdown of the YAP1 gene reversed these carcinogenic effects of E2 and inhibited E2-induced autophagy. Lastly, we showed that YAP1 is highly expressed and hypomethylated in human BC tissues and that increased YAP1 expression correlates negatively with DNMT3B expression but strongly associated with ER expression. Our data provide the basis for considering screening of YAP1 expression and its promoter methylation status in the diagnosis and prognosis of BC.

Health benefits attributed to 17α-estradiol, a lifespan-extending compound, are mediated through estrogen receptor α

Metabolic dysfunction underlies several chronic diseases, many of which are exacerbated by obesity. Dietary interventions can reverse metabolic declines and slow aging, although compliance issues remain paramount. 17α-estradiol treatment improves metabolic parameters and slows aging in male mice. The mechanisms by which 17α-estradiol elicits these benefits remain unresolved. Herein, we show that 17α-estradiol elicits similar genomic binding and transcriptional activation through estrogen receptor α (ERα) to that of 17β-estradiol. In addition, we show that the ablation of ERα completely attenuates the beneficial metabolic effects of 17α-E2 in male mice. Our findings suggest that 17α-E2 may act through the liver and hypothalamus to improve metabolic parameters in male mice. Lastly, we also determined that 17α-E2 improves metabolic parameters in male rats, thereby proving that the beneficial effects of 17α-E2 are not limited to mice. Collectively, these studies suggest ERα may be a drug target for mitigating chronic diseases in male mammals.

Assessing the Ability of Developmentally Precocious Estrogen Signaling to Recapitulate Ovarian Transcriptomes and Follicle Dynamics in Alligators from a Contaminated Lake

BACKGROUND

Concern has grown in recent decades over anthropogenic contaminants that interfere with the functioning of endocrine hormones. However, mechanisms connecting developmental processes to pathologies associated with endocrine-disrupting chemical (EDC) exposure are poorly understood in naturally exposed populations.

OBJECTIVES

We sought to a) characterize divergence in ovarian transcriptomic and follicular profiles between alligators originating from a historically EDC-contaminated site, Lake Apopka, and a reference site; b) test the ability of developmentally precocious estrogen exposure to recapitulate site-associated patterns of divergence; and c) test whether treatment with exogenous follicle-stimulating hormone (FSH) is capable of rescuing phenotypes associated with contaminant exposure and/or embryonic estrogen treatment.

METHODS

Alligators eggs were collected from a contaminated site and a reference site, and a subset of eggs from the reference site were treated with estradiol (E2) during embryonic development prior to gonadal differentiation. After hatching, alligators were raised under controlled laboratory settings for 5 months. Juveniles from both sites were divided and treated with exogenous FSH. Histological analyses and RNA-sequencing were conducted to characterize divergence in ovarian follicle dynamics and transcriptomes between sites, between reference and E2-treated animals, and between FSH-treated and nontreated animals.

RESULTS

We observed broad site-of-origin divergence in ovarian transcriptomes and reductions in ovarian follicle density between juvenile alligators from Lake Apopka and the reference site. Treating embryos from the reference site with E2 overwhelmingly recapitulated transcriptional and histological alterations observed in Lake Apopka juveniles. Ovarian phenotypes observed in Lake Apopka alligators or resulting from estrogen treatment were only partially rescued by treatment with exogenous FSH.

DISCUSSION

Recapitulation of ovarian abnormalities by precocious E2 revealed a relatively simple mechanism underlying contaminant-induced pathologies in a historical example of environmental endocrine disruption. Findings reported here support a model where the developmental timing of estrogen signaling has the potential to permanently alter ovarian organization and function. https://doi.org/10.1289/EHP6627.

β-catenin activates TGF-β-induced epithelial-mesenchymal transition in adenomyosis

Adenomyosis is defined as the presence of ectopic nests of endometrial glands and stroma within the myometrium. Adenomyosis is a common cause of dysmenorrhea, menorrhagia, and chronic pelvic pain but is often underdiagnosed. Despite its prevalence and severity of symptoms, its pathogenesis and etiology are poorly understood. Our previous study showed that aberrant activation of β-catenin results in adenomyosis through epithelial–mesenchymal transition. Using transcriptomic and ChIP-seq analysis, we identified activation of TGF-β signaling in the uteri of mutant mice that expressed dominant stabilized β-catenin in the uterus. There was a strong positive correlation between β-catenin and TGF-β2 proteins in women with adenomyosis. Furthermore, treatment with pirfenidone, a TGF-β inhibitor, increased E-cadherin expression and reduced cell invasiveness in Ishikawa cells with nuclear β-catenin. Our results suggest that β-catenin activates TGF-β-induced epithelial–mesenchymal transition in adenomyosis. This finding describes the molecular pathogenesis of adenomyosis and the use of TGF-β as a potential therapeutic target for adenomyosis.

A regulatory link between two proteins involved in the progression of a debilitating uterine condition highlights a potential therapeutic target. Adenomyosis involves the invasion of cells from the inner lining of the uterus (the endometrium) into the uterine muscle wall (the myometrium), resulting in heavy, prolonged periods and chronic pain. The aberrent activation of a protein called β-catenin triggers adenomyosis, but the precise mechanisms are unclear. A team led by Jung-Ho Shin at the Korea University Medical Center, Seoul, South Korea, and Jae-Wook Jeong, Michigan State University, Grand Rapids, USA, used sequencing techniques in mice and human tissue samples to identify the pathways governed by β-catenin in adenomyosis. They found that the Tgf-β2 gene is directly regulated by β-catenin in the uterus. TGF-β2 levels were elevated in human adenomyosis lesions, suggesting the protein could be a therapeutic target.

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

Estrogen receptors (ER) are activated by the steroid hormone 17β-estradiol. Estrogen receptor alpha (ER-α) forms a regulatory network in mammary epithelial cells and in breast cancer with the transcription factors FOXA1 and GATA3. GATA3 is one of the most frequently mutated genes in breast cancer and is capable of specifying chromatin localization of FOXA1 and ER-α. How GATA3 mutations found in breast cancer impact genomic localization of ER-α and the transcriptional network downstream of ER-α and FOXA1 remains unclear. Here, we investigate the function of a recurrent patient-derived GATA3 mutation (R330fs) on this regulatory network. Genomic analysis indicates that the R330fs mutant can disrupt localization of ER-α and FOXA1. Loci co-bound by all three factors are enriched for genes integral to mammary gland development as well as epithelial cell biology. This gene set is differentially regulated in GATA3 mutant cells in culture and in tumors bearing similar mutations in vivo. The altered distribution of ER-α and FOXA1 in GATA3-mutant cells is associated with altered chromatin architecture, which leads to differential gene expression. These results suggest an active role for GATA3 zinc finger 2 mutants in ER-α positive breast tumors.

Estrogen-ERα signaling and DNA hypomethylation co-regulate expression of stem cell protein PIWIL1 in ERα-positive endometrial cancer cells

Background

We previously identified PIWIL1 as an oncogene involved in endometrial carcinogenesis. However, the mechanism of Piwil1 mediated regulation of tumorigenesis remains poorly understood.

Methods

The expression levels of target genes in endometrial cancer cells were detected by quantitative reverse transcription-PCR (RT-qPCR) and western blotting. Up- or down-regulation of ERα or PIWIL1 was achieved by transient transfection with expressing plasmids or short hairpin RNA (shRNA). Dual-luciferase reporter assays and chromatin immunoprecipitation (ChIP) were used to demonstrate the ERα bound to the half estrogen response element (half-ERE) located in PIWIL1 promoter. The expression of PIWIL1 and ERα in endometrial carcinoma tissues were investigated using immunohistochemistry and RT-qPCR. The proliferation ability of cancer cells were evaluated by MTT. Methylation status of the PIWIL1 promoter was detected by bisulfite sequencing PCR (BSP).

Results

In the present study, we found that PIWIL1 mediated E₂-stimulated cancer cell proliferation. In ERα-positive endometrial cancer cells, we demonstrated that estrogen-ERα signaling significantly up-regulated the expression of PIWIL1, which was mediated by binding of the ERα onto the PIWIL1 promoter. Furthermore, we found that a half-ERE in the PIWIL1 promoter was essential for ERα binding. The PIWIL1 promoter was hypomethylated in ERα-positive endometrial cancer cells. Treatment with 5-aza-deoxycytidine (5-aza-dC) could up-regulate PIWIL1 expression.

Conclusions

These findings uncover a novel molecular mechanism by which estrogen-ERα signaling and DNA hypomethylation co-regulate PIWIL1 expression. These findings provide novel insights into the hormonal regulation of PIWIL1 in endometrial cancer and the PIWIL1’s role in estrogen-stimulated endometrial carcinogenesis.

GATA2 and Progesterone Receptor Interaction in Endometrial Stromal Cells Undergoing Decidualization

The transcription factor GATA2 is important for endometrial stromal cell decidualization in early pregnancy. Progesterone receptor (PGR) is also critical during decidualization but its interaction with GATA2 in regulating genes and pathways necessary for decidualization in human endometrium are unclear. RNA-sequencing (RNA-seq) was performed to compare gene expression profiles (n = 3), and chromatin immunoprecipitation followed by sequencing (ChIP-seq) using an antibody against GATA2 (n = 2) was performed to examine binding to target genes in human endometrial stromal cells undergoing in vitro decidualization (IVD including estrogen, progestin, and 3',5'-cyclic AMP analogue) or vehicle treatment. We identified 1232 differentially expressed genes (DEGs) in IVD vs vehicle. GATA2 cistrome in IVD-treated cells was enriched with motifs for GATA, ATF, and JUN, and gene ontology analysis of GATA2 cistrome revealed pathways that regulate cholesterol storage, p38 mitogen-activated protein kinase, and the c-Jun N-terminal kinase cascades. Integration of RNA-seq and ChIP-seq data revealed that the PGR motif is highly enriched at GATA2 binding regions surrounding upregulated genes in IVD-treated cells. The integration of a mined public PGR cistrome in IVD-treated human endometrial cells with our GATA2 cistrome showed that GATA2 binding was significantly enhanced at PGR-binding regions in IVD vs vehicle. Interrogating 2 separate ChIP-seq data sets together with RNA-seq revealed integration of GATA2 and PGR action to coregulate biologic processes during decidualization of human endometrial stromal cells, specifically via WNT activation and stem cell differentiation pathways. These findings reveal the key pathways that are coactivated by GATA2 and PGR that may be therapeutic targets for supporting implantation and early pregnancy.

The Estrogen-Responsive Transcriptome of Female Secondary Sexual Traits in the Gulf Pipefish

Sexual dimorphism often results from hormonally regulated trait differences between the sexes. In sex-role-reversed vertebrates, females often have ornaments used in mating competition that are expected to be under hormonal control. Males of the sex-role-reversed Gulf pipefish (Syngnathus scovelli) develop female-typical traits when they are exposed to estrogens. We aimed to identify genes whose expression levels changed during the development and maintenance of female-specific ornaments. We performed RNA-sequencing on skin and muscle tissue in male Gulf pipefish with and without exposure to estrogen to investigate the transcriptome of the sexually dimorphic ornament of vertical iridescent bands found in females and estrogen-exposed males. We further compared differential gene expression patterns between males and females to generate a list of genes putatively involved in the female secondary sex traits of bands and body depth. A detailed analysis of estrogen-receptor binding sites demonstrates that estrogen-regulated genes tend to have nearby cis-regulatory elements. Our results identified a number of genes that differed between the sexes and confirmed that many of these were estrogen-responsive. These estrogen-regulated genes may be involved in the arrangement of chromatophores for color patterning, as well as in the growth of muscles to achieve the greater body depth typical of females in this species. In addition, anaerobic respiration and adipose tissue could be involved in the rigors of female courtship and mating competition. Overall, this study generates a number of interesting hypotheses regarding the genetic basis of a female ornament in a sex-role-reversed pipefish.

Annotation and cluster analysis of long noncoding RNA linked to male sex and estrogen in cancers

The sex difference in cancer occurrence is a consistent finding in cancer epidemiology. Several solid tumors, including lung cancer, colorectal cancer, hepatic carcinoma, and renal carcinoma, are generally more common in males. Although sexual dimorphism is attributed to hormonal or behavioral differences, evidence for the function of lncRNA is lacking in sex-specific cancers. We show here that LINC00263 is one of the most dysregulated lncRNAs in lung adenocarcinomas and is upregulated in lung adenocarcinoma, colorectal cancer, and renal carcinoma, especially in male patients compared to females. LINC00263 functions as an oncogene by promoting translocation of p65 into the nucleus to activate the NF-κB-signaling pathway through interaction with IKKα in the cytoplasm. The expression of LINC00263 is strongly correlated with ESR1, and it is decreased after treatment with estrogen. Ligand-activated ER could inhibit the function of LINC00263 by inhibiting NF-κB from cytoplasmic translocation into the nucleus. The inhibitory effect of estrogen on LINC00263 indicates its differential expression in male and female patients. Our findings indicate that LINC00263 is linked to male sex and estrogen as an oncogene, and these findings might help in the exploration of the mechanisms of differential gene regulation in sex-specific cancers.

Ret Receptor Has Distinct Alterations and Functions in Breast Cancer

Ret receptor tyrosine kinase is a proto-oncogene that participates in development of various cancers. Several independent studies have recently identified Ret as a key player in breast cancer. Although Ret overexpression and function have been under investigation, mainly in estrogen receptor positive breast cancer, a more comprehensive analysis of the impact of recurring Ret alterations in breast cancer is needed. This review consolidates the current knowledge of Ret alterations and their potential effects in breast cancer. We discuss and integrate data on Ret changes in different breast cancer subtypes and potential function in progression, as well as the participation of distinct Ret network signaling partners in these processes. We propose that it will be essential to define a shared molecular feature of tumors with alteration in Ret receptor, be this at the genetic level or via overexpression in order to design effective therapies to target the Ret pathway. Here we review experimental evidence from basic research and pre-clinical studies concentrating on Ret alterations as potential biomarkers for recurrence, and we discuss the possibility that targeting the Ret pathway might in the future become a treatment for breast cancer.

From Influenza Virus Infections to Lupus: Synchronous Estrogen Receptor α and RNA Polymerase II Binding Within the Immunoglobulin Heavy Chain Locus

Males and females respond to pathogens differently and exhibit significantly different frequencies of autoimmune disease. For example, vaccinated adult females control influenza virus better than males, but females suffer systemic lupus erythematosus at a 9:1 frequency compared to males. Numerous explanations have been offered for these sex differences, but most have involved indirect mechanisms by which estrogen, a nuclear hormone, modifies cell barriers or immunity. In search of a direct mechanism, we examined the binding of estrogen receptor α (ERα), a class I nuclear hormone receptor, to the immunoglobulin heavy chain locus. Here, we show that in purified murine B cells, ERα and RNA polymerase II (RNA Pol II) exhibit extraordinarily similar DNA binding patterns. We further demonstrate that ERα preferentially binds adenosine-cytidine (AC)-repeats in the immunoglobulin heavy chain locus when supplemental estrogen is added to purified, lipopolysaccharide-activated B cells. Based on these and previous data, we hypothesize that (i) estrogen guides the binding of ERα and its RNA Pol II partner within the locus, which in turn instructs sterile transcription and class switch recombination (CSR), (ii) ERα binding to AC-repeats modifies the DNA architecture and loops associated with CSR, and (iii) by these mechanisms, estrogen instructs antibody expression. By targeting ERα-DNA interactions in the immunoglobulin heavy chain locus, clinicians may ultimately enhance antibody responses in the context of infectious diseases and reduce antibody responses in the context of allergic or autoimmune reactions.

Sex steroids as mediators of phenotypic integration, genetic correlations, and evolutionary transitions

In recent decades, endocrinologists have increasingly adopted evolutionary methods and perspectives to characterize the evolution of the vertebrate endocrine system and leverage it as a model for developing and testing evolutionary theories. This review summarizes recent research on sex steroids (androgens and estrogens) to illustrate three ways in which a detailed understanding of the molecular and cellular architecture of hormonally mediated gene expression can enhance our understanding of general evolutionary principles. By virtue of their massively pleiotropic effects on the expression of genes and phenotypes, sex steroids and their receptors can (1) structure the patterns of phenotypic variance and covariance that are available to natural selection, (2) alter the underlying genetic correlations that determine a population's evolutionary response to selection, and (3) facilitate evolutionary transitions in fitness-related phenotypes via subtle regulatory shifts in underlying tissues and genes. These principles are illustrated by the author's research on testosterone and sexual dimorphism in lizards, and by recent examples drawn from other vertebrate systems. Mechanistically, these examples call attention to the importance of evolutionary changes in (1) androgen- and estrogen-mediated gene expression, (2) androgen and estrogen receptor expression, and (3) the distribution of androgen and estrogen response elements in target genes throughout the genome. A central theme to emerge from this review is that the rapidly increasing availability of genomic and transcriptomic data from non-model organisms places evolutionary endocrinologist in an excellent position to address the hormonal regulation of the key evolutionary interface between genes and phenotypes.

Complex sex-biased antibody responses: estrogen receptors bind estrogen response elements centered within immunoglobulin heavy chain gene enhancers

Nuclear hormone receptors including the estrogen receptor (ERα) and the retinoic acid receptor regulate a plethora of biological functions including reproduction, circulation and immunity. To understand how estrogen and other nuclear hormones influence antibody production, we characterized total serum antibody isotypes in female and male mice of C57BL/6J, BALB/cJ and C3H/HeJ mouse strains. Antibody levels were higher in females compared to males in all strains and there was a female preference for IgG2b production. Sex-biased patterns were influenced by vitamin levels, and by antigen specificity toward influenza virus or pneumococcus antigens. To help explain sex biases, we examined the direct effects of estrogen on immunoglobulin heavy chain sterile transcript production among purified, lipopolysaccharide-stimulated B cells. Supplemental estrogen in B-cell cultures significantly increased immunoglobulin heavy chain sterile transcripts. Chromatin immunoprecipitation analyses of activated B cells identified significant ERα binding to estrogen response elements (EREs) centered within enhancer elements of the immunoglobulin heavy chain locus, including the Eµ enhancer and hypersensitive site 1,2 (HS1,2) in the 3' regulatory region. The ERE in HS1,2 was conserved across animal species, and in humans marked a site of polymorphism associated with the estrogen-augmented autoimmune disease, lupus. Taken together, the results highlight: (i) the important targets of ERα in regulatory regions of the immunoglobulin heavy chain locus that influence antibody production, and (ii) the complexity of mechanisms by which estrogen instructs sex-biased antibody production profiles.

The Impact of Skeletal Muscle ERα on Mitochondrial Function and Metabolic Health

The incidence of chronic disease is elevated in women after menopause. Increased expression of ESR1 (the gene that encodes the estrogen receptor alpha, ERα) in muscle is highly associated with metabolic health and insulin sensitivity. Moreover, reduced muscle expression levels of ESR1 are observed in women, men, and animals presenting clinical features of the metabolic syndrome (MetSyn). Considering that metabolic dysfunction elevates chronic disease risk, including type 2 diabetes, heart disease, and certain cancers, treatment strategies to combat metabolic dysfunction and associated pathologies are desperately needed. This review will provide published work supporting a critical and protective role for skeletal muscle ERα in the regulation of mitochondrial function, metabolic homeostasis, and insulin action. We will provide evidence that muscle-selective targeting of ERα may be effective for the preservation of mitochondrial and metabolic health. Collectively published findings support a compelling role for ERα in the control of muscle metabolism via its regulation of mitochondrial function and quality control. Studies identifying ERα-regulated pathways essential for disease prevention will lay the important foundation for the design of novel therapeutics to improve metabolic health of women while limiting secondary complications that have historically plagued traditional hormone replacement interventions.

Cooperativity of co-factor NR2F2 with Pioneer Factors GATA3, FOXA1 in promoting ERα function

Estrogen receptor α (ERα) drives growth in the majority of human breast cancers by binding to regulatory elements and inducing transcriptional events that promote tumor growth. ERα binding activity largely depends on access to binding sites on chromatin, which is facilitated in part by Pioneer Factors (PFs). Transcription factors operate in complexes through thousands of genomic binding sites in a combinatorial fashion to control the expression of genes. However, the extent of crosstalk and cooperation between ERα pioneer factors and more collaborative transcription factors in breast cancer still remains to be elucidated systematically. Methods: Here, we determined the genomic binding information of 40 transcription-related factors and histone modifications with ChIP-seq in ENCODE and integrated it with other genomic information (RNA-seq, ATAC-seq, Gene microarray, 450k methylation chip, GRO-seq), forming a multi-dimension network to illuminate ERα associated transcription. Results: We show that transcription factor, NR2F2 binds to most sites independently of estrogen. Perturbation of NR2F2 expression decreases ERα DNA binding, chromatin openning, and estrogen-dependent cell growth. In the genome-wide analysis, we show that most binding events of NR2F2 and known pioneer factors FOXA1, GATA3 occur together, covering 85% of the ERα binding sites. Regions bound by all the three TFs appeared to be the most active, to have the strongest ERα binding and to be enriched for the super enhancers. Conclusions: The ERα binds to pre-accessible sites containing ERE elements bound by the three transcription factors (NR2F2, FOXA1 and GATA3).The three genes were also identified to correlate with decreased metastatic potential in patient cohorts and co-regulate each other. Together, our results suggest that NR2F2 is a cofactor with FOXA1 and GATA3 in ERα-mediated transcription.

The genomic landscape of estrogen receptor α binding sites in mouse mammary gland

Estrogen receptor α (ERα) is the major driving transcription factor in the mammary gland development as well as breast cancer initiation and progression. However, the genomic landscape of ERα binding sites in the normal mouse mammary gland has not been completely elucidated. Here, we mapped genome-wide ERα binding events by chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) in the mouse mammary gland in response to estradiol. We identified 6237 high confidence ERα binding sites in two biological replicates and showed that many of these were located at distal enhancer regions. Furthermore, we discovered 3686 unique genes in the mouse genome that recruit ER in response to estradiol. Interrogation of ER-DNA binding sites in ER-positive luminal epithelial cells showed that the ERE, PAX2, SF1, and AP1 motifs were highly enriched at distal enhancer regions. In addition, comprehensive transcriptome analysis by RNA-seq revealed that 493 genes are differentially regulated by acute treatment with estradiol in the mouse mammary gland in vivo. Through integration of RNA-seq and ERα ChIP-seq data, we uncovered a novel ERα targetome in mouse mammary epithelial cells. Taken together, our study has identified the genomic landscape of ERα binding events in mouse mammary epithelial cells. Furthermore, our study also highlights the cis-regulatory elements and cofactors that are involved in estrogen signaling and may contribute to ductal elongation in the normal mouse mammary gland.

Therapeutic Ligands Antagonize Estrogen Receptor Function by Impairing Its Mobility

Estrogen receptor-positive (ER⁺) breast cancers frequently remain dependent on ER signaling even after acquiring resistance to endocrine agents, prompting the development of optimized ER antagonists. Fulvestrant is unique among approved ER therapeutics due to its capacity for full ER antagonism, thought to be achieved through ER degradation. The clinical potential of fulvestrant is limited by poor physicochemical features, spurring attempts to generate ER degraders with improved drug-like properties. We show that optimization of ER degradation does not guarantee full ER antagonism in breast cancer cells; ER "degraders" exhibit a spectrum of transcriptional activities and anti-proliferative potential. Mechanistically, we find that fulvestrant-like antagonists suppress ER transcriptional activity not by ER elimination, but by markedly slowing the intra-nuclear mobility of ER. Increased ER turnover occurs as a consequence of ER immobilization. These findings provide proof-of-concept that small molecule perturbation of transcription factor mobility may enable therapeutic targeting of this challenging target class.

erefinder: Genome-wide detection of oestrogen response elements

Oestrogen response elements (EREs) are specific DNA sequences to which ligand-bound oestrogen receptors (ERs) physically bind, allowing them to act as transcription factors for target genes. Locating EREs and ER responsive regions is therefore a potentially important component of the study of oestrogen-regulated pathways. Here, we report the development of a novel software tool, erefinder, which conducts a genome-wide, sliding-window analysis of oestrogen receptor binding affinity. We demonstrate the effects of adjusting window size and highlight the program's general agreement with ChIP studies. We further provide two examples of how erefinder can be used for comparative approaches. erefinder can handle large input files, has settings to allow for broad and narrow searches, and provides the full output to allow for greater data manipulation. These features facilitate a wide range of hypothesis testing for researchers and make erefinder an excellent tool to aid in oestrogen-related research.

Estrogen receptor α regulates the expression of syndecan-1 in human breast carcinoma cells

Breast cancer (BC) is the primary cause of cancer-related mortality among women. Patients who express the estrogen receptor (ER), which mediates the tumorigenic effects of estrogens, respond to antihormonal therapy. Loss of ER expression or acquired resistance to E2 is associated with aggressive malignant phenotypes, which lead to relapse. These BC subtypes overexpress syndecan-1 (SDC1), a transmembrane heparan sulfate proteoglycan that mediates angiogenesis as well as the proliferation and invasiveness of cancer cells. We showed here that the activation of ER-alpha (ERα) by estrogens induces downregulation of SDC1 expression in ER(+) MCF7 cells but not in T47D cells. Loss of ERα expression, induced by RNA interference or a selective ER downregulator, led to subsequent SDC1 overexpression. E2-dependent downregulation of SDC1 expression required de novo protein synthesis and was antagonized by treatment with BAY 11-7085, an irreversible inhibitor of IκBα phosphorylation, which inhibits the activation of NFκB. Downregulation of SDC1 expression required ERα and activation of IKK, but was independent to downstream transcriptional regulators of NFκB. BAY 11-7085 prevented E2-mediated phosphorylation of ERα on Ser118, increasing its proteasomal degradation, suggesting that IKK stabilized E2-activated ERα, leading to subsequent downregulation of SDC1 expression. Our results showed that sustained ER signaling inhibits SDC1 expression. Such antagonism elucidates the inverse correlation between SDC1 and ER expression in ER(+) BC as well as the overexpression of SDC1 in hormone receptor-negative BC subtypes with the most aggressive phenotypes. These results identify SDC1 as an attractive therapeutic target for BC as well as for other endocrine-associated cancers.

Temporal dynamic reorganization of 3D chromatin architecture in hormone-induced breast cancer and endocrine resistance

Recent studies have demonstrated that chromatin architecture is linked to the progression of cancers. However, the roles of 3D structure and its dynamics in hormone-dependent breast cancer and endocrine resistance are largely unknown. Here we report the dynamics of 3D chromatin structure across a time course of estradiol (E2) stimulation in human estrogen receptor α (ERα)-positive breast cancer cells. We identified subsets of temporally highly dynamic compartments predominantly associated with active open chromatin and found that these highly dynamic compartments showed higher alteration in tamoxifen-resistant breast cancer cells. Remarkably, these compartments are characterized by active chromatin states, and enhanced ERα binding but decreased transcription factor CCCTC-binding factor (CTCF) binding. We finally identified a set of ERα-bound promoter-enhancer looping genes enclosed within altered domains that are enriched with cancer invasion, aggressiveness or metabolism signaling pathways. This large-scale analysis expands our understanding of high-order temporal chromatin reorganization underlying hormone-dependent breast cancer.

Comprehensive Genomic Analysis Reveals that the Pioneering Function of FOXA1 Is Independent of Hormonal Signaling

Considerable work has linked hormone receptors, such as estrogen receptor-alpha (ER), with the pioneer factor FOXA1. Altered FOXA1 levels contribute to endocrine-resistant breast cancer, where it maintains ER-chromatin interactions, even in contexts in which cells are refractory to ER-targeted drugs. A recent study controversially suggests that FOXA1 binding can be induced by hormonal pathways, including the estrogen-ER complex. We now show that the vast majority (>99%) of FOXA1 binding events are unaffected by steroid activation. A small number (<1%) of FOXA1 binding sites appear to be induced by estrogen, but these are created from chromatin interactions between ER binding sites and adjacent FOXA1 binding sites and do not represent genuine new FOXA1-pioneering elements. FOXA1 is therefore not regulated by estrogen and remains a bone fide pioneer factor that is entirely upstream of the ER complex.

Estrogen Signaling in Endometrial Cancer: a Key Oncogenic Pathway with Several Open Questions

Endometrial cancer is the most common gynecological cancer in the developed world, and it is one of the few cancer types that is becoming more prevalent and leading to more deaths in the USA each year. The majority of endometrial tumors are considered to be hormonally driven, where estrogen signaling through estrogen receptor α (ER) acts as an oncogenic signal. The major risk factors and some treatment options for endometrial cancer patients emphasize a key role for estrogen signaling in the disease. Despite the strong connections between estrogen signaling and endometrial cancer, important molecular aspects of ER function remain poorly understood; however, progress is being made in our understanding of estrogen signaling in endometrial cancer. Here, we discuss the evidence for the importance of estrogen signaling in endometrial cancer, details of the endometrial cancer-specific actions of ER, and open questions surrounding estrogen signaling in endometrial cancer.

The Phosphorylated Estrogen Receptor α (ER) Cistrome Identifies a Subset of Active Enhancers Enriched for Direct ER-DNA Binding and the Transcription Factor GRHL2

Posttranslational modifications are key regulators of protein function, providing cues that can alter protein interactions and cellular location. Phosphorylation of estrogen receptor α (ER) at serine 118 (pS118-ER) occurs in response to multiple stimuli and is involved in modulating ER-dependent gene transcription. While the cistrome of ER is well established, surprisingly little is understood about how phosphorylation impacts ER-DNA binding activity. To define the pS118-ER cistrome, chromatin immunoprecipitation sequencing was performed on pS118-ER and ER in MCF-7 cells treated with estrogen. pS118-ER occupied a subset of ER binding sites which were associated with an active enhancer mark, acetylated H3K27. Unlike ER, pS118-ER sites were enriched in GRHL2 DNA binding motifs, and estrogen treatment increased GRHL2 recruitment to sites occupied by pS118-ER. Additionally, pS118-ER occupancy sites showed greater enrichment of full-length estrogen response elements relative to ER sites. In an in vitro DNA binding array of genomic binding sites, pS118-ER was more commonly associated with direct DNA binding events than indirect binding events. These results indicate that phosphorylation of ER at serine 118 promotes direct DNA binding at active enhancers and is a distinguishing mark for associated transcription factor complexes on chromatin.