Integration into host target DNA (tDNA), a hallmark of retroviral replication, is mediated by the intasome, a multimer of integrase (IN) assembled on viral DNA (vDNA) ends. To ascertain aspects of tDNA recognition during integration, we have solved the 3.5 Å resolution cryo-EM structure of the mouse mammary tumor virus (MMTV) strand transfer complex (STC) intasome. The tDNA adopts an A-like conformation in the region encompassing the sites of vDNA joining, which exposes the sugar-phosphate backbone for IN-mediated strand transfer. Examination of existing retroviral STC structures revealed conservation of A-form tDNA in the analogous regions of these complexes. Furthermore, analyses of sequence preferences in genomic integration sites selectively targeted by six different retroviruses highlighted consistent propensity for A-philic sequences at the sites of vDNA joining. Our structure additionally revealed several novel MMTV IN-DNA interactions, as well as contacts seen in prior STC structures, including conserved Pro125 and Tyr149 residues interacting with tDNA. In infected cells, Pro125 substitutions impacted the global pattern of MMTV integration without significantly altering local base sequence preferences at vDNA insertion sites. Collectively, these data advance our understanding of retroviral intasome structure and function, as well as factors that influence patterns of vDNA integration in genomic DNA.

The genetic modification of microorganisms is conducive to the selection of high-yield producers of high-value-added chemicals, but a lack of genetic tools hinders the industrialization of most wild species. Therefore, it is crucial to develop host-independent gene editing tools that can be used for genetic manipulation-deprived strains. The Tn7-like transposon from Scytonema hofmanni has been shown to mediate homologous recombination-independent genomic integration after heterologous expression in Escherichia coli, but the integration efficiency of heterologous sequences larger than 5 kb remains suboptimal. Here, we constructed a versatile Cas12k-based genetic engineering toolkit (C12KGET) that can achieve genomic integration of fragments up to 10 kb in size with up to 100% efficiency in challenging strains. Using C12KGET, we achieved the first example of highly efficient genome editing in Sinorhizobium meliloti, which successfully solved the problem that industrial strains are difficult to genetically modify, and increased vitamin B12 production by 25%. In addition, Cas12k can be directly used for transcriptional regulation of genes with up to 92% efficiency due to its naturally inactivated nuclease domain. The C12KGET established in this study is a versatile and efficient marker-free tool for gene integration as well as transcriptional regulation that can be used for challenging strains with underdeveloped genetic toolkits.

Cold shock adaptability is a key survival skill of gut bacteria of warm-blooded animals. Escherichia coli cold shock responses are controlled by a complex multi-gene, timely-ordered transcriptional program. We investigated its underlying mechanisms. Having identified short-term, cold shock repressed genes, we show that their responsiveness is unrelated to their transcription factors or global regulators, while their single-cell protein numbers' variability increases after cold shock. We hypothesized that some cold shock repressed genes could be triggered by high propensity for transcription locking due to changes in DNA supercoiling (likely due to DNA relaxation caused by an overall reduction in negative supercoiling). Concomitantly, we found that nearly half of cold shock repressed genes are also highly responsive to gyrase inhibition (albeit most genes responsive to gyrase inhibition are not cold shock responsive). Further, their response strengths to cold shock and gyrase inhibition correlate. Meanwhile, under cold shock, nucleoid density increases, and gyrases and nucleoid become more colocalized. Moreover, the cellular energy decreases, which may hinder positive supercoils resolution. Overall, we conclude that sensitivity to diminished negative supercoiling is a core feature of E. coli's short-term, cold shock transcriptional program, and could be used to regulate the temperature sensitivity of synthetic circuits.

DNA methylation (5-methylcytosine (5mC)) is critical for genome stability and transcriptional regulation in mammals. The discovery that ten-eleven translocation (TET) proteins catalyze the oxidation of 5mC to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC) revolutionized our perspective on the complexity and regulation of DNA modifications. However, to what extent the regulatory functions of TET1 can be attributed to its catalytic activity remains unclear. Here, we use genome engineering and quantitative multi-omics approaches to dissect the precise catalytic vs. non-catalytic functions of TET1 in murine embryonic stem cells (mESCs). Our study identifies TET1 as an essential interaction hub for multiple chromatin modifying complexes and a global regulator of histone modifications. Strikingly, we find that the majority of transcriptional regulation depends on non-catalytic functions of TET1. In particular, we show that TET1 is critical for the establishment of H3K9me3 and H4K20me3 at endogenous retroviral elements (ERVs) and their silencing that is independent of its canonical role in DNA demethylation. Furthermore, we provide evidence that this repression of ERVs depends on the interaction between TET1 and SIN3A. In summary, we demonstrate that the non-catalytic functions of TET1 are critical for regulation of gene expression and the silencing of endogenous retroviruses in mESCs.

A series of multiple logic circuits based on a single biomolecular platform is constructed to perform nonarithmetic and arithmetic functions, including 4-to-2 encoder, 1-to-2 demultiplexer, 1-to-4 demultiplexer, and multi-input OR gate. The encoder to a DNA circuit is the equivalent of a sensory receptor to a reflex arc. They all function to encode information from outside the pathway (DNA circuit or reflex arc) into a form that subsequent pathways can recognize and utilize. Current molecular encoders are based on optical or electrical signals as outputs, while DNA circuits are based on DNA strands as transmission signals. The output of existing encoders cannot be recognized by subsequent DNA circuits. It is the first time the DNA-based encoder with DNA strands as outputs can be truly applied to the DNA circuit, enabling the application of DNA circuits in non-binary biological environments. Another novel feature of the designed system is that the developed nanodevices all have a simple structure, low leakage and low crosstalk, which allows them to implement higher-level encoders and demultiplexers easily. Our work is based on the idea of complex functionality in a simple form, which will also provide a new route for developing advanced molecular logic circuits.

Mitochondrial DNA has been investigated for nearly fifty years, but many aspects of the maintenance of this essential small genome remain unknown. Like any genome, mammalian mitochondrial DNA requires the function of topoisomerases to counter and regulate the topological tension arising during replication, transcription, segregation, and repair. However, the functions of the different mitochondrial topoisomerases are poorly understood. Here, we investigate the role of Topoisomerase 3α (Top3α) in mtDNA replication and transcription, providing evidence that this enzyme, previously reported to act in mtDNA segregation, also participates in mtDNA replication fork progression. Top3α knockdown caused replication fork stalling, increased mtDNA catenation and decreased mtDNA levels. Overexpression in contrast induced abundant double-strand breaks around the replication origin OH and abortion of early replication, while at the same time improving the resolution of mtDNA replication termination intermediates. Both Top3α knockdown and overexpression affected mitochondrial RNA transcription, leading to a decrease in steady-state levels of mitochondrial transcripts. Together, our results indicate that the mitochondrial isoform of Top3α is not only involved in mtDNA segregation, as reported previously, but also supports the progression of the replication fork. Mitochondrial Top3α is also influencing the progression of transcription, with its absence affecting downstream transcript levels.

The use of multiple drugs simultaneously targeting DNA is a promising strategy in cancer therapy for potentially overcoming single drug resistance. In support of this concept, we report that a combination of actinomycin D (ActD) and echinomycin (Echi), can interact in novel ways with native and mismatched DNA sequences, distinct from the structural effects produced by either drug alone. Changes in the former with GpC and CpG steps separated by a A:G or G:A mismatch or in a native DNA with canonical G:C and C:G base pairs, result in significant asymmetric backbone twists through staggered intercalation and base pair modulations. A wobble or Watson-Crick base pair at the two drug-binding interfaces can result in a single-stranded 'chair-shaped' DNA duplex with a straight helical axis. However, a novel sugar-edged hydrogen bonding geometry in the G:A mismatch leads to a 'curved-shaped' duplex. Two non-canonical G:C Hoogsteen base pairings produce a sharply kinked duplex in different forms and a four-way junction-like superstructure, respectively. Therefore, single base pair modulations on the two drug-binding interfaces could significantly affect global DNA structure. These structures thus provide a rationale for atypical DNA recognition via multiple DNA intercalators and a structural basis for the drugs' potential synergetic use.

The combination of genome-editing and epitope tagging provides a powerful strategy to study proteins with high affinity and specificity while preserving their physiological expression patterns. However, stably modifying endogenous genes in cells that do not allow for clonal selection has been challenging. Here, we present a simple and fast strategy to generate stable, endogenously tagged alleles in a non-transformed cell culture model. At the example of piwi in Drosophila ovarian somatic sheath cells, we show that this strategy enables the generation of an N-terminally tagged protein that emulates the expression level and subcellular localization of the wild type protein and forms functional Piwi-piRNA complexes. We further present a concise workflow to establish endogenously N-terminally and C-terminally tagged proteins, and knockout alleles through rapid selection of cell pools in fly and human models.

Nonsense-mediated mRNA decay (NMD) protects cells from the toxic and potentially dominant effects of truncated proteins. Targeting of mRNAs with early stop codons is mediated by the ribosome and spatiotemporally aligned with translation termination. Previously we identified a novel NMD intermediate: ribosomes stalled on cleaved stop codons, raising the possibility that NMD begins even prior to ribosome removal from the stop codon. Here we show that this intermediate is the result of mRNA cleavage by the endonuclease SMG-6. Our work supports a model in which ribosomes stall secondary to SMG-6 mRNA cleavage in Caenorhabditis elegans and humans, i.e. that the novel NMD intermediate occurs after a prior ribosome elicits NMD. Our genetic analysis of C. elegans' SMG-6 supports a central role for SMG-6 in metazoan NMD, and provides a context for evaluating its function in other metazoans.

Correct pre-mRNA processing in higher eukaryotes vastly depends on splice site recognition. Beyond conserved 5'ss and 3'ss motifs, splicing regulatory elements (SREs) play a pivotal role in this recognition process. Here, we present in silico designed sequences with arbitrary a priori prescribed splicing regulatory HEXplorer properties that can be concatenated to arbitrary length without changing their regulatory properties. We experimentally validated in silico predictions in a massively parallel splicing reporter assay on more than 3000 sequences and exemplarily identified some SRE binding proteins. Aiming at a unified 'functional splice site strength' encompassing both U1 snRNA complementarity and impact from neighboring SREs, we developed a novel RNA-seq based 5'ss usage landscape, mapping the competition of pairs of high confidence 5'ss and neighboring exonic GT sites along HBond and HEXplorer score coordinate axes on human fibroblast and endothelium transcriptome datasets. These RNA-seq data served as basis for a logistic 5'ss usage prediction model, which greatly improved discrimination between strong but unused exonic GT sites and annotated highly used 5'ss. Our 5'ss usage landscape offers a unified view on 5'ss and SRE neighborhood impact on splice site recognition, and may contribute to improved mutation assessment in human genetics.

Mitochondrial DNA (mtDNA) is prone to mutation in aging and over evolutionary time, yet the processes that regulate the accumulation of de novo mtDNA mutations and modulate mtDNA heteroplasmy are not fully elucidated. Mitochondria lack certain DNA repair processes, which could contribute to polymerase error-induced mutations and increase susceptibility to chemical-induced mtDNA mutagenesis. We conducted error-corrected, ultra-sensitive Duplex Sequencing to investigate the effects of two known nuclear genome mutagens, cadmium and Aflatoxin B1, on germline mtDNA mutagenesis in Caenorhabditis elegans. Detection of thousands of mtDNA mutations revealed pervasive heteroplasmy in C. elegans and that mtDNA mutagenesis is dominated by C:G → A:T mutations generally attributed to oxidative damage. However, there was no effect of either exposure on mtDNA mutation frequency, spectrum, or trinucleotide context signature despite a significant increase in nuclear mutation rate after aflatoxin B1 exposure. Mitophagy-deficient mutants pink-1 and dct-1 accumulated significantly higher levels of mtDNA damage compared to wild-type C. elegans after exposures. However, there were only small differences in mtDNA mutation frequency, spectrum, or trinucleotide context signature compared to wild-type after 3050 generations, across all treatments. These findings suggest mitochondria harbor additional previously uncharacterized mechanisms that regulate mtDNA mutational processes across generations.

