Dissecting the retinoid-induced differentiation of F9 embryonal stem cells by integrative genomics

Retinoids: Mechanisms of Action in Neuronal Cell Fate Acquisition

Neuronal differentiation has been shown to be directed by retinoid action during embryo development and has been exploited in various in vitro cell differentiation systems. In this review, we summarize the role of retinoids through the activation of their specific retinoic acid nuclear receptors during embryo development and also in a variety of in vitro strategies for neuronal differentiation, including recent efforts in driving cell specialization towards a range of neuronal subtypes and glial cells. Finally, we highlight the role of retinoic acid in recent protocols recapitulating nervous tissue complexity (cerebral organoids). Overall, we expect that this effort might pave the way for exploring the usage of specific synthetic retinoids for directing complex nervous tissue differentiation.

Deregulation of All-Trans Retinoic Acid Signaling and Development in Cancer

Cancer stem cells are the root cause of cancer, which, in essence, is a developmental disorder. All-trans retinoic acid (ATRA) signaling via ligand-activation of the retinoic acid receptors (RARs) plays a crucial role in tissue patterning and development during mammalian embryogenesis. In adults, active RARγ maintains the pool of hematopoietic stem cells, whereas active RARα drives myeloid cell differentiation. Various findings have revealed that ATRA signaling is deregulated in many cancers. The enzymes for ATRA synthesis are downregulated in colorectal, gastric, lung, and oropharyngeal cancers. ATRA levels within breast, ovarian, pancreatic, prostate, and renal cancer cells were lower than within their normal counterpart cells. The importance is that 0.24 nM ATRA activates RARγ (for stem cell stemness), whereas 100 times more is required to activate RARα (for differentiation). Moreover, RARγ is an oncogene regarding overexpression within colorectal, cholangiocarcinoma, hepatocellular, ovarian, pancreatic, and renal cancer cells. The microRNA (miR) 30a-5p downregulates expression of RARγ, and miR-30a/miR-30a-5p is a tumor suppressor for breast, colorectal, gastric, hepatocellular, lung, oropharyngeal, ovarian, pancreatic, prostate, and renal cancer. These complementary findings support the view that perturbations to ATRA signaling play a role in driving the abnormal behavior of cancer stem cells. Targeting ATRA synthesis and RARγ has provided promising approaches to eliminating cancer stem cells because such agents have been shown to drive cell death.

Targeting the Retinoic Acid Pathway to Eradicate Cancer Stem Cells

All-trans retinoic acid is a morphogen during embryogenesis and a teratogen. Cancer is an error of development, and the retinoic acid receptors (RAR) for all-trans retinoic acid play a role in cancer. Expression of the cytosolic aldehyde dehydrogenases, which mediate the last step to the synthesis of all-trans retinoic acid, is deregulated in various human cancers. Inhibiting these enzymes using a variety of agents reduced the proliferation of lung cancer cells, reduced the proliferation and induced apoptosis of ovarian, prostate, squamous, and uterine cancer cells, and sensitised breast, colorectal and ovarian cancer cells to chemotherapeutic agents. RARγ is an oncogene within some cases of AML, cholangiocarcinoma, colorectal cancer, clear cell renal cell carcinoma, hepatocellular carcinoma, pancreatic ductal adenocarcinoma, prostate cancer, and ovarian cancer. Pan-RAR and RARγ antagonist inhibition of the action of RARγ led to necroptosis of human prostate and pediatric brain tumour cancer stem cells. Treatment of hepatocellular carcinoma cells with the flavenoid acacetin, which interferes with the action of RARγ, decreased cell growth and induced apoptosis. Targeting the retinoic acid pathway is promising regarding the development of new drugs to eradicate cancer stem cells.

Antagonizing RARγ Drives Necroptosis of Cancer Stem Cells

There is a need for agents that eliminate cancer stem cells, which sustain cancer and are also largely responsible for disease relapse and metastasis. Conventional chemotherapeutics and radiotherapy are often highly effective against the bulk of cancer cells, which are proliferating, but spare cancer stem cells. Therapeutics that target cancer stem cells may also provide a bona fide cure for cancer. There are two rationales for targeting the retinoic acid receptor (RAR)γ. First, RARγ is expressed selectively within primitive cells. Second, RARγ is a putative oncogene for a number of human cancers, including cases of acute myeloid leukemia, cholangiocarcinoma, and colorectal, renal and hepatocellular carcinomas. Prostate cancer cells depend on active RARγ for their survival. Antagonizing all RARs caused necroptosis of prostate and breast cancer stem cell-like cells, and the cancer stem cells that gave rise to neurospheres from pediatric patients’ primitive neuroectodermal tumors and an astrocytoma. As tested for prostate cancer, antagonizing RARγ was sufficient to drive necroptosis. Achieving cancer-selectively is a longstanding paradigm for developing new treatments. The normal prostate epithelium was less sensitive to the RARγ antagonist and pan-RAR antagonist than prostate cancer cells, and fibroblasts and blood mononuclear cells were insensitive. The RARγ antagonist and pan-RAR antagonist are promising new cancer therapeutics.

Ligands and DNA in the allosteric control of retinoid receptors function

Retinoids are a family of compounds that include both vitamin A (all-trans retinol) and its naturally occurring metabolites such as retinoic acids (e.g. all-trans retinoic acid) as well as synthetic analogs. They are critically involved in the regulation of a wide variety of essential biological processes, such as embryogenesis and organogenesis, apoptosis, reproduction, vision, and the growth and differentiation of normal and neoplastic cells in vertebrates. The ability of these small molecules to control the expression of several hundred genes through binding to nuclear ligand-dependent transcription factors accounts for most of their functions. Three retinoic acid receptor (RARα,β,γ) and three retinoid X receptor (RXRα,β,γ) subtypes form a variety of RXR-RAR heterodimers that have been shown to mediate the pleiotropic effects of retinoids through the recruitment of high-molecular weight co-regulatory complexes to response-element DNA sequences found in the promoter region of their target genes. Hence, heterodimeric retinoid receptors are multidomain entities that respond to various incoming signals, such as ligand and DNA binding, by allosteric structural alterations which are the basis of further signal propagation. Here, we provide an overview of the current state of knowledge with regard to the structural mechanisms by which retinoids and DNA response elements act as allosteric effectors that may combine to finely tune RXR-RAR heterodimers activity.

The RARγ Oncogene: An Achilles Heel for Some Cancers

Cancer "stem cells" (CSCs) sustain the hierarchies of dividing cells that characterize cancer. The main causes of cancer-related mortality are metastatic disease and relapse, both of which originate primarily from CSCs, so their eradication may provide a bona fide curative strategy, though there maybe also the need to kill the bulk cancer cells. While classic anti-cancer chemotherapy is effective against the dividing progeny of CSCs, non-dividing or quiescent CSCs are often spared. Improved anti-cancer therapies therefore require approaches that target non-dividing CSCs, which must be underpinned by a better understanding of factors that permit these cells to maintain a stem cell-like state. During hematopoiesis, retinoic acid receptor (RAR) γ is selectively expressed by stem cells and their immediate progeny. It is overexpressed in, and is an oncogene for, many cancers including colorectal, renal and hepatocellular carcinoma, cholangiocarcinomas and some cases of acute myeloid leukemia that harbor RARγ fusion proteins. In vitro studies suggest that RARγ-selective and pan-RAR antagonists provoke the death of CSCs by necroptosis and point to antagonism of RARγ as a potential strategy to treat metastatic disease and relapse, and perhaps provide a cure for some cancers.

