Factorbook: an updated catalog of transcription factor motifs and candidate regulatory motif sites.

The human genome contains ∼2000 transcriptional regulatory proteins, including ∼1600 DNA-binding transcription factors (TFs) recognizing characteristic sequence motifs to exert regulatory effects on gene expression. The binding specificities of these factors have been profiled both in vitro, using techniques such as HT-SELEX, and in vivo, using techniques including ChIP-seq. We previously developed Factorbook, a TF-centric database of annotations, motifs, and integrative analyses based on ChIP-seq data from Phase II of the ENCODE Project. Here we present an update to Factorbook which significantly expands the breadth of cell type and TF coverage. The update includes an expanded motif catalog derived from thousands of ENCODE Phase II and III ChIP-seq experiments and HT-SELEX experiments; this motif catalog is integrated with the ENCODE registry of candidate cis-regulatory elements to annotate a comprehensive collection of genome-wide candidate TF binding sites. The database also offers novel tools for applying the motif models within machine learning frameworks and using these models for integrative analysis, including annotation of variants and disease and trait heritability. Factorbook is publicly available at www.factorbook.org; we will continue to expand the resource as ENCODE Phase IV data are released.

Annotation of chromatin states in 66 complete mouse epigenomes during development.

The morphologically and functionally distinct cell types of a multicellular organism are maintained by their unique epigenomes and gene expression programs. Phase III of the ENCODE Project profiled 66 mouse epigenomes across twelve tissues at daily intervals from embryonic day 11.5 to birth. Applying the ChromHMM algorithm to these epigenomes, we annotated eighteen chromatin states with characteristics of promoters, enhancers, transcribed regions, repressed regions, and quiescent regions. Our integrative analyses delineate the tissue specificity and developmental trajectory of the loci in these chromatin states. Approximately 0.3% of each epigenome is assigned to a bivalent chromatin state, which harbors both active marks and the repressive mark H3K27me3. Highly evolutionarily conserved, these loci are enriched in silencers bound by polycomb repressive complex proteins, and the transcription start sites of their silenced target genes. This collection of chromatin state assignments provides a useful resource for studying mammalian development.

Expanded encyclopaedias of DNA elements in the human and mouse genomes.

