Pronounced sequence specificity of the TET enzyme catalytic domain guides its cellular function

TET (ten-eleven translocation) enzymes catalyze the oxidation of 5-methylcytosine bases in DNA, thus driving active and passive DNA demethylation. Here, we report that the catalytic cores of mammalian TET enzymes favor CpGs embedded within bHLH and bZIP transcription factor binding sites, with 250-fold preference in vitro. Crystal structures and molecular dynamics calculations show that sequence preference is caused by intra-substrate interactions and CpG flanking sequence indirectly affecting enzyme conformation. TET sequence preferences are physiologically relevant as they explain the rates of DNA demethylation in TET-rescue experiments in culture and in vivo within the zygote and germline. Most and least favorable TET motifs represent DNA sites that are bound by methylation-sensitive immediate-early transcription factors and OCT4, respectively, illuminating TET function in transcriptional responses and pluripotency support. One-Sentence Summary: The catalytic domains of the enzymes that facilitate passive and drive active DNA demethylation have intrinsic sequence preferences that target DNA demethylation to bHLH and bZIP transcription factor binding sites.

Database and Statistical Analyses of Transcription Factor Binding Sites in the Non-Coding Control Region of JC Virus

JC virus (JCV), as an archetype, establishes a lifelong latent or persistent infection in many healthy individuals. In immunocompromised patients, prototype JCV with variable mutations in the non-coding control region (NCCR) causes progressive multifocal leukoencephalopathy (PML), a severe demyelinating disease. This study was conducted to create a database of NCCR sequences annotated with transcription factor binding sites (TFBSs) and statistically analyze the mutational pattern of the JCV NCCR. JCV NCCRs were extracted from >1000 sequences registered in GenBank, and TFBSs within each NCCR were identified by computer simulation, followed by examination of their prevalence, multiplicity, and location by statistical analyses. In the NCCRs of the prototype JCV, the limited types of TFBSs, which are mainly present in regions D through F of archetype JCV, were significantly reduced. By contrast, modeling count data revealed that several TFBSs located in regions C and E tended to overlap in the prototype NCCRs. Based on data from the BioGPS database, genes encoding transcription factors that bind to these TFBSs were expressed not only in the brain but also in the peripheral sites. The database and NCCR patterns obtained in this study could be a suitable platform for analyzing JCV mutations and pathogenicity.

Global properties of regulatory sequences are predicted by transcription factor recognition mechanisms

Background Mammalian genomes contain millions of putative regulatory sequences, which are delineated by binding of multiple transcription factors. The degree to which spacing and orientation constraints among transcription factor binding sites contribute to the recognition and identity of regulatory sequence is an unresolved but important question that impacts our understanding of genome function and evolution. Global mechanisms that underlie phenomena including the size of regulatory sequences, their uniqueness, and their evolutionary turnover remain poorly described. Results Here, we ask whether models incorporating different degrees of spacing and orientation constraints among transcription factor binding sites are broadly consistent with several global properties of regulatory sequence. These properties include length, sequence diversity, turnover rate, and dominance of specific TFs in regulatory site identity and cell type specification. Models with and without spacing and orientation constraints are generally consistent with all observed properties of regulatory sequence, and with regulatory sequences being fundamentally small (~ 1 nucleosome). Uniqueness of regulatory regions and their rapid evolutionary turnover are expected under all models examined. An intriguing issue we identify is that the complexity of eukaryotic regulatory sites must scale with the number of active transcription factors, in order to accomplish observed specificity. Conclusions Models of transcription factor binding with or without spacing and orientation constraints predict that regulatory sequences should be fundamentally short, unique, and turn over rapidly. We posit that the existence of master regulators may be, in part, a consequence of evolutionary pressure to limit the complexity and increase evolvability of regulatory sites.

Regulatory SNPs: Altered Transcription Factor Binding Sites Implicated in Complex Traits and Diseases

The vast majority of the genetic variants (mainly SNPs) associated with various human traits and diseases map to a noncoding part of the genome and are enriched in its regulatory compartment, suggesting that many causal variants may affect gene expression. The leading mechanism of action of these SNPs consists in the alterations in the transcription factor binding via creation or disruption of transcription factor binding sites (TFBSs) or some change in the affinity of these regulatory proteins to their cognate sites. In this review, we first focus on the history of the discovery of regulatory SNPs (rSNPs) and systematized description of the existing methodical approaches to their study. Then, we brief the recent comprehensive examples of rSNPs studied from the discovery of the changes in the TFBS sequence as a result of a nucleotide substitution to identification of its effect on the target gene expression and, eventually, to phenotype. We also describe state-of-the-art genome-wide approaches to identification of regulatory variants, including both making molecular sense of genome-wide association studies (GWAS) and the alternative approaches the primary goal of which is to determine the functionality of genetic variants. Among these approaches, special attention is paid to expression quantitative trait loci (eQTLs) analysis and the search for allele-specific events in RNA-seq (ASE events) as well as in ChIP-seq, DNase-seq, and ATAC-seq (ASB events) data.

Differential gene regulation in DAPT-treated Hydra reveals candidate direct notch signalling targets

In Hydra, Notch inhibition causes defects in head patterning and prevents differentiation of proliferating nematocyte progenitor cells into mature nematocytes. To understand the molecular mechanisms by which the Notch pathway regulates these processes we performed RNAseq and identified genes that are differentially regulated in response to 48 hours of treating the animals with the Notch-inhibitor DAPT. To identify candidate direct regulators of Notch-signalling, we profiled gene expression changes that occur during subsequent restoration of Notch-activity and performed promoter analyses to identify RBPJ transcription factor binding sites in the regulatory regions of Notch-responsive genes. Interrogating the available single cell sequencing data set revealed gene expression patterns of Notch-regulated Hydra genes. By these analyses a comprehensive picture of the molecular pathways regulated by Notch signalling in head patterning and in interstitial cell differentiation in Hydra emerged. As prime candidates for direct Notch-target genes, in addition to HyHes, we suggest Sp5 and HyAlx. They rapidly recovered their expression levels after DAPT removal and possess Notch-responsive RBPJ transcription factor binding sites in their regulatory regions.