G-quadruplexes (G4s) are a set of stable secondary structures that form within guanine-rich regions of single-stranded nucleic acids that pose challenges for DNA maintenance. The G-rich DNA sequence at telomeres has a propensity to form G4s of various topologies. The human protein complexes Replication Protein A (RPA) and CTC1-STN1-TEN1 (CST) are implicated in managing G4s at telomeres, leading to DNA unfolding and allowing telomere replication to proceed. Here, we use fluorescence anisotropy equilibrium binding measurements to determine the ability of these proteins to bind various telomeric G4s. We find that the ability of CST to specifically bind G-rich ssDNA is substantially inhibited by the presence of G4s. In contrast, RPA tightly binds telomeric G4s, showing negligible changes in affinity for G4 structure compared to linear ssDNAs. Using a mutagenesis strategy, we found that RPA DNA-binding domains work together for G4 binding, and simultaneous disruption of these domains reduces the affinity of RPA for G4 ssDNA. The relative inability of CST to disrupt G4s, combined with the greater cellular abundance of RPA, suggests that RPA could act as a primary protein complex responsible for resolving G4s at telomeres.

Phosphorothioates (PS) have proven their effectiveness in the area of therapeutic oligonucleotides with applications spanning from cancer treatment to neurodegenerative disorders. Initially, PS substitution was introduced for the antisense oligonucleotides (PS ASOs) because it confers an increased nuclease resistance meanwhile ameliorates cellular uptake and in-vivo bioavailability. Thus, PS oligonucleotides have been elevated to a fundamental asset in the realm of gene silencing therapeutic methodologies. But, despite their wide use, little is known on the possibly different structural changes PS-substitutions may provoke in DNA·RNA hybrids. Additionally, scarce information and significant controversy exists on the role of phosphorothioate chirality in modulating PS properties. Here, through comprehensive computational investigations and experimental measurements, we shed light on the impact of PS chirality in DNA-based antisense oligonucleotides; how the different phosphorothioate diastereomers impact DNA topology, stability and flexibility to ultimately disclose pro-Sp S and pro-Rp S roles at the catalytic core of DNA Exonuclease and Human Ribonuclease H; two major obstacles in ASOs-based therapies. Altogether, our results provide full-atom and mechanistic insights on the structural aberrations PS-substitutions provoke and explain the origin of nuclease resistance PS-linkages confer to DNA·RNA hybrids; crucial information to improve current ASOs-based therapies.

