Enhancer switching in cell lineage priming is linked to eRNA, Brg1's AT-hook, and SWI/SNF recruitment.

RNA transcribed from enhancers, i.e., eRNA, has been suggested to directly activate transcription by recruiting transcription factors and co-activators. Although there have been specific examples of eRNA functioning in this way, it is not clear how general this may be. We find that the AT-hook of SWI/SNF preferentially binds RNA and, as part of the esBAF complex, associates with eRNA transcribed from intronic and intergenic regions. Our data suggest that SWI/SNF is globally recruited in cis by eRNA to cell-type-specific enhancers, representative of two distinct stages that mimic early mammalian development, and not at enhancers that are shared between the two stages. In this manner, SWI/SNF facilitates recruitment and/or activation of MLL3/4, p300/CBP, and Mediator to stage-specific enhancers and super-enhancers that regulate the transcription of metabolic and cell lineage priming-related genes. These findings highlight a connection between ATP-dependent chromatin remodeling and eRNA in cell identity and typical- and super-enhancer activation.

Systematic discovery of Xist RNA binding proteins.

Noncoding RNAs (ncRNAs) function with associated proteins to effect complex structural and regulatory outcomes. To reveal the composition and dynamics of specific noncoding RNA-protein complexes (RNPs) in vivo, we developed comprehensive identification of RNA binding proteins by mass spectrometry (ChIRP-MS). ChIRP-MS analysis of four ncRNAs captures key protein interactors, including a U1-specific link to the 3' RNA processing machinery. Xist, an essential lncRNA for X chromosome inactivation (XCI), interacts with 81 proteins from chromatin modification, nuclear matrix, and RNA remodeling pathways. The Xist RNA-protein particle assembles in two steps coupled with the transition from pluripotency to differentiation. Specific interactors include HnrnpK, which participates in Xist-mediated gene silencing and histone modifications but not Xist localization, and Drosophila Split ends homolog Spen, which interacts via the A-repeat domain of Xist and is required for gene silencing. Thus, Xist lncRNA engages with proteins in a modular and developmentally controlled manner to coordinate chromatin spreading and silencing.

Oncogenic transcription factors and neogenes: New opportunities for cancer immunotherapy?

Vibert et al. (2022) demonstrate that oncogenic transcription factor fusion proteins bind otherwise silent genomic regions, producing RNAs that can be spliced, exported, and translated. These "neogenes" represent possible targets for immunotherapy and may even be universal byproducts of altered transcription in cancer.

Site-specific silencing of regulatory elements as a mechanism of X inactivation.

The inactive X chromosome's (Xi) physical territory is microscopically devoid of transcriptional hallmarks and enriched in silencing-associated modifications. How these microscopic signatures relate to specific Xi sequences is unknown. Therefore, we profiled Xi gene expression and chromatin states at high resolution via allele-specific sequencing in mouse trophoblast stem cells. Most notably, X-inactivated transcription start sites harbored distinct epigenetic signatures relative to surrounding Xi DNA. These sites displayed H3-lysine27-trimethylation enrichment and DNaseI hypersensitivity, similar to autosomal Polycomb targets, yet excluded Pol II and other transcriptional hallmarks, similar to nontranscribed genes. CTCF bound X-inactivated and escaping genes, irrespective of measured chromatin boundaries. Escape from X inactivation occurred within, and X inactivation was maintained exterior to, the area encompassed by Xist in cells subject to imprinted and random X inactivation. The data support a model whereby inactivation of specific regulatory elements, rather than a simple chromosome-wide separation from transcription machinery, governs gene silencing over the Xi.

lncRNA-Induced Spread of Polycomb Controlled by Genome Architecture, RNA Abundance, and CpG Island DNA.

