Insulin Resistance in Women Correlates with Chromatin Histone Lysine Acetylation, Inflammatory Signaling, and Accelerated Aging.

BACKGROUND: Epigenetic changes link medical, social, and environmental factors with cardiovascular and kidney disease and, more recently, with cancer. The mechanistic link between metabolic health and epigenetic changes is only starting to be investigated. In our in vitro and in vivo studies, we performed a broad analysis of the link between hyperinsulinemia and chromatin acetylation; our top "hit" was chromatin opening at H3K9ac. METHODS: Building on our published preclinical studies, here, we performed a detailed analysis of the link between insulin resistance, chromatin acetylation, and inflammation using an initial test set of 28 women and validation sets of 245, 22, and 53 women. RESULTS: ChIP-seq identified chromatin acetylation and opening at the genes coding for TNFα and IL6 in insulin-resistant women. Pathway analysis identified inflammatory response genes, NFκB/TNFα-signaling, reactome cytokine signaling, innate immunity, and senescence. Consistent with this finding, flow cytometry identified increased senescent circulating peripheral T-cells. DNA methylation analysis identified evidence of accelerated aging in insulin-resistant vs. metabolically healthy women. CONCLUSIONS: This study shows that insulin-resistant women have increased chromatin acetylation/opening, inflammation, and, perhaps, accelerated aging. Given the role that inflammation plays in cancer initiation and progression, these studies provide a potential mechanistic link between insulin resistance and cancer.

Loss of epigenetic suppression of retrotransposons with oncogenic potential in aging mammary luminal epithelial cells.

A primary function of DNA methylation in mammalian genomes is to repress transposable elements (TEs). The widespread methylation loss that is commonly observed in cancer cells results in the loss of epigenetic repression of TEs. The aging process is similarly characterized by changes to the methylome. However, the impact of these epigenomic alterations on TE silencing and the functional consequences of this have remained unclear. To assess the epigenetic regulation of TEs in aging, we profiled DNA methylation in human mammary luminal epithelial cells (LEps)-a key cell lineage implicated in age-related breast cancers-from younger and older women. We report here that several TE subfamilies function as regulatory elements in normal LEps, and a subset of these display consistent methylation changes with age. Methylation changes at these TEs occurred at lineage-specific transcription factor binding sites, consistent with loss of lineage specificity. Whereas TEs mainly showed methylation loss, CpG islands (CGIs) that are targets of the Polycomb repressive complex 2 (PRC2) show a gain of methylation in aging cells. Many TEs with methylation loss in aging LEps have evidence of regulatory activity in breast cancer samples. We furthermore show that methylation changes at TEs impact the regulation of genes associated with luminal breast cancers. These results indicate that aging leads to DNA methylation changes at TEs that undermine the maintenance of lineage specificity, potentially increasing susceptibility to breast cancer.

Inhibition of DNMT1 methyltransferase activity via glucose-regulated O-GlcNAcylation alters the epigenome.

The DNA methyltransferase activity of DNMT1 is vital for genomic maintenance of DNA methylation. We report here that DNMT1 function is regulated by O-GlcNAcylation, a protein modification that is sensitive to glucose levels, and that elevated O-GlcNAcylation of DNMT1 from high glucose environment leads to alterations to the epigenome. Using mass spectrometry and complementary alanine mutation experiments, we identified S878 as the major residue that is O-GlcNAcylated on human DNMT1. Functional studies in human and mouse cells further revealed that O-GlcNAcylation of DNMT1-S878 results in an inhibition of methyltransferase activity, resulting in a general loss of DNA methylation that preferentially occurs at partially methylated domains (PMDs). This loss of methylation corresponds with an increase in DNA damage and apoptosis. These results establish O-GlcNAcylation of DNMT1 as a mechanism through which the epigenome is regulated by glucose metabolism and implicates a role for glycosylation of DNMT1 in metabolic diseases characterized by hyperglycemia.