RNA-binding proteins (RBPs) form highly diverse and dynamic ribonucleoprotein complexes, whose functions determine the molecular fate of the bound RNA. In the model organism Sacchromyces cerevisiae, the number of proteins identified as RBPs has greatly increased over the last decade. However, the cellular function of most of these novel RBPs remains largely unexplored. We used mass spectrometry-based quantitative proteomics to systematically identify protein-protein interactions (PPIs) and RNA-dependent interactions (RDIs) to create a novel dataset for 40 RBPs that are associated with the mRNA life cycle. Domain, functional and pathway enrichment analyses revealed an over-representation of RNA functionalities among the enriched interactors. Using our extensive PPI and RDI networks, we revealed putative new members of RNA-associated pathways, and highlighted potential new roles for several RBPs. Our RBP interactome resource is available through an online interactive platform as a community tool to guide further in-depth functional studies and RBP network analysis (https://www.butterlab.org/RINE).

Mosaic mutations can be used to track cell ancestries and reconstruct high-resolution lineage trees during cancer progression and during development, starting from the first cell divisions of the zygote. However, this approach requires sampling and analyzing the genomes of multiple cells, which can be redundant in lineage representation, limiting the scalability of the approach. We describe a strategy for cost- and time-efficient lineage reconstruction using clonal induced pluripotent stem cell lines from human skin fibroblasts. The approach leverages shallow sequencing coverage to assess the clonality of the lines, clusters redundant lines and sums their coverage to accurately discover mutations in the corresponding lineages. Only a fraction of lines needs to be sequenced to high coverage. We demonstrate the effectiveness of this approach for reconstructing lineage trees during development and in hematologic malignancies. We discuss and propose an optimal experimental design for reconstructing lineage trees.

CST (CTC1-STN1-TEN1) is a telomere associated complex that binds ssDNA and is required for multiple steps in telomere replication, including termination of G-strand extension by telomerase and synthesis of the complementary C-strand. CST contains seven OB-folds which appear to mediate CST function by modulating CST binding to ssDNA and the ability of CST to recruit or engage partner proteins. However, the mechanism whereby CST achieves its various functions remains unclear. To address the mechanism, we generated a series of CTC1 mutants and studied their effect on CST binding to ssDNA and their ability to rescue CST function in CTC1-/- cells. We identified the OB-B domain as a key determinant of telomerase termination but not C-strand synthesis. CTC1-ΔB expression rescued C-strand fill-in, prevented telomeric DNA damage signaling and growth arrest. However, it caused progressive telomere elongation and the accumulation of telomerase at telomeres, indicating an inability to limit telomerase action. The CTC1-ΔB mutation greatly reduced CST-TPP1 interaction but only modestly affected ssDNA binding. OB-B point mutations also weakened TPP1 association, with the deficiency in TPP1 interaction tracking with an inability to limit telomerase action. Overall, our results indicate that CTC1-TPP1 interaction plays a key role in telomerase termination.

Unveiling the nucleic acid binding sites of a protein helps reveal its regulatory functions in vivo. Current methods encode protein sites from the handcrafted features of their local neighbors and recognize them via a classification, which are limited in expressive ability. Here, we present GeoBind, a geometric deep learning method for predicting nucleic binding sites on protein surface in a segmentation manner. GeoBind takes the whole point clouds of protein surface as input and learns the high-level representation based on the aggregation of their neighbors in local reference frames. Testing GeoBind on benchmark datasets, we demonstrate GeoBind is superior to state-of-the-art predictors. Specific case studies are performed to show the powerful ability of GeoBind to explore molecular surfaces when deciphering proteins with multimer formation. To show the versatility of GeoBind, we further extend GeoBind to five other types of ligand binding sites prediction tasks and achieve competitive performances.

The long non-coding RNA EPR is expressed in epithelial tissues, binds to chromatin and controls distinct biological activities in mouse mammary gland cells. Because of its high expression in the intestine, in this study we have generated a colon-specific conditional targeted deletion (EPR cKO) to evaluate EPR in vivo functions in mice. EPR cKO mice display epithelium hyperproliferation, impaired mucus production and secretion, as well as inflammatory infiltration in the proximal portion of the large intestine. RNA sequencing analysis reveals a rearrangement of the colon crypt transcriptome with strong reduction of goblet cell-specific factors including those involved in the synthesis, assembly, transport and control of mucus proteins. Further, colon mucosa integrity and permeability are impaired in EPR cKO mice, and this results in higher susceptibility to dextran sodium sulfate (DSS)-induced colitis and tumor formation. Human EPR is down-regulated in human cancer cell lines as well as in human cancers, and overexpression of EPR in a colon cancer cell line results in enhanced expression of pro-apoptotic genes. Mechanistically, we show that EPR directly interacts with select genes involved in mucus metabolism whose expression is reduced in EPR cKO mice and that EPR deletion causes tridimensional chromatin organization changes.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of coronavirus disease 2019 (COVID-19). The NSP15 endoribonuclease enzyme, known as NendoU, is highly conserved and plays a critical role in the ability of the virus to evade the immune system. NendoU is a promising target for the development of new antiviral drugs. However, the complexity of the enzyme's structure and kinetics, along with the broad range of recognition sequences and lack of structural complexes, hampers the development of inhibitors. Here, we performed enzymatic characterization of NendoU in its monomeric and hexameric form, showing that hexamers are allosteric enzymes with a positive cooperative index, and with no influence of manganese on enzymatic activity. Through combining cryo-electron microscopy at different pHs, X-ray crystallography and biochemical and structural analysis, we showed that NendoU can shift between open and closed forms, which probably correspond to active and inactive states, respectively. We also explored the possibility of NendoU assembling into larger supramolecular structures and proposed a mechanism for allosteric regulation. In addition, we conducted a large fragment screening campaign against NendoU and identified several new allosteric sites that could be targeted for the development of new inhibitors. Overall, our findings provide insights into the complex structure and function of NendoU and offer new opportunities for the development of inhibitors.

Targeted in vivo hypermutation accelerates directed evolution of proteins through concurrent DNA diversification and selection. Although systems employing a fusion protein of a nucleobase deaminase and T7 RNA polymerase present gene-specific targeting, their mutational spectra have been limited to exclusive or dominant C:G→T:A mutations. Here we describe eMutaT7transition, a new gene-specific hypermutation system, that installs all transition mutations (C:G→T:A and A:T→G:C) at comparable frequencies. By using two mutator proteins in which two efficient deaminases, PmCDA1 and TadA-8e, are separately fused to T7 RNA polymerase, we obtained similar numbers of C:G→T:A and A:T→G:C substitutions at a sufficiently high frequency (∼6.7 substitutions in 1.3 kb gene during 80-h in vivo mutagenesis). Through eMutaT7transition-mediated TEM-1 evolution for antibiotic resistance, we generated many mutations found in clinical isolates. Overall, with a high mutation frequency and wider mutational spectrum, eMutaT7transition is a potential first-line method for gene-specific in vivo hypermutation.

A large number of transcription factors have been shown to bind and interact with mitotic chromosomes, which may promote the efficient reactivation of transcriptional programs following cell division. Although the DNA-binding domain (DBD) contributes strongly to TF behavior, the mitotic behaviors of TFs from the same DBD family may vary. To define the mechanisms governing TF behavior during mitosis in mouse embryonic stem cells, we examined two related TFs: Heat Shock Factor 1 and 2 (HSF1 and HSF2). We found that HSF2 maintains site-specific binding genome-wide during mitosis, whereas HSF1 binding is somewhat decreased. Surprisingly, live-cell imaging shows that both factors appear excluded from mitotic chromosomes to the same degree, and are similarly more dynamic in mitosis than in interphase. Exclusion from mitotic DNA is not due to extrinsic factors like nuclear import and export mechanisms. Rather, we found that the HSF DBDs can coat mitotic chromosomes, and that HSF2 DBD is able to establish site-specific binding. These data further confirm that site-specific binding and chromosome coating are independent properties, and that for some TFs, mitotic behavior is largely determined by the non-DBD regions.

Apurinic/apyrimidinic (AP) sites are DNA lesions created under normal growth conditions that result in cytotoxicity, replication-blocks, and mutations. AP sites are susceptible to β-elimination and are liable to be converted to DNA strand breaks. HMCES (5-hydroxymethylcytosine binding, ES cell specific) protein interacts with AP sites in single stranded (ss) DNA exposed at DNA replication forks to generate a stable thiazolidine protein-DNA crosslink and protect cells against AP site toxicity. The crosslinked HMCES is resolved by proteasome-mediated degradation; however, it is unclear how HMCES-crosslinked ssDNA and the resulting proteasome-degraded HMCES adducts are processed and repaired. Here, we describe methods for the preparation of thiazolidine adduct-containing oligonucleotides and determination of their structure. We demonstrate that the HMCES-crosslink is a strong replication blocking adduct and that protease-digested HMCES adducts block DNA replication to a similar extent as AP sites. Moreover, we show that the human AP endonuclease APE1 incises DNA 5' to the protease-digested HMCES adduct. Interestingly, while HMCES-ssDNA crosslinks are stable, the crosslink is reversed upon the formation of dsDNA, possibly due to a catalytic reverse reaction. Our results shed new light on damage tolerance and repair pathways for HMCES-DNA crosslinks in human cells.

How multi-segmented double-stranded RNA (dsRNA) viruses correctly incorporate their genomes into their capsids remains unclear for many viruses, including Bluetongue virus (BTV), a Reoviridae member, with a genome of 10 segments. To address this, we used an RNA-cross-linking and peptide-fingerprinting assay (RCAP) to identify RNA binding sites of the inner capsid protein VP3, the viral polymerase VP1 and the capping enzyme VP4. Using a combination of mutagenesis, reverse genetics, recombinant proteins and in vitro assembly, we validated the importance of these regions in virus infectivity. Further, to identify which RNA segments and sequences interact with these proteins, we used viral photo-activatable ribonucleoside crosslinking (vPAR-CL) which revealed that the larger RNA segments (S1-S4) and the smallest segment (S10) have more interactions with viral proteins than the other smaller segments. Additionally, using a sequence enrichment analysis we identified an RNA motif of nine bases that is shared by the larger segments. The importance of this motif for virus replication was confirmed by mutagenesis followed by virus recovery. We further demonstrated that these approaches could be applied to a related Reoviridae member, rotavirus (RV), which has human epidemic impact, offering the possibility of novel intervention strategies for a human pathogen.

Most cell type-specific genes are regulated by the interaction of enhancers with their promoters. The identification of enhancers is not trivial as enhancers are diverse in their characteristics and dynamic in their interaction partners. We present Esearch3D, a new method that exploits network theory approaches to identify active enhancers. Our work is based on the fact that enhancers act as a source of regulatory information to increase the rate of transcription of their target genes and that the flow of this information is mediated by the folding of chromatin in the three-dimensional (3D) nuclear space between the enhancer and the target gene promoter. Esearch3D reverse engineers this flow of information to calculate the likelihood of enhancer activity in intergenic regions by propagating the transcription levels of genes across 3D genome networks. Regions predicted to have high enhancer activity are shown to be enriched in annotations indicative of enhancer activity. These include: enhancer-associated histone marks, bidirectional CAGE-seq, STARR-seq, P300, RNA polymerase II and expression quantitative trait loci (eQTLs). Esearch3D leverages the relationship between chromatin architecture and transcription, allowing the prediction of active enhancers and an understanding of the complex underpinnings of regulatory networks. The method is available at: https://github.com/InfOmics/Esearch3D and https://doi.org/10.5281/zenodo.7737123.

Gene expression changes are orchestrated by transcription factors (TFs), which bind to DNA to regulate gene expression. It remains surprisingly difficult to predict basic features of the transcriptional process, including in vivo TF occupancy. Existing thermodynamic models of TF function are often not concordant with experimental measurements, suggesting undiscovered biology. Here, we analyzed one of the most well-studied TFs, the yeast zinc cluster Gal4, constructed a Shea-Ackers thermodynamic model to describe its binding, and compared the results of this model to experimentally measured Gal4p binding in vivo. We found that at many promoters, the model predicted no Gal4p binding, yet substantial binding was observed. These outlier promoters lacked canonical binding motifs, and subsequent investigation revealed Gal4p binds unexpectedly to DNA sequences with high densities of its half site (CGG). We confirmed this novel mode of binding through multiple experimental and computational paradigms; we also found most other zinc cluster TFs we tested frequently utilize this binding mode, at 27% of their targets on average. Together, these results demonstrate a novel mode of binding where zinc clusters, the largest class of TFs in yeast, bind DNA sequences with high densities of half sites.

Prokaryotic Argonaute proteins (pAgos) are homologs of eukaryotic Argonautes (eAgos) and are also thought to play a role in cell defense against invaders. However, pAgos are much more diverse than eAgos and little is known about their functional activities and target specificities in vivo. Here, we describe five pAgos from mesophilic bacteria that act as programmable DNA endonucleases and analyze their ability to target chromosomal and invader DNA. In vitro, the analyzed proteins use small guide DNAs for precise cleavage of single-stranded DNA at a wide range of temperatures. Upon their expression in Escherichia coli, all five pAgos are loaded with small DNAs preferentially produced from plasmids and chromosomal regions of replication termination. One of the tested pAgos, EmaAgo from Exiguobacterium marinum, can induce DNA interference between homologous sequences resulting in targeted processing of multicopy plasmid and genomic elements. EmaAgo also protects bacteria from bacteriophage infection, by loading phage-derived guide DNAs and decreasing phage DNA content and phage titers. Thus, the ability of pAgos to target multicopy elements may be crucial for their protective function. The wide spectrum of pAgo activities suggests that they may have diverse functions in vivo and paves the way for their use in biotechnology.

Ribosome biogenesis occurs co-transcriptionally and entails rRNA folding, ribosomal protein binding, rRNA processing, and rRNA modification. In most bacteria, the 16S, 23S and 5S rRNAs are co-transcribed, often with one or more tRNAs. Transcription involves a modified RNA polymerase, called the antitermination complex, which forms in response to cis-acting elements (boxB, boxA and boxC) in the nascent pre-rRNA. Sequences flanking the rRNAs are complementary and form long helices known as leader-trailer helices. Here, we employed an orthogonal translation system to interrogate the functional roles of these RNA elements in 30S subunit biogenesis in Escherichia coli. Mutations that disrupt the leader-trailer helix caused complete loss of translation activity, indicating that this helix is absolutely essential for active subunit formation in the cell. Mutations of boxA also reduced translation activity, but by only 2- to 3-fold, suggesting a smaller role for the antitermination complex. Similarly modest drops in activity were seen upon deletion of either or both of two leader helices, termed here hA and hB. Interestingly, subunits formed in the absence of these leader features exhibited defects in translational fidelity. These data suggest that the antitermination complex and precursor RNA elements help to ensure quality control during ribosome biogenesis.

Conversely to canonical splicing, back-splicing connects the upstream 3' splice site (SS) with a downstream 5'SS and generates exonic circular RNAs (circRNAs) that are widely identified and have regulatory functions in eukaryotic gene expression. However, sex-specific back-splicing in Drosophila has not been investigated and its regulation remains unclear. Here, we performed multiple RNA analyses of a variety sex-specific Drosophila samples and identified over ten thousand circular RNAs, in which hundreds are sex-differentially and -specifically back-spliced. Intriguingly, we found that expression of SXL, an RNA-binding protein encoded by Sex-lethal (Sxl), the master Drosophila sex-determination gene that is only spliced into functional proteins in females, promoted back-splicing of many female-differential circRNAs in the male S2 cells, whereas expression of a SXL mutant (SXLRRM) did not promote those events. Using a monoclonal antibody, we further obtained the transcriptome-wide RNA-binding sites of SXL through PAR-CLIP. After splicing assay of mini-genes with mutations in the SXL-binding sites, we revealed that SXL-binding on flanking exons and introns of pre-mRNAs facilitates back-splicing, whereas SXL-binding on the circRNA exons inhibits back-splicing. This study provides strong evidence that SXL has a regulatory role in back-splicing to generate sex-specific and -differential circRNAs, as well as in the initiation of sex-determination cascade through canonical forward-splicing.

Liquid-liquid phase separation (LLPS) plays a critical role in regulating gene transcription via the formation of transcriptional condensates. However, LLPS has not been reported to be engineered as a tool to activate endogenous gene expression in mammalian cells or in vivo. Here, we developed a droplet-forming CRISPR (clustered regularly interspaced short palindromic repeats) gene activation system (DropCRISPRa) to activate transcription with high efficiency via combining the CRISPR-SunTag system with FETIDR-AD fusion proteins, which contain an N-terminal intrinsically disordered region (IDR) of a FET protein (FUS or TAF15) and a transcription activation domain (AD, VP64/P65/VPR). In this system, the FETIDR-AD fusion protein formed phase separation condensates at the target sites, which could recruit endogenous BRD4 and RNA polymerase II with an S2 phosphorylated C-terminal domain (CTD) to enhance transcription elongation. IDR-FUS9Y>S and IDR-FUSG156E, two mutants with deficient and aberrant phase separation respectively, confirmed that appropriate phase separation was required for efficient gene activation. Further, the DropCRISPRa system was compatible with a broad set of CRISPR-associated (Cas) proteins and ADs, including dLbCas12a, dAsCas12a, dSpCas9 and the miniature dUnCas12f1, and VP64, P65 and VPR. Finally, the DropCRISPRa system could activate target genes in mice. Therefore, this study provides a robust tool to activate gene expression for foundational research and potential therapeutics.

Mutational signatures discerned in cancer genomes, in aging tissues and in cells exposed to toxic agents, reflect complex processes underlying transformation of cells from normal to dysfunctional. Due to its ubiquitous and chronic nature, redox stress contributions to cellular makeover remain equivocal. The deciphering of a new mutational signature of an environmentally-relevant oxidizing agent, potassium bromate, in yeast single strand DNA uncovered a surprising heterogeneity in the mutational signatures of oxidizing agents. NMR-based analysis of molecular outcomes of redox stress revealed profound dissimilarities in metabolic landscapes following exposure to hydrogen peroxide versus potassium bromate. The predominance of G to T substitutions in the mutational spectra distinguished potassium bromate from hydrogen peroxide and paraquat and mirrored the observed metabolic changes. We attributed these changes to the generation of uncommon oxidizing species in a reaction with thiol-containing antioxidants; a nearly total depletion of intracellular glutathione and a paradoxical augmentation of potassium bromate mutagenicity and toxicity by antioxidants. Our study provides the framework for understanding multidimensional processes triggered by agents collectively known as oxidants. Detection of increased mutational loads associated with potassium bromate-related mutational motifs in human tumors may be clinically relevant as a biomarker of this distinct type of redox stress.

Bacteria use a diverse arsenal of anti-phage immune systems, including CRISPR-Cas and restriction enzymes. Recent advances in anti-phage system discovery and annotation tools have unearthed many unique systems, often encoded in horizontally transferred defense islands, which can be horizontally transferred. Here, we developed Hidden Markov Models (HMMs) for defense systems and queried microbial genomes on the NCBI database. Out of the 30 species with >200 completely sequenced genomes, our analysis found Pseudomonas aeruginosa exhibits the greatest diversity of anti-phage systems, as measured by Shannon entropy. Using network analysis to identify the common neighbors of anti-phage systems, we identified two core defense hotspot loci (cDHS1 and cDHS2). cDHS1 is up to 224 kb (median: 26 kb) with varied arrangements of more than 30 distinct immune systems across isolates, while cDHS2 has 24 distinct systems (median: 6 kb). Both cDHS regions are occupied in a majority of P. aeruginosa isolates. Most cDHS genes are of unknown function potentially representing new anti-phage systems, which we validated by identifying a novel anti-phage system (Shango) commonly encoded in cDHS1. Identifying core genes flanking immune islands could simplify immune system discovery and may represent popular landing spots for diverse MGEs carrying anti-phage systems.

The marine thermophilic archaeon Nanoarchaeum equitans possesses a monomeric primase encompassing the conserved domains of the small catalytic and the large regulatory subunits of archaeoeukaryotic heterodimeric primases in one protein chain. The recombinant protein primes on templates containing a triplet with a central thymidine, thus displaying a pronounced sequence specificity typically observed with bacterial type primases only. The N. equitans primase (NEQ395) is a highly active primase enzyme synthesizing short RNA primers. Termination occurs preferentially at about nine nucleotides, as determined by HPLC analysis and confirmed with mass spectrometry. Possibly, the compact monomeric primase NEQ395 represents the minimal archaeoeukaryotic primase and could serve as a functional and structural model of the heterodimeric archaeoeukaryotic primases, whose study is hindered by engagement in protein assemblies and rather low activity.

Facioscapulohumeral muscular dystrophy (FSHD) is one of the most prevalent neuromuscular disorders. The disease is linked to copy number reduction and/or epigenetic alterations of the D4Z4 macrosatellite on chromosome 4q35 and associated with aberrant gain of expression of the transcription factor DUX4, which triggers a pro-apoptotic transcriptional program leading to muscle wasting. As today, no cure or therapeutic option is available to FSHD patients. Given its centrality in FSHD, blocking DUX4 expression with small molecule drugs is an attractive option. We previously showed that the long non protein-coding RNA DBE-T is required for aberrant DUX4 expression in FSHD. Using affinity purification followed by proteomics, here we identified the chromatin remodeling protein WDR5 as a novel DBE-T interactor and a key player required for the biological activity of the lncRNA. We found that WDR5 is required for the expression of DUX4 and its targets in primary FSHD muscle cells. Moreover, targeting WDR5 rescues both cell viability and myogenic differentiation of FSHD patient cells. Notably, comparable results were obtained by pharmacological inhibition of WDR5. Importantly, WDR5 targeting was safe to healthy donor muscle cells. Our results support a pivotal role of WDR5 in the activation of DUX4 expression identifying a druggable target for an innovative therapeutic approach for FSHD.

Accumulation of DNA damage resulting from reactive oxygen species was proposed to cause neurological and degenerative disease in patients, deficient in nucleotide excision repair (NER) or its transcription-coupled subpathway (TC-NER). Here, we assessed the requirement of TC-NER for the repair of specific types of oxidatively generated DNA modifications. We incorporated synthetic 5',8-cyclo-2'-deoxypurine nucleotides (cyclo-dA, cyclo-dG) and thymine glycol (Tg) into an EGFP reporter gene to measure transcription-blocking potentials of these modifications in human cells. Using null mutants, we further identified the relevant DNA repair components by a host cell reactivation approach. The results indicated that NTHL1-initiated base excision repair is by far the most efficient pathway for Tg. Moreover, Tg was efficiently bypassed during transcription, which effectively rules out TC-NER as an alternative repair mechanism. In a sharp contrast, both cyclopurine lesions robustly blocked transcription and were repaired by NER, wherein the specific TC-NER components CSB/ERCC6 and CSA/ERCC8 were as essential as XPA. Instead, repair of classical NER substrates, cyclobutane pyrimidine dimer and N-(deoxyguanosin-8-yl)-2-acetylaminofluorene, occurred even when TC-NER was disrupted. The strict requirement of TC-NER highlights cyclo-dA and cyclo-dG as candidate damage types, accountable for cytotoxic and degenerative responses in individuals affected by genetic defects in this pathway.

In enteric bacteria, several small RNAs (sRNAs) including MicC employ endoribonuclease RNase E to stimulate target RNA decay. A current model proposes that interaction of the sRNA 5' monophosphate (5'P) with the N-terminal sensing pocket of RNase E allosterically activates cleavage of the base-paired target in the active site. In vivo evidence supporting this model is lacking. Here, we engineered a genetic tool allowing us to generate 5' monophosphorylated sRNAs of choice in a controllable manner in the cell. Four sRNAs were tested and none performed better in target destabilization when 5' monophosphorylated. MicC retains full activity even when RNase E is defective in 5'P sensing, whereas regulation is lost upon removal of its scaffolding domain. Interestingly, sRNAs MicC and RyhB that originate with a 5' triphosphate group are dramatically destabilized when 5' monophosphorylated, but stable when in 5' triphosphorylated form. In contrast, the processing-derived sRNAs CpxQ and SroC, which carry 5'P groups naturally, are highly stable. Thus, the 5' phosphorylation state determines stability of naturally triphosphorylated sRNAs, but plays no major role for target RNA destabilization in vivo. In contrast, the RNase E C-terminal half is crucial for MicC-mediated ompD decay, suggesting that interaction with Hfq is mandatory.