Long noncoding RNAs (lncRNAs) cause Polycomb repressive complexes (PRCs) to spread over broad regions of the mammalian genome. We report that in mouse trophoblast stem cells, the Airn and Kcnq1ot1 lncRNAs induce PRC-dependent chromatin modifications over multi-megabase domains. Throughout the Airn-targeted domain, the extent of PRC-dependent modification correlated with intra-nuclear distance to the Airn locus, preexisting genome architecture, and the abundance of Airn itself. Specific CpG islands (CGIs) displayed characteristics indicating that they nucleate the spread of PRCs upon exposure to Airn. Chromatin environments surrounding Xist, Airn, and Kcnq1ot1 suggest common mechanisms of PRC engagement and spreading. Our data indicate that lncRNA potency can be tightly linked to lncRNA abundance and that within lncRNA-targeted domains, PRCs are recruited to CGIs via lncRNA-independent mechanisms. We propose that CGIs that autonomously recruit PRCs interact with lncRNAs and their associated proteins through three-dimensional space to nucleate the spread of PRCs in lncRNA-targeted domains.

SAFB associates with nascent RNAs and can promote gene expression in mouse embryonic stem cells.

Scaffold attachment factor B (SAFB) is a conserved RNA-binding protein that is essential for early mammalian development. However, the functions of SAFB in mouse embryonic stem cells (ESCs) have not been characterized. Using RNA immunoprecipitation followed by RNA-seq (RIP-seq), we examined the RNAs associated with SAFB in wild-type and SAFB/SAFB2 double-knockout ESCs. SAFB predominantly associated with introns of protein-coding genes through purine-rich motifs. The transcript most enriched in SAFB association was the lncRNA Malat1, which also contains a purine-rich region in its 5' end. Knockout of SAFB/SAFB2 led to differential expression of approximately 1000 genes associated with multiple biological processes, including apoptosis, cell division, and cell migration. Knockout of SAFB/SAFB2 also led to splicing changes in a set of genes that were largely distinct from those that exhibited changes in expression level. The spliced and nascent transcripts of many genes whose expression levels were positively regulated by SAFB also associated with high levels of SAFB, implying that SAFB binding promotes their expression. Reintroduction of SAFB into double-knockout cells restored gene expression toward wild-type levels, an effect again observable at the level of spliced and nascent transcripts. Proteomics analysis revealed a significant enrichment of nuclear speckle-associated and RS domain-containing proteins among SAFB interactors. Neither Xist nor Polycomb functions were dramatically altered in SAFB/2 knockout ESCs. Our findings suggest that among other potential functions in ESCs, SAFB promotes the expression of certain genes through its ability to bind nascent RNA.

Connecting microRNA genes to the core transcriptional regulatory circuitry of embryonic stem cells.

MicroRNAs (miRNAs) are crucial for normal embryonic stem (ES) cell self-renewal and cellular differentiation, but how miRNA gene expression is controlled by the key transcriptional regulators of ES cells has not been established. We describe here the transcriptional regulatory circuitry of ES cells that incorporates protein-coding and miRNA genes based on high-resolution ChIP-seq data, systematic identification of miRNA promoters, and quantitative sequencing of short transcripts in multiple cell types. We find that the key ES cell transcription factors are associated with promoters for miRNAs that are preferentially expressed in ES cells and with promoters for a set of silent miRNA genes. This silent set of miRNA genes is co-occupied by Polycomb group proteins in ES cells and shows tissue-specific expression in differentiated cells. These data reveal how key ES cell transcription factors promote the ES cell miRNA expression program and integrate miRNAs into the regulatory circuitry controlling ES cell identity.

Elements at the 5' end of Xist harbor SPEN-independent transcriptional antiterminator activity.

The Xist lncRNA requires Repeat A, a conserved RNA element located in its 5' end, to induce gene silencing during X-chromosome inactivation. Intriguingly, Repeat A is also required for production of Xist. While silencing by Repeat A requires the protein SPEN, how Repeat A promotes Xist production remains unclear. We report that in mouse embryonic stem cells, expression of a transgene comprising the first two kilobases of Xist (Xist-2kb) causes transcriptional readthrough of downstream polyadenylation sequences. Readthrough required Repeat A and the ∼750 nucleotides downstream, did not require SPEN, and was attenuated by splicing. Despite associating with SPEN and chromatin, Xist-2kb did not robustly silence transcription, whereas a 5.5-kb Xist transgene robustly silenced transcription and read through its polyadenylation sequence. Longer, spliced Xist transgenes also induced robust silencing yet terminated efficiently. Thus, in contexts examined here, Xist requires sequence elements beyond its first two kilobases to robustly silence transcription, and the 5' end of Xist harbors SPEN-independent transcriptional antiterminator activity that can repress proximal cleavage and polyadenylation. In endogenous contexts, this antiterminator activity may help produce full-length Xist RNA while rendering the Xist locus resistant to silencing by the same repressive complexes that the lncRNA recruits to other genes.

