Higher-Order Inter-chromosomal Hubs Shape 3D Genome Organization in the Nucleus.

Eukaryotic genomes are packaged into a 3-dimensional structure in the nucleus. Current methods for studying genome-wide structure are based on proximity ligation. However, this approach can fail to detect known structures, such as interactions with nuclear bodies, because these DNA regions can be too far apart to directly ligate. Accordingly, our overall understanding of genome organization remains incomplete. Here, we develop split-pool recognition of interactions by tag extension (SPRITE), a method that enables genome-wide detection of higher-order interactions within the nucleus. Using SPRITE, we recapitulate known structures identified by proximity ligation and identify additional interactions occurring across larger distances, including two hubs of inter-chromosomal interactions that are arranged around the nucleolus and nuclear speckles. We show that a substantial fraction of the genome exhibits preferential organization relative to these nuclear bodies. Our results generate a global model whereby nuclear bodies act as inter-chromosomal hubs that shape the overall packaging of DNA in the nucleus.

RNA-RNA interactions enable specific targeting of noncoding RNAs to nascent Pre-mRNAs and chromatin sites.

Intermolecular RNA-RNA interactions are used by many noncoding RNAs (ncRNAs) to achieve their diverse functions. To identify these contacts, we developed a method based on RNA antisense purification to systematically map RNA-RNA interactions (RAP-RNA) and applied it to investigate two ncRNAs implicated in RNA processing: U1 small nuclear RNA, a component of the spliceosome, and Malat1, a large ncRNA that localizes to nuclear speckles. U1 and Malat1 interact with nascent transcripts through distinct targeting mechanisms. Using differential crosslinking, we confirmed that U1 directly hybridizes to 5' splice sites and 5' splice site motifs throughout introns and found that Malat1 interacts with pre-mRNAs indirectly through protein intermediates. Interactions with nascent pre-mRNAs cause U1 and Malat1 to localize proximally to chromatin at active genes, demonstrating that ncRNAs can use RNA-RNA interactions to target specific pre-mRNAs and genomic sites. RAP-RNA is sensitive to lower abundance RNAs as well, making it generally applicable for investigating ncRNAs.

Transcriptional and epigenetic dynamics during specification of human embryonic stem cells.

Differentiation of human embryonic stem cells (hESCs) provides a unique opportunity to study the regulatory mechanisms that facilitate cellular transitions in a human context. To that end, we performed comprehensive transcriptional and epigenetic profiling of populations derived through directed differentiation of hESCs representing each of the three embryonic germ layers. Integration of whole-genome bisulfite sequencing, chromatin immunoprecipitation sequencing, and RNA sequencing reveals unique events associated with specification toward each lineage. Lineage-specific dynamic alterations in DNA methylation and H3K4me1 are evident at putative distal regulatory elements that are frequently bound by pluripotency factors in the undifferentiated hESCs. In addition, we identified germ-layer-specific H3K27me3 enrichment at sites exhibiting high DNA methylation in the undifferentiated state. A better understanding of these initial specification events will facilitate identification of deficiencies in current approaches, leading to more faithful differentiation strategies as well as providing insights into the rewiring of human regulatory programs during cellular transitions.

Pooled CRISPR screens with imaging on microraft arrays reveals stress granule-regulatory factors.

Genetic screens using pooled CRISPR-based approaches are scalable and inexpensive, but restricted to standard readouts, including survival, proliferation and sortable markers. However, many biologically relevant cell states involve cellular and subcellular changes that are only accessible by microscopic visualization, and are currently impossible to screen with pooled methods. Here we combine pooled CRISPR-Cas9 screening with microraft array technology and high-content imaging to screen image-based phenotypes (CRaft-ID; CRISPR-based microRaft followed by guide RNA identification). By isolating microrafts that contain genetic clones harboring individual guide RNAs (gRNA), we identify RNA-binding proteins (RBPs) that influence the formation of stress granules, the punctate protein-RNA assemblies that form during stress. To automate hit identification, we developed a machine-learning model trained on nuclear morphology to remove unhealthy cells or imaging artifacts. In doing so, we identified and validated previously uncharacterized RBPs that modulate stress granule abundance, highlighting the applicability of our approach to facilitate image-based pooled CRISPR screens.

The Xist lncRNA interacts directly with SHARP to silence transcription through HDAC3.