The DNA-glycosylase OGG1 oversees the detection and clearance of the 7,8-dihydro-8-oxoguanine (8-oxoG), which is the most frequent form of oxidized base in the genome. This lesion is deeply buried within the double-helix and its detection requires careful inspection of the bases by OGG1 via a mechanism that remains only partially understood. By analyzing OGG1 dynamics in the nucleus of living human cells, we demonstrate that the glycosylase constantly samples the DNA by rapidly alternating between diffusion within the nucleoplasm and short transits on the DNA. This sampling process, that we find to be tightly regulated by the conserved residue G245, is crucial for the rapid recruitment of OGG1 at oxidative lesions induced by laser micro-irradiation. Furthermore, we show that residues Y203, N149 and N150, while being all involved in early stages of 8-oxoG probing by OGG1 based on previous structural data, differentially regulate the sampling of the DNA and recruitment to oxidative lesions.

The Ccr4-Not complex is a conserved multi protein complex with diverse roles in the mRNA life cycle. Recently we determined that the Not1 and Not4 subunits of Ccr4-Not inversely regulate mRNA solubility and thereby impact dynamics of co-translation events. One mRNA whose solubility is limited by Not4 is MMF1 encoding a mitochondrial matrix protein. In this work we uncover a mechanism that limits MMF1 overexpression and depends upon its co-translational targeting to the mitochondria. We have named this mechanism Mito-ENCay. This mechanism relies on Not4 promoting ribosome pausing during MMF1 translation, and hence the co-translational docking of the MMF1 mRNA to mitochondria via the mitochondrial targeting sequence of the Mmf1 nascent chain, the Egd1 chaperone, the Om14 mitochondrial outer membrane protein and the co-translational import machinery. Besides co-translational Mitochondrial targeting, Mito-ENCay depends upon Egd1 ubiquitination by Not4, the Caf130 subunit of the Ccr4-Not complex, the mitochondrial outer membrane protein Cis1, autophagy and no-go-decay.

Transposon Tn4430 belongs to a widespread family of bacterial transposons, the Tn3 family, which plays a prevalent role in the dissemination of antibiotic resistance among pathogens. Despite recent data on the structural architecture of the transposition complex, the molecular mechanisms underlying the replicative transposition of these elements are still poorly understood. Here, we use force-distance curve-based atomic force microscopy to probe the binding of the TnpA transposase of Tn4430 to DNA molecules containing one or two transposon ends and to extract the thermodynamic and kinetic parameters of transposition complex assembly. Comparing wild-type TnpA with previously isolated deregulated TnpA mutants supports a stepwise pathway for transposition complex formation and activation during which TnpA first binds as a dimer to a single transposon end and then undergoes a structural transition that enables it to bind the second end cooperatively and to become activated for transposition catalysis, the latter step occurring at a much faster rate for the TnpA mutants. Our study thus provides an unprecedented approach to probe the dynamic of a complex DNA processing machinery at the single-particle level.

TbMex67 is the major mRNA export factor known to date in trypanosomes, forming part of the docking platform within the nuclear pore. To explore its role in co-transcriptional mRNA export, recently reported in Trypanosoma brucei, pulse labelling of nascent RNAs with 5-ethynyl uridine (5-EU) was performed with cells depleted of TbMex67 and complemented with a dominant-negative mutant (TbMex67-DN). RNA polymerase (Pol) II transcription was unaffected, but the procyclin loci, which encode mRNAs transcribed by Pol I from internal sites on chromosomes 6 and 10, showed increased levels of 5-EU incorporation. This was due to Pol I readthrough transcription, which proceeded beyond the procyclin and procyclin-associated genes up to the Pol II transcription start site on the opposite strand. Complementation by TbMex67-DN also increased Pol I-dependent formation of R-loops and γ-histone 2A foci. The DN mutant exhibited reduced nuclear localisation and binding to chromatin compared to wild-type TbMex67. Together with its interaction with chromatin remodelling factor TbRRM1 and Pol II, and transcription-dependent association of Pol II with nucleoporins, our findings support a role for TbMex67 in connecting transcription and export in T. brucei. In addition, TbMex67 stalls readthrough by Pol I in specific contexts, thereby limiting R-loop formation and replication stress.

