Search

The design of high-affinity, RNA-binding ligands has proven very challenging. This is due to the unique structural properties of RNA, often characterized by polar surfaces and high flexibility. In addition, the frequent lack of well-defined binding pockets complicates the development of small molecule binders. This has triggered the search for alternative scaffolds of intermediate size. Among these, peptide-derived molecules represent appealing entities as they can mimic structural features also present in RNA-binding proteins. However, the application of peptidic RNA-targeting ligands is hampered by a lack of design principles and their inherently low bio-stability. Here, the structure-based design of constrained α-helical peptides derived from the viral suppressor of RNA silencing, TAV2b, is described. We observe that the introduction of two inter-side chain crosslinks provides peptides with increased α-helicity and protease stability. One of these modified peptides (B3) shows high affinity for double-stranded RNA structures including a palindromic siRNA as well as microRNA-21 and its precursor pre-miR-21. Notably, B3 binding to pre-miR-21 inhibits Dicer processing in a biochemical assay. As a further characteristic this peptide also exhibits cellular entry. Our findings show that constrained peptides can efficiently mimic RNA-binding proteins rendering them potentially useful for the design of bioactive RNA-targeting ligands.

Human fetal globin (γ-globin) genes are developmentally silenced after birth, and reactivation of γ-globin expression in adulthood ameliorates symptoms of hemoglobin disorders, such as sickle cell disease (SCD) and β-thalassemia. However, the mechanisms by which γ-globin expression is precisely regulated are still incompletely understood. Here, we found that NonO (non-POU domain-containing octamer-binding protein) interacted directly with SOX6, and repressed the expression of γ-globin gene in human erythroid cells. We showed that NonO bound to the octamer binding motif, ATGCAAAT, of the γ-globin proximal promoter, resulting in inhibition of γ-globin transcription. Depletion of NonO resulted in significant activation of γ-globin expression in K562, HUDEP-2, and primary human erythroid progenitor cells. To confirm the role of NonO in vivo, we further generated a conditional knockout of NonO by using IFN-inducible Mx1-Cre transgenic mice. We found that induced NonO deletion reactivated murine embryonic globin and human γ-globin gene expression in adult β-YAC mice, suggesting a conserved role for NonO during mammalian evolution. Thus, our data indicate that NonO acts as a novel transcriptional repressor of γ-globin gene expression through direct promoter binding, and is essential for γ-globin gene silencing.

Chromosome stability is primarily determined by telomere length. TRF1 is the core subunit of shelterin that plays a critical role in telomere organization and replication. However, the dynamics of TRF1 in scenarios of telomere-processing activities remain elusive. Using single-molecule magnetic tweezers, we here investigated the dynamics of TRF1 upon organizing a human telomere and the protein-DNA interactions at a moving telomeric fork. We first developed a method to obtain telomeres from human cells for directly measuring the telomere length by single-molecule force spectroscopy. Next, we examined the compaction and decompaction of a telomere by TRF1 dimers. TRF1 dissociates from a compacted telomere with heterogenous loops in ∼20 s. We also found a negative correlation between the number of telomeric loops and loop sizes. We further characterized the dynamics of TRF1 at a telomeric DNA fork. With binding energies of 11 kBT, TRF1 can modulate the forward and backward steps of DNA fork movements by 2–9 s at a critical force of F1/2, temporarily maintaining the telomeric fork open. Our results shed light on the mechanisms of how TRF1 organizes human telomeres and facilitates the efficient replication of telomeric DNA. Our work will help future research on the chemical biology of telomeres and shelterin-targeted drug discovery., Graphical Abstract Graphical AbstractWe use single-molecule magnetic tweezers to precisely measure the absolute length of single human telomeres and reveal the dynamics of TRF1 upon organizing a telomere.

The microprocessor complex cleaves the primary transcript of microRNA (pri-miRNA) to initiate miRNA maturation. Microprocessor is known to consist of RNase III DROSHA and dsRNA-binding DGCR8. Here, we identify Enhancer of Rudimentary Homolog (ERH) as a new component of Microprocessor. Through a crystal structure and biochemical experiments, we reveal that ERH uses its hydrophobic groove to bind to a conserved region in the N-terminus of DGCR8, in a 2:2 stoichiometry. Knock-down of ERH or deletion of the DGCR8 N-terminus results in a reduced processing of suboptimal pri-miRNAs in polycistronic miRNA clusters. ERH increases the processing of suboptimal pri-miR-451 in a manner dependent on its neighboring pri-miR-144. Thus, the ERH dimer may mediate ‘cluster assistance’ in which Microprocessor is loaded onto a poor substrate with help from a high-affinity substrate in the same cluster. Our study reveals a role of ERH in the miRNA biogenesis pathway.

