Background

Motif enrichment analysis (MEA) identifies over-represented transcription factor binding (TF) motifs in the DNA sequence of regulatory regions, enabling researchers to infer which transcription factors can regulate transcriptional response to a stimulus, or identify sequence features found near a target protein in a ChIP-seq experiment. Score-based MEA determines motifs enriched in regions exhibiting extreme differences in regulatory activity, but existing methods do not control for biases in GC content or dinucleotide composition. This lack of control for sequence bias, such as those often found in CpG islands, can obscure the enrichment of biologically relevant motifs.

Results

We developed Motif Enrichment In Ranked Lists of Peaks (MEIRLOP), a novel MEA method that determines enrichment of TF binding motifs in a list of scored regulatory regions, while controlling for sequence bias. In this study, we compare MEIRLOP against other MEA methods in identifying binding motifs found enriched in differentially active regulatory regions after interferon-beta stimulus, finding that using logistic regression and covariates improves the ability to call enrichment of ISGF3 binding motifs from differential acetylation ChIP-seq data compared to other methods. Our method achieves similar or better performance compared to other methods when quantifying the enrichment of TF binding motifs from ENCODE TF ChIP-seq datasets. We also demonstrate how MEIRLOP is broadly applicable to the analysis of numerous types of NGS assays and experimental designs.

Conclusions

Our results demonstrate the importance of controlling for sequence bias when accurately identifying enriched DNA sequence motifs using score-based MEA. MEIRLOP is available for download from https://github.com/npdeloss/meirlop under the MIT license.

How transcription affects genome 3D organization is not well understood. We found that during influenza A (IAV) infection, rampant transcription rapidly reorganizes host cell chromatin interactions. These changes occur at the ends of highly transcribed genes, where global inhibition of transcription termination by IAV NS1 protein causes readthrough transcription for hundreds of kilobases. In these readthrough regions, elongating RNA polymerase II disrupts chromatin interactions by inducing cohesin displacement from CTCF sites, leading to locus decompaction. Readthrough transcription into heterochromatin regions switches them from the inert (B) to the permissive (A) chromatin compartment and enables transcription factor binding. Data from non-viral transcription stimuli show that transcription similarly affects cohesin-mediated chromatin contacts within gene bodies. Conversely, inhibition of transcription elongation allows cohesin to accumulate at previously transcribed intragenic CTCF sites and to mediate chromatin looping and compaction. Our data indicate that transcription elongation by RNA polymerase II remodels genome 3D architecture.

Cells regulate gene expression in response to salient external stimuli. In neurons, depolarization leads to the expression of inducible transcription factors (ITFs) that direct subsequent gene regulation. Depolarization encodes both a neuron's action potential (AP) output and synaptic inputs, via excitatory postsynaptic potentials (EPSPs). However, it is unclear if distinct types of electrical activity can be transformed by an ITF into distinct modes of genomic regulation. Here, we show that APs and EPSPs in mouse hippocampal neurons trigger two spatially segregated and molecularly distinct induction mechanisms that lead to the expression of the ITF NPAS4. These two pathways culminate in the formation of stimulus-specific NPAS4 heterodimers that exhibit distinct DNA binding patterns. Thus, NPAS4 differentially communicates increases in a neuron's spiking output and synaptic inputs to the nucleus, enabling gene regulation to be tailored to the type of depolarizing activity along the somato-dendritic axis of a neuron.

Tissue environment plays a powerful role in establishing and maintaining the distinct phenotypes of resident macrophages, but the underlying molecular mechanisms remain poorly understood. Here, we characterized transcriptomic and epigenetic changes in repopulating liver macrophages following acute Kupffer cell depletion as a means to infer signaling pathways and transcription factors that promote Kupffer cell differentiation. We obtained evidence that combinatorial interactions of the Notch ligand DLL4 and transforming growth factor-b (TGF-β) family ligands produced by sinusoidal endothelial cells and endogenous LXR ligands were required for the induction and maintenance of Kupffer cell identity. DLL4 regulation of the Notch transcriptional effector RBPJ activated poised enhancers to rapidly induce LXRα and other Kupffer cell lineage-determining factors. These factors in turn reprogrammed the repopulating liver macrophage enhancer landscape to converge on that of the original resident Kupffer cells. Collectively, these findings provide a framework for understanding how macrophage progenitor cells acquire tissue-specific phenotypes.