FOXC1 Downregulates Nanog Expression by Recruiting HDAC2 to Its Promoter in F9 Cells Treated by Retinoic Acid

FOXC1, a transcription factor involved in cell differentiation and embryogenesis, is demonstrated to be a negative regulator of Nanog in this study. FOXC1 is up-regulated in retinoic acid-induced differentiation of F9 Embryonal Carcinoma (EC) cells; furthermore, FOXC1 specifically inhibits the core pluripotency factor Nanog by binding to the proximal promoter. Overexpression of FOXC1 in F9 or knockdown in 3T3 results in the down-regulation or up-regulation of Nanog mRNA and proteins, respectively. In order to explain the mechanism by which FOXC1 inhibits Nanog expression, we identified the co-repressor HDAC2 from the FOXC1 interactome. FOXC1 recruits HDAC2 to Nanog promoter to decrease H3K27ac enrichment, resulting in transcription inhibition of Nanog. To the best of our knowledge, this is the first report that FOXC1 is involved in the epigenetic regulation of gene expression.

Structural basis for DNA recognition and allosteric control of the retinoic acid receptors RAR-RXR

Retinoic acid receptors (RARs) as a functional heterodimer with retinoid X receptors (RXRs), bind a diverse series of RA-response elements (RAREs) in regulated genes. Among them, the non-canonical DR0 elements are bound by RXR–RAR with comparable affinities to DR5 elements but DR0 elements do not act transcriptionally as independent RAREs. In this work, we present structural insights for the recognition of DR5 and DR0 elements by RXR–RAR heterodimer using x-ray crystallography, small angle x-ray scattering, and hydrogen/deuterium exchange coupled to mass spectrometry. We solved the crystal structures of RXR–RAR DNA-binding domain in complex with the Rarb2 DR5 and RXR–RXR DNA-binding domain in complex with Hoxb13 DR0. While cooperative binding was observed on DR5, the two molecules bound non-cooperatively on DR0 on opposite sides of the DNA. In addition, our data unveil the structural organization and dynamics of the multi-domain RXR–RAR DNA complexes providing evidence for DNA-dependent allosteric communication between domains. Differential binding modes between DR0 and DR5 were observed leading to differences in conformation and structural dynamics of the multi-domain RXR–RAR DNA complexes. These results reveal that the topological organization of the RAR binding element confer regulatory information by modulating the overall topology and structural dynamics of the RXR–RAR heterodimers.

DNA recognition by retinoic acid nuclear receptors

Retinoic acid receptors (RARs) heterodimerize with retinoid X receptors (RXRs) to regulate gene expression. The heterodimer recognizes the genome via a large and diverse repertoire of DNA response elements. Assessing the binding mode of RAR and RXR with various DNA response elements is important for understanding how they select their binding site and how DNA sequence and topology allosterically regulate RAR function. A number of complementary assays are often employed for analysis of the binding mode. To biochemically and structurally characterize RAR and RXR-DNA complexes, we describe how to express and purify RAR and RXR-DNA binding domains (DBDs) and multidomain constructs. We also describe the use of electrospray ionization mass spectrometry (ESI MS) and isothermal titration calorimetry (ITC) that give information about stoichiometry and binding affinity, as well as our approaches for co-crystallization of RAR and RXR DBDs with DNA.

Retinoic Acid-Regulated Target Genes During Development: Integrative Genomics Analysis

Retinoic acid (RA), a major natural active metabolite of vitamin A (VA) is well known to play critical roles in embryonic development. The effects of RA are mediated by nuclear receptors (RARs), which regulate the expression of gene batteries involved in cell growth and differentiation. Since the early 1990s several laboratories have focused on understanding how RA-regulated genes and RAR binding sites operate by studying the differentiation of embryonal carcinoma cells and embryonic stem cells. The development of hybridization-based microarray technology and high performance software analysis programs has allowed the characterization of thousands of RA-regulated genes. During the two last decades, publication of the genome sequence of various organisms has allowed advances in massive parallel sequencing and bioinformatics analysis of genome-wide data sets. These new generation sequencing (NGS) technologies have revolutionized the field by providing a global integrated picture of RA-regulated gene networks and the regulatory programs involved in cell fate decisions during embryonal carcinoma and embryonic stem cells differentiation. Now the challenge is to reconstruct the RA-regulated gene networks at the single cell level during the development of specialized embryonic tissues.

Regulation of RXR-RAR Heterodimers by RXR- and RAR-Specific Ligands and Their Combinations

The three subtypes (α, β, and γ) of the retinoic acid receptor (RAR) are ligand-dependent transcription factors that mediate retinoic acid signaling by forming heterodimers with the retinoid X receptor (RXR). Heterodimers are functional units that bind ligands (retinoids), transcriptional co-regulators and DNA, to regulate gene networks controlling cell growth, differentiation, and death. Using biochemical, crystallographic, and cellular approaches, we have set out to explore the spectrum of possibilities to regulate RXR-RAR heterodimer-dependent transcription through various pharmacological classes of RAR- and RXR- specific ligands, alone or in combination. We reveal the molecular details by which these compounds direct specificity and functionality of RXR-RAR heterodimers. Among these ligands, we have reevaluated and improved the molecular and structural definition of compounds CD2665, Ro41-5253, LE135, or LG100754, highlighting novel functional features of these molecules. Our analysis reveals a model of RXR-RAR heterodimer action in which each subunit retains its intrinsic properties in terms of ligand and co-regulator binding. However, their interplay upon the combined action of RAR- and RXR-ligands allows for the fine tuning of heterodimer activity. It also stresses the importance of accurate ligand characterization to use synthetic selective retinoids appropriately and avoid data misinterpretations.

Inferring Regulatory Programs Governing Region Specificity of Neuroepithelial Stem Cells during Early Hindbrain and Spinal Cord Development

Neuroepithelial stem cells (NSC) from different anatomical regions of the embryonic neural tube's rostrocaudal axis can differentiate into diverse central nervous system tissues, but the transcriptional regulatory networks governing these processes are incompletely understood. Here, we measure region-specific NSC gene expression along the rostrocaudal axis in a human pluripotent stem cell model of early central nervous system development over a 72-h time course, spanning the hindbrain to cervical spinal cord. We introduce Escarole, a probabilistic clustering algorithm for non-stationary time series, and combine it with prior-based regulatory network inference to identify genes that are regulated dynamically and predict their upstream regulators. We identify known regulators of patterning and neural development, including the HOX genes, and predict a direct regulatory connection between the transcription factor POU3F2 and target gene STMN2. We demonstrate that POU3F2 is required for expression of STMN2, suggesting that this regulatory connection is important for region specificity of NSCs.

Znf703 is a novel RA target in the neural plate border

Znf703 is an RAR- and Wnt-inducible transcription factor that exhibits a complex expression pattern in the developing embryo: Znf703 mRNA is found in the early circumblastoporal ring, then later throughout the neural plate and its border, and subsequently in the mid/hindbrain and somites. We show that Znf703 has a different and separable function in early mesoderm versus neural crest and placode development. Independent of its early knockdown phenotype on Gdf3 and Wnt8, Znf703 disrupts patterning of distinct neural crest migratory streams normally delineated by Sox10, Twist, and Foxd3 and inhibits otocyst formation and otic expression of Sox10 and Eya1. Furthermore, Znf703 promotes massive overgrowth of SOX2+ cells, disrupting the SoxB1 balance at the neural plate border. Despite prominent expression in other neural plate border-derived cranial and sensory domains, Znf703 is selectively absent from the otocyst, suggesting that Znf703 must be specifically cleared or down-regulated for proper otic development. We show that mutation of the putative Groucho-repression domain does not ameliorate Znf703 effects on mesoderm, neural crest, and placodes. We instead provide evidence that Znf703 requires the Buttonhead domain for transcriptional repression.