The human and mouse genomes contain instructions that specify RNAs and proteins and govern the timing, magnitude, and cellular context of their production. To better delineate these elements, phase III of the Encyclopedia of DNA Elements (ENCODE) Project has expanded analysis of the cell and tissue repertoires of RNA transcription, chromatin structure and modification, DNA methylation, chromatin looping, and occupancy by transcription factors and RNA-binding proteins. Here we summarize these efforts, which have produced 5,992 new experimental datasets, including systematic determinations across mouse fetal development. All data are available through the ENCODE data portal (https://www.encodeproject.org), including phase II ENCODE1 and Roadmap Epigenomics2 data. We have developed a registry of 926,535 human and 339,815 mouse candidate cis-regulatory elements, covering 7.9 and 3.4% of their respective genomes, by integrating selected datatypes associated with gene regulation, and constructed a web-based server (SCREEN; http://screen.encodeproject.org) to provide flexible, user-defined access to this resource. Collectively, the ENCODE data and registry provide an expansive resource for the scientific community to build a better understanding of the organization and function of the human and mouse genomes.

A curated benchmark of enhancer-gene interactions for evaluating enhancer-target gene prediction methods.

BACKGROUND: Many genome-wide collections of candidate cis-regulatory elements (cCREs) have been defined using genomic and epigenomic data, but it remains a major challenge to connect these elements to their target genes. RESULTS: To facilitate the development of computational methods for predicting target genes, we develop a Benchmark of candidate Enhancer-Gene Interactions (BENGI) by integrating the recently developed Registry of cCREs with experimentally derived genomic interactions. We use BENGI to test several published computational methods for linking enhancers with genes, including signal correlation and the TargetFinder and PEP supervised learning methods. We find that while TargetFinder is the best-performing method, it is only modestly better than a baseline distance method for most benchmark datasets when trained and tested with the same cell type and that TargetFinder often does not outperform the distance method when applied across cell types. CONCLUSIONS: Our results suggest that current computational methods need to be improved and that BENGI presents a useful framework for method development and testing.

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

Enhancers are distal cis-regulatory elements that modulate gene expression. They are depleted of nucleosomes and enriched in specific histone modifications; thus, calling DNase-seq and histone mark ChIP-seq peaks can predict enhancers. We evaluated nine peak-calling algorithms for predicting enhancers validated by transgenic mouse assays. DNase and H3K27ac peaks were consistently more predictive than H3K4me1/2/3 and H3K9ac peaks. DFilter and Hotspot2 were the best DNase peak callers, while HOMER, MUSIC, MACS2, DFilter and F-seq were the best H3K27ac peak callers. We observed that the differential DNase or H3K27ac signals between two distant tissues increased the area under the precision-recall curve (PR-AUC) of DNase peaks by 17.5-166.7% and that of H3K27ac peaks by 7.1-22.2%. We further improved this differential signal method using multiple contrast tissues. Evaluated using a blind test, the differential H3K27ac signal method substantially improved PR-AUC from 0.48 to 0.75 for predicting heart enhancers. We further validated our approach using postnatal retina and cerebral cortex enhancers identified by massively parallel reporter assays, and observed improvements for both tissues. In summary, we compared nine peak callers and devised a superior method for predicting tissue-specific mouse developmental enhancers by reranking the called peaks.

Where fMRI and electrophysiology agree to disagree: corticothalamic and striatal activity patterns in the WAG/Rij rat.

The relationship between neuronal activity and hemodynamic changes plays a central role in functional neuroimaging. Under normal conditions and in neurological disorders such as epilepsy, it is commonly assumed that increased functional magnetic resonance imaging (fMRI) signals reflect increased neuronal activity and that fMRI decreases represent neuronal activity decreases. Recent work suggests that these assumptions usually hold true in the cerebral cortex. However, less is known about the basis of fMRI signals from subcortical structures such as the thalamus and basal ganglia. We used WAG/Rij rats (Wistar albino Glaxo rats of Rijswijk), an established animal model of human absence epilepsy, to perform fMRI studies with blood oxygen level-dependent and cerebral blood volume (CBV) contrasts at 9.4 tesla, as well as laser Doppler cerebral blood flow (CBF), local field potential (LFP), and multiunit activity (MUA) recordings. We found that, during spike-wave discharges, the somatosensory cortex and thalamus showed increased fMRI, CBV, CBF, LFP, and MUA signals. However, the caudate-putamen showed fMRI, CBV, and CBF decreases despite increases in LFP and MUA signals. Similarly, during normal whisker stimulation, the cortex and thalamus showed increases in CBF and MUA, whereas the caudate-putamen showed decreased CBF with increased MUA. These findings suggest that neuroimaging-related signals and electrophysiology tend to agree in the cortex and thalamus but disagree in the caudate-putamen. These opposite changes in vascular and electrical activity indicate that caution should be applied when interpreting fMRI signals in both health and disease from the caudate-putamen, as well as possibly from other subcortical structures.

Impaired consciousness in temporal lobe seizures: role of cortical slow activity.

Impaired consciousness requires altered cortical function. This can occur either directly from disorders that impair widespread bilateral regions of the cortex or indirectly through effects on subcortical arousal systems. It has therefore long been puzzling why focal temporal lobe seizures so often impair consciousness. Early work suggested that altered consciousness may occur with bilateral or dominant temporal lobe seizure involvement. However, other bilateral temporal lobe disorders do not impair consciousness. More recent work supports a 'network inhibition hypothesis' in which temporal lobe seizures disrupt brainstem-diencephalic arousal systems, leading indirectly to depressed cortical function and impaired consciousness. Indeed, prior studies show subcortical involvement in temporal lobe seizures and bilateral frontoparietal slow wave activity on intracranial electroencephalography. However, the relationships between frontoparietal slow waves and impaired consciousness and between cortical slowing and fast seizure activity have not been directly investigated. We analysed intracranial electroencephalography recordings during 63 partial seizures in 26 patients with surgically confirmed mesial temporal lobe epilepsy. Behavioural responsiveness was determined based on blinded review of video during seizures and classified as impaired (complex-partial seizures) or unimpaired (simple-partial seizures). We observed significantly increased delta-range 1-2 Hz slow wave activity in the bilateral frontal and parietal neocortices during complex-partial compared with simple-partial seizures. In addition, we confirmed prior work suggesting that propagation of unilateral mesial temporal fast seizure activity to the bilateral temporal lobes was significantly greater in complex-partial than in simple-partial seizures. Interestingly, we found that the signal power of frontoparietal slow wave activity was significantly correlated with the temporal lobe fast seizure activity in each hemisphere. Finally, we observed that complex-partial seizures were somewhat more common with onset in the language-dominant temporal lobe. These findings provide direct evidence for cortical dysfunction in the form of bilateral frontoparietal slow waves associated with impaired consciousness in temporal lobe seizures. We hypothesize that bilateral temporal lobe seizures may exert a powerful inhibitory effect on subcortical arousal systems. Further investigations will be needed to fully determine the role of cortical-subcortical networks in ictal neocortical dysfunction and may reveal treatments to prevent this important negative consequence of temporal lobe epilepsy.

A prospective study of loss of consciousness in epilepsy using virtual reality driving simulation and other video games.

Patients with epilepsy are at risk of traffic accidents when they have seizures while driving. However, driving is an essential part of normal daily life in many communities, and depriving patients of driving privileges can have profound consequences for their economic and social well-being. In the current study, we collected ictal performance data from a driving simulator and two other video games in patients undergoing continuous video/EEG monitoring. We captured 22 seizures in 13 patients and found that driving impairment during seizures differed in terms of both magnitude and character, depending on the seizure type. Our study documents the feasibility of a prospective study of driving and other behaviors during seizures through the use of computer-based tasks. This methodology may be applied to further describe differential driving impairment in specific types of seizures and to gain data on anatomical networks disrupted in seizures that impair consciousness and driving safety.

Focal BOLD fMRI changes in bicuculline-induced tonic-clonic seizures in the rat.

Generalized tonic-clonic seizures cause widespread physiological changes throughout the cerebral cortex and subcortical structures in the brain. Using combined blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) at 9.4 T and electroencephalography (EEG), these changes can be characterized with high spatiotemporal resolution. We studied BOLD changes in anesthetized Wistar rats during bicuculline-induced tonic-clonic seizures. Bicuculline, a GABA(A) receptor antagonist, was injected systemically and seizure activity was observed on EEG as high-amplitude, high-frequency polyspike discharges followed by clonic paroxysmal activity of lower frequency, with mean electrographic seizure duration of 349 s. Our aim was to characterize the spatial localization, direction, and timing of BOLD signal changes during the pre-ictal, ictal and post-ictal periods. Group analysis was performed across seizures using paired t-maps of BOLD signal superimposed on high-resolution anatomical images. Regional analysis was then performed using volumes of interest to quantify BOLD timecourses. In the pre-ictal period we found focal BOLD increases in specific areas of somatosensory cortex (S1, S2) and thalamus several seconds before seizure onset. During seizures we observed BOLD increases in cortex, brainstem and thalamus and BOLD decreases in the hippocampus. The largest ictal BOLD increases remained in the focal regions of somatosensory cortex showing pre-ictal increases. During the post-ictal period we observed widespread BOLD decreases. These findings support a model in which "generalized" tonic-clonic seizures begin with focal changes before electrographic seizure onset, which progress to non-uniform changes during seizures, possibly shedding light on the etiology and pathophysiology of similar seizures in humans.