X chromosome inactivation (XCI) serves as a paradigm for RNA-mediated regulation of gene expression, wherein the long non-coding RNA XIST spreads across the X chromosome in cis to mediate gene silencing chromosome-wide. In female naive human pluripotent stem cells (hPSCs), XIST is in a dispersed configuration, and XCI does not occur, raising questions about XIST's function. We found that XIST spreads across the X chromosome and induces dampening of X-linked gene expression in naive hPSCs. Surprisingly, XIST also targets specific autosomal regions, where it induces repressive chromatin changes and gene expression dampening. Thereby, XIST equalizes X-linked gene dosage between male and female cells while inducing differences in autosomes. The dispersed Xist configuration and autosomal localization also occur transiently during XCI initiation in mouse PSCs. Together, our study identifies XIST as the regulator of X chromosome dampening, uncovers an evolutionarily conserved trans-acting role of XIST/Xist, and reveals a correlation between XIST/Xist dispersal and autosomal targeting.

To regulate gene expression, the macromolecular components of the mammalian interphase nucleus are spatially organized into a myriad of functional compartments. Over the past decade, increasingly sophisticated genomics, microscopy, and functional approaches have probed this organization in unprecedented detail. These investigations have linked chromatin-associated noncoding RNAs to specific nuclear compartments and uncovered mechanisms by which these RNAs establish such domains. In this review, we focus on the long non-coding RNA Xist and summarize new evidence demonstrating the significance of chromatin reconfiguration in creating the inactive X-chromosome compartment. Differences in chromatin compaction correlate with distinct levels of gene repression on the X-chromosome, potentially explaining how human XIST can induce chromosome-wide dampening and silencing of gene expression at different stages of human development.

Skin pigmentation is paused after sun exposure; however, the mechanism behind this pausing is unknown. In this study, we found that the UVB-induced DNA repair system, led by the ataxia telangiectasia mutated (ATM) protein kinase, represses MITF transcriptional activity of pigmentation genes while placing MITF in DNA repair mode, thus directly inhibiting pigment production. Phosphoproteomics analysis revealed ATM to be the most significantly enriched pathway among all UVB-induced DNA repair systems. ATM inhibition in mouse or human skin, either genetically or chemically, induces pigmentation. Upon UVB exposure, MITF transcriptional activation is blocked owing to ATM-dependent phosphorylation of MITF on S414, which modifies MITF activity and interactome toward DNA repair, including binding to TRIM28 and RBBP4. Accordingly, MITF genome occupancy is enriched in sites of high DNA damage that are likely repaired. This suggests that ATM harnesses the pigmentation key activator for the necessary rapid, efficient DNA repair, thus optimizing the chances of the cell surviving. Data are available from ProteomeXchange with the identifier PXD041121.

Human germ cell development is a highly regulated process beginning soon after embryo implantation with the specification of primordial germ cells (PGCs) and ending in adulthood with the differentiation of gametes. Here, we show that fibroblast growth factor receptor 3 (FGFR3) is expressed by human PGCs during the first and second trimester, becoming repressed as PGCs differentiate into primordial oocytes. Using fluorescence-activated cell sorting (FACS) with antibodies that recognize FGFR3 followed by single-cell RNA sequencing, we show that isolating FGFR3-positive cells enriches for human PGCs. Taken together, FGFR3 could be used in future studies as a strategy to identify maturing PGCs in vitro.

The lncRNA Xist forms ∼50 diffraction-limited foci to transcriptionally silence one X chromosome. How this small number of RNA foci and interacting proteins regulate a much larger number of X-linked genes is unknown. We show that Xist foci are locally confined, contain ∼2 RNA molecules, and nucleate supramolecular complexes (SMACs) that include many copies of the critical silencing protein SPEN. Aggregation and exchange of SMAC proteins generate local protein gradients that regulate broad, proximal chromatin regions. Partitioning of numerous SPEN molecules into SMACs is mediated by their intrinsically disordered regions and essential for transcriptional repression. Polycomb deposition via SMACs induces chromatin compaction and the increase in SMACs density around genes, which propagates silencing across the X chromosome. Our findings introduce a mechanism for functional nuclear compartmentalization whereby crowding of transcriptional and architectural regulators enables the silencing of many target genes by few RNA molecules.