The underlying mechanism of transcriptional co-repressor ETO2 during early erythropoiesis and hemoglobin switching is unclear. We find that absence of ETO2 in mice interferes with down-regulation of PU.1 and GATA2 in the fetal liver, impeding a key step required for commitment to erythroid maturation. In human β-globin transgenic Eto2 null mice and in human CD34+ erythroid progenitor cells with reduced ETO2, loss of ETO2 results in ineffective silencing of embryonic/fetal globin gene expression, impeding hemoglobin switching during erythroid differentiation. ETO2 occupancy genome-wide occurs virtually exclusively at LDB1-complex binding sites in enhancers and ETO2 loss leads to increased enhancer activity and expression of target genes. ETO2 recruits the NuRD nucleosome remodeling and deacetylation complex to regulate histone acetylation and nucleosome occupancy in the β-globin locus control region and γ-globin gene. Loss of ETO2 elevates LDB1, MED1 and Pol II in the locus and facilitates fetal γ-globin/LCR looping and γ-globin transcription. Absence of the ETO2 hydrophobic heptad repeat region impairs ETO2-NuRD interaction and function in antagonizing γ-globin/LCR looping. Our results reveal a pivotal role for ETO2 in erythropoiesis and globin gene switching through its repressive role in the LDB1 complex, affecting the transcription factor and epigenetic environment and ultimately restructuring chromatin organization.

Extracellular RNAs participate in intercellular communication, and are being studied as promising minimally invasive diagnostic markers. Several studies in recent years showed that tRNA halves and distinct Y RNA fragments are abundant in the extracellular space, including in biofluids. While their regulatory and diagnostic potential has gained a substantial amount of attention, the biogenesis of these extracellular RNA fragments remains largely unexplored. Here, we demonstrate that these fragments are produced by RNase 1, a highly active secreted nuclease. We use RNA sequencing to investigate the effect of a null mutation of RNase 1 on the levels of tRNA halves and Y RNA fragments in the extracellular environment of cultured human cells. We complement and extend our RNA sequencing results with northern blots, showing that tRNAs and Y RNAs in the non-vesicular extracellular compartment are released from cells as full-length precursors and are subsequently cleaved to distinct fragments. In support of these results, formation of tRNA halves is recapitulated by recombinant human RNase 1 in our in vitro assay. These findings assign a novel function for RNase 1, and position it as a strong candidate for generation of tRNA halves and Y RNA fragments in biofluids.

Precise genomic modification using prime editing (PE) holds enormous potential for research and clinical applications. In this study, we generated all-in-one prime editing (PEA1) constructs that carry all the components required for PE, along with a selection marker. We tested these constructs (with selection) in HEK293T, K562, HeLa and mouse embryonic stem (ES) cells. We discovered that PE efficiency in HEK293T cells was much higher than previously observed, reaching up to 95% (mean 67%). The efficiency in K562 and HeLa cells, however, remained low. To improve PE efficiency in K562 and HeLa, we generated a nuclease prime editor and tested this system in these cell lines as well as mouse ES cells. PE-nuclease greatly increased prime editing initiation, however, installation of the intended edits was often accompanied by extra insertions derived from the repair template. Finally, we show that zygotic injection of the nuclease prime editor can generate correct modifications in mouse fetuses with up to 100% efficiency.

The identification of transcription factor binding sites and cis-regulatory motifs is a frontier whereupon the rules governing protein–DNA binding are being revealed. Here, we developed a new method (DEep Sequence and Shape mOtif or DESSO) for cis-regulatory motif prediction using deep neural networks and the binomial distribution model. DESSO outperformed existing tools, including DeepBind, in predicting motifs in 690 human ENCODE ChIP-sequencing datasets. Furthermore, the deep-learning framework of DESSO expanded motif discovery beyond the state-of-the-art by allowing the identification of known and new protein–protein–DNA tethering interactions in human transcription factors (TFs). Specifically, 61 putative tethering interactions were identified among the 100 TFs expressed in the K562 cell line. In this work, the power of DESSO was further expanded by integrating the detection of DNA shape features. We found that shape information has strong predictive power for TF–DNA binding and provides new putative shape motif information for human TFs. Thus, DESSO improves in the identification and structural analysis of TF binding sites, by integrating the complexities of DNA binding into a deep-learning framework.

We have investigated the structure of the most compact 30-nm chromatin fibres by modelling those with 2-start or 1-start crossed-linker organisations. Using an iterative procedure we obtained possible structural solutions for fibres of the highest possible compaction permitted by physical constraints, including the helical repeat of linker DNA. We find that this procedure predicts a quantized nucleosome repeat length (NRL) and that only fibres with longer NRLs (≥197 bp) can more likely adopt the 1-start organisation. The transition from 2-start to 1-start fibres is consistent with reported differing binding modes of the linker histone. We also calculate that in 1-start fibres the DNA constrains more torsion (as writhe) than 2-start fibres with the same NRL and that the maximum constraint obtained is in accord with previous experimental results. We posit that the coiling of the fibre is driven by overtwisting of linker DNA which, in the most compact forms - for example, in echinoderm sperm and avian erythrocytes - could adopt a helical repeat of ∼10 bp/turn. We argue that in vivo the total twist of linker DNA could be modulated by interaction with other abundant chromatin-associated proteins and by epigenetic modifications of the C-terminal tail of linker histones.

