RNA-DNA triplexes: molecular mechanisms and functional relevance.

Interactions of RNA with DNA are principles of gene expression control that have recently gained considerable attention. Among RNA-DNA interactions are R-loops and RNA-DNA hybrid G-quadruplexes, as well as RNA-DNA triplexes. It is proposed that RNA-DNA triplexes guide RNA-associated regulatory proteins to specific genomic locations, influencing transcription and epigenetic decision making. Although triplex formation initially was considered solely an in vitro event, recent progress in computational, biochemical, and biophysical methods support in vivo functionality with relevance for gene expression control. Here, we review the central methodology and biology of triplexes, outline paradigms required for triplex function, and provide examples of physiologically important triplex-forming long non-coding RNAs.

NCoR1 limits angiogenic capacity by altering Notch signaling.

Corepressors negatively regulate gene expression by chromatin compaction. Targeted regulation of gene expression could provide a means to control endothelial cell phenotype. We hypothesize that by targeting corepressor proteins, endothelial angiogenic function can be improved. To study this, the expression and function of nuclear corepressors in human umbilical vein endothelial cells (HUVEC) and in murine organ culture was studied. RNA-seq revealed that nuclear receptor corepressor 1 (NCoR1), silencing mediator of retinoid and thyroid hormone receptors (SMRT) and repressor element-1 silencing transcription factor (REST) are the highest expressed corepressors in HUVECs. Knockout and knockdown strategies demonstrated that the depletion of NCoR1 increased the angiogenic capacity of endothelial cells, whereas depletion of SMRT or REST did not. Interestingly, the effect was VEGF signaling independent. NCoR1 depletion significantly upregulated angiogenesis-associated genes, especially tip cell genes, including ESM1, DLL4 and NOTCH4, as observed by RNA- and ATAC-seq. Confrontation assays comparing cells with and without NCoR1-deficiency revealed that loss of NCoR1 promotes a tip-cell position during spheroid sprouting. Moreover, a proximity ligation assay identified NCoR1 as a direct binding partner of the Notch-signaling-related transcription factor RBPJk. Luciferase assays showed that siRNA-mediated knockdown of NCOR1 promotes RBPJk activity. Furthermore, NCoR1 depletion prompts upregulation of several elements in the Notch signaling cascade. Downregulation of NOTCH4, but not NOTCH1, prevented the positive effect of NCOR1 knockdown on spheroid outgrowth. Collectively, these data indicate that decreasing NCOR1 expression is an attractive approach to promote angiogenic function.

Acute injury to the mouse carotid artery provokes a distinct healing response.

Treatment of vascular stenosis with angioplasty results in acute vascular damage, which may lead to restenosis. Owing to the highly complex cellularity of blood vessels, the healing response following this damage is incompletely understood. To gain further insight into this process, scRNA-seq of mouse carotid tissue after wire injury was performed. Stages of acute inflammation, resolution and remodeling were recapitulated in these data. To identify cell types which give rise to neointima, analyses focused on smooth muscle cell and fibroblast populations, and included data integration with scRNA-seq data from myocardial infarction and atherosclerosis datasets. Following carotid injury, a subpopulation of smooth muscle cells which also arises during atherosclerosis and myocardial infarction was identified. So-called stem cell/endothelial cell/monocyte (SEM) cells are candidates for repopulating injured vessels, and were amongst the most proliferative cell clusters following wire-injury of the carotid artery. Importantly, SEM cells exhibit specific transcriptional profiles which could be therapeutically targeted. SEM cell gene expression patterns could also be detected in bulk RNA-sequencing of neointimal tissue isolated from injured carotid vessels by laser capture microdissection. These data indicate that phenotypic plasticity of smooth muscle cells is highly important to the progression of lumen loss following acute carotid injury. Interference with SEM cell formation could be an innovative approach to combat development of restenosis.

Computational Methods to Study DNA:DNA:RNA Triplex Formation by lncRNAs.

