TherMos: Estimating protein-DNA binding energies from in vivo binding profiles

Accurately characterizing transcription factor (TF)-DNA affinity is a central goal of regulatory genomics. Although thermodynamics provides the most natural language for describing the continuous range of TF-DNA affinity, traditional motif discovery algorithms focus instead on classification paradigms that aim to discriminate 'bound' and 'unbound' sequences. Moreover, these algorithms do not directly model the distribution of tags in ChIP-seq data. Here, we present a new algorithm named Thermodynamic Modeling of ChIP-seq (TherMos), which directly estimates a positionspecific binding energy matrix (PSEM) from ChIPseq/exo tag profiles. In cross-validation tests on seven genome-wide TF-DNA binding profiles, one of which we generated via ChIP-seq on a complex developing tissue, TherMos predicted quantitative TF-DNA binding with greater accuracy than five well-known algorithms. We experimentally validated TherMos binding energy models for Klf4 and Esrrb, using a novel protocol to measure PSEMs in vitro. Strikingly, our measurements revealed strong nonadditivity at multiple positions within the two PSEMs. Among the algorithms tested, only TherMos was able to model the entire binding energy landscape of Klf4 and Esrrb. Our study reveals new insights into the energetics of TF-DNA binding in vivo and provides an accurate first-principles approach to binding energy inference from ChIP-seq and ChIP-exo data. © 2013 The Author(s).Link_to_subscribed_fulltex

The transcription factor SOX6 contributes to the developmental origins of obesity by promoting adipogenesis

10.1242/dev.131573Development (Cambridge, England)1436950-961GUSTO (Growing up towards Healthy Outcomes

Ultraconservation identifies a small subset of extremely constrained developmental enhancers

While experimental studies have suggested that non-coding ultraconserved DNA elements are central nodes in the regulatory circuitry that specifies mammalian embryonic development, the possible functional relevance of their>200bp of perfect sequence conservation between human-mouse-rat remains obscure 1,2. Here we have compared the in vivo enhancer activity of a genome-wide set of 231 non-exonic sequences with ultraconserved cores to that of 206 sequences that are under equivalently severe human-rodent constraint (ultra-like), but lack perfect sequence conservation. In transgenic mouse assays, 50percent of the ultraconserved and 50percent of the ultra-like conserved elements reproducibly functioned as tissue-specific enhancers at embryonic day 11.5. In this in vivo assay, we observed that ultraconserved enhancers and constrained non-ultraconserved enhancers targeted expression to a similar spectrum of tissues with a particular enrichment in the developing central nervous system. A human genome-wide comparative screen uncovered ~;;2,600 non-coding elements that evolved under ultra-like human-rodent constraint and are similarly enriched near transcriptional regulators and developmental genes as the much smaller number of ultraconserved elements. These data indicate that ultraconserved elements possessing absolute human-rodent sequence conservation are not distinct from other non-coding elements that are under comparable purifying selection in mammals and suggest they are principal constituents of the cis-regulatory framework of mammalian development

Intrinsic mutant HTT-mediated defects in oligodendroglia cause myelination deficits and behavioral abnormalities in Huntington disease

White matter abnormalities are a nearly universal pathological feature of neurodegenerative disorders including Huntington disease (HD). A long-held assumption is that this white matter pathology is simply a secondary outcome of the progressive neuronal loss that manifests with advancing disease. Using a mouse model of HD, here we show that white matter and myelination abnormalities are an early disease feature appearing before the manifestation of any behavioral abnormalities or neuronal loss. We further show that selective inactivation of mutant huntingtin (mHTT) in the NG2+ oligodendrocyte progenitor cell population prevented myelin abnormalities and certain behavioral deficits in HD mice. Strikingly, the improvements in behavioral outcomes were seen despite the continued expression of mHTT in nonoligodendroglial cells including neurons, astrocytes, and microglia. Using RNA-seq and ChIP-seq analyses, we implicate a pathogenic mechanism that involves enhancement of polycomb repressive complex 2 (PRC2) activity by mHTT in the intrinsic oligodendroglial dysfunction and myelination deficits observed in HD. Our findings challenge the long-held dogma regarding the etiology of white matter pathology in HD and highlight the contribution of epigenetic mechanisms to the observed intrinsic oligodendroglial dysfunction. Our results further suggest that ameliorating white matter pathology and oligodendroglial dysfunction may be beneficial for HD.Agency for Science, Technology and Research (A*STAR)Accepted versionWe thank members of the M.A.P. laboratory for helpful discussions and comments. C.F.B. is supported by a Singapore International Graduate Award from the Agency for Science, Technology and Research (A*STAR). M.A.P. is supported by grants from A*STAR and the National University of Singapore

The complete sequence of human chromosome 5

Chromosome 5 is one of the largest human chromosomes yet has one of the lowest gene densities. This is partially explained by numerous gene-poor regions that display a remarkable degree of noncoding and syntenic conservation with non-mammalian vertebrates, suggesting they are functionally constrained. In total, we compiled 177.7 million base pairs of highly accurate finished sequence containing 923 manually curated protein-encoding genes including the protocadherin and interleukin gene families and the first complete versions of each of the large chromosome 5 specific internal duplications. These duplications are very recent evolutionary events and play a likely mechanistic role, since deletions of these regions are the cause of debilitating disorders including spinal muscular atrophy (SMA)

Eye-View: A retinal image capture system

This paper presents a low cost but high resolution retinal image acquisition system of the human eye. The images acquired by a CMOS image sensor are communicated through the Universal Serial Bus (USB) interface to a personal computer for viewing and further processing. The image acquisition time was estimated to be 2.5 seconds. This system can also be used in telemedicine applications

Transposable elements that have recently been mobile in the human genome

10.1186/s12864-021-08085-0BMC Genomics22178

Integrative model of genomic factors for determining binding site selection by estrogen receptor‐α

A major question in transcription factor (TF) biology is why a TF binds to only a small fraction of motif eligible binding sites in the genome. Using the estrogen receptor-α as a model system, we sought to explicitly define parameters that determine TF-binding site selection. By examining 12 genetic and epigenetic parameters, we find that an energetically favorable estrogen response element (ERE) motif sequence, co-occupancy by the TF FOXA1, the presence of the H3K4me1 mark and an open chromatin configuration in the pre-ligand state provide specificity for ER binding. These factors can model estrogen-induced ER binding with high accuracy (ROC-AUC=0.95 and 0.88 using different genomic backgrounds). Moreover, when assessed in another estrogen-responsive cell line, this model was highly predictive for ERα binding (ROC-AUC=0.86). Variance in binding site selection between MCF-7 and T47D resides in sites with suboptimal ERE motifs, but modulated by the chromatin configuration. These results suggest a definable interplay between sequence motifs and local chromatin in selecting TF binding