High-throughput biological technologies (e.g. ChIP-seq, RNA-seq and single-cell RNA-seq) rapidly accelerate the accumulation of genome-wide omics data in diverse interrelated biological scenarios (e.g. cells, tissues and conditions). Integration and differential analysis are two common paradigms for exploring and analyzing such data. However, current integrative methods usually ignore the differential part, and typical differential analysis methods either fail to identify combinatorial patterns of difference or require matched dimensions of the data. Here, we propose a flexible framework CSMF to combine them into one paradigm to simultaneously reveal Common and Specific patterns via Matrix Factorization from data generated under interrelated biological scenarios. We demonstrate the effectiveness of CSMF with four representative applications including pairwise ChIP-seq data describing the chromatin modification map between K562 and Huvec cell lines; pairwise RNA-seq data representing the expression profiles of two different cancers; RNA-seq data of three breast cancer subtypes; and single-cell RNA-seq data of human embryonic stem cell differentiation at six time points. Extensive analysis yields novel insights into hidden combinatorial patterns in these multi-modal data. Results demonstrate that CSMF is a powerful tool to uncover common and specific patterns with significant biological implications from data of interrelated biological scenarios.

Sickle cell disease (SCD) is a monogenic disorder that affects millions worldwide. Allogeneic hematopoietic stem cell transplantation is the only available cure. Here, we demonstrate the use of CRISPR/Cas9 and a short single-stranded oligonucleotide template to correct the sickle mutation in the β-globin gene in hematopoietic stem and progenitor cells (HSPCs) from peripheral blood or bone marrow of patients with SCD, with 24.5 ± 7.6% efficiency without selection. Erythrocytes derived from gene-edited cells showed a marked reduction of sickle cells, with the level of normal hemoglobin (HbA) increased to 25.3 ± 13.9%. Gene-corrected SCD HSPCs retained the ability to engraft when transplanted into non-obese diabetic (NOD)-SCID-gamma (NSG) mice with detectable levels of gene correction 16–19 weeks post-transplantation. We show that, by using a high-fidelity SpyCas9 that maintained the same level of on-target gene modification, the off-target effects including chromosomal rearrangements were significantly reduced. Taken together, our results demonstrate efficient gene correction of the sickle mutation in both peripheral blood and bone marrow-derived SCD HSPCs, a significant reduction in sickling of red blood cells, engraftment of gene-edited SCD HSPCs in vivo and the importance of reducing off-target effects; all are essential for moving genome editing based SCD treatment into clinical practice.

Type V CRISPR–Cas12a systems provide an alternate nuclease platform to Cas9, with potential advantages for specific genome editing applications. Here we describe improvements to the Cas12a system that facilitate efficient targeted mutagenesis in mammalian cells and zebrafish embryos. We show that engineered variants of Cas12a with two different nuclear localization sequences (NLS) on the C terminus provide increased editing efficiency in mammalian cells. Additionally, we find that pre-crRNAs comprising a full-length direct repeat (full-DR-crRNA) sequence with specific stem-loop G-C base substitutions exhibit increased editing efficiencies compared with the standard mature crRNA framework. Finally, we demonstrate in zebrafish embryos that the improved LbCas12a and FnoCas12a nucleases in combination with these modified crRNAs display high mutagenesis efficiencies and low toxicity when delivered as ribonucleoprotein complexes at high concentration. Together, these results define a set of enhanced Cas12a components with broad utility in vertebrate systems.

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is used to identify genome-wide DNA regions bound by proteins. Several sources of variation can affect the reproducibility of a particular ChIP-seq assay, which can lead to a misinterpretation of where the protein under investigation binds to the genome in a particular cell type. Given one ChIP-seq experiment with replicates, binding sites not observed in all the replicates will usually be interpreted as noise and discarded. However, the recent discovery of high-occupancy target (HOT) regions suggests that there are regions where binding of multiple transcription factors can be identified. To investigate these regions, we developed a reproducibility score and a method that identifies cell-specific variable regions in ChIP-seq data by integrating replicated ChIP-seq experiments for multiple protein targets on a particular cell type. Using our method, we found variable regions in human cell lines K562, GM12878, HepG2, MCF-7, and in mouse embryonic stem cells, defined as protein binding regions with non-reproducible results across replicated experiments. These variable-occupancy target (VOT) regions are CG dinucleotide rich, and show enrichment at promoters and R-loops. They overlap significantly with HOT regions, but are not blacklisted regions producing non-specific binding ChIP-seq peaks. Interestingly, among various genomic features, DNA accessibility is a better predictor of VOTs than CpG islands or epigenetic marks. Our method can be useful to point to such regions along the genome in a given cell type of interest, to improve the downstream interpretative analysis before follow up experiments.

High-throughput sequencing of RNAs isolated by crosslinking immunoprecipitation (HITS-CLIP, also called CLIP-Seq) has been used to map global RNA–protein interactions. However, a critical caveat of HITS-CLIP results is that they contain non-linear background noise—different extent of non-specific interactions caused by individual transcript abundance—that has been inconsiderately normalized, resulting in sacrifice of sensitivity. To properly deconvolute RNA–protein interactions, we have implemented CLIPick, a flexible peak calling pipeline for analyzing HITS-CLIP data, which statistically determines the signal-to-noise ratio for each transcript based on the expression-dependent background simulation. Comprising of streamlined Python modules with an easy-to-use standalone graphical user interface, CLIPick robustly identifies significant peaks and quantitatively defines footprint regions within which RNA–protein interactions were occurred. CLIPick outperforms other peak callers in accuracy and sensitivity, selecting the largest number of peaks particularly in lowly expressed transcripts where such marginal signals are hard to discriminate. Specifically, the application of CLIPick to Argonaute (Ago) HITS-CLIP data were sensitive enough to uncover extended features of microRNA target sites, and these sites were experimentally validated. CLIPick enables to resolve critical interactions in a wide spectrum of transcript levels and extends the scope of HITS-CLIP analysis. CLIPick is available at: http://clip.korea.ac.kr/clipick/

