Ryan K. Dale

A compendium of RNA-binding motifs for decoding gene regulation

RNA-binding proteins are key regulators of gene expression, yet only a small fraction have been functionally characterized. Here we report a systematic analysis of the RNA motifs recognized by RNA-binding proteins, encompassing 205 distinct genes from 24 diverse eukaryotes. The sequence specificities of RNA-binding proteins display deep evolutionary conservation, and the recognition preferences for a large fraction of metazoan RNA-binding proteins can thus be inferred from their RNA-binding domain sequence. The motifs that we identify in vitro correlate well with in vivo RNA-binding data. Moreover, we can associate them with distinct functional roles in diverse types of post-transcriptional regulation, enabling new insights into the functions of RNA-binding proteins both in normal physiology and in human disease. These data provide an unprecedented overview of RNA-binding proteins and their targets, and constitute an invaluable resource for determining post-transcriptional regulatory mechanisms in eukaryotes.

Cell type specificity of chromatin organization mediated by CTCF and cohesin

CTCF sites are abundant in the genomes of diverse species but their function is enigmatic. We used chromosome conformation capture to determine long-range interactions among CTCF/cohesin sites over 2 Mb on human chromosome 11 encompassing the β-globin locus and flanking olfactory receptor genes. Although CTCF occupies these sites in both erythroid K562 cells and fibroblast 293T cells, the long-range interaction frequencies among the sites are highly cell type specific, revealing a more densely clustered organization in the absence of globin gene activity. Both CTCF and cohesins are required for the cell-type-specific chromatin conformation. Furthermore, loss of the organizational loops in K562 cells through reduction of CTCF with shRNA results in acquisition of repressive histone marks in the globin locus and reduces globin gene expression whereas silent flanking olfactory receptor genes are unaffected. These results support a genome-wide role for CTCF/cohesin sites through loop formation that both influences transcription and contributes to cell-type-specific chromatin organization and function.

Pybedtools: a flexible Python library for manipulating genomic datasets and annotations

Summary: pybedtools is a flexible Python software library for manipulating and exploring genomic datasets in many common formats. It provides an intuitive Python interface that extends upon the popular BEDTools genome arithmetic tools. The library is well documented and efficient, and allows researchers to quickly develop simple, yet powerful scripts that enable complex genomic analyses. Availability: pybedtools is maintained under the GPL license. Stable versions of pybedtools as well as documentation are available on the Python Package Index at http://pypi.python.org/pypi/pybedtools. Contact: dalerr@niddk.nih.gov; arq5x@virginia.edu Supplementary Information: Supplementary data are available at Bioinformatics online.

RNAi-independent role for Argonaute2 in CTCF/CP190 chromatin insulator function

A major role of the RNAi pathway in Schizosaccharomyces pombe is to nucleate heterochromatin, but it remains unclear whether this mechanism is conserved. To address this question in Drosophila, we performed genome-wide localization of Argonaute2 (AGO2) by chromatin immunoprecipitation (ChIP)-seq in two different embryonic cell lines and found that AGO2 localizes to euchromatin but not heterochromatin. This localization pattern is further supported by immunofluorescence staining of polytene chromosomes and cell lines, and these studies also indicate that a substantial fraction of AGO2 resides in the nucleus. Intriguingly, AGO2 colocalizes extensively with CTCF/CP190 chromatin insulators but not with genomic regions corresponding to endogenous siRNA production. Moreover, AGO2, but not its catalytic activity or Dicer-2, is required for CTCF/CP190-dependent Fab-8 insulator function. AGO2 interacts physically with CTCF and CP190, and depletion of either CTCF or CP190 results in genome-wide loss of AGO2 chromatin association. Finally, mutation of CTCF, CP190, or AGO2 leads to reduction of chromosomal looping interactions, thereby altering gene expression. We propose that RNAi-independent recruitment of AGO2 to chromatin by insulator proteins promotes the definition of transcriptional domains throughout the genome.

Mediation of Drosophila autosomal dosage effects and compensation by network interactions

BackgroundGene dosage change is a mild perturbation that is a valuable tool for pathway reconstruction in Drosophila. While it is often assumed that reducing gene dose by half leads to two-fold less expression, there is partial autosomal dosage compensation in Drosophila, which may be mediated by feedback or buffering in expression networks.ResultsWe profiled expression in engineered flies where gene dose was reduced from two to one. While expression of most one-dose genes was reduced, the gene-specific dose responses were heterogeneous. Expression of two-dose genes that are first-degree neighbors of one-dose genes in novel network models also changed, and the directionality of change depended on the response of one-dose genes.ConclusionsOur data indicate that expression perturbation propagates in network space. Autosomal compensation, or the lack thereof, is a gene-specific response, largely mediated by interactions with the rest of the transcriptome.

Ldb1-nucleated transcription complexes function as primary mediators of global erythroid gene activation