Many essential cellular processes rely on substrate rotation or translocation by a multi-subunit, ring-type NTPase. A large number of double-stranded DNA viruses, including tailed bacteriophages and herpes viruses, use a homomeric ring ATPase to processively translocate viral genomic DNA into procapsids during assembly. Our current understanding of viral DNA packaging comes from three archetypal bacteriophage systems: cos, pac and phi29. Detailed mechanistic understanding exists for pac and phi29, but not for cos. Here, we reconstituted in vitro a cos packaging system based on bacteriophage HK97 and provided a detailed biochemical and structural description. We used a photobleaching-based, single-molecule assay to determine the stoichiometry of the DNA-translocating ATPase large terminase. Crystal structures of the large terminase and DNA-recruiting small terminase, a first for a biochemically defined cos system, reveal mechanistic similarities between cos and pac systems. At the same time, mutational and biochemical analyses indicate a new regulatory mechanism for ATPase multimerization and coordination in the HK97 system. This work therefore establishes a framework for studying the evolutionary relationships between ATP-dependent DNA translocation machineries in double-stranded DNA viruses.

CRISPR-Cas systems are bacterial defense systems for fighting against invaders such as bacteriophages and mobile genetic elements. To escape destruction by these bacterial immune systems, phages have co-evolved multiple anti-CRISPR (Acr) proteins, which inhibit CRISPR-Cas function. Many acr genes form an operon with genes encoding transcriptional regulators, called anti-CRISPR-associated (Aca) proteins. Aca10 is the most recently discovered Aca family that is encoded within an operon containing acrIC7 and acrIC6 in Pseudomonas citronellolis. Here, we report the high-resolution crystal structure of an Aca10 protein to unveil the molecular basis of transcriptional repressor role of Aca10 in the acrIC7-acrIC6-aca10 operon. We identified that Aca10 forms a dimer in solution, which is critical for binding specific DNA. We also showed that Aca10 directly recognizes a 21 bp palindromic sequence in the promoter of the acr operon. Finally, we revealed that R44 of Aca10 is a critical residue involved in the DNA binding, which likely results in a high degree of DNA bending.

Despite advances in bacterial genome engineering, delivery of large synthetic constructs remains challenging in practice. In this study, we propose a straightforward and robust approach for the markerless integration of DNA fragments encoding whole metabolic pathways into the genome. This approach relies on the replacement of a counterselection marker with cargo DNA cassettes via λRed recombineering. We employed a counterselection strategy involving a genetic circuit based on the CI repressor of λ phage. Our design ensures elimination of most spontaneous mutants, and thus provides a counterselection stringency close to the maximum possible. We improved the efficiency of integrating long PCR-generated cassettes by exploiting the Ocr antirestriction function of T7 phage, which completely prevents degradation of unmethylated DNA by restriction endonucleases in wild-type bacteria. The employment of highly restrictive counterselection and ocr-assisted λRed recombineering allowed markerless integration of operon-sized cassettes into arbitrary genomic loci of four enterobacterial species with an efficiency of 50-100%. In the case of Escherichia coli, our strategy ensures simple combination of markerless mutations in a single strain via P1 transduction. Overall, the proposed approach can serve as a general tool for synthetic biology and metabolic engineering in a range of bacterial hosts.

Staphylococcus aureus, a human opportunist pathogen, adjusts its metabolism to cope with iron deprivation within the host. We investigated the potential role of small non-coding RNAs (sRNAs) in dictating this process. A single sRNA, named here IsrR, emerged from a competition assay with tagged-mutant libraries as being required during iron starvation. IsrR is iron-repressed and predicted to target mRNAs expressing iron-containing enzymes. Among them, we demonstrated that IsrR down-regulates the translation of mRNAs of enzymes that catalyze anaerobic nitrate respiration. The IsrR sequence reveals three single-stranded C-rich regions (CRRs). Mutational and structural analysis indicated a differential contribution of these CRRs according to targets. We also report that IsrR is required for full lethality of S. aureus in a mouse septicemia model, underscoring its role as a major contributor to the iron-sparing response for bacterial survival during infection. IsrR is conserved among staphylococci, but it is not ortholog to the proteobacterial sRNA RyhB, nor to other characterized sRNAs down-regulating mRNAs of iron-containing enzymes. Remarkably, these distinct sRNAs regulate common targets, illustrating that RNA-based regulation provides optimal evolutionary solutions to improve bacterial fitness when iron is scarce.

Identification of metabolite caps including FAD on the 5' end of RNA has uncovered a previously unforeseen intersection between cellular metabolism and gene expression. To understand the function of FAD caps in cellular physiology, we characterised the proteins interacting with FAD caps in budding yeast. Here we demonstrate that highly conserved 5'-3' exoribonucleases, Xrn1 and Rat1, physically interact with the RNA 5' FAD cap and both possess FAD cap decapping (deFADding) activity and subsequently degrade the resulting RNA. Xrn1 deFADding activity was also evident in human cells indicating its evolutionary conservation. Furthermore, we report that the recently identified bacterial 5'-3' exoribonuclease RNase AM also possesses deFADding activity that can degrade FAD-capped RNAs in vitro and in Escherichia coli cells. To gain a molecular understanding of the deFADding reaction, an RNase AM crystal structure with three manganese ions coordinated by a sulfate molecule and the active site amino acids was generated that provided details underlying hydrolysis of the FAD cap. Our findings reveal a general propensity for 5'-3' exoribonucleases to hydrolyse and degrade RNAs with 5' end noncanonical caps in addition to their well characterized 5' monophosphate RNA substrates indicating an intrinsic property of 5'-3' exoribonucleases.

Recent in vitro reconstitution analyses have proven that the physical interaction between the exosome core and MTR4 helicase, which promotes the exosome activity, is maintained by either MPP6 or RRP6. However, knowledge regarding the function of MPP6 with respect to in vivo exosome activity remains scarce. Here, we demonstrate a facilitative function of MPP6 that composes a specific part of MTR4-dependent substrate decay by the human exosome. Using RNA polymerase II-transcribed poly(A)+ substrate accumulation as an indicator of a perturbed exosome, we found functional redundancy between RRP6 and MPP6 in the decay of these poly(A)+ transcripts. MTR4 binding to the exosome core via MPP6 was essential for MPP6 to exert its redundancy with RRP6. However, at least for the decay of our identified exosome substrates, MTR4 recruitment by MPP6 was not functionally equivalent to recruitment by RRP6. Genome-wide classification of substrates based on their sensitivity to each exosome component revealed that MPP6 deals with a specific range of substrates and highlights the importance of MTR4 for their decay. Considering recent findings of competitive binding to the exosome between auxiliary complexes, our results suggest that the MPP6-incorporated MTR4-exosome complex is one of the multiple alternative complexes rather than the prevailing one.

PLEs (phage-inducible chromosomal island-like elements) are phage parasites integrated into the chromosome of epidemic Vibrio cholerae. In response to infection by its viral host ICP1, PLE excises, replicates and hijacks ICP1 structural components for transduction. Through an unknown mechanism, PLE prevents ICP1 from transitioning to rolling circle replication (RCR), a prerequisite for efficient packaging of the viral genome. Here, we characterize a PLE-encoded nuclease, NixI, that blocks phage development likely by nicking ICP1's genome as it transitions to RCR. NixI-dependent cleavage sites appear in ICP1's genome during infection of PLE(+) V. cholerae. Purified NixI demonstrates in vitro nuclease activity specifically for sites in ICP1's genome and we identify a motif that is necessary for NixI-mediated cleavage. Importantly, NixI is sufficient to limit ICP1 genome replication and eliminate progeny production, representing the most inhibitory PLE-encoded mechanism revealed to date. We identify distant NixI homologs in an expanded family of putative phage parasites in vibrios that lack nucleotide homology to PLEs but nonetheless share genomic synteny with PLEs. More generally, our results reveal a previously unknown mechanism deployed by phage parasites to limit packaging of their viral hosts' genome and highlight the prominent role of nuclease effectors as weapons in the arms race between antagonizing genomes.

Many bacteria use CRISPR-Cas systems to defend against invasive mobile genetic elements (MGEs). In response, MGEs have developed strategies to resist CRISPR-Cas, including the use of anti-CRISPR (Acr) proteins. Known acr genes may be followed in an operon by a putative regulatory Acr-associated gene (aca), suggesting the importance of regulation. Although ten families of helix-turn-helix (HTH) motif containing Aca proteins have been identified (Aca1-10), only three have been tested and shown to be transcriptional repressors of acr-aca expression. The AcrIIA1 protein (a Cas9 inhibitor) also contains a functionally similar HTH containing repressor domain. Here, we identified and analysed Aca and AcrIIA1 homologs across all bacterial genomes. Using HMM models we found aca-like genes are widely distributed in bacteria, both with and without known acr genes. The putative promoter regions of acr-aca operons were analysed and members of each family of bacterial Aca tested for regulatory function. For each Aca family, we predicted a conserved inverted repeat binding site within a core promoter. Promoters containing these sites directed reporter expression in E. coli and were repressed by the cognate Aca protein. These data demonstrate that acr repression by Aca proteins is widely conserved in nature.

The bacterial FtsK motor harvests energy from ATP to translocate double-stranded DNA during cell division. Here, we probe the molecular mechanisms underlying coordinated DNA translocation in FtsK by performing long timescale simulations of its hexameric assembly and individual subunits. From these simulations we predict signaling pathways that connect the ATPase active site to DNA-gripping residues, which allows the motor to coordinate its translocation activity with its ATPase activity. Additionally, we utilize well-tempered metadynamics simulations to compute free-energy landscapes that elucidate the extended-to-compact transition involved in force generation. We show that nucleotide binding promotes a compact conformation of a motor subunit, whereas the apo subunit is flexible. Together, our results support a mechanism whereby each ATP-bound subunit of the motor conforms to the helical pitch of DNA, and ATP hydrolysis/product release causes a subunit to lose grip of DNA. By ordinally engaging and disengaging with DNA, the FtsK motor unidirectionally translocates DNA.

G-rich sequences found at multiple sites throughout all genomes may form secondary structures called G-quadruplexes (G4), which act as roadblocks for molecular motors. Among the enzymes thought to process these structures, the Pif1 DNA helicase is considered as an archetypical G4-resolvase and its absence has been linked to G4-related genomic instabilities in yeast. Here we developed a single-molecule assay to observe Pif1 opening a DNA duplex and resolving the G4 in real time. In support of former enzymological studies, we show that the helicase reduces the lifetime of G4 from hours to seconds. However, we observe that in the presence of a G4, Pif1 exhibits a strong strand switching behavior, which can lead to Pif1 escaping G4 resolution, depending on the structural context surrounding the substrate. This behavior is also detected in the presence of other roadblocks (LNA or RNA). We propose that the efficiency of Pif1 to remove a roadblock (G4 or other) is affected by its strand switching behavior and depends on the context surrounding the obstacle. We discuss how this switching behavior may explain several aspects of Pif1 substrate preference and affect its activity as a G4 resolvase in vivo.

