Imprinted X chromosome inactivation at the gamete-to-embryo transition.

In mammals, dosage compensation involves two parallel processes: (1) X inactivation, which equalizes X chromosome dosage between males and females, and (2) X hyperactivation, which upregulates the active X for X-autosome balance. The field currently favors models whereby dosage compensation initiates "de novo" during mouse development. Here, we develop "So-Smart-seq" to revisit the question and interrogate a comprehensive transcriptome including noncoding genes and repeats in mice. Intriguingly, de novo silencing pertains only to a subset of Xp genes. Evolutionarily older genes and repetitive elements demonstrate constitutive Xp silencing, adopt distinct signatures, and do not require Xist to initiate silencing. We trace Xp silencing backward in developmental time to meiotic sex chromosome inactivation in the male germ line and observe that Xm hyperactivation is timed to Xp silencing on a gene-by-gene basis. Thus, during the gamete-to-embryo transition, older Xp genes are transmitted in a "pre-inactivated" state. These findings have implications for the evolution of imprinting.

Targeting Xist with compounds that disrupt RNA structure and X inactivation.

Although more than 98% of the human genome is non-coding1, nearly all of the drugs on the market target one of about 700 disease-related proteins. The historical reluctance to invest in non-coding RNA stems partly from requirements for drug targets to adopt a single stable conformation2. Most RNAs can adopt several conformations of similar stabilities. RNA structures also remain challenging to determine3. Nonetheless, an increasing number of diseases are now being attributed to non-coding RNA4 and the ability to target them would vastly expand the chemical space for drug development. Here we devise a screening strategy and identify small molecules that bind the non-coding RNA prototype Xist5. The X1 compound has drug-like properties and binds specifically the RepA motif6 of Xist in vitro and in vivo. Small-angle X-ray scattering analysis reveals that RepA can adopt multiple conformations but favours one structure in solution. X1 binding reduces the conformational space of RepA, displaces cognate interacting protein factors (PRC2 and SPEN), suppresses histone H3K27 trimethylation, and blocks initiation of X-chromosome inactivation. X1 inhibits cell differentiation and growth in a female-specific manner. Thus, RNA can be systematically targeted by drug-like compounds that disrupt RNA structure and epigenetic function.

Jpx RNA regulates CTCF anchor site selection and formation of chromosome loops.

Chromosome loops shift dynamically during development, homeostasis, and disease. CCCTC-binding factor (CTCF) is known to anchor loops and construct 3D genomes, but how anchor sites are selected is not yet understood. Here, we unveil Jpx RNA as a determinant of anchor selectivity. Jpx RNA targets thousands of genomic sites, preferentially binding promoters of active genes. Depleting Jpx RNA causes ectopic CTCF binding, massive shifts in chromosome looping, and downregulation of >700 Jpx target genes. Without Jpx, thousands of lost loops are replaced by de novo loops anchored by ectopic CTCF sites. Although Jpx controls CTCF binding on a genome-wide basis, it acts selectively at the subset of developmentally sensitive CTCF sites. Specifically, Jpx targets low-affinity CTCF motifs and displaces CTCF protein through competitive inhibition. We conclude that Jpx acts as a CTCF release factor and shapes the 3D genome by regulating anchor site usage.

Denaturing cross-linking immunoprecipitation to identify footprints for RNA-binding proteins.

The isolation of protein-RNA complexes in the "denaturing cross-linked RNA immunoprecipitation" (dCLIP) protocol is based on biotin-tagging proteins of interest, UV cross-linking RNA to protein in vivo, RNase protection assay, and isolating RNA-protein complexes under denaturing conditions over a streptavidin column. Insofar as conventional antibody-based CLIP assays have been challenging to apply to Polycomb complexes, dCLIP has been applied successfully and yields small RNA footprints from which de novo motif analysis can be performed to identify RNA binding motifs. For complete details on the use and execution of this protocol, please refer to Rosenberg et al. (2017).

A disproportionate impact of G9a methyltransferase deficiency on the X chromosome.

