Topoisomerase 3β (TOP3B) and TDRD3 form a dual-activity topoisomerase complex that interacts with FMRP and can change the topology of both DNA and RNA. Here, we investigated the post-transcriptional influence of TOP3B and associated proteins on mRNA translation and turnover. First, we discovered that in human HCT116 colon cancer cells, knock-out (KO) of TOP3B had similar effects on mRNA turnover and translation as did TDRD3-KO, while FMRP-KO resulted in rather distinct effects, indicating that TOP3B had stronger coordination with TDRD3 than FMRP in mRNA regulation. Second, we identified TOP3B-bound mRNAs in HCT116 cells; we found that while TOP3B did not directly influence the stability or translation of most TOP3B target mRNAs, it stabilized a subset of target mRNAs but had a more complex effect on translation-enhancing for some mRNAs whereas reducing for others. Interestingly, a point mutation that specifically disrupted TOP3B catalytic activity only partially recapitulated the effects of TOP3B-KO on mRNA stability and translation, suggesting that the impact of TOP3B on target mRNAs is partly linked to its ability to change topology of mRNAs. Collectively, our data suggest that TOP3B-TDRD3 can regulate mRNA translation and turnover by mechanisms that are dependent and independent of topoisomerase activity.

Break-induced replication (BIR) proceeds via a migrating D-loop for hundreds of kilobases and is highly mutagenic. Previous studies identified long single-stranded (ss) nascent DNA that accumulates during leading strand synthesis to be a target for DNA damage and a primary source of BIR-induced mutagenesis. Here, we describe a new important source of mutagenic ssDNA formed during BIR: the ssDNA template for leading strand BIR synthesis formed during D-loop migration. Specifically, we demonstrate that this D-loop bottom template strand (D-BTS) is susceptible to APOBEC3A (A3A)-induced DNA lesions leading to mutations associated with BIR. Also, we demonstrate that BIR-associated ssDNA promotes an additional type of genetic instability: replication slippage between microhomologies stimulated by inverted DNA repeats. Based on our results we propose that these events are stimulated by both known sources of ssDNA formed during BIR, nascent DNA formed by leading strand synthesis, and the D-BTS that we describe here. Together we report a new source of mutagenesis during BIR that may also be shared by other homologous recombination pathways driven by D-loop repair synthesis.

The Caenorhabditis elegans model has greatly contributed to the understanding of the role of G-quadruplexes in genomic instability. The GGCTTA repeats of the C. elegans telomeres resemble the GGGTTA repeats of the human telomeres. However, the comparison of telomeric sequences (Homo sapiens, Tetrahymena, Oxytricha, Bombyx mori and Giardia) revealed that small changes in these repeats can drastically change the topology of the folded G-quadruplex. In the present work we determined the structure adopted by the C. elegans telomeric sequence d[GG(CTTAGG)3]. The investigated C. elegans telomeric sequence is shown to fold into an intramolecular two G-tetrads basket type G-quadruplex structure that includes a C-T base pair in the diagonal loop. This work sheds light on the telomeric structure of the widely used C. elegans animal model.

Preclinical mechanistic studies have pointed towards RNA interference-mediated off-target effects as a major driver of hepatotoxicity for GalNAc-siRNA conjugates. Here, we demonstrate that a single glycol nucleic acid or 2'-5'-RNA modification can substantially reduce small interfering RNA (siRNA) seed-mediated binding to off-target transcripts while maintaining on-target activity. In siRNAs with established hepatotoxicity driven by off-target effects, these novel designs with seed-pairing destabilization, termed enhanced stabilization chemistry plus (ESC+), demonstrated a substantially improved therapeutic window in rats. In contrast, siRNAs thermally destabilized to a similar extent by the incorporation of multiple DNA nucleotides in the seed region showed little to no improvement in rat safety suggesting that factors in addition to global thermodynamics play a role in off-target mitigation. We utilized the ESC+ strategy to improve the safety of ALN-HBV, which exhibited dose-dependent, transient and asymptomatic alanine aminotransferase elevations in healthy volunteers. The redesigned ALN-HBV02 (VIR-2218) showed improved specificity with comparable on-target activity and the program was reintroduced into clinical development.

Nitric oxide (NO) is a key player in numerous physiological processes. Excessive NO induces DNA damage, but how plants respond to this damage remains unclear. We screened and identified an Arabidopsis NO hypersensitive mutant and found it to be allelic to TEBICHI/POLQ, encoding DNA polymerase θ. The teb mutant plants were preferentially sensitive to NO- and its derivative peroxynitrite-induced DNA damage and subsequent double-strand breaks (DSBs). Inactivation of TEB caused the accumulation of spontaneous DSBs largely attributed to endogenous NO and was synergistic to DSB repair pathway mutations with respect to growth. These effects were manifested in the presence of NO-inducing agents and relieved by NO scavengers. NO induced G2/M cell cycle arrest in the teb mutant, indicative of stalled replication forks. Genetic analyses indicate that Polθ is required for translesion DNA synthesis across NO-induced lesions, but not oxidation-induced lesions. Whole-genome sequencing revealed that Polθ bypasses NO-induced base adducts in an error-free manner and generates mutations characteristic of Polθ-mediated end joining. Our experimental data collectively suggests that Polθ plays dual roles in protecting plants from NO-induced DNA damage. Since Polθ is conserved in higher eukaryotes, mammalian Polθ may also be required for balancing NO physiological signaling and genotoxicity.

The E3 ubiquitin ligase TRIM25 is a key factor in the innate immune response to RNA viruses. TRIM25 has been shown to play a role in the retinoic-acid-inducible gene-1 (RIG-I) pathway, which triggers expression of type 1 interferons upon viral infection. We and others have shown that TRIM25 is an RNA-binding protein; however, the role of TRIM25 RNA-binding in the innate immune response to RNA viruses is unclear. Here, we demonstrate that influenza A virus (IAV A/PR/8/34_NS1(R38A/K41A)) infection is inhibited by TRIM25. Surprisingly, previously identified RNA-binding deficient mutant TRIM25ΔRBD and E3 ubiquitin ligase mutant TRIM25ΔRING, which lack E3 ubiquitin ligase activity, still inhibited IAV replication. Furthermore, we show that in human-derived cultured cells, activation of the RIG-I/interferon type 1 pathway mediated by either an IAV-derived 5'-triphosphate RNA or by IAV itself does not require TRIM25 activity. Additionally, we present new evidence that instead of TRIM25 directly inhibiting IAV transcription it binds and destabilizes IAV mRNAs. Finally, we show that direct tethering of TRIM25 to RNA is sufficient to downregulate the targeted RNA. In summary, our results uncover a potential mechanism that TRIM25 uses to inhibit IAV infection and regulate RNA metabolism.

Although reverse-transcriptase (RT) enzymes are critical reagents for research and biotechnology, their mechanical properties are not well understood. In particular, we know little about their relative speed and response to structural obstacles in the template. Commercial retroviral RTs stop at many positions along mixed sequence templates, resulting in truncated cDNA products that complicate downstream analysis. By contrast, group II intron-encoded RTs appear to copy long RNAs with high processivity and minimal stops. However, their speed, consistency and pausing behavior have not been explored. Here, we analyze RT velocity as the enzyme moves through heterogeneous sequences and structures that are embedded within a long noncoding RNA transcript. We observe that heterogeneities in the template are highly disruptive to primer extension by retroviral RTs. However, sequence composition and template structure have negligible effects on behavior of group II intron RTs, such as MarathonRT (MRT). Indeed, MRT copies long RNAs in a single pass, and displays synchronized primer extension at a constant speed of 25 nt/sec. In addition, it passes through stable RNA structural motifs without perturbation of velocity. Taken together, the results demonstrate that consistent, robust translocative behavior is a hallmark of group II intron-encoded RTs, some of which operate at high velocity.

Designer RNA scaffolds constitute a promising tool for synthetic biology, as they can be genetically expressed to perform specific functions in vivo such as scaffolding enzymatic cascades and regulating gene expression through CRISPR-dCas9 applications. RNA origami is a recently developed RNA design approach that allows construction of large RNA nanostructures that can position aptamer motifs to spatially organize other molecules, including proteins. However, it is still not fully understood how positioning multiple aptamers on a scaffold and the orientation of a scaffold affects functional properties. Here, we investigate fusions of single-guide RNAs and RNA origami scaffolds (termed sgRNAO) capable of recruiting activating domains for control of gene expression in yeast. Using MS2 and PP7 as orthogonal protein-binding aptamers, we observe a gradual increase in transcriptional activation for up to four aptamers. We demonstrate that different aptamer positions on a scaffold and scaffold orientation affect transcriptional activation. Finally, sgRNAOs are used to regulate expression of enzymes of the violacein biosynthesis pathway to control metabolic flux. The integration of RNA origami nanostructures at promoter sites achieved here, can in the future be expanded by the addition of functional motifs such as riboswitches, ribozymes and sensor elements to allow for complex gene regulation.

Erratum in Nucleic Acids Res. 2022 Nov 28;50(21):12598.

Dimension reduction and (spatial) clustering is usually performed sequentially; however, the low-dimensional embeddings estimated in the dimension-reduction step may not be relevant to the class labels inferred in the clustering step. We therefore developed a computation method, Dimension-Reduction Spatial-Clustering (DR-SC), that can simultaneously perform dimension reduction and (spatial) clustering within a unified framework. Joint analysis by DR-SC produces accurate (spatial) clustering results and ensures the effective extraction of biologically informative low-dimensional features. DR-SC is applicable to spatial clustering in spatial transcriptomics that characterizes the spatial organization of the tissue by segregating it into multiple tissue structures. Here, DR-SC relies on a latent hidden Markov random field model to encourage the spatial smoothness of the detected spatial cluster boundaries. Underlying DR-SC is an efficient expectation-maximization algorithm based on an iterative conditional mode. As such, DR-SC is scalable to large sample sizes and can optimize the spatial smoothness parameter in a data-driven manner. With comprehensive simulations and real data applications, we show that DR-SC outperforms existing clustering and spatial clustering methods: it extracts more biologically relevant features than conventional dimension reduction methods, improves clustering performance, and offers improved trajectory inference and visualization for downstream trajectory inference analyses.

The study and manipulation of T cell receptors (TCRs) is central to multiple fields across basic and translational immunology research. Produced by V(D)J recombination, TCRs are often only recorded in the literature and data repositories as a combination of their V and J gene symbols, plus their hypervariable CDR3 amino acid sequence. However, numerous applications require full-length coding nucleotide sequences. Here we present Stitchr, a software tool developed to specifically address this limitation. Given minimal V/J/CDR3 information, Stitchr produces complete coding sequences representing a fully spliced TCR cDNA. Due to its modular design, Stitchr can be used for TCR engineering using either published germline or novel/modified variable and constant region sequences. Sequences produced by Stitchr were validated by synthesizing and transducing TCR sequences into Jurkat cells, recapitulating the expected antigen specificity of the parental TCR. Using a companion script, Thimble, we demonstrate that Stitchr can process a million TCRs in under ten minutes using a standard desktop personal computer. By systematizing the production and modification of TCR sequences, we propose that Stitchr will increase the speed, repeatability, and reproducibility of TCR research. Stitchr is available on GitHub.

Protein dynamics involving higher-energy sparsely populated conformational substates are frequently critical for protein function. This study describes the dynamics of the homodimer (p50)2 of the p50 Rel homology region (RHR) of the transcription factor NF-κB, using 13C relaxation dispersion experiments with specifically (13C, 1H)-labeled methyl groups of Ile (δ), Leu and Val. Free (p50)2 is highly dynamic in solution, showing μs-ms relaxation dispersion consistent with exchange between the ground state and higher energy substates. These fluctuations propagate from the DNA-binding loops through the core of the domain. The motions are damped in the presence of κB DNA, but the NMR spectra of the DNA complexes reveal multiple local conformations of the p50 RHR homodimer bound to certain κB DNA sequences. Varying the length and sequence of κB DNA revealed two factors that promote a single bound conformation for the complex: the length of the κB site in the duplex and a symmetrical sequence of guanine nucleotides at both ends of the recognition motif. The dynamic nature of the DNA-binding loops, together with the multiple bound conformations of p50 RHR with certain κB sites, is consistent with variations in the transcriptional activity of the p50 homodimer with different κB sequences.

