Multifaceted Transcriptional Network of Estrogen-Related Receptor Alpha in Health and Disease.

Estrogen-related receptors (ERRα, ß and γ in mammals) are orphan members of the nuclear receptor superfamily acting as transcription factors. ERRs are expressed in several cell types and they display various functions in normal and pathological contexts. Amongst others, they are notably involved in bone homeostasis, energy metabolism and cancer progression. In contrast to other nuclear receptors, the activities of the ERRs are apparently not controlled by a natural ligand but they rely on other means such as the availability of transcriptional co-regulators. Here we focus on ERRα and review the variety of co-regulators that have been identified by various means for this receptor and their reported target genes. ERRα cooperates with distinct co-regulators to control the expression of distinct sets of target genes. This exemplifies the combinatorial specificity of transcriptional regulation that induces discrete cellular phenotypes depending on the selected coregulator. We finally propose an integrated view of the ERRα transcriptional network.

ERRα induces H3K9 demethylation by LSD1 to promote cell invasion.

Lysine Specific Demethylase 1 (LSD1) removes mono- and dimethyl groups from lysine 4 of histone H3 (H3K4) or H3K9, resulting in repressive or activating (respectively) transcriptional histone marks. The mechanisms that control the balance between these two antagonist activities are not understood. We here show that LSD1 and the orphan nuclear receptor estrogen-related receptor α (ERRα) display commonly activated genes. Transcriptional activation by LSD1 and ERRα involves H3K9 demethylation at the transcriptional start site (TSS). Strikingly, ERRα is sufficient to induce LSD1 to demethylate H3K9 in vitro. The relevance of this mechanism is highlighted by functional data. LSD1 and ERRα coregulate several target genes involved in cell migration, including the MMP1 matrix metallo-protease, also activated through H3K9 demethylation at the TSS. Depletion of LSD1 or ERRα reduces the cellular capacity to invade the extracellular matrix, a phenomenon that is rescued by MMP1 reexpression. Altogether our results identify a regulatory network involving a direct switch in the biochemical activities of a histone demethylase, leading to increased cell invasion.

Computational identification of potential transcriptional regulators of TGF-ß1 in human atherosclerotic arteries.

TGF-ß is protective in atherosclerosis but deleterious in metastatic cancers. Our aim was to determine whether TGF-ß transcriptional regulation is tissue-specific in early atherosclerosis. The computational methods included 5 steps: (i) from microarray data of human atherosclerotic carotid tissue, to identify the 10 best co-expressed genes with TGFB1 (TGFB1 gene cluster), (ii) to choose the 11 proximal promoters, (iii) to predict the TFBS shared by the promoters, (iv) to identify the common TFs co-expressed with the TGFB1 gene cluster, and (v) to compare the common TFs in the early lesions to those identified in advanced atherosclerotic lesions and in various cancers. Our results show that EGR1, SP1 and KLF6 could be responsible for TGFB1 basal expression, KLF6 appearing specific to atherosclerotic lesions. Among the TFs co-expressed with the gene cluster, transcriptional activators (SLC2A4RG, MAZ) and repressors (ZBTB7A, PATZ1, ZNF263) could be involved in the fine-tuning of TGFB1 expression in atherosclerosis.

NRF2 Shortage in Human Skin Fibroblasts Dysregulates Matrisome Gene Expression and Affects Collagen Fibrillogenesis.

NRF2 is a master regulator of the antioxidative response that was recently proposed as a potential regulator of extracellular matrix (ECM) gene expression. Fibroblasts are major ECM producers in all connective tissues, including the dermis. A better understanding of NRF2-mediated ECM regulation in skin fibroblasts is thus of great interest for skin homeostasis maintenance and aging protection. In this study, we investigate the impact of NRF2 downregulation on matrisome gene expression and ECM deposits in human primary dermal fibroblasts. RNA-sequencing‒based transcriptome analysis of NRF2 silenced dermal fibroblasts shows that ECM genes are the most regulated gene sets, highlighting the relevance of the NRF2-mediated matrisome program in these cells. Using complementary light and electron microscopy methods, we show that NRF2 deprivation in dermal fibroblasts results in reduced collagen I biosynthesis and impacts collagen fibril deposition. Moreover, we identify ZNF469, a putative transcriptional regulator of collagen biosynthesis, as a target of NRF2. Both ZNF469 silenced fibroblasts and fibroblasts derived from Brittle Corneal Syndrome patients carrying variants in ZNF469 gene show reduced collagen I gene expression. Our study shows that NRF2 orchestrates matrisome expression in human skin fibroblasts through direct or indirect transcriptional mechanisms that could be prioritized to target dermal ECM homeostasis in health and disease.

