Protocol for bulk and single-nuclei chromatin accessibility quantification in mouse liver tissue.

The accessibility of different chromatin regions to transcription factors and other DNA-binding proteins is a critical determinant of cell function. Here, we detail a modified assay for transposase-accessible chromatin sequencing (ATAC-seq) protocol which measures chromatin accessibility genome wide. We describe nuclei isolation, tagmentation, PCR amplification, and pre- and post-sequencing quality control. Our protocol is optimized for the liver, a tissue where nuclei isolation requires distinct steps. We provide two detailed vignettes: one for bulk ATAC-seq and another for single-nuclei ATAC-seq.

Dynamic chromatin accessibility during nutritional iron overload reveals a BMP6-independent induction of cell cycle genes.

Iron is essential to organism physiology as it participates in numerous biological processes including oxygen transport, respiration, and erythropoiesis. Although iron is critical to physiology, excess iron is toxic to cells and tissues due to generation of reactive oxygen species. Therefore, well-kept iron homeostasis is a mainstay of proper cell and organ function. Iron overload disorders, caused by nutritional or genetic factors, contribute to many pathologies such as diabetes, non-alcoholic steatohepatitis and hepatocellular carcinoma. The liver is not only vulnerable to the effects of iron overload, it is also the major organ controlling iron homeostasis. During iron overload, Bone Morphogenic Protein (BMP) levels increase and initiate a hepatic response aimed at lowering iron levels. The transcriptional effects of iron overload are not well-characterized and the underlining enhancer regulation is uncharted. Here, we profiled the liver's transcriptome and chromatin accessibility following nutritional iron overload. We found marked changes in gene expression and enhancer accessibility following iron overload. Surprisingly, 16% of genes induced following iron overload participate in propagating the cell cycle. Induction of cell cycle genes was independent of BMP. Genome-wide enhancer landscape profiling revealed hundreds of enhancers with altered activity following iron overload. Characterization of transcription factor motifs and footprints in iron-regulated enhancers showed a role for the Activator Protein 1 (AP-1) transcription factor in promoting cell cycle-related transcription. In summary, we found that the transcriptional program at play during iron overload is bifurcated in which BMP signaling controls iron homeostasis genes while an AP-1-driven program controls cell cycle genes.

Hormone-controlled cooperative binding of transcription factors drives synergistic induction of fasting-regulated genes.

During fasting, hepatocytes produce glucose in response to hormonal signals. Glucagon and glucocorticoids are principal fasting hormones that cooperate in regulating glucose production via gluconeogenesis. However, how these hormone signals are integrated and interpreted to a biological output is unknown. Here, we use genome-wide profiling of gene expression, enhancer dynamics and transcription factor (TF) binding in primary mouse hepatocytes to uncover the mode of cooperation between glucagon and glucocorticoids. We found that compared to a single treatment with each hormone, a dual treatment directs hepatocytes to a pro-gluconeogenic gene program by synergistically inducing gluconeogenic genes. The cooperative mechanism driving synergistic gene expression is based on 'assisted loading' whereby a glucagon-activated TF (cAMP responsive element binding protein; CREB) leads to enhancer activation which facilitates binding of the glucocorticoid receptor (GR) upon glucocorticoid stimulation. Glucagon does not only activate single enhancers but also activates enhancer clusters, thereby assisting the loading of GR also across enhancer units within the cluster. In summary, we show that cells integrate extracellular signals by an enhancer-specific mechanism: one hormone-activated TF activates enhancers, thereby assisting the loading of a TF stimulated by a second hormone, leading to synergistic gene induction and a tailored transcriptional response to fasting.

The adipokine FABP4 is a key regulator of neonatal glucose homeostasis.