Genome editing of human induced pluripotent stem cells (iPSCs) is instrumental for functional genomics, disease modeling, and regenerative medicine. However, low editing efficiency has hampered the applications of CRISPR–Cas9 technology in creating knockin (KI) or knockout (KO) iPSC lines, which is largely due to massive cell death after electroporation with editing plasmids. Here, we report that the transient delivery of BCL-XL increases iPSC survival by ∼10-fold after plasmid transfection, leading to a 20- to 100-fold increase in homology-directed repair (HDR) KI efficiency and a 5-fold increase in non-homologous end joining (NHEJ) KO efficiency. Treatment with a BCL inhibitor ABT-263 further improves HDR efficiency by 70% and KO efficiency by 40%. The increased genome editing efficiency is attributed to higher expressions of Cas9 and sgRNA in surviving cells after electroporation. HDR or NHEJ efficiency reaches 95% with dual editing followed by selection of cells with HDR insertion of a selective gene. Moreover, KO efficiency of 100% can be achieved in a bulk population of cells with biallelic HDR KO followed by double selection, abrogating the necessity for single cell cloning. Taken together, these simple yet highly efficient editing strategies provide useful tools for applications ranging from manipulating human iPSC genomes to creating gene-modified animal models.

International audience; During cell division, maintenance of chromatin features from the parental genome requires their proper establishment on its newly synthetized copy. The loss of epigenetic marks within heterochromatin, typically enriched in repetitive elements, endangers genome stability and permits chromosomal rearrangements via recombination. However, how histone modifications associated with heterochromatin are maintained across mitosis remains poorly understood. KAP1 is known to act as a scaffold for a repressor complex that mediates local heterochromatin formation, and was previously demonstrated to play an important role during DNA repair. Accordingly, we investigated a putative role for this protein in the replication of heterochromatic regions. We first found that KAP1 associates with several DNA replication factors including PCNA, MCM3 and MCM6. We then observed that these interactions are promoted by KAP1 phosphorylation on serine 473 during S phase. Finally, we could demonstrate that KAP1 forms a complex with PCNA and the histone-lysine methyltransferase Suv39h1 to reinstate heterochromatin after DNA replication.

The two paralogous zinc finger factors CTCF and CTCFL differ in expression such that CTCF is ubiquitously expressed, whereas CTCFL is found during spermatogenesis and in some cancer types in addition to other cell types. Both factors share the highly conserved DNA binding domain and are bound to DNA sequences with an identical consensus. In contrast, both factors differ substantially in the number of bound sites in the genome. Here, we addressed the molecular features for this binding specificity. In contrast to CTCF we found CTCFL highly enriched at ‘open’ chromatin marked by H3K27 acetylation, H3K4 di- and trimethylation, H3K79 dimethylation and H3K9 acetylation plus the histone variant H2A.Z. CTCFL is enriched at transcriptional start sites and regions bound by transcription factors. Consequently, genes deregulated by CTCFL are highly cell specific. In addition to a chromatin-driven choice of binding sites, we determined nucleotide positions critical for DNA binding by CTCFL, but not by CTCF.

The PolyC binding proteins (PCBPs) impact alternative splicing of a subset of mammalian genes that are enriched in basic cellular functions. Here, we focus our analysis on PCBP-controlled cassette exon-splicing within the cell cycle control regulator cyclin-dependent kinase-2 (CDK2) transcript. We demonstrate that PCBP binding to a C-rich polypyrimidine tract (PPT) preceding exon 5 of the CDK2 transcript enhances cassette exon inclusion. This splice enhancement is U2AF65-independent and predominantly reflects actions of the PCBP1 isoform. Remarkably, PCBPs’ control of CDK2 ex5 splicing has evolved subsequent to mammalian divergence via conversion of constitutive exon 5 inclusion in the mouse CDK2 transcript to PCBP-responsive exon 5 alternative splicing in humans. Importantly, exclusion of exon 5 from the hCDK2 transcript dramatically represses the expression of CDK2 protein with a corresponding perturbation in cell cycle kinetics. These data highlight a recently evolved post-transcriptional pathway in primate species with the potential to modulate cell cycle control.

CRISPR systems have emerged as transformative tools for altering genomes in living cells with unprecedented ease, inspiring keen interest in increasing their specificity for perfectly matched targets. We have developed a novel approach for improving specificity by incorporating chemical modifications in guide RNAs (gRNAs) at specific sites in their DNA recognition sequence (‘guide sequence’) and systematically evaluating their on-target and off-target activities in biochemical DNA cleavage assays and cell-based assays. Our results show that a chemical modification (2′-O-methyl-3′-phosphonoacetate, or ‘MP’) incorporated at select sites in the ribose-phosphate backbone of gRNAs can dramatically reduce off-target cleavage activities while maintaining high on-target performance, as demonstrated in clinically relevant genes. These findings reveal a unique method for enhancing specificity by chemically modifying the guide sequence in gRNAs. Our approach introduces a versatile tool for augmenting the performance of CRISPR systems for research, industrial and therapeutic applications.