The in vivo role for RNase H1 in mammalian mitochondria has been much debated. Loss of RNase H1 is embryonic lethal and to further study its role in mtDNA expression we characterized a conditional knockout of Rnaseh1 in mouse heart. We report that RNase H1 is essential for processing of RNA primers to allow site-specific initiation of mtDNA replication. Without RNase H1, the RNA:DNA hybrids at the replication origins are not processed and mtDNA replication is initiated at non-canonical sites and becomes impaired. Importantly, RNase H1 is also needed for replication completion and in its absence linear deleted mtDNA molecules extending between the two origins of mtDNA replication are formed accompanied by mtDNA depletion. The steady-state levels of mitochondrial transcripts follow the levels of mtDNA, and RNA processing is not altered in the absence of RNase H1. Finally, we report the first patient with a homozygous pathogenic mutation in the hybrid-binding domain of RNase H1 causing impaired mtDNA replication. In contrast to catalytically inactive variants of RNase H1, this mutant version has enhanced enzyme activity but shows impaired primer formation. This finding shows that the RNase H1 activity must be strictly controlled to allow proper regulation of mtDNA replication.

Defining the impact of missense mutations on the recognition of DNA motifs is highly dependent on bioinformatic tools that define DNA binding elements. However, classical motif analysis tools remain limited in their capacity to identify subtle changes in complex binding motifs between distinct conditions. To overcome this limitation, we developed a new tool, MoMotif, that facilitates a sensitive identification, at the single base-pair resolution, of complex, or subtle, alterations to core binding motifs, discerned from ChIP-seq data. We employed MoMotif to define the previously uncharacterized recognition motif of CTCF zinc-finger 1 (ZF1), and to further define the impact of CTCF ZF1 mutation on its association with chromatin. Mutations of CTCF ZF1 are exclusive to breast cancer and are associated with metastasis and therapeutic resistance, but the underlying mechanisms are unclear. Using MoMotif, we identified an extension of the CTCF core binding motif, necessitating a functional ZF1 to bind appropriately. Using a combination of ChIP-Seq and RNA-Seq, we discover that the inability to bind this extended motif drives an altered transcriptional program associated with the oncogenic phenotypes observed clinically. Our study demonstrates that MoMotif is a powerful new tool for comparative ChIP-seq analysis and characterising DNA-protein contacts.

The lung is a complex organ with various cell types having distinct roles. Antisense oligonucleotides (ASOs) have been studied in the lung, but it has been challenging to determine their effectiveness in each cell type due to the lack of appropriate analytical methods. We employed three distinct approaches to study silencing efficacy within different cell types. First, we used lineage markers to identify cell types in flow cytometry, and simultaneously measured ASO-induced silencing of cell-surface proteins CD47 or CD98. Second, we applied single-cell RNA sequencing (scRNA-seq) to measure silencing efficacy in distinct cell types; to the best of our knowledge, this is the first time scRNA-seq has been applied to measure the efficacy of oligonucleotide therapeutics. In both approaches, fibroblasts were the most susceptible to locally delivered ASOs, with significant silencing also in endothelial cells. Third, we confirmed that the robust silencing in fibroblasts is broadly applicable by silencing two targets expressed mainly in fibroblasts, Mfap4 and Adam33. Across independent approaches, we demonstrate that intratracheally administered LNA gapmer ASOs robustly induce gene silencing in lung fibroblasts. ASO-induced gene silencing in fibroblasts was durable, lasting 4-8 weeks after a single dose. Thus, lung fibroblasts are well aligned with ASOs as therapeutics.

The transcription factor SOX9 is activated at the onset of endothelial-to-mesenchymal transition (EndMT) during embryonic development and in pathological conditions. Its roles in regulating these processes, however, are not clear. Using human umbilical vein endothelial cells (HUVECs) as an EndMT model, we show that SOX9 expression alone is sufficient to activate mesenchymal genes and steer endothelial cells towards a mesenchymal fate. By genome-wide mapping of the chromatin landscape, we show that SOX9 displays features of a pioneer transcription factor, such as opening of chromatin and leading to deposition of active histone modifications at silent chromatin regions, guided by SOX dimer motifs and H2A.Z enrichment. We further observe highly transient and dynamic SOX9 binding, possibly promoted through its eviction by histone phosphorylation. However, while SOX9 binding is dynamic, changes in the chromatin landscape and cell fate induced by SOX9 are persistent. Finally, our analysis of single-cell chromatin accessibility indicates that SOX9 opens chromatin to drive EndMT in atherosclerotic lesions in vivo. This study provides new insight into key molecular functions of SOX9 and mechanisms of EndMT and highlights the crucial developmental role of SOX9 and relevance to human disease.

Long interspersed nuclear element 1 (L1) parasitized most vertebrates and constitutes ∼20% of the human genome. It encodes ORF1p and ORF2p which form an L1-ribonucleoprotein (RNP) with their encoding transcript that is copied into genomic DNA (retrotransposition). ORF1p binds single-stranded nucleic acid (ssNA) and exhibits NA chaperone activity. All vertebrate ORF1ps contain a coiled coil (CC) domain and we previously showed that a CC-retrotransposition null mutant prevented formation of stably bound ORF1p complexes on ssNA. Here, we compared CC variants using our recently improved method that measures ORF1p binding to ssDNA at different forces. Bound proteins decrease ssDNA contour length and at low force, retrotransposition-competent ORF1ps (111p and m14p) exhibit two shortening phases: the first is rapid, coincident with ORF1p binding; the second is slower, consistent with formation of tightly compacted complexes by NA-bound ORF1p. In contrast, two retrotransposition-null CC variants (151p and m15p) did not attain the second tightly compacted state. The C-terminal half of the ORF1p trimer (not the CC) contains the residues that mediate NA-binding. Our demonstrating that the CC governs the ability of NA-bound retrotransposition-competent trimers to form tightly compacted complexes reveals the biochemical phenotype of these coiled coil mutants.

Genome wide association studies provide statistical measures of gene-trait associations that reveal how genetic variation influences phenotypes. This study develops an unsupervised dimensionality reduction method called UnTANGLeD (Unsupervised Trait Analysis of Networks from Gene Level Data) which organizes 16,849 genes into discrete gene programs by measuring the statistical association between genetic variants and 1,393 diverse complex traits. UnTANGLeD reveals 173 gene clusters enriched for protein-protein interactions and highly distinct biological processes governing development, signalling, disease, and homeostasis. We identify diverse gene networks with robust interactions but not associated with known biological processes. Analysis of independent disease traits shows that UnTANGLeD gene clusters are conserved across all complex traits, providing a simple and powerful framework to predict novel gene candidates and programs influencing orthogonal disease phenotypes. Collectively, this study demonstrates that gene programs co-ordinately orchestrating cell functions can be identified without reliance on prior knowledge, providing a method for use in functional annotation, hypothesis generation, machine learning and prediction algorithms, and the interpretation of diverse genomic data.

Topologically associated domains (TADs) are crucial chromatin structural units. Evidence has illustrated that RNA-chromatin and RNA-RNA spatial interactions, so-called RNA-associated interactions (RAIs), may be associated with TAD-like domains (TLDs). To decode hierarchical TLDs from RAIs, we proposed SuperTLD, a domain detection algorithm incorporating imputation. We applied SuperTLD on four RAI data sets and compared TLDs with the TADs identified from the corresponding Hi-C datasets. The TLDs and TADs share a moderate similarity of hierarchies ≥ 0.5312 and the finest structures ≥ 0.8295. Comparison between boundaries and domains further demonstrated the novelty of TLDs. Enrichment analysis of epigenetic characteristics illustrated that the novel TLDs exhibit an enriched CTCF by 0.6245 fold change and H3 histone marks enriched within domains. GO analysis on the TLD novel boundaries exhibited enriched diverse terms, revealing TLDs' formation mechanism related closely to gene regulation.