Long non-coding RNAs (lncRNAs) impact cell function via numerous mechanisms. In the nucleus, interactions between lncRNAs and DNA and the consequent formation of non-canonical nucleic acid structures seems to be particularly relevant. Along with interactions between single-stranded RNA (ssRNA) and single-stranded DNA (ssDNA), such as R-loops, ssRNA can also interact with double-stranded DNA (dsDNA) to form DNA:DNA:RNA triplexes. A major challenge in the study of DNA:DNA:RNA triplexes is the identification of the precise RNA component interacting with specific regions of the dsDNA. As this is a crucial step towards understanding lncRNA function, there exist several computational methods designed to predict these sequences. This review summarises the recent progress in the prediction of triplex formation and highlights important DNA:DNA:RNA triplexes. In particular, different prediction tools (Triplexator, LongTarget, TRIPLEXES, Triplex Domain Finder, TriplexFFP, TriplexAligner and Fasim-LongTarget) will be discussed and their use exemplified by selected lncRNAs, whose DNA:DNA:RNA triplex forming potential was validated experimentally. Collectively, these tools revealed that DNA:DNA:RNA triplexes are likely to be numerous and make important contributions to gene expression regulation.

The endothelial-enriched lncRNA LINC00607 mediates angiogenic function.

Long non-coding RNAs (lncRNAs) can act as regulatory RNAs which, by altering the expression of target genes, impact on the cellular phenotype and cardiovascular disease development. Endothelial lncRNAs and their vascular functions are largely undefined. Deep RNA-Seq and FANTOM5 CAGE analysis revealed the lncRNA LINC00607 to be highly enriched in human endothelial cells. LINC00607 was induced in response to hypoxia, arteriosclerosis regression in non-human primates, post-atherosclerotic cultured endothelial cells from patients and also in response to propranolol used to induce regression of human arteriovenous malformations. siRNA knockdown or CRISPR/Cas9 knockout of LINC00607 attenuated VEGF-A-induced angiogenic sprouting. LINC00607 knockout in endothelial cells also integrated less into newly formed vascular networks in an in vivo assay in SCID mice. Overexpression of LINC00607 in CRISPR knockout cells restored normal endothelial function. RNA- and ATAC-Seq after LINC00607 knockout revealed changes in the transcription of endothelial gene sets linked to the endothelial phenotype and in chromatin accessibility around ERG-binding sites. Mechanistically, LINC00607 interacted with the SWI/SNF chromatin remodeling protein BRG1. CRISPR/Cas9-mediated knockout of BRG1 in HUVEC followed by CUT&RUN revealed that BRG1 is required to secure a stable chromatin state, mainly on ERG-binding sites. In conclusion, LINC00607 is an endothelial-enriched lncRNA that maintains ERG target gene transcription by interacting with the chromatin remodeler BRG1 to ultimately mediate angiogenesis.

Endothelial deletion of the cytochrome P450 reductase leads to cardiac remodelling.

The cytochrome P450 reductase (POR) transfers electrons to all microsomal cytochrome P450 enzymes (CYP450) thereby driving their activity. In the vascular system, the POR/CYP450 system has been linked to the production of epoxyeicosatrienoic acids (EETs) but also to the generation of reactive oxygen species. In cardiac myocytes (CMs), EETs have been shown to modulate the cardiac function and have cardioprotective effects. The functional importance of the endothelial POR/CYP450 system in the heart is unclear and was studied here using endothelial cell-specific, inducible knockout mice of POR (ecPOR-/-). RNA sequencing of murine cardiac cells revealed a cell type-specific expression of different CYP450 homologues. Cardiac endothelial cells mainly expressed members of the CYP2 family which produces EETs, and of the CYP4 family that generates omega fatty acids. Tamoxifen-induced endothelial deletion of POR in mice led to cardiac remodelling under basal conditions, as shown by an increase in heart weight to body weight ratio and an increased CM area as compared to control animals. Endothelial deletion of POR was associated with a significant increase in endothelial genes linked to protein synthesis with no changes in genes of the oxidative stress response. CM of ecPOR-/- mice exhibited attenuated expression of genes linked to mitochondrial function and an increase in genes related to cardiac myocyte contractility. In a model of pressure overload (transverse aortic constriction, TAC with O-rings), ecPOR-/- mice exhibited an accelerated reduction in cardiac output (CO) and stroke volume (SV) as compared to control mice. These results suggest that loss of endothelial POR along with a reduction in EETs leads to an increase in vascular stiffness and loss in cardioprotection, resulting in cardiac remodelling.

Long non-coding RNA PCAT19 safeguards DNA in quiescent endothelial cells by preventing uncontrolled phosphorylation of RPA2.