The long non-coding RNA Xist exploits numerous effector proteins to progressively induce gene silencing across the X chromosome and form the inactive X (Xi)-compartment. The mechanism underlying formation of the chromosome-wide Xi-compartment is poorly understood. Here, we find that formation of the Xi-compartment is induced by ∼50 locally confined granules, where two Xist RNA molecules nucleate s upra- m olecular c omplexes (SMCs) of interacting proteins. Xist-SMCs are transient structures that concentrate rapidly recycling proteins in the X by increasing protein binding affinity. We find that gene silencing originates at Xist-SMCs and propagates across the entire chromosome over time, achieved by Polycomb-mediated coalescence of chromatin regions and aggregation, via its intrinsically disordered domains, of the critical silencing factor SPEN. Our results suggest a new model for X chromosome inactivation, in which Xist RNA induces macromolecular crowding of heterochromatinizing proteins near distinct sites which ultimately increases their density throughout the chromosome. This mechanism enables deterministic gene silencing without the need for Xist ribonucleoprotein complex-chromatin interactions at each target gene.

X-chromosome dosage compensation in female placental mammals is achieved by X-chromosome inactivation (XCI). Human pre-implantation embryos are an exception, in which dosage compensation occurs by X-chromosome dampening (XCD). Here, we examined whether XCD extends to human prenatal germ cells given their similarities to naive pluripotent cells. We found that female human primordial germ cells (hPGCs) display reduced X-linked gene expression before entering meiosis. Moreover, in hPGCs, both X chromosomes are active and express the long non-coding RNAs X active coating transcript (XACT) and X inactive specific transcript (XIST)-the master regulator of XCI-which are silenced after entry into meiosis. We find that XACT is a hPGC marker, describe XCD associated with XIST expression in hPGCs and suggest that XCD evolved in humans to regulate X-linked genes in pre-implantation embryos and PGCs. Furthermore, we found a unique mechanism of X-chromosome regulation in human primordial oocytes. Therefore, future studies of human germline development must consider the sexually dimorphic X-chromosome dosage compensation mechanisms in the prenatal germline.

Skin sun exposure induces two protection programs: stress responses and pigmentation, the former within minutes and the latter only hours afterward. Although serving the same physiological purpose, it is not known whether and how these programs are coordinated. Here, we report that UVB exposure every other day induces significantly more skin pigmentation than the higher frequency of daily exposure, without an associated increase in stress responses. Using mathematical modeling and empirical studies, we show that the melanocyte master regulator, MITF, serves to synchronize stress responses and pigmentation and, furthermore, functions as a UV-protection timer via damped oscillatory dynamics, thereby conferring a trade-off between the two programs. MITF oscillations are controlled by multiple negative regulatory loops, one at the transcriptional level involving HIF1α and another post-transcriptional loop involving microRNA-148a. These findings support trait linkage between the two skin protection programs, which, we speculate, arose during furless skin evolution to minimize skin damage.

Ectopic expression of lineage master regulators induces transdifferentiation. Whether cell fate transitions can be induced during various developmental stages has not been systemically examined. Here we discover that amongst different developmental stages, mouse embryonic stem cells (mESCs) are resistant to cell fate conversion induced by the melanocyte lineage master regulator MITF. By generating a transgenic system we exhibit that in mESCs, the pluripotency master regulator Oct4, counteracts pro-differentiation induced by Mitf by physical interference with MITF transcriptional activity. We further demonstrate that mESCs must be released from Oct4-maintained pluripotency prior to ectopically induced differentiation. Moreover, Oct4 induction in various differentiated cells represses their lineage identity in vivo. Alongside, chromatin architecture combined with ChIP-seq analysis suggest that Oct4 competes with various lineage master regulators for binding promoters and enhancers. Our analysis reveals pluripotency and transdifferentiation regulatory principles and could open new opportunities in the field of regenerative medicine.

