Differential analysis of chromatin accessibility and histone modifications for predicting mouse developmental enhancers

Atypical Modes of CTCF Binding Facilitate Tissue-Specific and Neuronal Activity-Dependent Gene Expression States

Multivalent binding of CTCF to variable DNA sequences is thought to underlie its ability to mediate diverse cellular functions. CTCF typically binds a 20 base-pair consensus DNA sequence, but the full diversity of CTCF binding sites (CBS) within the genome has not been interrogated. We assessed CTCF occupancy in cultured cortical neurons and observed surprisingly that ~ 22% of CBS lack the consensus CTCF motif. We report here that sequence diversity at most of these atypical CBS involves degeneracy at specific nucleotide positions within the consensus CTCF motif, which likely affect the binding of CTCF zinc fingers 6 and 7. This mode of atypical CTCF binding defines most CBS at gene promoters, as well as CBS that are dynamically altered during neural differentiation and following neuronal stimulation, revealing how atypical CTCF binding could influence gene activity. Dynamic CBS are distributed both within and outside loop anchors and TAD boundaries, suggesting both looping-dependent and independent roles for CTCF. Finally, we describe a second mode of atypical CTCF binding to DNA sequences that are completely unrelated to the consensus CTCF motif, which are enriched within the bodies of tissue-specific genes. These tissue-specific atypical CBS are also enriched in H3K27ac, which marks cis-regulatory elements within chromatin, including enhancers. Overall, these results indicate how atypical CBS could dynamically regulate gene activity patterns during differentiation, development, and in response to environmental cues.

The genetic regulatory architecture and epigenomic basis for age-related changes in rattlesnake venom

Developmental phenotypic changes can evolve under selection imposed by age- and size-related ecological differences. Many of these changes occur through programmed alterations to gene expression patterns, but the molecular mechanisms and gene-regulatory networks underlying these adaptive changes remain poorly understood. Many venomous snakes, including the eastern diamondback rattlesnake (Crotalus adamanteus), undergo correlated changes in diet and venom expression as snakes grow larger with age, providing models for identifying mechanisms of timed expression changes that underlie adaptive life history traits. By combining a highly contiguous, chromosome-level genome assembly with measures of expression, chromatin accessibility, and histone modifications, we identified cis-regulatory elements and trans-regulatory factors controlling venom ontogeny in the venom glands of C. adamanteus. Ontogenetic expression changes were significantly correlated with epigenomic changes within genes, immediately adjacent to genes (e.g., promoters), and more distant from genes (e.g., enhancers). We identified 37 candidate transcription factors (TFs), with the vast majority being up-regulated in adults. The ontogenetic change is largely driven by an increase in the expression of TFs associated with growth signaling, transcriptional activation, and circadian rhythm/biological timing systems in adults with corresponding epigenomic changes near the differentially expressed venom genes. However, both expression activation and repression contributed to the composition of both adult and juvenile venoms, demonstrating the complexity and potential evolvability of gene regulation for this trait. Overall, given that age-based trait variation is common across the tree of life, we provide a framework for understanding gene-regulatory-network-driven life-history evolution more broadly.

Identification of somatic mutation-driven enhancers and their clinical utility in breast cancer

Somatic mutations contribute to cancer development by altering the activity of enhancers. In the study, a total of 135 mutation-driven enhancers, which displayed significant chromatin accessibility changes, were identified as candidate risk factors for breast cancer (BRCA). Furthermore, we identified four mutation-driven enhancers as independent prognostic factors for BRCA subtypes. In Her2 subtype, enhancer G > C mutation was associated with poorer prognosis through influencing its potential target genes FBXW9, TRIR, and WDR83. We identified aminoglutethimide and quinpirole as candidate drugs targeting the mutated enhancer. In normal subtype, enhancer G > A mutation was associated with poorer prognosis through influencing its target genes ALOX15B, LINC00324, and MPDU1. We identified eight candidate drugs such as erastin, colforsin, and STOCK1N-35874 targeting the mutated enhancer. Our findings suggest that somatic mutations contribute to breast cancer subtype progression by altering enhancer activity, which could be potential candidates for cancer therapy.

ASH2L regulates postnatal neurogenesis through Onecut2-mediated inhibition of TGF-β signaling pathway