ERRα coordinates actin and focal adhesion dynamics.

Cell migration depends on the dynamic organisation of the actin cytoskeleton and assembly and disassembly of focal adhesions (FAs). However, the precise mechanisms coordinating these processes remain poorly understood. We previously identified the oestrogen-related receptor α (ERRα) as a major regulator of cell migration. Here, we show that loss of ERRα leads to abnormal accumulation of actin filaments that is associated with an increased level of inactive form of the actin-depolymerising factor cofilin. We further show that ERRα depletion decreases cell adhesion and results in defective FA formation and turnover. Interestingly, specific inhibition of the RhoA-ROCK-LIMK-cofilin pathway rescues the actin polymerisation defects resulting from ERRα silencing, but not cell adhesion. Instead, we found that MAP4K4 is a direct target of ERRα and down-regulation of its activity rescues cell adhesion and FA formation in the ERRα-depleted cells. Altogether, our results highlight a crucial role of ERRα in coordinating the dynamic of actin network and FAs through the independent regulation of the RhoA and MAP4K4 pathways.

Computational identification of new potential transcriptional partners of ERRα in breast cancer cells: specific partners for specific targets.

Estrogen related receptors are orphan members of the nuclear receptor superfamily acting as transcription factors (TFs). In contrast to classical nuclear receptors, the activities of the ERRs are not controlled by a natural ligand. Regulation of their activities thus relies on availability of transcriptional co-regulators. In this paper, we focus on ERRα, whose involvement in cancer progression has been broadly demonstrated. We propose a new approach to identify potential co-activators, starting from previously identified ERRα-activated genes in a breast cancer (BC) cell line. Considering mRNA gene expression from two sets of human BC cells as major endpoint, we used sparse partial least squares modeling to uncover new transcriptional regulators associated with ERRα. Among them, DDX21, MYBBP1A, NFKB1, and SETD7 are functionally relevant in MDA-MB-231 cells, specifically activating the expression of subsets of ERRα-activated genes. We studied SET7 in more details and showed its co-localization with ERRα and its ERRα-dependent transcriptional and phenotypic effects. Our results thus demonstrate the ability of a modeling approach to identify new transcriptional partners from gene expression. Finally, experimental results show that ERRα cooperates with distinct co-regulators to control the expression of distinct sets of target genes, thus reinforcing the combinatorial specificity of transcription.

MEN1 silencing triggers the dysregulation of mTORC1 and MYC pathways in ER+ breast cancer cells.

Menin, encoded by the MEN1 gene, has been identified as a critical factor regulating ESR1 transcription, playing an oncogenic role in ER+ breast cancer (BC) cells. Here, we further dissected the consequences of menin inactivation in ER+ BC cells by focusing on factors within two major pathways involved in BC, mTOR and MYC. MEN1 silencing in MCF7 and T-47D resulted in an increase in phosphor-p70S6K1, phosphor-p85S6K1 and phosphor-4EBP1 expression. The use of an AKT inhibitor inhibited the activation of S6K1 and S6RP triggered by MEN1 knockdown (KD). Moreover, MEN1 silencing in ER+ BC cells led to increased formation of the eIF4E and 4G complex. Clinical studies showed that patients with menin-low breast cancer receiving tamoxifen plus everolimus displayed a trend toward better overall survival. Importantly, MEN1 KD in MCF7 and T-47D cells led to reduced MYC expression. ChIP analysis demonstrated that menin bound not only to the MYC promoter but also to its 5' enhancer. Furthermore, E2-treated MEN1 KD MCF7 cells displayed a decrease in MYC activation, suggesting its role in estrogen-mediated MYC transcription. Finally, expression data mining in tumors revealed a correlation between the expression of MEN1 mRNA and that of several mTORC1 components and targets and a significant inverse correlation between MEN1 and two MYC inhibitory factors, MYCBP2 and MYCT1, in ER+ BC. The current work thus highlights altered mTORC1 and MYC pathways after menin inactivation in ER+ BC cells, providing insight into the crosstalk between menin, mTORC1 and MYC in ER+ BC.

