TNFα increases the degradation of pyruvate dehydrogenase kinase 4 by the Lon protease to support proinflammatory genes.

The endothelium is a major target of the proinflammatory cytokine, tumor necrosis factor alpha (TNFα). Exposure of endothelial cells (EC) to proinflammatory stimuli leads to an increase in mitochondrial metabolism; however, the function and regulation of elevated mitochondrial metabolism in EC in response to proinflammatory cytokines remain unclear. Studies using high-resolution metabolomics and 13C-glucose and 13C-glutamine labeling flux techniques showed that pyruvate dehydrogenase activity (PDH) and oxidative tricarboxylic acid cycle (TCA) flux are elevated in human umbilical vein ECs in response to overnight (16 h) treatment with TNFα (10 ng/mL). Mechanistic studies indicated that TNFα mediated these metabolic changes via mitochondrial-specific protein degradation of pyruvate dehydrogenase kinase 4 (PDK4, inhibitor of PDH) by the Lon protease via an NF-κB-dependent mechanism. Using RNA sequencing following siRNA-mediated knockdown of the catalytically active subunit of PDH, PDHE1α (PDHA1 gene), we show that PDH flux controls the transcription of approximately one-third of the genes that are up-regulated by TNFα stimulation. Notably, TNFα-induced PDH flux regulates a unique signature of proinflammatory mediators (cytokines and chemokines) but not inducible adhesion molecules. Metabolomics and ChIP sequencing for acetylated modification on lysine 27 of histone 3 (H3K27ac) showed that TNFα-induced PDH flux promotes histone acetylation of specific gene loci via citrate accumulation and ATP-citrate lyase-mediated generation of acetyl CoA. Together, these results uncover a mechanism by which TNFα signaling increases oxidative TCA flux of glucose to support TNFα-induced gene transcription through extramitochondrial acetyl CoA generation and histone acetylation.

Systematic analysis of naturally occurring insertions and deletions that alter transcription factor spacing identifies tolerant and sensitive transcription factor pairs.

Regulation of gene expression requires the combinatorial binding of sequence-specific transcription factors (TFs) at promoters and enhancers. Prior studies showed that alterations in the spacing between TF binding sites can influence promoter and enhancer activity. However, the relative importance of TF spacing alterations resulting from naturally occurring insertions and deletions (InDels) has not been systematically analyzed. To address this question, we first characterized the genome-wide spacing relationships of 73 TFs in human K562 cells as determined by ChIP-seq (chromatin immunoprecipitation sequencing). We found a dominant pattern of a relaxed range of spacing between collaborative factors, including 45 TFs exclusively exhibiting relaxed spacing with their binding partners. Next, we exploited millions of InDels provided by genetically diverse mouse strains and human individuals to investigate the effects of altered spacing on TF binding and local histone acetylation. These analyses suggested that spacing alterations resulting from naturally occurring InDels are generally tolerated in comparison to genetic variants directly affecting TF binding sites. To experimentally validate this prediction, we introduced synthetic spacing alterations between PU.1 and C/EBPß binding sites at six endogenous genomic loci in a macrophage cell line. Remarkably, collaborative binding of PU.1 and C/EBPß at these locations tolerated changes in spacing ranging from 5 bp increase to >30 bp decrease. Collectively, these findings have implications for understanding mechanisms underlying enhancer selection and for the interpretation of non-coding genetic variation.

Breast tumor stiffness instructs bone metastasis via maintenance of mechanical conditioning.

While the immediate and transitory response of breast cancer cells to pathological stiffness in their native microenvironment has been well explored, it remains unclear how stiffness-induced phenotypes are maintained over time after cancer cell dissemination in vivo. Here, we show that fibrotic-like matrix stiffness promotes distinct metastatic phenotypes in cancer cells, which are preserved after transition to softer microenvironments, such as bone marrow. Using differential gene expression analysis of stiffness-responsive breast cancer cells, we establish a multigenic score of mechanical conditioning (MeCo) and find that it is associated with bone metastasis in patients with breast cancer. The maintenance of mechanical conditioning is regulated by RUNX2, an osteogenic transcription factor, established driver of bone metastasis, and mitotic bookmarker that preserves chromatin accessibility at target gene loci. Using genetic and functional approaches, we demonstrate that mechanical conditioning maintenance can be simulated, repressed, or extended, with corresponding changes in bone metastatic potential.