Compared with the dual binding components in a binary interaction, the third component of a ternary interaction often serves as modulator or regulator in biochemical processes. Here, we presented a programmable ternary interaction strategy based on the natural DNA triplex structure. With the DNA triplex-based ternary interaction, we have successfully demonstrated controllable hierarchical assemblies from nanometer scale synthetic DNA nanostructure units to micrometer scale live bacteria. A selective signaling system responsive to orthogonal nucleic acid signals via ternary interaction was also demonstrated. This assembly method could further enrich the diversified design schemes of DNA nanotechnology.

tRNA fragmentation is an evolutionarily conserved molecular phenomenon. tRNA-derived small RNAs (tsRNAs) have been associated with many cellular processes, including improved survival during stress conditions. Here, we have revisited accepted experimental paradigms for modeling oxidative stress resulting in tRNA fragmentation. Various cell culture models were exposed to oxidative stressors followed by determining cell viability, the production of specific tsRNAs and stress granule formation. These experiments revealed that exposure to stress parameters commonly used to induce tRNA fragmentation negatively affected cell viability after stress removal. Quantification of specific tsRNA species in cells responding to experimental stress and in cells that were transfected with synthetic tsRNAs indicated that neither physiological nor non-physiological copy numbers of tsRNAs induced the formation of stress granules. Furthermore, the increased presence of tsRNA species in culture medium collected from stressed cells indicated that cells suffering from experimental stress exposure gave rise to stable extracellular tsRNAs. These findings suggest a need to modify current experimental stress paradigms in order to allow separating the function of tRNA fragmentation during the acute stress response from tRNA fragmentation as a consequence of ongoing cell death, which will have major implications for the current perception of the biological function of stress-induced tsRNAs.

Translational control of maternal mRNAs generates spatial and temporal patterns of protein expression necessary to begin animal development. Translational repression of unlocalized nanos (nos) mRNA in late-stage Drosophila oocytes by the hnRNP F/H homolog, Glorund (Glo), is important for embryonic body patterning. While previous work has suggested that repression occurs at both the translation initiation and elongation phases, the molecular mechanism by which Glo regulates nos translation remains elusive. Here, we have identified the Drosophila fragile X mental retardation protein, dFMRP, as a Glo interaction partner with links to the translational machinery. Using an oocyte-based in vitro translation system, we confirmed that Glo regulates both initiation and elongation of a nos translational reporter and showed that dFMRP specifically represses translation elongation and promotes ribosome stalling. Furthermore, we combined mutational analysis and in vivo and in vitro binding assays to show that Glo's qRRM2 domain specifically and directly interacts with dFMRP. Our findings suggest that Glo regulates nos translation elongation by recruiting dFMRP and that Glo's RNA-binding domains can also function as protein-protein interaction interfaces critical for its regulatory functions. Additionally, they reveal a mechanism for targeting dFMRP to specific transcripts.

Gliomas are one of the most common and lethal brain tumors among adults. One process that contributes to glioma progression and recurrence is the epithelial to mesenchymal transition (EMT). EMT is regulated by a set of defined transcription factors which tightly regulate this process, among them is the basic helix-loop-helix family member, TWIST1. Here we show that TWIST1 is methylated on lysine-33 at chromatin by SETD6, a methyltransferase with expression levels correlating with poor survival in glioma patients. RNA-seq analysis in U251 glioma cells suggested that both SETD6 and TWIST1 regulate cell adhesion and migration processes. We further show that TWIST1 methylation attenuates the expression of the long-non-coding RNA, LINC-PINT, thereby promoting EMT in glioma. Mechanistically, TWIST1 methylation represses the transcription of LINC-PINT by increasing the occupancy of EZH2 and the catalysis of the repressive H3K27me3 mark at the LINC-PINT locus. Under un-methylated conditions, TWIST1 dissociates from the LINC-PINT locus, allowing the expression of LINC-PINT which leads to increased cell adhesion and decreased cell migration. Together, our findings unravel a new mechanistic dimension for selective expression of LINC-PINT mediated by TWIST1 methylation.

The standard analysis pipeline for single-cell RNA-seq data consists of sequential steps initiated by clustering the cells. An innate limitation of this pipeline is that an imperfect clustering result can irreversibly affect the succeeding steps. For example, there can be cell types not well distinguished by clustering because they largely share the global structure, such as the anterior primitive streak and mid primitive streak cells. If one searches differentially expressed genes (DEGs) solely based on clustering, marker genes for distinguishing these types will be missed. Moreover, clustering depends on many parameters and can often be subjective to manual decisions. To overcome these limitations, we propose MarcoPolo, a method that identifies informative DEGs independently of prior clustering. MarcoPolo sorts out genes by evaluating if the distributions are bimodal, if similar expression patterns are observed in other genes, and if the expressing cells are proximal in a low-dimensional space. Using real datasets with FACS-purified cell labels, we demonstrate that MarcoPolo recovers marker genes better than competing methods. Notably, MarcoPolo finds key genes that can distinguish cell types that are not distinguishable by the standard clustering. MarcoPolo is built in a convenient software package that provides analysis results in an HTML file.

G-quadruplex structure (G4) is a type of DNA secondary structure that widely exists in the genomes of many organisms. G4s are believed to participate in multiple biological processes. Acyl-CoA binding protein (ACBP), a ubiquitously expressed and highly conserved protein in eukaryotic cells, plays important roles in lipid metabolism by transporting and protecting acyl-CoA esters. Here, we report the functional identification of a G4 in the promoter of the ACBP gene in silkworm and human cancer cells. We found that G4 exists as a conserved element in the promoters of ACBP genes in invertebrates and vertebrates. The BmACBP G4 bound with G4-binding protein LARK regulated BmACBP transcription, which was blocked by the G4 stabilizer pyridostatin (PDS) and G4 antisense oligonucleotides. PDS treatment with fifth instar silkworm larvae decreased the BmACBP expression and triacylglycerides (TAG) level, resulting in reductions in fat body mass, body size and weight and growth and metamorphic rates. PDS treatment and knocking out of the HsACBP G4 in human hepatic adenocarcinoma HepG2 cells inhibited the expression of HsACBP and decreased the TAG level and cell proliferation. Altogether, our findings suggest that G4 of the ACBP genes is involved in regulation of lipid metabolism processes in invertebrates and vertebrates.

Spermatogenesis is precisely controlled by sophisticated gene expression programs and is driven by epigenetic reprogramming, including histone modification alterations and histone-to-protamine transition. Nuclear receptor binding SET domain protein 2 (Nsd2) is the predominant histone methyltransferase catalyzing H3K36me2 and its role in male germ cell development remains elusive. Here, we report that NSD2 protein is abundant in spermatogenic cells. Conditional loss of Nsd2 in postnatal germ cells impaired fertility owing to apoptosis of spermatocytes and aberrant spermiogenesis. Nsd2 deficiency results in dysregulation of thousands of genes and remarkable reduction of both H3K36me2 and H3K36me3 in spermatogenic cells, with H3K36me2 occupancy correlating positively with expression of germline genes. Nsd2 deficiency leads to H4K16ac elevation in spermatogenic cells, probably through interaction between NSD2 and PSMA8, which regulates acetylated histone degradation. We further reveal that Nsd2 deficiency impairs EP300-induced H4K5/8ac, recognized by BRDT to mediate the eviction of histones. Accordingly, histones are largely retained in Nsd2-deficient spermatozoa. In addition, Nsd2 deficiency enhances expression of protamine genes, leading to increased protamine proteins in Nsd2-deficient spermatozoa. Our findings thus reveal a previously unappreciated role of the Nsd2-dependent chromatin remodeling during spermatogenesis and provide clues to the molecular mechanisms in epigenetic abnormalities impacting male reproductive health.

A G-rich sequence was designed to allow folding into either a stable parallel or hybrid-type topology. With the parent sequence featuring coexisting species, various related sequences with single and double mutations and with a shortened central propeller loop affected the topological equilibrium. Two simple modifications, likewise introduced separately to all sequences, were employed to lock folds into one of the topologies without noticeable structural alterations. The unique combination of sequence mutations, high-resolution NMR structural information, and the thermodynamic stability for both topological competitors identified critical loop residue interactions. In contrast to first loop residues, which are mostly disordered and exposed to solvent in both propeller and lateral loops bridging a narrow groove, the last loop residue in a lateral three-nucleotide loop is engaged in stabilizing stacking interactions. The propensity of single-nucleotide loops to favor all-parallel topologies by enforcing a propeller-like conformation of an additional longer loop is shown to result from their preference in linking two outer tetrads of the same tetrad polarity. Taken together, the present studies contribute to a better structural and thermodynamic understanding of delicate loop interactions in genomic and artificially designed quadruplexes, e.g. when employed as therapeutics or in other biotechnological applications.

Two-component systems (TCS) and small RNAs (sRNA) are widespread regulators that participate in the response and the adaptation of bacteria to their environments. TCSs and sRNAs mostly act at the transcriptional and post-transcriptional levels, respectively, and can be found integrated in regulatory circuits, where TCSs control sRNAs transcription and/or sRNAs post-transcriptionally regulate TCSs synthesis. In response to nitrate and nitrite, the paralogous NarQ-NarP and NarX-NarL TCSs regulate the expression of genes involved in anaerobic respiration of these alternative electron acceptors to oxygen. In addition to the previously reported repression of NarP synthesis by the SdsN137 sRNA, we show here that RprA, another Hfq-dependent sRNA, also negatively controls narP. Interestingly, the repression of narP by RprA actually relies on two independent mechanisms of control. The first is via the direct pairing of the central region of RprA to the narP translation initiation region and presumably occurs at the translation initiation level. In contrast, the second requires only the very 5' end of the narP mRNA, which is targeted, most likely indirectly, by the full-length or the shorter, processed, form of RprA. In addition, our results raise the possibility of a direct role of Hfq in narP control, further illustrating the diversity of post-transcriptional regulation mechanisms in the synthesis of TCSs.

Ribonucleotides can be incorporated into DNA during replication by the replicative DNA polymerases. These aberrant DNA subunits are efficiently recognized and removed by Ribonucleotide Excision Repair, which is initiated by the heterotrimeric enzyme RNase H2. While RNase H2 is essential in higher eukaryotes, the yeast Saccharomyces cerevisiae can survive without RNase H2 enzyme, although the genome undergoes mutation, recombination and other genome instability events at an increased rate. Although RNase H2 can be considered as a protector of the genome from the deleterious events that can ensue from recognition and removal of embedded ribonucleotides, under conditions of high ribonucleotide incorporation and retention in the genome in a RNase H2-negative strain, sudden introduction of active RNase H2 causes massive DNA breaks and genome instability in a condition which we term 'ribodysgenesis'. The DNA breaks and genome instability arise solely from RNase H2 cleavage directed to the ribonucleotide-containing genome. Survivors of ribodysgenesis have massive loss of heterozygosity events stemming from recombinogenic lesions on the ribonucleotide-containing DNA, with increases of over 1000X from wild-type. DNA breaks are produced over one to two divisions and subsequently cells adapt to RNase H2 and ribonucleotides in the genome and grow with normal levels of genome instability.

The robustness and sensitivity of gene networks to environmental changes is critical for cell survival. How gene networks produce specific, chronologically ordered responses to genome-wide perturbations, while robustly maintaining homeostasis, remains an open question. We analysed if short- and mid-term genome-wide responses to shifts in RNA polymerase (RNAP) concentration are influenced by the known topology and logic of the transcription factor network (TFN) of Escherichia coli. We found that, at the gene cohort level, the magnitude of the single-gene, mid-term transcriptional responses to changes in RNAP concentration can be explained by the absolute difference between the gene's numbers of activating and repressing input transcription factors (TFs). Interestingly, this difference is strongly positively correlated with the number of input TFs of the gene. Meanwhile, short-term responses showed only weak influence from the TFN. Our results suggest that the global topological traits of the TFN of E. coli shape which gene cohorts respond to genome-wide stresses.