During pregnancy, fetal glucose production is suppressed, with rapid activation immediately postpartum. Fatty acid-binding protein 4 (FABP4) was recently demonstrated as a regulator of hepatic glucose production and systemic metabolism in animal models. Here, we studied the role of FABP4 in regulating neonatal glucose hemostasis. Serum samples were collected from pregnant women with normoglycemia or gestational diabetes at term, from the umbilical circulation, and from the newborns within 6 hours of life. The level of FABP4 was higher in the fetal versus maternal circulation, with a further rise in neonates after birth of approximately 3-fold. Neonatal FABP4 inversely correlated with blood glucose, with an approximately 10-fold increase of FABP4 in hypoglycemic neonates. When studied in mice, blood glucose of 12-hour-old WT, Fabp4-/+, and Fabp4-/- littermate mice was 59 ± 13 mg/dL, 50 ± 11 mg/dL, and 43 ± 11 mg/dL, respectively. Similar to our observations in humans, FABP4 levels in WT mouse neonates were approximately 8-fold higher compared with those in adult mice. RNA sequencing of the neonatal liver suggested altered expression of multiple glucagon-regulated pathways in Fabp4-/- mice. Indeed, Fabp4-/- liver glycogen was inappropriately intact, despite a marked hypoglycemia, with rapid restoration of normoglycemia upon injection of recombinant FABP4. Our data suggest an important biological role for the adipokine FABP4 in the orchestrated regulation of postnatal glucose metabolism.

Fasting Hormones Synergistically Induce Amino Acid Catabolism Genes to Promote Gluconeogenesis.

BACKGROUND & AIMS: Gluconeogenesis from amino acids (AAs) maintains glucose homeostasis during fasting. Although glucagon is known to regulate AA catabolism, the contribution of other hormones to it and the scope of transcriptional regulation dictating AA catabolism are unknown. We explored the role of the fasting hormones glucagon and glucocorticoids in transcriptional regulation of AA catabolism genes and AA-dependent gluconeogenesis. METHODS: We tested the RNA expression of AA catabolism genes and glucose production in primary mouse hepatocytes treated with fasting hormones (glucagon, corticosterone) and feeding hormones (insulin, fibroblast growth factor 19). We analyzed genomic data of chromatin accessibility and chromatin immunoprecipitation in mice and primary mouse hepatocytes. We performed chromatin immunoprecipitation in livers of fasted mice to show binding of cAMP responsive element binding protein (CREB) and the glucocorticoid receptor (GR). RESULTS: Fasting induced the expression of 31 genes with various roles in AA catabolism. Of them, 15 were synergistically induced by co-treatment of glucagon and corticosterone. Synergistic gene expression relied on the activity of both CREB and GR and was abolished by treatment with either insulin or fibroblast growth factor 19. Enhancers adjacent to synergistically induced genes became more accessible and were bound by CREB and GR on fasting. Akin to the gene expression pattern, gluconeogenesis from AAs was synergistically induced by glucagon and corticosterone in a CREB- and GR-dependent manner. CONCLUSIONS: Transcriptional regulation of AA catabolism genes during fasting is widespread and is driven by glucagon (via CREB) and corticosterone (via GR). Glucose production in hepatocytes is also synergistically augmented, showing that glucagon alone is insufficient in fully activating gluconeogenesis.

Protocol for Primary Mouse Hepatocyte Isolation.

Primary hepatocytes are a vital tool in various biomedical research disciplines, serving as an ex vivo model for liver physiology. Obtaining high yields of viable primary mouse hepatocytes is technically challenging, limiting their use. Here, we present an improved protocol based on the classic two-step collagenase perfusion technique. The liver is washed by perfusion, hepatocytes are dissociated by collagenase, separated from other cells, and cultured. This protocol was optimized to significantly reduce procedure duration and improve hepatocyte yield and viability.

The Three Ds of Transcription Activation by Glucagon: Direct, Delayed, and Dynamic.