Despite recent advances in inferring cellular dynamics using single-cell RNA-seq data, existing trajectory inference (TI) methods face difficulty in accurately reconstructing the cell-state manifold and cell-fate plasticity for complex topologies. Here, we present MARGARET (https://github.com/Zafar-Lab/Margaret) for inferring single-cell trajectory and fate mapping for diverse dynamic cellular processes. MARGARET reconstructs complex trajectory topologies using a deep unsupervised metric learning and a graph-partitioning approach based on a novel connectivity measure, automatically detects terminal cell states, and generalizes the quantification of fate plasticity for complex topologies. On a diverse benchmark consisting of synthetic and real datasets, MARGARET outperformed state-of-the-art methods in recovering global topology and cell pseudotime ordering. For human hematopoiesis, MARGARET accurately identified all major lineages and associated gene expression trends and helped identify transitional progenitors associated with key branching events. For embryoid body differentiation, MARGARET identified novel transitional populations that were validated by bulk sequencing and functionally characterized different precursor populations in the mesoderm lineage. For colon differentiation, MARGARET characterized the lineage for BEST4/OTOP2 cells and the heterogeneity in goblet cell lineage in the colon under normal and inflamed ulcerative colitis conditions. Finally, we demonstrated that MARGARET can scale to large scRNA-seq datasets consisting of ∼ millions of cells.

PCI domain proteins play important roles in post-transcriptional gene regulation. In the TREX-2 complex, PCI domain-containing Sac3 and Thp1 proteins and accessory Sem1 protein form a ternary complex required for mRNA nuclear export. In contrast, structurally related Thp3-Csn12-Sem1 complex mediates pre-mRNA splicing. In this study, we determined the structure of yeast Thp3186-470-Csn12-Sem1 ternary complex at 2.9 Å resolution. Both Thp3 and Csn12 structures have a typical PCI structural fold, characterized by a stack of α-helices capped by a C-terminal winged-helix (WH) domain. The overall structure of Thp3186-470-Csn12-Sem1 complex has an inverted V-shape with Thp3 and Csn12 forming the two sides. A fishhook-shaped Sem1 makes extensive contacts on Csn12 to stabilize its conformation. The overall structure of Thp3186-470-Csn12-Sem1 complex resembles the previously reported Sac3-Thp1-Sem1 complex, but also has significant structural differences. The C-terminal WH domains of Thp3 and Csn12 form a continuous surface to bind different forms of nucleic acids with micromolar affinity. Mutation of the basic residues in the WH domains of Thp3 and Csn12 affects nucleic acid binding in vitro and mRNA splicing in vivo. The Thp3-Csn12-Sem1 structure provides a foundation for further exploring the structural elements required for its specific recruitment to spliceosome for pre-mRNA splicing.

Bacterial adaptation is largely shaped by horizontal gene transfer, xenogeneic silencing mediated by lineage-specific DNA bridgers (H-NS, Lsr2, MvaT and Rok), and various anti-silencing mechanisms. No xenogeneic silencing DNA bridger is known for α-proteobacteria, from which mitochondria evolved. By investigating α-proteobacterium Sinorhizobium fredii, a facultative legume microsymbiont, here we report the conserved zinc-finger bearing MucR as a novel xenogeneic silencing DNA bridger. Self-association mediated by its N-terminal domain (NTD) is required for DNA-MucR-DNA bridging complex formation, maximizing MucR stability, transcriptional silencing, and efficient symbiosis in legume nodules. Essential roles of NTD, CTD (C-terminal DNA-binding domain), or full-length MucR in symbiosis can be replaced by non-homologous NTD, CTD, or full-length protein of H-NS from γ-proteobacterium Escherichia coli, while NTD rather than CTD of Lsr2 from Gram-positive Mycobacterium tuberculosis can replace the corresponding domain of MucR in symbiosis. Chromatin immunoprecipitation sequencing reveals similar recruitment profiles of H-NS, MucR and various functional chimeric xenogeneic silencers across the multipartite genome of S. fredii, i.e. preferring AT-rich genomic islands and symbiosis plasmid with key symbiosis genes as shared targets. Collectively, the convergently evolved DNA bridger MucR predisposed α-proteobacteria to integrate AT-rich foreign DNA including symbiosis genes, horizontal transfer of which is strongly selected in nature.

Tools for synthetically controlling gene expression are a cornerstone of genetic engineering. CRISPRi and CRISPRa technologies have been applied extensively for programmable modulation of gene transcription, but there are few such tools for targeted modulation of protein translation rates. Here, we employ CRISPR-Cas13 as a programmable activator of translation. We develop a novel variant of the catalytically-deactivated Cas13d enzyme dCasRx by fusing it to translation initiation factor IF3. We demonstrate dCasRx-IF3's ability to enhance expression 21.3-fold above dCasRx when both are targeted to the start of the 5' untranslated region of mRNA encoding red fluorescent protein in Escherichia coli. Activation of translation is location-dependent, and we show dCasRx-IF3 represses translation when targeted to the ribosomal binding site, rather than enhancing it. We provide evidence that dCasRx-IF3 targeting enhances mRNA stability relative to dCasRx, providing mechanistic insights into how this new tool functions to enhance gene expression. We also demonstrate targeted upregulation of native LacZ 2.6-fold, showing dCasRx-IF3's ability to enhance expression of endogenous genes. dCasRx-IF3 requires no additional host modification to influence gene expression. This work outlines a novel approach, CRISPR-RNAa, for post-transcriptional control of translation to activate gene expression.

Information processing functions are essential for organisms to perceive and react to their complex environment, and for humans to analyze and rationalize them. While our brain is extraordinary at processing complex information, winner-take-all, as a type of biased competition is one of the simplest models of lateral inhibition and competition among biological neurons. It has been implemented as DNA-based neural networks, for example, to mimic pattern recognition. However, the utility of DNA-based computation in information processing for real biotechnological applications remains to be demonstrated. In this paper, a biased competition method for nonlinear manipulation and analysis of mixtures of DNA sequences was developed. Unlike conventional biological experiments, selected species were not directly subjected to analysis. Instead, parallel computation among a myriad of different DNA sequences was carried out to reduce the information entropy. The method could be used for various oligonucleotide-encoded libraries, as we have demonstrated its application in decoding and data analysis for selection experiments with DNA-encoded chemical libraries against protein targets.

Regulation of gene expression is a key factor influencing the success of antimicrobial resistance determinants. A variety of determinants conferring resistance against aminoglycosides (Ag) are commonly found in clinically relevant bacteria, but whether their expression is regulated or not is controversial. The expression of several Ag resistance genes has been reported to be controlled by a riboswitch mechanism encoded in a conserved sequence. Yet this sequence corresponds to the integration site of an integron, a genetic platform that recruits genes of different functions, making the presence of such a riboswitch counterintuitive. We provide, for the first time, experimental evidence against the existence of such Ag-sensing riboswitch. We first tried to reproduce the induction of the well characterized aacA5 gene using its native genetic environment, but were unsuccessful. We then broadened our approach and analyzed the inducibility of all AgR genes encoded in integrons against a variety of antibiotics. We could not observe biologically relevant induction rates for any gene in the presence of several aminoglycosides. Instead, unrelated antibiotics produced mild but consistently higher increases in expression, that were the result of pleiotropic effects. Our findings rule out the riboswitch control of aminoglycoside resistance genes in integrons.

SRRM2 is a nuclear-speckle marker containing multiple disordered domains, whose dysfunction is associated with several human diseases. Using mainly EGFP-SRRM2 knock-in HEK293T cells, we show that SRRM2 forms biomolecular condensates satisfying most hallmarks of liquid-liquid phase separation, including spherical shape, dynamic rearrangement, coalescence and concentration dependence supported by in vitro experiments. Live-cell imaging shows that SRRM2 organizes nuclear speckles along the cell cycle. As bona-fide splicing factor present in spliceosome structures, SRRM2 deficiency induces skipping of cassette exons with short introns and weak splice sites, tending to change large protein domains. In THP-1 myeloid-like cells, SRRM2 depletion compromises cell viability, upregulates differentiation markers, and sensitizes cells to anti-leukemia drugs. SRRM2 induces a FES splice isoform that attenuates innate inflammatory responses, and MUC1 isoforms that undergo shedding with oncogenic properties. We conclude that SRRM2 acts as a scaffold to organize nuclear speckles, regulating alternative splicing in innate immunity and cell homeostasis.

Replication of the circular bacterial chromosome is initiated from a locus oriC with the aid of an essential protein DnaA. One approach to identify factors acting to prevent aberrant oriC-independent replication initiation in Escherichia coli has been that to obtain mutants which survive loss of DnaA. Here, we show that a ΔrecD mutation, associated with attenuation of RecBCD's DNA double strand end-resection activity, provokes abnormal replication and rescues ΔdnaA lethality in two situations: (i) in absence of 5'-3' single-strand DNA exonuclease RecJ, or (ii) when multiple two-ended DNA double strand breaks (DSBs) are generated either by I-SceI endonucleolytic cleavages or by radiomimetic agents phleomycin or bleomycin. One-ended DSBs in the ΔrecD mutant did not rescue ΔdnaA lethality. With two-ended DSBs in the ΔrecD strain, ΔdnaA viability was retained even after linearization of the chromosome. Data from genome-wide DNA copy number determinations in ΔdnaA-rescued cells lead us to propose a model that nuclease-mediated DNA resection activity of RecBCD is critical for prevention of a σ-mode of rolling-circle over-replication when convergent replication forks merge and fuse, as may be expected to occur during normal replication at the chromosomal terminus region or during repair of two-ended DSBs following 'ends-in' replication.

Although both the p53 and forkhead box (FOX) family proteins are key transcription factors associated with cancer progression, their direct relationship is unknown. Here, we found that FOX family proteins bind to the non-canonical homotypic cluster of the p53 promoter region (TP53). Analysis of crystal structures of FOX proteins (FOXL2 and FOXA1) bound to the p53 homotypic cluster indicated that they interact with a 2:1 stoichiometry accommodated by FOX-induced DNA allostery. In particular, FOX proteins exhibited distinct dimerization patterns in recognition of the same p53-DNA; dimer formation of FOXA1 involved protein-protein interaction, but FOXL2 did not. Biochemical and biological functional analyses confirmed the cooperative binding of FOX proteins to the TP53 promoter for the transcriptional activation of TP53. In addition, up-regulation of TP53 was necessary for FOX proteins to exhibit anti-proliferative activity in cancer cells. These analyses reveal the presence of a discrete characteristic within FOX family proteins in which FOX proteins regulate the transcription activity of the p53 tumor suppressor via cooperative binding to the TP53 promoter in alternative dimer configurations.

Update of bioRxiv. 2021 Apr 21;:

CGG repeat expansions in the FMR1 5'UTR cause the neurodegenerative disease Fragile X-associated tremor/ataxia syndrome (FXTAS). These repeats form stable RNA secondary structures that support aberrant translation in the absence of an AUG start codon (RAN translation), producing aggregate-prone peptides that accumulate within intranuclear neuronal inclusions and contribute to neurotoxicity. Here, we show that the most abundant RAN translation product, FMRpolyG, is markedly less toxic when generated from a construct with a non-repetitive alternating codon sequence in place of the CGG repeat. While exploring the mechanism of this differential toxicity, we observed a +1 translational frameshift within the CGG repeat from the arginine to glycine reading frame. Frameshifts occurred within the first few translated repeats and were triggered predominantly by RNA sequence and structural features. Short chimeric R/G peptides form aggregates distinct from those formed by either pure arginine or glycine, and these chimeras induce toxicity in cultured rodent neurons. Together, this work suggests that CGG repeats support translational frameshifting and that chimeric RAN translated peptides may contribute to CGG repeat-associated toxicity in FXTAS and related disorders.

Correct B cell identity at each stage of cellular differentiation during B lymphocyte development is critically dependent on a tightly controlled epigenomic landscape. We previously identified HDAC7 as an essential regulator of early B cell development and its absence leads to a drastic block at the pro-B to pre-B cell transition. More recently, we demonstrated that HDAC7 loss in pro-B-ALL in infants associates with a worse prognosis. Here we delineate the molecular mechanisms by which HDAC7 modulates early B cell development. We find that HDAC7 deficiency drives global chromatin de-condensation, histone marks deposition and deregulates other epigenetic regulators and mobile elements. Specifically, the absence of HDAC7 induces TET2 expression, which promotes DNA 5-hydroxymethylation and chromatin de-condensation. HDAC7 deficiency also results in the aberrant expression of microRNAs and LINE-1 transposable elements. These findings shed light on the mechanisms by which HDAC7 loss or misregulation may lead to B cell-based hematological malignancies.

Alternative pre-mRNA splicing is key to proteome diversity; however, the biological roles of alternative splicing (AS) in signaling pathways remain elusive. Here, we focus on TEA domain transcription factor 1 (TEAD1), a YAP binding factor in the Hippo signaling pathway. Public database analyses showed that expression of YAP-TEAD target genes negatively correlated with the expression of a TEAD1 isoform lacking exon 6 (TEAD1ΔE6) but did not correlate with overall TEAD1 expression. We confirmed that the transcriptional activity and oncogenic properties of the full-length TEAD1 isoform were greater than those of TEAD1ΔE6, with the difference in transcription related to YAP interaction. Furthermore, we showed that RNA-binding Fox-1 homolog 2 (RBFOX2) promoted the inclusion of TEAD1 exon 6 via binding to the conserved GCAUG element in the downstream intron. These results suggest a regulatory mechanism of RBFOX2-mediated TEAD1 AS and provide insight into AS-specific modulation of signaling pathways.

RNA polymerase II (RNAPII) is emerging as an important factor in DNA damage responses, but how it responds to genotoxic stress is not fully understood. We have developed a rapid and sensitive flow cytometry method to study chromatin binding of RNAPII in individual human cells through the cell cycle. Indicating enhanced transcription initiation at early timepoints, levels of RNAPII were increased at 15-30min after UV-induced DNA damage. This was particularly evident for the S5 phosphorylated form of RNAPII (pRNAPII S5), which is typically associated with promoter proximal pausing. Furthermore, degradation of pRNAPII S5 frequently occurs, as its levels on chromatin were strongly enhanced by the proteasome inhibitor MG132 with and without UV. Remarkably, inhibiting pause release with 5,6-dichloro-1-beta-ribo-furanosyl benzimidazole (DRB) further promoted UV-induced degradation of pRNAPII S5, suggesting enhanced initiation may lead to a phenomenon of 'promoter proximal crowding' resulting in premature termination via degradation of RNAPII. Moreover, pRNAPII S2 levels on chromatin were more stable in S phase of the cell cycle 2h after UV, indicating cell cycle specific effects. Altogether our results demonstrate a useful new method and suggest that degradation of promoter proximal RNAPII plays an unanticipated large role both during normal transcription and after UV.

There is an increased interest in the determination of RNA structures in vivo as it is now possible to probe them in a high-throughput manner, e.g. using SHAPE protocols. By now, there exist a range of computational methods that integrate experimental SHAPE-probing evidence into computational RNA secondary structure prediction. The state-of-the-art in this field is currently provided by computational methods that employ the minimum-free energy strategy for prediction RNA secondary structures with SHAPE-probing evidence. These methods, however, rely on the assumption that transcripts in vivo fold into the thermodynamically most stable configuration and ignore evolutionary evidence for conserved RNA structure features. We here present a new computational method, ShapeSorter, that predicts RNA structure features without employing the thermodynamic strategy. Instead, ShapeSorter employs a fully probabilistic framework to identify RNA structure features that are supported by evolutionary and SHAPE-probing evidence. Our method can capture RNA structure heterogeneity, pseudo-knotted RNA structures as well as transient and mutually exclusive RNA structure features. Moreover, it estimates P-values for the predicted RNA structure features which allows for easy filtering and ranking. We investigate the merits of our method in a comprehensive performance benchmarking and conclude that ShapeSorter has a significantly superior performance for predicting base-pairs than the existing state-of-the-art methods.

Noncoding, structured 5'-untranslated regions (5'-UTRs) of bacterial messenger RNAs (mRNAs) can control translation efficiency by forming structures that either recruit or repel the ribosome. Here we exploit a 5'-UTR embedded preQ1-sensing, pseudoknotted translational riboswitch to probe how binding of a small ligand controls recruitment of the bacterial ribosome to the partially overlapping Shine-Dalgarno (SD) sequence. Combining single-molecule fluorescence microscopy with mutational analyses, we find that the stability of 30S ribosomal subunit binding is inversely correlated with the free energy needed to unfold the 5'-UTR during mRNA accommodation into the mRNA binding cleft. Ligand binding to the riboswitch stabilizes the structure to both antagonize 30S recruitment and accelerate 30S dissociation. Proximity of the 5'-UTR and stability of the SD:anti-SD interaction both play important roles in modulating the initial 30S-mRNA interaction. Finally, depletion of small ribosomal subunit protein S1, known to help resolve structured 5'-UTRs, further increases the energetic penalty for mRNA accommodation. The resulting model of rapid standby site exploration followed by gated non-equilibrium unfolding of the 5'-UTR during accommodation provides a mechanistic understanding of how translation efficiency is governed by riboswitches and other dynamic structure motifs embedded upstream of the translation initiation site of bacterial mRNAs.