The ability of neural stem/progenitor cells (NSPCs) to proliferate and differentiate is required through different stages of neurogenesis. Disturbance in the regulation of neurogenesis causes many neurological diseases, such as intellectual disability, autism, and schizophrenia. However, the intrinsic mechanisms of this regulation in neurogenesis remain poorly understood. Here, we report that Ash2l (Absent, small or homeotic discs-like 2), one core component of a multimeric histone methyltransferase complex, is essential for NSPC fate determination during postnatal neurogenesis. Deletion of Ash2l in NSPCs impairs their capacity for proliferation and differentiation, leading to simplified dendritic arbors in adult-born hippocampal neurons and deficits in cognitive abilities. RNA sequencing data reveal that Ash2l primarily regulates cell fate specification and neuron commitment. Furthermore, we identified Onecut2, a major downstream target of ASH2L characterized by bivalent histone modifications, and demonstrated that constitutive expression of Onecut2 restores defective proliferation and differentiation of NSPCs in adult Ash2l-deficient mice. Importantly, we identified that Onecut2 modulates TGF-β signaling in NSPCs and that treatment with a TGF-β inhibitor rectifies the phenotype of Ash2l-deficient NSPCs. Collectively, our findings reveal the ASH2L-Onecut2-TGF-β signaling axis that mediates postnatal neurogenesis to maintain proper forebrain function.

Toward a comprehensive catalog of regulatory elements

Regulatory elements are the genomic regions that interact with transcription factors to control cell-type-specific gene expression in different cellular environments. A precise and complete catalog of functional elements encoded by the human genome is key to understanding mammalian gene regulation. Here, we review the current state of regulatory element annotation. We first provide an overview of assays for characterizing functional elements, including genome, epigenome, transcriptome, three-dimensional chromatin interaction, and functional validation assays. We then discuss computational methods for defining regulatory elements, including peak-calling and other statistical modeling methods. Finally, we introduce several high-quality lists of regulatory element annotations and suggest potential future directions.

Relating enhancer genetic variation across mammals to complex phenotypes using machine learning

Protein-coding differences between species often fail to explain phenotypic diversity, suggesting the involvement of genomic elements that regulate gene expression such as enhancers. Identifying associations between enhancers and phenotypes is challenging because enhancer activity can be tissue-dependent and functionally conserved despite low sequence conservation. We developed the Tissue-Aware Conservation Inference Toolkit (TACIT) to associate candidate enhancers with species' phenotypes using predictions from machine learning models trained on specific tissues. Applying TACIT to associate motor cortex and parvalbumin-positive interneuron enhancers with neurological phenotypes revealed dozens of enhancer-phenotype associations, including brain size-associated enhancers that interact with genes implicated in microcephaly or macrocephaly. TACIT provides a foundation for identifying enhancers associated with the evolution of any convergently evolved phenotype in any large group of species with aligned genomes.

Matching queried single-cell open-chromatin profiles to large pools of single-cell transcriptomes and epigenomes for reference supported analysis

The true benefits of large single-cell transcriptome and epigenome data sets can be realized only with the development of new approaches and search tools for annotating individual cells. Matching a single-cell epigenome profile to a large pool of reference cells remains a major challenge. Here, we present scEpiSearch, which enables searching, comparison, and independent classification of single-cell open-chromatin profiles against a large reference of single-cell expression and open-chromatin data sets. Across performance benchmarks, scEpiSearch outperformed multiple methods in accuracy of search and low-dimensional coembedding of single-cell profiles, irrespective of platforms and species. Here we also demonstrate the unconventional utilities of scEpiSearch by applying it on single-cell epigenome profiles of K562 cells and samples from patients with acute leukaemia to reveal different aspects of their heterogeneity, multipotent behavior, and dedifferentiated states. Applying scEpiSearch on our single-cell open-chromatin profiles from embryonic stem cells (ESCs), we identified ESC subpopulations with more activity and poising for endoplasmic reticulum stress and unfolded protein response. Thus, scEpiSearch solves the nontrivial problem of amalgamating information from a large pool of single cells to identify and study the regulatory states of cells using their single-cell epigenomes.