Nonlinear sequence similarity between the Xist and Rsx long noncoding RNAs suggests shared functions of tandem repeat domains.

The marsupial inactive X chromosome expresses a long noncoding RNA (lncRNA) called Rsx that has been proposed to be the functional analog of eutherian Xist Despite the possibility that Xist and Rsx encode related functions, the two lncRNAs harbor no linear sequence similarity. However, both lncRNAs harbor domains of tandemly repeated sequence. In Xist, these repeat domains are known to be critical for function. Using k-mer based comparison, we show that the repeat domains of Xist and Rsx unexpectedly partition into two major clusters that each harbor substantial levels of nonlinear sequence similarity. Xist Repeats B, C, and D were most similar to each other and to Rsx Repeat 1, whereas Xist Repeats A and E were most similar to each other and to Rsx Repeats 2, 3, and 4. Similarities at the level of k-mers corresponded to domain-specific enrichment of protein-binding motifs. Within individual domains, protein-binding motifs were often enriched to extreme levels. Our data support the hypothesis that Xist and Rsx encode similar functions through different spatial arrangements of functionally analogous protein-binding domains. We propose that the two clusters of repeat domains in Xist and Rsx function in part to cooperatively recruit PRC1 and PRC2 to chromatin. The physical manner in which these domains engage with protein cofactors may be just as critical to the function of the domains as the protein cofactors themselves. The general approaches we outline in this report should prove useful in the study of any set of RNAs.

A piggyBac-based toolkit for inducible genome editing in mammalian cells.

We describe the development and application of a novel series of vectors that facilitate CRISPR-Cas9-mediated genome editing in mammalian cells, which we call CRISPR-Bac. CRISPR-Bac leverages the piggyBac transposon to randomly insert CRISPR-Cas9 components into mammalian genomes. In CRISPR-Bac, a single piggyBac cargo vector containing a doxycycline-inducible Cas9 or catalytically dead Cas9 (dCas9) variant and a gene conferring resistance to Hygromycin B is cotransfected with a plasmid expressing the piggyBac transposase. A second cargo vector, expressing a single-guide RNA (sgRNA) of interest, the reverse-tetracycline TransActivator (rtTA), and a gene conferring resistance to G418, is also cotransfected. Subsequent selection on Hygromycin B and G418 generates polyclonal cell populations that stably express Cas9, rtTA, and the sgRNA(s) of interest. We show that CRISPR-Bac can be used to knock down proteins of interest, to create targeted genetic deletions with high efficiency, and to activate or repress transcription of protein-coding genes and an imprinted long noncoding RNA. The ratio of sgRNA-to-Cas9-to-transposase can be adjusted in transfections to alter the average number of cargo insertions into the genome. sgRNAs targeting multiple genes can be inserted in a single transfection. CRISPR-Bac is a versatile platform for genome editing that simplifies the generation of mammalian cells that stably express the CRISPR-Cas9 machinery.

RETRACTED: A 5' fragment of Xist can sequester RNA produced from adjacent genes on chromatin.

Xist requires Repeat-A, a protein-binding module in its first two kilobases (2kb), to repress transcription. We report that when expressed as a standalone transcript in mouse embryonic stem cells (ESCs), the first 2kb of Xist (Xist-2kb) does not induce transcriptional silencing. Instead, Xist-2kb sequesters RNA produced from adjacent genes on chromatin. Sequestration does not spread beyond adjacent genes, requires the same sequence elements in Repeat-A that full-length Xist requires to repress transcription and can be induced by lncRNAs with similar sequence composition to Xist-2kb. We did not detect sequestration by full-length Xist, but we did detect it by mutant forms of Xist with attenuated transcriptional silencing capability. Xist-2kb associated with SPEN, a Repeat-A binding protein required for Xist-induced transcriptional silencing, but SPEN was not necessary for sequestration. Thus, when expressed in mouse ESCs, a 5' fragment of Xist that contains Repeat-A sequesters RNA from adjacent genes on chromatin and associates with the silencing factor SPEN, but it does not induce transcriptional silencing. Instead, Xist-induced transcriptional silencing requires synergy between Repeat-A and additional sequence elements in Xist. We propose that sequestration is mechanistically related to the Repeat-A dependent stabilization and tethering of Xist near actively transcribed regions of chromatin.