Ethanol promotes differentiation of embryonic stem cells through retinoic acid receptor-γ

Ethanol (EtOH) is a teratogen, but its teratogenic mechanisms are not fully understood. The alcohol form of vitamin A (retinol/ROL) can be oxidized to all-trans-retinoic acid (RA), which plays a critical role in stem cell differentiation and development. Using an embryonic stem cell (ESC) model to analyze EtOH's effects on differentiation, we show here that EtOH and acetaldehyde, but not acetate, increase differentiation-associated mRNA levels, and that EtOH decreases pluripotency-related mRNAs. Using reporter assays, ChIP assays, and retinoic acid receptor-γ (RARγ) knockout ESC lines generated by CRISPR/Cas9 and homologous recombination, we demonstrate that EtOH signals via RARγ binding to RA response elements (RAREs) in differentiation-associated gene promoters or enhancers. We also report that EtOH-mediated increases in homeobox A1 (Hoxa1) and cytochrome P450 family 26 subfamily A member 1 (Cyp26a1) transcripts, direct RA target genes, require the expression of the RA-synthesizing enzyme, aldehyde dehydrogenase 1 family member A2 (Aldh1a2), suggesting that EtOH-mediated induction of Hoxa1 and Cyp26a1 requires ROL from the serum. As shown with CRISPR/Cas9 knockout lines, the retinol dehydrogenase gene Rdh10 and a functional RARE in the ROL transporter stimulated by retinoic acid 6 (Stra6) gene are required for EtOH induction of Hoxa1 and Cyp26a1 We conclude that EtOH stimulates stem cell differentiation by increasing the influx and metabolism of ROL for downstream RARγ-dependent transcription. In stem cells, EtOH may shift cell fate decisions to alter developmental outcomes by increasing endogenous ROL/RA signaling via increased Stra6 expression and ROL oxidation.

Metabolic profile and differentiation potential of extraembryonic endoderm-like cells

Glucose metabolism has a crucial role for providing substrates required to generate ATP and regulate the epigenetic landscape. We reported that F9 embryonal carcinoma stem-like cells require cytosolic reactive oxygen species to differentiate into extraembryonic endoderm; however, mitochondrial sources were not examined. To extend these studies, we examined the metabolic profile of early and late-passage F9 cells, and show that their ability to differentiate is similar, even though each population has dramatically different metabolic profiles. Differentiated early-passage cells relied on glycolysis, while differentiated late-passage cells transitioned towards oxidative phosphorylation (OXPHOS). Unexpectedly, electron transport chain protein stoichiometry was disrupted in differentiated late-passage cells, whereas genes encoding mitofusion 1 and 2, which promote mitochondrial fusion and favor OXPHOS, were upregulated in differentiated early-passage cells. Despite this, early-passage cells cultured under conditions to promote glycolysis showed enhanced differentiation, whereas promoting OXPHOS in late-passage cells showed a similar trend. Further analysis revealed that the distinct metabolic profiles seen between the two populations is largely associated with changes in genomic integrity, linking metabolism to passage number. Together, these results indicate that passaging has no effect on the potential for F9 cells to differentiate into extraembryonic endoderm; however, it does impact their metabolic profile. Thus, it is imperative to determine the molecular and metabolic status of a stem cell population before considering its utility as a therapeutic tool for regenerative medicine.

RXR heterodimers orchestrate transcriptional control of neurogenesis and cell fate specification

Retinoid X Receptors (RXRs) are unique and enigmatic members of the nuclear receptor (NR) family with extensive and complex biological functions in cellular differentiation. On the one hand, RXRs through permissive heterodimerization with other NRs are able to integrate multiple lipid signaling pathways and are believed to play a central role to coordinate the development of the central nervous system. On the other hand, RXRs may have heterodimer-independent functions as well. Therefore, a more RXR-centric analysis is warranted to identify its genomic binding sites and regulated gene networks, which are orchestrating the earliest events in neuronal differentiation. Recently developed genome-wide approaches allow systematic analyses of the RXR-driven neural differentiation. Here we applied next generation sequencing-based methodology to track the dynamic redistribution of the RXR cistrome along the path of embryonic stem cell to glutamatergic neuron differentiation. We identified Retinoic Acid Receptor (RAR) and Liver X Receptor (LXR) as dominant heterodimeric partners of RXR in these cellular stages. Our data presented here characterize the RAR:RXR and LXR:RXR-mediated transcriptional program in embryonic stem cells, neural progenitors and terminally differentiated neurons. Considering the growing evidence for dysregulated RXR-mediated signaling in neurodegenerative disorders, such as Alzheimer's Disease or Amyotrophic Lateral Sclerosis, the data presented here will be also a valuable resource for the field of neuro(patho)biology.

Combinatorial knockout of RARα, RARβ, and RARγ completely abrogates transcriptional responses to retinoic acid in murine embryonic stem cells

All-trans-retinoic acid (RA), a potent inducer of cellular differentiation, functions as a ligand for retinoic acid receptors (RARα, β, and γ). RARs are activated by ligand binding, which induces transcription of direct genomic targets. However, whether embryonic stem cells respond to RA through routes that do not involve RARs is unknown. Here, we used CRISPR technology to introduce biallelic frameshift mutations in RARα, RARβ, and RARγ, thereby abrogating all RAR functions in murine embryonic stem cells. We then evaluated RA-responsiveness of the RAR-null cells using RNA-Seq transcriptome analysis. We found that the RAR-null cells display no changes in transcripts in response to RA, demonstrating that the RARs are essential for the regulation of all transcripts in murine embryonic stem cells in response to RA. Our key finding, that in embryonic stem cells the transcriptional effects of RA all depend on RARs, addresses a long-standing topic of discussion in the field of retinoic acid signaling.

Protagonist or antagonist? The complex roles of retinoids in the regulation of hematopoietic stem cells and their specification from pluripotent stem cells

Hematopoietic stem cells (HSCs) are multipotent cells responsible for the maintenance of the hematopoietic system throughout life. Dysregulation of the balance in HSC self-renewal, death, and differentiation can have serious consequences such as myelodysplastic syndromes or leukemia. All-trans retinoic acid (ATRA), the biologically active metabolite of vitamin A/RA, has been shown to have pleiotropic effects on hematopoietic cells, enhancing HSC self-renewal while also increasing differentiation of more mature progenitors. Furthermore, ATRA has been shown to have key roles in regulating the specification and formation of hematopoietic cells from pluripotent stem cells including embryonic stem cells (ESCs) and induced pluripotent stem cells (iPSCs). Here, we summarize the known roles of vitamin A and RA receptors in the regulation of hematopoiesis from HSCs, ES, and iPSCs.

The ancestral retinoic acid receptor was a low-affinity sensor triggering neuronal differentiation