In healthy vessels, endothelial cells maintain a stable, differentiated, and growth-arrested phenotype for years. Upon injury, a rapid phenotypic switch facilitates proliferation to restore tissue perfusion. Here we report the identification of the endothelial cell-enriched long non-coding RNA (lncRNA) PCAT19, which contributes to the proliferative switch and acts as a safeguard for the endothelial genome. PCAT19 is enriched in confluent, quiescent endothelial cells and binds to the full replication protein A (RPA) complex in a DNA damage- and cell-cycle-related manner. Our results suggest that PCAT19 limits the phosphorylation of RPA2, primarily on the serine 33 (S33) residue, and thereby facilitates an appropriate DNA damage response while slowing cell cycle progression. Reduction in PCAT19 levels in response to either loss of cell contacts or knockdown promotes endothelial proliferation and angiogenesis. Collectively, PCAT19 acts as a dynamic guardian of the endothelial genome and facilitates rapid switching from quiescence to proliferation.

HIF1α-AS1 is a DNA:DNA:RNA triplex-forming lncRNA interacting with the HUSH complex.

DNA:DNA:RNA triplexes that are formed through Hoogsteen base-pairing of the RNA in the major groove of the DNA duplex have been observed in vitro, but the extent to which these interactions occur in cells and how they impact cellular functions remains elusive. Using a combination of bioinformatic techniques, RNA/DNA pulldown and biophysical studies, we set out to identify functionally important DNA:DNA:RNA triplex-forming long non-coding RNAs (lncRNA) in human endothelial cells. The lncRNA HIF1α-AS1 was retrieved as a top hit. Endogenous HIF1α-AS1 reduces the expression of numerous genes, including EPH Receptor A2 and Adrenomedullin through DNA:DNA:RNA triplex formation by acting as an adapter for the repressive human silencing hub complex (HUSH). Moreover, the oxygen-sensitive HIF1α-AS1 is down-regulated in pulmonary hypertension and loss-of-function approaches not only result in gene de-repression but also enhance angiogenic capacity. As exemplified here with HIF1α-AS1, DNA:DNA:RNA triplex formation is a functionally important mechanism of trans-acting gene expression control.

A universal model of RNA.DNA:DNA triplex formation accurately predicts genome-wide RNA-DNA interactions.

RNA.DNA:DNA triple helix (triplex) formation is a form of RNA-DNA interaction which regulates gene expression but is difficult to study experimentally in vivo. This makes accurate computational prediction of such interactions highly important in the field of RNA research. Current predictive methods use canonical Hoogsteen base pairing rules, which whilst biophysically valid, may not reflect the plastic nature of cell biology. Here, we present the first optimization approach to learn a probabilistic model describing RNA-DNA interactions directly from motifs derived from triplex sequencing data. We find that there are several stable interaction codes, including Hoogsteen base pairing and novel RNA-DNA base pairings, which agree with in vitro measurements. We implemented these findings in TriplexAligner, a program that uses the determined interaction codes to predict triplex binding. TriplexAligner predicts RNA-DNA interactions identified in all-to-all sequencing data more accurately than all previously published tools in human and mouse and also predicts previously studied triplex interactions with known regulatory functions. We further validated a novel triplex interaction using biophysical experiments. Our work is an important step towards better understanding of triplex formation and allows genome-wide analyses of RNA-DNA interactions.

Epoxyeicosatrienoic Acid and Prostanoid Crosstalk at the Receptor and Intracellular Signaling Levels to Maintain Vascular Tone.

Epoxyeicosatrienoic acids (EETs) are signaling lipids produced by the cytochrome P450-(CYP450)-mediated epoxygenation of arachidonic acid. EETs have numerous biological effects on the vascular system, but aspects including their species specificity make their effects on vascular tone controversial. CYP450 enzymes require the 450-reductase (POR) for their activity. We set out to determine the contribution of endothelial CYP450 to murine vascular function using isolated aortic ring preparations from tamoxifen-inducible endothelial cell-specific POR knockout mice (ecPOR-/-). Constrictor responses to phenylephrine were similar between control (CTR) and ecPOR-/- mice. Contrastingly, sensitivity to the thromboxane receptor agonist U46619 and prostaglandin E2 (PGE2) was increased following the deletion of POR. Ex vivo incubation with a non-hydrolyzable EET (14,15-EE-8(Z)-E, EEZE) reversed the increased sensitivity to U46619 to the levels of CTR. EETs had no effect on vascular tone in phenylephrine-preconstricted vessels, but dilated vessels contracted with U46619 or PGE2. As U46619 acts through RhoA-dependent kinase, this system was analyzed. The deletion of POR affected the expression of genes in this pathway and the inhibition of Rho-GTPase with SAR407899 decreased sensitivity to U46619. These data suggest that EET and prostanoid crosstalk at the receptor level and that lack of EET production sensitizes vessels to vasoconstriction via the induction of the Rho kinase system.