Enzymatic ligation is a popular method in DNA nanotechnology for structural enforcement. When employed as stability switch for chosen components, ligation can be applied to induce DNA nanostructure reconfiguration. In this study, we investigate the reinforcement effect of ligation on addressable DNA nanostructures assembled entirely from short synthetic strands as the basis of structural reconfiguration. A careful calibration of ligation efficiency is performed on structures with programmable nicks. Systematic investigation using comparative agarose gel electrophoresis enables quantitative assessment of enhanced survivability with ligation treatment on a number of unique structures. The solid ligation performance sets up the foundation for the ligation-based structural reconfiguration. With the capability of switching base pairing status between permanent and transient (ON and OFF) by a simple round of enzymatic treatment, ligation induced reconfiguration can be engineered for DNA nanostructures accordingly.

Accumulation of DNA-RNA hybrids in the form of R-loops can result in replication-transcription conflict that leads to the formation of DNA double strand breaks (DSBs). Using null mutants for the two Caenorhabditis elegans genes encoding for RNaseH1 and RNaseH2, we identify novel effects of R-loop accumulation in the germline. R-loop accumulation leads, as expected, to replication stress, followed by the formation of DSBs. A subset of these DSBs are irreparable. However, unlike irreparable DSBs generated in other systems, which trigger permanent cell cycle arrest, germline irreparable DSBs are propagated to oocytes. Despite DNA damage checkpoint activation in the stem cell niche, the signaling cannot be sustained and nuclei with irreparable DNA damage progress into meiosis. Moreover, unlike other forms of DNA damage that increase germline apoptosis, R-loop-generated DSBs remain undetected by the apoptotic checkpoint. This coincides with attenuation of ATM/ATR signaling in mid-to-late meiotic prophase I. These data altogether indicate that in the germline, DSBs that are generated by R-loops can lead to irreparable DSBs that evade cellular machineries designed for damage recognition. These studies implicate germline R-loops as an especially dangerous driver of germline mutagenesis.

SUMO proteins are important regulators of many key cellular functions in part through their ability to form interactions with other proteins containing SUMO interacting motifs (SIMs). One characteristic feature of all SUMO proteins is the presence of a highly divergent intrinsically disordered region at their N-terminus. In this study, we examine the role of this N-terminal region of SUMO proteins in SUMO-SIM interactions required for the formation of nuclear bodies by the promyelocytic leukemia (PML) protein (PML-NBs). We demonstrate that the N-terminal region of SUMO1 functions in a paralog specific manner as an auto-inhibition domain by blocking its binding to the phosphorylated SIMs of PML and Daxx. Interestingly, we find that this auto-inhibition in SUMO1 is relieved by zinc, and structurally show that zinc stabilizes the complex between SUMO1 and a phospho-mimetic form of the SIM of PML. In addition, we demonstrate that increasing cellular zinc levels enhances PML-NB formation in senescent cells. Taken together, these results provide important insights into a paralog specific function of SUMO1, and suggest that zinc levels could play a crucial role in regulating SUMO1-SIM interactions required for PML-NB formation and function.

Update of bioRxiv. 2021 Jun 15;:

Translation of SARS-CoV-2-encoded mRNAs by the host ribosomes is essential for its propagation. Following infection, the early expressed viral protein NSP1 binds the ribosome, represses translation, and induces mRNA degradation, while the host elicits an anti-viral response. The mechanisms enabling viral mRNAs to escape this multifaceted repression remain obscure. Here we show that expression of NSP1 leads to destabilization of multi-exon cellular mRNAs, while intron-less transcripts, such as viral mRNAs and anti-viral interferon genes, remain relatively stable. We identified a conserved and precisely located cap-proximal RNA element devoid of guanosines that confers resistance to NSP1-mediated translation inhibition. Importantly, the primary sequence rather than the secondary structure is critical for protection. We further show that the genomic 5'UTR of SARS-CoV-2 drives cap-independent translation and promotes expression of NSP1 in an eIF4E-independent and Torin1-resistant manner. Upon expression, NSP1 further enhances cap-independent translation. However, the sub-genomic 5'UTRs are highly sensitive to eIF4E availability, rendering viral propagation partially sensitive to Torin1. We conclude that the combined NSP1-mediated degradation of spliced mRNAs and translation inhibition of single-exon genes, along with the unique features present in the viral 5'UTRs, ensure robust expression of viral mRNAs. These features can be exploited as potential therapeutic targets.

RNase P is a ribonucleoprotein (RNP) that catalyzes removal of the 5' leader from precursor tRNAs in all domains of life. A recent cryo-EM study of Methanocaldococcus jannaschii (Mja) RNase P produced a model at 4.6-Å resolution in a dimeric configuration, with each holoenzyme monomer containing one RNase P RNA (RPR) and one copy each of five RNase P proteins (RPPs; POP5, RPP30, RPP21, RPP29, L7Ae). Here, we used native mass spectrometry (MS), mass photometry (MP), and biochemical experiments that (i) validate the oligomeric state of the Mja RNase P holoenzyme in vitro, (ii) find a different stoichiometry for each holoenzyme monomer with up to two copies of L7Ae, and (iii) assess whether both L7Ae copies are necessary for optimal cleavage activity. By mutating all kink-turns in the RPR, we made the discovery that abolishing the canonical L7Ae-RPR interactions was not detrimental for RNase P assembly and function due to the redundancy provided by protein-protein interactions between L7Ae and other RPPs. Our results provide new insights into the architecture and evolution of RNase P, and highlight the utility of native MS and MP in integrated structural biology approaches that seek to augment the information obtained from low/medium-resolution cryo-EM models.

Translocation of messenger RNA (mRNA) and transfer RNA (tRNA) substrates through the ribosome during protein synthesis, an exemplar of directional molecular movement in biology, entails a complex interplay of conformational, compositional, and chemical changes. The molecular determinants of early translocation steps have been investigated rigorously. However, the elements enabling the ribosome to complete translocation and reset for subsequent protein synthesis reactions remain poorly understood. Here, we have combined molecular simulations with single-molecule fluorescence resonance energy transfer imaging to gain insights into the rate-limiting events of the translocation mechanism. We find that diffusive motions of the ribosomal small subunit head domain to hyper-swivelled positions, governed by universally conserved rRNA, can maneuver the mRNA and tRNAs to their fully translocated positions. Subsequent engagement of peptidyl-tRNA and disengagement of deacyl-tRNA from mRNA, within their respective small subunit binding sites, facilitate the ribosome resetting mechanism after translocation has occurred to enable protein synthesis to resume.

Cell-cell interactions are the fundamental building blocks of tissue organization and multicellular life. We developed Neighbor-seq, a method to identify and annotate the architecture of direct cell-cell interactions and relevant ligand-receptor signaling from the undissociated cell fractions in massively parallel single cell sequencing data. Neighbor-seq accurately identifies microanatomical features of diverse tissue types such as the small intestinal epithelium, terminal respiratory tract, and splenic white pulp. It also captures the differing topologies of cancer-immune-stromal cell communications in pancreatic and skin tumors, which are consistent with the patterns observed in spatial transcriptomic data. Neighbor-seq is fast and scalable. It draws inferences from routine single-cell data and does not require prior knowledge about sample cell-types or multiplets. Neighbor-seq provides a framework to study the organ-level cellular interactome in health and disease, bridging the gap between single-cell and spatial transcriptomics.

Tandem repeats of guanine-rich sequences in RNA often form thermodynamically stable four-stranded RNA structures. Such RNA G-quadruplexes have long been considered to be linked to essential biological processes, yet their physiological significance in cells remains unclear. Here, we report a approach that permits the detection of RNA G-quadruplex structures that modulate protein translation in mammalian cells. The approach combines antibody arrays and RGB-1, a small molecule that selectively stabilizes RNA G-quadruplex structures. Analysis of the protein and mRNA products of 84 cancer-related human genes identified Nectin-4 and CapG as G-quadruplex-controlled genes whose mRNAs harbor non-canonical G-quadruplex structures on their 5'UTR region. Further investigations revealed that the RNA G-quadruplex of CapG exhibits a structural polymorphism, suggesting a possible mechanism that ensures the translation repression in a KCl concentration range of 25-100 mM. The approach described in the present study sets the stage for further discoveries of RNA G-quadruplexes.

The kinetics of DNA hybridization are fundamental to biological processes and DNA-based technologies. However, the precise physical mechanisms that determine why different DNA sequences hybridize at different rates are not well understood. Secondary structure is one predictable factor that influences hybridization rates but is not sufficient on its own to fully explain the observed sequence-dependent variance. In this context, we measured hybridization rates of 43 different DNA sequences that are not predicted to form secondary structure and present a parsimonious physically justified model to quantify our observations. Accounting only for the combinatorics of complementary nucleating interactions and their sequence-dependent stability, the model achieves good correlation with experiment with only two free parameters. Our results indicate that greater repetition of Watson-Crick pairs increases the number of initial states able to proceed to full hybridization, with the stability of those pairings dictating the likelihood of such progression, thus providing new insight into the physical factors underpinning DNA hybridization rates.

The RNA recognition motif (RRM) occurs widely in RNA-binding proteins, but does not always by itself support full binding. For example, it is known that binding of SL1 RNA to the protein U1-70K in the U1 spliceosomal particle is reduced when a region flanking the RRM is truncated. How the RRM flanking regions that together with the RRM make up an 'extended RRM' (eRRM) contribute to complex stability and structural organization is unknown. We study the U1-70K eRRM bound to SL1 RNA by thermal dissociation and laser temperature jump kinetics; long-time molecular dynamics simulations interpret the experiments with atomistic resolution. Truncation of the helix flanking the RRM on its N-terminal side, 'N-helix,' strongly reduces overall binding, which is further weakened under higher salt and temperature conditions. Truncating the disordered region flanking the RRM on the C-terminal side, 'C-IDR', affects the local binding site. Surprisingly, all-atom simulations show that protein truncation enhances base stacking interactions in the binding site and leaves the overall number of hydrogen bonds intact. Instead, the flanking regions of the eRRM act in a distributed fashion via collective interactions with the RNA when external stresses such as temperature or high salt mimicking osmotic imbalance are applied.