Hnf1b renal expression directed by a distal enhancer responsive to Pax8

Xenopus provides a simple and efficient model system to study nephrogenesis and explore the mechanisms causing renal developmental defects in human. Hnf1b (hepatocyte nuclear factor 1 homeobox b), a gene whose mutations are the most commonly identified genetic cause of developmental kidney disease, is required for the acquisition of a proximo-intermediate nephron segment in Xenopus as well as in mouse. Genetic networks involved in Hnf1b expression during kidney development remain poorly understood. We decided to explore the transcriptional regulation of Hnf1b in the developing Xenopus pronephros and mammalian renal cells. Using phylogenetic footprinting, we identified an evolutionary conserved sequence (CNS1) located several kilobases (kb) upstream the Hnf1b transcription start and harboring epigenomic marks characteristics of a distal enhancer in embryonic and adult renal cells in mammals. By means of functional expression assays in Xenopus and mammalian renal cell lines we showed that CNS1 displays enhancer activity in renal tissue. Using CRISPR/cas9 editing in Xenopus tropicalis, we demonstrated the in vivo functional relevance of CNS1 in driving hnf1b expression in the pronephros. We further showed the importance of Pax8-CNS1 interaction for CNS1 enhancer activity allowing us to conclude that Hnf1b is a direct target of Pax8. Our work identified for the first time a Hnf1b renal specific enhancer and may open important perspectives into the diagnosis for congenital kidney anomalies in human, as well as modeling HNF1B-related diseases.

Acetate supplementation restores cognitive deficits caused by ARID1A haploinsufficiency in excitatory neurons

Mutations in AT-rich interactive domain-containing protein 1A (ARID1A) cause Coffin-Siris syndrome (CSS), a rare genetic disorder that results in mild to severe intellectual disabilities. However, the biological role of ARID1A in the brain remains unclear. In this study, we report that the haploinsufficiency of ARID1A in excitatory neurons causes cognitive impairment and defects in hippocampal synaptic transmission and dendritic morphology in mice. Similarly, human embryonic stem cell-derived excitatory neurons with deleted ARID1A exhibit fewer dendritic branches and spines, and abnormal electrophysiological activity. Importantly, supplementation of acetate, an epigenetic metabolite, can ameliorate the morphological and electrophysiological deficits observed in mice with Arid1a haploinsufficiency, as well as in ARID1A-null human excitatory neurons. Mechanistically, transcriptomic and ChIP-seq analyses demonstrate that acetate supplementation can increase the levels of H3K27 acetylation at the promoters of key regulatory genes associated with neural development and synaptic transmission. Collectively, these findings support the essential roles of ARID1A in the excitatory neurons and cognition and suggest that acetate supplementation could be a potential therapeutic intervention for CSS.

Differential enrichment of H3K9me3 in intrahepatic cholangiocarcinoma

Background

Intrahepatic cholangiocarcinoma (ICC) is a malignant tumor, which poses a serious threat to human health. Histone 3 lysine 9 trimethylation (H3K9me3) is a post-translational modification involved in regulating a broad range of biological processes and has been considered as potential therapeutic target in types of cancer. However, there is limited research on investigating profiles of histone modification H3K9me3 in ICC patients.

Methods

In this study, we applied the ChIP-seq technique to investigate the effect of H3K9me3 on ICC. Anti-H3K9me3 antibody was used for ChIP-seq in ICC (RBE cell lines) and HIBEpic (normal cell lines). MACS2 (peak-calling tools) was then used to identify the peaks recorded in RBE and HIBEpic cell lines. Gene expression, mutation and clinical data were downloaded from TCGA and cBioPortal databases.

Results

H3K9me3 exhibited abnormal methylation and influenced the process of abnormal gene expression in patients suffering from ICC. The Wnt/β-Catenin signaling pathway (also known as simply the WNT signaling pathway) was enriched in H3K9me3-regulated genes.

Conclusions

We are the first to report that H3K9me3 may play an important role in the progression of ICC. It promotes the understanding of epigenetic molecular mechanisms for ICC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01338-1.

Ranking reprogramming factors for cell differentiation