Many long non-coding RNAs (lncRNAs) affect gene expression, but the mechanisms by which they act are still largely unknown. One of the best-studied lncRNAs is Xist, which is required for transcriptional silencing of one X chromosome during development in female mammals. Despite extensive efforts to define the mechanism of Xist-mediated transcriptional silencing, we still do not know any proteins required for this role. The main challenge is that there are currently no methods to comprehensively define the proteins that directly interact with a lncRNA in the cell. Here we develop a method to purify a lncRNA from cells and identify proteins interacting with it directly using quantitative mass spectrometry. We identify ten proteins that specifically associate with Xist, three of these proteins--SHARP, SAF-A and LBR--are required for Xist-mediated transcriptional silencing. We show that SHARP, which interacts with the SMRT co-repressor that activates HDAC3, is not only essential for silencing, but is also required for the exclusion of RNA polymerase II (Pol II) from the inactive X. Both SMRT and HDAC3 are also required for silencing and Pol II exclusion. In addition to silencing transcription, SHARP and HDAC3 are required for Xist-mediated recruitment of the polycomb repressive complex 2 (PRC2) across the X chromosome. Our results suggest that Xist silences transcription by directly interacting with SHARP, recruiting SMRT, activating HDAC3, and deacetylating histones to exclude Pol II across the X chromosome.

Genetic and epigenetic fine mapping of causal autoimmune disease variants.

Genome-wide association studies have identified loci underlying human diseases, but the causal nucleotide changes and mechanisms remain largely unknown. Here we developed a fine-mapping algorithm to identify candidate causal variants for 21 autoimmune diseases from genotyping data. We integrated these predictions with transcription and cis-regulatory element annotations, derived by mapping RNA and chromatin in primary immune cells, including resting and stimulated CD4(+) T-cell subsets, regulatory T cells, CD8(+) T cells, B cells, and monocytes. We find that ∼90% of causal variants are non-coding, with ∼60% mapping to immune-cell enhancers, many of which gain histone acetylation and transcribe enhancer-associated RNA upon immune stimulation. Causal variants tend to occur near binding sites for master regulators of immune differentiation and stimulus-dependent gene activation, but only 10-20% directly alter recognizable transcription factor binding motifs. Rather, most non-coding risk variants, including those that alter gene expression, affect non-canonical sequence determinants not well-explained by current gene regulatory models.

Genome-wide screen identifies pathways that govern GAA/TTC repeat fragility and expansions in dividing and nondividing yeast cells.

Triplex structure-forming GAA/TTC repeats pose a dual threat to the eukaryotic genome integrity. Their potential to expand can lead to gene inactivation, the cause of Friedreich's ataxia disease in humans. In model systems, long GAA/TTC tracts also act as chromosomal fragile sites that can trigger gross chromosomal rearrangements. The mechanisms that regulate the metabolism of GAA/TTC repeats are poorly understood. We have developed an experimental system in the yeast Saccharomyces cerevisiae that allows us to systematically identify genes crucial for maintaining the repeat stability. Two major groups of mutants defective in DNA replication or transcription initiation are found to be prone to fragility and large-scale expansions. We demonstrate that problems imposed by the repeats during DNA replication in actively dividing cells and during transcription initiation in nondividing cells can culminate in genome instability. We propose that similar mechanisms can mediate detrimental metabolism of GAA/TTC tracts in human cells.

Robust transcriptome-wide discovery of RNA-binding protein binding sites with enhanced CLIP (eCLIP).