mRNA 5' cap recognition by eIF4F is a key element of eukaryotic translational control. Kinetic differences in eIF4F-mRNA interactions have long been proposed to mediate translation-efficiency differences between mRNAs, and recent transcriptome-wide studies have revealed significant heterogeneity in eIF4F engagement with differentially-translated mRNAs. However, detailed kinetic information exists only for eIF4F interactions with short model RNAs. We developed and applied single-molecule fluorescence approaches to directly observe real-time Saccharomyces cerevisiae eIF4F subunit interactions with full-length polyadenylated mRNAs. We found that eIF4E-mRNA association rates linearly anticorrelate with mRNA length. eIF4G-mRNA interaction accelerates eIF4E-mRNA association in proportion to mRNA length, as does an eIF4F-independent activity of eIF4A, though cap-proximal secondary structure still plays an important role in defining the final association rates. eIF4F-mRNA interactions remained dominated by effects of eIF4G, but were modulated to different extents for different mRNAs by the presence of eIF4A and ATP. We also found that eIF4A-catalyzed ATP hydrolysis ejects eIF4E, and likely eIF4E•eIF4G from the mRNA after initial eIF4F•mRNA complex formation, suggesting a mechanism to prepare the mRNA 5' end for ribosome recruitment. Our results support a role for mRNA-specific, factor-driven eIF4F association rates in kinetically controlling translation.

The transcriptional coactivator YAP is emerging as a master regulator of cell growth. In the liver, YAP activity is linked to hepatomegaly, regeneration, dedifferentiation, and aggressive tumor growth. Here we present genomic studies to address how YAP may elicit such profound biological changes in murine models. YAP bound the genome in a TEAD-dependent manner, either at loci constitutively occupied by TEAD or by pioneering enhancers, which comprised a fraction of HNF4a/FOXA-bound embryonic enhancers active during embryonic development but silent in the adult. YAP triggered transcription on promoters by recruiting BRD4, enhancing H3K122 acetylation, and promoting RNApol2 loading and pause-release. YAP also repressed HNF4a target genes by binding to their promoters and enhancers, thus preventing RNApol2 pause-release. YAP activation led to the induction of hepatocyte proliferation, accompanied by tissue remodeling, characterized by polarized macrophages, exhausted T-lymphocytes and dedifferentiation of endothelial cells into proliferative progenitors. Overall, these analyses suggest that YAP is a master regulator of liver function that reshapes the enhancer landscape to control transcription of genes involved in metabolism, proliferation, and inflammation, subverts lineage specification programs by antagonizing HNF4a and modulating the immune infiltrate and the vascular architecture of the liver.

Streptomyces coelicolor (Sc) is a model organism of actinobacteria to study morphological differentiation and production of bioactive metabolites. Sc zinc uptake regulator (Zur) affects both processes by controlling zinc homeostasis. It activates transcription by binding to palindromic Zur-box sequences upstream of -35 elements. Here we deciphered the molecular mechanism by which ScZur interacts with promoter DNA and Sc RNA polymerase (RNAP) by cryo-EM structures and biochemical assays. The ScZur-DNA structures reveal a sequential and cooperative binding of three ScZur dimers surrounding a Zur-box spaced 8 nt upstream from a -35 element. The ScRNAPσHrdB-Zur-DNA structures define protein-protein and protein-DNA interactions involved in the principal housekeeping σHrdB-dependent transcription initiation from a noncanonical promoter with a -10 element lacking the critical adenine residue at position -11 and a TTGCCC -35 element deviating from the canonical TTGACA motif. ScZur interacts with the C-terminal domain of ScRNAP α subunit (αCTD) in a complex structure trapped in an active conformation. Key ScZur-αCTD interfacial residues accounting for ScZur-dependent transcription activation were confirmed by mutational studies. Together, our structural and biochemical results provide a comprehensive model for transcription activation of Zur family regulators.

DNA damage response pathways rely extensively on nuclease activity to process DNA intermediates. Exonuclease 1 (EXO1) is a pleiotropic evolutionary conserved DNA exonuclease involved in various DNA repair pathways, replication, antibody diversification, and meiosis. But, whether EXO1 facilitates these DNA metabolic processes through its enzymatic or scaffolding functions remains unclear. Here, we dissect the contribution of EXO1 enzymatic versus scaffolding activity by comparing Exo1DA/DA mice expressing a proven nuclease-dead mutant form of EXO1 to entirely EXO1-deficient Exo1-/- and EXO1 wild type Exo1+/+ mice. We show that Exo1DA/DA and Exo1-/- mice are compromised in canonical DNA repair processing, suggesting that the EXO1 enzymatic role is important for error-free DNA mismatch and double-strand break repair pathways. However, in non-canonical repair pathways, EXO1 appears to have a more nuanced function. Next-generation sequencing of heavy chain V region in B cells showed the mutation spectra of Exo1DA/DA mice to be intermediate between Exo1+/+ and Exo1-/- mice, suggesting that both catalytic and scaffolding roles of EXO1 are important for somatic hypermutation. Similarly, while overall class switch recombination in Exo1DA/DA and Exo1-/- mice was comparably defective, switch junction analysis suggests that EXO1 might fulfill an additional scaffolding function downstream of class switching. In contrast to Exo1-/- mice that are infertile, meiosis progressed normally in Exo1DA/DA and Exo1+/+ cohorts, indicating that a structural but not the nuclease function of EXO1 is critical for meiosis. However, both Exo1DA/DA and Exo1-/- mice displayed similar mortality and cancer predisposition profiles. Taken together, these data demonstrate that EXO1 has both scaffolding and enzymatic functions in distinct DNA repair processes and suggest a more composite and intricate role for EXO1 in DNA metabolic processes and disease.

The RNA programmed non-specific (trans) nuclease activity of CRISPR-Cas Type V and VI systems has opened a new era in the field of nucleic acid-based detection. Here, we report on the enhancement of trans-cleavage activity of Cas12a enzymes using hairpin DNA sequences as FRET-based reporters. We discover faster rate of trans-cleavage activity of Cas12a due to its improved affinity (Km) for hairpin DNA structures, and provide mechanistic insights of our findings through Molecular Dynamics simulations. Using hairpin DNA probes we significantly enhance FRET-based signal transduction compared to the widely used linear single stranded DNA reporters. Our signal transduction enables faster detection of clinically relevant double stranded DNA targets with improved sensitivity and specificity either in the presence or in the absence of an upstream pre-amplification step.

Autophagy is a catabolic pathway that maintains cellular homeostasis under various stress conditions, including conditions of nutrient deprivation. To elevate autophagic flux to a sufficient level under stress conditions, transcriptional activation of autophagy genes occurs to replenish autophagy components. Thus, the transcriptional and epigenetic control of the genes regulating autophagy is essential for cellular homeostasis. Here, we applied integrated transcriptomic and epigenomic profiling to reveal the roles of plant homeodomain finger protein 20 (PHF20), which is an epigenetic reader possessing methyl binding activity, in controlling the expression of autophagy genes. Phf20 deficiency led to impaired autophagic flux and autophagy gene expression under glucose starvation. Interestingly, the genome-wide characterization of chromatin states by Assay for Transposase-Accessible Chromatin (ATAC)-sequencing revealed that the PHF20-dependent chromatin remodelling occurs in enhancers that are co-occupied by dimethylated lysine 36 on histone H3 (H3K36me2). Importantly, the recognition of H3K36me2 by PHF20 was found to be highly correlated with increased levels of H3K4me1/2 at the enhancer regions. Collectively, these results indicate that PHF20 regulates autophagy genes through enhancer activation via H3K36me2 recognition as an epigenetic reader. Our findings emphasize the importance of nuclear events in the regulation of autophagy.

G-quadruplex (G4) transitions play integral roles in regulating biological functions and can be modified by ligands. However, little is known about G4 transitions. Herein, we reveal distinct pathways of a platinum(II) compound Pt-phen converting parallel-stranded MYC G4 to a hybrid-type structure. Three NMR structures, 1:1 5'-end binding, 1:1 3'-end binding and 2:1 Pt-phen-MYC G4 complexes, were determined by NMR. We find that Pt-phen drives G4 transition at a low ratio. Under physiological 100 mM K+ conditions, a significant stable hydrogen-bonded T:T:A triad is formed at 3'-end of hybrid-type Myc1234, and consequently, Pt-phen first binds the 5'-end to form a 1:1 5'-end binding complex and then disrupts the 3' T:T:A triad and binds 3'-end to form a 2:1 complex with more Pt-phen. Remarkably, the G4 transition pathway is different in 5 mM K+ with Pt-phen first binding the 3'-end and then the 5'-end. 'Edgewise-loop and flanking/ligand/G-tetrad' sandwich structure formation and terminal T:T:A triad stabilization play decisive roles in advancing and altering transition pathways. Our work is the first to elucidate the molecular structures of G4 transitions driven by a small molecule. The ligand-driven G4 transition is a dynamic process that includes a quick G4 transition and multiple complexes formation.

Trypanosomatids cause the neglected tropical diseases, sleeping sickness, Chagas disease and the leishmaniases. Studies on these lethal parasites would be further facilitated by new and improved genetic technologies. Scalable precision editing methods, for example, could be used to improve our understanding of potential mutations associated with drug resistance, a current priority given that several new anti-trypanosomal drugs, with known targets, are currently in clinical development. We report the development of a simple oligo targeting method for rapid and precise editing of priority drug targets in otherwise wild type trypanosomatids. In Trypanosoma brucei, approx. 50-b single-stranded oligodeoxynucleotides were optimal, multiple base edits could be incorporated, and editing efficiency was substantially increased when mismatch repair was suppressed. Resistance-associated edits were introduced in T. brucei cyclin dependent kinase 12 (CRK12, L482F) or cleavage and polyadenylation specificity factor 3 (N232H), in the Trypanosoma cruzi proteasome β5 subunit (G208S), or in Leishmania donovani CRK12 (G572D). We further implemented oligo targeting for site saturation mutagenesis, targeting codon G492 in T. brucei CRK12. This approach, combined with amplicon sequencing for codon variant scoring, revealed fourteen resistance conferring G492 edits encoding six distinct amino acids. The outputs confirm on-target drug activity, reveal a variety of resistance-associated mutations, and facilitate rapid assessment of potential impacts on drug efficacy.

Signal transduction pathways often involve transcription factors that promote activation of defined target gene sets. The transcription factor RBPJ is the central player in Notch signaling and either forms an activator complex with the Notch intracellular domain (NICD) or a repressor complex with corepressors like KYOT2/FHL1. The balance between these two antagonizing RBPJ-complexes depends on the activation state of the Notch receptor regulated by cell-to-cell interaction, ligand binding and proteolytic cleavage events. Here, we depleted RBPJ in mature T-cells lacking active Notch signaling and performed RNA-Seq, ChIP-Seq and ATAC-seq analyses. RBPJ depletion leads to upregulation of many Notch target genes. Ectopic expression of NICD1 activates several Notch target genes and enhances RBPJ occupancy. Based on gene expression changes and RBPJ occupancy we define four different clusters, either RBPJ- and/or Notch-regulated genes. Importantly, we identify early (Hes1 and Hey1) and late Notch-responsive genes (IL2ra). Similarly, to RBPJ depletion, interfering with transcriptional repression by squelching with cofactor KYOT2/FHL1, leads to upregulation of Notch target genes. Taken together, RBPJ is not only an essential part of the Notch co-activator complex but also functions as a repressor in a Notch-independent manner.

Coactivator complexes regulate chromatin accessibility and transcription. SAGA (Spt-Ada-Gcn5 Acetyltransferase) is an evolutionary conserved coactivator complex. The core module scaffolds the entire SAGA complex and adopts a histone octamer-like structure, which consists of six histone-fold domain (HFD)-containing proteins forming three histone-fold (HF) pairs, to which the double HFD-containing SUPT3H adds one HF pair. Spt3, the yeast ortholog of SUPT3H, interacts genetically and biochemically with the TATA binding protein (TBP) and contributes to global RNA polymerase II (Pol II) transcription. Here we demonstrate that (i) SAGA purified from human U2OS or mouse embryonic stem cells (mESC) can assemble without SUPT3H, (ii) SUPT3H is not essential for mESC survival, but required for their growth and self-renewal, and (iii) the loss of SUPT3H from mammalian cells affects the transcription of only a specific subset of genes. Accordingly, in the absence of SUPT3H no major change in TBP accumulation at gene promoters was observed. Thus, SUPT3H is not required for the assembly of SAGA, TBP recruitment, or overall Pol II transcription, but plays a role in mESC growth and self-renewal. Our data further suggest that yeast and mammalian SAGA complexes contribute to transcription regulation by distinct mechanisms.