MicroRNAs (miRNAs) bind to complementary target RNAs and regulate their gene expression post-transcriptionally. These non-coding regulatory RNAs become functional after loading into Argonaute (AGO) proteins to form the effector complexes. Humans have four AGO proteins, AGO1, AGO2, AGO3 and AGO4, which share a high sequence identity. Since most miRNAs are found across the four AGOs, it has been thought that they work redundantly, and AGO2 has been heavily studied as the exemplified human paralog. Nevertheless, an increasing number of studies have found that the other paralogs play unique roles in various biological processes and diseases. In the last decade, the structural study of the four AGOs has provided the field with solid structural bases. This review exploits the completed structural catalog to describe common features and differences in target specificity across the four AGOs.

The non-structural protein 3 helicase (NS3h) is a multifunctional protein that is critical in RNA replication and other stages in the flavivirus life cycle. NS3h uses energy from ATP hydrolysis to translocate along single stranded nucleic acid and to unwind double stranded RNA. Here we present a detailed mechanistic analysis of the product release stage in the catalytic cycle of the dengue virus (DENV) NS3h. This study is based on a combined experimental and computational approach of product-inhibition studies and free energy calculations. Our results support a model in which the catalytic cycle of ATP hydrolysis proceeds through an ordered sequential mechanism that includes a ternary complex intermediate (NS3h-Pi-ADP), which evolves releasing the first product, phosphate (Pi), and subsequently ADP. Our results indicate that in the product release stage of the DENV NS3h a novel open-loop conformation plays an important role that may be conserved in NS3 proteins of other flaviviruses as well.

Long noncoding RNAs (lncRNAs) and microRNAs (miRNAs) modulate gene expression programs in physiology and disease. Here, we report a noncoding RNA regulatory network that modulates myoblast fusion into multinucleated myotubes, a process that occurs during muscle development and muscle regeneration after injury. In early stages of human myogenesis, the levels of lncRNA OIP5-AS1 increased, while the levels of miR-7 decreased. Moreover, OIP5-AS1 bound and induced miR-7 decay via target RNA-directed miRNA decay; accordingly, loss of OIP5-AS1 attenuated, while antagonizing miR-7 accelerated, myotube formation. We found that the OIP5-AS1-mediated miR-7 degradation promoted myoblast fusion, as it derepressed the miR-7 target MYMX mRNA, which encodes the fusogenic protein myomixer (MYMX). Remarkably, an oligonucleotide site blocker interfered with the OIP5-AS1-directed miR-7 degradation, allowing miR-7 to accumulate, lowering MYMX production and suppressing myotube formation. These results highlight a mechanism whereby lncRNA OIP5-AS1-mediated miR-7 decay promotes myotube formation by stimulating a myogenic fusion program.