Single-Cell Epigenomics and Functional Fine-Mapping of Atherosclerosis GWAS Loci.

[Figure: see text].

Systems Genetics in Human Endothelial Cells Identifies Non-coding Variants Modifying Enhancers, Expression, and Complex Disease Traits.

The identification of causal variants and mechanisms underlying complex disease traits in humans is important for the progress of human disease genetics; this requires finding strategies to detect functional regulatory variants in disease-relevant cell types. To achieve this, we collected genetic and transcriptomic data from the aortic endothelial cells of up to 157 donors and four epigenomic phenotypes in up to 44 human donors representing individuals of both sexes and three major ancestries. We found thousands of expression quantitative trait loci (eQTLs) at all ranges of effect sizes not detected by the Gene-Tissue Expression Project (GTEx) in human tissues, showing that novel biological relationships unique to endothelial cells (ECs) are enriched in this dataset. Epigenetic profiling enabled discovery of over 3,000 regulatory elements whose activity is modulated by genetic variants that most frequently mutated ETS, AP-1, and NF-kB binding motifs, implicating these motifs as governors of EC regulation. Using CRISPR interference (CRISPRi), allele-specific reporter assays, and chromatin conformation capture, we validated candidate enhancer variants located up to 750 kb from their target genes, VEGFC, FGD6, and KIF26B. Regulatory SNPs identified were enriched in coronary artery disease (CAD) loci, and this result has specific implications for PECAM-1, FES, and AXL. We also found significant roles for EC regulatory variants in modifying the traits pulse pressure, blood protein levels, and monocyte count. Lastly, we present two unlinked SNPs in the promoter of MFAP2 that exhibit pleiotropic effects on human disease traits. Together, this supports the possibility that genetic predisposition for complex disease is manifested through the endothelium.

Transcriptomic and epigenetic mechanisms underlying myeloid diversity in the lung.

The lung is inhabited by resident alveolar and interstitial macrophages as well as monocytic cells that survey lung tissues. Each cell type plays distinct functional roles under homeostatic and inflammatory conditions, but mechanisms establishing their molecular identities and functional potential remain poorly understood. In the present study, systematic evaluation of transcriptomes and open chromatin of alveolar macrophages (AMs), interstitial macrophages (IMs) and lung monocytes from two mouse strains enabled inference of common and cell-specific transcriptional regulators. We provide evidence that these factors drive selection of regulatory landscapes that specify distinct phenotypes of AMs and IMs and entrain qualitatively different responses to toll-like receptor 4 signaling in vivo. These studies reveal a striking divergence in a fundamental innate immune response pathway in AMs and establish a framework for further understanding macrophage diversity in the lung.

MMARGE: Motif Mutation Analysis for Regulatory Genomic Elements.

Cell-specific patterns of gene expression are determined by combinatorial actions of sequence-specific transcription factors at cis-regulatory elements. Studies indicate that relatively simple combinations of lineage-determining transcription factors (LDTFs) play dominant roles in the selection of enhancers that establish cell identities and functions. LDTFs require collaborative interactions with additional transcription factors to mediate enhancer function, but the identities of these factors are often unknown. We have shown that natural genetic variation between individuals has great utility for discovering collaborative transcription factors. Here, we introduce MMARGE (Motif Mutation Analysis of Regulatory Genomic Elements), the first publicly available suite of software tools that integrates genome-wide genetic variation with epigenetic data to identify collaborative transcription factor pairs. MMARGE is optimized to work with chromatin accessibility assays (such as ATAC-seq or DNase I hypersensitivity), as well as transcription factor binding data collected by ChIP-seq. Herein, we provide investigators with rationale for each step in the MMARGE pipeline and key differences for analysis of datasets with different experimental designs. We demonstrate the utility of MMARGE using mouse peritoneal macrophages, liver cells, and human lymphoblastoid cells. MMARGE provides a powerful tool to identify combinations of cell type-specific transcription factors while simultaneously interpreting functional effects of non-coding genetic variation.

Analysis of Genetically Diverse Macrophages Reveals Local and Domain-wide Mechanisms that Control Transcription Factor Binding and Function.

