Advances in experimental and computational techniques enable the exploration of large and complex RNA 3D structures. These, in turn, reveal previously unstudied properties and motifs not characteristic for small molecules with simple architectures. Examples include entanglements of structural elements in RNA molecules and knot-like folds discovered, among others, in the genomes of RNA viruses. Recently, we presented the first classification of entanglements, determined by their topology and the type of entangled structural elements. Here, we introduce RNAspider - a web server to automatically identify, classify, and visualize primary and higher-order entanglements in RNA tertiary structures. The program applies to evaluate RNA 3D models obtained experimentally or by computational prediction. It supports the analysis of uncommon topologies in the pseudoknotted RNA structures. RNAspider is implemented as a publicly available tool with a user-friendly interface and can be freely accessed at https://rnaspider.cs.put.poznan.pl/.

ChIP-Atlas (https://chip-atlas.org) is a web service providing both GUI- and API-based data-mining tools to reveal the architecture of the transcription regulatory landscape. ChIP-Atlas is powered by comprehensively integrating all data sets from high-throughput ChIP-seq and DNase-seq, a method for profiling chromatin regions accessible to DNase. In this update, we further collected all the ATAC-seq and whole-genome bisulfite-seq data for six model organisms (human, mouse, rat, fruit fly, nematode, and budding yeast) with the latest genome assemblies. These together with ChIP-seq data can be visualized with the Peak Browser tool and a genome browser to explore the epigenomic landscape of a query genomic locus, such as its chromatin accessibility, DNA methylation status, and protein-genome interactions. This epigenomic landscape can also be characterized for multiple genes and genomic loci by querying with the Enrichment Analysis tool, which, for example, revealed that inflammatory bowel disease-associated SNPs are the most significantly hypo-methylated in neutrophils. Therefore, ChIP-Atlas provides a panoramic view of the whole epigenomic landscape. All datasets are free to download via either a simple button on the web page or an API.

DAVID is a popular bioinformatics resource system including a web server and web service for functional annotation and enrichment analyses of gene lists. It consists of a comprehensive knowledgebase and a set of functional analysis tools. Here, we report all updates made in 2021. The DAVID Gene system was rebuilt to gain coverage of more organisms, which increased the taxonomy coverage from 17 399 to 55 464. All existing annotation types have been updated, if available, based on the new DAVID Gene system. Compared with the last version, the number of gene-term records for most annotation types within the updated Knowledgebase have significantly increased. Moreover, we have incorporated new annotations in the Knowledgebase including small molecule-gene interactions from PubChem, drug-gene interactions from DrugBank, tissue expression information from the Human Protein Atlas, disease information from DisGeNET, and pathways from WikiPathways and PathBank. Eight of ten subgroups split from Uniprot Keyword annotation were assigned to specific types. Finally, we added a species parameter for uploading a list of gene symbols to minimize the ambiguity between species, which increases the efficiency of the list upload and eliminates confusion for users. These current updates have significantly expanded the Knowledgebase and enhanced the discovery power of DAVID.

Upon SARS-CoV-2 infection, viral intermediates specifically activate the IFN response through MDA5-mediated sensing and accordingly induce ADAR1 p150 expression, which might lead to viral A-to-I RNA editing. Here, we developed an RNA virus-specific editing identification pipeline, surveyed 7622 RNA-seq data from diverse types of samples infected with SARS-CoV-2, and constructed an atlas of A-to-I RNA editing sites in SARS-CoV-2. We found that A-to-I editing was dynamically regulated, varied between tissue and cell types, and was correlated with the intensity of innate immune response. On average, 91 editing events were deposited at viral dsRNA intermediates per sample. Moreover, editing hotspots were observed, including recoding sites in the spike gene that affect viral infectivity and antigenicity. Finally, we provided evidence that RNA editing accelerated SARS-CoV-2 evolution in humans during the epidemic. Our study highlights the ability of SARS-CoV-2 to hijack components of the host antiviral machinery to edit its genome and fuel its evolution, and also provides a framework and resource for studying viral RNA editing.

Mobile element insertions (MEIs) are a major class of structural variants (SVs) and have been linked to many human genetic disorders, including hemophilia, neurofibromatosis, and various cancers. However, human MEI resources from large-scale genome sequencing are still lacking compared to those for SNPs and SVs. Here, we report a comprehensive map of 36 699 non-reference MEIs constructed from 5675 genomes, comprising 2998 Chinese samples (∼26.2×, NyuWa) and 2677 samples from the 1000 Genomes Project (∼7.4×, 1KGP). We discovered that LINE-1 insertions were highly enriched in centromere regions, implying the role of chromosome context in retroelement insertion. After functional annotation, we estimated that MEIs are responsible for about 9.3% of all protein-truncating events per genome. Finally, we built a companion database named HMEID for public use. This resource represents the latest and largest genomewide study on MEIs and will have broad utility for exploration of human MEI findings.

G-quadruplexes (GQs) are non-canonical DNA structures composed of stacks of stabilized G-tetrads. GQs play an important role in a variety of biological processes and may form at telomeres and oncogene promoters among other genomic locations. Here, we investigate nine variants of telomeric DNA from Tetrahymena thermophila with the repeat (TTGGGG)n. Biophysical data indicate that the sequences fold into stable four-tetrad GQs which adopt multiple conformations according to native PAGE. Excitingly, we solved the crystal structure of two variants, TET25 and TET26. The two variants differ by the presence of a 3'-T yet adopt different GQ conformations. TET25 forms a hybrid [3 + 1] GQ and exhibits a rare 5'-top snapback feature. Consequently, TET25 contains four loops: three lateral (TT, TT, and GTT) and one propeller (TT). TET26 folds into a parallel GQ with three TT propeller loops. To the best of our knowledge, TET25 and TET26 are the first reported hybrid and parallel four-tetrad unimolecular GQ structures. The results presented here expand the repertoire of available GQ structures and provide insight into the intricacy and plasticity of the 3D architecture adopted by telomeric repeats from T. thermophila and GQs in general.

Recognition of mammalian mitochondrial promoters requires the concerted action of mitochondrial RNA polymerase (mtRNAP) and transcription initiation factors TFAM and TFB2M. In this work, we found that transcript slippage results in heterogeneity of the human mitochondrial transcripts in vivo and in vitro. This allowed us to correctly interpret the RNAseq data, identify the bona fide transcription start sites (TSS), and assign mitochondrial promoters for > 50% of mammalian species and some other vertebrates. The divergent structure of the mammalian promoters reveals previously unappreciated aspects of mtDNA evolution. The correct assignment of TSS also enabled us to establish the precise register of the DNA in the initiation complex and permitted investigation of the sequence-specific protein-DNA interactions. We determined the molecular basis of promoter recognition by mtRNAP and TFB2M, which cooperatively recognize bases near TSS in a species-specific manner. Our findings reveal a role of mitochondrial transcription machinery in mitonuclear coevolution and speciation.

Two families of DNA glycosylases (YtkR2/AlkD, AlkZ/YcaQ) have been found to remove bulky and crosslinking DNA adducts produced by bacterial natural products. Whether DNA glycosylases eliminate other types of damage formed by structurally diverse antibiotics is unknown. Here, we identify four DNA glycosylases-TxnU2, TxnU4, LldU1 and LldU5-important for biosynthesis of the aromatic polyketide antibiotics trioxacarcin A (TXNA) and LL-D49194 (LLD), and show that the enzymes provide self-resistance to the producing strains by excising the intercalated guanine adducts of TXNA and LLD. These enzymes are highly specific for TXNA/LLD-DNA lesions and have no activity toward other, less stable alkylguanines as previously described for YtkR2/AlkD and AlkZ/YcaQ. Similarly, TXNA-DNA adducts are not excised by other alkylpurine DNA glycosylases. TxnU4 and LldU1 possess unique active site motifs that provide an explanation for their tight substrate specificity. Moreover, we show that abasic (AP) sites generated from TxnU4 excision of intercalated TXNA-DNA adducts are incised by AP endonuclease less efficiently than those formed by 7mG excision. This work characterizes a distinct class of DNA glycosylase acting on intercalated DNA adducts and furthers our understanding of specific DNA repair self-resistance activities within antibiotic producers of structurally diverse, highly functionalized DNA damaging agents.

The Sleeping Beauty (SB) transposon system is a popular tool for genome engineering, but random integration into the genome carries a certain genotoxic risk in therapeutic applications. Here we investigate the role of amino acids H187, P247 and K248 in target site selection of the SB transposase. Structural modeling implicates these three amino acids located in positions analogous to amino acids with established functions in target site selection in retroviral integrases and transposases. Saturation mutagenesis of these residues in the SB transposase yielded variants with altered target site selection properties. Transposon integration profiling of several mutants reveals increased specificity of integrations into palindromic AT repeat target sequences in genomic regions characterized by high DNA bendability. The H187V and K248R mutants redirect integrations away from exons, transcriptional regulatory elements and nucleosomal DNA in the human genome, suggesting enhanced safety and thus utility of these SB variants in gene therapy applications.

Many cellular processes occur out of equilibrium. This includes site-specific unwinding in supercoiled DNA, which may play an important role in gene regulation. Here, we use the Convex Lens-induced Confinement (CLiC) single-molecule microscopy platform to study these processes with high-throughput and without artificial constraints on molecular structures or interactions. We use two model DNA plasmid systems, pFLIP-FUSE and pUC19, to study the dynamics of supercoiling-induced secondary structural transitions after perturbations away from equilibrium. We find that structural transitions can be slow, leading to long-lived structural states whose kinetics depend on the duration and direction of perturbation. Our findings highlight the importance of out-of-equilibrium studies when characterizing the complex structural dynamics of DNA and understanding the mechanisms of gene regulation.

Specific sequence features of the protospacer and protospacer-adjacent motif (PAM) are critical for efficient cleavage by CRISPR-Cas9, but current knowledge is largely derived from single-guide RNA (sgRNA) systems assessed in cultured cells. In this study, we sought to determine gRNA sequence features of a more native CRISPR-Cas9 ribonucleoprotein (RNP) complex with dual-guide RNAs (dgRNAs) composed of crRNA and tracrRNA, which has been used increasingly in recent CRISPR-Cas9 applications, particularly in zebrafish. Using both wild-type and HiFi SpCas9, we determined on-target cleavage efficiencies of 51 crRNAs in zebrafish embryos by assessing indel occurrence. Statistical analysis of these data identified novel position-specific mononucleotide features relevant to cleavage efficiencies throughout the protospacer sequence that may be unique to CRISPR-Cas9 RNPs pre-assembled with perfectly matched gRNAs. Overall features for wild-type Cas9 resembled those for HiFi Cas9, but specific differences were also observed. Mutational analysis of mononucleotide features confirmed their relevance to cleavage efficiencies. Moreover, the mononucleotide feature-based score, CRISPR-kp, correlated well with efficiencies of gRNAs reported in previous zebrafish RNP injection experiments, as well as independently tested crRNAs only in RNP format, but not with Cas9 mRNA co-injection. These findings will facilitate design of gRNA/crRNAs in genome editing applications, especially when using pre-assembled RNPs.

Ssd1, a conserved fungal RNA-binding protein, is important in stress responses, cell division and virulence. Ssd1 is closely related to Dis3L2 of the RNase II family of nucleases, but lacks catalytic activity and likely suppresses translation of bound mRNAs. Previous studies identified RNA motifs enriched in Ssd1-associated transcripts, yet the sequence requirements for Ssd1 binding are not defined. Here, we identify precise binding sites of Ssd1 on RNA using in vivo cross-linking and cDNA analysis. These sites are enriched in 5' untranslated regions of a subset of mRNAs encoding cell wall proteins. We identified a conserved bipartite motif that binds Ssd1 with high affinity in vitro. Active RNase II enzymes have a characteristic, internal RNA binding path; the Ssd1 crystal structure at 1.9 Å resolution shows that remnants of regulatory sequences block this path. Instead, RNA binding activity has relocated to a conserved patch on the surface of the protein. Structure-guided mutations of this surface prevent Ssd1 from binding RNA in vitro and phenocopy Ssd1 deletion in vivo. These studies provide a new framework for understanding the function of a pleiotropic post-transcriptional regulator of gene expression and give insights into the evolution of regulatory and binding elements in the RNase II family.

Human papillomavirus type 16 (HPV16) E1 and E6 proteins are produced from mRNAs with retained introns, but it has been unclear how these mRNAs are generated. Here, we report that hnRNP D act as a splicing inhibitor of HPV16 E1/E2- and E6/E7-mRNAs thereby generating intron-containing E1- and E6-mRNAs, respectively. N- and C-termini of hnRNP D contributed to HPV16 mRNA splicing control differently. HnRNP D interacted with the components of splicing machinery and with HPV16 RNA to exert its inhibitory function. As a result, the cytoplasmic levels of intron-retained HPV16 mRNAs were increased in the presence of hnRNP D. Association of hnRNP D with HPV16 mRNAs in the cytoplasm was observed, and this may correlate with unexpected inhibition of HPV16 E1- and E6-mRNA translation. Notably, hnRNP D40 interacted with HPV16 mRNAs in an HPV16-driven tonsillar cancer cell line and in HPV16-immortalized human keratinocytes. Furthermore, knockdown of hnRNP D in HPV16-driven cervical cancer cells enhanced production of the HPV16 E7 oncoprotein. Our results suggest that hnRNP D plays significant roles in the regulation of HPV gene expression and HPV-associated cancer development.

Comment in 10.1093/nar/gkab519.

Anti-CRISPR (Acr) proteins are encoded by many mobile genetic elements (MGEs) such as phages and plasmids to combat CRISPR-Cas adaptive immune systems employed by prokaryotes, which provide powerful tools for CRISPR-Cas-based applications. Here, we discovered nine distinct type II-A anti-CRISPR (AcrIIA24-32) families from Streptococcus MGEs and found that most Acrs can potently inhibit type II-A Cas9 orthologs from Streptococcus (SpyCas9, St1Cas9 or St3Cas9) in bacterial and human cells. Among these Acrs, AcrIIA26, AcrIIA27, AcrIIA30 and AcrIIA31 are able to block Cas9 binding to DNA, while AcrIIA24 abrogates DNA cleavage by Cas9. Notably, AcrIIA25.1 and AcrIIA32.1 can inhibit both DNA binding and DNA cleavage activities of SpyCas9, exhibiting unique anti-CRISPR characteristics. Importantly, we developed several chemically inducible anti-CRISPR variants based on AcrIIA25.1 and AcrIIA32.1 by comprising hybrids of Acr protein and the 4-hydroxytamoxifen-responsive intein, which enabled post-translational control of CRISPR-Cas9-mediated genome editing in human cells. Taken together, our work expands the diversity of type II-A anti-CRISPR families and the toolbox of Acr proteins for the chemically inducible control of Cas9-based applications.

The Spt4-Spt5 complex is conserved and essential RNA polymerase elongation factor. To investigate the role of the Spt4-Spt5 complex in non-coding transcription during development, we used the unicellular model Paramecium tetraurelia. In this organism harboring both germline and somatic nuclei, massive transcription of the entire germline genome takes place during meiosis. This phenomenon starts a series of events mediated by different classes of non-coding RNAs that control developmentally programmed DNA elimination. We focused our study on Spt4, a small zinc-finger protein encoded in P. tetraurelia by two genes expressed constitutively and two genes expressed during meiosis. SPT4 genes are not essential in vegetative growth, but they are indispensable for sexual reproduction, even though genes from both expression families show functional redundancy. Silencing of the SPT4 genes resulted in the absence of double-stranded ncRNAs and reduced levels of scnRNAs - 25 nt-long sRNAs produced from these double-stranded precursors in the germline nucleus. Moreover, we observed that the presence of a germline-specific Spt4-Spt5m complex is necessary for transfer of the scnRNA-binding PIWI protein between the germline and somatic nucleus. Our study establishes that Spt4, together with Spt5m, is essential for expression of the germline genome and necessary for developmental genome rearrangements.

Thymine dimers are a major mutagenic photoproduct induced by UV radiation. While they have been the subject of extensive theoretical and experimental investigations, questions of how DNA supercoiling affects local defect properties, or, conversely, how the presence of such defects changes global supercoiled structure, are largely unexplored. Here, we introduce a model of thymine dimers in the oxDNA forcefield, parametrized by comparison to melting experiments and structural measurements of the thymine dimer induced bend angle. We performed extensive molecular dynamics simulations of double-stranded DNA as a function of external twist and force. Compared to undamaged DNA, the presence of a thymine dimer lowers the supercoiling densities at which plectonemes and bubbles occur. For biologically relevant supercoiling densities and forces, thymine dimers can preferentially segregate to the tips of the plectonemes, where they enhance the probability of a localized tip-bubble. This mechanism increases the probability of highly bent and denatured states at the thymine dimer site, which may facilitate repair enzyme binding. Thymine dimer-induced tip-bubbles also pin plectonemes, which may help repair enzymes to locate damage. We hypothesize that the interplay of supercoiling and local defects plays an important role for a wider set of DNA damage repair systems.

Accelerated evolution of any portion of the genome is of significant interest, potentially signaling positive selection of phenotypic traits and adaptation. Accelerated evolution remains understudied for structured RNAs, despite the fact that an RNA's structure is often key to its function. RNA structures are typically characterized by compensatory (structure-preserving) basepair changes that are unexpected given the underlying sequence variation, i.e., they have evolved through negative selection on structure. We address the question of how fast the primary sequence of an RNA can change through evolution while conserving its structure. Specifically, we consider predicted and known structures in vertebrate genomes. After careful control of false discovery rates, we obtain 13 de novo structures (and three known Rfam structures) that we predict to have rapidly evolving sequences-defined as structures where the primary sequences of human and mouse have diverged at least twice as fast (1.5 times for Rfam) as nearby neutrally evolving sequences. Two of the three known structures function in translation inhibition related to infection and immune response. We conclude that rapid sequence divergence does not preclude RNA structure conservation in vertebrates, although these events are relatively rare.

The combined analysis of haplotype panels with phenotype clinical cohorts is a common approach to explore the genetic architecture of human diseases. However, genetic studies are mainly based on single nucleotide variants (SNVs) and small insertions and deletions (indels). Here, we contribute to fill this gap by generating a dense haplotype map focused on the identification, characterization, and phasing of structural variants (SVs). By integrating multiple variant identification methods and Logistic Regression Models (LRMs), we present a catalogue of 35 431 441 variants, including 89 178 SVs (≥50 bp), 30 325 064 SNVs and 5 017 199 indels, across 785 Illumina high coverage (30x) whole-genomes from the Iberian GCAT Cohort, containing a median of 3.52M SNVs, 606 336 indels and 6393 SVs per individual. The haplotype panel is able to impute up to 14 360 728 SNVs/indels and 23 179 SVs, showing a 2.7-fold increase for SVs compared with available genetic variation panels. The value of this panel for SVs analysis is shown through an imputed rare Alu element located in a new locus associated with Mononeuritis of lower limb, a rare neuromuscular disease. This study represents the first deep characterization of genetic variation within the Iberian population and the first operational haplotype panel to systematically include the SVs into genome-wide genetic studies.

Ribosome assembly is an essential process that is linked to human congenital diseases and tumorigenesis. While great progress has been made in deciphering mechanisms governing ribosome biogenesis in eukaryotes, an inventory of factors that support ribosome synthesis in human cells is still missing, in particular regarding the maturation of the large 60S subunit. Here, we performed a genome-wide RNAi screen using an imaging-based, single cell assay to unravel the cellular machinery promoting 60S subunit assembly in human cells. Our screen identified a group of 310 high confidence factors. These highlight the conservation of the process across eukaryotes and reveal the intricate connectivity of 60S subunit maturation with other key cellular processes, including splicing, translation, protein degradation, chromatin organization and transcription. Intriguingly, we also identified a cluster of hits comprising metabolic enzymes of the polyamine synthesis pathway. We demonstrate that polyamines, which have long been used as buffer additives to support ribosome assembly in vitro, are required for 60S maturation in living cells. Perturbation of polyamine metabolism results in early defects in 60S but not 40S subunit maturation. Collectively, our data reveal a novel function for polyamines in living cells and provide a rich source for future studies on ribosome synthesis.

DNA replication occurring in S-phase is critical for the maintenance of the cell fate from one generation to the next, and requires the duplication of epigenetic information. The integrity of the epigenome is, in part, insured by the recycling of parental histones and de novo deposition of newly synthesized histones. While the histone variants have revealed important functions in epigenetic regulations, the deposition in chromatin during S-phase of newly synthesized histone variants remains unclear. The identification of histone variants of H3 and unique features of Physarum polycephalum provides a powerful system for investigating de novo deposition of newly synthesized histones by tracking the incorporation of exogenous histones within cells. The analyses revealed that the rate of deposition of H3.1 and H3.3 is anticorrelated as S-phase progresses, H3.3 is predominately produced and utilized in early S and dropped throughout S-phase, while H3.1 behaved in the opposite way. Disturbing the expression of H3 variants by siRNAs revealed mutual compensation of histone transcripts. Interestingly, the incorporation of pre-formed constrained histone complexes showed that tetramers of H3/H4 are more efficiently utilized by the cell than dimers. These results support the model whereby the histone variant distribution is established upon replication and new histone deposition.

Patient-derived tumor organoids (PDOs) have emerged as a reliable in vitro model for drug discovery. However, RNA sequencing-based analysis of PDOs treated with drugs has not been realized in a high-throughput format due to the limited quantity of organoids. Here, we translated a newly developed pooled RNA-seq methodology onto a superhydrophobic microwell array chip to realize an assay of genome-wide RNA output unified with phenotypic data (Grouped-seq). Over 10-fold reduction of sample and reagent consumption together with a new ligation-based barcode synthesis method lowers the cost to ∼$2 per RNA-seq sample. Patient-derived colorectal cancer (CRC) organoids with a number of 10 organoids per microwell were treated with four anti-CRC drugs across eight doses and analyzed by the Grouped-seq. Using a phenotype-assisted pathway enrichment analysis (PAPEA) method, the mechanism of actions of the drugs were correctly derived, illustrating the great potential of Grouped-seq for pharmacological screening with tumor organoids.

Comment in 10.1093/nar/gkac105.

Parasitic helminths infecting humans are highly prevalent infecting ∼2 billion people worldwide, causing inflammatory responses, malnutrition and anemia that are the primary cause of morbidity. In addition, helminth infections of cattle have a significant economic impact on livestock production, milk yield and fertility. The etiological agents of helminth infections are mainly Nematodes (roundworms) and Platyhelminths (flatworms). G-quadruplexes (G4) are unusual nucleic acid structures formed by G-rich sequences that can be recognized by specific G4 ligands. Here we used the G4Hunter Web Tool to identify and compare potential G4 sequences (PQS) in the nuclear and mitochondrial genomes of various helminths to identify G4 ligand targets. PQS are nonrandomly distributed in these genomes and often located in the proximity of genes. Unexpectedly, a Nematode, Ascaris lumbricoides, was found to be highly enriched in stable PQS. This species can tolerate high-stability G4 structures, which are not counter selected at all, in stark contrast to most other species. We experimentally confirmed G4 formation for sequences found in four different parasitic helminths. Small molecules able to selectively recognize G4 were found to bind to Schistosoma mansoni G4 motifs. Two of these ligands demonstrated potent activity both against larval and adult stages of this parasite.

Comment in 10.1093/nar/gkab1027.

The autosomal recessive genome instability disorder Ataxia-telangiectasia, caused by mutations in ATM kinase, is characterized by the progressive loss of cerebellar neurons. We find that DNA damage associated with ATM loss results in dysfunctional behaviour of human microglia, immune cells of the central nervous system. Microglial dysfunction is mediated by the pro-inflammatory RELB/p52 non-canonical NF-κB transcriptional pathway and leads to excessive phagocytic clearance of neuronal material. Activation of the RELB/p52 pathway in ATM-deficient microglia is driven by persistent DNA damage and is dependent on the NIK kinase. Activation of non-canonical NF-κB signalling is also observed in cerebellar microglia of individuals with Ataxia-telangiectasia. These results provide insights into the underlying mechanisms of aberrant microglial behaviour in ATM deficiency, potentially contributing to neurodegeneration in Ataxia-telangiectasia.

CMTR1 (cap methyltransferase 1) catalyses methylation of the first transcribed nucleotide of RNAPII transcripts (N1 2'-O-Me), creating part of the mammalian RNA cap structure. In addition to marking RNA as self, N1 2'-O-Me has ill-defined roles in RNA expression and translation. Here, we investigated the gene specificity of CMTR1 and its impact on RNA expression in embryonic stem cells. Using chromatin immunoprecipitation, CMTR1 was found to bind to transcription start sites (TSS) correlating with RNAPII levels, predominantly binding at histone genes and ribosomal protein (RP) genes. Repression of CMTR1 expression resulted in repression of RNAPII binding at the TSS and repression of RNA expression, particularly of histone and RP genes. In correlation with regulation of histones and RP genes, CMTR1 repression resulted in repression of translation and induction of DNA replication stress and damage. Indicating a direct role for CMTR1 in transcription, addition of recombinant CMTR1 to purified nuclei increased transcription of the histone and RP genes. CMTR1 was found to be upregulated during neural differentiation and there was an enhanced requirement for CMTR1 for gene expression and proliferation during this process. We highlight the distinct roles of the cap methyltransferases RNMT and CMTR1 in target gene expression and differentiation.

Serine integrases are emerging as one of the most powerful biological tools for synthetic biology. They have been widely used across genome engineering and genetic circuit design. However, developing serine integrase-based tools for directly/precisely manipulating synthetic biobricks is still missing. Here, we report SYMBIOSIS, a versatile method that can robustly manipulate DNA parts in vivo and in vitro. First, we propose a 'keys match locks' model to demonstrate that three orthogonal serine integrases are able to irreversibly and stably switch on seven synthetic biobricks with high accuracy in vivo. Then, we demonstrate that purified integrases can facilitate the assembly of 'donor' and 'acceptor' plasmids in vitro to construct composite plasmids. Finally, we use SYMBIOSIS to assemble different chromoprotein genes and create novel colored Escherichia coli. We anticipate that our SYMBIOSIS strategy will accelerate synthetic biobrick manipulation, genetic circuit design and multiple plasmid assembly for synthetic biology with broad potential applications.

Cyclic guanosine monophosphate-adenosine monophosphate synthase (cGAS) is activated in cells with defective DNA damage repair and signaling (DDR) factors, but a direct role for DDR factors in regulating cGAS activation in response to micronuclear DNA is still poorly understood. Here, we provide novel evidence that Nijmegen breakage syndrome 1 (NBS1) protein, a well-studied DNA double-strand break (DSB) sensor-in coordination with Ataxia Telangiectasia Mutated (ATM), a protein kinase, and Carboxy-terminal binding protein 1 interacting protein (CtIP), a DNA end resection factor-functions as an upstream regulator that prevents cGAS from binding micronuclear DNA. When NBS1 binds to micronuclear DNA via its fork-head-associated domain, it recruits CtIP and ATM via its N- and C-terminal domains, respectively. Subsequently, ATM stabilizes NBS1's interaction with micronuclear DNA, and CtIP converts DSB ends into single-strand DNA ends; these two key events prevent cGAS from binding micronuclear DNA. Additionally, by using a cGAS tripartite system, we show that cells lacking NBS1 not only recruit cGAS to a major fraction of micronuclear DNA but also activate cGAS in response to these micronuclear DNA. Collectively, our results underscore how NBS1 and its binding partners prevent cGAS from binding micronuclear DNA, in addition to their classical functions in DDR signaling.

Regulated transcription termination provides an efficient and responsive means to control gene expression. In bacteria, rho-independent termination occurs through the formation of an intrinsic RNA terminator loop, which disrupts the RNA polymerase elongation complex, resulting in its dissociation from the DNA template. Bacteria have a number of pathways for overriding termination, one of which is the formation of mutually exclusive RNA motifs. ANTAR domains are a class of antiterminator that bind and stabilize dual hexaloop RNA motifs within the nascent RNA chain to prevent terminator loop formation. We have determined the structures of the dimeric ANTAR domain protein EutV, from Enterococcus faecialis, in the absence of and in complex with the dual hexaloop RNA target. The structures illustrate conformational changes that occur upon RNA binding and reveal that the molecular interactions between the ANTAR domains and RNA are restricted to a single hexaloop of the motif. An ANTAR domain dimer must contact each hexaloop of the dual hexaloop motif individually to prevent termination in eubacteria. Our findings thereby redefine the minimal ANTAR domain binding motif to a single hexaloop and revise the current model for ANTAR-mediated antitermination. These insights will inform and facilitate the discovery of novel ANTAR domain RNA targets.

The use of complex biological molecules to solve computational problems is an emerging field at the interface between biology and computer science. There are two main categories in which biological molecules, especially DNA, are investigated as alternatives to silicon-based computer technologies. One is to use DNA as a storage medium, and the other is to use DNA for computing. Both strategies come with certain constraints. In the current study, we present a novel approach derived from chaos game representation for DNA to generate DNA code words that fulfill user-defined constraints, namely GC content, homopolymers, and undesired motifs, and thus, can be used to build codes for reliable DNA storage systems.

Multi-omics integration is key to fully understand complex biological processes in an holistic manner. Furthermore, multi-omics combined with new longitudinal experimental design can unreveal dynamic relationships between omics layers and identify key players or interactions in system development or complex phenotypes. However, integration methods have to address various experimental designs and do not guarantee interpretable biological results. The new challenge of multi-omics integration is to solve interpretation and unlock the hidden knowledge within the multi-omics data. In this paper, we go beyond integration and propose a generic approach to face the interpretation problem. From multi-omics longitudinal data, this approach builds and explores hybrid multi-omics networks composed of both inferred and known relationships within and between omics layers. With smart node labelling and propagation analysis, this approach predicts regulation mechanisms and multi-omics functional modules. We applied the method on 3 case studies with various multi-omics designs and identified new multi-layer interactions involved in key biological functions that could not be revealed with single omics analysis. Moreover, we highlighted interplay in the kinetics that could help identify novel biological mechanisms. This method is available as an R package netOmics to readily suit any application.

The function and identity of a cell is shaped by transcription factors controlling transcriptional networks, and further shaped by RNA binding proteins controlling post-transcriptional networks. To overcome limitations inherent to analysis of sparse single-cell post-transcriptional data, we leverage the invariant Caenorhabditis elegans cell lineage, isolating thousands of identical neuron types from thousands of isogenic individuals. The resulting deep transcriptomes facilitate splicing network analysis due to increased sequencing depth and uniformity. We focus on mechanosensory touch-neuron splicing regulated by MEC-8/RBPMS. We identify a small MEC-8-regulated network, where MEC-8 establishes touch-neuron isoforms differing from default isoforms found in other cells. MEC-8 establishes the canonical long mec-2/Stomatin isoform in touch neurons, but surprisingly the non-canonical short isoform predominates in other neurons, including olfactory neurons, and mec-2 is required for olfaction. Forced endogenous isoform-specific expression reveals that the short isoform functions in olfaction but not mechanosensation. The long isoform is functional in both processes. Remarkably, restoring the long isoform completely rescues mec-8 mutant mechanosensation, indicating a single MEC-8 touch-neuron target is phenotypically relevant. Within the long isoform we identify a cassette exon further diversifying mec-2 into long/extra-long isoforms. Neither is sufficient for mechanosensation. Both are simultaneously required, likely functioning as heteromers to mediate mechanosensation.

Comment in 10.1093/nar/gkac108.

The tumour suppressor SLX4 plays multiple roles in the maintenance of genome stability, acting as a scaffold for structure-specific endonucleases and other DNA repair proteins. It directly interacts with the mismatch repair (MMR) protein MSH2 but the significance of this interaction remained unknown until recent findings showing that MutSβ (MSH2-MSH3) stimulates in vitro the SLX4-dependent Holliday junction resolvase activity. Here, we characterize the mode of interaction between SLX4 and MSH2, which relies on an MSH2-interacting peptide (SHIP box) that drives interaction of SLX4 with both MutSβ and MutSα (MSH2-MSH6). While we show that this MSH2 binding domain is dispensable for the well-established role of SLX4 in interstrand crosslink repair, we find that it mediates inhibition of MutSα-dependent MMR by SLX4, unravelling an unanticipated function of SLX4.

The histone acetyltransferase p300 (also known as KAT3B) is a general transcriptional coactivator that introduces the H3K27ac mark on enhancers triggering their activation and gene transcription. Genome-wide screenings demonstrated that a large fraction of long non-coding RNAs (lncRNAs) plays a role in cellular processes and organ development although the underlying molecular mechanisms remain largely unclear (1,2). We found 122 lncRNAs that interacts directly with p300. In depth analysis of one of these, lncSmad7, is required to maintain ESC self-renewal and it interacts to the C-terminal domain of p300. lncSmad7 also contains predicted RNA-DNA Hoogsteen forming base pairing. Combined Chromatin Isolation by RNA precipitation followed by sequencing (ChIRP-seq) together with CRISPR/Cas9 mutagenesis of the target sites demonstrate that lncSmad7 binds and recruits p300 to enhancers in trans, to trigger enhancer acetylation and transcriptional activation of its target genes. Thus, these results unveil a new mechanism by which p300 is recruited to the genome.

Non-coding variants have long been recognized as important contributors to common disease risks, but with the expansion of clinical whole genome sequencing, examples of rare, high-impact non-coding variants are also accumulating. Despite recent advances in the study of regulatory elements and the availability of specialized data collections, the systematic annotation of non-coding variants from genome sequencing remains challenging. Here, we propose a new framework for the prioritization of non-coding regulatory variants that integrates information about regulatory regions with prediction scores and HPO-based prioritization. Firstly, we created a comprehensive collection of annotations for regulatory regions including a database of 2.4 million regulatory elements (GREEN-DB) annotated with controlled gene(s), tissue(s) and associated phenotype(s) where available. Secondly, we calculated a variation constraint metric and showed that constrained regulatory regions associate with disease-associated genes and essential genes from mouse knock-outs. Thirdly, we compared 19 non-coding impact prediction scores providing suggestions for variant prioritization. Finally, we developed a VCF annotation tool (GREEN-VARAN) that can integrate all these elements to annotate variants for their potential regulatory impact. In our evaluation, we show that GREEN-DB can capture previously published disease-associated non-coding variants as well as identify additional candidate disease genes in trio analyses.

The mammalian high mobility group protein AT-hook 2 (HMGA2) houses three motifs that preferentially bind short stretches of AT-rich DNA regions. These DNA binding motifs, known as 'AT-hooks', are traditionally characterized as being unstructured. Upon binding to AT-rich DNA, they form ordered assemblies. It is this disordered-to-ordered transition that has implicated HMGA2 as a protein actively involved in many biological processes, with abnormal HMGA expression linked to a variety of health problems including diabetes, obesity, and oncogenesis. In the current work, the solution binding dynamics of the three 'AT-hook' peptides (ATHPs) with AT-rich DNA hairpin substrates were studied using DNA UV melting studies, fluorescence spectroscopy, native ion mobility spectrometry-mass spectrometry (IMS-MS), solution isothermal titration calorimetry (ITC) and molecular modeling. Results showed that the ATHPs bind to the DNA to form a single, 1:1 and 2:1, 'key-locked' conformational ensemble. The molecular models showed that 1:1 and 2:1 complex formation is driven by the capacity of the ATHPs to bind to the minor and major grooves of the AT-rich DNA oligomers. Complementary solution ITC results confirmed that the 2:1 stoichiometry of ATHP: DNA is originated under native conditions in solution.

Biogenesis of spliceosomal small nuclear ribonucleoproteins (snRNPs) and their recycling after splicing require numerous assembly/recycling factors whose modes of action are often poorly understood. The intrinsically disordered TSSC4 protein has been identified as a nuclear-localized U5 snRNP and U4/U6-U5 tri-snRNP assembly/recycling factor, but how TSSC4's intrinsic disorder supports TSSC4 functions remains unknown. Using diverse interaction assays and cryogenic electron microscopy-based structural analysis, we show that TSSC4 employs four conserved, non-contiguous regions to bind the PRPF8 Jab1/MPN domain and the SNRNP200 helicase at functionally important sites. It thereby inhibits SNRNP200 helicase activity, spatially aligns the proteins, coordinates formation of a U5 sub-module and transiently blocks premature interaction of SNRNP200 with at least three other spliceosomal factors. Guided by the structure, we designed a TSSC4 variant that lacks stable binding to the PRPF8 Jab1/MPN domain or SNRNP200 in vitro. Comparative immunoprecipitation/mass spectrometry from HEK293 nuclear extract revealed distinct interaction profiles of wild type TSSC4 and the variant deficient in PRPF8/SNRNP200 binding with snRNP proteins, other spliceosomal proteins as well as snRNP assembly/recycling factors and chaperones. Our findings elucidate molecular strategies employed by an intrinsically disordered protein to promote snRNP assembly, and suggest multiple TSSC4-dependent stages during snRNP assembly/recycling.

The Isw1b chromatin-remodeling complex is specifically recruited to gene bodies to help retain pre-existing histones during transcription by RNA polymerase II. Recruitment is dependent on H3K36 methylation and the Isw1b subunit Ioc4, which contains an N-terminal PWWP domain. Here, we present the crystal structure of the Ioc4-PWWP domain, including a detailed functional characterization of the domain on its own as well as in the context of full-length Ioc4 and the Isw1b remodeler. The Ioc4-PWWP domain preferentially binds H3K36me3-containing nucleosomes. Its ability to bind DNA is required for nucleosome binding. It is also furthered by the unique insertion motif present in Ioc4-PWWP. The ability to bind H3K36me3 and DNA promotes the interaction of full-length Ioc4 with nucleosomes in vitro and they are necessary for its recruitment to gene bodies in vivo. Furthermore, a fully functional Ioc4-PWWP domain promotes efficient remodeling by Isw1b and the maintenance of ordered chromatin in vivo, thereby preventing the production of non-coding RNAs.

In ribosomal translation, peptidyl transfer occurs between P-site peptidyl-tRNA and A-site aminoacyl-tRNA, followed by translocation of the resulting P-site deacylated-tRNA and A-site peptidyl-tRNA to E and P site, respectively, mediated by EF-G. Here, we report that mistranslocation of P-site peptidyl-tRNA and A-site aminoacyl-tRNA toward E and A site occurs when high concentration of EF-G triggers the migration of two tRNAs prior to completion of peptidyl transfer. Consecutive incorporation of less reactive amino acids, such as Pro and d-Ala, makes peptidyl transfer inefficient and thus induces the mistranslocation event. Consequently, the E-site peptidyl-tRNA drops off from ribosome to give a truncated peptide lacking the C-terminal region. The P-site aminoacyl-tRNA allows for reinitiation of translation upon accommodation of a new aminoacyl-tRNA at A site, leading to synthesis of a truncated peptide lacking the N-terminal region, which we call the 'reinitiated peptide'. We also revealed that such a drop-off-reinitiation event can be alleviated by EF-P that promotes peptidyl transfer of Pro. Moreover, this event takes place both in vitro and in cell, showing that reinitiated peptides during protein synthesis could be accumulated in this pathway in cells.

Some proteins, like the lac repressor (LacI), mediate long-range loops that alter DNA topology and create torsional barriers. During transcription, RNA polymerase generates supercoiling that may facilitate passage through such barriers. We monitored E. coli RNA polymerase progress along templates in conditions that prevented, or favored, 400 bp LacI-mediated DNA looping. Tethered particle motion measurements revealed that RNA polymerase paused longer at unlooped LacI obstacles or those barring entry to a loop than those barring exit from the loop. Enhanced dissociation of a LacI roadblock by the positive supercoiling generated ahead of a transcribing RNA polymerase within a torsion-constrained DNA loop may be responsible for this reduction in pause time. In support of this idea, RNA polymerase transcribed 6-fold more slowly through looped DNA and paused at LacI obstacles for 66% less time on positively supercoiled compared to relaxed templates, especially under increased tension (torque). Positive supercoiling propagating ahead of polymerase facilitated elongation along topologically complex, protein-coated templates.

In leiomyosarcoma class IIa HDACs (histone deacetylases) bind MEF2 and convert these transcription factors into repressors to sustain proliferation. Disruption of this complex with small molecules should antagonize cancer growth. NKL54, a PAOA (pimeloylanilide o-aminoanilide) derivative, binds a hydrophobic groove of MEF2, which is used as a docking site by class IIa HDACs. However, NKL54 could also act as HDAC inhibitor (HDACI). Therefore, it is unclear which activity is predominant. Here, we show that NKL54 and similar derivatives are unable to release MEF2 from binding to class IIa HDACs. Comparative transcriptomic analysis classifies these molecules as HDACIs strongly related to SAHA/vorinostat. Low expressed genes are upregulated by HDACIs, while abundant genes are repressed. This transcriptional resetting correlates with a reorganization of H3K27 acetylation around the transcription start site (TSS). Among the upregulated genes there are several BH3-only family members, thus explaining the induction of apoptosis. Moreover, NKL54 triggers the upregulation of MEF2 and the downregulation of class IIa HDACs. NKL54 also increases the binding of MEF2D to promoters of genes that are upregulated after treatment. In summary, although NKL54 cannot outcompete MEF2 from binding to class IIa HDACs, it supports MEF2-dependent transcription through several actions, including potentiation of chromatin binding.

Reverse transcriptases (RTs) are found in different systems including group II introns, Diversity Generating Retroelements (DGRs), retrons, CRISPR-Cas systems, and Abortive Infection (Abi) systems in prokaryotes. Different classes of RTs can play different roles, such as template switching and mobility in group II introns, spacer acquisition in CRISPR-Cas systems, mutagenic retrohoming in DGRs, programmed cell suicide in Abi systems, and recently discovered phage defense in retrons. While some classes of RTs have been studied extensively, others remain to be characterized. There is a lack of computational tools for identifying and characterizing various classes of RTs. In this study, we built a tool (called myRT) for identification and classification of prokaryotic RTs. In addition, our tool provides information about the genomic neighborhood of each RT, providing potential functional clues. We applied our tool to predict RTs in all complete and draft bacterial genomes, and created a collection that can be used for exploration of putative RTs and their associated protein domains. Application of myRT to metagenomes showed that gut metagenomes encode proportionally more RTs related to DGRs, outnumbering retron-related RTs, as compared to the collection of reference genomes. MyRT is both available as a standalone software (https://github.com/mgtools/myRT) and also through a website (https://omics.informatics.indiana.edu/myRT/).

CUGBP Elav-like family member 1 (CELF1), an RNA-binding protein (RBP), plays important roles in the pathogenesis of diseases such as myotonic dystrophy, liver fibrosis and cancers. However, targeting CELF1 is still a challenge, as RBPs are considered largely undruggable. Here, we discovered that compound 27 disrupted CELF1-RNA binding via structure-based virtual screening and biochemical assays. Compound 27 binds directly to CELF1 and competes with RNA for binding to CELF1. Compound 27 promotes IFN-γ secretion and suppresses TGF-β1-induced hepatic stellate cell (HSC) activation by inhibiting CELF1-mediated IFN-γ mRNA decay. In vivo, compound 27 attenuates CCl4-induced murine liver fibrosis. Furthermore, the structure-activity relationship analysis was performed and compound 841, a derivative of compound 27, was identified as a selective CELF1 inhibitor. In conclusion, targeting CELF1 RNA-binding activity with small molecules was achieved, which provides a novel strategy for treating liver fibrosis and other CELF1-mediated diseases.

Selection of the appropriate DNA double-strand break (DSB) repair pathway is decisive for genetic stability. It is proposed to act according to two steps: 1-canonical nonhomologous end-joining (C-NHEJ) versus resection that generates single-stranded DNA (ssDNA) stretches; 2-on ssDNA, gene conversion (GC) versus nonconservative single-strand annealing (SSA) or alternative end-joining (A-EJ). Here, we addressed the mechanisms by which RAD51 regulates this second step, preventing nonconservative repair in human cells. Silencing RAD51 or BRCA2 stimulated both SSA and A-EJ, but not C-NHEJ, validating the two-step model. Three different RAD51 dominant-negative forms (DN-RAD51s) repressed GC and stimulated SSA/A-EJ. However, a fourth DN-RAD51 repressed SSA/A-EJ, although it efficiently represses GC. In living cells, the three DN-RAD51s that stimulate SSA/A-EJ failed to load efficiently onto damaged chromatin and inhibited the binding of endogenous RAD51, while the fourth DN-RAD51, which inhibits SSA/A-EJ, efficiently loads on damaged chromatin. Therefore, the binding of RAD51 to DNA, rather than its ability to promote GC, is required for SSA/A-EJ inhibition by RAD51. We showed that RAD51 did not limit resection of endonuclease-induced DSBs, but prevented spontaneous and RAD52-induced annealing of complementary ssDNA in vitro. Therefore, RAD51 controls the selection of the DSB repair pathway, protecting genome integrity from nonconservative DSB repair through ssDNA occupancy, independently of the promotion of CG.