The glucocorticoid receptor (GR) binds as a homodimer to genomic response elements, which have particular sequence and shape characteristics. Here we show that the nucleotides directly flanking the core-binding site, differ depending on the strength of GR-dependent activation of nearby genes. Our study indicates that these flanking nucleotides change the three-dimensional structure of the DNA-binding site, the DNA-binding domain of GR and the quaternary structure of the dimeric complex. Functional studies in a defined genomic context show that sequence-induced changes in GR activity cannot be explained by differences in GR occupancy. Rather, mutating the dimerization interface mitigates DNA-induced changes in both activity and structure, arguing for a role of DNA-induced structural changes in modulating GR activity. Together, our study shows that DNA sequence identity of genomic binding sites modulates GR activity downstream of binding, which may play a role in achieving regulatory specificity towards individual target genes.

Transcription factors (TFs) have to find their binding sites, which are distributed throughout the genome. Facilitated diffusion is currently the most widely accepted model for this search process. Based on this model the TF alternates between one-dimensional sliding along the DNA, and three-dimensional bulk diffusion. In this view, the non-specific associations between the proteins and the DNA play a major role in the search dynamics. However, little is known about how the DNA properties around the motif contribute to the search. Accumulating evidence showing that TF binding sites are embedded within a unique environment, specific to each TF, leads to the hypothesis that the search process is facilitated by favorable DNA features that help to improve the search efficiency. Here, we review the field and present the hypothesis that TF-DNA recognition is dictated not only by the motif, but is also influenced by the environment in which the motif resides.

Transcriptional regulation requires the binding of transcription factors (TFs) to short sequence-specific DNA motifs, usually located at the gene regulatory regions. Interestingly, based on a vast amount of data accumulated from genomic assays, it has been shown that only a small fraction of all potential binding sites containing the consensus motif of a given TF actually bind the protein. Recent in vitro binding assays, which exclude the effects of the cellular environment, also demonstrate selective TF binding. An intriguing conjecture is that the surroundings of cognate binding sites have unique characteristics that distinguish them from other sequences containing a similar motif that are not bound by the TF. To test this hypothesis, we conducted a comprehensive analysis of the sequence and DNA shape features surrounding the core-binding sites of 239 and 56 TFs extracted from in vitro HT-SELEX binding assays and in vivo ChIP-seq data, respectively. Comparing the nucleotide content of the regions around the TF-bound sites to the counterpart unbound regions containing the same consensus motifs revealed significant differences that extend far beyond the core-binding site. Specifically, the environment of the bound motifs demonstrated unique sequence compositions, DNA shape features, and overall high similarity to the core-binding motif. Notably, the regions around the binding sites of TFs that belong to the same TF families exhibited similar features, with high agreement between the in vitro and in vivo data sets. We propose that these unique features assist in guiding TFs to their cognate binding sites.

Protein–DNA recognition is a critical component of gene regulatory processes but the underlying molecular mechanisms are not yet completely understood. Whereas the DNA binding preferences of transcription factors (TFs) are commonly described using nucleotide sequences, the 3D DNA structure is recognized by proteins and is crucial for achieving binding specificity. However, the ability to analyze DNA shape in a high-throughput manner made it only recently feasible to integrate structural information into studies of protein–DNA binding. Here we focused on the homeodomain family of TFs and analyzed the DNA shape of thousands of their DNA binding sites, investigating the covariation between the protein sequence and the sequence and shape of their DNA targets. We found distinct homeodomain regions that were more correlated with either the nucleotide sequence or the DNA shape of their preferred binding sites, demonstrating different readout mechanisms through which homeodomains attain DNA binding specificity. We identified specific homeodomain residues that likely play key roles in DNA recognition via shape readout. Finally, we showed that adding DNA shape information when characterizing binding sites improved the prediction accuracy of homeodomain binding specificities. Taken together, our findings indicate that DNA shape information can generally provide new mechanistic insights into TF binding.