Non-coding genetic variation is a major driver of phenotypic diversity and allows the investigation of mechanisms that control gene expression. Here, we systematically investigated the effects of >50 million variations from five strains of mice on mRNA, nascent transcription, transcription start sites, and transcription factor binding in resting and activated macrophages. We observed substantial differences associated with distinct molecular pathways. Evaluating genetic variation provided evidence for roles of ∼100 TFs in shaping lineage-determining factor binding. Unexpectedly, a substantial fraction of strain-specific factor binding could not be explained by local mutations. Integration of genomic features with chromatin interaction data provided evidence for hundreds of connected cis-regulatory domains associated with differences in transcription factor binding and gene expression. This system and the >250 datasets establish a substantial new resource for investigation of how genetic variation affects cellular phenotypes.

Deleting an Nr4a1 Super-Enhancer Subdomain Ablates Ly6Clow Monocytes while Preserving Macrophage Gene Function.

Mononuclear phagocytes are a heterogeneous family that occupy all tissues and assume numerous roles to support tissue function and systemic homeostasis. Our ability to dissect the roles of individual subsets is limited by a lack of technologies that ablate gene function within specific mononuclear phagocyte sub-populations. Using Nr4a1-dependent Ly6Clow monocytes, we present a proof-of-principle approach that addresses these limitations. Combining ChIP-seq and molecular approaches we identified a single, conserved, sub-domain within the Nr4a1 enhancer that was essential for Ly6Clow monocyte development. Mice lacking this enhancer lacked Ly6Clow monocytes but retained Nr4a1 gene expression in macrophages during steady state and in response to LPS. Because Nr4a1 regulates inflammatory gene expression and differentiation of Ly6Clow monocytes, decoupling these processes allows Ly6Clow monocytes to be studied independently.

Exploiting genomics and natural genetic variation to decode macrophage enhancers.

The mammalian genome contains on the order of a million enhancer-like regions that are required to establish the identities and functions of specific cell types. Here, we review recent studies in immune cells that have provided insight into the mechanisms that selectively activate certain enhancers in response to cell lineage and environmental signals. We describe a working model wherein distinct classes of transcription factors define the repertoire of active enhancers in macrophages through collaborative and hierarchical interactions, and discuss important challenges to this model, specifically providing examples from T cells. We conclude by discussing the use of natural genetic variation as a powerful approach for decoding transcription factor combinations that play dominant roles in establishing the enhancer landscapes, and the potential that these insights have for advancing our understanding of the molecular causes of human disease.

MAFG is a transcriptional repressor of bile acid synthesis and metabolism.

Specific bile acids are potent signaling molecules that modulate metabolic pathways affecting lipid, glucose and bile acid homeostasis, and the microbiota. Bile acids are synthesized from cholesterol in the liver, and the key enzymes involved in bile acid synthesis (Cyp7a1, Cyp8b1) are regulated transcriptionally by the nuclear receptor FXR. We have identified an FXR-regulated pathway upstream of a transcriptional repressor that controls multiple bile acid metabolism genes. We identify MafG as an FXR target gene and show that hepatic MAFG overexpression represses genes of the bile acid synthetic pathway and modifies the biliary bile acid composition. In contrast, loss-of-function studies using MafG(+/-) mice causes de-repression of the same genes with concordant changes in biliary bile acid levels. Finally, we identify functional MafG response elements in bile acid metabolism genes using ChIP-seq analysis. Our studies identify a molecular mechanism for the complex feedback regulation of bile acid synthesis controlled by FXR.

Environment drives selection and function of enhancers controlling tissue-specific macrophage identities.

Macrophages reside in essentially all tissues of the body and play key roles in innate and adaptive immune responses. Distinct populations of tissue macrophages also acquire context-specific functions that are important for normal tissue homeostasis. To investigate mechanisms responsible for tissue-specific functions, we analyzed the transcriptomes and enhancer landscapes of brain microglia and resident macrophages of the peritoneal cavity. In addition, we exploited natural genetic variation as a genome-wide "mutagenesis" strategy to identify DNA recognition motifs for transcription factors that promote common or subset-specific binding of the macrophage lineage-determining factor PU.1. We find that distinct tissue environments drive divergent programs of gene expression by differentially activating a common enhancer repertoire and by inducing the expression of divergent secondary transcription factors that collaborate with PU.1 to establish tissue-specific enhancers. These findings provide insights into molecular mechanisms by which tissue environment influences macrophage phenotypes that are likely to be broadly applicable to other cell types.