Upon lowered blood glucose occurring during fasting, glucagon is secreted from pancreatic islets, exerting various metabolic effects to normalize glucose levels. A considerable portion of these effects is mediated by glucagon-activated transcription factors (TFs) in liver. Glucagon directly activates several TFs via immediate cyclic adenosine monophosphate (cAMP)- and calcium-dependent signaling events. Among these TFs, cAMP response element-binding protein (CREB) is a major factor. CREB recruits histone-modifying enzymes and cooperates with other TFs on the chromatin template to increase the rate of gene transcription. In addition to direct signal transduction, the transcriptional effects of glucagon are also influenced by dynamic TF cross talk. Specifically, assisted loading of one TF by a companion TF leads to increased binding and activity. Lastly, transcriptional regulation by glucagon is also exerted by TF cascades by which a primary TF induces the gene expression of secondary TFs that bring about their activity a few hours after the initial glucagon signal. This mechanism of a delayed response may be instrumental in establishing the temporal organization of the fasting response by which distinct metabolic events separate early from prolonged fasting. In this mini-review, we summarize recent advances and critical discoveries in glucagon-dependent gene regulation with a focus on direct TF activation, dynamic TF cross talk, and TF cascades.

Dynamic enhancer function in the chromatin context.

Enhancers serve as critical regulatory elements in higher eukaryotic cells. The characterization of enhancer function has evolved primarily from genome-wide methodologies, including chromatin immunoprecipitation (ChIP-seq), DNase-I hypersensitivity (DNase-seq), digital genomic footprinting (DGF), and the chromosome conformation capture techniques (3C, 4C, and Hi-C). These population-based assays average signals across millions of cells and lead to enhancer models characterized by static and sequential binding. More recently, fluorescent microscopy techniques, including fluorescence recovery after photobleaching, fluorescence correlation spectroscopy, and single molecule tracking (SMT), reveal a highly dynamic binding behavior for these factors in live cells. Furthermore, a refined analysis of genomic footprinting suggests that many transcription factors leave minimal or no footprints in chromatin, even when present and active in a given cell type. In this study, we review the implications of these new approaches for an accurate understanding of enhancer function in real time. In vivo SMT, in particular, has recently evolved as a promising methodology to probe enhancer function in live cells. Integration of findings from the many approaches now employed in the study of enhancer function suggest a highly dynamic view for the action of enhancer activating factors, viewed on a time scale of milliseconds to seconds, rather than minutes to hours. WIREs Syst Biol Med 2018, 10:e1390. doi: 10.1002/wsbm.1390 This article is categorized under: Analytical and Computational Methods > Computational Methods Laboratory Methods and Technologies > Genetic/Genomic Methods Laboratory Methods and Technologies > Imaging.

Synergistic gene expression during the acute phase response is characterized by transcription factor assisted loading.

The cytokines interleukin 1ß and 6 (IL-1ß, IL-6) mediate the acute phase response (APR). In liver, they regulate the secretion of acute phase proteins. Using RNA-seq in primary hepatocytes, we show that these cytokines regulate transcription in a bifurcated manner, leading to both synergistic and antagonistic gene expression. By mapping changes in enhancer landscape and transcription factor occupancy (using ChIP-seq), we show that synergistic gene induction is achieved by assisted loading of STAT3 on chromatin by NF-κB. With IL-6 treatment alone, STAT3 does not efficiently bind 20% of its coordinated binding sites. In the presence of IL-1ß, NF-κB is activated, binds a subset of enhancers and primes their activity, as evidenced by increasing H3K27ac. This facilitates STAT3 binding and synergistic gene expression. Our findings reveal an enhancer-specific crosstalk whereby NF-κB enables STAT3 binding at some enhancers while perturbing it at others. This model reconciles seemingly contradictory reports of NF-κB-STAT3 crosstalk.

Bivariate Genomic Footprinting Detects Changes in Transcription Factor Activity.

In response to activating signals, transcription factors (TFs) bind DNA and regulate gene expression. TF binding can be measured by protection of the bound sequence from DNase digestion (i.e., footprint). Here, we report that 80% of TF binding motifs do not show a measurable footprint, partly because of a variable cleavage pattern within the motif sequence. To more faithfully portray the effect of TFs on chromatin, we developed an algorithm that captures two TF-dependent effects on chromatin accessibility: footprinting and motif-flanking accessibility. The algorithm, termed bivariate genomic footprinting (BaGFoot), efficiently detects TF activity. BaGFoot is robust to different accessibility assays (DNase-seq, ATAC-seq), all examined peak-calling programs, and a variety of cut bias correction approaches. BaGFoot reliably predicts TF binding and provides valuable information regarding the TFs affecting chromatin accessibility in various biological systems and following various biological events, including in cases where an absolute footprint cannot be determined.