Retinoic acid (RA) is an important intercellular signaling molecule in vertebrate development, with a well-established role in the regulation of hox genes during hindbrain patterning and in neurogenesis. However, the evolutionary origin of the RA signaling pathway remains elusive. To elucidate the evolution of the RA signaling system, we characterized RA metabolism and signaling in the marine annelid Platynereis dumerilii, a powerful model for evolution, development, and neurobiology. Binding assays and crystal structure analyses show that the annelid retinoic acid receptor (RAR) binds RA and activates transcription just as vertebrate RARs, yet with a different ligand-binding pocket and lower binding affinity, suggesting a permissive rather than instructive role of RA signaling. RAR knockdown and RA treatment of swimming annelid larvae further reveal that the RA signal is locally received in the medial neuroectoderm, where it controls neurogenesis and axon outgrowth, whereas the spatial colinear hox gene expression in the neuroectoderm remains unaffected. These findings suggest that one early role of the new RAR in bilaterian evolution was to control the spatially restricted onset of motor and interneuron differentiation in the developing ventral nerve cord and to indicate that the regulation of hox-controlled anterior-posterior patterning arose only at the base of the chordates, concomitant with a high-affinity RAR needed for the interpretation of a complex RA gradient.

Inhibition of Starvation-Triggered Endoplasmic Reticulum Stress, Autophagy, and Apoptosis in ARPE-19 Cells by Taurine through Modulating the Expression of Calpain-1 and Calpain-2

Age-related macular degeneration (AMD) is a complex disease with multiple initiators and pathways that converge on death for retinal pigment epithelial (RPE) cells. In this study, effects of taurine on calpains, autophagy, endoplasmic reticulum (ER) stress, and apoptosis in ARPE-19 cells (a human RPE cell line) were investigated. We first confirmed that autophagy, ER stress and apoptosis in ARPE-19 cells were induced by Earle’s balanced salt solution (EBSS) through starvation to induce RPE metabolic stress. Secondly, inhibition of ER stress by 4-phenyl butyric acid (4-PBA) alleviated autophagy and apoptosis, and suppression of autophagy by 3-methyl adenine (3-MA) reduced the cell apoptosis, but the ER stress was minimally affected. Thirdly, the apoptosis, ER stress and autophagy were inhibited by gene silencing of calpain-2 and overexpression of calpain-1, respectively. Finally, taurine suppressed both the changes of the important upstream regulators (calpain-1 and calpain-2) and the activation of ER stress, autophagy and apoptosis, and taurine had protective effects on the survival of ARPE-19 cells. Collectively, this data indicate that taurine inhibits starvation-triggered endoplasmic reticulum stress, autophagy, and apoptosis in ARPE-19 cells by modulating the expression of calpain-1 and calpain-2.

Epimetheus - a multi-profile normalizer for epigenomic sequencing data

Background

Exponentially increasing numbers of NGS-based epigenomic datasets in public repositories like GEO constitute an enormous source of information that is invaluable for integrative and comparative studies of gene regulatory mechanisms. One of today’s challenges for such studies is to identify functionally informative local and global patterns of chromatin states in order to describe the regulatory impact of the epigenome in normal cell physiology and in case of pathological aberrations. Critically, the most preferred Chromatin ImmunoPrecipitation-Sequencing (ChIP-Seq) is inherently prone to significant variability between assays, which poses significant challenge on comparative studies. One challenge concerns data normalization to adjust sequencing depth variation.

Results

Currently existing tools either apply linear scaling corrections and/or are restricted to specific genomic regions, which can be prone to biases. To overcome these restrictions without any external biases, we developed Epimetheus, a genome-wide quantile-based multi-profile normalization tool for histone modification data and related datasets.

Conclusions

Epimetheus has been successfully used to normalize epigenomics data in previous studies on X inactivation in breast cancer and in integrative studies of neuronal cell fate acquisition and tumorigenic transformation; Epimetheus is freely available to the scientific community.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1655-3) contains supplementary material, which is available to authorized users.

Mechanisms Regulating Stemness and Differentiation in Embryonal Carcinoma Cells

Just over ten years have passed since the seminal Takahashi-Yamanaka paper, and while most attention nowadays is on induced, embryonic, and cancer stem cells, much of the pioneering work arose from studies with embryonal carcinoma cells (ECCs) derived from teratocarcinomas. This original work was broad in scope, but eventually led the way for us to focus on the components involved in the gene regulation of stemness and differentiation. As the name implies, ECCs are malignant in nature, yet maintain the ability to differentiate into the 3 germ layers and extraembryonic tissues, as well as behave normally when reintroduced into a healthy blastocyst. Retinoic acid signaling has been thoroughly interrogated in ECCs, especially in the F9 and P19 murine cell models, and while we have touched on this aspect, this review purposely highlights how some key transcription factors regulate pluripotency and cell stemness prior to this signaling. Another major focus is on the epigenetic regulation of ECCs and stem cells, and, towards that end, this review closes on what we see as a new frontier in combating aging and human disease, namely, how cellular metabolism shapes the epigenetic landscape and hence the pluripotency of all stem cells.

Allosteric Regulation in the Ligand Binding Domain of Retinoic Acid Receptorγ

Retinoic acid (RA) plays key roles in cell differentiation and growth arrest through nuclear retinoic acid receptors (RARs), which are ligand-dependent transcription factors. While the main trigger of RAR activation is the binding of RA, phosphorylation of the receptors has also emerged as an important regulatory signal. Phosphorylation of the RARγ N-terminal domain (NTD) is known to play a functional role in neuronal differentiation. In this work, we investigated the phosphorylation of RARγ ligand binding domain (LBD), and present evidence that the phosphorylation status of the LBD affects the phosphorylation of the NTD region. We solved the X-ray structure of a phospho-mimetic mutant of the LBD (RARγ S371E), which we used in molecular dynamics simulations to characterize the consequences of the S371E mutation on the RARγ structural dynamics. Combined with simulations of the wild-type LBD, we show that the conformational equilibria of LBD salt bridges (notably R387-D340) are affected by the S371E mutation, which likely affects the recruitment of the kinase complex that phosphorylates the NTD. The molecular dynamics simulations also showed that a conservative mutation in this salt bridge (R387K) affects the dynamics of the LBD without inducing large conformational changes. Finally, cellular assays showed that the phosphorylation of the NTD of RARγ is differentially regulated by retinoic acid in RARγWT and in the S371N, S371E and R387K mutants. This multidisciplinary work highlights an allosteric coupling between phosphorylations of the LBD and the NTD of RARγ and supports the importance of structural dynamics involving electrostatic interactions in the regulation of RARs activity.

Mechanisms of retinoic acid signaling during cardiogenesis

Substantial experimental and epidemiological data have highlighted the interplay between nutritional and genetic factors in the development of congenital heart defects. Retinoic acid (RA), a derivative of vitamin A, plays a key role during vertebrate development including the formation of the heart. Retinoids bind to RA and retinoid X receptors (RARs and RXRs) which then regulate tissue-specific genes. Here, we will focus on the roles of RA signaling and receptors in gene regulation during cardiogenesis, and the consequence of deregulated retinoid signaling on heart formation and congenital heart defects.

Reconstructed cell fate-regulatory programs in stem cells reveal hierarchies and key factors of neurogenesis

Cell lineages, which shape the body architecture and specify cell functions, derive from the integration of a plethora of cell intrinsic and extrinsic signals. These signals trigger a multiplicity of decisions at several levels to modulate the activity of dynamic gene regulatory networks (GRNs), which ensure both general and cell-specific functions within a given lineage, thereby establishing cell fates. Significant knowledge about these events and the involved key drivers comes from homogeneous cell differentiation models. Even a single chemical trigger, such as the morphogen all-trans retinoic acid (RA), can induce the complex network of gene-regulatory decisions that matures a stem/precursor cell to a particular step within a given lineage. Here we have dissected the GRNs involved in the RA-induced neuronal or endodermal cell fate specification by integrating dynamic RXRA binding, chromatin accessibility, epigenetic promoter epigenetic status, and the transcriptional activity inferred from RNA polymerase II mapping and transcription profiling. Our data reveal how RA induces a network of transcription factors (TFs), which direct the temporal organization of cognate GRNs, thereby driving neuronal/endodermal cell fate specification. Modeling signal transduction propagation using the reconstructed GRNs indicated critical TFs for neuronal cell fate specification, which were confirmed by CRISPR/Cas9-mediated genome editing. Overall, this study demonstrates that a systems view of cell fate specification combined with computational signal transduction models provides the necessary insight in cellular plasticity for cell fate engineering. The present integrated approach can be used to monitor the in vitro capacity of (engineered) cells/tissues to establish cell lineages for regenerative medicine.