Loss of Endothelial Cytochrome P450 Reductase Induces Vascular Dysfunction in Mice.

BACKGROUND: POR (cytochrome P450 reductase) provides electrons for the catalytic activity of the CYP (cytochrome P450) monooxygenases. CYPs are dual-function enzymes as they generate protective vasoactive mediators derived from polyunsaturated fatty acids but also reactive oxygen species. It is not known in which conditions the endothelial POR/CYP system is beneficial versus deleterious. Here, the activity of all CYP enzymes was eliminated in the vascular endothelium to examine its impact on vascular function. METHODS: An endothelial-specific, tamoxifen-inducible POR knockout mouse (ecPOR-/-) was generated. Vascular function was studied by organ chamber experiments. eNOS (endothelial nitric oxide synthase) activity was accessed by heavy arginine/citrulline LC-MS/MS detection and phosphorylation of serine1177 in aortic rings. CYP-derived epoxyeicosatrienoic acids and prostanoids were measured by LC-MS/MS. Gene expression of aorta and endothelial cells was profiled by RNA sequencing. Blood pressure was measured by telemetry. RESULTS: Acetylcholine-induced endothelium-dependent relaxation was attenuated in isolated vessels of ecPOR-/- as compared with control mice. Additionally, ecPOR-/- mice had attenuated eNOS activity and eNOS/AKT phosphorylation. POR deletion reduced endothelial stores of CYP-derived epoxyeicosatrienoic acids but increased vascular prostanoids. This phenomenon was paralleled by the induction of genes implicated in eicosanoid generation. In response to Ang II (angiotensin II) infusion, blood pressure increased significantly more in ecPOR-/- mice. Importantly, the cyclooxygenase inhibitor Naproxen selectively lowered the Ang II-induced hypertension in ecPOR-/- mice. CONCLUSIONS: POR expression in endothelial cells maintains eNOS activity and its loss results in an overactivation of the vasoconstrictor prostanoid system. Through these mechanisms, loss of endothelial POR induces vascular dysfunction and hypertension.

Nuclear receptor activation shapes spatial genome organization essential for gene expression control: lessons learned from the vitamin D receptor.

Spatial genome organization is tightly controlled by several regulatory mechanisms and is essential for gene expression control. Nuclear receptors are ligand-activated transcription factors that modulate physiological and pathophysiological processes and are primary pharmacological targets. DNA binding of the important loop-forming insulator protein CCCTC-binding factor (CTCF) was modulated by 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3). We performed CTCF HiChIP assays to produce the first genome-wide dataset of CTCF long-range interactions in 1,25(OH)2D3-treated cells, and to determine whether dynamic changes of spatial chromatin interactions are essential for fine-tuning of nuclear receptor signaling. We detected changes in 3D chromatin organization upon vitamin D receptor (VDR) activation at 3.1% of all observed CTCF interactions. VDR binding was enriched at both differential loop anchors and within differential loops. Differential loops were observed in several putative functional roles including TAD border formation, promoter-enhancer looping, and establishment of VDR-responsive insulated neighborhoods. Vitamin D target genes were enriched in differential loops and at their anchors. Secondary vitamin D effects related to dynamic chromatin domain changes were linked to location of downstream transcription factors in differential loops. CRISPR interference and loop anchor deletion experiments confirmed the functional relevance of nuclear receptor ligand-induced adjustments of the chromatin 3D structure for gene expression regulation.

DNA topoisomerase inhibition with the HIF inhibitor acriflavine promotes transcription of lncRNAs in endothelial cells.