Transcription factor over-expression is a proven method for reprogramming cells to a desired cell type for regenerative medicine and therapeutic discovery. However, a general method for the identification of reprogramming factors to create an arbitrary cell type is an open problem. Here we examine the success rate of methods and data for differentiation by testing the ability of nine computational methods (CellNet, GarNet, EBseq, AME, DREME, HOMER, KMAC, diffTF and DeepAccess) to discover and rank candidate factors for eight target cell types with known reprogramming solutions. We compare methods that use gene expression, biological networks and chromatin accessibility data, and comprehensively test parameter and preprocessing of input data to optimize performance. We find the best factor identification methods can identify an average of 50-60% of reprogramming factors within the top ten candidates, and methods that use chromatin accessibility perform the best. Among the chromatin accessibility methods, complex methods DeepAccess and diffTF have higher correlation with the ranked significance of transcription factor candidates within reprogramming protocols for differentiation. We provide evidence that AME and diffTF are optimal methods for transcription factor recovery that will allow for systematic prioritization of transcription factor candidates to aid in the design of new reprogramming protocols.

Investigating chromatin accessibility during development and differentiation by ATAC-sequencing to guide the identification of cis-regulatory elements

Mapping accessible chromatin across time scales can give insights into its dynamic nature, for example during cellular differentiation and tissue or organism development. Analysis of such data can be utilised to identify functional cis-regulatory elements (CRE) and transcription factor binding sites and, when combined with transcriptomics, can reveal gene regulatory networks (GRNs) of expressed genes. Chromatin accessibility mapping is a powerful approach and can be performed using ATAC-sequencing (ATAC-seq), whereby Tn5 transposase inserts sequencing adaptors into genomic DNA to identify differentially accessible regions of chromatin in different cell populations. It requires low sample input and can be performed and analysed relatively quickly compared with other methods. The data generated from ATAC-seq, along with other genomic approaches, can help uncover chromatin packaging and potential cis-regulatory elements that may be responsible for gene expression. Here, we describe the ATAC-seq approach and give examples from mainly vertebrate embryonic development, where such datasets have identified the highly dynamic nature of chromatin, with differing landscapes between cellular precursors for different lineages.

Analysis of long and short enhancers in melanoma cell states

Understanding how enhancers drive cell-type specificity and efficiently identifying them is essential for the development of innovative therapeutic strategies. In melanoma, the melanocytic (MEL) and the mesenchymal-like (MES) states present themselves with different responses to therapy, making the identification of specific enhancers highly relevant. Using massively parallel reporter assays (MPRAs) in a panel of patient-derived melanoma lines (MM lines), we set to identify and decipher melanoma enhancers by first focusing on regions with state-specific H3K27 acetylation close to differentially expressed genes. An in-depth evaluation of those regions was then pursued by investigating the activity of overlapping ATAC-seq peaks along with a full tiling of the acetylated regions with 190 bp sequences. Activity was observed in more than 60% of the selected regions, and we were able to precisely locate the active enhancers within ATAC-seq peaks. Comparison of sequence content with activity, using the deep learning model DeepMEL2, revealed that AP-1 alone is responsible for the MES enhancer activity. In contrast, SOX10 and MITF both influence MEL enhancer function with SOX10 being required to achieve high levels of activity. Overall, our MPRAs shed light on the relationship between long and short sequences in terms of their sequence content, enhancer activity, and specificity across melanoma cell states.

Coordinated Gene Expression and Chromatin Regulation during Hydra Head Regeneration

The cnidarian model organism Hydra has long been studied for its remarkable ability to regenerate its head, which is controlled by a head organizer located near the hypostome. The canonical Wnt pathway plays a central role in head organizer function during regeneration and during bud formation, which is the asexual mode of reproduction in Hydra. However, it is unclear how shared the developmental programs of head organizer genesis are in budding and regeneration. Time-series analysis of gene expression changes during head regeneration and budding revealed a set of 298 differentially expressed genes during the 48-h head regeneration and 72-h budding time courses. In order to understand the regulatory elements controlling Hydra head regeneration, we first identified 27,137 open-chromatin elements that are open in one or more sections of the organism body or regenerating tissue. We used histone modification ChIP-seq to identify 9,998 candidate proximal promoter and 3,018 candidate enhancer-like regions respectively. We show that a subset of these regulatory elements is dynamically remodeled during head regeneration and identify a set of transcription factor motifs that are enriched in the enhancer regions activated during head regeneration. Our results show that Hydra displays complex gene regulatory structures of developmentally dynamic enhancers, which suggests that the evolution of complex developmental enhancers predates the split of cnidarians and bilaterians.