Control of VEGF-A transcriptional programs by pausing and genomic compartmentalization.

Vascular endothelial growth factor A (VEGF-A) is a master regulator of angiogenesis, vascular development and function. In this study we investigated the transcriptional regulation of VEGF-A-responsive genes in primary human aortic endothelial cells (HAECs) and human umbilical vein endothelial cells (HUVECs) using genome-wide global run-on sequencing (GRO-Seq). We demonstrate that half of VEGF-A-regulated gene promoters are characterized by a transcriptionally competent paused RNA polymerase II (Pol II). We show that transition into productive elongation is a major mechanism of gene activation of virtually all VEGF-regulated genes, whereas only ∼40% of the genes are induced at the level of initiation. In addition, we report a comprehensive chromatin interaction map generated in HUVECs using tethered conformation capture (TCC) and characterize chromatin interactions in relation to transcriptional activity. We demonstrate that sites of active transcription are more likely to engage in chromatin looping and cell type-specific transcriptional activity reflects the boundaries of chromatin interactions. Furthermore, we identify large chromatin compartments with a tendency to be coordinately transcribed upon VEGF-A stimulation. We provide evidence that these compartments are enriched for clusters of regulatory regions such as super-enhancers and for disease-associated single nucleotide polymorphisms (SNPs). Collectively, these findings provide new insights into mechanisms behind VEGF-A-regulated transcriptional programs in endothelial cells.

25-Hydroxycholesterol activates the integrated stress response to reprogram transcription and translation in macrophages.

25-Hydroxycholesterol (25OHC) is an enzymatically derived oxidation product of cholesterol that modulates lipid metabolism and immunity. 25OHC is synthesized in response to interferons and exerts broad antiviral activity by as yet poorly characterized mechanisms. To gain further insights into the basis for antiviral activity, we evaluated time-dependent responses of the macrophage lipidome and transcriptome to 25OHC treatment. In addition to altering specific aspects of cholesterol and sphingolipid metabolism, we found that 25OHC activates integrated stress response (ISR) genes and reprograms protein translation. Effects of 25OHC on ISR gene expression were independent of liver X receptors and sterol-response element-binding proteins and instead primarily resulted from activation of the GCN2/eIF2α/ATF4 branch of the ISR pathway. These studies reveal that 25OHC activates the integrated stress response, which may contribute to its antiviral activity.

Remodeling of the enhancer landscape during macrophage activation is coupled to enhancer transcription.

Recent studies suggest a hierarchical model in which lineage-determining factors act in a collaborative manner to select and prime cell-specific enhancers, thereby enabling signal-dependent transcription factors to bind and function in a cell-type-specific manner. Consistent with this model, TLR4 signaling primarily regulates macrophage gene expression through a pre-existing enhancer landscape. However, TLR4 signaling also induces priming of ∼3,000 enhancer-like regions de novo, enabling visualization of intermediates in enhancer selection and activation. Unexpectedly, we find that enhancer transcription precedes local mono- and dimethylation of histone H3 lysine 4 (H3K4me1/2). H3K4 methylation at de novo enhancers is primarily dependent on the histone methyltransferases Mll1, Mll2/4, and Mll3 and is significantly reduced by inhibition of RNA polymerase II elongation. Collectively, these findings suggest an essential role of enhancer transcription in H3K4me1/2 deposition at de novo enhancers that is independent of potential functions of the resulting eRNA transcripts.

Identification of CAD candidate genes in GWAS loci and their expression in vascular cells.

Recent genome-wide association studies (GWAS) have identified 35 loci that significantly associate with coronary artery disease (CAD) susceptibility. The majority of the genes represented in these loci have not previously been studied in the context of atherosclerosis. To characterize the roles of these candidate genes in the vessel wall, we determined their expression levels in endothelial, smooth muscle, and macrophage cells isolated from healthy, prelesioned, and lesioned mouse aortas. We also performed expression quantitative locus (eQTL) mapping of these genes in human endothelial cells under control and proatherogenic conditions. Of the 57 genes studied, 31 were differentially expressed in one or more cell types in disease state in mice, and the expression levels of 8 were significantly associated with the CAD SNPs in human cells, 7 of which were also differentially expressed in mice. By integrating human and mouse results, we predict that PPAP2B, GALNT4, MAPKAPK5, TCTN1, SRR, SNF8, and ICAM1 play a causal role in the susceptibility to atherosclerosis through a role in the vasculature. Additionally, we highlight the genetic complexity of a subset of CAD loci through the differential expression of multiple candidate genes per locus and the involvement of genes that lie outside linkage disequilibrium blocks.