The action of cis-regulatory elements with either activation or repression functions underpins the precise regulation of gene expression during normal development and cell differentiation. Gene activation by the combined activities of promoters and distal enhancers has been extensively studied in normal and pathological contexts. In sharp contrast, gene repression by cis-acting silencers, defined as genetic elements that negatively regulate gene transcription in a position-independent fashion, is less well understood. Here, we repurpose the STARR-seq approach as a novel high-throughput reporter strategy to quantitatively assess silencer activity in mammals. We assessed silencer activity from DNase hypersensitive I sites in a mouse T cell line. Identified silencers were associated with either repressive or active chromatin marks and enriched for binding motifs of known transcriptional repressors. CRISPR-mediated genomic deletions validated the repressive function of distinct silencers involved in the repression of non-T cell genes and genes regulated during T cell differentiation. Finally, we unravel an association of silencer activity with short tandem repeats, highlighting the role of repetitive elements in silencer activity. Our results provide a general strategy for genome-wide identification and characterization of silencer elements.

Recycling and de-novo deposition of histones during DNA replication is a critical challenge faced by eukaryotic cells and is coordinated by histone chaperones. Spermatogenesis is highly regulated sophisticated process necessitating not only histone modification but loading of testis specific histone variants. Here, we show that Germ Cell Nuclear Acidic protein (GCNA), a germ cell specific protein in adult mice, can bind histones and purified GCNA exhibits histone chaperone activity. GCNA associates with the DNA replication machinery and supports progression through S-phase in murine undifferentiated spermatogonia (USGs). Whilst GCNA is dispensable for embryonic germ cell development, it is required for the maintenance of the USG pool and for long-term production of sperm. Our work describes the role of a germ cell specific histone chaperone in USGs maintenance in mice. These findings provide a mechanistic basis for the male infertility observed in patients carrying GCNA mutations.

Long INterspersed Element 1 (LINE-1 or L1) acts as a major remodeling force in genome regulation and evolution. Accumulating evidence shows that virus infection impacts L1 expression, potentially impacting host antiviral response and diseases. The underlying regulation mechanism is unclear. Epstein-Barr virus (EBV), a double-stranded DNA virus linked to B-cell and epithelial malignancies, is known to have viral-host genome interaction, resulting in transcriptional rewiring in EBV-associated gastric cancer (EBVaGC). By analyzing publicly available datasets from the Gene Expression Omnibus (GEO), we found that EBVaGC has L1 transcriptional repression compared with EBV-negative gastric cancer (EBVnGC). More specifically, retrotransposition-associated young and full-length L1s (FL-L1s) were among the most repressed L1s. Epigenetic alterations, especially increased H3K9me3, were observed on FL-L1s. H3K9me3 deposition was potentially attributed to increased TASOR expression, a key component of the human silencing hub (HUSH) complex for H3K9 trimethylation. The 4C- and HiC-seq data indicated that the viral DNA interacted in the proximity of the TASOR enhancer, strengthening the loop formation between the TASOR enhancer and its promoter. These results indicated that EBV infection is associated with increased H3K9me3 deposition, leading to L1 repression. This study uncovers a regulation mechanism of L1 expression by chromatin topology remodeling associated with viral-host genome interaction in EBVaGC.

Although splicing occurs largely co-transcriptionally, the order by which introns are removed does not necessarily follow the order in which they are transcribed. Whereas several genomic features are known to influence whether or not an intron is spliced before its downstream neighbor, multiple questions related to adjacent introns' splicing order (AISO) remain unanswered. Here, we present Insplico, the first standalone software for quantifying AISO that works with both short and long read sequencing technologies. We first demonstrate its applicability and effectiveness using simulated reads and by recapitulating previously reported AISO patterns, which unveiled overlooked biases associated with long read sequencing. We next show that AISO around individual exons is remarkably constant across cell and tissue types and even upon major spliceosomal disruption, and it is evolutionarily conserved between human and mouse brains. We also establish a set of universal features associated with AISO patterns across various animal and plant species. Finally, we used Insplico to investigate AISO in the context of tissue-specific exons, particularly focusing on SRRM4-dependent microexons. We found that the majority of such microexons have non-canonical AISO, in which the downstream intron is spliced first, and we suggest two potential modes of SRRM4 regulation of microexons related to their AISO and various splicing-related features. Insplico is available on gitlab.com/aghr/insplico.

Because DNA double-strand breaks (DSBs) greatly threaten genomic integrity, effective DNA damage sensing and repair are essential for cellular survival in all organisms. However, DSB repair mainly occurs during interphase and is repressed during mitosis. Here, we show that, unlike mitotic cells, oocytes can repair DSBs during meiosis I through microtubule-dependent chromosomal recruitment of the CIP2A-MDC1-TOPBP1 complex from spindle poles. After DSB induction, we observed spindle shrinkage and stabilization, as well as BRCA1 and 53BP1 recruitment to chromosomes and subsequent DSB repair during meiosis I. Moreover, p-MDC1 and p-TOPBP1 were recruited from spindle poles to chromosomes in a CIP2A-dependent manner. This pole-to-chromosome relocation of the CIP2A-MDC1-TOPBP1 complex was impaired not only by depolymerizing microtubules but also by depleting CENP-A or HEC1, indicating that the kinetochore/centromere serves as a structural hub for microtubule-dependent transport of the CIP2A-MDC1-TOPBP1 complex. Mechanistically, DSB-induced CIP2A-MDC1-TOPBP1 relocation is regulated by PLK1 but not by ATM activity. Our data provide new insights into the critical crosstalk between chromosomes and spindle microtubules in response to DNA damage to maintain genomic stability during oocyte meiosis.

UV-damaged DNA-binding protein (UV-DDB) is a heterodimeric protein, consisting of DDB1 and DDB2 subunits, that works to recognize DNA lesions induced by UV damage during global genome nucleotide excision repair (GG-NER). Our laboratory previously discovered a non-canonical role for UV-DDB in the processing of 8-oxoG, by stimulating 8-oxoG glycosylase, OGG1, activity 3-fold, MUTYH activity 4-5-fold, and APE1 (apurinic/apyrimidinic endonuclease 1) activity 8-fold. 5-hydroxymethyl-deoxyuridine (5-hmdU) is an important oxidation product of thymidine which is removed by single-strand selective monofunctional DNA glycosylase (SMUG1). Biochemical experiments with purified proteins indicated that UV-DDB stimulates the excision activity of SMUG1 on several substrates by 4-5-fold. Electrophoretic mobility shift assays indicated that UV-DDB displaced SMUG1 from abasic site products. Single-molecule analysis revealed that UV-DDB decreases the half-life of SMUG1 on DNA by ∼8-fold. Immunofluorescence experiments demonstrated that cellular treatment with 5-hmdU (5 μM for 15 min), which is incorporated into DNA during replication, produces discrete foci of DDB2-mCherry, which co-localize with SMUG1-GFP. Proximity ligation assays supported a transient interaction between SMUG1 and DDB2 in cells. Poly(ADP)-ribose accumulated after 5-hmdU treatment, which was abrogated with SMUG1 and DDB2 knockdown. These data support a novel role for UV-DDB in the processing of the oxidized base, 5-hmdU.

The Paf1 complex (Paf1C) is a conserved transcription elongation factor that regulates transcription elongation efficiency, facilitates co-transcriptional histone modifications, and impacts molecular processes linked to RNA synthesis, such as polyA site selection. Coupling of the activities of Paf1C to transcription elongation requires its association with RNA polymerase II (Pol II). Mutational studies in yeast identified Paf1C subunits Cdc73 and Rtf1 as important mediators of Paf1C association with Pol II on active genes. While the interaction between Rtf1 and the general elongation factor Spt5 is relatively well-understood, the interactions involving Cdc73 have not been fully elucidated. Using a site-specific protein cross-linking strategy in yeast cells, we identified direct interactions between Cdc73 and two components of the Pol II elongation complex, the elongation factor Spt6 and the largest subunit of Pol II. Both of these interactions require the tandem SH2 domain of Spt6. We also show that Cdc73 and Spt6 can interact in vitro and that rapid depletion of Spt6 dissociates Paf1 from chromatin, altering patterns of Paf1C-dependent histone modifications genome-wide. These results reveal interactions between Cdc73 and the Pol II elongation complex and identify Spt6 as a key factor contributing to the occupancy of Paf1C at active genes in Saccharomyces cerevisiae.

Endogenous retroviruses (ERVs) have been proposed as a driving force for the evolution of the mammalian placenta, however, the contribution of ERVs to placental development and the underlying regulatory mechanism remain largely elusive. A key process of placental development is the formation of multinucleated syncytiotrophoblasts (STBs) in direct contact with maternal blood, through which constitutes the maternal-fetal interface critical for nutrient allocation, hormone production and immunological modulation during pregnancy. We delineate that ERVs profoundly rewire the transcriptional program of trophoblast syncytialization. Here, we first determined the dynamic landscape of bivalent ERV-derived enhancers with dual occupancy of H3K27ac and H3K9me3 in human trophoblast stem cells (hTSCs). We further demonstrated that enhancers overlapping several ERV families tend to exhibit increased H3K27ac and reduced H3K9me3 occupancy in STBs relative to hTSCs. Particularly, bivalent enhancers derived from the Simiiformes-specific MER50 transposons were linked to a cluster of genes important for STB formation. Importantly, deletions of MER50 elements adjacent to several STB genes, including MFSD2A and TNFAIP2, significantly attenuated their expression concomitant to compromised syncytium formation. Together, we propose that ERV-derived enhancers, MER50 specifically, fine-tune the transcriptional networks accounting for human trophoblast syncytialization, which sheds light on a novel ERV-mediated regulatory mechanism underlying placental development.

Intrinsically disordered regions (IDRs) direct transcription factors (TFs) towards selected genomic occurrences of their binding motif, as exemplified by budding yeast's Msn2. However, the sequence basis of IDR-directed TF binding selectivity remains unknown. To reveal this sequence grammar, we analyze the genomic localizations of >100 designed IDR mutants, each carrying up to 122 mutations within this 567-AA region. Our data points at multivalent interactions, carried by hydrophobic-mostly aliphatic-residues dispersed within a disordered environment and independent of linear sequence motifs, as the key determinants of Msn2 genomic localization. The implications of our results for the mechanistic basis of IDR-based TF binding preferences are discussed.