Blockade of trans PD-L1 interaction with CD80 augments antitumor immunity.

PD-L1 has two receptors: PD-1 and CD80. Previous reports assumed that PD-L1 and CD80 interacted in trans, but recent reports showed that only cis PD-L1/CD80 interactions existed, and prevention of cis PD-L1/CD80 interactions on antigen-presenting cells (APCs) reduced antitumor immunity via augmenting PD-L1/PD-1 and CD80/CTLA4 interactions between T and APCs. Here, using tumor-bearing mice capable of cis and trans or trans only PD-L1/CD80 interactions, we show that trans PD-L1/CD80 interactions do exist between tumor and T cells, and the effects of trans PD-L1/CD80 interactions require tumor cell expression of MHC-I and T cell expression of CD28. The blockade of PD-L1/CD80 interactions in mice with both cis and trans interactions or with only trans interactions augments antitumor immunity by expanding IFN-γ-producing CD8+ T cells and IFN-γ-dependent NOS2-expressing tumor-associated macrophages. Our studies indicate that although cis and trans PD-L1/CD80 interactions may have opposite effects on antitumor immunity, the net effect of blocking PD-L1/CD80 interactions in vivo augments CD8+ T cell-mediated antitumor immunity.

Epigenetic dosage identifies two major and functionally distinct ß cell subtypes.

The mechanisms that specify and stabilize cell subtypes remain poorly understood. Here, we identify two major subtypes of pancreatic ß cells based on histone mark heterogeneity (ßHI and ßLO). ßHI cells exhibit ∼4-fold higher levels of H3K27me3, distinct chromatin organization and compaction, and a specific transcriptional pattern. ßHI and ßLO cells also differ in size, morphology, cytosolic and nuclear ultrastructure, epigenomes, cell surface marker expression, and function, and can be FACS separated into CD24+ and CD24- fractions. Functionally, ßHI cells have increased mitochondrial mass, activity, and insulin secretion in vivo and ex vivo. Partial loss of function indicates that H3K27me3 dosage regulates ßHI/ßLO ratio in vivo, suggesting that control of ß cell subtype identity and ratio is at least partially uncoupled. Both subtypes are conserved in humans, with ßHI cells enriched in humans with type 2 diabetes. Thus, epigenetic dosage is a novel regulator of cell subtype specification and identifies two functionally distinct ß cell subtypes.

Epigenetic Dysregulation of KCNK9 Imprinting and Triple-Negative Breast Cancer.

Genomic imprinting is an inherited form of parent-of-origin specific epigenetic gene regulation that is dysregulated by poor prenatal nutrition and environmental toxins. KCNK9 encodes for TASK3, a pH-regulated potassium channel membrane protein that is overexpressed in 40% of breast cancer. However, KCNK9 gene amplification accounts for increased expression in <10% of these breast cancers. Here, we showed that KCNK9 is imprinted in breast tissue and identified a differentially methylated region (DMR) controlling its imprint status. Hypomethylation at the DMR, coupled with biallelic expression of KCNK9, occurred in 63% of triple-negative breast cancers (TNBC). The association between hypomethylation and TNBC status was highly significant in African-Americans (p = 0.006), but not in Caucasians (p = 0.70). KCNK9 hypomethylation was also found in non-cancerous tissue from 77% of women at high-risk of developing breast cancer. Functional studies demonstrated that the KCNK9 gene product, TASK3, regulates mitochondrial membrane potential and apoptosis-sensitivity. In TNBC cells and non-cancerous mammary epithelial cells from high-risk women, hypomethylation of the KCNK9 DMR predicts for increased TASK3 expression and mitochondrial membrane potential (p < 0.001). This is the first identification of the KCNK9 DMR in mammary epithelial cells and demonstration that its hypomethylation in breast cancer is associated with increases in both mitochondrial membrane potential and apoptosis resistance. The high frequency of hypomethylation of the KCNK9 DMR in TNBC and non-cancerous breast tissue from high-risk women provides evidence that hypomethylation of the KNCK9 DMR/TASK3 overexpression may serve as a marker of risk and a target for prevention of TNBC, particularly in African American women.