Homologs of the mutagenic Escherichia coli DNA polymerase V (pol V) are encoded by numerous pathogens and mobile elements. We have used Rum pol (RumA'2B), from the integrative conjugative element (ICE), R391, as a model mobile element-encoded polymerase (MEPol). The highly mutagenic Rum pol is transferred horizontally into a variety of recipient cells, including many pathogens. Moving between species, it is unclear if Rum pol can function on its own or requires activation by host factors. Here, we show that Rum pol biochemical activity requires the formation of a physical mutasomal complex, Rum Mut, containing RumA'2B-RecA-ATP, with RecA being donated by each recipient bacteria. For R391, Rum Mut specific activities in vitro and mutagenesis rates in vivo depend on the phylogenetic distance of host-cell RecA from E. coli RecA. Rum pol is a highly conserved and effective mobile catalyst of rapid evolution, with the potential to generate a broad mutational landscape that could serve to ensure bacterial adaptation in antibiotic-rich environments leading to the establishment of antibiotic resistance.

In Escherichia coli, transcription-translation coupling is mediated by NusG. Although chloroplasts are descendants of endosymbiotic prokaryotes, the mechanism underlying this coupling in chloroplasts remains unclear. Here, we report transcription-translation coupling through AtNusG in chloroplasts. AtNusG is localized in chloroplast nucleoids and is closely associated with the chloroplast PEP complex by interacting with its essential component PAP9. It also comigrates with chloroplast ribosomes and interacts with their two components PRPS5 (uS5c) and PRPS10 (uS10c). These data suggest that the transcription and translation machineries are coupled in chloroplasts. In the atnusg mutant, the accumulation of chloroplast-encoded photosynthetic gene transcripts, such as psbA, psbB, psbC and psbD, was not obviously changed, but that of their proteins was clearly decreased. Chloroplast polysomic analysis indicated that the decrease in these proteins was due to the reduced efficiency of their translation in this mutant, leading to reduced photosynthetic efficiency and enhanced sensitivity to cold stress. These data indicate that AtNusG-mediated coupling between transcription and translation in chloroplasts ensures the rapid establishment of photosynthetic capacity for plant growth and the response to environmental changes. Therefore, our study reveals a conserved mechanism of transcription-translation coupling between chloroplasts and E. coli, which perhaps represents a regulatory mechanism of chloroplast gene expression. This study provides insights into the underlying mechanisms of chloroplast gene expression in higher plants.

Bacteria and bacteriophages co-evolve in a constant arms race, wherein one tries and finds newer ways to overcome the other. Phage resistance poses a great threat to the development of phage therapy. Hence, it is both essential and important to understand the mechanism of phage resistance in bacteria. First identified in Mycobacterium smegmatis, the gene mpr, upon overexpression, confers resistance against D29 mycobacteriophage. Presently, the mechanism behind phage resistance by mpr is poorly understood. Here we show that Mpr is a membrane-bound DNA exonuclease, which digests DNA in a non-specific manner independent of the sequence, and shares no sequence or structural similarity with any known nuclease. Exonuclease activity of mpr provides resistance against phage infection, but the role of mpr may very well go beyond just phage resistance. Our experiments show that mpr plays a crucial role in the appearance of mutant colonies (phage resistant strains). However, the molecular mechanism behind the emergence of these mutant/resistant colonies is yet to be understood. Nevertheless, it appears that mpr is involved in the survival and evolution of M. smegmatis against phage. A similar mechanism may be present in other organisms, which requires further exploration.

In nucleotide excision repair (NER), the xeroderma pigmentosum D helicase (XPD) scans DNA searching for bulky lesions, stalls when encountering such damage to verify its presence, and allows repair to proceed. Structural studies have shown XPD bound to its single-stranded DNA substrate, but molecular and dynamic characterization of how XPD translocates on undamaged DNA and how it stalls to verify lesions remains poorly understood. Here, we have performed extensive all-atom MD simulations of human XPD bound to undamaged and damaged ssDNA, containing a mutagenic pyrimidine (6-4) pyrimidone UV photoproduct (6-4PP), near the XPD pore entrance. We characterize how XPD responds to the presence of the DNA lesion, delineating the atomistic-scale mechanism that it utilizes to discriminate between damaged and undamaged nucleotides. We identify key amino acid residues, including FeS residues R112, R196, H135, K128, Arch residues E377 and R380, and ATPase lobe 1 residues 215-221, that are involved in damage verification and show how movements of Arch and ATPase lobe 1 domains relative to the FeS domain modulate these interactions. These structural and dynamic molecular depictions of XPD helicase activity with unmodified DNA and its inhibition by the lesion elucidate how the lesion is verified by inducing XPD stalling.

Antiviral innate immunity represents the first defense against invading viruses and is key to control viral infections, including SARS-CoV-2. Body temperature is an omnipresent variable but was neglected when addressing host defense mechanisms and susceptibility to SARS-CoV-2 infection. Here, we show that increasing temperature in a 1.5°C window, between 36.5 and 38°C, strongly increases the expression of genes in two branches of antiviral immunity, nitric oxide production and type I interferon response. We show that alternative splicing coupled to nonsense-mediated decay decreases STAT2 expression in colder conditions and suggest that increased STAT2 expression at elevated temperature induces the expression of diverse antiviral genes and SARS-CoV-2 restriction factors. This cascade is activated in a remarkably narrow temperature range below febrile temperature, which reflects individual, circadian and age-dependent variation. We suggest that decreased body temperature with aging contributes to reduced expression of antiviral genes in older individuals. Using cell culture and in vivo models, we show that higher body temperature correlates with reduced SARS-CoV-2 replication, which may affect the different vulnerability of children versus seniors toward severe SARS-CoV-2 infection. Altogether, our data connect body temperature and pre-mRNA processing to provide new mechanistic insight into the regulation of antiviral innate immunity.

The retrovirus human immunodeficiency virus-1 (HIV-1) is the causative agent of AIDS. Although treatment of HIV/AIDS with antiretroviral therapy provides suppression of viremia, latent reservoirs of integrated proviruses preclude cure by current antiviral treatments. Understanding the mechanisms of host-viral interactions may elucidate new treatment strategies. Here, we performed a CRISPR/Cas9 transcriptional activation screen using a high-complexity, genome-wide sgRNA library to identify cellular factors that inhibit HIV-1 infection of human CD4+ T cells. MT4 cells were transduced with a CRISPR/Cas9 sgRNA library and infected with nef-deficient HIV-1NL4-3 expressing ganciclovir-sensitive thymidine kinase, thus enabling selection of HIV-1-resistant cells for analysis of enriched sgRNAs. After validation of screen hits, multiple host factors essential for HIV-1 infection were identified, including SET (SET nuclear proto-oncogene) and ANP32A (acidic nuclear phosphoprotein 32A, PP32A), which together form a histone acetylase inhibitor complex. Using multiple human cell lines and peripheral blood mononuclear cells (PBMCs) from healthy donors and HIV-1-infected individuals, we demonstrate that SET depletion increased HIV-1 infectivity by augmenting DNA integration without significantly changing sites of integration. Conversely, SET overexpression decreased HIV-1 integration and infectivity. SET protein expression was significantly reduced in PBMCs from HIV-1-infected individuals and was downregulated by HIV-1 infection of healthy donor cells in vitro. Notably, HIV-1-induced downregulation of SET could be alleviated by inhibition of the protease granzyme A. Altogether, we have identified cellular inhibitors of HIV-1 infection on a genome-wide scale, which affords new insight into host-virus interactions and may provide new strategies for HIV-1 treatment.

DNA-encoded library (DEL) technology is a powerful tool for small molecule identification in drug discovery, yet the reported DEL selection strategies were applied primarily on protein targets in either purified form or in cellular context. To expand the application of this technology, we employed DEL selection on an RNA target HIV-1 TAR (trans-acting responsive region), but found that the majority of signals were resulted from false positive DNA-RNA binding. We thus developed an optimized selection strategy utilizing RNA patches and competitive elution to minimize unwanted DNA binding, followed by k-mer analysis and motif search to differentiate false positive signal. This optimized strategy resulted in a very clean background in a DEL selection against Escherichia coli FMN Riboswitch, and the enriched compounds were determined with double digit nanomolar binding affinity, as well as similar potency in functional FMN competition assay. These results demonstrated the feasibility of small molecule identification against RNA targets using DEL selection. The developed experimental and computational strategy provided a promising opportunity for RNA ligand screening and expanded the application of DEL selection to a much wider context in drug discovery.

Human ribosomes have long been thought to be uniform factories with little regulatory function. Accumulating evidence emphasizes the heterogeneity of ribosomal protein (RP) expression in specific cellular functions and development. However, a systematic understanding of functional relevance of RPs is lacking. Here, we surveyed translational and transcriptional changes after individual knockdown of 75 RPs, 44 from the large subunit (60S) and 31 from the small subunit (40S), by Ribo-seq and RNA-seq analyses. Deficiency of individual RPs altered specific subsets of genes transcriptionally and translationally. RP genes were under cotranslational regulation upon ribosomal stress, and deficiency of the 60S RPs and the 40S RPs had opposite effects. RP deficiency altered the expression of genes related to eight major functional classes, including the cell cycle, cellular metabolism, signal transduction and development. 60S RP deficiency led to greater inhibitory effects on cell growth than did 40S RP deficiency, through P53 signaling. Particularly, we showed that eS8/RPS8 deficiency stimulated apoptosis while eL13/RPL13 or eL18/RPL18 deficiency promoted senescence. We also validated the phenotypic impacts of uL5/RPL11 and eL15/RPL15 deficiency on retina development and angiogenesis, respectively. Overall, our study provides a valuable resource for and novel insights into ribosome regulation in cellular activities, development and diseases.

Flexible regions in biomolecular complexes, although crucial to understanding structure-function relationships, are often unclear in high-resolution crystal structures. In this study, we showed that single-molecule techniques, in combination with computational modeling, can characterize dynamic conformations not resolved by high-resolution structure determination methods. Taking two Pif1 helicases (ScPif1 and BsPif1) as model systems, we found that, besides a few tightly bound nucleotides, adjacent solvent-exposed nucleotides interact dynamically with the helicase surfaces. The whole nucleotide segment possessed curved conformations and covered the two RecA-like domains of the helicases, which are essential for the inch-worm mechanism. The synergetic approach reveals that the interactions between the exposed nucleotides and the helicases could be reduced by large stretching forces or electrostatically shielded with high-concentration salt, subsequently resulting in reduced translocation rates of the helicases. The dynamic interactions between the exposed nucleotides and the helicases underlay the force- and salt-dependences of their enzymatic activities. The present single-molecule based approach complements high-resolution structural methods in deciphering the molecular mechanisms of the helicases.

Riboswitches are gene regulatory elements located in untranslated mRNA regions. They bind inducer molecules with high affinity and specificity. Cyclic-di-nucleotide-sensing riboswitches are major regulators of genes for the environment, membranes and motility (GEMM) of bacteria. Up to now, structural probing assays or crystal structures have provided insight into the interaction between cyclic-di-nucleotides and their corresponding riboswitches. ITC analysis, NMR analysis and computational modeling allowed us to gain a detailed understanding of the gene regulation mechanisms for the Cd1 (Clostridium difficile) and for the pilM (Geobacter metallireducens) riboswitches and their respective di-nucleotides c-di-GMP and c-GAMP. Binding capability showed a 25 nucleotide (nt) long window for pilM and a 61 nt window for Cd1. Within this window, binding affinities ranged from 35 μM to 0.25 μM spanning two orders of magnitude for Cd1 and pilM showing a strong dependence on competing riboswitch folds. Experimental results were incorporated into a Markov simulation to further our understanding of the transcriptional folding pathways of riboswitches. Our model showed the ability to predict riboswitch gene regulation and its dependence on transcription speed, pausing and ligand concentration.

DICER1 syndrome is a cancer pre-disposition disorder caused by mutations that disrupt the function of DICER1 in miRNA processing. Studying the molecular, cellular and oncogenic effects of these mutations can reveal novel mechanisms that control cell homeostasis and tumor biology. Here, we conduct the first analysis of pathogenic DICER1 syndrome allele from the DICER1 3'UTR. We find that the DICER1 syndrome allele, rs1252940486, abolishes interaction with the PUMILIO RNA binding protein with the DICER1 3'UTR, resulting in the degradation of the DICER1 mRNA by AUF1. This single mutational event leads to diminished DICER1 mRNA and protein levels, and widespread reprogramming of miRNA networks. The in-depth characterization of the rs1252940486 DICER1 allele, reveals important post-transcriptional regulatory events that control DICER1 levels.

Aha1 is a co-chaperone of heat shock protein 90 (HSP90), and it stimulates the ATPase activity of HSP90 to promote the folding of its client proteins. By employing ascorbate peroxidase (APEX)-based proximity labeling and proteomic analysis, we identified over 30 proteins exhibiting diminished abundances in the proximity proteome of HSP90 in HEK293T cells upon genetic depletion of Aha1. Dicer1 is a top-ranked protein, and we confirmed its interactions with HSP90 and Aha1 by immunoprecipitation followed by western blot analysis. Genetic depletion of Aha1 and pharmacological inhibition of HSP90 both led to reduced levels of Dicer1 protein. Additionally, HSP90 and Aha1 bind preferentially to newly translated Dicer1. Reconstitution of Aha1-depleted cells with wild-type Aha1 substantially rescued Dicer1 protein level, and a lower level of restoration was observed for complementation with the HSP90-binding-defective Aha1-E67K, whereas an Aha1 mutant lacking the first 20 amino acids-which abolishes its chaperone activity-failed to rescue Dicer1 protein level. Moreover, knockdown of Aha1 and inhibition of HSP90 led to diminished levels of mature microRNAs (miRNAs), but not their corresponding primary miRNAs. Together, we uncovered a novel mechanism of HSP90 and Aha1 in regulating the miRNA pathway through promoting the folding of Dicer1 protein, and we also demonstrated that Aha1 modulates this process by acting as an autonomous chaperone and a co-chaperone for HSP90.

We analysed DNA methylation data from 30 datasets comprising 3474 individuals, 19 tissues and 8 ethnicities at CpGs covered by the Illumina450K array. We identified 4143 hypervariable CpGs ('hvCpGs') with methylation in the top 5% most variable sites across multiple tissues and ethnicities. hvCpG methylation was influenced but not determined by genetic variation, and was not linked to probe reliability, epigenetic drift, age, sex or cell heterogeneity effects. hvCpG methylation tended to covary across tissues derived from different germ-layers and hvCpGs were enriched for proximity to ERV1 and ERVK retrovirus elements. hvCpGs were also enriched for loci previously associated with periconceptional environment, parent-of-origin-specific methylation, and distinctive methylation signatures in monozygotic twins. Together, these properties position hvCpGs as strong candidates for studying how stochastic and/or environmentally influenced DNA methylation states which are established in the early embryo and maintained stably thereafter can influence life-long health and disease.

The rapid transcriptional response to the transcription factor, glucocorticoid receptor (GR), including gene activation or repression, is mediated by the spatial association of genes with multiple GR binding sites (GBSs) over large genomic distances. However, only a minority of the GBSs have independent GR-mediated activating capacity, and GBSs with independent repressive activity were rarely reported. To understand the positive and negative effects of GR we mapped the regulatory environment of its gene targets. We show that the chromatin interaction networks of GR-activated and repressed genes are spatially separated and vary in the features and configuration of their GBS and other non-GBS regulatory elements. The convergence of the KLF4 pathway in GR-activated domains and the STAT6 pathway in GR-repressed domains, impose opposite transcriptional effects to GR, independent of hormone application. Moreover, the ROR and Rev-erb transcription factors serve as positive and negative regulators, respectively, of GR-mediated gene activation. We found that the spatial crosstalk between GBSs and non-GBSs provides a physical platform for sequestering the Ep300 co-activator from non-GR regulatory loci in both GR-activated and -repressed gene compartments. While this allows rapid gene repression, Ep300 recruitment to GBSs is productive specifically in the activated compartments, thus providing the basis for gene induction.

We report the discovery and functional characterization of a new bacterial tRNA species. The tRNA-Asp-AUC, from a fast-growing desert streptomycete, decodes GAU codons. In the absence of queuosine tRNA anticodon modification in streptomycetes, the new tRNA circumvents inefficient wobble base-pairing during translation. The tRNA, which is constitutively expressed, greatly enhances synthesis of 4 different antibiotics in the model mesophilic species Streptomyces coelicolor, including the product of a so-called cryptic pathway, and increases yields of medically-important antibiotics in other species. This can be rationalised due to increased expression of both pleiotropic and pathway-specific transcriptional activators of antibiotic biosynthesis whose genes generally possess one or more GAT codons; the frequency of this codon in these gene sets is significantly higher than the average for streptomycete genes. In addition, the tRNA enhances production of cobalamin, a precursor of S-adenosyl methionine, itself an essential cofactor for synthesis of many antibiotics. The results establish a new paradigm of inefficient wobble base-pairing involving GAU codons as an evolved strategy to regulate gene expression and, in particular, antibiotic biosynthesis. Circumventing this by expression of the new cognate tRNA offers a generic strategy to increase antibiotic yields and to expand the repertoire of much-needed new bioactive metabolites produced by these valuable bacteria.

Mature microRNAs are bound by a member of the Argonaute (Ago1-4) protein family, forming the core of the RNA-induced silencing complex (RISC). Association of RISC with target mRNAs results in ribonucleoprotein (RNP) assembly involved in translational silencing or RNA degradation. Yet, the dynamics of RNP assembly and its underlying functional implications are unknown. Here, we have characterized the role of the RNA-binding protein Staufen2, a candidate Ago interactor, in RNP assembly. Staufen2 depletion resulted in the upregulation of Ago1/2 and the RISC effector proteins Ddx6 and Dcp1a. This upregulation was accompanied by the displacement of Ago1/2 from processing bodies, large RNPs implicated in RNA storage, and subsequent association of Ago2 with polysomes. In parallel, Staufen2 deficiency decreased global translation and increased dendritic branching. As the observed phenotypes can be rescued by Ago1/2 knockdown, we propose a working model in which both Staufen2 and Ago proteins depend on each other and contribute to neuronal homeostasis.

The nucleolus is a subnuclear membraneless compartment intimately involved in ribosomal RNA synthesis, ribosome biogenesis and stress response. Multiple optogenetic devices have been developed to manipulate nuclear protein import and export, but molecular tools tailored for remote control over selective targeting or partitioning of cargo proteins into subnuclear compartments capable of phase separation are still limited. Here, we report a set of single-component photoinducible nucleolus-targeting tools, designated pNUTs, to enable rapid and reversible nucleoplasm-to-nucleolus shuttling, with the half-lives ranging from milliseconds to minutes. pNUTs allow both global protein infiltration into nucleoli and local delivery of cargoes into the outermost layer of the nucleolus, the granular component. When coupled with the amyotrophic lateral sclerosis (ALS)-associated C9ORF72 proline/arginine-rich dipeptide repeats, pNUTs allow us to photomanipulate poly-proline-arginine nucleolar localization, perturb nucleolar protein nucleophosmin 1 and suppress nascent protein synthesis. pNUTs thus expand the optogenetic toolbox by permitting light-controllable interrogation of nucleolar functions and precise induction of ALS-associated toxicity in cellular models.

Currently, gapmer antisense oligonucleotide (ASO) therapeutics are under clinical development for the treatment of various diseases, including previously intractable human disorders; however, they have the potential to induce hepatotoxicity. Although several groups have reported the reduced hepatotoxicity of gapmer ASOs following chemical modifications of sugar residues or internucleotide linkages, only few studies have described nucleobase modifications to reduce hepatotoxicity. In this study, we introduced single or multiple combinations of 17 nucleobase derivatives, including four novel derivatives, into hepatotoxic locked nucleic acid gapmer ASOs and examined their effects on hepatotoxicity. The results demonstrated successful identification of chemical modifications that strongly reduced the hepatotoxicity of gapmer ASOs. This approach expands the ability to design gapmer ASOs with optimal therapeutic profiles.

G-quadruplexes (G4s) are well known non-canonical DNA secondary structures that can form in human cells. Most of the tools available to investigate G4-biology rely on small molecule ligands that stabilise these structures. However, the development of probes that disrupt G4s is equally important to study their biology. In this study, we investigated the disruption of G4s using Locked Nucleic Acids (LNA) as invader probes. We demonstrated that strategic positioning of LNA-modifications within short oligonucleotides (10 nts.) can significantly accelerate the rate of G4-disruption. Single-molecule experiments revealed that short LNA-probes can promote disruption of G4s with mechanical stability sufficient to stall polymerases. We corroborated this using a single-step extension assay, revealing that short LNA-probes can relieve replication dependent polymerase-stalling at G4 sites. We further demonstrated the potential of such LNA-based probes to study G4-biology in cells. By using a dual-luciferase assay, we found that short LNA probes can enhance the expression of c-KIT to levels similar to those observed when the c-KIT promoter is mutated to prevent the formation of the c-KIT1 G4. Collectively, our data suggest a potential use of rationally designed LNA-modified oligonucleotides as an accessible chemical-biology tool for disrupting individual G4s and interrogating their biological functions in cells.

DNA supercoiling acts as a global transcriptional regulator in bacteria, but the promoter sequence or structural determinants controlling its effect remain unclear. It was previously proposed to modulate the torsional angle between the -10 and -35 hexamers, and thereby regulate the formation of the closed-complex depending on the length of the 'spacer' between them. Here, we develop a thermodynamic model of this notion based on DNA elasticity, providing quantitative and parameter-free predictions of the relative activation of promoters containing a short versus long spacer when the DNA supercoiling level is varied. The model is tested through an analysis of in vitro and in vivo expression assays of mutant promoters with variable spacer lengths, confirming its accuracy for spacers ranging from 15 to 19 nucleotides, except those of 16 nucleotides where other regulatory mechanisms likely overcome the effect of this specific step. An analysis at the whole-genome scale in Escherichia coli then demonstrates a significant effect of the spacer length on the genomic expression after transient or inheritable superhelical variations, validating the model's predictions. Altogether, this study shows an example of mechanical constraints associated to promoter binding by RNA Polymerase underpinning a basal and global regulatory mechanism.

The nuanced mechanisms driving primordial germ cells (PGC) specification remain incompletely understood since genome-wide transcriptional regulation in developing PGCs has previously only been defined indirectly. Here, using SLAMseq analysis, we determined genome-wide transcription rates during the differentiation of embryonic stem cells (ESCs) to form epiblast-like (EpiLC) cells and ultimately PGC-like cells (PGCLCs). This revealed thousands of genes undergoing bursts of transcriptional induction and rapid shut-off not detectable by RNAseq analysis. Our SLAMseq datasets also allowed us to infer RNA turnover rates, which revealed thousands of mRNAs stabilized and destabilized during PGCLC specification. mRNAs tend to be unstable in ESCs and then are progressively stabilized as they differentiate. For some classes of genes, mRNA turnover regulation collaborates with transcriptional regulation, but these processes oppose each other in a surprisingly high frequency of genes. To test whether regulated mRNA turnover has a physiological role in PGC development, we examined three genes that we found were regulated by RNA turnover: Sox2, Klf2 and Ccne1. Circumvention of their regulated RNA turnover severely impaired the ESC-to-EpiLC and EpiLC-to-PGCLC transitions. Our study demonstrates the functional importance of regulated RNA stability in germline development and provides a roadmap of transcriptional and post-transcriptional regulation during germline specification.

Histone H3 Lysine 9 (H3K9) methylation, a characteristic mark of heterochromatin, is progressively implemented during development to contribute to cell fate restriction as differentiation proceeds. Accordingly, in undifferentiated and pluripotent mouse Embryonic Stem (ES) cells the global levels of H3K9 methylation are rather low and increase only upon differentiation. How global H3K9 methylation levels are coupled with the loss of pluripotency remains largely unknown. Here, we identify SUV39H1, a major H3K9 di- and tri-methylase, as an indirect target of the pluripotency network of Transcription Factors (TFs). We find that pluripotency TFs, principally OCT4, activate the expression of Suv39h1as, an antisense long non-coding RNA to Suv39h1. In turn, Suv39h1as downregulates Suv39h1 transcription in cis via a mechanism involving the modulation of the chromatin status of the locus. The targeted deletion of the Suv39h1as promoter region triggers increased SUV39H1 expression and H3K9me2 and H3K9me3 levels, affecting all heterochromatic regions, particularly peri-centromeric major satellites and retrotransposons. This increase in heterochromatinization efficiency leads to accelerated and more efficient commitment into differentiation. We report, therefore, a simple genetic circuitry coupling the genetic control of pluripotency with the global efficiency of H3K9 methylation associated with a major cell fate restriction, the irreversible loss of pluripotency.

R-loops are three-stranded nucleic acid structures formed from the hybridization of RNA and DNA. While the pathological consequences of R-loops have been well-studied to date, the locations, classes, and dynamics of physiological R-loops remain poorly understood. R-loop mapping studies provide insight into R-loop dynamics, but their findings are challenging to generalize. This is due to the narrow biological scope of individual studies, the limitations of each mapping modality, and, in some cases, poor data quality. In this study, we reprocessed 810 R-loop mapping datasets from a wide array of biological conditions and mapping modalities. From this data resource, we developed an accurate R-loop data quality control method, and we reveal the extent of poor-quality data within previously published studies. We then identified a set of high-confidence R-loop mapping samples and used them to define consensus R-loop sites called 'R-loop regions' (RL regions). In the process, we identified a stark divergence between RL regions detected by S9.6 and dRNH-based mapping methods, particularly with respect to R-loop size, location, and colocalization with RNA binding factors. Taken together, this work provides a much-needed method to assess R-loop data quality and offers novel context regarding the differences between dRNH- and S9.6-based R-loop mapping approaches.

As of today, the majority of environmental microorganisms remain uncultured and is therefore referred to as 'microbial dark matter' (MDM). Hence, genomic insights into these organisms are limited to cultivation-independent approaches such as single-cell- and metagenomics. However, without access to cultured representatives for verifying correct taxon-assignments, MDM genomes may cause potentially misleading conclusions based on misclassified or contaminant contigs, thereby obfuscating our view on the uncultured microbial majority. Moreover, gradual database contaminations by past genome submissions can cause error propagations which affect present as well as future comparative genome analyses. Consequently, strict contamination detection and filtering need to be applied, especially in the case of uncultured MDM genomes. Current genome reporting standards, however, emphasize completeness over purity and the de facto gold standard genome assessment tool, checkM, discriminates against uncultured taxa and fragmented genomes. To tackle these issues, we present a novel contig classification, screening, and filtering workflow and corresponding open-source python implementation called MDMcleaner, which was tested and compared to other tools on mock and real datasets. MDMcleaner revealed substantial contaminations overlooked by current screening approaches and sensitively detects misattributed contigs in both novel genomes and the underlying reference databases, thereby greatly improving our view on 'microbial dark matter'.

The dynamic rearrangement of chromatin is critical for gene regulation, but mapping both the spatial organization of chromatin and its dynamics remains a challenge. Many structural conformations are too small to be resolved via conventional fluorescence microscopy and the long acquisition time of super-resolution photoactivated localization microscopy (PALM) precludes the structural characterization of chromatin below the optical diffraction limit in living cells due to chromatin motion. Here we develop a correlative conventional fluorescence and PALM imaging approach to quantitatively map time-averaged chromatin structure and dynamics below the optical diffraction limit in living cells. By assigning localizations to a locus as it moves, we reliably discriminate between bound and unbound dCas9 molecules, whose mobilities overlap. Our approach accounts for changes in DNA mobility and relates local chromatin motion to larger scale domain movement. In our experimental system, we show that compacted telomeres move faster and have a higher density of bound dCas9 molecules, but the relative motion of those molecules is more restricted than in less compacted telomeres. Correlative conventional and PALM imaging therefore improves the ability to analyze the mobility and time-averaged nanoscopic structural features of locus specific chromatin with single molecule sensitivity and yields unprecedented insights across length and time scales.

Single-molecule techniques such as optical tweezers and fluorescence imaging are powerful tools for probing the biophysics of DNA and DNA-protein interactions. The application of these methods requires efficient approaches for creating designed DNA structures with labels for binding to a surface or microscopic beads. In this paper, we develop a simple and fast technique for making a diverse range of such DNA constructs by combining PCR amplicons and synthetic oligonucleotides using golden gate assembly rules. We demonstrate high yield fabrication of torsionally-constrained duplex DNA up to 10 kbp in length and a variety of DNA hairpin structures. We also show how tethering to a cross-linked antibody substrate significantly enhances measurement lifetime under high force. This rapid and adaptable fabrication method streamlines the assembly of DNA constructs for single molecule biophysics.

Post-transcriptional modifications can impact the stability and functionality of many different classes of RNA molecules and are an especially important aspect of tRNA regulation. It is hypothesized that cells can orchestrate rapid responses to changing environmental conditions by adjusting the specific types and levels of tRNA modifications. We uncovered strong evidence in support of this tRNA global regulation hypothesis by examining effects of the well-conserved tRNA modifying enzyme MiaA in extraintestinal pathogenic Escherichia coli (ExPEC), a major cause of urinary tract and bloodstream infections. MiaA mediates the prenylation of adenosine-37 within tRNAs that decode UNN codons, and we found it to be crucial to the fitness and virulence of ExPEC. MiaA levels shifted in response to stress via a post-transcriptional mechanism, resulting in marked changes in the amounts of fully modified MiaA substrates. Both ablation and forced overproduction of MiaA stimulated translational frameshifting and profoundly altered the ExPEC proteome, with variable effects attributable to UNN content, changes in the catalytic activity of MiaA, or availability of metabolic precursors. Cumulatively, these data indicate that balanced input from MiaA is critical for optimizing cellular responses, with MiaA acting much like a rheostat that can be used to realign global protein expression patterns.

Expression of therapeutically important proteins has benefited dramatically from the advent of chemically modified mRNAs that feature decreased lability and immunogenicity. This had a momentous effect on the rapid development of COVID-19 mRNA vaccines. Incorporation of the naturally occurring pseudouridine (Ψ) or N1-methyl-pseudouridine (N1mΨ) into in vitro transcribed mRNAs prevents the activation of unwanted immune responses by blocking eIF2α phosphorylation, which inhibits translation. Here, we report that Ψs in luciferase (Luc) mRNA exacerbate translation pausing in nuclease-untreated rabbit reticulocyte lysate (uRRL) and promote the formation of high-order-ribosome structures. The major deceleration of elongation occurs at the Ψ-rich nucleotides 1294-1326 of Ψ-Luc mRNA and results in premature termination of translation. The impairment of translation is mainly due to the shortage of membranous components. Supplementing uRRL with canine microsomal membranes (CMMs) relaxes the impediments to ribosome movement, resolves collided ribosomes, and greatly enhances full-size luciferase production. CMMs also strongly stimulated an extremely inefficient translation of N1mΨ-Luc mRNA in uRRL. Evidence is presented that translational pausing can promote membrane recruitment of polysomes with nascent polypeptides that lack a signal sequence. Our results highlight an underappreciated role of membrane binding to polysomes in the prevention of ribosome collision and premature release of nascent polypeptides.

The ribose 2'-hydroxyl is the key chemical difference between RNA and DNA and primary source of their divergent structural and functional characteristics. Macromolecular X-ray diffraction experiments typically do not reveal the positions of hydrogen atoms. Thus, standard crystallography cannot determine 2'-OH orientation (H2'-C2'-O2'-HO2' torsion angle) and its potential roles in sculpting the RNA backbone and the expansive fold space. Here, we report the first neutron crystal structure of an RNA, the Escherichia coli rRNA Sarcin-Ricin Loop (SRL). 2'-OD orientations were established for all 27 residues and revealed O-D bonds pointing toward backbone (O3', 13 observations), nucleobase (11) or sugar (3). Most riboses in the SRL stem region show a 2'-OD backbone-orientation. GAGA-tetraloop riboses display a 2'-OD base-orientation. An atypical C2'-endo sugar pucker is strictly correlated with a 2'-OD sugar-orientation. Neutrons reveal the strong preference of the 2'-OH to donate in H-bonds and that 2'-OH orientation affects both backbone geometry and ribose pucker. We discuss 2'-OH and water molecule orientations in the SRL neutron structure and compare with results from a solution phase 10 μs MD simulation. We demonstrate that joint cryo-neutron/X-ray crystallography offers an all-in-one approach to determine the complete structural properties of RNA, i.e. geometry, conformation, protonation state and hydration structure.