Transcription factor assisted loading and enhancer dynamics dictate the hepatic fasting response.

Fasting elicits transcriptional programs in hepatocytes leading to glucose and ketone production. This transcriptional program is regulated by many transcription factors (TFs). To understand how this complex network regulates the metabolic response to fasting, we aimed at isolating the enhancers and TFs dictating it. Measuring chromatin accessibility revealed that fasting massively reorganizes liver chromatin, exposing numerous fasting-induced enhancers. By utilizing computational methods in combination with dissecting enhancer features and TF cistromes, we implicated four key TFs regulating the fasting response: glucocorticoid receptor (GR), cAMP responsive element binding protein 1 (CREB1), peroxisome proliferator activated receptor alpha (PPARA), and CCAAT/enhancer binding protein beta (CEBPB). These TFs regulate fuel production by two distinctly operating modules, each controlling a separate metabolic pathway. The gluconeogenic module operates through assisted loading, whereby GR doubles the number of sites occupied by CREB1 as well as enhances CREB1 binding intensity and increases accessibility of CREB1 binding sites. Importantly, this GR-assisted CREB1 binding was enhancer-selective and did not affect all CREB1-bound enhancers. Single-molecule tracking revealed that GR increases the number and DNA residence time of a portion of chromatin-bound CREB1 molecules. These events collectively result in rapid synergistic gene expression and higher hepatic glucose production. Conversely, the ketogenic module operates via a GR-induced TF cascade, whereby PPARA levels are increased following GR activation, facilitating gradual enhancer maturation next to PPARA target genes and delayed ketogenic gene expression. Our findings reveal a complex network of enhancers and TFs that dynamically cooperate to restore homeostasis upon fasting.

TM7SF3, a novel p53-regulated homeostatic factor, attenuates cellular stress and the subsequent induction of the unfolded protein response.

Earlier reported small interfering RNA (siRNA) high-throughput screens, identified seven-transmembrane superfamily member 3 (TM7SF3) as a novel inhibitor of pancreatic ß-cell death. Here we show that TM7SF3 maintains protein homeostasis and promotes cell survival through attenuation of ER stress. Overexpression of TM7SF3 inhibits caspase 3/7 activation. In contrast, siRNA-mediated silencing of TM7SF3 accelerates ER stress and activation of the unfolded protein response (UPR). This involves inhibitory phosphorylation of eukaryotic translation initiation factor 2α activity and increased expression of activating transcription factor-3 (ATF3), ATF4 and C/EBP homologous protein, followed by induction of apoptosis. This process is observed both in human pancreatic islets and in a number of cell lines. Some of the effects of TM7SF3 silencing are evident both under basal conditions, in otherwise untreated cells, as well as under different stress conditions induced by thapsigargin, tunicamycin or a mixture of pro-inflammatory cytokines (tumor necrosis factor alpha, interleukin-1 beta and interferon gamma). Notably, TM7SF3 is a downstream target of p53: activation of p53 by Nutlin increases TM7SF3 expression in a time-dependent manner, although silencing of p53 abrogates this effect. Furthermore, p53 is found in physical association with the TM7SF3 promoter. Interestingly, silencing of TM7SF3 promotes p53 activity, suggesting the existence of a negative-feedback loop, whereby p53 promotes expression of TM7SF3 that acts to restrict p53 activity. Our findings implicate TM7SF3 as a novel p53-regulated pro-survival homeostatic factor that attenuates the development of cellular stress and the subsequent induction of the UPR.

Steroid Receptors Reprogram FoxA1 Occupancy through Dynamic Chromatin Transitions.