Conformation capture technologies measure frequencies of interactions between chromatin regions. However, understanding gene-regulation require knowledge of detailed spatial structures of heterogeneous chromatin in cells. Here we describe the nC-SAC (n-Constrained-Self Avoiding Chromatin) method that transforms experimental interaction frequencies into 3D ensembles of chromatin chains. nC-SAC first distinguishes specific from non-specific interaction frequencies, then generates 3D chromatin ensembles using identified specific interactions as spatial constraints. Application to α-globin locus shows that these constraints (∼20%) drive the formation of ∼99% all experimentally captured interactions, in which ∼30% additional to the imposed constraints is found to be specific. Many novel specific spatial contacts not captured by experiments are also predicted. A subset, of which independent ChIA-PET data are available, is validated to be RNAPII-, CTCF-, and RAD21-mediated. Their positioning in the architectural context of imposed specific interactions from nC-SAC is highly important. Our results also suggest the presence of a many-body structural unit involving α-globin gene, its enhancers, and POL3RK gene for regulating the expression of α-globin in silent cells.

Mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene cause cystic fibrosis (CF), but are not good predictors of lung phenotype. Genome-wide association studies (GWAS) previously identified additional genomic sites associated with CF lung disease severity. One of these, at chromosome 11p13, is an intergenic region between Ets homologous factor (EHF) and Apaf-1 interacting protein (APIP). Our goal was to determine the functional significance of this region, which being intergenic is probably regulatory. To identify cis-acting elements, we used DNase-seq and H3K4me1 and H3K27Ac ChIP-seq to map open and active chromatin respectively, in lung epithelial cells. Two elements showed strong enhancer activity for the promoters of EHF and the 5′ adjacent gene E47 like ETS transcription factor 5 (ELF5) in reporter gene assays. No enhancers of the APIP promoter were found. Circular chromosome conformation capture (4C-seq) identified direct physical interactions of elements within 11p13. This confirmed the enhancer-promoter associations, identified additional interacting elements and defined topologically associating domain (TAD) boundaries, enriched for CCCTC-binding factor (CTCF). No strong interactions were observed with the APIP promoter, which lies outside the main TAD encompassing the GWAS signal. These results focus attention on the role of EHF in modifying CF lung disease severity.

Chromatin structure affects DNA replication patterns, but the role of specific chromatin modifiers in regulating the replication process is yet unclear. We report that phosphorylation of the human SIRT1 deacetylase on Threonine 530 (T530-pSIRT1) modulates DNA synthesis. T530-pSIRT1 associates with replication origins and inhibits replication from a group of ‘dormant’ potential replication origins, which initiate replication only when cells are subject to replication stress. Although both active and dormant origins bind T530-pSIRT1, active origins are distinguished from dormant origins by their unique association with an open chromatin mark, histone H3 methylated on lysine 4. SIRT1 phosphorylation also facilitates replication fork elongation. SIRT1 T530 phosphorylation is essential to prevent DNA breakage upon replication stress and cells harboring SIRT1 that cannot be phosphorylated exhibit a high prevalence of extrachromosomal elements, hallmarks of perturbed replication. These observations suggest that SIRT1 phosphorylation modulates the distribution of replication initiation events to insure genomic stability.

TET1 oxidizes methylated cytosine into 5-hydroxymethylcytosine (5hmC), resulting in regulation of DNA methylation and gene expression. Full length TET1 (TET1FL) has a CXXC domain that binds to unmethylated CpG islands (CGIs). This CXXC domain allows TET1 to protect CGIs from aberrant methylation, but it also limits its ability to regulate genes outside of CGIs. Here, we report a novel isoform of TET1 (TET1ALT) that has a unique transcription start site from an alternate promoter in intron 2, yielding a protein with a unique translation start site. Importantly, TET1ALT lacks the CXXC domain but retains the catalytic domain. TET1ALT is repressed in embryonic stem cells (ESCs) but becomes activated in embryonic and adult tissues while TET1FL is expressed in ESCs, but repressed in adult tissues. Overexpression of TET1ALT shows production of 5hmC with distinct (and weaker) effects on DNA methylation or gene expression when compared to TET1FL. TET1ALT is aberrantly activated in multiple cancer types including breast, uterine and glioblastoma, and TET1 activation is associated with a worse overall survival in breast, uterine and ovarian cancers. Our data suggest that the predominantly activated isoform of TET1 in cancer cells does not protect from CGI methylation and likely mediates dynamic site-specific demethylation outside of CGIs.

The transcription factor c-Myb is involved in early differentiation and proliferation of haematopoietic cells, where it operates as a regulator of self-renewal and multi-lineage differentiation. Deregulated c-Myb plays critical roles in leukaemias and other human cancers. Due to its role as a master regulator, we hypothesized it might function as a pioneer transcription factor. Our approach to test this was to analyse a mutant of c-Myb, D152V, previously reported to cause haematopoietic defects in mice by an unknown mechanism. Our transcriptome data from K562 cells indicates that this mutation specifically affects c-Myb's ability to regulate genes involved in differentiation, causing failure in c-Myb's ability to block differentiation. Furthermore, we see a major effect of this mutation in assays where chromatin opening is involved. We show that each repeat in the minimal DNA-binding domain of c-Myb binds to histones and that D152V disrupts histone binding of the third repeat. ATAC-seq data indicates this mutation impairs the ability of c-Myb to cause chromatin opening at specific sites. Taken together, our findings support that c-Myb acts as a pioneer factor and show that D152V impairs this function. The D152V mutant is the first mutant of a transcription factor specifically destroying pioneer factor function.