Genome-wide Analysis of RARβ Transcriptional Targets in Mouse Striatum Links Retinoic Acid Signaling with Huntington's Disease and Other Neurodegenerative Disorders

Retinoic acid (RA) signaling through retinoic acid receptors (RARs), known for its multiple developmental functions, emerged more recently as an important regulator of adult brain physiology. How RAR-mediated regulation is achieved is poorly known, partly due to the paucity of information on critical target genes in the brain. Also, it is not clear how reduced RA signaling may contribute to pathophysiology of diverse neuropsychiatric disorders. We report the first genome-wide analysis of RAR transcriptional targets in the brain. Using chromatin immunoprecipitation followed by high-throughput sequencing and transcriptomic analysis of RARβ-null mutant mice, we identified genomic targets of RARβ in the striatum. Characterization of RARβ transcriptional targets in the mouse striatum points to mechanisms through which RAR may control brain functions and display neuroprotective activity. Namely, our data indicate with statistical significance (FDR 0.1) a strong contribution of RARβ in controlling neurotransmission, energy metabolism, and transcription, with a particular involvement of G-protein coupled receptor (p = 5.0e^-5), cAMP (p = 4.5e^-4), and calcium signaling (p = 3.4e^-3). Many identified RARβ target genes related to these pathways have been implicated in Alzheimer's, Parkinson's, and Huntington's disease (HD), raising the possibility that compromised RA signaling in the striatum may be a mechanistic link explaining the similar affective and cognitive symptoms in these diseases. The RARβ transcriptional targets were particularly enriched for transcripts affected in HD. Using the R6/2 transgenic mouse model of HD, we show that partial sequestration of RARβ in huntingtin protein aggregates may account for reduced RA signaling reported in HD.

Phosphoproteome and Transcriptome of RA-Responsive and RA-Resistant Breast Cancer Cell Lines

Retinoic acid (RA), the main active vitamin A metabolite, controls multiple biological processes such as cell proliferation and differentiation through genomic programs and kinase cascades activation. Due to these properties, RA has proven anti-cancer capacity. Several breast cancer cells respond to the antiproliferative effects of RA, while others are RA-resistant. However, the overall signaling and transcriptional pathways that are altered in such cells have not been elucidated. Here, in a large-scale analysis of the phosphoproteins and in a genome-wide analysis of the RA-regulated genes, we compared two human breast cancer cell lines, a RA-responsive one, the MCF7 cell line, and a RA-resistant one, the BT474 cell line, which depicts several alterations of the “kinome”. Using high-resolution nano-LC-LTQ-Orbitrap mass spectrometry associated to phosphopeptide enrichment, we found that several proteins involved in signaling and in transcription, are differentially phosphorylated before and after RA addition. The paradigm of these proteins is the RA receptor α (RARα), which was phosphorylated in MCF7 cells but not in BT474 cells after RA addition. The panel of the RA-regulated genes was also different. Overall our results indicate that RA resistance might correlate with the deregulation of the phosphoproteome with consequences on gene expression.

cDREM: inferring dynamic combinatorial gene regulation

Genes are often combinatorially regulated by multiple transcription factors (TFs). Such combinatorial regulation plays an important role in development and facilitates the ability of cells to respond to different stresses. While a number of approaches have utilized sequence and ChIP-based datasets to study combinational regulation, these have often ignored the combinational logic and the dynamics associated with such regulation. Here we present cDREM, a new method for reconstructing dynamic models of combinatorial regulation. cDREM integrates time series gene expression data with (static) protein interaction data. The method is based on a hidden Markov model and utilizes the sparse group Lasso to identify small subsets of combinatorially active TFs, their time of activation, and the logical function they implement. We tested cDREM on yeast and human data sets. Using yeast we show that the predicted combinatorial sets agree with other high throughput genomic datasets and improve upon prior methods developed to infer combinatorial regulation. Applying cDREM to study human response to flu, we were able to identify several combinatorial TF sets, some of which were known to regulate immune response while others represent novel combinations of important TFs.

Pan-cancer analyses of the nuclear receptor superfamily

Nuclear receptors (NR) act as an integrated conduit for environmental and hormonal signals to govern genomic responses, which relate to cell fate decisions. We review how their integrated actions with each other, shared co-factors and other transcription factors are disrupted in cancer. Steroid hormone nuclear receptors are oncogenic drivers in breast and prostate cancer and blockade of signaling is a major therapeutic goal. By contrast to blockade of receptors, in other cancers enhanced receptor function is attractive, as illustrated initially with targeting of retinoic acid receptors in leukemia. In the post-genomic era large consortia, such as The Cancer Genome Atlas, have developed a remarkable volume of genomic data with which to examine multiple aspects of nuclear receptor status in a pan-cancer manner. Therefore to extend the review of NR function we have also undertaken bioinformatics analyses of NR expression in over 3000 tumors, spread across six different tumor types (bladder, breast, colon, head and neck, liver and prostate). Specifically, to ask how the NR expression was distorted (altered expression, mutation and CNV) we have applied bootstrapping approaches to simulate data for comparison, and also compared these NR findings to 12 other transcription factor families. Nuclear receptors were uniquely and uniformly downregulated across all six tumor types, more than predicted by chance. These approaches also revealed that each tumor type had a specific NR expression profile but these were most similar between breast and prostate cancer. Some NRs were down-regulated in at least five tumor types (e.g. NR3C2/MR and NR5A2/LRH-1)) whereas others were uniquely down-regulated in one tumor (e.g. NR1B3/RARG). The downregulation was not driven by copy number variation or mutation and epigenetic mechanisms maybe responsible for the altered nuclear receptor expression.

Integrative genomic analysis in K562 chronic myelogenous leukemia cells reveals that proximal NCOR1 binding positively regulates genes that govern erythroid differentiation and Imatinib sensitivity

To define the functions of NCOR1 we developed an integrative analysis that combined ENCODE and NCI-60 data, followed by in vitro validation. NCOR1 and H3K9me3 ChIP-Seq, FAIRE-seq and DNA CpG methylation interactions were related to gene expression using bootstrapping approaches. Most NCOR1 combinations (24/44) were associated with significantly elevated level expression of protein coding genes and only very few combinations related to gene repression. DAVID's biological process annotation revealed that elevated gene expression was uniquely associated with acetylation and ETS binding. A matrix of gene and drug interactions built on NCI-60 data identified that Imatinib significantly targeted the NCOR1 governed transcriptome. Stable knockdown of NCOR1 in K562 cells slowed growth and significantly repressed genes associated with NCOR1 cistrome, again, with the GO terms acetylation and ETS binding, and significantly dampened sensitivity to Imatinib-induced erythroid differentiation. Mining public microarray data revealed that NCOR1-targeted genes were significantly enriched in Imatinib response gene signatures in cell lines and chronic myelogenous leukemia (CML) patients. These approaches integrated cistrome, transcriptome and drug sensitivity relationships to reveal that NCOR1 function is surprisingly most associated with elevated gene expression, and that these targets, both in CML cell lines and patients, associate with sensitivity to Imatinib.