Sequence features of retrotransposons allow for epigenetic variability.

Transposable elements (TEs) are mobile genetic elements that make up a large fraction of mammalian genomes. While select TEs have been co-opted in host genomes to have function, the majority of these elements are epigenetically silenced by DNA methylation in somatic cells. However, some TEs in mice, including the Intracisternal A-particle (IAP) subfamily of retrotransposons, have been shown to display interindividual variation in DNA methylation. Recent work has revealed that IAP sequence differences and strain-specific KRAB zinc finger proteins (KZFPs) may influence the methylation state of these IAPs. However, the mechanisms underlying the establishment and maintenance of interindividual variability in DNA methylation still remain unclear. Here, we report that sequence content and genomic context influence the likelihood that IAPs become variably methylated. IAPs that differ from consensus IAP sequences have altered KZFP recruitment that can lead to decreased KAP1 recruitment when in proximity of constitutively expressed genes. These variably methylated loci have a high CpG density, similar to CpG islands, and can be bound by ZF-CxxC proteins, providing a potential mechanism to maintain this permissive chromatin environment and protect from DNA methylation. These observations indicate that variably methylated IAPs escape silencing through both attenuation of KZFP binding and recognition by ZF-CxxC proteins to maintain a hypomethylated state.

Dietary glutamine supplementation suppresses epigenetically-activated oncogenic pathways to inhibit melanoma tumour growth.

Tumour cells adapt to nutrient deprivation in vivo, yet strategies targeting the nutrient poor microenvironment remain unexplored. In melanoma, tumour cells often experience low glutamine levels, which promote cell dedifferentiation. Here, we show that dietary glutamine supplementation significantly inhibits melanoma tumour growth, prolongs survival in a transgenic melanoma mouse model, and increases sensitivity to a BRAF inhibitor. Metabolomic analysis reveals that dietary uptake of glutamine effectively increases the concentration of glutamine in tumours and its downstream metabolite, αKG, without increasing biosynthetic intermediates necessary for cell proliferation. Mechanistically, we find that glutamine supplementation uniformly alters the transcriptome in tumours. Our data further demonstrate that increase in intra-tumoural αKG concentration drives hypomethylation of H3K4me3, thereby suppressing epigenetically-activated oncogenic pathways in melanoma. Therefore, our findings provide evidence that glutamine supplementation can serve as a potential dietary intervention to block melanoma tumour growth and sensitize tumours to targeted therapy via epigenetic reprogramming.

Hyperinsulinemia promotes aberrant histone acetylation in triple-negative breast cancer.

BACKGROUND: Hyperinsulinemia, the presence of excess insulin relative to glucose in the blood, is considered to be a poor prognostic indicator for patients with triple-negative breast cancer (TNBC). mTOR, a downstream effector of insulin, enhances mitochondrial biogenesis and activity, thereby increasing acetyl-CoA precursors. Increased acetyl-CoA can, in turn, be utilized by nuclear acetyltransferases for histone acetylation, a critical feature of genome regulation. While signaling pathways downstream of insulin have been established for sometime, the effect of insulin on chromatin remains unclear. We hypothesized that hyperinsulinemia-induced metabolic changes lead to genome-wide changes in histone acetylation in TNBC. RESULTS: MDA-MB-231 cells were xenografted into hyperinsulinemic and wild-type mice. Tumors in the hyperinsulinemic mice displayed elevated levels of histone acetylation compared to tumors in normal insulin conditions. We show that insulin treatment in vitro leads to global increase in chromatin-associated histone acetylation, in particular at H3K9, through the PI3K/AKT/mTOR pathway. Genome-wide analyses revealed that most promoter regions have an increase in histone acetylation upon insulin treatment. In addition, insulin induces higher levels of reactive oxygen species and DNA damage foci in cells. CONCLUSIONS: These results demonstrate the impact of hyperinsulinemia on altered gene regulation through chromatin and the importance of targeting hyperinsulinemia-induced processes that lead to chromatin dysfunction in TNBC.