The estrogen receptor (ER), glucocorticoid receptor (GR), and forkhead box protein 1 (FoxA1) are significant factors in breast cancer progression. FoxA1 has been implicated in establishing ER-binding patterns though its unique ability to serve as a pioneer factor. However, the molecular interplay between ER, GR, and FoxA1 requires further investigation. Here we show that ER and GR both have the ability to alter the genomic distribution of the FoxA1 pioneer factor. Single-molecule tracking experiments in live cells reveal a highly dynamic interaction of FoxA1 with chromatin in vivo. Furthermore, the FoxA1 factor is not associated with detectable footprints at its binding sites throughout the genome. These findings support a model wherein interactions between transcription factors and pioneer factors are highly dynamic. Moreover, at a subset of genomic sites, the role of pioneer can be reversed, with the steroid receptors serving to enhance binding of FoxA1.

Transcriptional and Chromatin Regulation during Fasting - The Genomic Era.

An elaborate metabolic response to fasting is orchestrated by the liver and is heavily reliant on transcriptional regulation. In response to hormones (glucagon, glucocorticoids) many transcription factors (TFs) are activated and regulate various genes involved in metabolic pathways aimed at restoring homeostasis: gluconeogenesis, fatty acid oxidation, ketogenesis, and amino acid shuttling. We summarize recent discoveries regarding fasting-related TFs with an emphasis on genome-wide binding patterns. Collectively, the findings we discuss reveal a large degree of cooperation between TFs during fasting that occurs at motif-rich DNA sites bound by a combination of TFs. These new findings implicate transcriptional and chromatin regulation as major determinants of the response to fasting and unravels the complex, multi-TF nature of this response.

p53 promotes the expression of gluconeogenesis-related genes and enhances hepatic glucose production.

BACKGROUND: The p53 tumor suppressor protein is a transcription factor that initiates transcriptional programs aimed at inhibiting carcinogenesis. p53 represses metabolic pathways that support tumor development (such as glycolysis and the pentose phosphate pathway (PPP)) and enhances metabolic pathways that are considered counter-tumorigenic such as fatty acid oxidation. FINDINGS: In an attempt to comprehensively define metabolic pathways regulated by p53, we performed two consecutive high-throughput analyses in human liver-derived cells with varying p53 statuses. A gene expression microarray screen followed by constraint-based modeling (CBM) predicting metabolic changes imposed by the transcriptomic changes suggested a role for p53 in enhancing gluconeogenesis (de novo synthesis of glucose). Examining glucogenic gene expression revealed a p53-dependent induction of genes involved in both gluconeogenesis (G6PC, PCK2) and in supplying glucogenic precursors (glycerol kinase (GK), aquaporin 3 (AQP3), aquaporin 9 (AQP9) and glutamic-oxaloacetic transaminase 1 (GOT1)). Accordingly, p53 augmented hepatic glucose production (HGP) in both human liver cells and primary mouse hepatocytes. CONCLUSIONS: These findings portray p53 as a novel regulator of glucose production. By facilitating glucose export, p53 may prevent it from being shunted to pro-cancerous pathways such as glycolysis and the PPP. Thus, our findings suggest a metabolic pathway through which p53 may inhibit tumorigenesis.

'Cancer associated fibroblasts'--more than meets the eye.

Cancer associated fibroblasts (CAFs) are a subpopulation of cells that reside within the tumor microenvironment and promotes the transformation process by encouraging tumor growth, angiogenesis, inflammation, and metastasis. CAF-specific proteins serve as both prognostic markers and targets for anticancer drugs. With the growing interest in CAFs, several controversial issues have been raised, including the genomic landscape of these cells, the identity of specific markers, and their cell of origin. Here, we tackle these debated issues and put forward a new definition for 'CAF' as a cell 'state' rather than a cell type. We hope this conceptualization can resolve the ongoing discrepancies revolving around CAF research and aid in designing better anti-cancer treatment strategies.

Mutant p53 attenuates the anti-tumorigenic activity of fibroblasts-secreted interferon beta.

