Gene regulation shapes the evolution of phenotypic diversity. We investigated the evolution of liver promoters and enhancers in six primate species using ChIP-seq (H3K27ac and H3K4me1) to profile cis-regulatory elements (CREs) and using RNAseq to characterize gene expression in the same individuals. To quantify regulatory divergence, we compared CRE activity across species by testing differential ChIP-seq read depths directly measured for orthologous sequences. We show that the primate regulatory landscape is largely conserved across the lineage, with 63% of the tested human liver CREs showing similar activity across species. Conserved CRE function is associated with sequence conservation, proximity to coding genes, cell-type specificity, and transcription factor binding. Newly evolved CREs are enriched in immune response and neurodevelopmental functions. We further demonstrate that conserved CREs bind master regulators, suggesting that while CREs contribute to species adaptation to the environment, core functions remain intact. Newly evolved CREs are enriched in young transposable elements (TEs), including Long-Terminal-Repeats (LTRs) and SINE-VNTR-Alus (SVAs), that significantly affect gene expression. Conversely, only 16% of conserved CREs overlap TEs. We tested the cis-regulatory activity of 69 TE subfamilies by luciferase reporter assays, spanning all major TE classes, and showed that 95.6% of tested TEs can function as either transcriptional activators or repressors. In conclusion, we demonstrated the critical role of TEs in primate gene regulation and illustrated potential mechanisms underlying evolutionary divergence among the primate species through the noncoding genome.

Pancreatic ductal adenocarcinoma (PDA) has a poor prognosis, and new strategies for prevention and treatment are urgently needed. We previously reported that histone H4 acetylation is elevated in pancreatic acinar cells harboring Kras mutations prior to the appearance of premalignant lesions. Because acetyl-CoA abundance regulates global histone acetylation, we hypothesized that altered acetyl-CoA metabolism might contribute to metabolic or epigenetic alterations that promote tumorigenesis. We found that acetyl-CoA abundance is elevated in KRAS -mutant acinar cells and that its use in the mevalonate pathway supports acinar-to-ductal metaplasia (ADM). Pancreas-specifi c loss of the acetyl-CoA-producing enzyme ATP-citrate lyase (ACLY) accordingly suppresses ADM and tumor formation. In PDA cells, growth factors promote AKT-ACLY signaling and histone acetylation, and both cell proliferation and tumor growth can be suppressed by concurrent BET inhibition and statin treatment. Thus, KRAS-driven metabolic alterations promote acinar cell plasticity and tumor development, and targeting acetyl-CoA-dependent processes exerts anticancer effects. SIGNIFICANCE :Pancreatic cancer is among the deadliest of human malignancies. We identify a key role for the metabolic enzyme ACLY, which produces acetyl-CoA, in pancreatic carcinogenesis. The data suggest that acetyl-CoA use for histone acetylation and in the mevalonate pathway facilitates cell plasticity and proliferation, suggesting potential to target these pathways. Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): A. Carrer, S. Trefely,

A molecular analysis was carried out on the European hermit beetles (the Osmoderma eremita species complex) to explore their genetic diversification and the robustness of previous morphologically based taxonomic arrangements. Complete sequences of mtDNA cytochrome C oxidase I gene were obtained from 26 individuals. Mean levels of interspecific sequence divergence ranged from 0.044 to 0.186. The results indicate a clear-cut distinction between two clades. The first one includes the W-European O. eremita Scopoli 1763; and the two Italian endemic taxa Osmoderma italicum Sparacio, 2000 and Osmoderma cristinae Sparacio, 1994; from southern peninsular Italy and Sicily, respectively. The second one includes the widespread E-European Osmoderma barnabita; and the southern Balcanic Osmoderma lassallei Baraud and Tauzin, 1991 from Greece and European Turkey. Within the two clades, molecular and morphological data well support a specific rank for O. lassallei and O. barnabita on one side, and for O. eremita and O. cristinae on the other side, while the taxonomic position of O. italicum, more closely related to O. eremita, is still uncertain. Current geographical distribution, interspecific genetic diversification, and very low levels of intraspecific genetic divergence in western European populations of O. eremita sensu stricto are hypothesized to be the result of multiple speciation events (mainly occurred in refugial forest areas of the Italian and Balkan peninsulas and Sicily before and during the Pleistocene glacial peaks), followed by fast post-glacial northward and westward expansion of some species.

BackgroundTransposable elements (TE) are an important source of evolutionary novelty in gene regulation. However, the mechanisms by which TEs contribute to gene expression are largely uncharacterized.ResultsHere, we leverage Roadmap and GTEx data to investigate the association of TEs with active and repressed chromatin in 24 tissues. We find 112 human TE families enriched in active regions of the genome across tissues. Short Interspersed Nuclear Elements (SINEs) and DNA transposons are the most frequently enriched classes, while Long Terminal Repeat Retrotransposons (LTRs) are often enriched in a tissue-specific manner. We report across-tissue variability in TE enrichment in active regions. Genes with consistent expression across tissues are less likely to be associated with TE insertions. TE presence in repressed regions similarly follows tissue-specific patterns. Moreover, different TE classes correlate with different repressive marks: LTRs and Long Interspersed Nuclear Elements (LINEs) are overrepresented in regions marked by H3K9me3, while the other TEs are more likely to overlap regions with H3K27me3. Young TEs are typically enriched in repressed regions and depleted in active regions. We detect multiple instances of TEs that are enriched in tissue-specific active regulatory regions. Such TEs contain binding sites for transcription factors that are master regulators for the given tissue. These TEs are enriched in intronic enhancers, and their tissue-specific enrichment correlates with tissue-specific variations in the expression of the nearest genes.ConclusionsWe provide an integrated overview of the contribution of TEs to human gene regulation. Expanding previous analyses, we demonstrate that TEs can potentially contribute to the turnover of regulatory sequences in a tissue-specific fashion.Electronic supplementary materialThe online version of this article (10.1186/s12864-018-4850-3) contains supplementary material, which is available to authorized users.