Normal erythropoiesis requires the precise regulation of gene expression patterns, and transcription cofactors play a vital role in this process. Deregulation of cofactors has emerged as a key mechanism contributing to erythroid disorders. Through gene expression profiling, we found HES6 as an abundant cofactor expressed at gene level during human erythropoiesis. HES6 physically interacted with GATA1 and influenced the interaction of GATA1 with FOG1. Knockdown of HES6 impaired human erythropoiesis by decreasing GATA1 expression. Chromatin immunoprecipitation and RNA sequencing revealed a rich set of HES6- and GATA1-co-regulated genes involved in erythroid-related pathways. We also discovered a positive feedback loop composed of HES6, GATA1 and STAT1 in the regulation of erythropoiesis. Notably, erythropoietin (EPO) stimulation led to up-regulation of these loop components. Increased expression levels of loop components were observed in CD34+ cells of polycythemia vera patients. Interference by either HES6 knockdown or inhibition of STAT1 activity suppressed proliferation of erythroid cells with the JAK2V617F mutation. We further explored the impact of HES6 on polycythemia vera phenotypes in mice. The identification of the HES6-GATA1 regulatory loop and its regulation by EPO provides novel insights into human erythropoiesis regulated by EPO/EPOR and a potential therapeutic target for the management of polycythemia vera.

Besides entrapping sister chromatids, cohesin drives other high-order chromosomal structural dynamics like looping, compartmentalization and condensation. ESCO2 acetylates a subset of cohesin so that cohesion must be established and only be established between nascent sister chromatids. How this process is precisely achieved remains unknown. Here, we report that GSK3 family kinases provide higher hierarchical control through an ESCO2 regulator, CRL4MMS22L. GSK3s phosphorylate Thr105 in MMS22L, resulting in homo-dimerization of CRL4MMS22L and ESCO2 during S phase as evidenced by single-molecule spectroscopy and several biochemical approaches. A single phospho-mimicking mutation on MMS22L (T105D) is sufficient to mediate their dimerization and rescue the cohesion defects caused by GSK3 or MMS22L depletion, whereas non-phosphorylable T105A exerts dominant-negative effects even in wildtype cells. Through cell fractionation and time-course measurements, we show that GSK3s facilitate the timely chromatin association of MMS22L and ESCO2 and subsequently SMC3 acetylation. The necessity of ESCO2 dimerization implicates symmetric control of cohesion establishment in eukaryotes.

Inter-chromosomal interactions play a crucial role in genome organization, yet the organizational principles remain elusive. Here, we introduce a novel computational method to systematically characterize inter-chromosomal interactions using in situ Hi-C results from various cell types. Our method successfully identifies two apparently hub-like inter-chromosomal contacts associated with nuclear speckles and nucleoli, respectively. Interestingly, we discover that nuclear speckle-associated inter-chromosomal interactions are highly cell-type invariant with a marked enrichment of cell-type common super-enhancers (CSEs). Validation using DNA Oligopaint fluorescence in situ hybridization (FISH) shows a strong but probabilistic interaction behavior between nuclear speckles and CSE-harboring genomic regions. Strikingly, we find that the likelihood of speckle-CSE associations can accurately predict two experimentally measured inter-chromosomal contacts from Hi-C and Oligopaint DNA FISH. Our probabilistic establishment model well describes the hub-like structure observed at the population level as a cumulative effect of summing individual stochastic chromatin-speckle interactions. Lastly, we observe that CSEs are highly co-occupied by MAZ binding and MAZ depletion leads to significant disorganization of speckle-associated inter-chromosomal contacts. Taken together, our results propose a simple organizational principle of inter-chromosomal interactions mediated by MAZ-occupied CSEs.

Structural analysis of RNA is an important and versatile tool to investigate the function of this type of molecules in the cell as well as in vitro. Several robust and reliable procedures are available, relying on chemical modification inducing RT stops or nucleotide misincorporations during reverse transcription. Others are based on cleavage reactions and RT stop signals. However, these methods address only one side of the RT stop or misincorporation position. Here, we describe Led-Seq, a new approach based on lead-induced cleavage of unpaired RNA positions, where both resulting cleavage products are investigated. The RNA fragments carrying 2', 3'-cyclic phosphate or 5'-OH ends are selectively ligated to oligonucleotide adapters by specific RNA ligases. In a deep sequencing analysis, the cleavage sites are identified as ligation positions, avoiding possible false positive signals based on premature RT stops. With a benchmark set of transcripts in Escherichia coli, we show that Led-Seq is an improved and reliable approach based on metal ion-induced phosphodiester hydrolysis to investigate RNA structures in vivo.

Despite the need in various applications, accurate quantification of nucleic acids still remains a challenge. The widely-used qPCR has reduced accuracy at ultralow template concentration and is susceptible to nonspecific amplifications. The more recently developed dPCR is costly and cannot handle high-concentration samples. We combine the strengths of qPCR and dPCR by performing PCR in silicon-based microfluidic chips and demonstrate high quantification accuracy in a large concentration range. Importantly, at low template concentration, we observe on-site PCR (osPCR), where only certain sites of the channel show amplification. The sites have almost identical ct values, showing osPCR is a quasi-single molecule phenomenon. Using osPCR, we can measure both the ct values and the absolute concentration of templates in the same reaction. Additionally, osPCR enables identification of each template molecule, allowing removal of nonspecific amplification during quantification and greatly improving quantification accuracy. We develop sectioning algorithm that improves the signal amplitude and demonstrate improved detection of COVID in patient samples.

In bacteria, release of newly synthesized proteins from ribosomes during translation termination is catalyzed by class-I release factors (RFs) RF1 or RF2, reading UAA and UAG or UAA and UGA codons, respectively. Class-I RFs are recycled from the post-termination ribosome by a class-II RF, the GTPase RF3, which accelerates ribosome intersubunit rotation and class-I RF dissociation. How conformational states of the ribosome are coupled to the binding and dissociation of the RFs remains unclear and the importance of ribosome-catalyzed guanine nucleotide exchange on RF3 for RF3 recycling in vivo has been disputed. Here, we profile these molecular events using a single-molecule fluorescence assay to clarify the timings of RF3 binding and ribosome intersubunit rotation that trigger class-I RF dissociation, GTP hydrolysis, and RF3 dissociation. These findings in conjunction with quantitative modeling of intracellular termination flows reveal rapid ribosome-dependent guanine nucleotide exchange to be crucial for RF3 action in vivo.

The single-stranded DNA (ssDNA) binding protein complex RPA plays a critical role in promoting DNA replication and multiple DNA repair pathways. However, how RPA is regulated to achieve its functions precisely in these processes remains elusive. Here, we found that proper acetylation and deacetylation of RPA are required to regulate RPA function in promoting high-fidelity DNA replication and repair. We show that yeast RPA is acetylated on multiple conserved lysines by the acetyltransferase NuA4 upon DNA damage. Mimicking constitutive RPA acetylation or blocking its acetylation causes spontaneous mutations with the signature of micro-homology-mediated large deletions or insertions. In parallel, improper RPA acetylation/deacetylation impairs DNA double-strand break (DSB) repair by the accurate gene conversion or break-induced replication while increasing the error-prone repair by single-strand annealing or alternative end joining. Mechanistically, we show that proper acetylation and deacetylation of RPA ensure its normal nuclear localization and ssDNA binding ability. Importantly, mutation of the equivalent residues in human RPA1 also impairs RPA binding on ssDNA, leading to attenuated RAD51 loading and homologous recombination repair. Thus, timely RPA acetylation and deacetylation likely represent a conserved mechanism promoting high-fidelity replication and repair while discriminating the error-prone repair mechanisms in eukaryotes.

DNA translocases, such as RNA polymerases, inevitably collide with nucleosomes on eukaryotic chromatin. Upon these collisions, histone chaperones are suggested to facilitate nucleosome disassembly and re-assembly. In this study, by performing in vitro transcription assays and molecular simulations, we found that partial unwrapping of a nucleosome by an RNA polymerase dramatically facilitates an H2A/H2B dimer dismantling from the nucleosome by Nucleosome Assembly Protein 1 (Nap1). Furthermore, the results uncovered molecular mechanisms of Nap1 functions in which the highly acidic C-terminal flexible tails of Nap1 contribute to the H2A/H2B binding by associating with the binding interface buried and not accessible to Nap1 globular domains, supporting the penetrating fuzzy binding mechanism seemingly shared across various histone chaperones. These findings have broad implications for the mechanisms by which histone chaperones process nucleosomes upon collisions with translocases in transcription, histone recycling and nucleosomal DNA repair.

In this study, we specifically visualized DNA molecules at their AT base pairs after in vitro phage ejection. Our AT-specific visualization revealed that either end of the DNA molecule could be ejected first with a nearly 50% probability. This observation challenges the generally accepted theory of Last In First Out (LIFO), which states that the end of the phage λ DNA that enters the capsid last during phage packaging is the first to be ejected, and that both ends of the DNA are unable to move within the extremely condensed phage capsid. To support our observations, we conducted computer simulations that revealed that both ends of the DNA molecule are randomized, resulting in the observed near 50% probability. Additionally, we found that the length of the ejected DNA by LIFO was consistently longer than that by First In First Out (FIFO) during in vitro phage ejection. Our simulations attributed this difference in length to the stiffness difference of the remaining DNA within the phage capsid. In conclusion, this study demonstrates that a DNA molecule within an extremely dense phage capsid exhibits a degree of mobility, allowing it to switch ends during ejection.

Quantifying the nucleotide preferences of DNA binding proteins is essential to understanding how transcription factors (TFs) interact with their targets in the genome. High-throughput in vitro binding assays have been used to identify the inherent DNA binding preferences of TFs in a controlled environment isolated from confounding factors such as genome accessibility, DNA methylation, and TF binding cooperativity. Unfortunately, many of the most common approaches for measuring binding preferences are not sensitive enough for the study of moderate-to-low affinity binding sites, and are unable to detect small-scale differences between closely related homologs. The Forkhead box (FOX) family of TFs is known to play a crucial role in regulating a variety of key processes from proliferation and development to tumor suppression and aging. By using the high-sequencing depth SELEX-seq approach to study all four FOX homologs in Saccharomyces cerevisiae, we have been able to precisely quantify the contribution and importance of nucleotide positions all along an extended binding site. Essential to this process was the alignment of our SELEX-seq reads to a set of candidate core sequences determined using a recently developed tool for the alignment of enriched k-mers and a newly developed approach for the reprioritization of candidate cores.