Transforming growth factor-ß1 inhibits interleukin-1ß-induced expression of inflammatory genes and Cathepsin S activity in human vascular smooth muscle cells.

OBJECTIVE AND DESIGN: This study investigated the opposite mechanisms by which IL-1ß and TGF-ß1 modulated the inflammatory and migratory phenotypes in cultured human intimal vascular smooth muscle cells vSMCs. MATERIALS AND TREATMENT: Primary human vSMCs, obtained from twelve hypertensive patients who underwent carotid endarterectomy, were incubated for 24 hours with either 40 pM TGF-ß1, or 1 nmol/L IL-1ß, or their combination in presence or absence of anti-TGF-ß neutralizing antibody. METHODS: The expression levels of matrix metalloproteases and their inhibitors, and the elastolytic enzyme cathepsin S (CTSS) and its inhibitor cystatin C were evaluated with RT-PCR. CTSS activity was measured by fluorometry. RESULTS: TGF-ß1 reversed IL-1ß-induced expression of iNOS, CXCL6, IL1R1, MMP12, and CTSS, while upregulated TIMP2 expression. Furthermore, anti-TGF-ß neutralizing antibody abrogated TGF-ß effects. Combination with IL-1ß and TGF-ß1 induced the expression of IL1α, IL1ß, IL1R1, and CTSS, but suppressed CST3 expression. CTSS expression in the combination treatment was higher than that of cells treated with anti-TGF-ß antibodies alone. Moreover, IL-1ß-induced CTSS enzymatic activity was reduced when human vSMCs were co-treated with TGF-ß, whereas this reduction was abrogated by anti-TGF-ß neutralizing antibody. CONCLUSION: TGF-ß1 abrogated IL-1ß-induced expression of inflammatory genes and elastolytic activity in cultured human vSMCs. Thus, TGF-ß1 can play a crucial role in impairing IL-1ß-induced vascular inflammation and damage involved in the etiology of cardiovascular diseases.

A Cell-Based Method to Detect Agonist and Antagonist Activities of Endocrine-Disrupting Chemicals on GPER.

Endocrine-disrupting chemicals (EDCs) are exogenous compounds that impact endogenous hormonal systems, resulting in adverse health effects. These chemicals can exert their actions by interfering with several pathways. Simple biological systems to determine whether EDCs act positively or negatively on a given receptor are often lacking. Here we describe a low-to-middle throughput method to screen the agonist/antagonist potential of EDCs specifically on the GPER membrane estrogen receptor. Application of this assay to 23 candidate EDCs from different chemical families reveals the existence of six agonists and six antagonists.

Functional meta-analysis of double connectivity in gene coexpression networks in mammals.

In functional genomics, the high-throughput methods such as microarrays 1) allow analysis of the relationships between genes considering them as elements of a network and 2) lead to biological interpretations thanks to Gene Ontology. But up to now it has not been possible to find relationships between the functions and the connectivity of the genes in coexpression networks. To achieve this aim, we have defined a double connectivity for each gene by the numbers of its significant negative and positive correlations with the other genes within a given biological condition, or group. Here, based on the analysis of 1,260 DNA microarrays, we show that this double connectivity clearly separates two types of genes, those with a predominantly strong negative connectivity, hub- genes, and those with a predominantly strong positive connectivity, hub+ genes. Interestingly, the hub+ genes concerned transcription factors more often than by chance and, similarly, for the hub- genes concerning miRNA predicted targets. Furthermore, a meta-analysis of GO annotations carried out on 67 groups in humans and rats shows that these two types of genes correspond to a functional biological duality. The hub- genes were mainly involved in basic functions common to all eukaryote cells, whereas the hub+ genes were mainly involved in specialized functions related to cell differentiation and communication. The separation and the biological role of these hub- and hub+ genes provide a powerful new tool for a better understanding of the control and regulation of the key genes involved in cellular differentiation and physiopathological conditions.