It is well established that transcription factors (TFs) identify their binding sites through direct contacts with unique chemical groups of the base pairs mainly in the major groove. A second type of mechanism, which has been relatively less studied, is the readout of the DNA shape, in which the protein recognizes the three-dimensional DNA structure (Rohs et al., 2010). Here, we focused on the homeodomain family of TFs and analyzed the DNA shape of thousands of sequences in order to study the correlation between the amino acid sequence of homeodomains and the nucleotide sequence and the shape of their DNA binding sites. We have found regions in the homeodomains that are significantly correlated with the sequence or with the shape of their preferred binding sites, demonstrating the role of different homeodomain regions in attaining binding specificity through different modes of recognition. Next, we predicted specific residues in homeodomains which likely play an important role in DNA recognition through DNA shape attributes. Furthermore, we show that adding DNA shape information to the characterization of TF binding sites can improve predictions of homeodomain binding specificity. Finally, our work indicates that DNA shape information can provide new mechanistic insights into TF binding.

Protein-DNA binding is mediated by the recognition of the chemical signatures of the DNA bases and the 3D shape of the DNA molecule. Because DNA shape is a consequence of sequence, it is difficult to dissociate these modes of recognition. Here, we tease them apart in the context of Hox-DNA binding by mutating residues that, in a co-crystal structure, only recognize DNA shape. Complexes made with these mutants lose the preference to bind sequences with specific DNA shape features. Introducing shape-recognizing residues from one Hox protein to another swapped binding specificities in vitro and gene regulation in vivo. Statistical machine learning revealed that the accuracy of binding specificity predictions improves by adding shape features to a model that only depends on sequence, and feature selection identified shape features important for recognition. Thus, shape readout is a direct and independent component of binding site selection by Hox proteins.

DNA sequence is a major determinant of the binding specificity of transcription factors (TFs) for their genomic targets. However, eukaryotic cells often express, at the same time, TFs with highly similar DNA binding motifs but distinct in vivo targets. Currently, it is not well understood how TFs with seemingly identical DNA motifs achieve unique specificities in vivo. Here, we used custom protein-binding microarrays to analyze TF specificity for putative binding sites in their genomic sequence context. Using yeast TFs Cbf1 and Tye7 as our case studies, we found that binding sites of these bHLH TFs (i.e., E-boxes) are bound differently in vitro and in vivo, depending on their genomic context. Computational analyses suggest that nucleotides outside E-box binding sites contribute to specificity by influencing the three-dimensional structure of DNA binding sites. Thus, the local shape of target sites might play a widespread role in achieving regulatory specificity within TF families.

Transcription factor binding sites (TFBSs) are most commonly characterized by the nucleotide preferences at each position of the DNA target. Whereas these sequence motifs are quite accurate descriptions of DNA binding specificities of transcription factors (TFs), proteins recognize DNA as a three-dimensional object. DNA structural features refine the description of TF binding specificities and provide mechanistic insights into protein–DNA recognition. Existing motif databases contain extensive nucleotide sequences identified in binding experiments based on their selection by a TF. To utilize DNA shape information when analysing the DNA binding specificities of TFs, we developed a new tool, the TFBSshape database (available at http://rohslab.cmb.usc.edu/TFBSshape/), for calculating DNA structural features from nucleotide sequences provided by motif databases. The TFBSshape database can be used to generate heat maps and quantitative data for DNA structural features (i.e., minor groove width, roll, propeller twist and helix twist) for 739 TF datasets from 23 different species derived from the motif databases JASPAR and UniPROBE. As demonstrated for the basic helix-loop-helix and homeodomain TF families, our TFBSshape database can be used to compare, qualitatively and quantitatively, the DNA binding specificities of closely related TFs and, thus, uncover differential DNA binding specificities that are not apparent from nucleotide sequence alone.

We present a method and web server for predicting DNA structural features in a high-throughput (HT) manner for massive sequence data. This approach provides the framework for the integration of DNA sequence and shape analyses in genome-wide studies. The HT methodology uses a sliding-window approach to mine DNA structural information obtained from Monte Carlo simulations. It requires only nucleotide sequence as input and instantly predicts multiple structural features of DNA (minor groove width, roll, propeller twist and helix twist). The results of rigorous validations of the HT predictions based on DNA structures solved by X-ray crystallography and NMR spectroscopy, hydroxyl radical cleavage data, statistical analysis and cross-validation, and molecular dynamics simulations provide strong confidence in this approach. The DNAshape web server is freely available at http://rohslab.cmb.usc.edu/DNAshape/.