Dynamic protein gradients are exploited for the spatial organization and segregation of replicated chromosomes. However, mechanisms of protein gradient formation and how that spatially organizes chromosomes remain poorly understood. Here, we have determined the kinetic principles of subcellular localizations of ParA2 ATPase, an essential spatial regulator of chromosome 2 segregation in the multichromosome bacterium, Vibrio cholerae. We found that ParA2 gradients self-organize in V. cholerae cells into dynamic pole-to-pole oscillations. We examined the ParA2 ATPase cycle and ParA2 interactions with ParB2 and DNA. In vitro, ParA2-ATP dimers undergo a rate-limiting conformational switch, catalysed by DNA to achieve DNA-binding competence. This active ParA2 state loads onto DNA cooperatively as higher order oligomers. Our results indicate that the midcell localization of ParB2-parS2 complexes stimulate ATP hydrolysis and ParA2 release from the nucleoid, generating an asymmetric ParA2 gradient with maximal concentration toward the poles. This rapid dissociation coupled with slow nucleotide exchange and conformational switch provides for a temporal lag that allows the redistribution of ParA2 to the opposite pole for nucleoid reattachment. Based on our data, we propose a 'Tug-of-war' model that uses dynamic oscillations of ParA2 to spatially regulate symmetric segregation and positioning of bacterial chromosomes.

Classic promoter mutagenesis strategies can be used to study how proximal promoter regions regulate the expression of particular genes of interest. This is a laborious process, in which the smallest sub-region of the promoter still capable of recapitulating expression in an ectopic setting is first identified, followed by targeted mutation of putative transcription factor binding sites. Massively parallel reporter assays such as survey of regulatory elements (SuRE) provide an alternative way to study millions of promoter fragments in parallel. Here we show how a generalized linear model (GLM) can be used to transform genome-scale SuRE data into a high-resolution genomic track that quantifies the contribution of local sequence to promoter activity. This coefficient track helps identify regulatory elements and can be used to predict promoter activity of any sub-region in the genome. It thus allows in silico dissection of any promoter in the human genome to be performed. We developed a web application, available at cissector.nki.nl, that lets researchers easily perform this analysis as a starting point for their research into any promoter of interest.

Aptamers are ligand-binding RNA or DNA molecules and have been widely examined as biosensors, diagnostic tools, and therapeutic agents. The application of aptamers as biosensors commonly requires an expression platform to produce a signal to report the aptamer-ligand binding event. Traditionally, aptamer selection and expression platform integration are two independent steps and the aptamer selection requires the immobilization of either the aptamer or the ligand. These drawbacks can be easily overcome through the selection of allosteric DNAzymes (aptazymes). Herein, we used the technique of Expression-SELEX developed in our laboratory to select for aptazymes that can be specifically activated by low concentrations of l-phenylalanine. We chose a previous DNA-cleaving DNAzyme known as II-R1 as the expression platform for its low cleavage rate and used stringent selection conditions to drive the selection of high-performance aptazyme candidates. Three aptazymes were chosen for detailed characterization and these DNAzymes were found to exhibit a dissociation constant for l-phenylalanine as low as 4.8 μM, a catalytic rate constant improvement as high as 20 000-fold in the presence of l-phenylalanine, and the ability to discriminate against closely related l-phenylalanine analogs including d-phenylalanine. This work has established the Expression-SELEX as an effective SELEX method to enrich high-quality ligand-responsive aptazymes.

Understanding the effects of genetic variation in gene regulatory elements is crucial to interpreting genome function. This is particularly pertinent for the hundreds of thousands of disease-associated variants identified by GWAS, which frequently sit within gene regulatory elements but whose functional effects are often unknown. Current methods are limited in their scalability and ability to assay regulatory variants in their endogenous context, independently of other tightly linked variants. Here, we present a new medium-throughput screening system: genome engineering based interrogation of enhancers assay for transposase accessible chromatin (GenIE-ATAC), that measures the effect of individual variants on chromatin accessibility in their endogenous genomic and chromatin context. We employ this assay to screen for the effects of regulatory variants in human induced pluripotent stem cells, validating a subset of causal variants, and extend our software package (rgenie) to analyse these new data. We demonstrate that this methodology can be used to understand the impact of defined deletions and point mutations within transcription factor binding sites. We thus establish GenIE-ATAC as a method to screen for the effect of gene regulatory element variation, allowing identification and prioritisation of causal variants from GWAS for functional follow-up and understanding the mechanisms of regulatory element function.

Single-stranded DNA (ssDNA) gapped regions are common intermediates in DNA transactions. Using a new non-denaturing bisulfite treatment combined with ChIP-seq, abbreviated 'ssGap-seq', we explore RecA and SSB binding to ssDNA on a genomic scale in E. coli in a wide range of genetic backgrounds. Some results are expected. During log phase growth, RecA and SSB assembly profiles coincide globally, concentrated on the lagging strand and enhanced after UV irradiation. Unexpected results also abound. Near the terminus, RecA binding is favored over SSB, binding patterns change in the absence of RecG, and the absence of XerD results in massive RecA assembly. RecA may substitute for the absence of XerCD to resolve chromosome dimers. A RecA loading pathway may exist that is independent of RecBCD and RecFOR. Two prominent and focused peaks of RecA binding revealed a pair of 222 bp and GC-rich repeats, equidistant from dif and flanking the Ter domain. The repeats, here named RRS for replication risk sequence, trigger a genomically programmed generation of post-replication gaps that may play a special role in relieving topological stress during replication termination and chromosome segregation. As demonstrated here, ssGap-seq provides a new window on previously inaccessible aspects of ssDNA metabolism.

Many transcription factors (TFs) localize in nuclear clusters of locally increased concentrations, but how TF clustering is regulated and how it influences gene expression is not well understood. Here, we use quantitative microscopy in living cells to study the regulation and function of clustering of the budding yeast TF Gal4 in its endogenous context. Our results show that Gal4 forms clusters that overlap with the GAL loci. Cluster number, density and size are regulated in different growth conditions by the Gal4-inhibitor Gal80 and Gal4 concentration. Gal4 truncation mutants reveal that Gal4 clustering is facilitated by, but does not completely depend on DNA binding and intrinsically disordered regions. Moreover, we discover that clustering acts as a double-edged sword: self-interactions aid TF recruitment to target genes, but recruited Gal4 molecules that are not DNA-bound do not contribute to, and may even inhibit, transcription activation. We propose that cells need to balance the different effects of TF clustering on target search and transcription activation to facilitate proper gene expression.

Accumulating evidence suggests that posttranscriptional control of gene expression, including RNA splicing, transport, modification, translation and degradation, primarily relies on RNA binding proteins (RBPs). However, the functions of many RBPs remain understudied. Here, we characterized the function of a novel RBP, Proline-Rich Coiled-coil 2B (PRRC2B). Through photoactivatable ribonucleoside-enhanced crosslinking and immunoprecipitation and sequencing (PAR-CLIP-seq), we identified transcriptome-wide CU- or GA-rich PRRC2B binding sites near the translation initiation codon on a specific cohort of mRNAs in HEK293T cells. These mRNAs, including oncogenes and cell cycle regulators such as CCND2 (cyclin D2), exhibited decreased translation upon PRRC2B knockdown as revealed by polysome-associated RNA-seq, resulting in reduced G1/S phase transition and cell proliferation. Antisense oligonucleotides blocking PRRC2B interactions with CCND2 mRNA decreased its translation, thus inhibiting G1/S transition and cell proliferation. Mechanistically, PRRC2B interactome analysis revealed RNA-independent interactions with eukaryotic translation initiation factors 3 (eIF3) and 4G2 (eIF4G2). The interaction with translation initiation factors is essential for PRRC2B function since the eIF3/eIF4G2-interacting defective mutant, unlike wild-type PRRC2B, failed to rescue the translation deficiency or cell proliferation inhibition caused by PRRC2B knockdown. Altogether, our findings reveal that PRRC2B is essential for efficiently translating specific proteins required for cell cycle progression and cell proliferation.

The compaction of mitochondrial DNA (mtDNA) is regulated by architectural HMG-box proteins whose limited cross-species similarity suggests diverse underlying mechanisms. Viability of Candida albicans, a human antibiotic-resistant mucosal pathogen, is compromised by altering mtDNA regulators. Among them, there is the mtDNA maintenance factor Gcf1p, which differs in sequence and structure from its human and Saccharomyces cerevisiae counterparts, TFAM and Abf2p. Our crystallographic, biophysical, biochemical and computational analysis showed that Gcf1p forms dynamic protein/DNA multimers by a combined action of an N-terminal unstructured tail and a long helix. Furthermore, an HMG-box domain canonically binds the minor groove and dramatically bends the DNA while, unprecedentedly, a second HMG-box binds the major groove without imposing distortions. This architectural protein thus uses its multiple domains to bridge co-aligned DNA segments without altering the DNA topology, revealing a new mechanism of mtDNA condensation.

ANKLE1 is a nuclease that provides a final opportunity to process unresolved junctions in DNA that would otherwise create chromosomal linkages blocking cell division. It is a GIY-YIG nuclease. We have expressed an active domain of human ANKLE1 containing the GIY-YIG nuclease domain in bacteria, that is monomeric in solution and when bound to a DNA Y-junction, and unilaterally cleaves a cruciform junction. Using an AlphaFold model of the enzyme we identify the key active residues, and show that mutation of each leads to impairment of activity. There are two components in the catalytic mechanism. Cleavage rate is pH dependent, corresponding to a pKa of 6.9, suggesting an involvement of the conserved histidine in proton transfer. The reaction rate depends on the nature of the divalent cation, likely bound by glutamate and asparagine side chains, and is log-linear with the metal ion pKa. We propose that the reaction is subject to general acid-base catalysis, using a combination of tyrosine and histidine acting as general base and water directly coordinated to the metal ion as general acid. The reaction is temperature dependent; activation energy Ea = 37 kcal mol-1, suggesting that cleavage is coupled to opening of DNA in the transition state.

Streptococcus pyogenes Cas9 (SpCas9) nuclease exhibits considerable position-dependent sequence preferences. The reason behind these preferences is not well understood and is difficult to rationalise, since the protein establishes interactions with the target-spacer duplex in a sequence-independent manner. We revealed here that intramolecular interactions within the single guide RNA (sgRNA), between the spacer and the scaffold, cause most of these preferences. By using in cellulo and in vitro SpCas9 activity assays with systematically designed spacer and scaffold sequences and by analysing activity data from a large SpCas9 sequence library, we show that some long (>8 nucleotides) spacer motifs, that are complementary to the RAR unit of the scaffold, interfere with sgRNA loading, and that some motifs of more than 4 nucleotides, that are complementary to the SL1 unit, inhibit DNA binding and cleavage. Furthermore, we show that intramolecular interactions are present in the majority of the inactive sgRNA sequences of the library, suggesting that they are the most important intrinsic determinants of the activity of the SpCas9 ribonucleoprotein complex. We also found that in pegRNAs, sequences at the 3' extension of the sgRNA that are complementary to the SL2 unit are also inhibitory to prime editing, but not to the nuclease activity of SpCas9.