Increased insulin-stimulated expression of arterial angiotensinogen and angiotensin type 1 receptor in patients with type 2 diabetes mellitus and atheroma.

OBJECTIVE: Because inhibition of the renin-angiotensin system (RAS) reduces the onset of type 2 diabetes (T2D) and prevents atherosclerosis, we investigated the expression of RAS in the arterial wall of T2D and nondiabetic (CTR) patients. METHODS AND RESULTS: mRNA and protein levels of angiotensinogen (AGT), angiotensin-converting enzyme (ACE) and AT1 receptor (AT1R) were determined in carotid atheroma plaque, nearby macroscopically intact tissue (MIT), and in vascular smooth muscle cells (VSMCs) before and after insulin stimulation from 21 T2D and 22 CTR patients. AGT and ACE mRNA and their protein levels were 2- to 3-fold higher in atheroma and in MIT of T2D patients. VSMCs from T2D patients had respectively 2.5- and 5-fold higher AGT and AT1R mRNA and protein contents. Insulin induced an increase in AGT and AT1R mRNA with similar ED50. These responses were blocked by PD98059, an inhibitor of MAP-kinase in the two groups whereas wortmannin, an inhibitor of PI3-kinase, partially prevented the response in CTR patients. Phosphorylated ERK1-2 was 4-fold higher in MIT from T2D than from CTR patients. CONCLUSIONS: The arterial RAS is upregulated in T2D patients, which can be partly explained by an hyperactivation of the ERK1-2 pathway by insulin.

LSD1-ERRα complex requires NRF1 to positively regulate transcription and cell invasion.

Lysine-specific demethylase 1 (LSD1) exerts dual effects on histone H3, promoting transcriptional repression via Lys4 (H3K4) demethylation or transcriptional activation through Lys9 (H3K9) demethylation. These activities are often exerted at transcriptional start sites (TSSs) and depend on the type of enhancer-bound transcription factor (TFs) with which LSD1 interacts. In particular, the Estrogen-Receptor Related α (ERRα) TF interacts with LSD1 and switches its activities toward H3K9 demethylation, resulting in transcriptional activation of a set of common target genes. However, how are the LSD1-TF and, in particular LSD1-ERRα, complexes determined to act at TSSs is not understood. Here we show that promoter-bound nuclear respiratory factor 1 (NRF1), but not ERRα, is essential to LSD1 recruitment at the TSSs of positive LSD1-ERRα targets. In contrast to ERRα, NRF1 does not impact on the nature of LSD1 enzymatic activity. We propose a three factor model, in which the LSD1 histone modifier requires a TSS tethering factor (NRF1) as well as an activity inducer (ERRα) to transcriptionally activate common targets. The relevance of this common network is illustrated by functional data, showing that all three factors are required for cell invasion in an MMP1 (Matrix MetalloProtease 1)-dependent manner, the expression of which is regulated by NRF1/LSD1/ERRα-mediated H3K9me2 demethylation.

ERRα protein is stabilized by LSD1 in a demethylation-independent manner.

The LSD1 histone demethylase is highly expressed in breast tumors where it constitutes a factor of poor prognosis and promotes traits of cancer aggressiveness such as cell invasiveness. Recent work has shown that the Estrogen-Related Receptor α (ERRα) induces LSD1 to demethylate the Lys 9 of histone H3. This results in the transcriptional activation of a number of common target genes, several of which being involved in cellular invasion. High expression of ERRα protein is also a factor of poor prognosis in breast tumors. Here we show that, independently of its demethylase activities, LSD1 protects ERRα from ubiquitination, resulting in overexpression of the latter protein. Our data also suggests that the elevation of LSD1 mRNA and protein in breast cancer (as compared to normal tissue) may be a key event to increase ERRα protein, independently of its corresponding mRNA.

Transcriptional alterations in the left ventricle of three hypertensive rat models.