Arginine starvation kills tumor cells through aspartate exhaustion and mitochondrial dysfunction.

Defective arginine synthesis, due to the silencing of argininosuccinate synthase 1 (ASS1), is a common metabolic vulnerability in cancer, known as arginine auxotrophy. Understanding how arginine depletion kills arginine-auxotrophic cancer cells will facilitate the development of anti-cancer therapeutic strategies. Here we show that depletion of extracellular arginine in arginine-auxotrophic cancer cells causes mitochondrial distress and transcriptional reprogramming. Mechanistically, arginine starvation induces asparagine synthetase (ASNS), depleting these cancer cells of aspartate, and disrupting their malate-aspartate shuttle. Supplementation of aspartate, depletion of mitochondria, and knockdown of ASNS all protect the arginine-starved cells, establishing the causal effects of aspartate depletion and mitochondrial dysfunction on the arginine starvation-induced cell death. Furthermore, dietary arginine restriction reduced tumor growth in a xenograft model of ASS1-deficient breast cancer. Our data challenge the view that ASNS promotes homeostasis, arguing instead that ASNS-induced aspartate depletion promotes cytotoxicity, which can be exploited for anti-cancer therapies.

LTRs activated by Epstein-Barr virus-induced transformation of B cells alter the transcriptome.

Endogenous retroviruses (ERVs) are ancient viral elements that have accumulated in the genome through retrotransposition events. Although they have lost their ability to transpose, many of the long terminal repeats (LTRs) that originally flanked full-length ERVs maintain the ability to regulate transcription. While these elements are typically repressed in somatic cells, they can function as transcriptional enhancers and promoters when this repression is lost. Epstein-Barr virus (EBV), which transforms primary B cells into continuously proliferating cells, is a tumor virus associated with lymphomas. We report here that transformation of primary B cells by EBV leads to genome-wide activation of LTR enhancers and promoters. The activation of LTRs coincides with local DNA hypomethylation and binding by transcription factors such as RUNX3, EBF1, and EBNA2. The set of activated LTRs is unique to transformed B cells compared with other cell lines known to have activated LTRs. Furthermore, we found that LTR activation impacts the B cell transcriptome by up-regulating transcripts driven by cryptic LTR promoters. These transcripts include genes important to oncogenesis of Hodgkin lymphoma and other cancers, such as HUWE1/HECTH9 These data suggest that the activation of LTRs by EBV-induced transformation is important to the pathology of EBV-associated cancers. Altogether, our results indicate that EBV-induced transformation of B cells alters endogenous retroviral element activity, thereby impacting host gene regulatory networks and oncogenic potential.

Diabetes Mellitus-Induced Long Noncoding RNA Dnm3os Regulates Macrophage Functions and Inflammation via Nuclear Mechanisms.