G9a is a histone methyltransferase responsible for the dimethylation of histone H3 at lysine 9 (H3K9me2). G9a plays key roles in transcriptional silencing of developmentally regulated genes, but its role in X-chromosome inactivation (XCI) has been under debate. Here, we uncover a female-specific function of G9a and demonstrate that deleting G9a has a disproportionate impact on the X chromosome relative to the rest of the genome. G9a deficiency causes a failure of XCI and female-specific hypersensitivity to drug inhibition of H3K9me2. We show that G9a interacts with Tsix and Xist RNAs, and that competitive inhibition of the G9a-RNA interaction recapitulates the XCI defect. During XCI, Xist recruits G9a to silence X-linked genes on the future inactive X. In parallel on the future Xa, Tsix recruits G9a to silence Xist in cis Thus, RNA tethers G9a for allele-specific targeting of the H3K9me2 modification and the G9a-RNA interaction is essential for XCI.

Motif-driven interactions between RNA and PRC2 are rheostats that regulate transcription elongation.

Although polycomb repressive complex 2 (PRC2) is now recognized as an RNA-binding complex, the full range of binding motifs and why PRC2-RNA complexes often associate with active genes have not been elucidated. Here, we identify high-affinity RNA motifs whose mutations weaken PRC2 binding and attenuate its repressive function in mouse embryonic stem cells. Interactions occur at promoter-proximal regions and frequently coincide with pausing of RNA polymerase II (POL-II). Surprisingly, while PRC2-associated nascent transcripts are highly expressed, ablating PRC2 further upregulates expression via loss of pausing and enhanced transcription elongation. Thus, PRC2-nascent RNA complexes operate as rheostats to fine-tune transcription by regulating transitions between pausing and elongation, explaining why PRC2-RNA complexes frequently occur within active genes. Nascent RNA also targets PRC2 in cis and downregulates neighboring genes. We propose a unifying model in which RNA specifically recruits PRC2 to repress genes through POL-II pausing and, more classically, trimethylation of histone H3 at Lys27.

Denaturing CLIP, dCLIP, Pipeline Identifies Discrete RNA Footprints on Chromatin-Associated Proteins and Reveals that CBX7 Targets 3' UTRs to Regulate mRNA Expression.

Interaction networks between chromatin complexes and long noncoding RNAs have become a recurrent theme in epigenetic regulation. However, technical limitations have precluded identification of RNA binding motifs for chromatin-associated proteins. Here, we add a denaturation step to UV-crosslink RNA immunoprecipitation (dCLIP) and apply dCLIP to mouse and human chromobox homolog 7 (CBX7), an RNA binding subunit of Polycomb repressive complex 1 (PRC1). In both species, CBX7 predominantly binds 3' UTRs of messenger RNAs. CBX7 binds with a median RNA "footprint" of 171-183 nucleotides, the small size of which facilitates motif identification by bioinformatics. We find four families of consensus RNA motifs in mouse, and independent analysis of human CBX7 dCLIP data identifies similar motifs. Their mutation abolishes CBX7 binding in vitro. Pharmacological intervention with antisense oligonucleotides paradoxically increases CBX7 binding and enhances gene expression. These data support the utility of dCLIP and reveal an unexpected functional interaction between CBX7 and the 3' UTRs of mRNA.

TERRA RNA Antagonizes ATRX and Protects Telomeres.

Through an integration of genomic and proteomic approaches to advance understanding of long noncoding RNAs, we investigate the function of the telomeric transcript, TERRA. By identifying thousands of TERRA target sites in the mouse genome, we demonstrate that TERRA can bind both in cis to telomeres and in trans to genic targets. We then define a large network of interacting proteins, including epigenetic factors, telomeric proteins, and the RNA helicase, ATRX. TERRA and ATRX share hundreds of target genes and are functionally antagonistic at these loci: whereas TERRA activates, ATRX represses gene expression. At telomeres, TERRA competes with telomeric DNA for ATRX binding, suppresses ATRX localization, and ensures telomeric stability. Depleting TERRA increases telomerase activity and induces telomeric pathologies, including formation of telomere-induced DNA damage foci and loss or duplication of telomeric sequences. We conclude that TERRA functions as an epigenomic modulator in trans and as an essential regulator of telomeres in cis.

PAR-TERRA directs homologous sex chromosome pairing.