Topoisomerases facilitate transcription of long genes linked to autism.

Topoisomerases are expressed throughout the developing and adult brain and are mutated in some individuals with autism spectrum disorder (ASD). However, how topoisomerases are mechanistically connected to ASD is unknown. Here we find that topotecan, a topoisomerase 1 (TOP1) inhibitor, dose-dependently reduces the expression of extremely long genes in mouse and human neurons, including nearly all genes that are longer than 200 kilobases. Expression of long genes is also reduced after knockdown of Top1 or Top2b in neurons, highlighting that both enzymes are required for full expression of long genes. By mapping RNA polymerase II density genome-wide in neurons, we found that this length-dependent effect on gene expression was due to impaired transcription elongation. Interestingly, many high-confidence ASD candidate genes are exceptionally long and were reduced in expression after TOP1 inhibition. Our findings suggest that chemicals and genetic mutations that impair topoisomerases could commonly contribute to ASD and other neurodevelopmental disorders.

SHAPE reveals transcript-wide interactions, complex structural domains, and protein interactions across the Xist lncRNA in living cells.

The 18-kb Xist long noncoding RNA (lncRNA) is essential for X-chromosome inactivation during female eutherian mammalian development. Global structural architecture, cell-induced conformational changes, and protein-RNA interactions within Xist are poorly understood. We used selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP) to examine these features of Xist at single-nucleotide resolution both in living cells and ex vivo. The Xist RNA forms complex well-defined secondary structure domains and the cellular environment strongly modulates the RNA structure, via motifs spanning one-half of all Xist nucleotides. The Xist RNA structure modulates protein interactions in cells via multiple mechanisms. For example, repeat-containing elements adopt accessible and dynamic structures that function as landing pads for protein cofactors. Structured RNA motifs create interaction domains for specific proteins and also sequester other motifs, such that only a subset of potential binding sites forms stable interactions. This work creates a broad quantitative framework for understanding structure-function interrelationships for Xist and other lncRNAs in cells.

fourSig: a method for determining chromosomal interactions in 4C-Seq data.

The ability to correlate chromosome conformation and gene expression gives a great deal of information regarding the strategies used by a cell to properly regulate gene activity. 4C-Seq is a relatively new and increasingly popular technology where the set of genomic interactions generated by a single point in the genome can be determined. 4C-Seq experiments generate large, complicated data sets and it is imperative that signal is properly distinguished from noise. Currently, there are a limited number of methods for analyzing 4C-Seq data. Here, we present a new method, fourSig, which in addition to being precise and simple to use also includes a new feature that prioritizes detected interactions. Our results demonstrate the efficacy of fourSig with previously published and novel 4C-Seq data sets and show that our significance prioritization correlates with the ability to reproducibly detect interactions among replicates.

Detection of RNA-Protein Interactions in Living Cells with SHAPE.

SHAPE-MaP is unique among RNA structure probing strategies in that it both measures flexibility at single-nucleotide resolution and quantifies the uncertainties in these measurements. We report a straightforward analytical framework that incorporates these uncertainties to allow detection of RNA structural differences between any two states, and we use it here to detect RNA-protein interactions in healthy mouse trophoblast stem cells. We validate this approach by analysis of three model cytoplasmic and nuclear ribonucleoprotein complexes, in 2 min in-cell probing experiments. In contrast, data produced by alternative in-cell SHAPE probing methods correlate poorly (r = 0.2) with those generated by SHAPE-MaP and do not yield accurate signals for RNA-protein interactions. We then examine RNA-protein and RNA-substrate interactions in the RNase MRP complex and, by comparing in-cell interaction sites with disease-associated mutations, characterize these noncoding mutations in terms of molecular phenotype. Together, these results reveal that SHAPE-MaP can define true interaction sites and infer RNA functions under native cellular conditions with limited preexisting knowledge of the proteins or RNAs involved.