Objective- Macrophages play key roles in inflammation and diabetic vascular complications. Emerging evidence implicates long noncoding RNAs in inflammation, but their role in macrophage dysfunction associated with inflammatory diabetic complications is unclear and was therefore investigated in this study. Approach and Results- RNA-sequencing and real-time quantitative PCR demonstrated that a long noncoding RNA Dnm3os (dynamin 3 opposite strand) is upregulated in bone marrow-derived macrophages from type 2 diabetic db/db mice, diet-induced insulin-resistant mice, and diabetic ApoE-/- mice, as well as in monocytes from type 2 diabetic patients relative to controls. Diabetic conditions (high glucose and palmitic acid) induced Dnm3os in mouse and human macrophages. Promoter reporter analysis and chromatin immunoprecipitation assays demonstrated that diabetic conditions induce Dnm3os via NF-κB activation. RNA fluorescence in situ hybridization and real-time quantitative PCRs of subcellular fractions demonstrated nuclear localization and chromatin enrichment of Dnm3os in macrophages. Stable overexpression of Dnm3os in macrophages altered global histone modifications and upregulated inflammation and immune response genes and phagocytosis. Conversely, RNAi-mediated knockdown of Dnm3os attenuated these responses. RNA pull-down assays with macrophage nuclear lysates identified nucleolin and ILF-2 (interleukin enhancer-binding factor 2) as protein binding partners of Dnm3os, which was further confirmed by RNA fluorescence in situ hybridization immunofluorescence. Furthermore, nucleolin levels were decreased in diabetic conditions, and its knockdown enhanced Dnm3os-induced inflammatory gene expression and histone H3K9-acetylation at their promoters. Conclusions- These results demonstrate novel mechanisms involving upregulation of long noncoding RNA Dnm3os, disruption of its interaction with nucleolin, and epigenetic modifications at target genes that promote macrophage inflammatory phenotype in diabetes mellitus. The data could lead to long noncoding RNA-based therapies for inflammatory diabetes mellitus complications.

Estrogens and selective estrogen receptor modulators differentially antagonize Runx2 in ST2 mesenchymal progenitor cells.

Estrogens attenuate bone turnover by inhibiting both osteoclasts and osteoblasts, in part through antagonizing Runx2. Apparently conflicting, stimulatory effects in osteoblast lineage cells, however, sway the balance between bone resorption and bone formation in favor of the latter. Consistent with this dualism, 17ß-estradiol (E2) both stimulates and inhibits Runx2 in a locus-specific manner, and here we provide evidence for such locus-specific regulation of Runx2 by E2 in vivo. We also demonstrate dual, negative and positive, regulation of Runx2-driven alkaline phosphatase (ALP) activity by increasing E2 concentrations in ST2 osteoblast progenitor cells. We further compared the effects of E2 to those of the Selective Estrogen Receptor Modulators (SERMs) raloxifene (ral) and lasofoxifene (las) and the phytoestrogen puerarin. We found that E2 at the physiological concentrations of 0.1-1 nM, as well as ral and las, but not puerarin, antagonize Runx2-driven ALP activity. At ≥10 nM, E2 and puerarin, but not ral or las, stimulate ALP relative to the activity measured at 0.1-1 nM. Contrasting the difference between E2 and SERMs in ST2 cells, they all shared a similar dose-response profile when inhibiting pre-osteoclast proliferation. That ral and las poorly mimic the locus- and concentration-dependent effects of E2 in mesenchymal progenitor cells may help explain their limited clinical efficacy.

JMJD1B Demethylates H4R3me2s and H3K9me2 to Facilitate Gene Expression for Development of Hematopoietic Stem and Progenitor Cells.

The arginine methylation status of histones dynamically changes during many cellular processes, including hematopoietic stem/progenitor cell (HSPC) development. The arginine methyltransferases and the readers that transduce the histone codes have been defined. However, whether arginine demethylation actively occurs in cells and what enzyme demethylates the methylarginine residues during various cellular processes are unknown. We report that JMJD1B, previously identified as a lysine demethylase for H3K9me2, mediates arginine demethylation of H4R3me2s and its intermediate, H4R3me1. We show that demethylation of H4R3me2s and H3K9me2s in promoter regions is correlated with active gene expression. Furthermore, knockout of JMJD1B blocks demethylation of H4R3me2s and/or H3K9me2 at distinct clusters of genes and impairs the activation of genes important for HSPC differentiation and development. Consequently, JMJD1B-/- mice show defects in hematopoiesis. Altogether, our study demonstrates that demethylase-mediated active arginine demethylation process exists in eukaryotes and that JMJD1B demethylates both H4R3me2s and H3K9me2 for epigenetic programming during hematopoiesis.