Coordination of mitochondrial and nuclear processes is key to the cellular health; however, very little is known about the molecular mechanisms regulating nuclear-mitochondrial crosstalk. Here, we report a novel molecular mechanism controlling the shuttling of CREB (cAMP response element-binding protein) protein complex between mitochondria and nucleoplasm. We show that a previously unknown protein, herein termed as Jig, functions as a tissue-specific and developmental timing-specific coregulator in the CREB pathway. Our results demonstrate that Jig shuttles between mitochondria and nucleoplasm, interacts with CrebA protein and controls its delivery to the nucleus, thus triggering CREB-dependent transcription in nuclear chromatin and mitochondria. Ablating the expression of Jig prevents CrebA from localizing to the nucleoplasm, affecting mitochondrial functioning and morphology and leads to Drosophila developmental arrest at the early third instar larval stage. Together, these results implicate Jig as an essential mediator of nuclear and mitochondrial processes. We also found that Jig belongs to a family of nine similar proteins, each of which has its own tissue- and time-specific expression profile. Thus, our results are the first to describe the molecular mechanism regulating nuclear and mitochondrial processes in a tissue- and time-specific manner.

In response to oxidative stress cells reprogram gene expression to enhance levels of antioxidant enzymes and promote survival. In Saccharomyces cerevisiae the polysome-interacting La-related proteins (LARPs) Slf1 and Sro9 aid adaptation of protein synthesis during stress by undetermined means. To gain insight in their mechanisms of action in stress responses, we determined LARP mRNA binding positions in stressed and unstressed cells. Both proteins bind within coding regions of stress-regulated antioxidant enzyme and other highly translated mRNAs in both optimal and stressed conditions. LARP interaction sites are framed and enriched with ribosome footprints suggesting ribosome-LARP-mRNA complexes are identified. Although stress-induced translation of antioxidant enzyme mRNAs is attenuated in slf1Δ, these mRNAs remain on polysomes. Focusing further on Slf1, we find it binds to both monosomes and disomes following RNase treatment. slf1Δ reduces disome enrichment during stress and alters programmed ribosome frameshifting rates. We propose that Slf1 is a ribosome-associated translational modulator that stabilises stalled/collided ribosomes, prevents ribosome frameshifting and so promotes translation of a set of highly-translated mRNAs that together facilitate cell survival and adaptation to stress.

Mitochondrial U-indel RNA editing in kinetoplastid protozoa is directed by trans-acting gRNAs and mediated by a holoenzyme with associated factors. Here, we examine the function of the holoenzyme-associated KREH1 RNA helicase in U-indel editing. We show that KREH1 knockout (KO) impairs editing of a small subset of mRNAs. Overexpression of helicase-dead mutants results in expanded impairment of editing across multiple transcripts, suggesting the existence of enzymes that can compensate for KREH1 in KO cells. In depth analysis of editing defects using quantitative RT-PCR and high-throughput sequencing reveals compromised editing initiation and progression in both KREH1-KO and mutant-expressing cells. In addition, these cells exhibit a distinct defect in the earliest stages of editing in which the initiator gRNA is bypassed, and a small number of editing events takes place just outside this region. Wild type KREH1 and a helicase-dead KREH1 mutant interact similarly with RNA and holoenzyme, and overexpression of both similarly disorders holoenzyme homeostasis. Thus, our data support a model in which KREH1 RNA helicase activity facilitates remodeling of initiator gRNA-mRNA duplexes to permit accurate utilization of initiating gRNAs on multiple transcripts.

Saccharomyces cerevisiae DNA polymerase IV (Pol4) like its homolog, human DNA polymerase lambda (Polλ), is involved in Non-Homologous End-Joining and Microhomology-Mediated Repair. Using genetic analysis, we identified an additional role of Pol4 also in homology-directed DNA repair, specifically in Rad52-dependent/Rad51-independent direct-repeat recombination. Our results reveal that the requirement for Pol4 in repeat recombination was suppressed by the absence of Rad51, suggesting that Pol4 counteracts the Rad51 inhibition of Rad52-mediated repeat recombination events. Using purified proteins and model substrates, we reconstituted in vitro reactions emulating DNA synthesis during direct-repeat recombination and show that Rad51 directly inhibits Polδ DNA synthesis. Interestingly, although Pol4 was not capable of performing extensive DNA synthesis by itself, it aided Polδ in overcoming the DNA synthesis inhibition by Rad51. In addition, Pol4 dependency and stimulation of Polδ DNA synthesis in the presence of Rad51 occurred in reactions containing Rad52 and RPA where DNA strand-annealing was necessary. Mechanistically, yeast Pol4 displaces Rad51 from ssDNA independent of DNA synthesis. Together our in vitro and in vivo data suggest that Rad51 suppresses Rad52-dependent/Rad51-independent direct-repeat recombination by binding to the primer-template and that Rad51 removal by Pol4 is critical for strand-annealing dependent DNA synthesis.

Methods for cell clustering and gene expression from single-cell RNA sequencing (scRNA-seq) data are essential for biological interpretation of cell processes. Here, we present TRIAGE-Cluster which uses genome-wide epigenetic data from diverse bio-samples to identify genes demarcating cell diversity in scRNA-seq data. By integrating patterns of repressive chromatin deposited across diverse cell types with weighted density estimation, TRIAGE-Cluster determines cell type clusters in a 2D UMAP space. We then present TRIAGE-ParseR, a machine learning method which evaluates gene expression rank lists to define gene groups governing the identity and function of cell types. We demonstrate the utility of this two-step approach using atlases of in vivo and in vitro cell diversification and organogenesis. We also provide a web accessible dashboard for analysis and download of data and software. Collectively, genome-wide epigenetic repression provides a versatile strategy to define cell diversity and study gene regulation of scRNA-seq data.

Photochemical dimerization of adjacent pyrimidines is fundamental to the creation of mutagenic hotspots caused by ultraviolet light. Distribution of the resulting lesions (cyclobutane pyrimidine dimers, CPDs) is already known to be highly variable in cells, and in vitro models have implicated DNA conformation as a major basis for this observation. Past efforts have primarily focused on mechanisms that influence CPD formation and have rarely considered contributions of CPD reversion. However, reversion is competitive under the standard conditions of 254 nm irradiation as illustrated in this report based on the dynamic response of CPDs to changes in DNA conformation. A periodic profile of CPDs was recreated in DNA held in a bent conformation by λ repressor. After linearization of this DNA, the CPD profile relaxed to its characteristic uniform distribution over a similar time of irradiation to that required to generate the initial profile. Similarly, when a T tract was released from a bent conformation, its CPD profile converted under further irradiation to that consistent with a linear T tract. This interconversion of CPDs indicates that both its formation and reversion exert control on CPD populations long before photo-steady-state conditions are achieved and suggests that the dominant sites of CPDs will evolve as DNA conformation changes in response to natural cellular processes.

Universal Minicircle Sequence binding proteins (UMSBPs) are CCHC-type zinc-finger proteins that bind the single-stranded G-rich UMS sequence, conserved at the replication origins of minicircles in the kinetoplast DNA, the mitochondrial genome of kinetoplastids. Trypanosoma brucei UMSBP2 has been recently shown to colocalize with telomeres and to play an essential role in chromosome end protection. Here we report that TbUMSBP2 decondenses in vitro DNA molecules, which were condensed by core histones H2B, H4 or linker histone H1. DNA decondensation is mediated via protein-protein interactions between TbUMSBP2 and these histones, independently of its previously described DNA binding activity. Silencing of the TbUMSBP2 gene resulted in a significant decrease in the disassembly of nucleosomes in T. brucei chromatin, a phenotype that could be reverted, by supplementing the knockdown cells with TbUMSBP2. Transcriptome analysis revealed that silencing of TbUMSBP2 affects the expression of multiple genes in T. brucei, with a most significant effect on the upregulation of the subtelomeric variant surface glycoproteins (VSG) genes, which mediate the antigenic variation in African trypanosomes. These observations suggest that UMSBP2 is a chromatin remodeling protein that functions in the regulation of gene expression and plays a role in the control of antigenic variation in T. brucei.

DNA double-strand break (DSB) repair via homologous recombination is initiated by end resection. The extent of DNA end resection determines the choice of the DSB repair pathway. Nucleases for end resection have been extensively studied. However, it is still unclear how the potential DNA structures generated by the initial short resection by MRE11-RAD50-NBS1 are recognized and recruit proteins, such as EXO1, to DSB sites to facilitate long-range resection. We found that the MSH2-MSH3 mismatch repair complex is recruited to DSB sites through interaction with the chromatin remodeling protein SMARCAD1. MSH2-MSH3 facilitates the recruitment of EXO1 for long-range resection and enhances its enzymatic activity. MSH2-MSH3 also inhibits access of POLθ, which promotes polymerase theta-mediated end-joining (TMEJ). Collectively, we present a direct role of MSH2-MSH3 in the initial stages of DSB repair by promoting end resection and influencing the DSB repair pathway by favoring homologous recombination over TMEJ.

Cellular totipotency is critical for whole-organism generation, yet how totipotency is established remains poorly illustrated. Abundant transposable elements (TEs) are activated in totipotent cells, which is critical for embryonic totipotency. Here, we show that the histone chaperone RBBP4, but not its homolog RBBP7, is indispensable for maintaining the identity of mouse embryonic stem cells (mESCs). Auxin-induced degradation of RBBP4, but not RBBP7, reprograms mESCs to the totipotent 2C-like cells. Also, loss of RBBP4 enhances transition from mESCs to trophoblast cells. Mechanistically, RBBP4 binds to the endogenous retroviruses (ERVs) and functions as an upstream regulator by recruiting G9a to deposit H3K9me2 on ERVL elements, and recruiting KAP1 to deposit H3K9me3 on ERV1/ERVK elements, respectively. Moreover, RBBP4 facilitates the maintenance of nucleosome occupancy at the ERVK and ERVL sites within heterochromatin regions through the chromatin remodeler CHD4. RBBP4 depletion leads to the loss of the heterochromatin marks and activation of TEs and 2C genes. Together, our findings illustrate that RBBP4 is required for heterochromatin assembly and is a critical barrier for inducing cell fate transition from pluripotency to totipotency.