Crosslink repair depends on the Fanconi anemia pathway and translesion synthesis polymerases that replicate over unhooked crosslinks. Translesion synthesis is regulated via ubiquitination of PCNA, and independently via translesion synthesis polymerase REV1. The division of labor between PCNA-ubiquitination and REV1 in interstrand crosslink repair is unclear. Inhibition of either of these pathways has been proposed as a strategy to increase cytotoxicity of platinating agents in cancer treatment. Here, we defined the importance of PCNA-ubiquitination and REV1 for DNA in mammalian ICL repair. In mice, loss of PCNA-ubiquitination, but not REV1, resulted in germ cell defects and hypersensitivity to cisplatin. Loss of PCNA-ubiquitination, but not REV1 sensitized mammalian cancer cell lines to cisplatin. We identify polymerase Kappa as essential in tolerating DNA damage-induced lesions, in particular cisplatin lesions. Polk-deficient tumors were controlled by cisplatin treatment and it significantly delayed tumor outgrowth and increased overall survival of tumor bearing mice. Our results indicate that PCNA-ubiquitination and REV1 play distinct roles in DNA damage tolerance. Moreover, our results highlight POLK as a critical TLS polymerase in tolerating multiple genotoxic lesions, including cisplatin lesions. The relative frequent loss of Polk in cancers indicates an exploitable vulnerability for precision cancer medicine.