Role of phospholipid oxidation products in atherosclerosis.

There is increasing clinical evidence that phospholipid oxidation products (Ox-PL) play a role in atherosclerosis. This review focuses on the mechanisms by which Ox-PL interact with endothelial cells, monocyte/macrophages, platelets, smooth muscle cells, and HDL to promote atherogenesis. In the past few years major progress has been made in identifying these mechanisms. It has been recognized that Ox-PL promote phenotypic changes in these cell types that have long-term consequences for the vessel wall. Individual Ox-PL responsible for specific cellular effects have been identified. A model of the configuration of bioactive truncated Ox-PL within membranes has been developed that demonstrates that the oxidized fatty acid moiety protrudes into the aqueous phase, rendering it accessible for receptor recognition. Receptors and signaling pathways for individual Ox-PL species are now determined and receptor independent signaling pathways identified. The effects of Ox-PL are mediated both by gene regulation and transcription independent processes. It has now become apparent that Ox-PL affects multiple genes and pathways, some of which are proatherogenic and some are protective. However, at concentrations that are likely present in the vessel wall in atherosclerotic lesions, the effects promote atherogenesis. There have also been new insights on enzymes that metabolize Ox-PL and the significance of these enzymes for atherosclerosis. With the knowledge we now have of the regulation and effects of Ox-PL in different vascular cell types, it should be possible to design experiments to test the role of specific Ox-PL on the development of atherosclerosis.

Metalloproteinase processing of HBEGF is a proximal event in the response of human aortic endothelial cells to oxidized phospholipids.

OBJECTIVE: Atherosclerosis is a chronic inflammatory disease initiated by monocyte recruitment and retention in the vessel wall. An important mediator of monocyte endothelial interaction is the chemokine interleukin (IL)-8. The oxidation products of phospholipids, including oxidized 1-palmitoyl-2-arachidonyl-sn-glycerol-3-phosphocholine (Ox-PAPC), accumulate in atherosclerotic lesions and strongly induce IL-8 in human aortic endothelial cells (HAECs). The goal of this study was to identify the proximal events leading to induction of IL-8 by Ox-PAPC in vascular endothelial cells. METHODS AND RESULTS: In a systems genetics analysis of HAECs isolated from 96 different human donors, we showed that heparin-binding EGF-like growth factor (HBEGF) transcript levels are strongly correlated to IL-8 induction by Ox-PAPC. The silencing and overexpression of HBEGF in HAECs confirmed the role of HBEGF in regulating IL-8 expression. HBEGF has been shown to be stored in an inactive form and activation is dependent on processing by a dysintegrin and metalloproteinases (ADAM) to a form that can activate the epidermal growth factor (EGF) receptor. Ox-PAPC was shown to rapidly induce HBEGF processing and EGF receptor activation in HAECs. Using siRNA we identified 3 ADAMs that regulate IL-8 induction and directly demonstrated that Ox-PAPC increases ADAM activity in the cells using a substrate cleavage assay. We provide evidence for one mechanism of Ox-PAPC activation of ADAM involving covalent binding of Ox-PAPC to cysteine on ADAM. Free thiol cysteine analogs showed inhibition of IL-8 induction by Ox-PAPC, and both a cysteine analog and a cell surface thiol blocker strongly inhibited ADAM activity induction by Ox-PAPC. Using microarray analyses, we determined that this ADAM pathway may regulate at least 30% of genes induced by Ox-PAPC in HAECs. CONCLUSIONS: This study is the first report demonstrating a role for the ADAM-HBEGF-EGF receptor axis in Ox-PAPC induction of IL-8 in HAECs. These studies highlight a role for specific ADAMs as initiators of Ox-PAPC action and provide evidence for a role of covalent interaction of Ox-PAPC in activation of ADAMs.