A latent pro-survival function for the mir-290-295 cluster in mouse embryonic stem cells.

MicroRNAs (miRNAs) post-transcriptionally regulate the expression of thousands of distinct mRNAs. While some regulatory interactions help to maintain basal cellular functions, others are likely relevant in more specific settings, such as response to stress. Here we describe such a role for the mir-290-295 cluster, the dominant miRNA cluster in mouse embryonic stem cells (mESCs). Examination of a target list generated from bioinformatic prediction, as well as expression data following miRNA loss, revealed strong enrichment for apoptotic regulators, two of which we validated directly: Caspase 2, the most highly conserved mammalian caspase, and Ei24, a p53 transcriptional target. Consistent with these predictions, mESCs lacking miRNAs were more likely to initiate apoptosis following genotoxic exposure to gamma irradiation or doxorubicin. Knockdown of either candidate partially rescued this pro-apoptotic phenotype, as did transfection of members of the mir-290-295 cluster. These findings were recapitulated in a specific mir-290-295 deletion line, confirming that they reflect miRNA functions at physiological levels. In contrast to the basal regulatory roles previously identified, the pro-survival phenotype shown here may be most relevant to stressful gestations, where pro-oxidant metabolic states induce DNA damage. Similarly, this cluster may mediate chemotherapeutic resistance in a neoplastic context, making it a useful clinical target.

Improved functions for non-linear sequence comparison using SEEKR.

SEquence Evaluation through k -mer Representation (SEEKR) is a method of sequence comparison that utilizes sequence substrings called k -mers to quantify non-linear similarity between nucleic acid species. We describe the development of new functions within SEEKR that enable end-users to estimate p-values that ascribe statistical significance to SEEKR-derived similarities as well as visualize different aspects of k -mer similarity. We apply the new functions to identify chromatin-enriched long noncoding RNAs (lncRNAs) that harbor XIST -like sequence fragments and show that several of these fragments are bound by XIST -associated proteins. We also highlight the best practice of using RNA-Seq data to evaluate support for lncRNA annotations prior to their in-depth study in cell types of interest.

A monoclonal antibody raised against human EZH2 cross-reacts with the RNA-binding protein SAFB.

The Polycomb Repressive Complex 2 (PRC2) is a conserved enzyme that tri-methylates Lysine 27 on Histone 3 (H3K27me3) to promote gene silencing. PRC2 is remarkably responsive to the expression of certain long noncoding RNAs (lncRNAs). In the most notable example, PRC2 is recruited to the X-chromosome shortly after expression of the lncRNA Xist begins during X-chromosome inactivation. However, the mechanisms by which lncRNAs recruit PRC2 to chromatin are not yet clear. We report that a broadly used rabbit monoclonal antibody raised against human EZH2, a catalytic subunit of PRC2, cross-reacts with an RNA-binding protein called Scaffold Attachment Factor B (SAFB) in mouse embryonic stem cells (ESCs) under buffer conditions that are commonly used for chromatin immunoprecipitation (ChIP). Knockout of EZH2 in ESCs demonstrated that the antibody is specific for EZH2 by western blot (no cross-reactivity). Likewise, comparison to previously published datasets confirmed that the antibody recovers PRC2-bound sites by ChIP-Seq. However, RNA-IP from formaldehyde-crosslinked ESCs using ChIP wash conditions recovers distinct peaks of RNA association that co-localize with peaks of SAFB and whose enrichment disappears upon knockout of SAFB but not EZH2. IP and mass spectrometry-based proteomics in wild-type and EZH2 knockout ESCs confirm that the EZH2 antibody recovers SAFB in an EZH2-independent manner. Our data highlight the importance of orthogonal assays when studying interactions between chromatin-modifying enzymes and RNA.