Summary AT-rich interactive domain-containing protein 1A and 1B (ARID1A, ARID1B) are mutually exclusive subunits of the chromatin remodeler SWI/SNF. ARID1A is the most frequently mutated chromatin regulator across all cancers, and ovarian clear cell carcinoma (OCCC) carries the highest prevalence of ARID1A mutations (~57%). Despite evidence implicating ARID1A in tumorigenesis, the mechanism remains elusive. Here, we demonstrate that in OCCC ARID1A binds active regulatory elements. Depletion of ARID1A represses RNA Polymerase II (RNAPII) transcription, but results in modest changes to accessibility. Specifically, pausing of RNAPII is severely impaired after loss of ARID1A. Compromised pausing results in transcriptional dysregulation of active genes, which is compensated by upregulation of ARID1B. However, a subset of ARID1A-dependent genes is not rescued by ARID1B, including many p53 and estrogen receptor (ESR1) targets. Our results provide new insight on ARID1A-mediated tumorigenesis and unveil new functions of SWI/SNF in modulating RNAPII dynamics.

The control of cell fate is an epigenetic process initiated by transcription factors (TFs) that recognize DNA motifs and recruit activator complexes and transcriptional machineries to chromatin. Lineage specificity is thought to be provided solely by TF-motif pairing, while the recruited activators are passive. Here, we show that INTS13, a subunit of the Integrator complex, operates as monocytic/macrophagic differentiation factor. Integrator is a general activator of transcription at coding genes and is required for eRNA maturation. Here, we show that INTS13 functions as an independent sub-module and targets enhancers through Early Growth Response (EGR1/2) TFs and their co-factor NAB2. INTS13 binds poised monocytic enhancers eliciting chromatin looping and activation. Independent depletion of INTS13, EGR1, or NAB2 impairs monocytic differentiation of cell lines and primary human progenitors. Our data demonstrate that Integrator is not functionally homogeneous and has TF-specific regulatory potential, revealing a new enhancer regulatory axis that controls myeloid differentiation.

Results of a molecular analysis on the European hermit beetles (the <em>Osmoderma</em> <em>eremita</em> species-complex), recently published in a companion paper, are shortly dis- cussed and commented. These results indicate a clear-cut distinction between two clades. The first one includes the W-European <em>O. eremita</em> Scopoli, 1763, and the two Italian endemic taxa <em>O. italicum</em> Sparacio, 2000 and <em>O. cristinae</em> Sparacio, 1994, from southern peninsular Italy and Sicily, respectively. The second one includes the widespread E-European<em> O. barnabita</em> Motschulsky, 1845 (<strong>nom. resurr</strong>.), and the southern Balcanic <em>O. lassallei</em> Baraud &amp; Tauzin, 1991 from Greece and European Turkey. Within the two clades, molecular data well support a specific rank for <em>O. lassallei</em> and <em>O. barn­abita</em> on one side, and of <em>O.</em> <em>eremita</em> and <em>O. cristinae</em> on the other side, while the taxonomic position of <em>O. italicum</em>, more closely related with <em>O. eremita</em>, is still uncertain, waiting for analysis of additional specimens of this very rare taxon. Current geographical distribution, interspecific genetic diversification, and relatively low levels of intraspecific genetic divergence in <em>O. eremita</em> sensu stricto, are hypothesized to be the result of multiple speciation events (mainly occurred in refugial forest areas of the Italian and Balkan peninsulas and Sicily before and during the Pleistocene glacial peaks), followed by fast post-glacial northward and westward expansion of some species. The need of further genetic data on the rare and threatened hermit beetle species and the importance of more detailed information on their distribution ranges are emphasized and discussed, in order to plan conservation strategies in the near future. An updated worldwide checklist of the species of the genus <em>Osmoderma</em> is finally presented.

Monocytes and monocyte-derived macrophages originate through a multistep differentiation process. First, hematopoietic stem cells generate lineage-restricted progenitors that eventually develop into peripheral, postmitotic monocytes. Second, blood-circulating monocytes undergo differentiation into macrophages, which are specialized phagocytic cells capable of tissue infiltration. While monocytes mediate some level of inflammation and cell toxicity, macrophages boast the widest set of defense mechanisms against pathogens and elicit robust inflammatory responses. Here, we analyze the molecular determinants of monocytic and macrophagic commitment by profiling the EGR1 transcription factor. EGR1 is essential for monopoiesis and binds enhancers that regulate monocytic developmental genes such as CSF1R. However, differentiating macrophages present a very different EGR1 binding pattern. We identify novel binding sites of EGR1 at a large set of inflammatory enhancers, even in the absence of its binding motif. We show that EGR1 repressive activity results in suppression of inflammatory genes and is mediated by the NuRD corepressor complex.