In mammals, homologous chromosomes rarely pair outside meiosis. One exception is the X chromosome, which transiently pairs during X-chromosome inactivation (XCI). How two chromosomes find each other in 3D space is not known. Here, we reveal a required interaction between the X-inactivation center (Xic) and the telomere in mouse embryonic stem (ES) cells. The subtelomeric, pseudoautosomal regions (PARs) of the two sex chromosomes (X and Y) also undergo pairing in both female and male cells. PARs transcribe a class of telomeric RNA, dubbed PAR-TERRA, which accounts for a vast majority of all TERRA transcripts. PAR-TERRA binds throughout the genome, including to the PAR and Xic. During X-chromosome pairing, PAR-TERRA anchors the Xic to the PAR, creating a 'tetrad' of pairwise homologous interactions (Xic-Xic, PAR-PAR, and Xic-PAR). Xic pairing occurs within the tetrad. Depleting PAR-TERRA abrogates pairing and blocks initiation of XCI, whereas autosomal PAR-TERRA induces ectopic pairing. We propose a 'constrained diffusion model' in which PAR-TERRA creates an interaction hub to guide Xic homology searching during XCI.

A high-throughput small molecule screen identifies synergism between DNA methylation and Aurora kinase pathways for X reactivation.

X-chromosome inactivation is a mechanism of dosage compensation in which one of the two X chromosomes in female mammals is transcriptionally silenced. Once established, silencing of the inactive X (Xi) is robust and difficult to reverse pharmacologically. However, the Xi is a reservoir of >1,000 functional genes that could be potentially tapped to treat X-linked disease. To identify compounds that could reactivate the Xi, here we screened ∼367,000 small molecules in an automated high-content screen using an Xi-linked GFP reporter in mouse fibroblasts. Given the robust nature of silencing, we sensitized the screen by "priming" cells with the DNA methyltransferase inhibitor, 5-aza-2'-deoxycytidine (5azadC). Compounds that elicited GFP activity include VX680, MLN8237, and 5azadC, which are known to target the Aurora kinase and DNA methylation pathways. We demonstrate that the combinations of VX680 and 5azadC, as well as MLN8237 and 5azadC, synergistically up-regulate genes on the Xi. Thus, our work identifies a synergism between the DNA methylation and Aurora kinase pathways as being one of interest for possible pharmacological reactivation of the Xi.

Chromosomes. A comprehensive Xist interactome reveals cohesin repulsion and an RNA-directed chromosome conformation.

The inactive X chromosome (Xi) serves as a model to understand gene silencing on a global scale. Here, we perform "identification of direct RNA interacting proteins" (iDRiP) to isolate a comprehensive protein interactome for Xist, an RNA required for Xi silencing. We discover multiple classes of interactors-including cohesins, condensins, topoisomerases, RNA helicases, chromatin remodelers, and modifiers-that synergistically repress Xi transcription. Inhibiting two or three interactors destabilizes silencing. Although Xist attracts some interactors, it repels architectural factors. Xist evicts cohesins from the Xi and directs an Xi-specific chromosome conformation. Upon deleting Xist, the Xi acquires the cohesin-binding and chromosomal architecture of the active X. Our study unveils many layers of Xi repression and demonstrates a central role for RNA in the topological organization of mammalian chromosomes.

Locus-specific targeting to the X chromosome revealed by the RNA interactome of CTCF.

CTCF is a master regulator that plays important roles in genome architecture and gene expression. How CTCF is recruited in a locus-specific manner is not fully understood. Evidence from epigenetic processes, such as X chromosome inactivation (XCI), indicates that CTCF associates functionally with RNA. Using genome-wide approaches to investigate the relationship between its RNA interactome and epigenomic landscape, here we report that CTCF binds thousands of transcripts in mouse embryonic stem cells, many in close proximity to CTCF's genomic binding sites. CTCF is a specific and high-affinity RNA-binding protein (Kd < 1 nM). During XCI, CTCF differentially binds the active and inactive X chromosomes and interacts directly with Tsix, Xite, and Xist RNAs. Tsix and Xite RNAs target CTCF to the X inactivation center, thereby inducing homologous X chromosome pairing. Our work elucidates one mechanism by which CTCF is recruited in a locus-specific manner and implicates CTCF-RNA interactions in long-range chromosomal interactions.

High-resolution Xist binding maps reveal two-step spreading during X-chromosome inactivation.