The transcription factor hypoxia-inducible factor 1 (HIF1) is an important driver of cancer and is therefore an attractive drug target. Acriflavine (ACF) has been suggested to inhibit HIF1, but its mechanism of action is unknown. Here we investigated the interaction of ACF with DNA and long non-coding RNAs (lncRNAs) and its function in human endothelial cells. ACF promoted apoptosis and reduced proliferation, network formation, and angiogenic capacity. It also induced changes in gene expression, as determined by RNA sequencing (RNA-seq), which could not be attributed to specific inhibition of HIF1. A similar response was observed in murine lung endothelial cells. Although ACF increased and decreased a similar number of protein-coding genes, lncRNAs were preferentially upregulated under normoxic and hypoxic conditions. An assay for transposase accessibility with subsequent DNA sequencing (ATAC-seq) demonstrated that ACF induced strong changes in chromatin accessibility at lncRNA promoters. Immunofluorescence showed displacement of DNA:RNA hybrids. Such effects might be due to ACF-mediated topoisomerase inhibition, which was indeed the case, as reflected by DNA unwinding assays. Comparison with other acridine derivatives and topoisomerase inhibitors suggested that the specific function of ACF is an effect of acridinium-class compounds. This study demonstrates that ACF inhibits topoisomerases rather than HIF specifically and that it elicits a unique expression response of lncRNAs.

Reactive Oxygen Species Differentially Modulate the Metabolic and Transcriptomic Response of Endothelial Cells.

Reactive oxygen species (ROS) are important mediators of both physiological and pathophysiological signal transduction in the cardiovascular system. The effects of ROS on cellular processes depend on the concentration, localization, and duration of exposure. Cellular stress response mechanisms have evolved to mitigate the negative effects of acute oxidative stress. In this study, we investigate the short-term and long-term metabolic and transcriptomic response of human umbilical vein endothelial cells (HUVEC) to different types and concentrations of ROS. To generate intracellular H2O2, we utilized a lentiviral chemogenetic approach for overexpression of human D-amino acid oxidase (DAO). DAO converts D-amino acids into their corresponding imino acids and H2O2. HUVEC stably overexpressing DAO (DAO-HUVEC) were exposed to D-alanine (3 mM), exogenous H2O2 (10 µM or 300 µM), or menadione (5 µM) for various timepoints and subjected to global untargeted metabolomics (LC-MS/MS) and RNAseq by MACE (Massive analysis of cDNA ends). A total of 300 µM H2O2 led to pronounced changes on both the metabolic and transcriptomic level. In particular, metabolites linked to redox homeostasis, energy-generating pathways, and nucleotide metabolism were significantly altered. Furthermore, 300 µM H2O2 affected genes related to the p53 pathway and cell cycle. In comparison, the effects of menadione and DAO-derived H2O2 mainly occurred at gene expression level. Collectively, all types of ROS led to subtle changes in the expression of ribosomal genes. Our results show that different types and concentration of ROS lead to a different metabolic and transcriptomic response in endothelial cells.

The hydrogen-peroxide producing NADPH oxidase 4 does not limit neointima development after vascular injury in mice.

OBJECTIVE: The NADPH oxidase Nox4 is an important source of H2O2. Nox4-derived H2O2 limits vascular inflammation and promotes smooth muscle differentiation. On this basis, the role of Nox4 for restenosis development was determined in the mouse carotid artery injury model. METHODS AND RESULTS: Genetic deletion of Nox4 by a tamoxifen-activated Cre-Lox-system did not impact on neointima formation in the carotid artery wire injury model. To understand this unexpected finding, time-resolved single-cell RNA-sequencing (scRNAseq) from injured carotid arteries of control mice and massive-analysis-of-cDNA-ends (MACE)-RNAseq from the neointima harvested by laser capture microdissection of control and Nox4 knockout mice was performed. This revealed that resting smooth muscle cells (SMCs) and fibroblasts exhibit high Nox4 expression, but that the proliferating de-differentiated SMCs, which give rise to the neointima, have low Nox4 expression. In line with this, the first weeks after injury, gene expression was unchanged between the carotid artery neointimas of control and Nox4 knockout mice. CONCLUSION: Upon vascular injury, Nox4 expression is transiently lost in the cells which comprise the neointima. NADPH oxidase 4 therefore does not interfere with restenosis development after wire-induced vascular injury.

A hierarchical regulatory network analysis of the vitamin D induced transcriptome reveals novel regulators and complete VDR dependency in monocytes.