Left ventricular hypertrophy (LVH) is commonly associated with hypertension and represents an independent cardiovascular risk factor. The aim of this study was to test the hypothesis that the cardiac overload related to hypertension is associated to a specific gene expression pattern independently of genetic background. Gene expression levels were obtained with microarrays for 15,866 transcripts from RNA of left ventricles from 12-wk-old rats of three hypertensive models [spontaneously hypertensive rat (SHR), Lyon hypertensive rat (LH), and heterozygous TGR(mRen2)27 rat] and their respective controls. More than 60% of the detected transcripts displayed significant changes between the three groups of normotensive rats, showing large interstrain variability. Expression data were analyzed with respect to hypertension, LVH, and chromosomal distribution. Only four genes had significantly modified expression in the three hypertensive models among which a single gene, coding for sialyltransferase 7A, was consistently overexpressed. Correlation analysis between expression data and left ventricular mass index (LVMI) over all rats identified a larger set of genes whose expression was continuously related with LVMI, including known genes associated with cardiac remodeling. Positioning the detected transcripts along the chromosomes pointed out high-density regions mostly located within blood pressure and cardiac mass quantitative trait loci. Although our study could not detect a unique reprogramming of cardiac cells involving specific genes at early stage of LVH, it allowed the identification of some genes associated with LVH regardless of genetic background. This study thus provides a set of potentially important genes contained within restricted chromosomal regions involved in cardiovascular diseases.

Risk profile in hypertension genesis: A five-year follow-up study.

BACKGROUND: Pathogenesis of primary hypertension remains unclear, because many heterogeneous factors (diet, physiologic, and psychological factors) are simultaneously involved. We have conducted an original analysis to study the influence of the combination of these factors on BP evolution. METHODS: Seven homogeneous clusters were constituted from 213 healthy normotensive subjects taking into account 10 variables. Those variables used to cluster homogeneous "risk profiles" are usually considered as potential risk factors for hypertension: age, body mass index, alcohol consumption, sodium/potassium urinary ratio, systolic blood pressure (SBP) and heart rate response to mental stress, baroreflex sensitivity (BRS), job demand, job latitude, and behavioral pattern (personality score). Five-year BP evolution (DeltaSBP or Deltadiastolic BP [DBP]) was compared between risk profiles. RESULTS: Four clusters of subjects representing about 50% of the population had a significantly higher 5-year DeltaSBP (>or=5 mm Hg) compared to the 5-year DeltaSBP of the two clusters in which SBP did not increase. These four clusters had a low BRS. Two profiles that group six unfavorable risk factors had the most detrimental 5-year DeltaSBP. Interestingly, perceived high job demand in a cluster of younger subjects with a high personality score and a low BRS had also a detrimental SBP evolution. CONCLUSIONS: The major interest of this study is to highlight that hypertension development is not univocal between subjects and that different combinations of factors could explain differential BP evolution in groups of subjects sharing the same risk profile. A lower BRS was a consistent predictor for detrimental 5-year BP evolution.

Is high job strain associated with hypertension genesis?

BACKGROUND: The aim of this analysis was to test, in a large sample of normotensive subjects, the short-term influence of job strain on the onset of hypertension. METHODS: According to the questionnaire of Karasek et al, job strain was divided into four modalities: (high strain, low strain, passive, and active) based on job demand (eg, the need to work hard and quickly) and job latitude (eg, control over skill use, time allocation, and organizational decisions) scores. High strain (HS) was defined by a high demand and a low job decision latitude. Individual data obtained in 926 (age 41 +/- 6 years) healthy normotensive or newly diagnosed hypertensive subjects were analyzed. Subjects participated in two prospective work site surveys designed to assess the influence of job strain on hypertension development. Relationships between job strain modalities and work site blood pressure (BP) levels were assessed using a general linear model. A complementary analysis using the the Pearson Phi coefficient (Z analysis) was implemented to explore nonlinear or scattered relationships between job strain and onset of hypertension. RESULTS: Systolic BP (SBP) was linearly related significantly to BMI and alcohol consumption, whereas diastolic BP (DBP) was related to age. The linear model did not find any relationship between SBP or DBP and job strain modalities. Using the Z analysis, development of systolic hypertension (SBP >140 mm Hg) was significantly associated with high job strain (P < .001). CONCLUSIONS: Our results suggest that there is no global relationship between job strain and BP levels. However our methodology revealed a significant association between job strain and work site BP in a predominantly male subgroup of newly diagnosed hypertensive subjects exposed to high job strain.

Transcriptomic Analysis Reveals Novel Transcription Factors Associated With Renin-Angiotensin-Aldosterone System in Human Atheroma.