One third of tumor suppressors are haploinsufficient transcriptional regulators, yet it remains unknown how a 50% reduction of a transcription factor is translated at the cis-regulatory level into a malignant transcriptional program. We studied CUX1, a haploinsufficient transcription factor that is recurrently mutated in hematopoietic and solid tumors. We determined CUX1 DNA-binding and target gene regulation in the wildtype and haploinsufficient states. CUX1 binds with transcriptional activators and cohesin at distal enhancers across three different human cell types. Haploinsufficiency of CUX1 altered the expression of a large number of genes, including cell cycle regulators, with concomitant increased cellular proliferation. Surprisingly, CUX1 occupancy decreased genome-wide in the haploinsufficient state, and binding site affinity did not correlate with differential gene expression. Instead, differentially expressed genes had multiple, low-affinity CUX1 binding sites, features of analog gene regulation. A machine-learning algorithm determined that chromatin accessibility, enhancer activity, and distance to the transcription start site are features of dose-sensitive CUX1 transcriptional regulation. Moreover, CUX1 is enriched at sites of DNA looping, as determined by Hi-C analysis, and these loops connect CUX1 to the promoters of regulated genes. We propose an analog model for haploinsufficient transcriptional deregulation mediated by higher order genome architecture.

The XIST RNA is a non-coding RNA that induces X chromosome inactivation (XCI). Unlike the mouse Xist RNA, how the human XIST RNA controls XCI in female cells is less well characterized, and its functional motifs remain unclear. To systematically decipher the XCI-involving elements of XIST RNA, 11 smaller XIST segments, including repeats A, D and E; human-specific repeat elements; the promoter; and non-repetitive exons, as well as the entire XIST gene, were homozygously deleted in K562 cells using the Cas9 nuclease and paired guide RNAs at high efficiencies, followed by high-throughput RNA sequencing and RNA fluorescence in situ hybridization experiments. Clones containing en bloc and promoter deletions that consistently displayed no XIST RNAs and a global up-regulation of X-linked genes confirmed that the deletion of XIST reactivates the inactive X chromosome. Systematic analyses of segmental deletions delineated that exon 5 harboring the non-repeat element is important for X-inactivation maintenance, whereas exons 2, 3 and 4 as well as the other repeats in exon 1 are less important, a different situation from that of mouse Xist. This Cas9-assisted dissection of XIST allowed us to understand the unique functional domains within the human XIST RNA.

textabstractNeuroglobin (NGB) is predominantly expressed in the brain and retina. Studies suggest that NGB exerts protective effects to neuronal cells and is implicated in reducing the severity of stroke and Alzheimer's disease. However, little is known about the mechanisms which regulate the cell type-specific expression of the gene. In this study, we hypothesized that distal regulatory elements (DREs) are involved in optimal expression of the NGB gene. By chromosome conformation capture we identified two novel DREs located -70 kb upstream and +100 kb downstream from the NGB gene. ENCODE database showed the presence of DNaseI hypersensitive and transcription factors binding sites in these regions. Further analyses using luciferase reporters and chromatin immunoprecipitation suggested that the -70 kb region upstream of the NGB gene contained a neuronalspecific enhancer and GATA transcription factor binding sites. Knockdown of GATA-2 caused NGB expression to drop dramatically, indicating GATA-2 as an essential transcription factor for the activation of NGB expression. The crucial role of the DRE in NGB expression activation was further confirmed by the drop in NGB level after CRISPR-mediated deletion of the DRE. Taken together, we show that the NGB gene is regulated by a cell type-specific loop formed between its promoter and the novel DRE.

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 3-dimensional folding of a 1 Mbp region of human chromosome 11 containing the {\beta}-globin genes by integrating looping interactions of the 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 {\beta}-globin locus control region (LCR) and genes. In these cells, the modeled {\beta}-globin domain folds into a globule with the LCR and the active globin genes on the periphery. By 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 positioning 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., Comment: Full article, including Supp. Mat., is available at Nucleic Acids Research, doi: 10.1093/nar/gks536

A number of established and investigational anticancer drugs slow the religation step of DNA topoisomerase I (topo I). These agents induce cytotoxicity by stabilizing topo I-DNA covalent complexes, which in turn interact with advancing replication forks or transcription complexes to generate lethal lesions. Despite the importance of topo I-DNA covalent complexes, it has been difficult to detect these lesions within intact cells and tumors. Here, we report development of a monoclonal antibody that specifically recognizes covalent topo I-DNA complexes, but not free topo I or DNA, by immunoblotting, immunofluorescence or flow cytometry. Utilizing this antibody, we demonstrate readily detectable topo I-DNA covalent complexes after treatment with camptothecins, indenoisoquinolines and cisplatin but not nucleoside analogues. Topotecan-induced topo I-DNA complexes peak at 15-30 min after drug addition and then decrease, whereas indotecan-induced complexes persist for at least 4 h. Interestingly, simultaneous staining for covalent topo I-DNA complexes, phospho-H2AX and Rad51 suggests that topotecan-induced DNA double-strand breaks occur at sites distinct from stabilized topo I-DNA covalent complexes. These studies not only provide new insight into the action of topo I-directed agents, but also illustrate a strategy that can be applied to study additional topoisomerases and their inhibitors in vitro and in vivo.