Genomic Marks Associated with Chromatin Compartments in the CTCF, RNAPII Loop and Genomic Windows

The nature of genome organization into two basic structural compartments is as yet undiscovered. However, it has been indicated to be a mechanism of gene expression regulation. Using the classification approach, we ranked genomic marks that hint at compartmentalization. We considered a broad range of marks, including GC content, histone modifications, DNA binding proteins, open chromatin, transcription and genome regulatory segmentation in GM12878 cells. Genomic marks were defined over CTCF or RNAPII loops, which are basic elements of genome 3D structure, and over 100 kb genomic windows. Experiments were carried out to empirically assess the whole set of features, as well as the individual features in classification of loops/windows, into compartment A or B. Using Monte Carlo Feature Selection and Analysis of Variance, we constructed a ranking of feature importance for classification. The best simple indicator of compartmentalization is DNase-seq open chromatin measurement for CTCF loops, H3K4me1 for RNAPII loops and H3K79me2 for genomic windows. Among DNA binding proteins, this is RUNX3 transcription factor for loops and RNAPII for genomic windows. Chromatin state prediction methods that indicate active elements like promoters, enhancers or heterochromatin enhance the prediction of loop segregation into compartments. However, H3K9me3, H4K20me1, H3K27me3 histone modifications and GC content poorly indicate compartments.

Mitochondrial dysfunction-induced H3K27 hyperacetylation perturbs enhancers in Parkinson's disease

Mitochondrial dysfunction is a major pathophysiological contributor to the progression of Parkinson’s disease (PD); however, whether it contributes to epigenetic dysregulation remains unknown. Here, we show that both chemically and genetically driven mitochondrial dysfunctions share a common mechanism of epigenetic dysregulation. Under both scenarios, lysine 27 acetylation of likely variant H3.3 (H3.3K27ac) increased in dopaminergic neuronal models of PD, thereby opening that region to active enhancer activity via H3K27ac. These vulnerable epigenomic loci represent potential transcription factor motifs for PD pathogenesis. We further confirmed that mitochondrial dysfunction induces H3K27ac in ex vivo and in vivo (MitoPark) neurodegenerative models of PD. Notably, the significantly increased H3K27ac in postmortem PD brains highlights the clinical relevance to the human PD population. Our results reveal an exciting mitochondrial dysfunction-metabolism-H3K27ac-transcriptome axis for PD pathogenesis. Collectively, the mechanistic insights link mitochondrial dysfunction to epigenetic dysregulation in dopaminergic degeneration and offer potential new epigenetic intervention strategies for PD.

PPARγ-Induced Global H3K27 Acetylation Maintains Osteo/Cementogenic Abilities of Periodontal Ligament Fibroblasts

The periodontal ligament is a soft connective tissue embedded between the alveolar bone and cementum, the surface hard tissue of teeth. Periodontal ligament fibroblasts (PDLF) actively express osteo/cementogenic genes, which contribute to periodontal tissue homeostasis. However, the key factors maintaining the osteo/cementogenic abilities of PDLF remain unclear. We herein demonstrated that PPARγ was expressed by in vivo periodontal ligament tissue and its distribution pattern correlated with alkaline phosphate enzyme activity. The knockdown of PPARγ markedly reduced the osteo/cementogenic abilities of PDLF in vitro, whereas PPARγ agonists exerted the opposite effects. PPARγ was required to maintain the acetylation status of H3K9 and H3K27, active chromatin markers, and the supplementation of acetyl-CoA, a donor of histone acetylation, restored PPARγ knockdown-induced decreases in the osteo/cementogenic abilities of PDLF. An RNA-seq/ChIP-seq combined analysis identified four osteogenic transcripts, RUNX2, SULF2, RCAN2, and RGMA, in the PPARγ-dependent active chromatin region marked by H3K27ac. Furthermore, RUNX2-binding sites were selectively enriched in the PPARγ-dependent active chromatin region. Collectively, these results identified PPARγ as the key transcriptional factor maintaining the osteo/cementogenic abilities of PDLF and revealed that global H3K27ac modifications play a role in the comprehensive osteo/cementogenic transcriptional alterations mediated by PPARγ.