Despite the well-known role of the renin-angiotensin-aldosterone system (RAAS) in atheroma, its global local organization is poorly understood. In this study, we used transcriptomic meta-analysis to reveal the local transcriptional organization and regulation of 37 extended RAAS (extRAAS) genes in atheroma. Expression analysis and hierarchical clustering were done on extRAAS genes in 32 paired early and advanced atherosclerotic lesions. Contrary to receptor-coding transcripts, multiple angiotensin-metabolizing enzymes showed higher expression in advance, in comparison to early lesions. Interestingly, similar results were obtained from GEO data sets containing human (n=839) and mouse (n=18) atherosclerotic samples, but different from normal human (n=11) arterial tissues. The expression and coordination patterns were then used to construct transcriptional maps of extRAAS, displaying favored pathways in atheroma. Three coexpression modules (M1, M2, and M3) with >80% reproducibility across human atheroma data sets were identified. M1 and M3 contained angiotensin-metabolizing enzymes transcripts, whereas M2 contained proatherogenic receptor-coding transcripts. Interestingly, M1 and M2 were negatively correlated. A total of 21 transcription factors with enriched binding sites in the promoters of coordinated genes were extracted, among which IRF5, MAX, and ETV5 showed significant positive correlations with M1, but negative correlations with M2. However, ETS1 and SMAD1 transcripts were positively correlated to receptor-coding genes in M2. Despite sharing some similarities in extRAAS organization with kidney and adipose, atheroma showed specific correlations between extRAAS and transcription factors. In conclusion, our transcriptional map helps in designing more efficient treatments for atherosclerosis. In addition, the identified transcription factors provide a basis for the discovery of atheroma-specific modulators of extRAAS.

miR-135a Inhibits the Invasion of Cancer Cells via Suppression of ERRα.

MicroRNA-135a (miR-135a) down-modulates parameters of cancer progression and its expression is decreased in metastatic breast cancers (as compared to non-metastatic tumors) as well as in prostate tumors relative to normal tissue. These expression and activity patterns are opposite to those of the Estrogen-Related Receptor α (ERRα), an orphan member of the nuclear receptor family. Indeed high expression of ERRα correlates with poor prognosis in breast and prostate cancers, and the receptor promotes various traits of cancer aggressiveness including cell invasion. Here we show that miR-135a down-regulates the expression of ERRα through specific sequences of its 3'UTR. As a consequence miR-135a also reduces the expression of downstream targets of ERRα. miR-135a also decreases cell invasive potential in an ERRα-dependent manner. Our results suggest that the decreased expression of miR-135a in metastatic tumors leads to elevated ERRα expression, resulting in increased cell invasion capacities.

Induction of tissue kallikrein in human carotid atheroma does not lead to kallikrein-kinins pathway activation.

OBJECTIVE: The impairment of the tissue kallikrein-kinin system (KKS) may result in atheroma development. To determine the involvement of KKS in pathophysiology of human atherosclerosis, we examined the expression of all components of this system as well as angiotensinogen (another tissue kallikrein (TK) substrate), at messenger ribonucleic acid (mRNA) and protein levels in the human carotid artery with and without atheroma. METHODS: mRNA levels were compared with semi-quantitative reverse transcriptase-polymerase chain reaction (RT-PCR) between atheroma plaque and intact tissue obtained during carotid endarterectomy in 15 patients. The cellular localization of the transcripts and proteins was analyzed with in situ hybridization and immunohistochemistry. TK activity was measured using chromogenic substrate. RESULTS: The kininogen mRNA was not detected in carotid wall. The TK mRNA was increased four-fold and TK activity 23-fold in atheroma plaque compared with intact tissue. No difference was observed for B1, B2 receptors, kallistatin, angiotensinogen and protein-kinase G type 1alpha (PK-G) mRNAs. The TK and angiotensinogen transcripts as well as kininogen and angiotensinogen proteins were present in both intimal and medial cells. The kininogen immunoreactivity was weaker in atheroma. CONCLUSIONS: All KKS components were synthesized in arterial wall except kininogen probably coming from plasma. The absence of PK-G mRNA down-regulation in atheroma suggests that the kallikrein induction does not lead to KKS activation.