DNA strand asymmetries can have a major effect on several biological functions, including replication, transcription and transcription factor binding. As such, DNA strand asymmetries and mutational strand bias can provide information about biological function. However, a versatile tool to explore this does not exist. Here, we present Asymmetron, a user-friendly computational tool that performs statistical analysis and visualizations for the evaluation of strand asymmetries. Asymmetron takes as input DNA features provided with strand annotation and outputs strand asymmetries for consecutive occurrences of a single DNA feature or between pairs of features. We illustrate the use of Asymmetron by identifying transcriptional and replicative strand asymmetries of germline structural variant breakpoints. We also show that the orientation of the binding sites of 45% of human transcription factors analyzed have a significant DNA strand bias in transcribed regions, that is also corroborated in ChIP-seq analyses, and is likely associated with transcription. In summary, we provide a novel tool to assess DNA strand asymmetries and show how it can be used to derive new insights across a variety of biological disciplines.

Eukaryotic transcriptomes are complex, involving thousands of overlapping transcripts. The interleaved nature of the transcriptomes limits our ability to identify regulatory regions, and in some cases can lead to misinterpretation of gene expression. To improve the understanding of the overlapping transcriptomes, we have developed an optimized method, TIF-Seq2, able to sequence simultaneously the 5′ and 3′ ends of individual RNA molecules at single-nucleotide resolution. We investigated the transcriptome of a well characterized human cell line (K562) and identified thousands of unannotated transcript isoforms. By focusing on transcripts which are challenging to be investigated with RNA-Seq, we accurately defined boundaries of lowly expressed unannotated and read-through transcripts putatively encoding fusion genes. We validated our results by targeted long-read sequencing and standard RNA-Seq for chronic myeloid leukaemia patient samples. Taking the advantage of TIF-Seq2, we explored transcription regulation among overlapping units and investigated their crosstalk. We show that most overlapping upstream transcripts use poly(A) sites within the first 2 kb of the downstream transcription units. Our work shows that, by paring the 5′ and 3′ end of each RNA, TIF-Seq2 can improve the annotation of complex genomes, facilitate accurate assignment of promoters to genes and easily identify transcriptionally fused genes.

To define the functions of NCOR1 we developed an integrative analysis that combined ENCODE and NCI-60 data, followed by in vitro validation. NCOR1 and H3K9me3 ChIP-Seq, FAIRE-seq and DNA CpG methylation interactions were related to gene expression using bootstrapping approaches. Most NCOR1 combinations (24/44) were associated with significantly elevated level expression of protein coding genes and only very few combinations related to gene repression. DAVID's biological process annotation revealed that elevated gene expression was uniquely associated with acetylation and ETS binding. A matrix of gene and drug interactions built on NCI-60 data identified that Imatinib significantly targeted the NCOR1 governed transcriptome. Stable knockdown of NCOR1 in K562 cells slowed growth and significantly repressed genes associated with NCOR1 cistrome, again, with the GO terms acetylation and ETS binding, and significantly dampened sensitivity to Imatinib-induced erythroid differentiation. Mining public microarray data revealed that NCOR1-targeted genes were significantly enriched in Imatinib response gene signatures in cell lines and chronic myelogenous leukemia (CML) patients. These approaches integrated cistrome, transcriptome and drug sensitivity relationships to reveal that NCOR1 function is surprisingly most associated with elevated gene expression, and that these targets, both in CML cell lines and patients, associate with sensitivity to Imatinib.

Polypyrimidine tract binding protein (PTBP1) is a widely expressed RNA binding protein that acts as a regulator of alternative splicing and of cytoplasmic mRNA functions. Vertebrates contain two closely-related paralogs with >75% amino acid sequence identity. Early replacement of PTBP1 by PTBP2 during neuronal differentiation causes a concerted set of splicing changes. By comparison, very little is known about the molecular functions or physiological roles of PTBP3, although its expression and conservation throughout the vertebrates suggest a role in haematopoietic cells. To begin to understand its functions we have characterized the mRNA and protein isoform repertoire of PTBP3. Combinatorial alternative splicing events at the 5' end of the gene allow for the generation of eight mRNA and three major protein isoforms. Individual mRNAs generate up to three protein isoforms via alternative translation initiation by re-initiation and leaky scanning using downstream AUG codons. The N-terminally truncated PTBP3 isoforms lack nuclear localization signals and/or most of the RRM1 domain and vary in their RNA binding properties and nuclear/cytoplasmic distribution, suggesting that PTBP3 may have major post-transcriptional cytoplasmic roles. Our findings set the stage for understanding the non-redundant physiological roles of PTBP3.