Transcription initiation is the first step in gene expression, and is therefore strongly regulated in all domains of life. The RNA polymerase (RNAP) first associates with the initiation factor $\sigma$ to form a holoenzyme, which binds, bends and opens the promoter in a succession of reversible states. These states are critical for transcription regulation, but remain poorly understood. Here, we addressed the mechanism of open complex formation by monitoring its assembly/disassembly kinetics on individual consensus lacUV5 promoters using high-throughput single-molecule magnetic tweezers. We probed the key protein-DNA interactions governing the open-complex formation and dissociation pathway by modulating the dynamics at different concentrations of monovalent salts and varying temperatures. Consistent with ensemble studies, we observed that RNAP-promoter open (RPO) complex is a stable, slowly reversible state that is preceded by a kinetically significant open intermediate (RPI), from which the holoenzyme dissociates. A strong anion concentration and type dependence indicates that the RPO stabilization may involve sequence-independent interactions between the DNA and the holoenzyme, driven by a non-Coulombic effect consistent with the non-template DNA strand interacting with $\sigma$ and the RNAP $\beta$ subunit. The temperature dependence provides the energy scale of open-complex formation and further supports the existence of additional intermediates.

Although the route to generate microRNAs (miRNAs) is often depicted as a linear series of sequential and constitutive cleavages, we now appreciate multiple alternative pathways as well as diverse strategies to modulate their processing and function. Here, we identify an unusually profound regulatory role of conserved loop sequences in vertebrate pre-mir-144, which are essential for its cleavage by the Dicer RNase III enzyme in human and zebrafish models. Our data indicate that pre-mir-144 dicing is positively regulated via its terminal loop, and involves the ILF3 complex (NF90 and its partner NF45/ILF2). We provide further evidence that this regulatory switch involves reshaping of the pre-mir-144 apical loop into a structure that is appropriate for Dicer cleavage. In light of our recent findings that mir-144 promotes the nuclear biogenesis of its neighbor mir-451, these data extend the complex hierarchy of nuclear and cytoplasmic regulatory events that can control the maturation of clustered miRNAs.

CRISPR-based precise gene-editing requires simultaneous delivery of multiple components into living cells, rapidly exceeding the cargo capacity of traditional viral vector systems. This challenge represents a major roadblock to genome engineering applications. Here we exploit the unmatched heterologous DNA cargo capacity of baculovirus to resolve this bottleneck in human cells. By encoding Cas9, sgRNA and Donor DNAs on a single, rapidly assembled baculoviral vector, we achieve with up to 30% efficacy whole-exon replacement in the intronic β-actin (ACTB) locus, including site-specific docking of very large DNA payloads. We use our approach to rescue wild-type podocin expression in steroid-resistant nephrotic syndrome (SRNS) patient derived podocytes. We demonstrate single baculovirus vectored delivery of single and multiplexed prime-editing toolkits, achieving up to 100% cleavage-free DNA search-and-replace interventions without detectable indels. Taken together, we provide a versatile delivery platform for single base to multi-gene level genome interventions, addressing the currently unmet need for a powerful delivery system accommodating current and future CRISPR technologies without the burden of limited cargo capacity.

Replication of the human genome initiates within broad zones of ∼150 kb. The extent to which firing of individual DNA replication origins within initiation zones is spatially stochastic or localised at defined sites remains a matter of debate. A thorough characterisation of the dynamic activation of origins within initiation zones is hampered by the lack of a high-resolution map of both their position and efficiency. To address this shortcoming, we describe a modification of initiation site sequencing (ini-seq), based on density substitution. Newly replicated DNA is rendered 'heavy-light' (HL) by incorporation of BrdUTP while unreplicated DNA remains 'light-light' (LL). Replicated HL-DNA is separated from unreplicated LL-DNA by equilibrium density gradient centrifugation, then both fractions are subjected to massive parallel sequencing. This allows precise mapping of 23,905 replication origins simultaneously with an assignment of a replication initiation efficiency score to each. We show that origin firing within early initiation zones is not randomly distributed. Rather, origins are arranged hierarchically with a set of very highly efficient origins marking zone boundaries. We propose that these origins explain much of the early firing activity arising within initiation zones, helping to unify the concept of replication initiation zones with the identification of discrete replication origin sites.

Labelling of oligonucleotides with dyes, targeting ligands, and other moieties has become ever more essential in life-sciences. Conventionally, modifications are introduced to oligonucleotides during solid phase synthesis by special phosphoramidites functionalised with a chemical handle or the desired functional group. In this work, we present a facile and inexpensive method to introduce modifications to oligonucleotides without the need for special phosphoramidites. Sulfonyl azides are applied to react with one or more selected phosphite intermediates during solid phase synthesis. We have prepared 11 sulfonyl azides with different chemical handles such as amine, azide, alkyne, and thiol, and we have further introduced functionalities such as pyrene, other dyes, photo-switchable azobenzenes, and a steroid. The method is compatible with current phosphoramidite-based automated oligonucleotide synthesis and serves as a simple alternative to the unstable and expensive special phosphoramidites currently used for conjugation to oligonucleotides.

The rise of antibiotic-resistant bacteria represents a major threat to global health, creating an urgent need to discover new antibiotics. Natural products derived from the genus Streptomyces represent a rich and diverse repertoire of chemical molecules from which new antibiotics are likely to be found. However, a major challenge is that the biosynthetic gene clusters (BGCs) responsible for natural product synthesis are often poorly expressed under laboratory culturing conditions, thus preventing the isolation and screening of novel chemicals. To address this, we describe a novel approach to activate silent BGCs through rewiring endogenous regulation using synthetic gene regulators based upon CRISPR-Cas. First, we refine CRISPR interference (CRISPRi) and create CRISPR activation (CRISPRa) systems that allow for highly programmable and effective gene repression and activation in Streptomyces. We then harness these tools to activate a silent BGC by perturbing its endogenous regulatory network. Together, this work advances the synthetic regulatory toolbox for Streptomyces and facilitates the programmable activation of silent BGCs for novel chemical discovery.

Erratum in Nucleic Acids Res. 2023 Jan 18;:

Polyglutamine (polyQ) diseases are a type of inherited neurodegenerative disorders caused by cytosine-adenine-guanine (CAG) trinucleotide expansion within the coding region of the disease-associated genes. We previously demonstrated that a pathogenic interaction between expanded CAG RNA and the nucleolin (NCL) protein triggers the nucleolar stress and neuronal cell death in polyQ diseases. However, mechanisms behind the molecular interaction remain unknown. Here, we report a 1.45 Å crystal structure of the r(CAG)5 oligo that comprises a full A'-form helical turn with widened grooves. Based on this structure, we simulated a model of r(CAG)5 RNA complexed with the RNA recognition motif 2 (RRM2) of NCL and identified NCL residues that are critical for its binding to CAG RNA. Combined with in vitro and in vivo site-directed mutagenesis studies, our model reveals that CAG RNA binds to NCL sites that are not important for other cellular functions like gene expression and rRNA synthesis regulation, indicating that toxic CAG RNA interferes with NCL functions by sequestering it. Accordingly, an NCL mutant that is aberrant in CAG RNA-binding could rescue RNA-induced cytotoxicity effectively. Taken together, our study provides new molecular insights into the pathogenic mechanism of polyQ diseases mediated by NCL-CAG RNA interaction.

Targeted chromosomal insertion of large genetic payloads in human cells leverages and broadens synthetic biology and genetic therapy efforts. Yet, obtaining large-scale gene knock-ins remains particularly challenging especially in hard-to-transfect stem and progenitor cells. Here, fully viral gene-deleted adenovector particles (AdVPs) are investigated as sources of optimized high-specificity CRISPR-Cas9 nucleases and donor DNA constructs tailored for targeted insertion of full-length dystrophin expression units (up to 14.8-kb) through homologous recombination (HR) or homology-mediated end joining (HMEJ). In muscle progenitor cells, donors prone to HMEJ yielded higher CRISPR-Cas9-dependent genome editing frequencies than HR donors, with values ranging between 6% and 34%. In contrast, AdVP transduction of HR and HMEJ substrates in induced pluripotent stem cells (iPSCs) resulted in similar CRISPR-Cas9-dependent genome editing levels. Notably, when compared to regular iPSCs, in p53 knockdown iPSCs, CRISPR-Cas9-dependent genome editing frequencies increased up to 6.7-fold specifically when transducing HMEJ donor constructs. Finally, single DNA molecule analysis by molecular combing confirmed that AdVP-based genome editing achieves long-term complementation of DMD-causing mutations through the site-specific insertion of full-length dystrophin expression units. In conclusion, AdVPs are a robust and flexible platform for installing large genomic edits in human cells and p53 inhibition fosters HMEJ-based genome editing in iPSCs.

Genomic integration techniques offer opportunities for generation of engineered microorganisms with improved or even entirely new functions but are currently limited by inability for efficient insertion of long genetic payloads due to multiplexing. Herein, using Shewanella oneidensis MR-1 as a model, we developed an optimized CRISPR-associated transposase from cyanobacteria Scytonema hofmanni (ShCAST system), which enables programmable, RNA-guided transposition of ultra-long DNA sequences (30 kb) onto bacterial chromosomes at ∼100% efficiency in a single orientation. In this system, a crRNA (CRISPR RNA) was used to target multicopy loci like insertion-sequence elements or combining I-SceI endonuclease, thereby allowing efficient single-step multiplexed or iterative DNA insertions. The engineered strain exhibited drastically improved substrate diversity and extracellular electron transfer ability, verifying the success of this system. Our work greatly expands the application range and flexibility of genetic engineering techniques and may be readily extended to other bacteria for better controlling various microbial processes.

Spliced leader trans-splicing is essential for gene expression in many eukaryotes. To elucidate the molecular mechanism of this process, we characterise the molecules associated with the Caenorhabditis elegans major spliced leader snRNP (SL1 snRNP), which donates the spliced leader that replaces the 5' untranslated region of most pre-mRNAs. Using a GFP-tagged version of the SL1 snRNP protein SNA-1 created by CRISPR-mediated genome engineering, we immunoprecipitate and identify RNAs and protein components by RIP-Seq and mass spectrometry. This reveals the composition of the SL1 snRNP and identifies associations with spliceosome components PRP-8 and PRP-19. Significantly, we identify a novel, nematode-specific protein required for SL1 trans-splicing, which we designate SNA-3. SNA-3 is an essential, nuclear protein with three NADAR domains whose function is unknown. Mutation of key residues in NADAR domains inactivates the protein, indicating that domain function is required for activity. SNA-3 interacts with the CBC-ARS2 complex and other factors involved in RNA metabolism, including SUT-1 protein, through RNA or protein-mediated contacts revealed by yeast two-hybrid assays, localisation studies and immunoprecipitations. Our data are compatible with a role for SNA-3 in coordinating trans-splicing with target pre-mRNA transcription or in the processing of the Y-branch product of the trans-splicing reaction.