The transcription factor vitamin D receptor (VDR) is the high affinity nuclear target of the biologically active form of vitamin D3 (1,25(OH)2D3). In order to identify pure genomic transcriptional effects of 1,25(OH)2D3, we used VDR cistrome, transcriptome and open chromatin data, obtained from the human monocytic cell line THP-1, for a novel hierarchical analysis applying three bioinformatics approaches. We predicted 75.6% of all early 1,25(OH)2D3-responding (2.5 or 4 h) and 57.4% of the late differentially expressed genes (24 h) to be primary VDR target genes. VDR knockout led to a complete loss of 1,25(OH)2D3-induced genome-wide gene regulation. Thus, there was no indication of any VDR-independent non-genomic actions of 1,25(OH)2D3 modulating its transcriptional response. Among the predicted primary VDR target genes, 47 were coding for transcription factors and thus may mediate secondary 1,25(OH)2D3 responses. CEBPA and ETS1 ChIP-seq data and RNA-seq following CEBPA knockdown were used to validate the predicted regulation of secondary vitamin D target genes by both transcription factors. In conclusion, a directional network containing 47 partly novel primary VDR target transcription factors describes secondary responses in a highly complex vitamin D signaling cascade. The central transcription factor VDR is indispensable for all transcriptome-wide effects of the nuclear hormone.

ZNF354C is a transcriptional repressor that inhibits endothelial angiogenic sprouting.

Zinc finger proteins (ZNF) are a large group of transcription factors with diverse functions. We recently discovered that endothelial cells harbour a specific mechanism to limit the action of ZNF354C, whose function in endothelial cells is unknown. Given that ZNF354C has so far only been studied in bone and tumour, its function was determined in endothelial cells. ZNF354C is expressed in vascular cells and localises to the nucleus and cytoplasm. Overexpression of ZNF354C in human endothelial cells results in a marked inhibition of endothelial sprouting. RNA-sequencing of human microvascular endothelial cells with and without overexpression of ZNF354C revealed that the protein is a potent transcriptional repressor. ZNF354C contains an active KRAB domain which mediates this suppression as shown by mutagenesis analysis. ZNF354C interacts with dsDNA, TRIM28 and histones, as observed by proximity ligation and immunoprecipitation. Moreover, chromatin immunoprecipitation revealed that the ZNF binds to specific endothelial-relevant target-gene promoters. ZNF354C suppresses these genes as shown by CRISPR/Cas knockout and RNAi. Inhibition of endothelial sprouting by ZNF354C is dependent on the amino acids DV and MLE of the KRAB domain. These results demonstrate that ZNF354C is a repressive transcription factor which acts through a KRAB domain to inhibit endothelial angiogenic sprouting.

The endocannabinoid anandamide has an anti-inflammatory effect on CCL2 expression in vascular smooth muscle cells.

Endocannabinoids are important lipid-signaling mediators. Both protective and deleterious effects of endocannabinoids in the cardiovascular system have been reported but the mechanistic basis for these contradicting observations is unclear. We set out to identify anti-inflammatory mechanisms of endocannabinoids in the murine aorta and in human vascular smooth muscle cells (hVSMC). In response to combined stimulation with cytokines, IL-1ß and TNFα, the murine aorta released several endocannabinoids, with anandamide (AEA) levels being the most significantly increased. AEA pretreatment had profound effects on cytokine-induced gene expression in hVSMC and murine aorta. As revealed by RNA-Seq analysis, the induction of a subset of 21 inflammatory target genes, including the important cytokine CCL2 was blocked by AEA. This effect was not mediated through AEA-dependent interference of the AP-1 or NF-κB pathways but rather through an epigenetic mechanism. In the presence of AEA, ATAC-Seq analysis and chromatin-immunoprecipitations revealed that CCL2 induction was blocked due to increased levels of H3K27me3 and a decrease of H3K27ac leading to compacted chromatin structure in the CCL2 promoter. These effects were mediated by recruitment of HDAC4 and the nuclear corepressor NCoR1 to the CCL2 promoter. This study therefore establishes a novel anti-inflammatory mechanism for the endogenous endocannabinoid AEA in vascular smooth muscle cells. Furthermore, this work provides a link between endogenous endocannabinoid signaling and epigenetic regulation.