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.

ATLAS: A database linking binding affinities with structures for wild-type and mutant TCR-pMHC complexes.

The ATLAS (Altered TCR Ligand Affinities and Structures) database (https://zlab.umassmed.edu/atlas/web/) is a manually curated repository containing the binding affinities for wild-type and mutant T cell receptors (TCRs) and their antigens, peptides presented by the major histocompatibility complex (pMHC). The database links experimentally measured binding affinities with the corresponding three dimensional (3D) structures for TCR-pMHC complexes. The user can browse and search affinities, structures, and experimental details for TCRs, peptides, and MHCs of interest. We expect this database to facilitate the development of next-generation protein design algorithms targeting TCR-pMHC interactions. ATLAS can be easily parsed using modeling software that builds protein structures for training and testing. As an example, we provide structural models for all mutant TCRs in ATLAS, built using the Rosetta program. Utilizing these structures, we report a correlation of 0.63 between experimentally measured changes in binding energies and our predicted changes. Proteins 2017; 85:908-916. © 2016 Wiley Periodicals, Inc.

Space-time wiring specificity supports direction selectivity in the retina.

How does the mammalian retina detect motion? This classic problem in visual neuroscience has remained unsolved for 50 years. In search of clues, here we reconstruct Off-type starburst amacrine cells (SACs) and bipolar cells (BCs) in serial electron microscopic images with help from EyeWire, an online community of 'citizen neuroscientists'. On the basis of quantitative analyses of contact area and branch depth in the retina, we find evidence that one BC type prefers to wire with a SAC dendrite near the SAC soma, whereas another BC type prefers to wire far from the soma. The near type is known to lag the far type in time of visual response. A mathematical model shows how such 'space-time wiring specificity' could endow SAC dendrites with receptive fields that are oriented in space-time and therefore respond selectively to stimuli that move in the outward direction from the soma.

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.

Simultaneous EEG, fMRI, and behavior in typical childhood absence seizures.

PURPOSE: Absence seizures cause transient impairment of consciousness. Typical absence seizures occur in children, and are accompanied by 3-4-Hz spike-wave discharges (SWDs) on electroencephalography (EEG). Prior EEG-functional magnetic resonance imaging (fMRI) studies of SWDs have shown a network of cortical and subcortical changes during these electrical events. However, fMRI during typical childhood absence seizures with confirmed impaired consciousness has not been previously investigated. METHODS: We performed EEG-fMRI with simultaneous behavioral testing in 37 children with typical childhood absence epilepsy (CAE). Attentional vigilance was evaluated by a continuous performance task (CPT), and simpler motor performance was evaluated by a repetitive tapping task (RTT). RESULTS: SWD episodes were obtained during fMRI scanning from 9 patients among the 37 studied. fMRI signal increases during SWDs were observed in the thalamus, frontal cortex, primary visual, auditory, somatosensory, and motor cortex, and fMRI decreases were seen in the lateral and medial parietal cortex, cingulate gyrus, and basal ganglia. Omission error rate (missed targets) with SWDs during fMRI was 81% on CPT and 39% on RTT. For those seizure epochs during which CPT performance was impaired, fMRI changes were seen in cortical and subcortical structures typically involved in SWDs, whereas minimal changes were observed for the few epochs during which performance was spared. DISCUSSION: These findings suggest that typical absence seizures involve a network of cortical-subcortical areas necessary for normal attention and primary information processing. Identification of this network may improve understanding of cognitive impairments in CAE, and may help guide development of new therapies for this disorder.

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.