Chromatin modifications in metabolic disease: Potential mediators of long-term disease risk.

Metabolic diseases such as obesity and diabetes are complex diseases resulting from multiple genetic and environmental factors, such as diet and activity levels. These factors are well known contributors to the development of metabolic diseases. One manner by which environmental factors can influence metabolic disease progression is through modifications to chromatin. These modifications can lead to altered gene regulatory programs, which alters disease risk. Furthermore, there is evidence that parents exposed to environmental factors can influence the metabolic health of offspring, especially if exposures are during intrauterine growth periods. In this review, we outline the evidence that chromatin modifications are associated with metabolic diseases, including diabetes and obesity. We also consider evidence that these chromatin modifications can lead to long-term disease risk and contribute to disease risk for future generations. This article is categorized under: Biological Mechanisms > Metabolism Developmental Biology > Developmental Processes in Health and Disease Physiology > Organismal Responses to Environment.

Regulation of angiotensin II actions by enhancers and super-enhancers in vascular smooth muscle cells.

Angiotensin II (AngII) promotes hypertension and atherosclerosis by activating growth-promoting and pro-inflammatory gene expression in vascular smooth muscle cells (VSMCs). Enhancers and super-enhancers (SEs) play critical roles in driving disease-associated gene expression. However, enhancers/SEs mediating VSMC dysfunction remain uncharacterized. Here, we show that AngII alters vascular enhancer and SE repertoires in cultured VSMCs in vitro, ex vivo, and in AngII-infused mice aortas in vivo. AngII-induced enhancers/SEs are enriched in binding sites for signal-dependent transcription factors and dependent on key signaling kinases. Moreover, CRISPR-Cas9-mediated deletion of candidate enhancers/SEs, targeting SEs with the bromodomain and extra-terminal domain inhibitor JQ1, or knockdown of overlapping long noncoding RNAs (lncRNAs) blocks AngII-induced genes associated with growth-factor signaling and atherosclerosis. Furthermore, JQ1 ameliorates AngII-induced hypertension, medial hypertrophy and inflammation in vivo in mice. These results demonstrate AngII-induced signals integrate enhancers/SEs and lncRNAs to increase expression of genes involved in VSMC dysfunction, and could uncover novel therapies.

Vertical sleeve gastrectomy reverses diet-induced gene-regulatory changes impacting lipid metabolism.

Vertical sleeve gastrectomy (VSG) produces sustainable weight loss, remission of type 2 diabetes (T2D), and improvement of nonalcoholic fatty liver disease (NAFLD). However, the molecular mechanisms underlying the metabolic benefits of VSG have remained elusive. According to our previous results, diet-induced obesity induces epigenetic modifications to chromatin in mouse liver. We demonstrate here that VSG in C57BL/6J wild-type male mice can reverse these chromatin modifications and thereby impact the expression of key metabolic genes. Genes involved in lipid metabolism, especially omega-6 fatty acid metabolism, are up-regulated in livers of mice after VSG while genes in inflammatory pathways are down-regulated after VSG. Consistent with gene expression changes, regulatory regions near genes involved in inflammatory response displayed decreased chromatin accessibility after VSG. Our results indicate that VSG induces global regulatory changes that impact hepatic inflammatory and lipid metabolic pathways, providing new insight into the mechanisms underlying the beneficial metabolic effects induced by VSG.

Transgenerational programming of longevity through E(z)-mediated histone H3K27 trimethylation in Drosophila.