AlpA positively regulates a programmed cell death pathway linked to the virulence of Pseudomonas aeruginosa by recognizing an AlpA binding element within the promoter, then binding RNA polymerase directly and allowing it to bypass an intrinsic terminator positioned downstream. Here, we report the single-particle cryo-electron microscopy structures of both an AlpA-loading complex and an AlpA-loaded complex. These structures indicate that the C-terminal helix-turn-helix motif of AlpA binds to the AlpA binding element and that the N-terminal segment of AlpA forms a narrow ring inside the RNA exit channel. AlpA was also revealed to render RNAP resistant to termination signals by prohibiting RNA hairpin formation in the RNA exit channel. Structural analysis predicted that AlpA, 21Q, λQ and 82Q share the same mechanism of transcription antitermination.

Erratum in Nucleic Acids Res. 2022 Sep 23;50(17):10200.

Nuclease-inactivated CRISPR/Cas-based (dCas-based) systems have emerged as powerful technologies to synthetically reshape the human epigenome and gene expression. Despite the increasing adoption of these platforms, their relative potencies and mechanistic differences are incompletely characterized, particularly at human enhancer-promoter pairs. Here, we systematically compared the most widely adopted dCas9-based transcriptional activators, as well as an activator consisting of dCas9 fused to the catalytic core of the human CBP protein, at human enhancer-promoter pairs. We find that these platforms display variable relative expression levels in different human cell types and that their transactivation efficacies vary based upon the effector domain, effector recruitment architecture, targeted locus and cell type. We also show that each dCas9-based activator can induce the production of enhancer RNAs (eRNAs) and that this eRNA induction is positively correlated with downstream mRNA expression from a cognate promoter. Additionally, we use dCas9-based activators to demonstrate that an intrinsic transcriptional and epigenetic reciprocity can exist between human enhancers and promoters and that enhancer-mediated tracking and engagement of a downstream promoter can be synthetically driven by targeting dCas9-based transcriptional activators to an enhancer. Collectively, our study provides new insights into the enhancer-mediated control of human gene expression and the use of dCas9-based activators.

RNA-binding proteins (RBPs) have been relatively overlooked in cancer research despite their contribution to virtually every cancer hallmark. Here, we use RNA interactome capture (RIC) to characterize the melanoma RBPome and uncover novel RBPs involved in melanoma progression. Comparison of RIC profiles of a non-tumoral versus a metastatic cell line revealed prevalent changes in RNA-binding capacities that were not associated with changes in RBP levels. Extensive functional validation of a selected group of 24 RBPs using five different in vitro assays unveiled unanticipated roles of RBPs in melanoma malignancy. As proof-of-principle we focused on PDIA6, an ER-lumen chaperone that displayed a novel RNA-binding activity. We show that PDIA6 is involved in metastatic progression, map its RNA-binding domain, and find that RNA binding is required for PDIA6 tumorigenic properties. These results exemplify how RIC technologies can be harnessed to uncover novel vulnerabilities of cancer cells.

Regnase-1 is an evolutionarily conserved endoribonuclease. It degrades diverse mRNAs important for many biological processes including immune homeostasis, development and cancer. There are two competing models of Regnase-1-mediated mRNA silencing. One model postulates that Regnase-1 works together with another RNA-binding protein, Roquin-1, which recruits Regnase-1 to specific mRNAs. The other model proposes that the two proteins function separately. Studying REGE-1, the Caenorhabditis elegans ortholog of Regnase-1, we have uncovered its functional relationship with RLE-1, the nematode counterpart of Roquin-1. While both proteins are essential for mRNA silencing, REGE-1 and RLE-1 appear to associate with target mRNA independently of each other. Thus, although the functional interdependence between REGE-1/Regnase-1 and RLE-1/Roquin-1 is conserved, the underlying mechanisms may display species-specific variation, providing a rare perspective on the evolution of this important post-transcriptional regulatory mechanism.

Genome binning has been essential for characterization of bacteria, archaea, and even eukaryotes from metagenomes. Yet, few approaches exist for viruses. We developed vRhyme, a fast and precise software for construction of viral metagenome-assembled genomes (vMAGs). vRhyme utilizes single- or multi-sample coverage effect size comparisons between scaffolds and employs supervised machine learning to identify nucleotide feature similarities, which are compiled into iterations of weighted networks and refined bins. To refine bins, vRhyme utilizes unique features of viral genomes, namely a protein redundancy scoring mechanism based on the observation that viruses seldom encode redundant genes. Using simulated viromes, we displayed superior performance of vRhyme compared to available binning tools in constructing more complete and uncontaminated vMAGs. When applied to 10,601 viral scaffolds from human skin, vRhyme advanced our understanding of resident viruses, highlighted by identification of a Herelleviridae vMAG comprised of 22 scaffolds, and another vMAG encoding a nitrate reductase metabolic gene, representing near-complete genomes post-binning. vRhyme will enable a convention of binning uncultivated viral genomes and has the potential to transform metagenome-based viral ecology.

Centromeres are the chromosomal loci essential for faithful chromosome segregation during cell division. Although centromeres are transcribed and produce non-coding RNAs (cenRNAs) that affect centromere function, we still lack a mechanistic understanding of how centromere transcription is regulated. Here, using a targeted RNA isoform sequencing approach, we identified the transcriptional landscape at and surrounding all centromeres in budding yeast. Overall, cenRNAs are derived from transcription readthrough of pericentromeric regions but rarely span the entire centromere and are a complex mixture of molecules that are heterogeneous in abundance, orientation, and sequence. While most pericentromeres are transcribed throughout the cell cycle, centromere accessibility to the transcription machinery is restricted to S-phase. This temporal restriction is dependent on Cbf1, a centromere-binding transcription factor, that we demonstrate acts locally as a transcriptional roadblock. Cbf1 deletion leads to an accumulation of cenRNAs at all phases of the cell cycle which correlates with increased chromosome mis-segregation that is partially rescued when the roadblock activity is restored. We propose that a Cbf1-mediated transcriptional roadblock protects yeast centromeres from untimely transcription to ensure genomic stability.

We recently reported that serine-arginine-rich (SR) protein-mediated pre-mRNA structural remodeling generates a pre-mRNA 3D structural scaffold that is stably recognized by the early spliceosomal components. However, the intermediate steps between the free pre-mRNA and the assembled early spliceosome are not yet characterized. By probing the early spliceosomal complexes in vitro and RNA-protein interactions in vivo, we show that the SR proteins bind the pre-mRNAs cooperatively generating a substrate that recruits U1 snRNP and U2AF65 in a splice signal-independent manner. Excess U1 snRNP selectively displaces some of the SR protein molecules from the pre-mRNA generating the substrate for splice signal-specific, sequential recognition by U1 snRNP, U2AF65 and U2AF35. Our work thus identifies a novel function of U1 snRNP in mammalian splicing substrate definition, explains the need for excess U1 snRNP compared to other U snRNPs in vivo, demonstrates how excess SR proteins could inhibit splicing, and provides a conceptual basis to examine if this mechanism of splicing substrate definition is employed by other splicing regulatory proteins.

Although originally described as transcriptional activator, SPI1/PU.1, a major player in haematopoiesis whose alterations are associated with haematological malignancies, has the ability to repress transcription. Here, we investigated the mechanisms underlying gene repression in the erythroid lineage, in which SPI1 exerts an oncogenic function by blocking differentiation. We show that SPI1 represses genes by binding active enhancers that are located in intergenic or gene body regions. HDAC1 acts as a cooperative mediator of SPI1-induced transcriptional repression by deacetylating SPI1-bound enhancers in a subset of genes, including those involved in erythroid differentiation. Enhancer deacetylation impacts on promoter acetylation, chromatin accessibility and RNA pol II occupancy. In addition to the activities of HDAC1, polycomb repressive complex 2 (PRC2) reinforces gene repression by depositing H3K27me3 at promoter sequences when SPI1 is located at enhancer sequences. Moreover, our study identified a synergistic relationship between PRC2 and HDAC1 complexes in mediating the transcriptional repression activity of SPI1, ultimately inducing synergistic adverse effects on leukaemic cell survival. Our results highlight the importance of the mechanism underlying transcriptional repression in leukemic cells, involving complex functional connections between SPI1 and the epigenetic regulators PRC2 and HDAC1.

Sustaining DNA damage response (DDR) signalling via retention of DDR factors at damaged sites is important for transmitting damage-sensing and repair signals. Herein, we found that DNA damage provoked the association of ribosomes with IRES region in lncRNA CTBP1-DT, which overcame the negative effect of upstream open reading frames (uORFs), and elicited the novel microprotein DNA damage-upregulated protein (DDUP) translation via a cap-independent translation mechanism. Activated ATR kinase-mediated phosphorylation of DDUP induced a drastic 'dense-to-loose' conformational change, which sustained the RAD18/RAD51C and RAD18/PCNA complex at damaged sites and initiated RAD51C-mediated homologous recombination and PCNA-mediated post-replication repair mechanisms. Importantly, treatment with ATR inhibitor abolished the effect of DDUP on chromatin retention of RAD51C and PCNA, thereby leading to hypersensitivity of cancer cells to DNA-damaging chemotherapeutics. Taken together, our results uncover a plausible mechanism underlying the DDR sustaining and might represent an attractive therapeutic strategy in improvement of DNA damage-based anticancer therapies.

Integrative and conjugative elements (ICEs) are important drivers of horizontal gene transfer in prokaryotes. They are responsible for antimicrobial resistance spread, a major current health concern. ICEs are initially processed by relaxases that recognize the binding site of oriT sequence and nick at a conserved nic site. The ICESt3/Tn916/ICEBs1 superfamily, which is widespread among Firmicutes, encodes uncanonical relaxases belonging to a recently identified family called MOBT. This family is related to the rolling circle replication initiators of the Rep_trans family. The nic site of these MOBT relaxases is conserved but their DNA binding site is still unknown. Here, we identified the bind site of RelSt3, the MOBT relaxase from ICESt3. Unexpectedly, we found this bind site distantly located from the nic site. We revealed that the binding of the RelSt3 N-terminal HTH domain is required for efficient nicking activity. We also deciphered the role of RelSt3 in the initial and final stages of DNA processing during conjugation. Especially, we demonstrated a strand transfer activity, and the formation of covalent DNA-relaxase intermediate for a MOBT relaxase.

Amino acid substitutions in the exonuclease domain of DNA polymerase ϵ (Polϵ) cause ultramutated tumors. Studies in model organisms suggested pathogenic mechanisms distinct from a simple loss of exonuclease. These mechanisms remain unclear for most recurrent Polϵ mutations. Particularly, the highly prevalent V411L variant remained a long-standing puzzle with no detectable mutator effect in yeast despite the unequivocal association with ultramutation in cancers. Using purified four-subunit yeast Polϵ, we assessed the consequences of substitutions mimicking human V411L, S459F, F367S, L424V and D275V. While the effects on exonuclease activity vary widely, all common cancer-associated variants have increased DNA polymerase activity. Notably, the analog of Polϵ-V411L is among the strongest polymerases, and structural analysis suggests defective polymerase-to-exonuclease site switching. We further show that the V411L analog produces a robust mutator phenotype in strains that lack mismatch repair, indicating a high rate of replication errors. Lastly, unlike wild-type and exonuclease-dead Polϵ, hyperactive variants efficiently synthesize DNA at low dNTP concentrations. We propose that this characteristic could promote cancer cell survival and preferential participation of mutator polymerases in replication during metabolic stress. Our results support the notion that polymerase fitness, rather than low fidelity alone, is an important determinant of variant pathogenicity.