Mutations in the p53 tumor suppressor protein are highly frequent in tumors and often endow cells with tumorigenic capacities. We sought to examine a possible role for mutant p53 in the cross-talk between cancer cells and their surrounding stroma, which is a crucial factor affecting tumor outcome. Here we present a novel model which enables individual monitoring of the response of cancer cells and stromal cells (fibroblasts) to co-culturing. We found that fibroblasts elicit the interferon beta (IFNß) pathway when in contact with cancer cells, thereby inhibiting their migration. Mutant p53 in the tumor was able to alleviate this response via SOCS1 mediated inhibition of STAT1 phosphorylation. IFNß on the other hand, reduced mutant p53 RNA levels by restricting its RNA stabilizer, WIG1. These data underscore mutant p53 oncogenic properties in the context of the tumor microenvironment and suggest that mutant p53 positive cancer patients might benefit from IFNß treatment.

Chemotherapeutic agents induce the expression and activity of their clearing enzyme CYP3A4 by activating p53.

Cytochrome P450 (P450) enzymes are abundantly expressed in the human liver where they hydroxylate organic substrates. In a microarray screen performed in human liver cells, we found a group of eleven P450 genes whose expression was induced by p53 (CYP3A4, CYP3A43, CYP3A5, CYP3A7, CYP4F2, CYP4F3, CYP4F11, CYP4F12, CYP19A1, CYP21A2 and CYP24A1). The mode of regulation of four representative genes (CYP3A4, CYP3A7, CYP4F2 and CYP4F3) was further characterized. The genes were induced in a p53-dependent manner in HepG2 and Huh6 cells (both are cancer-derived human liver cells) and in primary liver cells isolated from human donors. Furthermore, p53 was found to bind to p53-responsive elements in the genes' DNA-regulatory regions and to enhance their transcription in a reporter gene assay. Importantly, when p53 was activated following the administration of either of three different anticancer chemotherapeutic agents (cisplatin, etoposide or doxorubicin), it was able to induce CYP3A genes, which are the main factors in systemic clearance of these agents. Finally, the p53-dependent induction of P450 genes following either Nutlin or chemotherapy treatment led to enhanced P450 enzymatic activity. Thus, in addition to the well-established role of p53 at the tumor site, our data unravels a novel function of hepatic p53 in inducing P450 enzymes and position p53 as a major factor in the hepatic response to xenobiotic and metabolic signals. Importantly, this study reveals a novel pathway for the induction of CYP3As by their substrates through p53, warranting the need for careful consideration when designing systemically administered chemotherapeutic regimens.

Mutant p53R273H attenuates the expression of phase 2 detoxifying enzymes and promotes the survival of cells with high levels of reactive oxygen species.

Uncontrolled accumulation of reactive oxygen species (ROS) causes oxidative stress and induces harmful effects. Both high ROS levels and p53 mutations are frequent in human cancer. Mutant p53 forms are known to actively promote malignant growth. However, no mechanistic details are known about the contribution of mutant p53 to excessive ROS accumulation in cancer cells. Herein, we examine the effect of p53(R273H), a commonly occurring mutated p53 form, on the expression of phase 2 ROS-detoxifying enzymes and on the ability of cells to readopt a reducing environment after exposure to oxidative stress. Our data suggest that p53(R273H) mutant interferes with the normal response of human cells to oxidative stress. We show here that, upon oxidative stress, mutant p53(R273H) attenuates the activation and function of NF-E2-related factor 2 (NRF2), a transcription factor that induces the antioxidant response. This effect of mutant p53 is manifested by decreased expression of phase 2 detoxifying enzymes NQO1 and HO-1 and high ROS levels. These findings were observed in several human cancer cell lines, highlighting the general nature of this phenomenon. The failure of p53(R273H) mutant-expressing cells to restore a reducing oxidative environment was accompanied by increased survival, a known consequence of mutant p53 expression. These activities are attributable to mutant p53(R273H) gain of function and might underlie its well-documented oncogenic nature in human cancer.