Advances in genome sequencing have produced hundreds of thousands of bacterial genome sequences, many of which have integrated prophages derived from temperate bacteriophages. These prophages play key roles by influencing bacterial metabolism, pathogenicity, antibiotic resistance, and defense against viral attack. However, they vary considerably even among related bacterial strains, and they are challenging to identify computationally and to extract precisely for comparative genomic analyses. Here, we describe DEPhT, a multimodal tool for prophage discovery and extraction. It has three run modes that facilitate rapid screening of large numbers of bacterial genomes, precise extraction of prophage sequences, and prophage annotation. DEPhT uses genomic architectural features that discriminate between phage and bacterial sequences for efficient prophage discovery, and targeted homology searches for precise prophage extraction. DEPhT is designed for prophage discovery in Mycobacterium genomes but can be adapted broadly to other bacteria. We deploy DEPhT to demonstrate that prophages are prevalent in Mycobacterium strains but are absent not only from the few well-characterized Mycobacterium tuberculosis strains, but also are absent from all ∼30 000 sequenced M. tuberculosis strains.

Small molecule targeting of self-splicing RNAs like group I and II introns has been limited in part by the lack of a universal high-throughput screening platform for studies of splicing inhibition and kinetics. Here, we present the development of a molecular beacon assay for monitoring the accumulation of spliced exons during RNA splicing reactions. In this case, we applied it to the autocatalyzed reaction of the H.c.LSU group II intron found in the mitochondria of the pathogenic dimorphic fungus Histoplasma capsulatum. We find that a molecular beacon with the loop length of 18 nucleotides selectively recognizes ligated exons formed during self-splicing and exhibits high fluorescent signal upon binding of its target. We demonstrate that the fluorescent assay using molecular beacons can be successfully applied to kinetic characterization of the splicing reaction and determination of inhibition constants for small molecules. The results presented herein offer support for a molecular beacon approach to identifying small molecule inhibitors of intron splicing.

Recent technological advances have enabled the generation of large amounts of data consisting of RNA sequences and their functional activity. Here, we propose a method for extracting secondary structure features that affect the functional activity of RNA from sequence-activity data. Given pairs of RNA sequences and their corresponding bioactivity values, our method calculates position-specific structural features of the input RNA sequences, considering every possible secondary structure of each RNA. A Ridge regression model is trained using the structural features as feature vectors and the bioactivity values as response variables. Optimized model parameters indicate how secondary structure features affect bioactivity. We used our method to extract intramolecular structural features of bacterial translation initiation sites and self-cleaving ribozymes, and the intermolecular features between rRNAs and Shine-Dalgarno sequences and between U1 RNAs and splicing sites. We not only identified known structural features but also revealed more detailed insights into structure-activity relationships than previously reported. Importantly, the datasets we analyzed here were obtained from different experimental systems and differed in size, sequence length and similarity, and number of RNA molecules involved, demonstrating that our method is applicable to various types of data consisting of RNA sequences and bioactivity values.

Cellular senescence triggers various types of heterochromatin remodeling that contribute to aging. However, the age-related mechanisms that lead to these epigenetic alterations remain elusive. Here, we asked how two key aging hallmarks, telomere shortening and constitutive heterochromatin loss, are mechanistically connected during senescence. We show that, at the onset of senescence, pericentromeric heterochromatin is specifically dismantled consisting of chromatin decondensation, accumulation of DNA breakages, illegitimate recombination and loss of DNA. This process is caused by telomere shortening or genotoxic stress by a sequence of events starting from TP53-dependent downregulation of the telomere protective protein TRF2. The resulting loss of TRF2 at pericentromeres triggers DNA breaks activating ATM, which in turn leads to heterochromatin decondensation by releasing KAP1 and Lamin B1, recombination and satellite DNA excision found in the cytosol associated with cGAS. This TP53-TRF2 axis activates the interferon response and the formation of chromosome rearrangements when the cells escape the senescent growth arrest. Overall, these results reveal the role of TP53 as pericentromeric disassembler and define the basic principles of how a TP53-dependent senescence inducer hierarchically leads to selective pericentromeric dismantling through the downregulation of TRF2.

SAM and HD domain containing deoxynucleoside triphosphate triphosphohydrolase 1 (SAMHD1) is driven into its activated tetramer form by binding of GTP activator and dNTP activators/substrates. In addition, the inactive monomeric and dimeric forms of the enzyme bind to single-stranded (ss) nucleic acids. During DNA replication SAMHD1 can be phosphorylated by CDK1 and CDK2 at its C-terminal threonine 592 (pSAMHD1), localizing the enzyme to stalled replication forks (RFs) to promote their restart. Although phosphorylation has only a small effect on the dNTPase activity and ssDNA binding affinity of SAMHD1, perturbation of the native T592 by phosphorylation decreased the thermal stability of tetrameric SAMHD1 and accelerated tetramer dissociation in the absence and presence of ssDNA (∼15-fold). In addition, we found that ssDNA binds competitively with GTP to the A1 site. A full-length SAMHD1 cryo-EM structure revealed substantial dynamics in the C-terminal domain (which contains T592), which could be modulated by phosphorylation. We propose that T592 phosphorylation increases tetramer dynamics and allows invasion of ssDNA into the A1 site and the previously characterized DNA binding surface at the dimer-dimer interface. These features are consistent with rapid and regiospecific inactivation of pSAMHD1 dNTPase at RFs or other sites of free ssDNA in cells.

Processing bodies (P-bodies) are ribonucleoprotein granules that contain mRNAs, RNA-binding proteins and effectors of mRNA turnover. While P-bodies have been reported to contain translationally repressed mRNAs, a causative role for P-bodies in regulating mRNA decay has yet to be established. Enhancer of decapping protein 4 (EDC4) is a core P-body component that interacts with multiple mRNA decay factors, including the mRNA decapping (DCP2) and decay (XRN1) enzymes. EDC4 also associates with the RNA endonuclease MARF1, an interaction that antagonizes the decay of MARF1-targeted mRNAs. How EDC4 interacts with MARF1 and how it represses MARF1 activity is unclear. In this study, we show that human MARF1 and XRN1 interact with EDC4 using analogous conserved short linear motifs in a mutually exclusive manner. While the EDC4-MARF1 interaction is required for EDC4 to inhibit MARF1 activity, our data indicate that the interaction with EDC4 alone is not sufficient. Importantly, we show that P-body architecture plays a critical role in antagonizing MARF1-mediated mRNA decay. Taken together, our study suggests that P-bodies can directly regulate mRNA turnover by sequestering an mRNA decay enzyme and preventing it from interfacing with and degrading targeted mRNAs.

Acetaldehyde (AA), a by-product of ethanol metabolism, is acutely toxic due to its ability to react with various biological molecules including DNA and proteins, which can greatly impede key processes such as replication and transcription and lead to DNA damage. As such AA is classified as a group 1 carcinogen by the International Agency for Research on Cancer (IARC). Previous in vitro studies have shown that AA generates bulky adducts on DNA, with signature guanine-centered (GG→TT) mutations. However, due to its weak mutagenicity, short chemical half-life, and the absence of powerful genetic assays, there is considerable variability in reporting the mutagenic effects of AA in vivo. Here, we used an established yeast genetic reporter system and demonstrate that AA treatment is highly mutagenic to cells and leads to strand-biased mutations on guanines (G→T) at a high frequency on single stranded DNA (ssDNA). We further demonstrate that AA-derived mutations occur through lesion bypass on ssDNA by the translesion polymerase Polζ. Finally, we describe a unique mutation signature for AA, which we then identify in several whole-genome and -exome sequenced cancers, particularly those associated with alcohol consumption. Our study proposes a key mechanism underlying carcinogenesis by acetaldehyde-mutagenesis of single-stranded DNA.

Although previous studies have identified several autonomous pathway components that are required for the promotion of flowering, little is known about how these components cooperate. Here, we identified an autonomous pathway complex (AuPC) containing both known components (FLD, LD and SDG26) and previously unknown components (EFL2, EFL4 and APRF1). Loss-of-function mutations of all of these components result in increased FLC expression and delayed flowering. The delayed-flowering phenotype is independent of photoperiod and can be overcome by vernalization, confirming that the complex specifically functions in the autonomous pathway. Chromatin immunoprecipitation combined with sequencing indicated that, in the AuPC mutants, the histone modifications (H3Ac, H3K4me3 and H3K36me3) associated with transcriptional activation are increased, and the histone modification (H3K27me3) associated with transcriptional repression is reduced, suggesting that the AuPC suppresses FLC expression at least partially by regulating these histone modifications. Moreover, we found that the AuPC component SDG26 associates with FLC chromatin via a previously uncharacterized DNA-binding domain and regulates FLC expression and flowering time independently of its histone methyltransferase activity. Together, these results provide a framework for understanding the molecular mechanism by which the autonomous pathway regulates flowering time.

Ribosome serves as a universal molecular machine capable of synthesis of all the proteins in a cell. Small-molecule inhibitors, such as ribosome-targeting antibiotics, can compromise the catalytic versatility of the ribosome in a context-dependent fashion, preventing transpeptidation only between particular combinations of substrates. Classic peptidyl transferase center inhibitor chloramphenicol (CHL) fails to inhibit transpeptidation reaction when the incoming A site acceptor substrate is glycine, and the molecular basis for this phenomenon is unknown. Here, we present a set of high-resolution X-ray crystal structures that explain why CHL is unable to inhibit peptide bond formation between the incoming glycyl-tRNA and a nascent peptide that otherwise is conducive to the drug action. Our structures reveal that fully accommodated glycine residue can co-exist in the A site with the ribosome-bound CHL. Moreover, binding of CHL to a ribosome complex carrying glycyl-tRNA does not affect the positions of the reacting substrates, leaving the peptide bond formation reaction unperturbed. These data exemplify how small-molecule inhibitors can reshape the A-site amino acid binding pocket rendering it permissive only for specific amino acid residues and rejective for the other substrates extending our detailed understanding of the modes of action of ribosomal antibiotics.

EPR is a long non-coding RNA (lncRNA) that controls cell proliferation in mammary gland cells by regulating gene transcription. Here, we report on Mettl7a1 as a direct target of EPR. We show that EPR induces Mettl7a1 transcription by rewiring three-dimensional chromatin interactions at the Mettl7a1 locus. Our data indicate that METTL7A1 contributes to EPR-dependent inhibition of TGF-β signaling. METTL7A1 is absent in tumorigenic murine mammary gland cells and its human ortholog (METTL7A) is downregulated in breast cancers. Importantly, re-expression of METTL7A1 in 4T1 tumorigenic cells attenuates their transformation potential, with the putative methyltransferase activity of METTL7A1 being dispensable for its biological functions. We found that METTL7A1 localizes in the cytoplasm whereby it interacts with factors implicated in the early steps of mRNA translation, associates with ribosomes, and affects the levels of target proteins without altering mRNA abundance. Overall, our data indicates that METTL7A1-a transcriptional target of EPR-modulates translation of select transcripts.

5'-Adenylated oligonucleotides (AppOligos) are widely used for single-stranded DNA/RNA ligation in next-generation sequencing (NGS) applications such as microRNA (miRNA) profiling. The ligation between an AppOligo adapter and target molecules (such as miRNA) no longer requires ATP, thereby minimizing potential self-ligations and simplifying library preparation procedures. AppOligos can be produced by chemical synthesis or enzymatic modification. However, adenylation via chemical synthesis is inefficient and expensive, while enzymatic modification requires pre-phosphorylated substrate and additional purification. Here we cloned and characterized the Pfu RNA ligase encoded by the PF0353 gene in the hyperthermophilic archaea Pyrococcus furiosus. We further engineered fusion enzymes containing both Pfu RNA ligase and T4 polynucleotide kinase. One fusion enzyme, 8H-AP, was thermostable and can directly catalyze 5'-OH-terminated DNA substrates to adenylated products. The newly discovered Pfu RNA ligase and the engineered fusion enzyme may be useful tools for applications using AppOligos.