Environmental Enrichment Induces Epigenomic and Genome Organization Changes Relevant for Cognition

In early development, the environment triggers mnemonic epigenomic programs resulting in memory and learning experiences to confer cognitive phenotypes into adulthood. To uncover how environmental stimulation impacts the epigenome and genome organization, we used the paradigm of environmental enrichment (EE) in young mice constantly receiving novel stimulation. We profiled epigenome and chromatin architecture in whole cortex and sorted neurons by deep-sequencing techniques. Specifically, we studied chromatin accessibility, gene and protein regulation, and 3D genome conformation, combined with predicted enhancer and chromatin interactions. We identified increased chromatin accessibility, transcription factor binding including CTCF-mediated insulation, differential occupancy of H3K36me3 and H3K79me2, and changes in transcriptional programs required for neuronal development. EE stimuli led to local genome re-organization by inducing increased contacts between chromosomes 7 and 17 (inter-chromosomal). Our findings support the notion that EE-induced learning and memory processes are directly associated with the epigenome and genome organization.

Transcriptional profiling of microglia in the injured brain reveals distinct molecular features underlying neurodegeneration

Neurotrauma has been recognized as a risk factor for neurodegenerative diseases, and sex difference of the incidence and outcome of neurodegenerative diseases has long been recognized. Past studies suggest that microglia could play a versatile role in both health and disease. So far, the microglial mechanisms underlying neurodegeneration and potentially lead to sex-specific therapies are still very open. Here we applied whole transcriptome analysis of microglia acutely isolated at different timepoints after a cortical stab wound injury to gain insight into genes that might be dysregulated and transcriptionally different between males and females after cortical injury. We found that microglia displayed distinct temporal and sexual molecular signatures of transcriptome after cortical injury. Hypotheses and gene candidates that we presented in the present study could be worthy to be examined to explore the roles of microglia in neurotrauma and in sex-biased neurodegenerative diseases.

Integrated Transcriptome and Histone Modification Analysis Reveals NDV Infection Under Heat Stress Affects Bursa Development and Proliferation in Susceptible Chicken Line

Two environmental factors, Newcastle disease and heat stress, are concurrently negatively impacting poultry worldwide and warrant greater attention into developing genetic resistance within chickens. Using two genetically distinct and highly inbred layer lines, Fayoumi and Leghorn, we explored how different genetic backgrounds affect the bursal response to a treatment of simultaneous Newcastle disease virus (NDV) infection at 6 days postinfection (dpi) while under chronic heat stress. The bursa is a primary lymphoid organ within birds and is crucial for the development of B cells. We performed RNA-seq and ChIP-seq targeting histone modifications on bursa tissue. Differential gene expression revealed that Leghorn, compared to Fayoumi, had significant down-regulation in genes involved in cell proliferation, cell cycle, and cell division. Interestingly, we also found greater differences in histone modification levels in response to treatment in Leghorns than Fayoumis, and biological processes enriched in associated target genes of H3K27ac and H3K4me1 were similarly associated with cell cycle and receptor signaling of lymphocytes. Lastly, we found candidate variants between the two genetic lines within exons of differentially expressed genes and regulatory elements with differential histone modification enrichment between the lines, which provides a strong foundation for understanding the effects of genetic variation on NDV resistance under heat stress. This study provides further understanding of the cellular mechanisms affected by NDV infection under heat stress in chicken bursa and identified potential genes and regulatory regions that may be targets for developing genetic resistance within chickens.

Peax: Interactive Visual Pattern Search in Sequential Data Using Unsupervised Deep Representation Learning