As RNA-binding proteins (RBPs) play essential roles in cellular physiology by interacting with target RNA molecules, binding site identification by UV crosslinking and immunoprecipitation (CLIP) of ribonucleoprotein complexes is critical to understanding RBP function. However, current CLIP protocols are technically demanding and yield low-complexity libraries with high experimental failure rates. We have developed an enhanced CLIP (eCLIP) protocol that decreases requisite amplification by ∼1,000-fold, decreasing discarded PCR duplicate reads by ∼60% while maintaining single-nucleotide binding resolution. By simplifying the generation of paired IgG and size-matched input controls, eCLIP improves specificity in the discovery of authentic binding sites. We generated 102 eCLIP experiments for 73 diverse RBPs in HepG2 and K562 cells (available at https://www.encodeproject.org), demonstrating that eCLIP enables large-scale and robust profiling, with amplification and sample requirements similar to those of ChIP-seq. eCLIP enables integrative analysis of diverse RBPs to reveal factor-specific profiles, common artifacts for CLIP and RNA-centric perspectives on RBP activity.

Defining the 5΄ and 3΄ landscape of the Drosophila transcriptome with Exo-seq and RNaseH-seq.

Cells regulate biological responses in part through changes in transcription start sites (TSS) or cleavage and polyadenylation sites (PAS). To fully understand gene regulatory networks, it is therefore critical to accurately annotate cell type-specific TSS and PAS. Here we present a simple and straightforward approach for genome-wide annotation of 5΄- and 3΄-RNA ends. Our approach reliably discerns bona fide PAS from false PAS that arise due to internal poly(A) tracts, a common problem with current PAS annotation methods. We applied our methodology to study the impact of temperature on the Drosophila melanogaster head transcriptome. We found hundreds of previously unidentified TSS and PAS which revealed two interesting phenomena: first, genes with multiple PASs tend to harbor a motif near the most proximal PAS, which likely represents a new cleavage and polyadenylation signal. Second, motif analysis of promoters of genes affected by temperature suggested that boundary element association factor of 32 kDa (BEAF-32) and DREF mediates a transcriptional program at warm temperatures, a result we validated in a fly line where beaf-32 is downregulated. These results demonstrate the utility of a high-throughput platform for complete experimental and computational analysis of mRNA-ends to improve gene annotation.

Simultaneous generation of many RNA-seq libraries in a single reaction.

Although RNA-seq is a powerful tool, the considerable time and cost associated with library construction has limited its utilization for various applications. RNAtag-Seq, an approach to generate multiple RNA-seq libraries in a single reaction, lowers time and cost per sample, and it produces data on prokaryotic and eukaryotic samples that are comparable to those generated by traditional strand-specific RNA-seq approaches.

Variation in single-nucleotide sensitivity of eCLIP derived from reverse transcription conditions.

Crosslinking and immunoprecipitation (CLIP) followed by high-throughput sequencing identifies the binding sites of RNA binding proteins on RNAs. The covalent RNA-amino acid adducts produced by UV irradiation can cause premature reverse transcription termination and deletions (referred to as crosslink-induced mutation sites (CIMS)), which may decrease overall cDNA yield but are exploited in state-of-the-art CLIP methods to identify these crosslink sites at single-nucleotide resolution. Here, we show the ratio of both crosslinked base deletions and read-through versus termination are highly dependent on the identity of the reverse transcriptase enzyme as well as on buffer conditions used. AffinityScript and TGIRT showed a lack of deletion of the crosslinked base with other enzymes showing variable rates, indicating that utilization and interpretation of CIMS analysis requires knowledge of the reverse transcriptase enzyme used. Commonly used enzymes, including Superscript III and AffinityScript, show high termination rates in standard magnesium buffer conditions, but show a single base difference in the position of termination for TARDBP motifs. In contrast, manganese-containing buffer promoted read-through at the adduct site. These results validate the use of standard enzymes and also propose alternative enzyme and buffer choices for particularly challenging samples that contain extensive RNA adducts or other modifications that inhibit standard reverse transcription.

Large-scale expansions of Friedreich's ataxia GAA repeats in yeast.

Large-scale expansions of DNA repeats are implicated in numerous hereditary disorders in humans. We describe a yeast experimental system to analyze large-scale expansions of triplet GAA repeats responsible for the human disease Friedreich's ataxia. When GAA repeats were placed into an intron of the chimeric URA3 gene, their expansions caused gene inactivation, which was detected on the selective media. We found that the rates of expansions of GAA repeats increased exponentially with their lengths. These rates were only mildly dependent on the repeat's orientation within the replicon, whereas the repeat-mediated replication fork stalling was exquisitely orientation dependent. Expansion rates were significantly elevated upon inactivation of the replication fork stabilizers, Tof1 and Csm3, but decreased in the knockouts of postreplication DNA repair proteins, Rad6 and Rad5, and the DNA helicase Sgs1. We propose a model for large-scale repeat expansions based on template switching during replication fork progression through repetitive DNA.

Genome rearrangements caused by interstitial telomeric sequences in yeast.

Interstitial telomeric sequences (ITSs) are present in many eukaryotic genomes and are linked to genome instabilities and disease in humans. The mechanisms responsible for ITS-mediated genome instability are not understood in molecular detail. Here, we use a model Saccharomyces cerevisiae system to characterize genome instability mediated by yeast telomeric (Ytel) repeats embedded within an intron of a reporter gene inside a yeast chromosome. We observed a very high rate of small insertions and deletions within the repeats. We also found frequent gross chromosome rearrangements, including deletions, duplications, inversions, translocations, and formation of acentric minichromosomes. The inversions are a unique class of chromosome rearrangement involving an interaction between the ITS and the true telomere of the chromosome. Because we previously found that Ytel repeats cause strong replication fork stalling, we suggest that formation of double-stranded DNA breaks within the Ytel sequences might be responsible for these gross chromosome rearrangements.

Expansions, contractions, and fragility of the spinocerebellar ataxia type 10 pentanucleotide repeat in yeast.

Spinocerebellar ataxia 10 (SCA10) is an autosomal dominant disease caused by large-scale expansions of the (ATTCT)(n) repeat within an intron of the human ATXN10 gene. In contrast to other expandable repeats, this pentanucleotide repeat does not form stable intra- or interstranded DNA structures, being a DNA unwinding element instead. We analyzed the instability of the (ATTCT)(n) repeat in a yeast experimental system, where its expansions led to inactivation of the URA3 reporter gene. The inactivation was due to a dramatic decrease in the mRNA levels owing to premature transcription termination and RNA polyadenylation at the repeat. The rates of expansions strongly increased with the repeat's length, mimicking genetic anticipation in human pedigrees. A first round of genetic analysis showed that a functional TOF1 gene precludes, whereas a functional RAD5 gene promotes, expansions of the (ATTCT)(n) repeat. We hypothesize that repeat expansions could occur upon fortuitous template switching during DNA replication. The rate of repeat contractions was elevated in the Tof1 knockout strain, but it was not affected by the RAD5 gene. Supporting the notion of replication irregularities, we found that (ATTCT)(n) repeats also cause length-dependent chromosomal fragility in yeast. Repeat-mediated fragility was also affected by the Tof1 and Rad5 proteins, being reduced in their absence.

Partners in crime: Tbf1 and Vid22 promote expansions of long human telomeric repeats at an interstitial chromosome position in yeast.

In humans, telomeric repeats (TTAGGG)n are known to be present at internal chromosomal sites. These interstitial telomeric sequences (ITSs) are an important source of genomic instability, including repeat length polymorphism, but the molecular mechanisms responsible for this instability remain to be understood. Here, we studied the mechanisms responsible for expansions of human telomeric (Htel) repeats that were artificially inserted inside a yeast chromosome. We found that Htel repeats in an interstitial chromosome position are prone to expansions. The propensity of Htel repeats to expand depends on the presence of a complex of two yeast proteins: Tbf1 and Vid22. These two proteins are physically bound to an interstitial Htel repeat, and together they slow replication fork progression through it. We propose that slow progression of the replication fork through the protein complex formed by the Tbf1 and Vid22 partners at the Htel repeat cause DNA strand slippage, ultimately resulting in repeat expansions.

Transcriptome-wide identification of RNA-binding protein binding sites using seCLIP-seq.

Discovery of interaction sites between RNA-binding proteins (RBPs) and their RNA targets plays a critical role in enabling our understanding of how these RBPs control RNA processing and regulation. Cross-linking and immunoprecipitation (CLIP) provides a generalizable, transcriptome-wide method by which RBP/RNA complexes are purified and sequenced to identify sites of intermolecular contact. By simplifying technical challenges in prior CLIP methods and incorporating the generation of and quantitative comparison against size-matched input controls, the single-end enhanced CLIP (seCLIP) protocol allows for the profiling of these interactions with high resolution, efficiency and scalability. Here, we present a step-by-step guide to the seCLIP method, detailing critical steps and offering insights regarding troubleshooting and expected results while carrying out the ~4-d protocol. Furthermore, we describe a comprehensive bioinformatics pipeline that offers users the tools necessary to process two replicate datasets and identify reproducible and significant peaks for an RBP of interest in ~2 d.

Identification of the global miR-130a targetome reveals a role for TBL1XR1 in hematopoietic stem cell self-renewal and t(8;21) AML.

Gene expression profiling and proteome analysis of normal and malignant hematopoietic stem cells (HSCs) point to shared core stemness properties. However, discordance between mRNA and protein signatures highlights an important role for post-transcriptional regulation by microRNAs (miRNAs) in governing this critical nexus. Here, we identify miR-130a as a regulator of HSC self-renewal and differentiation. Enforced expression of miR-130a impairs B lymphoid differentiation and expands long-term HSCs. Integration of protein mass spectrometry and chimeric AGO2 crosslinking and immunoprecipitation (CLIP) identifies TBL1XR1 as a primary miR-130a target, whose loss of function phenocopies miR-130a overexpression. Moreover, we report that miR-130a is highly expressed in t(8;21) acute myeloid leukemia (AML), where it is critical for maintaining the oncogenic molecular program mediated by the AML1-ETO complex. Our study establishes that identification of the comprehensive miRNA targetome within primary cells enables discovery of genes and molecular networks underpinning stemness properties of normal and leukemic cells.

Microfluidic SlipChip device for multistep multiplexed biochemistry on a nanoliter scale.

We have developed a multistep microfluidic device that expands the current SlipChip capabilities by enabling multiple steps of droplet merging and multiplexing. Harnessing the interfacial energy between carrier and sample phases, this manually operated device accurately meters nanoliter volumes of reagents and transfers them into on-device reaction wells. Judiciously shaped microfeatures and surface-energy traps merge droplets in a parallel fashion. Wells can be tuned for different volumetric capacities and reagent types, including for pre-spotted reagents that allow for unique identification of original well contents even after their contents are pooled. We demonstrate the functionality of the multistep SlipChip by performing RNA transcript barcoding on-device for synthetic spiked-in standards and for biologically derived samples. This technology is a good candidate for a wide range of biological applications that require multiplexing of multistep reactions in nanoliter volumes, including single-cell analyses.

A Defective mRNA Cleavage and Polyadenylation Complex Facilitates Expansions of Transcribed (GAA)n Repeats Associated with Friedreich's Ataxia.

Expansions of microsatellite repeats are responsible for numerous hereditary diseases in humans, including myotonic dystrophy and Friedreich's ataxia. Whereas the length of an expandable repeat is the main factor determining disease inheritance, recent data point to genomic trans modifiers that can impact the likelihood of expansions and disease progression. Detection of these modifiers may lead to understanding and treating repeat expansion diseases. Here, we describe a method for the rapid, genome-wide identification of trans modifiers for repeat expansion in a yeast experimental system. Using this method, we found that missense mutations in the endoribonuclease subunit (Ysh1) of the mRNA cleavage and polyadenylation complex dramatically increase the rate of (GAA)n repeat expansions but only when they are actively transcribed. These expansions correlate with slower transcription elongation caused by the ysh1 mutation. These results reveal an interplay between RNA processing and repeat-mediated genome instability, confirming the validity of our approach.

Evolutionary analysis across mammals reveals distinct classes of long non-coding RNAs.

BACKGROUND: Recent advances in transcriptome sequencing have enabled the discovery of thousands of long non-coding RNAs (lncRNAs) across many species. Though several lncRNAs have been shown to play important roles in diverse biological processes, the functions and mechanisms of most lncRNAs remain unknown. Two significant obstacles lie between transcriptome sequencing and functional characterization of lncRNAs: identifying truly non-coding genes from de novo reconstructed transcriptomes, and prioritizing the hundreds of resulting putative lncRNAs for downstream experimental interrogation. RESULTS: We present slncky, a lncRNA discovery tool that produces a high-quality set of lncRNAs from RNA-sequencing data and further uses evolutionary constraint to prioritize lncRNAs that are likely to be functionally important. Our automated filtering pipeline is comparable to manual curation efforts and more sensitive than previously published computational approaches. Furthermore, we developed a sensitive alignment pipeline for aligning lncRNA loci and propose new evolutionary metrics relevant for analyzing sequence and transcript evolution. Our analysis reveals that evolutionary selection acts in several distinct patterns, and uncovers two notable classes of intergenic lncRNAs: one showing strong purifying selection on RNA sequence and another where constraint is restricted to the regulation but not the sequence of the transcript. CONCLUSION: Our results highlight that lncRNAs are not a homogenous class of molecules but rather a mixture of multiple functional classes with distinct biological mechanism and/or roles. Our novel comparative methods for lncRNAs reveals 233 constrained lncRNAs out of tens of thousands of currently annotated transcripts, which we make available through the slncky Evolution Browser.