Alternative end joining (alt-EJ) mechanisms, such as polymerase theta-mediated end joining, are increasingly recognized as important contributors to inaccurate double-strand break repair. We previously proposed an alt-EJ model whereby short DNA repeats near a double-strand break anneal to form secondary structures that prime limited DNA synthesis. The nascent DNA then pairs with microhomologous sequences on the other break end. This synthesis-dependent microhomology-mediated end joining (SD-MMEJ) explains many of the alt-EJ repair products recovered following I-SceI nuclease cutting in Drosophila. However, sequence-specific factors that influence SD-MMEJ repair remain to be fully characterized. Here, we expand the utility of the SD-MMEJ model through computational analysis of repair products at Cas9-induced double-strand breaks for 1100 different sequence contexts. We find evidence at single nucleotide resolution for sequence characteristics that drive successful SD-MMEJ repair. These include optimal primer repeat length, distance of repeats from the break, flexibility of DNA sequence between primer repeats, and positioning of microhomology templates relative to preferred primer repeats. In addition, we show that DNA polymerase theta is necessary for most SD-MMEJ repair at Cas9 breaks. The analysis described here includes a computational pipeline that can be utilized to characterize preferred mechanisms of alt-EJ repair in any sequence context.

Genome assembly has been benefited from long-read sequencing technologies with higher accuracy and higher continuity. However, most human genome assembly require large amount of DNAs from homogeneous cell lines without keeping cell heterogeneities, since cell heterogeneity could profoundly affect haplotype assembly results. Herein, using single-cell genome long-read sequencing technology (SMOOTH-seq), we have sequenced K562 and HG002 cells on PacBio HiFi and Oxford Nanopore Technologies (ONT) platforms and conducted de novo genome assembly. For the first time, we have completed the human genome assembly with high continuity (with NG50 of ∼2 Mb using 95 individual K562 cells) at single-cell levels, and explored the impact of different assemblers and sequencing strategies on genome assembly. With sequencing data from 30 diploid individual HG002 cells of relatively high genome coverage (average coverage ∼41.7%) on ONT platform, the NG50 can reach over 1.3 Mb. Furthermore, with the assembled genome from K562 single-cell dataset, more complete and accurate set of insertion events and complex structural variations could be identified. This study opened a new chapter on the practice of single-cell genome de novo assembly.

Stalling of the transcription elongation complex formed by DNA, RNA polymerase (RNAP) and RNA presents a serious obstacle to concurrent processes due to the extremely high stability of the DNA-bound polymerase. RapA, known to remove RNAP from DNA in an ATP-dependent fashion, was identified over 50 years ago as an abundant binding partner of RNAP; however, its mechanism of action remains unknown. Here, we use single-molecule magnetic trapping assays to characterize RapA activity and begin to specify its mechanism of action. We first show that stalled RNAP resides on DNA for times on the order of 106 seconds and that increasing positive torque on the DNA reduces this lifetime. Using stalled RNAP as a substrate we show that the RapA protein stimulates dissociation of stalled RNAP from positively supercoiled DNA but not negatively supercoiled DNA. We observe that RapA-dependent RNAP dissociation is torque-sensitive, is inhibited by GreB and depends on RNA length. We propose that stalled RNAP is dislodged from DNA by RapA via backtracking in a supercoiling- and torque-dependent manner, suggesting that RapA's activity on transcribing RNAP in vivo is responsible for resolving conflicts between converging polymerase molecular motors.

Autophagy, a catabolic process to remove unnecessary or dysfunctional organelles, is triggered by various signals including nutrient starvation. Depending on the types of the nutrient deficiency, diverse sensing mechanisms and signaling pathways orchestrate for transcriptional and epigenetic regulation of autophagy. However, our knowledge about nutrient type-specific transcriptional regulation during autophagy is limited. To understand nutrient type-dependent transcriptional mechanisms during autophagy, we performed single cell RNA sequencing (scRNAseq) in the mouse embryonic fibroblasts (MEFs) with or without glucose starvation (GS) as well as amino acid starvation (AAS). Trajectory analysis using scRNAseq identified sequential induction of potential transcriptional regulators for each condition. Gene regulatory rules inferred using TENET newly identified CCAAT/enhancer binding protein γ (C/EBPγ) as a regulator of autophagy in AAS, but not GS, condition, and knockdown experiment confirmed the TENET result. Cell biological and biochemical studies validated that activating transcription factor 4 (ATF4) is responsible for conferring specificity to C/EBPγ for the activation of autophagy genes under AAS, but not under GS condition. Together, our data identified C/EBPγ as a previously unidentified key regulator under AAS-induced autophagy.

Artemis nuclease and DNA-dependent protein kinase catalytic subunit (DNA-PKcs) are key components in nonhomologous DNA end joining (NHEJ), the major repair mechanism for double-strand DNA breaks. Artemis activation by DNA-PKcs resolves hairpin DNA ends formed during V(D)J recombination. Artemis deficiency disrupts development of adaptive immunity and leads to radiosensitive T- B- severe combined immunodeficiency (RS-SCID). An activated state of Artemis in complex with DNA-PK was solved by cryo-EM recently, which showed Artemis bound to the DNA. Here, we report that the pre-activated form (basal state) of the Artemis:DNA-PKcs complex is stable on an agarose-acrylamide gel system, and suitable for cryo-EM structural analysis. Structures show that the Artemis catalytic domain is dynamically positioned externally to DNA-PKcs prior to ABCDE autophosphorylation and show how both the catalytic and regulatory domains of Artemis interact with the N-HEAT and FAT domains of DNA-PKcs. We define a mutually exclusive binding site for Artemis and XRCC4 on DNA-PKcs and show that an XRCC4 peptide disrupts the Artemis:DNA-PKcs complex. All of the findings are useful in explaining how a hypomorphic L3062R missense mutation of DNA-PKcs could lead to insufficient Artemis activation, hence RS-SCID. Our results provide various target site candidates to design disruptors for Artemis:DNA-PKcs complex formation.

For more than three decades, Ionis Pharmaceutics has pursued the challenging mission of creating a new platform for drug discovery. To overcome the numerous challenges faced required the integration of innovation across many scientific areas, despite many disappointments and failures. The approaches implemented to create and maintain a scientific environment to achieve the mission demanded the rigorous practice of science over three decades. The approaches taken are discussed in this perspective.

SETDB1 is a key regulator of lineage-specific genes and endogenous retroviral elements (ERVs) through its deposition of repressive H3K9me3 mark. Apart from its H3K9me3 regulatory role, SETDB1 has seldom been studied in terms of its other potential regulatory roles. To investigate this, a genomic survey of SETDB1 binding in mouse embryonic stem cells across multiple libraries was conducted, leading to the unexpected discovery of regions bereft of common repressive histone marks (H3K9me3, H3K27me3). These regions were enriched with the CTCF motif that is often associated with the topological regulator Cohesin. Further profiling of these non-H3K9me3 regions led to the discovery of a cluster of non-repeat loci that were co-bound by SETDB1 and Cohesin. These regions, which we named DiSCs (domains involving SETDB1 and Cohesin) were seen to be proximal to the gene promoters involved in embryonic stem cell pluripotency and lineage development. Importantly, it was found that SETDB1-Cohesin co-regulate target gene expression and genome topology at these DiSCs. Depletion of SETDB1 led to localized dysregulation of Cohesin binding thereby locally disrupting topological structures. Dysregulated gene expression trends revealed the importance of this cluster in ES cell maintenance as well as at gene 'islands' that drive differentiation to other lineages. The 'unearthing' of the DiSCs thus unravels a unique topological and transcriptional axis of control regulated chiefly by SETDB1.

Prokaryotic Argonaute proteins (pAgos) widely participate in hosts to defend against the invasion of nucleic acids. Compared with the CRISPR-Cas system, which requires a specific motif on the target and can only use RNA as guide, pAgos exhibit precise endonuclease activity on any arbitrary target sequence and can use both RNA and DNA as guide, thus rendering great potential for genome editing applications. Hitherto, most in-depth studies on the structure-function relationship of pAgos were conducted on thermophilic ones, functioning at ∼60 to 100°C, whose structures were, however, determined experimentally at much lower temperatures (20-33°C). It remains unclear whether these low-temperature structures can represent the true conformations of the thermophilic pAgos under their physiological conditions. The present work studied three pAgos, PfAgo, TtAgo and CbAgo, whose physiological temperatures differ significantly (95, 75 and 37°C). By conducting thorough experimental and simulation studies, we found that thermophilic pAgos (PfAgo and TtAgo) adopt a loosely-packed structure with a partially-melted surface at the physiological temperatures, largely different from the compact crystalline structures determined at moderate temperatures. In contrast, the mesophilic pAgo (CbAgo) assumes a compact crystalline structure at its optimal function temperature. Such a partially-disrupted structure endows thermophilic pAgos with great flexibility both globally and locally at the catalytic sites, which is crucial for them to achieve high DNA-cleavage activity. To further prove this, we incubated thermophilic pAgos with urea to purposely disrupt their structures, and the resulting cleavage activity was significantly enhanced below the physiological temperature, even at human body temperature. Further testing of many thermophilic Agos present in various thermophilic prokaryotes demonstrated that their structures are generally disrupted under physiological conditions. Therefore, our findings suggest that the highly dynamical structure with a partially-melted surface, distinct from the low-temperature crystalline structure, could be a general strategy assumed by thermophilic pAgos to achieve the high DNA-cleavage activity.

The histone variant H3.3 is encoded by two distinct genes, H3f3a and H3f3b, exhibiting identical amino-acid sequence. H3.3 is required for spermatogenesis, but the molecular mechanism of its spermatogenic function remains obscure. Here, we have studied the role of each one of H3.3A and H3.3B proteins in spermatogenesis. We have generated transgenic conditional knock-out/knock-in (cKO/KI) epitope-tagged FLAG-FLAG-HA-H3.3B (H3.3BHA) and FLAG-FLAG-HA-H3.3A (H3.3AHA) mouse lines. We show that H3.3B, but not H3.3A, is required for spermatogenesis and male fertility. Analysis of the molecular mechanism unveils that the absence of H3.3B led to alterations in the meiotic/post-meiotic transition. Genome-wide RNA-seq reveals that the depletion of H3.3B in meiotic cells is associated with increased expression of the whole sex X and Y chromosomes as well as of both RLTR10B and RLTR10B2 retrotransposons. In contrast, the absence of H3.3B resulted in down-regulation of the expression of piRNA clusters. ChIP-seq experiments uncover that RLTR10B and RLTR10B2 retrotransposons, the whole sex chromosomes and the piRNA clusters are markedly enriched of H3.3. Taken together, our data dissect the molecular mechanism of H3.3B functions during spermatogenesis and demonstrate that H3.3B, depending on its chromatin localization, is involved in either up-regulation or down-regulation of expression of defined large chromatin regions.

DNA damage drives genetic mutations that underlie the development of cancer in humans. Multiple pathways have been described in mammalian cells which can repair this damage. However, most work to date has focused upon single lesions in DNA. We present here a combinatorial system which allows assembly of duplexes containing single or multiple types of damage by ligating together six oligonucleotides containing damaged or modified bases. The combinatorial system has dual fluorescent labels allowing examination of both strands simultaneously, in order to study interactions or competition between different DNA repair pathways. Using this system, we demonstrate how repair of oxidative damage in one DNA strand can convert a mispaired T:G deamination intermediate into a T:A mutation. We also demonstrate that slow repair of a T:G mispair, relative to a U:G mispair, by the human methyl-binding domain 4 DNA glycosylase provides a competitive advantage to competing repair pathways, and could explain why CpG dinucleotides are hotspots for C to T mutations in human tumors. Data is also presented that suggests repair of closely spaced lesions in opposing strands can be repaired by a combination of short and long-patch base excision repair and simultaneous repair of multiply damage sites can potentially lead to lethal double strand breaks.