RAR/RXR binding dynamics distinguish pluripotency from differentiation associated cis-regulatory elements

In mouse embryonic cells, ligand-activated retinoic acid receptors (RARs) play a key role in inhibiting pluripotency-maintaining genes and activating some major actors of cell differentiation. To investigate the mechanism underlying this dual regulation, we performed joint RAR/RXR ChIP-seq and mRNA-seq time series during the first 48 h of the RA-induced Primitive Endoderm (PrE) differentiation process in F9 embryonal carcinoma (EC) cells. We show here that this dual regulation is associated with RAR/RXR genomic redistribution during the differentiation process. In-depth analysis of RAR/RXR binding sites occupancy dynamics and composition show that in undifferentiated cells, RAR/RXR interact with genomic regions characterized by binding of pluripotency-associated factors and high prevalence of the non-canonical DR0-containing RA response element. By contrast, in differentiated cells, RAR/RXR bound regions are enriched in functional Sox17 binding sites and are characterized with a higher frequency of the canonical DR5 motif. Our data offer an unprecedentedly detailed view on the action of RA in triggering pluripotent cell differentiation and demonstrate that RAR/RXR action is mediated via two different sets of regulatory regions tightly associated with cell differentiation status.

Dissecting the calcium-induced differentiation of human primary keratinocytes stem cells by integrative and structural network analyses

The molecular details underlying the time-dependent assembly of protein complexes in cellular networks, such as those that occur during differentiation, are largely unexplored. Focusing on the calcium-induced differentiation of primary human keratinocytes as a model system for a major cellular reorganization process, we look at the expression of genes whose products are involved in manually-annotated protein complexes. Clustering analyses revealed only moderate co-expression of functionally related proteins during differentiation. However, when we looked at protein complexes, we found that the majority (55%) are composed of non-dynamic and dynamic gene products ('di-chromatic'), 19% are non-dynamic, and 26% only dynamic. Considering three-dimensional protein structures to predict steric interactions, we found that proteins encoded by dynamic genes frequently interact with a common non-dynamic protein in a mutually exclusive fashion. This suggests that during differentiation, complex assemblies may also change through variation in the abundance of proteins that compete for binding to common proteins as found in some cases for paralogous proteins. Considering the example of the TNF-α/NFκB signaling complex, we suggest that the same core complex can guide signals into diverse context-specific outputs by addition of time specific expressed subunits, while keeping other cellular functions constant. Thus, our analysis provides evidence that complex assembly with stable core components and competition could contribute to cell differentiation.

Nuclear and extranuclear effects of vitamin A

Vitamin A or retinol is a multifunctional vitamin that is essential at all stages of life from embryogenesis to adulthood. Up to now, it has been accepted that the effects of vitamin A are exerted by active metabolites, the major ones being 11-cis retinal for vision, and all trans-retinoic acid (RA) for cell growth and differentiation. Basically RA binds nuclear receptors, RARs, which regulate the expression of a battery of target genes in a ligand dependent manner. During the last decade, new scenarios have been discovered, providing a rationale for the understanding of other long-noted but not explained functions of retinol. These novel scenarios involve: (i) other nuclear receptors such as PPAR β/δ, which regulate the expression of other target genes with other functions; (ii) extranuclear and nontranscriptional effects, such as the activation of kinases, which phosphorylate RARs and other transcription factors, thus expanding the list of the RA-activated genes; (iii) finally, vitamin A is active per se and can work as a cytokine that regulates gene transcription by activating STRA6. New effects of vitamin A and RA are continuously being discovered in new fields, revealing new targets and new mechanisms thus improving the understanding the pleiotropicity of their effects.

Vitamin D receptor and RXR in the post-genomic era

Following the elucidation of the human genome and components of the epigenome, it is timely to revisit what is known of vitamin D receptor (VDR) function. Early transcriptomic studies using microarray approaches focused on the protein coding mRNA that were regulated by the VDR, usually following treatment with ligand. These studies quickly established the approximate size and surprising diversity of the VDR transcriptome, revealing it to be highly heterogenous and cell type and time dependent. Investigators also considered VDR regulation of non-protein coding RNA and again, cell and time dependency was observed. Attempts to integrate mRNA and miRNA regulation patterns are beginning to reveal patterns of co-regulation and interaction that allow for greater control of mRNA expression, and the capacity to govern more complex cellular events. Alternative splicing in the trasncriptome has emerged as a critical process in transcriptional control and there is evidence of the VDR interacting with components of the splicesome. ChIP-Seq approaches have proved to be pivotal to reveal the diversity of the VDR binding choices across cell types and following treatment, and have revealed that the majority of these are non-canonical in nature. The underlying causes driving the diversity of VDR binding choices remain enigmatic. Finally, genetic variation has emerged as important to impact the transcription factor affinity towards genomic binding sites, and recently the impact of this on VDR function has begun to be considered.

An alternative retinoic acid-responsive Stra6 promoter regulated in response to retinol deficiency

Cellular uptake of vitamin A (retinol) is essential for many biological functions. The Stra6 protein binds the serum retinol-binding protein, RBP4, and acts in conjunction with the enzyme lecithin:retinol acyltransferase to facilitate retinol uptake in some cell types. We show that in embryonic stem (ES) cells and in some tissues, the Stra6 gene encodes two distinct mRNAs transcribed from two different promoters. Whereas both are all-trans-retinoic acid (RA)-responsive in ES cells, the downstream promoter contains a half-site RA response element (RARE) and drives an ∼ 13-fold, RA-associated increase in luciferase reporter activity. We employed CRISPR-Cas9 genome editing to show that the endogenous RARE is required for RA-induced transcription of both Stra6 isoforms. We further demonstrate that in ES cells, 1) both RARγ and RXRα are present at the Stra6 RARE; 2) RA increases co-activator p300 (KAT3B) binding and histone H3 Lys-27 acetylation at both promoters; 3) RA decreases Suz12 levels and histone H3 Lys-27 trimethylation epigenetic marks at both promoters; and 4) these epigenetic changes are diminished in the absence of RARγ. In the brains of WT mice, both the longer and the shorter Stra6 transcript (Stra6L and Stra6S, respectively) are highly expressed, whereas these transcripts are found only at low levels in RARγ(-/-) mice. In the brains of vitamin A-deficient mice, both Stra6L and Stra6S levels are decreased. In contrast, in the vitamin A-deficient kidneys, the Stra6L levels are greatly increased, whereas Stra6S levels are decreased. Our data show that kidneys respond to retinol deficiency by differential Stra6 promoter usage, which may play a role in the retention of retinol when vitamin A is low.

Structural basis of natural promoter recognition by the retinoid X nuclear receptor

Retinoid X receptors (RXRs) act as homodimers or heterodimerisation partners of class II nuclear receptors. RXR homo- and heterodimers bind direct repeats of the half-site (A/G)G(G/T)TCA separated by 1 nucleotide (DR1). We present a structural characterization of RXR-DNA binding domain (DBD) homodimers on several natural DR1s and an idealized symmetric DR1. Homodimers displayed asymmetric binding, with critical high-affinity interactions accounting for the 3' positioning of RXR in heterodimers on DR1s. Differing half-site and spacer DNA sequence induce changes in RXR-DBD homodimer conformation notably in the dimerization interface such that natural DR1s are bound with higher affinity than an idealized symmetric DR1. Subtle changes in the consensus DR1 DNA sequence therefore specify binding affinity through altered RXR-DBD-DNA contacts and changes in DBD conformation suggesting a general model whereby preferential half-site recognition determines polarity of heterodimer binding to response elements.