G4C2 and G2C4 repeat expansions in chromosome 9 open reading frame 72 (C9orf72) are the most common cause of genetically defined amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD), or c9ALS/FTD. The gene is bidirectionally transcribed, producing G4C2 repeats [r(G4C2)exp] and G2C4 repeats [r(G2C4)exp]. The c9ALS/FTD repeat expansions are highly structured, and structural studies showed that r(G4C2)exp predominantly folds into a hairpin with a periodic array of 1 × 1 G/G internal loops and a G-quadruplex. A small molecule probe revealed that r(G4C2)exp also adopts a hairpin structure with 2 × 2 GG/GG internal loops. We studied the conformational dynamics adopted by 2 × 2 GG/GG loops using temperature replica exchange molecular dynamics (T-REMD) and further characterized the structure and underlying dynamics using traditional 2D NMR techniques. These studies showed that the loop's closing base pairs influence both structure and dynamics, particularly the configuration adopted around the glycosidic bond. Interestingly, r(G2C4) repeats, which fold into an array of 2 × 2 CC/CC internal loops, are not as dynamic. Collectively, these studies emphasize the unique sensitivity of r(G4C2)exp to small changes in stacking interactions, which is not observed in r(G2C4)exp, providing important considerations for further principles in structure-based drug design.

Transcription-associated cyclin-dependent kinases (CDKs) regulate the transcription cycle through sequential phosphorylation of RNA polymerase II (RNAPII). Herein, we report that dual inhibition of the highly homologous CDK12 and CDK13 impairs splicing of a subset of promoter-proximal introns characterized by weak 3' splice sites located at larger distance from the branchpoint. Nascent transcript analysis indicated that these introns are selectively retained upon pharmacological inhibition of CDK12/13 with respect to downstream introns of the same pre-mRNAs. Retention of these introns was also triggered by pladienolide B (PdB), an inhibitor of the U2 small nucelar ribonucleoprotein (snRNP) factor SF3B1 that recognizes the branchpoint. CDK12/13 activity promotes the interaction of SF3B1 with RNAPII phosphorylated on Ser2, and disruption of this interaction by treatment with the CDK12/13 inhibitor THZ531 impairs the association of SF3B1 with chromatin and its recruitment to the 3' splice site of these introns. Furthermore, by using suboptimal doses of THZ531 and PdB, we describe a synergic effect of these inhibitors on intron retention, cell cycle progression and cancer cell survival. These findings uncover a mechanism by which CDK12/13 couple RNA transcription and processing, and suggest that combined inhibition of these kinases and the spliceosome represents an exploitable anticancer approach.

DNA polymerases are essential for nucleic acid synthesis, cloning, sequencing and molecular diagnostics technologies. Conditional intein splicing is a powerful tool for controlling enzyme reactions. We have engineered a thermal switch into thermostable DNA polymerases from two structurally distinct polymerase families by inserting a thermally activated intein domain into a surface loop that is integral to the polymerase active site, thereby blocking DNA or RNA template access. The fusion proteins are inactive, but retain their structures, such that the intein excises during a heat pulse delivered at 70-80°C to generate spliced, active polymerases. This straightforward thermal activation step provides a highly effective, one-component 'hot-start' control of PCR reactions that enables accurate target amplification by minimizing unwanted by-products generated by off-target reactions. In one engineered enzyme, derived from Thermus aquaticus DNA polymerase, both DNA polymerase and reverse transcriptase activities are controlled by the intein, enabling single-reagent amplification of DNA and RNA under hot-start conditions. This engineered polymerase provides high-sensitivity detection for molecular diagnostics applications, amplifying 5-6 copies of the tested DNA and RNA targets with >95% certainty. The design principles used to engineer the inteins can be readily applied to construct other conditionally activated nucleic acid processing enzymes.

Dysfunction of the RNA-binding protein (RBP) FUS implicated in RNA metabolism can cause amyotrophic lateral sclerosis (ALS) and other neurodegenerative diseases. Mutations affecting FUS nuclear localization can drive RNA splicing defects and stimulate the formation of non-amyloid inclusions in affected neurons. However, the mechanism by which FUS mutations contribute to the development of ALS remains uncertain. Here we describe a pattern of RNA splicing changes in the dynamics of the continuous proteinopathy induced by mislocalized FUS. We show that the decrease in intron retention of FUS-associated transcripts represents the hallmark of the pathogenesis of ALS and is the earliest molecular event in the course of progression of the disease. As FUS aggregation increases, the pattern of RNA splicing changes, becoming more complex, including a decrease in the inclusion of neuron-specific microexons and induction of cryptic exon splicing due to the sequestration of additional RBPs into FUS aggregates. Crucially, the identified features of the pathological splicing pattern are also observed in ALS patients in both sporadic and familial cases. Our data provide evidence that both a loss of nuclear FUS function due to mislocalization and the subsequent cytoplasmic aggregation of mutant protein lead to the disruption of RNA splicing in a multistep fashion during FUS aggregation.

Gastrulation begins when the epiblast forms the primitive streak or becomes definitive ectoderm. During this lineage bifurcation, the DNA dioxygenase TET1 has bipartite functions in transcriptional activation and repression, but the mechanisms remain unclear. By converting mouse embryonic stem cells (ESCs) into neuroprogenitors, we defined how Tet1-/- cells switch from neuroectoderm fate to form mesoderm and endoderm. We identified the Wnt repressor Tcf7l1 as a TET1 target that suppresses Wnt/β-catenin and Nodal signalling. ESCs expressing catalytic dead TET1 retain neural potential but activate Nodal and subsequently Wnt/β-catenin pathways to generate also mesoderm and endoderm. At CpG-poor distal enhancers, TET1 maintains accessible chromatin at neuroectodermal loci independently of DNA demethylation. At CpG-rich promoters, DNA demethylation by TET1 affects the expression of bivalent genes. In ESCs, a non-catalytic TET1 cooperation with Polycomb represses primitive streak genes; post-lineage priming, the interaction becomes antagonistic at neuronal genes, when TET1's catalytic activity is further involved by repressing Wnt signalling. The convergence of repressive DNA and histone methylation does not inhibit neural induction in Tet1-deficient cells, but some DNA hypermethylated loci persist at genes with brain-specific functions. Our results reveal versatile switching of non-catalytic and catalytic TET1 activities based on genomic context, lineage and developmental stage.

Phosphorylation state-dependent interactions of the phosphoenolpyruvate (PEP):carbohydrate phosphotransferase system (PTS) components with transcription factors play a key role in carbon catabolite repression (CCR) by glucose in bacteria. Glucose inhibits the PTS-dependent transport of fructose and is preferred over fructose in Vibrio cholerae, but the mechanism is unknown. We have recently shown that, contrary to Escherichia coli, the fructose-dependent transcriptional regulator FruR acts as an activator of the fru operon in V. cholerae and binding of the FruR-fructose 1-phosphate (F1P) complex to an operator facilitates RNA polymerase (RNAP) binding to the fru promoter. Here we show that, in the presence of glucose, dephosphorylated HPr, a general PTS component, binds to FruR. Whereas HPr does not affect DNA-binding affinity of FruR, regardless of the presence of F1P, it prevents the FruR-F1P complex from facilitating the binding of RNAP to the fru promoter. Structural and biochemical analyses of the FruR-HPr complex identify key residues responsible for the V. cholerae-specific FruR-HPr interaction not observed in E. coli. Finally, we reveal how the dephosphorylated HPr interacts with FruR in V. cholerae, whereas the phosphorylated HPr binds to CcpA, which is a global regulator of CCR in Bacillus subtilis and shows structural similarity to FruR.

The human genome contains about 800 C2H2 zinc finger proteins (ZFPs), and most of them are composed of long arrays of zinc fingers. Standard ZFP recognition model asserts longer finger arrays should recognize longer DNA-binding sites. However, recent experimental efforts to identify in vivo ZFP binding sites contradict this assumption, with many exhibiting short motifs. Here we use ZFY, CTCF, ZIM3, and ZNF343 as examples to address three closely related questions: What are the reasons that impede current motif discovery methods? What are the functions of those seemingly unused fingers and how can we improve the motif discovery algorithms based on long ZFPs' biophysical properties? Using ZFY, we employed a variety of methods and find evidence for 'dependent recognition' where downstream fingers can recognize some previously undiscovered motifs only in the presence of an intact core site. For CTCF, high-throughput measurements revealed its upstream specificity profile depends on the strength of its core. Moreover, the binding strength of the upstream site modulates CTCF's sensitivity to different epigenetic modifications within the core, providing new insight into how the previously identified intellectual disability-causing and cancer-related mutant R567W disrupts upstream recognition and deregulates the epigenetic control by CTCF. Our results establish that, because of irregular motif structures, variable spacing and dependent recognition between sub-motifs, the specificities of long ZFPs are significantly underestimated, so we developed an algorithm, ModeMap, to infer the motifs and recognition models of ZIM3 and ZNF343, which facilitates high-confidence identification of specific binding sites, including repeats-derived elements. With revised concept, technique, and algorithm, we can discover the overlooked specificities and functions of those 'extra' fingers, and therefore decipher their broader roles in human biology and diseases.

The deubiquitinating enzyme Ataxin-3 (ATXN3) contains a polyglutamine (PolyQ) region, the expansion of which causes spinocerebellar ataxia type-3 (SCA3). ATXN3 has multiple functions, such as regulating transcription or controlling genomic stability after DNA damage. Here we report the role of ATXN3 in chromatin organization during unperturbed conditions, in a catalytic-independent manner. The lack of ATXN3 leads to abnormalities in nuclear and nucleolar morphology, alters DNA replication timing and increases transcription. Additionally, indicators of more open chromatin, such as increased mobility of histone H1, changes in epigenetic marks and higher sensitivity to micrococcal nuclease digestion were detected in the absence of ATXN3. Interestingly, the effects observed in cells lacking ATXN3 are epistatic to the inhibition or lack of the histone deacetylase 3 (HDAC3), an interaction partner of ATXN3. The absence of ATXN3 decreases the recruitment of endogenous HDAC3 to the chromatin, as well as the HDAC3 nuclear/cytoplasm ratio after HDAC3 overexpression, suggesting that ATXN3 controls the subcellular localization of HDAC3. Importantly, the overexpression of a PolyQ-expanded version of ATXN3 behaves as a null mutant, altering DNA replication parameters, epigenetic marks and the subcellular distribution of HDAC3, giving new insights into the molecular basis of the disease.