The Xist long noncoding RNA (lncRNA) is essential for X-chromosome inactivation (XCI), the process by which mammals compensate for unequal numbers of sex chromosomes. During XCI, Xist coats the future inactive X chromosome (Xi) and recruits Polycomb repressive complex 2 (PRC2) to the X-inactivation centre (Xic). How Xist spreads silencing on a 150-megabases scale is unclear. Here we generate high-resolution maps of Xist binding on the X chromosome across a developmental time course using CHART-seq. In female cells undergoing XCI de novo, Xist follows a two-step mechanism, initially targeting gene-rich islands before spreading to intervening gene-poor domains. Xist is depleted from genes that escape XCI but may concentrate near escapee boundaries. Xist binding is linearly proportional to PRC2 density and H3 lysine 27 trimethylation (H3K27me3), indicating co-migration of Xist and PRC2. Interestingly, when Xist is acutely stripped off from the Xi in post-XCI cells, Xist recovers quickly within both gene-rich and gene-poor domains on a timescale of hours instead of days, indicating a previously primed Xi chromatin state. We conclude that Xist spreading takes distinct stage-specific forms. During initial establishment, Xist follows a two-step mechanism, but during maintenance, Xist spreads rapidly to both gene-rich and gene-poor regions.

Isoform divergence of the filamin family of proteins.

The vertebrate filamin family (A, B, and C) is part of the spectrin family of actin cross-linking proteins. Family members share high sequence similarity (>64%) and have both common and isoform-distinct functionalities. To identify the basis for isoform-specific functionality, we perform an evolutionary trace of chordate filamin at the granularity of single residues. Our trace methodology is constrained to focus on neofunctionality by requiring that one isoform remain the ancestral type, whereas at least one isoform has an accepted mutation. We call divergence meeting these characteristics "class-distinctive." To obtain a temporal and spatial context for class-distinctive residues, we derive an all-atom model of full-length filamin A by homology modeling and joining individual domains. We map onto our model both conserved and class-distinctive residues along with the period (Teleostei, Amphibian, and Mammalian) in which they diverged. Our phylogenetic analysis suggests that filamins diverged from a common ancestral gene between urochordate and vertebrate lineages. Filamins also diverged the most just after gene duplication, in the Teleostei period, with filamin C remaining closest to ancestral filamin. At the residue level, domains with well-characterized interfaces, IgFLN 17 and IgFLN 21 (immunoglobulin, Ig), have diverged in potentially critical residues in their adhesion protein-binding interfaces, signifying that isoforms may bind or regulate ligand binding differentially. Similarly, isoform divergence in a region associated with F actin-binding regulation suggests that isoforms differentially regulate F-actin binding. In addition, we observe some class-distinctive residues in the vicinity of missense mutations that cause filamin A and B-associated skeletal disorders. Our analysis, utilizing both spatial and temporal granularity, has identified potentially important residues responsible for vertebrate filamin isoform-specific divergence-significantly in regions where few binding partners have been discovered to date- and suggests yet to be discovered filamin-binding partners and isoform-specific differential regulation with these binding partners.

N-terminal strands of filamin Ig domains act as a conformational switch under biological forces.

Conformational changes of filamin A under stress have been postulated to play crucial roles in signaling pathways of cell responses. Direct observation of conformational changes under stress is beyond the resolution of current experimental techniques. On the other hand, computational studies are mainly limited to either traditional molecular dynamics simulations of short durations and high forces or simulations of simplified models. Here we perform all-atom discrete molecular dynamics (DMD) simulations to study thermally and force-induced unfolding of filamin A. The high conformational sampling efficiency of DMD allows us to observe force-induced unfolding of filamin A Ig domains under physiological forces. The computationally identified critical unfolding forces agree well with experimental measurements. Despite a large heterogeneity in the population of force-induced intermediate states, we find a common initial unfolding intermediate in all the Ig domains of filamin, where the N-terminal strand unfolds. We also study the thermal unfolding of several filamin Ig-like domains. We find that thermally induced unfolding features an early-stage intermediate state similar to the one observed in force-induced unfolding and characterized by the N-terminal strand being unfurled. We propose that the N-terminal strand may act as a conformational switch that unfolds under physiological forces leading to exposure of cryptic binding sites, removal of native binding sites, and modulating the quaternary structure of domains.

Improved protein identification through the use of unstained gels.

In this work, a method for improved protein identification of low-abundance proteins using unstained gels, in combination with robotics and matrix-assisted laser desorption/ionization tandem mass spectrometry, has been developed and evaluated. Omitting the silver-staining process resulted in increased protein identification scores, an increase in the number of peptides observed in the MALDI mass spectrum, and improved quality of the tandem mass spectrometry data.