Retinoic acid receptors: from molecular mechanisms to cancer therapy

Retinoic acid (RA), the major bioactive metabolite of retinol or vitamin A, induces a spectrum of pleiotropic effects in cell growth and differentiation that are relevant for embryonic development and adult physiology. The RA activity is mediated primarily by members of the retinoic acid receptor (RAR) subfamily, namely RARα, RARβ and RARγ, which belong to the nuclear receptor (NR) superfamily of transcription factors. RARs form heterodimers with members of the retinoid X receptor (RXR) subfamily and act as ligand-regulated transcription factors through binding specific RA response elements (RAREs) located in target genes promoters. RARs also have non-genomic effects and activate kinase signaling pathways, which fine-tune the transcription of the RA target genes. The disruption of RA signaling pathways is thought to underlie the etiology of a number of hematological and non-hematological malignancies, including leukemias, skin cancer, head/neck cancer, lung cancer, breast cancer, ovarian cancer, prostate cancer, renal cell carcinoma, pancreatic cancer, liver cancer, glioblastoma and neuroblastoma. Of note, RA and its derivatives (retinoids) are employed as potential chemotherapeutic or chemopreventive agents because of their differentiation, anti-proliferative, pro-apoptotic, and anti-oxidant effects. In humans, retinoids reverse premalignant epithelial lesions, induce the differentiation of myeloid normal and leukemic cells, and prevent lung, liver, and breast cancer. Here, we provide an overview of the biochemical and molecular mechanisms that regulate the RA and retinoid signaling pathways. Moreover, mechanisms through which deregulation of RA signaling pathways ultimately impact on cancer are examined. Finally, the therapeutic effects of retinoids are reported.

High-wire act: the poised genome and cellular memory

Emerging evidence aided by genome-wide analysis of chromatin and transcriptional states has shed light on the mechanisms by which stem cells achieve cellular memory. The epigenetic and transcriptional plasticity governing stem cell behavior is highlighted by the identification of 'poised' genes, which permit cells to maintain readiness to undertake alternate developmental fates. This review focuses on two crucial mechanisms of gene poising: bivalent chromatin marks and RNA polymerase II stalling. We provide the context for these mechanisms by exploring the current consensus on the regulation of chromatin states, especially in quiescent adult stem cells, where poised genes are critical for recapitulating developmental choices, leading to regenerative function.

SMARTS: reconstructing disease response networks from multiple individuals using time series gene expression data

MOTIVATION

Current methods for reconstructing dynamic regulatory networks are focused on modeling a single response network using model organisms or cell lines. Unlike these models or cell lines, humans differ in their background expression profiles due to age, genetics and life factors. In addition, there are often differences in start and end times for time series human data and in the rate of progress based on the specific individual. Thus, new methods are required to integrate time series data from multiple individuals when modeling and constructing disease response networks.

RESULTS

We developed Scalable Models for the Analysis of Regulation from Time Series (SMARTS), a method integrating static and time series data from multiple individuals to reconstruct condition-specific response networks in an unsupervised way. Using probabilistic graphical models, SMARTS iterates between reconstructing different regulatory networks and assigning individuals to these networks, taking into account varying individual start times and response rates. These models can be used to group different sets of patients and to identify transcription factors that differentiate the observed responses between these groups. We applied SMARTS to analyze human response to influenza and mouse brain development. In both cases, it was able to greatly improve baseline groupings while identifying key relevant TFs that differ between the groups. Several of these groupings and TFs are known to regulate the relevant processes while others represent novel hypotheses regarding immune response and development.

AVAILABILITY AND IMPLEMENTATION

Software and Supplementary information are available at http://sb.cs.cmu.edu/smarts/.

CONTACT

zivbj@cs.cmu.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Evolution of retinoic acid receptors and retinoic acid signaling

Retinoic acid (RA) is a vitamin A-derived morphogen controlling important developmental processes in vertebrates, and more generally in chordates, including axial patterning and tissue formation and differentiation. In the embryo, endogenous RA levels are controlled by RA synthesizing and degrading enzymes and the RA signal is transduced by two retinoid receptors: the retinoic acid receptor (RAR) and the retinoid X receptor (RXR). Both RAR and RXR are members of the nuclear receptor superfamily of ligand-activated transcription factors and mainly act as heterodimers to activate the transcription of target genes in the presence of their ligand, all-trans RA. This signaling pathway was long thought to be a chordate innovation, however, recent findings of gene homologs involved in RA signaling in the genomes of a wide variety of non-chordate animals, including ambulacrarians (sea urchins and acorn worms) and lophotrochozoans (annelids and mollusks), challenged this traditional view and suggested that the RA signaling pathway might have a more ancient evolutionary origin than previously thought. In this chapter, we discuss the evolutionary history of the RA signaling pathway, and more particularly of the RARs, which might have experienced independent gene losses and duplications in different animal lineages. In sum, the available data reveal novel insights into the origin of the RA signaling pathway as well as into the evolutionary history of the RARs.

Phosphorylation of the retinoic acid receptor RARγ2 is crucial for the neuronal differentiation of mouse embryonic stem cells

Retinoic acid (RA) plays key roles in cell differentiation and growth arrest by activating nuclear RA receptors (RARs) (α, β and γ), which are ligand-dependent transcription factors. RARs are also phosphorylated in response to RA. Here, we investigated the in vivo relevance of the phosphorylation of RARs during RA-induced neuronal differentiation of mouse embryonic stem cells (mESCs). Using ESCs where the genes encoding each RAR subtype had been inactivated, and stable rescue lines expressing RARs mutated in phospho-acceptor sites, we show that RA-induced neuronal differentiation involves RARγ2 and requires RARγ2 phosphorylation. By gene expression profiling, we found that the phosphorylated form of RARγ2 regulates a small subset of genes through binding an unusual RA response element consisting of two direct repeats with a seven-base-pair spacer. These new findings suggest an important role for RARγ phosphorylation during cell differentiation and pave the way for further investigations during embryonic development.

Retinoic acid signaling and mouse embryonic stem cell differentiation: Cross talk between genomic and non-genomic effects of RA

Retinoic acid (RA), the active derivative of vitamin A, a fat-soluble vitamin, plays key roles in cell growth and differentiation by activating nuclear receptors, RARs (α, β and γ), which are ligand dependent regulators of transcription. The past years highlighted several novelties in the field that increased the complexity of RA effects. Indeed, in addition to its classical genomic effects, RA also has extranuclear and non-transcriptional effects. RA induces the rapid and transient activation of kinase cascades, which are integrated in the nucleus via the phosphorylation of RARs at a conserved serine residue located in the N-terminal domain and their coregulators. In order to investigate the relevance of RARs' phosphorylation in cell differentiation, mouse embryonic stem (mES) cells were used as a model. When treated with RA, these pluripotent cells give rise to neuronal cells. Cells invalidated for each RAR were generated as well as stable rescue lines expressing RARs mutated in phosphor acceptor sites. Such a strategy revealed that RA-induced neuronal differentiation involves the RARγ2 subtype and requires RARγ2 phosphorylation. Moreover, in gene expression profiling experiments, the phosphorylated form of RARγ2 was found to regulate a small subset of genes through binding a novel RA response element consisting of two direct repeats with a 7 base pair spacer. These new findings suggest an important role for RAR phosphorylation during cell differentiation, and pave the way for further investigations with other cell types and during embryonic development. This article is part of a Special Issue entitled Linking transcription to physiology in lipodomics.