We present Peax, a novel feature-based technique for interactive visual pattern search in sequential data, like time series or data mapped to a genome sequence. Visually searching for patterns by similarity is often challenging because of the large search space, the visual complexity of patterns, and the user's perception of similarity. For example, in genomics, researchers try to link patterns in multivariate sequential data to cellular or pathogenic processes, but a lack of ground truth and high variance makes automatic pattern detection unreliable. We have developed a convolutional autoencoder for unsupervised representation learning of regions in sequential data that can capture more visual details of complex patterns compared to existing similarity measures. Using this learned representation as features of the sequential data, our accompanying visual query system enables interactive feedback-driven adjustments of the pattern search to adapt to the users' perceived similarity. Using an active learning sampling strategy, Peax collects user-generated binary relevance feedback. This feedback is used to train a model for binary classification, to ultimately find other regions that exhibit patterns similar to the search target. We demonstrate Peax's features through a case study in genomics and report on a user study with eight domain experts to assess the usability and usefulness of Peax. Moreover, we evaluate the effectiveness of the learned feature representation for visual similarity search in two additional user studies. We find that our models retrieve significantly more similar patterns than other commonly used techniques.

MEDEA: analysis of transcription factor binding motifs in accessible chromatin

Deciphering the interplay between chromatin accessibility and transcription factor (TF) binding is fundamental to understanding transcriptional regulation, control of cellular states, and the establishment of new phenotypes. Recent genome-wide chromatin accessibility profiling studies have provided catalogs of putative open regions, where TFs can recognize their motifs and regulate gene expression programs. Here, we present motif enrichment in differential elements of accessibility (MEDEA), a computational tool that analyzes high-throughput chromatin accessibility genomic data to identify cell-type-specific accessible regions and lineage-specific motifs associated with TF binding therein. To benchmark MEDEA, we used a panel of reference cell lines profiled by ENCODE and curated by the ENCODE Project Consortium for the ENCODE-DREAM Challenge. By comparing results with RNA-seq data, ChIP-seq peaks, and DNase-seq footprints, we show that MEDEA improves the detection of motifs associated with known lineage specifiers. We then applied MEDEA to 610 ENCODE DNase-seq data sets, where it revealed significant motifs even when absolute enrichment was low and where it identified novel regulators, such as NRF1 in kidney development. Finally, we show that MEDEA performs well on both bulk and single-cell ATAC-seq data. MEDEA is publicly available as part of our Glossary-GENRE suite for motif enrichment analysis.

Developmental Transcriptional Enhancers: A Subtle Interplay between Accessibility and Activity: Considering Quantitative Accessibility Changes between Different Regulatory States of an Enhancer Deconvolutes the Complex Relationship between Accessibility and Activity

Measurements of open chromatin in specific cell types are widely used to infer the spatiotemporal activity of transcriptional enhancers. How reliable are these predictions? In this review, it is argued that the relationship between the accessibility and activity of an enhancer is insufficiently described by simply considering open versus closed chromatin, or active versus inactive enhancers. Instead, recent studies focusing on the quantitative nature of accessibility signal reveal subtle differences between active enhancers and their different inactive counterparts: the closed silenced state and the accessible primed and repressed states. While the open structure as such is not a specific indicator of enhancer activity, active enhancers display a higher degree of accessibility than the primed and repressed states. Molecular mechanisms that may account for these quantitative differences are discussed. A model that relates molecular events at an enhancer to changes in its activity and accessibility in a developing tissue is also proposed.

CRUP: a comprehensive framework to predict condition-specific regulatory units

We present the software Condition-specific Regulatory Units Prediction (CRUP) to infer from epigenetic marks a list of regulatory units consisting of dynamically changing enhancers with their target genes. The workflow consists of a novel pre-trained enhancer predictor that can be reliably applied across cell types and species, solely based on histone modification ChIP-seq data. Enhancers are subsequently assigned to different conditions and correlated with gene expression to derive regulatory units. We thoroughly test and then apply CRUP to a rheumatoid arthritis model, identifying enhancer-gene pairs comprising known disease genes as well as new candidate genes.

EOMES interacts with RUNX3 and BRG1 to promote innate memory cell formation through epigenetic reprogramming