Erythropoiesis is dependent on the lineage-specific transcription factors Gata1, Tal1, and Klf1. Several erythroid genes have been shown to require all 3 factors for their expression, suggesting that they function synergistically; however, there is little direct evidence for widespread cooperation. Gata1 and Tal1 can assemble within higher-order protein complexes (Ldb1 complexes) that include the adapter molecules Lmo2 and Ldb1. Ldb1 proteins are capable of coassociation, and long-range Ldb1-mediated oligomerization of enhancer- and promoter-bound Ldb1 complexes has been shown to be required for β-globin gene expression. In this study, we generated a genomewide map of Ldb1 complex binding sites that revealed widespread binding at erythroid genes and at known erythroid enhancer elements. Ldb1 complex binding sites frequently colocalized with Klf1 binding sites and with consensus binding motifs for other erythroid transcription factors. Transcriptomic analysis demonstrated a strong correlation between Ldb1 complex binding and Ldb1 dependency for gene expression and identified a large cohort of genes coregulated by Ldb1 complexes and Klf1. Together, these results provide a foundation for defining the mechanism and scope of Ldb1 complex activity during erythropoiesis.

Role of LDB1 in the transition from chromatin looping to transcription activation

Many questions remain about how close association of genes and distant enhancers occurs and how this is linked to transcription activation. In erythroid cells, lim domain binding 1 (LDB1) protein is recruited to the β-globin locus via LMO2 and is required for looping of the β-globin locus control region (LCR) to the active β-globin promoter. We show that the LDB1 dimerization domain (DD) is necessary and, when fused to LMO2, sufficient to completely restore LCR-promoter looping and transcription in LDB1-depleted cells. The looping function of the DD is unique and irreplaceable by heterologous DDs. Dissection of the DD revealed distinct functional properties of conserved subdomains. Notably, a conserved helical region (DD4/5) is dispensable for LDB1 dimerization and chromatin looping but essential for transcriptional activation. DD4/5 is required for the recruitment of the coregulators FOG1 and the nucleosome remodeling and deacetylating (NuRD) complex. Lack of DD4/5 alters histone acetylation and RNA polymerase II recruitment and results in failure of the locus to migrate to the nuclear interior, as normally occurs during erythroid maturation. These results uncouple enhancer-promoter looping from nuclear migration and transcription activation and reveal new roles for LDB1 in these processes.

Distinct Ldb1/NLI complexes orchestrate γ-globin repression and reactivation through ETO2 in human adult erythroid cells.

The Ldb1/GATA-1/TAL1/LMO2 complex mediates long-range interaction between the β-globin locus control region (LCR) and gene in adult mouse erythroid cells, but whether this complex mediates chromatin interactions at other developmental stages or in human cells is unknown. We investigated NLI (Ldb1 homolog) complex occupancy and chromatin conformation of the β-globin locus in human erythroid cells. In addition to the LCR, we found robust NLI complex occupancy at a site downstream of the (A)γ-globin gene within sequences of BGL3, an intergenic RNA transcript. In cells primarily transcribing β-globin, BGL3 is not transcribed and BGL3 sequences are occupied by NLI core complex members, together with corepressor ETO2 and by γ-globin repressor BCL11A. The LCR and β-globin gene establish proximity in these cells. In contrast, when γ-globin transcription is reactivated in these cells, ETO2 participation in the NLI complex at BGL3 is diminished, as is BCL11A occupancy, and both BGL3 and γ-globin are transcribed. In these cells, proximity between the BGL3/γ-globin region and the LCR is established. We conclude that alternative NLI complexes mediate γ-globin transcription or silencing through long-range LCR interactions involving an intergenic site of noncoding RNA transcription and that ETO2 is critical to this process.

Tissue-specific regulation of chromatin insulator function.

Chromatin insulators organize the genome into distinct transcriptional domains and contribute to cell type–specific chromatin organization. However, factors regulating tissue-specific insulator function have not yet been discovered. Here we identify the RNA recognition motif-containing protein Shep as a direct interactor of two individual components of the gypsy insulator complex in Drosophila. Mutation of shep improves gypsy-dependent enhancer blocking, indicating a role as a negative regulator of insulator activity. Unlike ubiquitously expressed core gypsy insulator proteins, Shep is highly expressed in the central nervous system (CNS) with lower expression in other tissues. We developed a novel, quantitative tissue-specific barrier assay to demonstrate that Shep functions as a negative regulator of insulator activity in the CNS but not in muscle tissue. Additionally, mutation of shep alters insulator complex nuclear localization in the CNS but has no effect in other tissues. Consistent with negative regulatory activity, ChIP–seq analysis of Shep in a CNS-derived cell line indicates substantial genome-wide colocalization with a single gypsy insulator component but limited overlap with intact insulator complexes. Taken together, these data reveal a novel, tissue-specific mode of regulation of a chromatin insulator.

CTCF-mediated transcriptional regulation through cell type-specific chromosome organization in the β-globin locus

The principles underlying the architectural landscape of chromatin beyond the nucleosome level in living cells remains largely unknown despite its potential to play a role in mammalian gene regulation. We investigated the three-dimensional folding of a 1 Mbp region of human chromosome 11 containing the β-globin genes by integrating looping interactions of the CCCTC-binding insulator protein CTCF determined comprehensively by chromosome conformation capture (3C) into a polymer model of chromatin. We find that CTCF-mediated cell type-specific interactions in erythroid cells are organized to favor contacts known to occur in vivo between the β-globin locus control region (LCR) and genes. In these cells, the modeled β-globin domain folds into a globule with the LCR and the active globin genes on the periphery. In contrast, in non-erythroid cells, the globule is less compact with few but dominant CTCF interactions driving the genes away from the LCR. This leads to a decrease in contact frequencies that can exceed 1000-fold depending on the stiffness of the chromatin and the exact position of the genes. Our findings show that an ensemble of CTCF contacts functionally affects spatial distances between control elements and target genes contributing to chromosomal organization required for transcription.