Characterization of genomic vitamin D receptor binding sites through chromatin looping and opening

The vitamin D receptor (VDR) is a transcription factor that mediates the genomic effects of 1α,25-dihydroxyvitamin D₃ (1,25(OH)₂D₃). Genome-wide there are several thousand binding sites and hundreds of primary 1,25(OH)₂D₃ target genes, but their functional relation is largely elusive. In this study, we used ChIA-PET data of the transcription factor CTCF in combination with VDR ChIP-seq data, in order to map chromatin domains containing VDR binding sites. In total, we found 1,599 such VDR containing chromatin domains and studied in THP-1 human monocytic leukemia cells four representatives of them. Our combined ChIP-seq and FAIRE-seq time course data showed that each of these four domains contained a master VDR binding site, where an increase of VDR binding pairs with 1,25(OH)₂D₃-promoted chromatin opening and the presence of a highly significant DR3-type sequence below the peak summit. These sites differed in their relative VDR binding but not in their kinetics, while other loci either had a weaker and delayed VDR association or could not be confirmed at all. All studied chromatin domains contained at least one primary 1,25(OH)₂D₃ target gene demonstrating a characteristic slope of mRNA increase, while neighboring genes responded delayed, if at all. In conclusion, the observation of ligand-inducible VDR binding and chromatin opening combined with a DR3-type sequence highlighted genome-wide 160 VDR loci that have within their chromatin domain a more than 4-fold increased likelihood to identify a primary 1,25(OH)₂D₃ target gene than in the vicinity of other genomic VDR binding sites.

Large-scale quality analysis of published ChIP-seq data

ChIP-seq has become the primary method for identifying in vivo protein-DNA interactions on a genome-wide scale, with nearly 800 publications involving the technique appearing in PubMed as of December 2012. Individually and in aggregate, these data are an important and information-rich resource. However, uncertainties about data quality confound their use by the wider research community. Recently, the Encyclopedia of DNA Elements (ENCODE) project developed and applied metrics to objectively measure ChIP-seq data quality. The ENCODE quality analysis was useful for flagging datasets for closer inspection, eliminating or replacing poor data, and for driving changes in experimental pipelines. There had been no similarly systematic quality analysis of the large and disparate body of published ChIP-seq profiles. Here, we report a uniform analysis of vertebrate transcription factor ChIP-seq datasets in the Gene Expression Omnibus (GEO) repository as of April 1, 2012. The majority (55%) of datasets scored as being highly successful, but a substantial minority (20%) were of apparently poor quality, and another ∼25% were of intermediate quality. We discuss how different uses of ChIP-seq data are affected by specific aspects of data quality, and we highlight exceptional instances for which the metric values should not be taken at face value. Unexpectedly, we discovered that a significant subset of control datasets (i.e., no immunoprecipitation and mock immunoprecipitation samples) display an enrichment structure similar to successful ChIP-seq data. This can, in turn, affect peak calling and data interpretation. Published datasets identified here as high-quality comprise a large group that users can draw on for large-scale integrated analysis. In the future, ChIP-seq quality assessment similar to that used here could guide experimentalists at early stages in a study, provide useful input in the publication process, and be used to stratify ChIP-seq data for different community-wide uses.

Genes involved in cell adhesion and signaling: a new repertoire of retinoic acid receptor target genes in mouse embryonic fibroblasts

Nuclear retinoic acid (RA) receptors (RARα, β and γ) are ligand-dependent transcription factors that regulate the expression of a battery of genes involved in cell differentiation and proliferation. They are also phosphoproteins and we previously showed the importance of their phosphorylation in their transcriptional activity. In the study reported here, we conducted a genome-wide analysis of the genes that are regulated by RARs in mouse embryonic fibroblasts (MEFs) by comparing wild-type MEFs to MEFs lacking the three RARs. We found that in the absence of RA, RARs control the expression of several gene transcripts associated with cell adhesion. Consequently the knockout MEFs are unable to adhere and to spread on substrates and they display a disrupted network of actin filaments, compared with the WT cells. In contrast, in the presence of the ligand, RARs control the expression of other genes involved in signaling and in RA metabolism. Taking advantage of rescue cell lines expressing the RARα or RARγ subtypes (either wild-type or mutated at the N-terminal phosphorylation sites) in the null background, we found that the expression of RA-target genes can be controlled either by a specific single RAR or by a combination of RAR isotypes, depending on the gene. We also selected genes that require the phosphorylation of the receptors for their regulation by RA. Our results increase the repertoire of genes that are regulated by RARs and highlight the complexity and diversity of the transcriptional programs regulated by RARs, depending on the gene.

History of retinoic acid receptors

The discovery of retinoic acid receptors arose from research into how vitamins are essential for life. Early studies indicated that Vitamin A was metabolized into an active factor, retinoic acid (RA), which regulates RNA and protein expression in cells. Each step forward in our understanding of retinoic acid in human health was accomplished by the development and application of new technologies. Development cDNA cloning techniques and discovery of nuclear receptors for steroid hormones provided the basis for identification of two classes of retinoic acid receptors, RARs and RXRs, each of which has three isoforms, α, β and ɣ. DNA manipulation and crystallographic studies revealed that the receptors contain discrete functional domains responsible for binding to DNA, ligands and cofactors. Ligand binding was shown to induce conformational changes in the receptors that cause release of corepressors and recruitment of coactivators to create functional complexes that are bound to consensus promoter DNA sequences called retinoic acid response elements (RAREs) and that cause opening of chromatin and transcription of adjacent genes. Homologous recombination technology allowed the development of mice lacking expression of retinoic acid receptors, individually or in various combinations, which demonstrated that the receptors exhibit vital, but redundant, functions in fetal development and in vision, reproduction, and other functions required for maintenance of adult life. More recent advancements in sequencing and proteomic technologies reveal the complexity of retinoic acid receptor involvement in cellular function through regulation of gene expression and kinase activity. Future directions will require systems biology approaches to decipher how these integrated networks affect human stem cells, health, and disease.

Integrative genomics to dissect retinoid functions

Retinoids and rexinoids, as all other ligands of the nuclear receptor (NR) family, act as ligand-regulated trans-acting transcription factors that bind to cis-acting DNA regulatory elements in the promoter regions of target genes (for reviews see [12, 22, 23, 26, 36]). Ligand binding modulates the communication functions of the receptor with the intracellular environment, which essentially entails receptor-protein and receptor-DNA or receptor-chromatin interactions. In this communication network, the receptor simultaneously serves as both intracellular sensor and regulator of cell/organ functions. Receptors are "intelligent" mediators of the information encoded in the chemical structure of a nuclear receptor ligand, as they interpret this information in the context of cellular identity and cell-physiological status and convert it into a dynamic chain of receptor-protein and receptor-DNA interactions. To process input and output information, they are composed of a modular structure with several domains that have evolved to exert particular molecular recognition functions. As detailed in other chapters in this volume, the main functional domains are the DNA-binding (DBD) and ligand-binding (LBD) [5-7, 38, 56, 71]. The LBD serves as a dual input-output information processor. Inputs, such as ligand binding or receptor phosphorylations, induce allosteric changes in receptor surfaces that serve as docking sites for outputs, such as subunits of transcription and epigenetic machineries or enzyme complexes. The complexity of input and output signals and their interdependencies is far from being understood.