Memory CD8⁺ T cells have the ability to provide lifelong immunity against pathogens. Although memory features generally arise after challenge with a foreign antigen, naïve CD8 single positive (SP) thymocytes may acquire phenotypic and functional characteristics of memory cells in response to cytokines such as interleukin-4. This process is associated with the induction of the T-box transcription factor Eomesodermin (EOMES). However, the underlying molecular mechanisms remain ill-defined. Using epigenomic profiling, we show that these innate memory CD8SP cells acquire only a portion of the active enhancer repertoire of conventional memory cells. This reprograming is secondary to EOMES recruitment, mostly to RUNX3-bound enhancers. Furthermore, EOMES is found within chromatin-associated complexes containing BRG1 and promotes the recruitment of this chromatin remodelling factor. Also, the in vivo acquisition of EOMES-dependent program is BRG1-dependent. In conclusion, our results support a strong epigenetic basis for the EOMES-driven establishment of CD8⁺ T cell innate memory program.

Widespread epigenetic changes to the enhancer landscape of mouse liver induced by a specific xenobiotic agonist ligand of the nuclear receptor CAR

CAR (Nr1i3), a liver nuclear receptor and xenobiotic sensor, induces drug, steroid and lipid metabolism and dysregulates genes linked to hepatocellular carcinogenesis, but its impact on the liver epigenome is poorly understood. TCPOBOP, a halogenated xenochemical and highly specific CAR agonist ligand, induces localized chromatin opening or closing at several thousand mouse liver genomic regions, discovered as differential DNase-hypersensitive sites (ΔDHS). Active enhancer and promoter histone marks induced by TCPOBOP were enriched at opening DHS and TCPOBOP-inducible genes. Enrichment of CAR binding and CAR motifs was seen at opening DHS and their inducible drug/lipid metabolism gene targets, and at many constitutively open DHS located nearby. TCPOBOP-responsive cell cycle and DNA replication genes co-dependent on MET/EGFR signaling for induction were also enriched for CAR binding. A subset of opening DHS and many closing DHS mapping to TCPOBOP-responsive target genes did not bind CAR, indicating an indirect mechanism for their changes in chromatin accessibility. TCPOBOP-responsive DHS were also enriched for induced binding of RXRA, CEBPA and CEBPB, and for motifs for liver-enriched factors that may contribute to liver-specific transcriptional responses to TCPOBOP exposure. These studies elucidate the enhancer landscape of TCPOBOP-exposed liver and the widespread epigenetic changes that are induced by both direct and indirect mechanisms linked to CAR activation. The global maps of thousands of environmental chemical-induced epigenetic changes described here constitute a rich resource for further research on xenochemical effects on liver chromatin states and the epigenome.

Identification of factors associated with duplicate rate in ChIP-seq data

Chromatin immunoprecipitation and sequencing (ChIP-seq) has been widely used to map DNA-binding proteins, histone proteins and their modifications. ChIP-seq data contains redundant reads termed duplicates, referring to those mapping to the same genomic location and strand. There are two main sources of duplicates: polymerase chain reaction (PCR) duplicates and natural duplicates. Unlike natural duplicates that represent true signals from sequencing of independent DNA templates, PCR duplicates are artifacts originating from sequencing of identical copies amplified from the same DNA template. In analysis, duplicates are removed from peak calling and signal quantification. Nevertheless, a significant portion of the duplicates is believed to represent true signals. Obviously, removing all duplicates will underestimate the signal level in peaks and impact the identification of signal changes across samples. Therefore, an in-depth evaluation of the impact from duplicate removal is needed. Using eight public ChIP-seq datasets from three narrow-peak and two broad-peak marks, we tried to understand the distribution of duplicates in the genome, the extent by which duplicate removal impacts peak calling and signal estimation, and the factors associated with duplicate level in peaks. The three PCR-free histone H3 lysine 4 trimethylation (H3K4me3) ChIP-seq data had about 40% duplicates and 97% of them were within peaks. For the other datasets generated with PCR amplification of ChIP DNA, as expected, the narrow-peak marks have a much higher proportion of duplicates than the broad-peak marks. We found that duplicates are enriched in peaks and largely represent true signals, more conspicuous in those with high confidence. Furthermore, duplicate level in peaks is strongly correlated with the target enrichment level estimated using nonredundant reads, which provides the basis to properly allocate duplicates between noise and signal. Our analysis supports the feasibility of retaining the portion of signal duplicates into downstream analysis, thus alleviating the limitation of complete deduplication.