Lentiviral vectors almost universally use heterologous internal promoters to express transgenes. One of the most commonly used promoter fragments is a 1.2-kb sequence from the human ubiquitin C (UBC) gene, encompassing the promoter, some enhancers, first exon, first intron and a small part of the second exon of UBC. Because splicing can occur after transcription of the vector genome during vector production, we investigated whether the intron within the UBC promoter fragment is faithfully transmitted to target cells. Genetic analysis revealed that more than 80% of proviral forms lack the intron of the UBC promoter. The human elongation factor 1 alpha (EEF1A1) promoter fragment intron was not lost during lentiviral packaging, and this difference between the UBC and EEF1A1 promoter introns was conferred by promoter exonic sequences. UBC promoter intron loss caused a 4-fold reduction in transgene expression. Movement of the expression cassette to the opposite strand prevented intron loss and restored full expression. This increase in expression was mostly due to non-classical enhancer activity within the intron, and movement of putative intronic enhancer sequences to multiple promoter-proximal sites actually repressed expression. Reversal of the UBC promoter also prevented intron loss and restored full expression in bidirectional lentiviral vectors.

Most retroviruses preferentially integrate into certain genomic locations and, as a result, their genome-wide integration patterns are non-random. We investigate the epigenetic landscape of integrated retroviral vectors and correlate it with the long-term stability of proviral transcription. Retroviral vectors derived from the avian sarcoma/leukosis virus expressing the GFP reporter were used to transduce the human myeloid lymphoblastoma cell line K562. Because of efficient silencing of avian retrovirus in mammalian cells, only ∼3% of established clones displayed stable proviral expression. We analyzed the vector integration sites in non-selected cells and in clones selected for the GFP expression. This selection led to overrepresentation of proviruses integrated in active transcription units, with particular accumulation in promoter-proximal areas. In parallel, we investigated the integration of vectors equipped with an anti-silencing CpG island core sequence. Such modification increased the frequency of stably expressing proviruses by one order. The modified vectors are also overrepresented in active transcription units, but stably expressed in distal parts of transcriptional units further away from promoters with marked accumulation in enhancers. These results suggest that integrated retroviruses subject to gradual epigenetic silencing during long-term cultivation. Among most genomic compartments, however, active promoters and enhancers protect the adjacent retroviruses from transcriptional silencing.

The initiation of mRNA translation from start codons other than AUG was previously believed to be rare and of relatively low impact. More recently, evidence has suggested that as much as half of all translation initiation utilizes non-AUG start codons, codons that deviate from AUG by a single base. Furthermore, non-AUG start codons have been shown to be involved in regulation of expression and disease etiology. Yet the ability to gauge expression based on the sequence of a translation initiation site (start codon and its flanking bases) has been limited. Here we have performed a comprehensive analysis of translation initiation sites that utilize non-AUG start codons. By combining genetic-reporter, cell-sorting, and high-throughput sequencing technologies, we have analyzed the expression associated with all possible variants of the -4 to +4 positions of non-AUG translation initiation site motifs. This complete motif analysis revealed that 1) with the right sequence context, certain non-AUG start codons can generate expression comparable to that of AUG start codons, 2) sequence context affects each non-AUG start codon differently, and 3) initiation at non-AUG start codons is highly sensitive to changes in the flanking sequences. Complete motif analysis has the potential to be a key tool for experimental and diagnostic genomics.

Proximal promoter regions (PPR) are heavily transcribed yielding different types of small RNAs. The act of transcription within PPRs might regulate downstream gene expression via transcriptional interference (TI). For analysis, we investigated capped and polyadenylated small RNA transcripts within PPRs of human RefSeq genes in eight different cell lines. Transcripts of our datasets overlapped with experimentally determined transcription factor binding sites (TFBS). For TFBSs intersected by these small RNA transcripts, we established negative correlation of sRNA expression levels and transcription factor (TF) DNA binding affinities; suggesting that the transcripts acted via TI. Accordingly, datasets were designated as TFbiTrs (TF-binding interfering transcripts). Expression of most TFbiTrs was restricted to certain cell lines. This facilitated the analysis of effects related to TFbiTr expression for the same RefSeq genes across cell lines. We consistently uncovered higher relative TF/DNA binding affinities and concomitantly higher expression levels for RefSeq genes in the absence of TFbiTrs. Analysis of corresponding chromatin landscapes supported these results. ChIA-PET revealed the participation of distal enhancers in TFbiTr transcription. Enhancers regulating TFbiTrs, in effect, act as repressors for corresponding downstream RefSeq genes. We demonstrate the significant impact of TI on gene expression using selected small RNA datasets.

The ubiquitously expressed transcription factor TFII-I exerts both positive and negative effects on transcription. Using biotinylation tagging technology and high-throughput sequencing, we determined sites of chromatin interactions for TFII-I in the human erythroleukemia cell line K562. This analysis revealed that TFII-I binds upstream of the transcription start site of expressed genes, both upstream and downstream of the transcription start site of repressed genes, and downstream of RNA polymerase II peaks at the ATF3 and other stress responsive genes. At the ATF3 gene, TFII-I binds immediately downstream of a Pol II peak located 5 kb upstream of exon 1. Induction of ATF3 expression increases transcription throughout the ATF3 gene locus which requires TFII-I and correlates with increased association of Pol II and Elongin A. Pull-down assays demonstrated that TFII-I interacts with Elongin A. Partial depletion of TFII-I expression caused a reduction in the association of Elongin A with and transcription of the DNMT1 and EFR3A genes without a decrease in Pol II recruitment. The data reveal different interaction patterns of TFII-I at active, repressed, or inducible genes, identify novel TFII-I interacting proteins, implicate TFII-I in the regulation of transcription elongation and provide insight into the role of TFII-I during the response to cellular stress.