The RAD9-RAD1-HUS1 (9-1-1) clamp forms one half of the DNA damage checkpoint system that signals the presence of substantial regions of single-stranded DNA arising from replication fork collapse or resection of DNA double strand breaks. Loaded at the 5'-recessed end of a dsDNA-ssDNA junction by the RAD17-RFC clamp loader complex, the phosphorylated C-terminal tail of the RAD9 subunit of 9-1-1 engages with the mediator scaffold TOPBP1 which in turn activates the ATR kinase, localised through the interaction of its constitutive partner ATRIP with RPA-coated ssDNA. Using cryogenic electron microscopy (cryoEM) we have determined the structure of a complex of the human RAD17-RFC clamp loader bound to human 9-1-1, engaged with a dsDNA-ssDNA junction. The structure answers the key questions of how RAD17 confers specificity for 9-1-1 over PCNA, and how the clamp loader specifically recognises the recessed 5' DNA end and fixes the orientation of 9-1-1 on the ssDNA.

Gene expression is controlled in part by post-translational modifications of core histones. Methylation of lysine 4 of histone H3 (H3K4), associated with open chromatin and gene transcription, is catalyzed by type 2 lysine methyltransferase complexes that require WDR5, RBBP5, ASH2L and DPY30 as core subunits. Ash2l is essential during embryogenesis and for maintaining adult tissues. To expand on the mechanistic understanding of Ash2l, we generated mouse embryo fibroblasts (MEFs) with conditional Ash2l alleles. Upon loss of Ash2l, methylation of H3K4 and gene expression were downregulated, which correlated with inhibition of proliferation and cell cycle progression. Moreover, we observed induction of senescence concomitant with a set of downregulated signature genes but independent of SASP. Many of the signature genes are FoxM1 responsive. Indeed, exogenous FOXM1 was sufficient to delay senescence. Thus, although the loss of Ash2l in MEFs has broad and complex consequences, a distinct set of downregulated genes promotes senescence.

The removal of RNA primers is essential for mitochondrial DNA (mtDNA) replication. Several nucleases have been implicated in RNA primer removal in human mitochondria, however, no conclusive mechanism has been elucidated. Here, we reconstituted minimal in vitro system capable of processing RNA primers into ligatable DNA ends. We show that human 5'-3' exonuclease, EXOG, plays a fundamental role in removal of the RNA primer. EXOG cleaves short and long RNA-containing flaps but also in cooperation with RNase H1, processes non-flap RNA-containing intermediates. Our data indicate that the enzymatic activity of both enzymes is necessary to process non-flap RNA-containing intermediates and that regardless of the pathway, EXOG-mediated RNA cleavage is necessary prior to ligation by DNA Ligase III. We also show that upregulation of EXOG levels in mitochondria increases ligation efficiency of RNA-containing substrates and discover physical interactions, both in vitro and in cellulo, between RNase H1 and EXOG, Pol γA, Pol γB and Lig III but not FEN1, which we demonstrate to be absent from mitochondria of human lung epithelial cells. Together, using human mtDNA replication enzymes, we reconstitute for the first time RNA primer removal reaction and propose a novel model for RNA primer processing in human mitochondria.

Chromatin insulators are DNA-protein complexes that can prevent the spread of repressive chromatin and block communication between enhancers and promoters to regulate gene expression. In Drosophila, the gypsy chromatin insulator complex consists of three core proteins: CP190, Su(Hw), and Mod(mdg4)67.2. These factors concentrate at nuclear foci termed insulator bodies, and changes in insulator body localization have been observed in mutants defective for insulator function. Here, we identified NURF301/E(bx), a nucleosome remodeling factor, as a novel regulator of gypsy insulator body localization through a high-throughput RNAi imaging screen. NURF301 promotes gypsy-dependent insulator barrier activity and physically interacts with gypsy insulator proteins. Using ChIP-seq, we found that NURF301 co-localizes with insulator proteins genome-wide, and NURF301 promotes chromatin association of Su(Hw) and CP190 at gypsy insulator binding sites. These effects correlate with NURF301-dependent nucleosome repositioning. At the same time, CP190 and Su(Hw) both facilitate recruitment of NURF301 to chromatin. Finally, Oligopaint FISH combined with immunofluorescence revealed that NURF301 promotes 3D contact between insulator bodies and gypsy insulator DNA binding sites, and NURF301 is required for proper nuclear positioning of gypsy binding sites. Our data provide new insights into how a nucleosome remodeling factor and insulator proteins cooperatively contribute to nuclear organization.

Histone modification is a key epigenetic mechanism for regulation of chromatin dynamics and gene expression. Deleted in breast cancer 1 (DBC1) has been shown to act as a negative regulator of epigenetic modifiers and as a co-activator for nuclear receptors and other transcription factors. However, little is known about the role of DBC1 in the regulation of histone modifications and chromatin landscapes. Here, we analyzed genome-wide profiles of active enhancer and promoter marks in colorectal cancer cells and report DBC1 as a critical positive regulator of histone epigenetic writers KMT2D (H3K4 methyltransferase) and p300 (histone acetyltransferase). DBC1 is required for establishing the landscape of active enhancers, for genome-wide chromatin binding and enhancer recruitment of KMT2D and p300, and for gene activation involved in colorectal cancer progression. DBC1 interacts directly with KMT2D and p300, and enhances KMT2D-mediated histone H3K4 methylation (H3K4me1/2/3) and p300-mediated H3 acetylation. Importantly, DBC1 contributes to super-enhancer formation and function by facilitating the recruitment of KMT2D and p300 and by enhancing their functional interaction and cooperative cross-talk. Our results highlight the critical role of DBC1 as a key positive regulator of KMT2D and p300, and provide insights into regulatory mechanisms underlying the interplay between the enhancer epigenomic writers in enhancer activation.

SMARCAL1, ZRANB3 and HLTF are required for the remodeling of replication forks upon stress to promote genome stability. RAD51, along with the RAD51 paralog complex, were also found to have recombination-independent functions in fork reversal, yet the underlying mechanisms remained unclear. Using reconstituted reactions, we build upon previous data to show that SMARCAL1, ZRANB3 and HLTF have unequal biochemical capacities, explaining why they have non-redundant functions. SMARCAL1 uniquely anneals RPA-coated ssDNA, which depends on its direct interaction with RPA, but not on ATP. SMARCAL1, along with ZRANB3, but not HLTF efficiently employ ATPase driven translocase activity to rezip RPA-covered bubbled DNA, which was proposed to mimic elements of fork reversal. In contrast, ZRANB3 and HLTF but not SMARCAL1 are efficient in branch migration that occurs downstream in fork remodeling. We also show that low concentrations of RAD51 and the RAD51 paralog complex, RAD51B-RAD51C-RAD51D-XRCC2 (BCDX2), directly stimulate the motor-driven activities of SMARCAL1 and ZRANB3 but not HLTF, and the interplay is underpinned by physical interactions. Our data provide a possible mechanism explaining previous cellular experiments implicating RAD51 and BCDX2 in fork reversal.

The development of RNA aptamers/fluorophores system is highly desirable for understanding the dynamic molecular biology of RNAs in vivo. Peppers-based imaging systems have been reported and applied for mRNA imaging in living cells. However, the need to insert corresponding RNA aptamer sequences into target RNAs and relatively low fluorescence signal limit its application in endogenous mRNA imaging. Herein, we remolded the original Pepper aptamer and developed a tandem array of inert Pepper (iPepper) fluorescence turn-on system. iPepper allows for efficient and selective imaging of diverse endogenous mRNA species in live cells with minimal agitation of the target mRNAs. We believe iPepper would significantly expand the applications of the aptamer/fluorophore system in endogenous mRNA imaging, and it has the potential to become a powerful tool for real-time studies in living cells and biological processing.

Spatial transcriptomics technologies have recently emerged as a powerful tool for measuring spatially resolved gene expression directly in tissues sections, revealing cell types and their dysfunction in unprecedented detail. However, spatial transcriptomics technologies are limited in their ability to separate transcriptionally similar cell types and can suffer further difficulties identifying cell types in slide regions where transcript capture is low. Here, we describe a conceptually novel methodology that can computationally integrate spatial transcriptomics data with cell-type-informative paired tissue images, obtained from, for example, the reverse side of the same tissue section, to improve inferences of tissue cell type composition in spatial transcriptomics data. The underlying statistical approach is generalizable to any spatial transcriptomics protocol where informative paired tissue images can be obtained. We demonstrate a use case leveraging cell-type-specific immunofluorescence markers obtained on mouse brain tissue sections and a use case for leveraging the output of AI annotated H&E tissue images, which we used to markedly improve the identification of clinically relevant immune cell infiltration in breast cancer tissue. Thus, combining spatial transcriptomics data with paired tissue images has the potential to improve the identification of cell types and hence to improve the applications of spatial transcriptomics that rely on accurate cell type identification.

Interpretation of non-coding genome remains an unsolved challenge in human genetics due to impracticality of exhaustively annotating biochemically active elements in all conditions. Deep learning based computational approaches emerge recently to help interpret non-coding regions. Here, we present LOGO (Language of Genome), a self-attention based contextualized pre-trained language model containing only two self-attention layers with 1 million parameters as a substantially light architecture that applies self-supervision techniques to learn bidirectional representations of the unlabelled human reference genome. LOGO is then fine-tuned for sequence labelling task, and further extended to variant prioritization task via a special input encoding scheme of alternative alleles followed by adding a convolutional module. Experiments show that LOGO achieves 15% absolute improvement for promoter identification and up to 4.5% absolute improvement for enhancer-promoter interaction prediction. LOGO exhibits state-of-the-art multi-task predictive power on thousands of chromatin features with only 3% parameterization benchmarking against the fully supervised model, DeepSEA and 1% parameterization against a recent BERT-based DNA language model. For allelic-effect prediction, locality introduced by one dimensional convolution shows improved sensitivity and specificity for prioritizing non-coding variants associated with human diseases. In addition, we apply LOGO to interpret type 2 diabetes (T2D) GWAS signals and infer underlying regulatory mechanisms. We make a conceptual analogy between natural language and human genome and demonstrate LOGO is an accurate, fast, scalable, and robust framework to interpret non-coding regions for global sequence labeling as well as for variant prioritization at base-resolution.

Update of bioRxiv. 2022 Mar 02;:

Non-CpG PS-ASOs can activate the innate immune system, leading to undesired outcomes. This response can vary-in part-as a function of 2'modifications and sequence. Here we investigated the molecular steps involved in the varied effects of PS-ASOs on the innate immune system. We found that pro-inflammatory PS-ASOs require TLR9 signaling based on the experimental systems used. However, the innate immunity of PS-ASOs does not correlate with their binding affinity with TLR9. Furthermore, the innate immune responses of pro-inflammatory PS-ASOs were reduced by coincubation with non-inflammatory PS-ASOs, suggesting that both pro-inflammatory and non-inflammatory PS-ASOs can interact with TLR9. We show that the kinetics of the PS-ASO innate immune responses can vary, which we speculate may be due to the existence of alternative PS-ASO binding sites on TLR9, leading to full, partial, or no activation of the pathway. In addition, we found that several extracellular proteins, including HMGB1, S100A8 and HRG, enhance the innate immune responses of PS-ASOs. Reduction of the binding affinity by reducing the PS content of PS-ASOs decreased innate immune responses, suggesting that PS-ASO-protein complexes may be sensed by TLR9. These findings thus provide critical information concerning how PS-ASOs can interact with and activate TLR9.