Transgenerational effects on health and development of early-life nutrition have gained increased attention recently. However, the underlying mechanisms of transgenerational transmission are only starting to emerge, with epigenetics as perhaps the most important mechanism. We recently reported the first animal model to study transgenerational programming of longevity after early-life dietary manipulations, enabling investigations to identify underlying epigenetic mechanisms. We report here that post-eclosion dietary manipulation (PDM) with a low-protein (LP) diet upregulates the protein level of E(z), an H3K27 specific methyltransferase, leading to higher levels of H3K27 trimethylation (H3K27me3). This PDM-mediated change in H3K27me3 corresponded with a shortened longevity of F0 flies as well as their F2 offspring. Specific RNAi-mediated post-eclosion knockdown of E(z) or pharmacological inhibition of its enzymatic function with EPZ-6438 in the F0 parents improved longevity while rendering H3K27me3 low across generations. Importantly, addition of EPZ-6438 to the LP diet fully alleviated the longevity-reducing effect of the LP PDM, supporting the increased level of E(z)-dependent H3K27me3 as the primary cause and immediate early-life period as the critical time to program longevity through epigenetic regulation. These observations establish E(z)-mediated H3K27me3 as one epigenetic mechanism underlying nutritional programming of longevity and support the use of EPZ-6438 to extend lifespan.

Regional glutamine deficiency in tumours promotes dedifferentiation through inhibition of histone demethylation.

Poorly organized tumour vasculature often results in areas of limited nutrient supply and hypoxia. Despite our understanding of solid tumour responses to hypoxia, how nutrient deprivation regionally affects tumour growth and therapeutic response is poorly understood. Here, we show that the core region of solid tumours displayed glutamine deficiency compared with other amino acids. Low glutamine in tumour core regions led to dramatic histone hypermethylation due to decreased α-ketoglutarate levels, a key cofactor for the Jumonji-domain-containing histone demethylases. Using patient-derived (V600E)BRAF melanoma cells, we found that low-glutamine-induced histone hypermethylation resulted in cancer cell dedifferentiation and resistance to BRAF inhibitor treatment, which was largely mediated by methylation on H3K27, as knockdown of the H3K27-specific demethylase KDM6B and the methyltransferase EZH2 respectively reproduced and attenuated the low-glutamine effects in vitro and in vivo. Thus, intratumoral regional variation in the nutritional microenvironment contributes to tumour heterogeneity and therapeutic response.

Chromatin variation associated with liver metabolism is mediated by transposable elements.

BACKGROUND: Functional regulatory regions in eukaryotic genomes are characterized by the disruption of nucleosomes leading to accessible chromatin. The modulation of chromatin accessibility is one of the key mediators of transcriptional regulation, and variation in chromatin accessibility across individuals has been linked to complex traits and disease susceptibility. While mechanisms responsible for chromatin variation across individuals have been investigated, the overwhelming majority of chromatin variation remains unexplained. Furthermore, the processes through which the variation of chromatin accessibility contributes to phenotypic diversity remain poorly understood. RESULTS: We profiled chromatin accessibility in liver from seven strains of mice with phenotypic diversity in response to a high-fat/high-sucrose (HF/HS) diet and identified reproducible chromatin variation across the individuals. We found that sites of variable chromatin accessibility were more likely to coincide with particular classes of transposable elements (TEs) than sites with common chromatin signatures. Evolutionarily younger long interspersed nuclear elements (LINEs) are particularly likely to harbor variable chromatin sites. These younger LINEs are enriched for binding sites of immune-associated transcription factors, whereas older LINEs are enriched for liver-specific transcription factors. Genomic region enrichment analysis indicates that variable chromatin sites at TEs may function to regulate liver metabolic pathways. CRISPR-Cas9 deletion of a number of variable chromatin sites at TEs altered expression of nearby metabolic genes. Finally, we show that polymorphism of TEs and differential DNA methylation at TEs can both influence chromatin variation. CONCLUSIONS: Our results demonstrate that specific classes of TEs show variable chromatin accessibility across strains of mice that display phenotypic diversity in response to a HF/HS diet. These results indicate that chromatin variation at TEs is an important contributor to phenotypic variation among populations.