DNA G4-structures from human c-MYC promoter and telomere are considered as important drug targets; however, the developing of small-molecule-based fluorescent binding ligands that are highly selective in targeting these G4-structures over other types of nucleic acids is challenging. We herein report a new approach of designing small molecules based on a non-selective thiazole orange scaffold to provide two-directional and multi-site interactions with flanking residues and loops of the G4-motif for better selectivity. The ligands are designed to establish multi-site interactions in the G4-binding pocket. This structural feature may render the molecules higher selectivity toward c-MYC G4s than other structures. The ligand-G4 interaction studied with 1H NMR may suggest a stacking interaction with the terminal G-tetrad. Moreover, the intracellular co-localization study with BG4 and cellular competition experiments with BRACO-19 may suggest that the binding targets of the ligands in cells are most probably G4-structures. Furthermore, the ligands that either preferentially bind to c-MYC promoter or telomeric G4s are able to downregulate markedly the c-MYC and hTERT gene expression in MCF-7 cells, and induce senescence and DNA damage to cancer cells. The in vivo antitumor activity of the ligands in MCF-7 tumor-bearing mice is also demonstrated.

MicroRNAs silence mRNAs by guiding the RISC complex. RISC assembly occurs following cleavage of pre-miRNAs by Dicer, assisted by TRBP or PACT, and the transfer of miRNAs to AGO proteins. The R2TP complex is an HSP90 co-chaperone involved in the assembly of ribonucleoprotein particles. Here, we show that the R2TP component RPAP3 binds TRBP but not PACT. The RPAP3-TPR1 domain interacts with the TRBP-dsRBD3, and the 1.5 Å resolution crystal structure of this complex identifies key residues involved in the interaction. Remarkably, binding of TRBP to RPAP3 or Dicer is mutually exclusive. Additionally, we found that AGO(1/2), TRBP and Dicer are all sensitive to HSP90 inhibition, and that TRBP sensitivity is increased in the absence of RPAP3. Finally, RPAP3 seems to impede miRNA activity, raising the possibility that the R2TP chaperone might sequester TRBP to regulate the miRNA pathway.

In rapidly growing cells, with recombinational DNA repair required often and a new replication fork passing every 20 min, the pace of RecA-mediated DNA strand exchange is potentially much too slow for bacterial DNA metabolism. The enigmatic RadD protein, a putative SF2 family helicase, exhibits no independent helicase activity on branched DNAs. Instead, RadD greatly accelerates RecA-mediated DNA strand exchange, functioning only when RecA protein is present. The RadD reaction requires the RadD ATPase activity, does not require an interaction with SSB, and may disassemble RecA filaments as it functions. We present RadD as a new class of enzyme, an accessory protein that accelerates DNA strand exchange, possibly with a helicase-like action, in a reaction that is entirely RecA-dependent. RadD is thus a DNA strand exchange (recombination) synergist whose primary function is to coordinate closely with and accelerate the DNA strand exchange reactions promoted by the RecA recombinase. Multiple observations indicate a uniquely close coordination of RadD with RecA function.

Erratum in Nucleic Acids Res. 2022 May 6;50(8):4799.

Bacterial exonuclease III (ExoIII) is a multifunctional enzyme that uses a single active site to perform two conspicuous activities: (i) apurinic/apyrimidinic (AP)-endonuclease and (ii) 3'→5' exonuclease activities. The AP endonuclease activity results in AP site incision, while the exonuclease activity results in the continuous excision of 3' terminal nucleobases to generate a partial duplex for recruiting the downstream DNA polymerase during the base excision repair process (BER). The key determinants of functional selection between the two activities are poorly understood. Here, we use a series of mutational analyses and single-molecule imaging to unravel the pivotal rules governing these endo- and exonuclease activities at the single amino acid level. An aromatic residue, either W212 or F213, recognizes AP sites to allow for the AP endonuclease activity, and the F213 residue also participates in the stabilization of the melted state of the 3' terminal nucleobases, leading to the catalytically competent state that activates the 3'→5' exonuclease activity. During exonucleolytic cleavage, the DNA substrate must be maintained as a B-form helix through a series of phosphate-stabilizing residues (R90, Y109, K121 and N153). Our work decouples the AP endonuclease and exonuclease activities of ExoIII and provides insights into how this multifunctional enzyme controls each function at the amino acid level.

A plethora of modified nucleotides extends the chemical and conformational space for natural occurring RNAs. tRNAs constitute the class of RNAs with the highest modification rate. The extensive modification modulates their overall stability, the fidelity and efficiency of translation. However, the impact of nucleotide modifications on the local structural dynamics is not well characterized. Here we show that the incorporation of the modified nucleotides in tRNAfMet from Escherichia coli leads to an increase in the local conformational dynamics, ultimately resulting in the stabilization of the overall tertiary structure. Through analysis of the local dynamics by NMR spectroscopic methods we find that, although the overall thermal stability of the tRNA is higher for the modified molecule, the conformational fluctuations on the local level are increased in comparison to an unmodified tRNA. In consequence, the melting of individual base pairs in the unmodified tRNA is determined by high entropic penalties compared to the modified. Further, we find that the modifications lead to a stabilization of long-range interactions harmonizing the stability of the tRNA's secondary and tertiary structure. Our results demonstrate that the increase in chemical space through introduction of modifications enables the population of otherwise inaccessible conformational substates.

In oligodendrocytes of the vertebrate central nervous system a complex network of transcriptional regulators is required to ensure correct and timely myelination of neuronal axons. Here we identify Zfp276, the only mammalian ZAD-domain containing zinc finger protein, as a transcriptional regulator of oligodendrocyte differentiation and central myelination downstream of Sox10. In the central nervous system, Zfp276 is exclusively expressed in mature oligodendrocytes. Oligodendroglial deletion of Zfp276 led to strongly reduced expression of myelin genes in the early postnatal mouse spinal cord. Retroviral overexpression of Zfp276 in cultured oligodendrocyte precursor cells induced precocious expression of maturation markers and myelin genes, further supporting its role in oligodendroglial differentiation. On the molecular level, Zfp276 directly binds to and represses Sox10-dependent gene regulatory regions of immaturity factors and functionally interacts with the transcriptional repressor Zeb2 to enable fast transition of oligodendrocytes to the myelinating stage.

N 6-Threonylcarbamoyladenosine (t6A) is a universal and pivotal tRNA modification. KEOPS in eukaryotes participates in its biogenesis, whose mutations are connected with Galloway-Mowat syndrome. However, the tRNA substrate selection mechanism by KEOPS and t6A modification function in mammalian cells remain unclear. Here, we confirmed that all ANN-decoding human cytoplasmic tRNAs harbor a t6A moiety. Using t6A modification systems from various eukaryotes, we proposed the possible coevolution of position 33 of initiator tRNAMet and modification enzymes. The role of the universal CCA end in t6A biogenesis varied among species. However, all KEOPSs critically depended on C32 and two base pairs in the D-stem. Knockdown of the catalytic subunit OSGEP in HEK293T cells had no effect on the steady-state abundance of cytoplasmic tRNAs but selectively inhibited tRNAIle aminoacylation. Combined with in vitro aminoacylation assays, we revealed that t6A functions as a tRNAIle isoacceptor-specific positive determinant for human cytoplasmic isoleucyl-tRNA synthetase (IARS1). t6A deficiency had divergent effects on decoding efficiency at ANN codons and promoted +1 frameshifting. Altogether, our results shed light on the tRNA recognition mechanism, revealing both commonality and diversity in substrate recognition by eukaryotic KEOPSs, and elucidated the critical role of t6A in tRNAIle aminoacylation and codon decoding in human cells.

Trypanosoma brucei causes human African trypanosomiasis and sequentially expresses distinct VSGs, its major surface antigen, to achieve host immune evasion. VSGs are monoallelically expressed from subtelomeric loci, and telomere proteins regulate VSG monoallelic expression and VSG switching. T. brucei telomerase is essential for telomere maintenance, but no regulators of telomerase have been identified. T. brucei appears to lack OB fold-containing telomere-specific ssDNA binding factors that are critical for coordinating telomere G- and C-strand syntheses in higher eukaryotes. We identify POLIE as a telomere protein essential for telomere integrity. POLIE-depleted cells have more frequent VSG gene conversion-mediated VSG switching and an increased amount of telomeric circles (T-circles), indicating that POLIE suppresses DNA recombination at the telomere/subtelomere. POLIE-depletion elongates telomere 3' overhangs dramatically, indicating that POLIE is essential for coordinating DNA syntheses of the two telomere strands. POLIE depletion increases the level of telomerase-dependent telomere G-strand extension, identifying POLIE as the first T. brucei telomere protein that suppresses telomerase. Furthermore, depletion of POLIE results in an elevated telomeric C-circle level, suggesting that the telomere C-strand experiences replication stress and that POLIE may promote telomere C-strand synthesis. Therefore, T. brucei uses a novel mechanism to coordinate the telomere G- and C-strand DNA syntheses.

Although there are several pathways to ensure that proteins are folded properly in the cell, little is known about the molecular mechanisms regulating histone folding and proteostasis. In this work, we identified that chaperone-mediated autophagy (CMA) is the main pathway involved in the degradation of newly synthesized histones H3 and H4. This degradation is finely regulated by the interplay between HSC70 and tNASP, two histone interacting proteins. tNASP stabilizes histone H3 levels by blocking the direct transport of histone H3 into lysosomes. We further demonstrate that CMA degrades unfolded histone H3. Thus, we reveal that CMA is the main degradation pathway involved in the quality control of histone biogenesis, evidencing an additional mechanism in the intricate network of histone cellular proteostasis.

Thousands of RNA species display nonuniform distribution within cells. However, quantification of the spatial patterns adopted by individual RNAs remains difficult, in part by a lack of quantitative tools for subcellular transcriptome analysis. In this study, we describe an RNA proximity labeling method that facilitates the quantification of subcellular RNA populations with high spatial specificity. This method, termed Halo-seq, pairs a light-activatable, radical generating small molecule with highly efficient Click chemistry to efficiently label and purify spatially defined RNA samples. We compared Halo-seq with previously reported similar methods and found that Halo-seq displayed a higher efficiency of RNA labeling, indicating that it is well suited to the investigation of small, precisely localized RNA populations. We then used Halo-seq to quantify nuclear, nucleolar and cytoplasmic transcriptomes, characterize their dynamic nature following perturbation, and identify RNA sequence features associated with their composition. Specifically, we found that RNAs containing AU-rich elements are relatively enriched in the nucleus. This enrichment becomes stronger upon treatment with the nuclear export inhibitor leptomycin B, both expanding the role of HuR in RNA export and generating a comprehensive set of transcripts whose export from the nucleus depends on HuR.

Most RNA footprinting approaches that require ribonuclease cleavage generate RNA fragments bearing a phosphate or cyclic phosphate group at their 3' end. Unfortunately, current library preparation protocols rely only on a 3' hydroxyl group for adaptor ligation or poly-A tailing. Here, we developed circAID-p-seq, a PCR-free library preparation for selective 3' phospho-RNA sequencing. As a proof of concept, we applied circAID-p-seq to ribosome profiling, which is based on sequencing of RNA fragments protected by ribosomes after endonuclease digestion. CircAID-p-seq, combined with the dedicated computational pipeline circAidMe, facilitates accurate, fast and highly efficient sequencing of phospho-RNA fragments from eukaryotic cells and tissues. We used circAID-p-seq to portray ribosome occupancy in transcripts, providing a versatile and PCR-free strategy to possibly unravel any endogenous 3'-phospho RNA molecules.

The A-repeat region of the lncRNA Xist is critical for X inactivation and harbors several N6-methyladenosine (m6A) modifications. How the m6A modification affects the conformation of the conserved AUCG tetraloop hairpin of the A-repeats and how it can be recognized by the YTHDC1 reader protein is unknown. Here, we report the NMR solution structure of the (m6A)UCG hairpin, which reveals that the m6A base extends 5' stacking of the A-form helical stem, resembling the unmethylated AUCG tetraloop. A crystal structure of YTHDC1 bound to the (m6A)UCG tetraloop shows that the (m6A)UC nucleotides are recognized by the YTH domain of YTHDC1 in a single-stranded conformation. The m6A base inserts into the aromatic cage and the U and C bases interact with a flanking charged surface region, resembling the recognition of single-stranded m6A RNA ligands. Notably, NMR and fluorescence quenching experiments show that the binding requires local unfolding of the upper stem region of the (m6A)UCG hairpin. Our data show that m6A can be readily accommodated in hairpin loop regions, but recognition by YTH readers requires local unfolding of flanking stem regions. This suggests how m6A modifications may regulate lncRNA function by modulating RNA structure.

Cytosine methylation at the 5-carbon position is an essential DNA epigenetic mark in many eukaryotic organisms. Although countless structural and functional studies of cytosine methylation have been reported, our understanding of how it influences the nucleosome assembly, structure, and dynamics remains obscure. Here, we investigate the effects of cytosine methylation at CpG sites on nucleosome dynamics and stability. By applying long molecular dynamics simulations on several microsecond time scale, we generate extensive atomistic conformational ensembles of full nucleosomes. Our results reveal that methylation induces pronounced changes in geometry for both linker and nucleosomal DNA, leading to a more curved, under-twisted DNA, narrowing the adjacent minor grooves, and shifting the population equilibrium of sugar-phosphate backbone geometry. These DNA conformational changes are associated with a considerable enhancement of interactions between methylated DNA and the histone octamer, doubling the number of contacts at some key arginines. H2A and H3 tails play important roles in these interactions, especially for DNA methylated nucleosomes. This, in turn, prevents a spontaneous DNA unwrapping of 3-4 helical turns for the methylated nucleosome with truncated histone tails, otherwise observed in the unmethylated system on several microseconds time scale.

Histone 3 lysine 4 trimethylation (H3K4me3) is an epigenetic mark found at gene promoters and CpG islands. H3K4me3 is essential for mammalian development, yet mechanisms underlying its genomic targeting are poorly understood. H3K4me3 methyltransferases SETD1B and MLL2 (KMT2B) are essential for oogenesis. We investigated changes in H3K4me3 in Setd1b conditional knockout (cKO) oocytes using ultra-low input ChIP-seq, with comparisons to DNA methylation and gene expression analyses. H3K4me3 was redistributed in Setd1b cKO oocytes showing losses at active gene promoters associated with downregulated gene expression. Remarkably, many regions also gained H3K4me3, in particular those that were DNA hypomethylated, transcriptionally inactive and CpG-rich, which are hallmarks of MLL2 targets. Consequently, loss of SETD1B disrupts the balance between MLL2 and de novo DNA methyltransferases in determining the epigenetic landscape during oogenesis. Our work reveals two distinct, complementary mechanisms of genomic targeting of H3K4me3 in oogenesis, with SETD1B linked to gene expression and MLL2 to CpG content.

Data alignment is one of the first key steps in single cell analysis for integrating multiple datasets and performing joint analysis across studies. Data alignment is challenging in extremely large datasets, however, as the major of the current single cell data alignment methods are not computationally efficient. Here, we present VIPCCA, a computational framework based on non-linear canonical correlation analysis for effective and scalable single cell data alignment. VIPCCA leverages both deep learning for effective single cell data modeling and variational inference for scalable computation, thus enabling powerful data alignment across multiple samples, multiple data platforms, and multiple data types. VIPCCA is accurate for a range of alignment tasks including alignment between single cell RNAseq and ATACseq datasets and can easily accommodate millions of cells, thereby providing researchers unique opportunities to tackle challenges emerging from large-scale single-cell atlas.

The CRISPR-Cas9 genome editing tool is used to study genomic variants and gene knockouts, and can be combined with transcriptomic analyses to measure the effects of such alterations on gene expression. But how can one be sure that differential gene expression is due to a successful intended edit and not to an off-target event, without performing an often resource-demanding genome-wide sequencing of the edited cell or strain? To address this question we developed CRISPRroots: CRISPR-Cas9-mediated edits with accompanying RNA-seq data assessed for on-target and off-target sites. Our method combines Cas9 and guide RNA binding properties, gene expression changes, and sequence variants between edited and non-edited cells to discover potential off-targets. Applied on seven public datasets, CRISPRroots identified critical off-target candidates that were overlooked in all of the corresponding previous studies. CRISPRroots is available via https://rth.dk/resources/crispr.

Accurately quantifying gene and isoform expression changes is essential to understanding cell functions, differentiation and disease. Sequencing full-length native RNAs using long-read direct RNA sequencing (DRS) has the potential to overcome many limitations of short and long-read sequencing methods that require RNA fragmentation, cDNA synthesis or PCR. However, there are a lack of tools specifically designed for DRS and its ability to identify differential expression in complex organisms is poorly characterised. We developed NanoCount for fast, accurate transcript isoform quantification in DRS and demonstrate it outperforms similar methods. Using synthetic controls and human SH-SY5Y cell differentiation into neuron-like cells, we show that DRS accurately quantifies RNA expression and identifies differential expression of genes and isoforms. Differential expression of 231 genes, 333 isoforms, plus 27 isoform switches were detected between undifferentiated and differentiated SH-SY5Y cells and samples clustered by differentiation state at the gene and isoform level. Genes upregulated in neuron-like cells were associated with neurogenesis. NanoCount quantification of thousands of novel isoforms discovered with DRS likewise enabled identification of their differential expression. Our results demonstrate enhanced DRS isoform quantification with NanoCount and establish the ability of DRS to identify biologically relevant differential expression of genes and isoforms.

It is currently unknown how many RNA transcripts are able to induce degradation of microRNAs (miRNA) via the mechanism known as target-directed miRNA degradation (TDMD). We developed TDMDfinder, a computational pipeline that identifies 'high confidence' TDMD interactions in the Human and Mouse transcriptomes by combining sequence alignment and feature selection approaches. Our predictions suggested that TDMD is widespread, with potentially every miRNA controlled by endogenous targets. We experimentally tested 37 TDMDfinder predictions, of which 17 showed TDMD effects as measured by RT-qPCR and small RNA sequencing, linking the miR-17, miR-19, miR-30, miR-221, miR-26 and miR-23 families to novel endogenous TDMDs. In some cases, TDMD was found to affect different members of the same miRNA family selectively. Features like complementarity to the miRNA 3' region, bulge size and hybridization energy appeared to be the main factors determining sensitivity. Computational analyses performed using the multiomic TCGA platform substantiated the involvement of many TDMD transcripts in human cancer and highlighted 36 highly significant interactions, suggesting TDMD as a new potential oncogenic mechanism. In conclusion, TDMDfinder provides the first inventory of bona fide human and mouse TDMDs. Available as a free webtool, TDMDfinder allows users to search for any TDMD interaction of interest by customizing its selection criteria.

Unlike their natural counterparts, synthetic genetic circuits are usually fragile in the face of environmental perturbations and genetic mutations. Several theoretical robust genetic circuits have been designed, but their performance under real-world conditions has not yet been carefully evaluated. Here, we designed and synthesized a new robust perfect adaptation circuit composed of two-node negative feedback coupling with linear positive feedback on the buffer node. As a key feature, the linear positive feedback was fine-tuned to evaluate its necessity. We found that the desired function was robustly achieved when genetic parameters were varied by systematically perturbing all interacting parts within the topology, and the necessity of the completeness of the topological structures was evaluated by destroying key circuit features. Furthermore, different environmental perturbances were imposed onto the circuit by changing growth rates, carbon metabolic strategies and even chassis cells, and the designed perfect adaptation function was still achieved under all conditions. The successful design of a robust perfect adaptation circuit indicated that the top-down design strategy is capable of predictably guiding bottom-up engineering for robust genetic circuits. This robust adaptation circuit could be integrated as a motif into more complex circuits to robustly implement more sophisticated and critical biological functions.

Transposable elements are an abundant source of transcription factor binding sites, and favorable genomic integration may lead to their recruitment by the host genome for gene regulatory functions. However, it is unclear how frequent co-option of transposable elements as regulatory elements is, to which regulatory programs they contribute and how they compare to regulatory elements devoid of transposable elements. Here, we report a transcription initiation-centric, in-depth characterization of the transposon-derived regulatory landscape of mouse embryonic stem cells. We demonstrate that a substantial number of transposable element insertions, in particular endogenous retroviral elements, are associated with open chromatin regions that are divergently transcribed into unstable RNAs in a cell-type specific manner, and that these elements contribute to a sizable proportion of active enhancers and gene promoters. We further show that transposon subfamilies contribute differently and distinctly to the pluripotency regulatory program through their repertoires of transcription factor binding site sequences, shedding light on the formation of regulatory programs and the origins of regulatory elements.

Transcription activator-like effectors (TALEs) are bacterial proteins with a programmable DNA-binding domain, which turned them into exceptional tools for biotechnology. TALEs contain a central array of consecutive 34 amino acid long repeats to bind DNA in a simple one-repeat-to-one-nucleotide manner. However, a few naturally occurring aberrant repeat variants break this strict binding mechanism, allowing for the recognition of an additional sequence with a -1 nucleotide frameshift. The limits and implications of this extended TALE binding mode are largely unexplored. Here, we analyse the complete diversity of natural and artificially engineered aberrant repeats for their impact on the DNA binding of TALEs. Surprisingly, TALEs with several aberrant repeats can loop out multiple repeats simultaneously without losing DNA-binding capacity. We also characterized members of the only natural TALE class harbouring two aberrant repeats and confirmed that their target is the major virulence factor OsSWEET13 from rice. In an aberrant TALE repeat, the position and nature of the amino acid sequence strongly influence its function. We explored the tolerance of TALE repeats towards alterations further and demonstrate that inserts as large as GFP can be tolerated without disrupting DNA binding. This illustrates the extraordinary DNA-binding capacity of TALEs and opens new uses in biotechnology.

The ribosome is a fundamental biomolecular complex that synthesizes proteins in cells. Nascent proteins emerge from the ribosome through a tunnel, where they may interact with the tunnel walls or small molecules such as antibiotics. These interactions can cause translational arrest with notable physiological consequences. Here, we studied the arrest caused by the regulatory peptide VemP, which is known to form α-helices inside the ribosome tunnel near the peptidyl transferase center under specific conditions. We used all-atom molecular dynamics simulations of the entire ribosome and circular dichroism spectroscopy to study the driving forces of helix formation and how VemP causes the translational arrest. To that aim, we compared VemP dynamics in the ribosome tunnel with its dynamics in solution. We show that the VemP peptide has a low helical propensity in water and that the propensity is higher in mixtures of water and trifluorethanol. We propose that helix formation within the ribosome is driven by the interactions of VemP with the tunnel and that a part of VemP acts as an anchor. This anchor might slow down VemP progression through the tunnel enabling α-helix formation, which causes the elongation arrest.

The nuclease ARTEMIS and the DNA-dependent protein kinase catalytic subunit (DNA-PKcs) are involved in the repair of physiological and pathogenic DNA double strand breaks. Both proteins are indispensable for the hairpin-opening activity in V(D)J recombination and therefore essential for the adaptive immune response. ARTEMIS and DNA-PKcs interact, however experimental evidence for in vivo significance is missing. We demonstrate that mutations abolishing this protein-protein interaction affect nuclease function. In DNA-PKcs, mutation L3062R impairs the physical interaction with ARTEMIS and was previously identified as pathogenic variant, resulting in radiosensitive severe combined immunodeficiency. In ARTEMIS, specific mutations in two conserved regions affect interaction with DNA-PKcs. In combination they impair V(D)J recombination activity, independent of ARTEMIS autoinhibitory self-interaction between the ARTEMIS C-terminus and the N-terminal nuclease domain. We describe small fragments from both proteins, capable of interaction with the corresponding full-length partner proteins: In DNA-PKcs 42 amino acids out of FAT region 2 (PKcs3041-3082) can mediate interaction with ARTEMIS. In the nuclease we have defined 26 amino acids (ARM378-403) as minimal DNA-PKcs interacting fragment. The exact mapping of the ARTEMIS:DNA-PKcs interaction may pave the way for the design of specific inhibitors targeting the repair of DNA double strand breaks.

The coexistence of DNA replication and transcription during S-phase requires their tight coordination to prevent harmful conflicts. While extensive research revealed important mechanisms for minimizing these conflicts and their consequences, little is known regarding how the replication and transcription machinery are coordinated in real-time. Here, we developed a live-cell imaging approach for the real-time monitoring of replisome progression and transcription dynamics during a transcription-replication encounter. We found a wave of partial transcriptional repression ahead of the moving replication fork, which may contribute to efficient fork progression through the transcribed gene. Real-time detection of conflicts revealed their negative impact on both processes, leading to fork stalling or slowdown as well as lower transcription levels during gene replication, with different trade-offs observed in defined subpopulations of cells. Our real-time measurements of transcription-replication encounters demonstrate how these processes can proceed simultaneously while maintaining genomic stability, and how conflicts can arise when coordination is impaired.

Erratum in Nucleic Acids Res. 2022 Feb 25;:

Mouse models have been engineered to reveal the biological mechanisms of human diseases based on an assumption. The assumption is that orthologous genes underlie conserved phenotypes across species. However, genetically modified mouse orthologs of human genes do not often recapitulate human disease phenotypes which might be due to the molecular evolution of phenotypic differences across species from the time of the last common ancestor. Here, we systematically investigated the evolutionary divergence of regulatory relationships between transcription factors (TFs) and target genes in functional modules, and found that the rewiring of gene regulatory networks (GRNs) contributes to the phenotypic discrepancies that occur between humans and mice. We confirmed that the rewired regulatory networks of orthologous genes contain a higher proportion of species-specific regulatory elements. Additionally, we verified that the divergence of target gene expression levels, which was triggered by network rewiring, could lead to phenotypic differences. Taken together, a careful consideration of evolutionary divergence in regulatory networks could be a novel strategy to understand the failure or success of mouse models to mimic human diseases. To help interpret mouse phenotypes in human disease studies, we provide quantitative comparisons of gene expression profiles on our website (http://sbi.postech.ac.kr/w/RN).

The second cell fate decision in the early stage of mammalian embryonic development is pivotal; however, the underlying molecular mechanism is largely unexplored. Here, we report that Prmt1 acts as an important regulator in primitive endoderm (PrE) formation. First, Prmt1 depletion promotes PrE gene expression in mouse embryonic stem cells (ESCs). Single-cell RNA sequencing and flow cytometry assays demonstrated that Prmt1 depletion in mESCs contributes to an emerging cluster, where PrE genes are upregulated significantly. Furthermore, the efficiency of extraembryonic endoderm stem cell induction increased in Prmt1-depleted ESCs. Second, the pluripotency factor Klf4 methylated at Arg396 by Prmt1 is required for recruitment of the repressive mSin3a/HDAC complex to silence PrE genes. Most importantly, an embryonic chimeric assay showed that Prmt1 inhibition and mutated Klf4 at Arg 396 induce the integration of mouse ESCs into the PrE lineage. Therefore, we reveal a regulatory mechanism for cell fate decisions centered on Prmt1-mediated Klf4 methylation.

Erratum in Nucleic Acids Res. 2022 Nov 11;50(20):11993.

CTR9 is the scaffold subunit in polymerase-associated factor complex (PAFc), a multifunctional complex employed in multiple steps of RNA Polymerase II (RNAPII)-mediated transcription. CTR9/PAFc is well known as an evolutionarily conserved elongation factor that regulates gene activation via coupling with histone modifications enzymes. However, little is known about its function to restrain repressive histone markers. Using inducible and stable CTR9 knockdown breast cancer cell lines, we discovered that the H3K27me3 levels are strictly controlled by CTR9. Quantitative profiling of histone modifications revealed a striking increase of H3K27me3 levels upon loss of CTR9. Moreover, loss of CTR9 leads to genome-wide expansion of H3K27me3, as well as increased recruitment of PRC2 on chromatin, which can be reversed by CTR9 restoration. Further, CTR9 depletion triggers a PRC2 subtype switch from the less active PRC2.2, to the more active PRC2.1 with higher methyltransferase activity. As a consequence, CTR9 depletion generates vulnerability that renders breast cancer cells hypersensitive to PRC2 inhibitors. Our findings that CTR9 demarcates PRC2-mediated H3K27me3 levels and genomic distribution provide a unique mechanism that explains the transition from transcriptionally active chromatin states to repressive chromatin states and sheds light on the biological functions of CTR9 in development and cancer.

Subdomain 5BSL3.2 of hepatitis C virus RNA lies at the core of a network of distal RNA-RNA contacts that connect the 5' and 3' regions of the viral genome and regulate the translation and replication stages of the viral cycle. Using small-angle X-ray scattering and NMR spectroscopy experiments, we have determined at low resolution the structural models of this subdomain and its distal complex with domain 3'X, located at the 3'-terminus of the viral RNA chain. 5BSL3.2 adopts a characteristic 'L' shape in solution, whereas the 5BSL3.2-3'X distal complex forms a highly unusual 'Y'-shaped kissing junction that blocks the dimer linkage sequence of domain 3'X and promotes translation. The structure of this complex may impede an effective association of the viral polymerase with 5BSL3.2 and 3'X to start negative-strand RNA synthesis, contributing to explain the likely mechanism used by these sequences to regulate viral replication and translation. In addition, sequence and shape features of 5BSL3.2 are present in functional RNA motifs of flaviviruses, suggesting conserved regulatory processes within the Flaviviridae family.

Elongation factor TFIIS (transcription factor IIS) is structurally and biochemically probably the best characterized elongation cofactor of RNA polymerase II. However, little is known about TFIIS regulation or its roles during stress responses. Here, we show that, although TFIIS seems unnecessary under optimal conditions in Arabidopsis, its absence renders plants supersensitive to heat; tfIIs mutants die even when exposed to sublethal high temperature. TFIIS activity is required for thermal adaptation throughout the whole life cycle of plants, ensuring both survival and reproductive success. By employing a transcriptome analysis, we unravel that the absence of TFIIS makes transcriptional reprogramming sluggish, and affects expression and alternative splicing pattern of hundreds of heat-regulated transcripts. Transcriptome changes indirectly cause proteotoxic stress and deterioration of cellular pathways, including photosynthesis, which finally leads to lethality. Contrary to expectations of being constantly present to support transcription, we show that TFIIS is dynamically regulated. TFIIS accumulation during heat occurs in evolutionary distant species, including the unicellular alga Chlamydomonas reinhardtii, dicot Brassica napus and monocot Hordeum vulgare, suggesting that the vital role of TFIIS in stress adaptation of plants is conserved.

MicroRNAs (miRNAs or miRs) are single-stranded, ∼22-nucleotide noncoding RNAs that regulate many cellular processes. While numerous miRNA quantification technologies are available, a recent analysis of 12 commercial platforms revealed high variations in reproducibility, sensitivity, accuracy, specificity and concordance within and/or between platforms. Here, we developed a universal hairpin primer (UHP) system that negates the use of miRNA-specific hairpin primers (MsHPs) for quantitative reverse transcription PCR (RT-qPCR)-based miRNA quantification. Specifically, we analyzed four UHPs that share the same hairpin structure but are anchored with two, three, four and six degenerate nucleotides at 3'-ends (namely UHP2, UHP3, UHP4 and UHP6), and found that the four UHPs yielded robust RT products and quantified miRNAs with high efficiency. UHP-based RT-qPCR miRNA quantification was not affected by long transcripts. By analyzing 14 miRNAs, we demonstrated that UHP4 closely mimicked MsHPs in miRNA quantification. Fine-tuning experiments identified an optimized UHP (OUHP) mix with a molar composition of UHP2:UHP4:UHP6 = 8:1:1, which closely recapitulated MsHPs in miRNA quantification. Using synthetic LET7 isomiRs, we demonstrated that the OUHP-based qPCR system exhibited high specificity and sensitivity. Collectively, our results demonstrate that the OUHP system can serve as a reliable and cost-effective surrogate of MsHPs for RT-qPCR-based miRNA quantification for basic research and precision medicine.

RNase H is involved in fundamental cellular processes and is responsible for removing the short stretch of RNA from Okazaki fragments and the long stretch of RNA from R-loops. Defects in RNase H lead to embryo lethality in mice and Aicardi-Goutieres syndrome in humans, suggesting the importance of RNase H. To date, RNase H is known to be a non-sequence-specific endonuclease, but it is not known whether it performs other functions on the structural variants of RNA:DNA hybrids. Here, we used Escherichia coli RNase H as a model, and examined its catalytic mechanism and its substrate recognition modes, using single-molecule FRET. We discovered that RNase H acts as a processive exoribonuclease on the 3' DNA overhang side but as a distributive non-sequence-specific endonuclease on the 5' DNA overhang side of RNA:DNA hybrids or on blunt-ended hybrids. The high affinity of previously unidentified double-stranded (ds) and single-stranded (ss) DNA junctions flanking RNA:DNA hybrids may help RNase H find the hybrid substrates in long genomic DNA. Our study provides new insights into the multifunctionality of RNase H, elucidating unprecedented roles of junctions and ssDNA overhang on RNA:DNA hybrids.

Bacteria and archaea use the CRISPR-Cas system to fend off invasions of bacteriophages and foreign plasmids. In response, bacteriophages encode anti-CRISPR (Acr) proteins that potently inhibit host Cas proteins to suppress CRISPR-mediated immunity. AcrIE4-F7, which was isolated from Pseudomonas citronellolis, is a fused form of AcrIE4 and AcrIF7 that inhibits both type I-E and type I-F CRISPR-Cas systems. Here, we determined the structure of AcrIE4-F7 and identified its Cas target proteins. The N-terminal AcrIE4 domain adopts a novel α-helical fold that targets the PAM interaction site of the type I-E Cas8e subunit. The C-terminal AcrIF7 domain exhibits an αβ fold like native AcrIF7, which disables target DNA recognition by the PAM interaction site in the type I-F Cas8f subunit. The two Acr domains are connected by a flexible linker that allows prompt docking onto their cognate Cas8 targets. Conserved negative charges in each Acr domain are required for interaction with their Cas8 targets. Our results illustrate a common mechanism by which AcrIE4-F7 inhibits divergent CRISPR-Cas types.

The shelterin protein complex is required for telomere protection in various eukaryotic organisms. In mammals, the shelterin subunit TRF2 is specialized in preventing ATM activation at telomeres and chromosome end fusion in somatic cells. Here, we demonstrate that the zebrafish ortholog of TRF2 (encoded by the terfa gene) is protecting against unwanted ATM activation genome-wide. The terfa-compromised fish develop a prominent and specific embryonic neurodevelopmental failure. The heterozygous fish survive to adulthood but exhibit a premature aging phenotype. The recovery from embryonic neurodevelopmental failure requires both ATM inhibition and transcriptional complementation of neural genes. Furthermore, restoring the expression of TRF2 in glial cells rescues the embryonic neurodevelopment phenotype. These results indicate that the shelterin subunit TRF2 evolved in zebrafish as a general factor of genome maintenance and transcriptional regulation that is required for proper neurodevelopment and normal aging. These findings uncover how TRF2 links development to aging by separate functions in gene expression regulation and genome stability control.

Staphylococcus aureus is a notorious and globally distributed pathogenic bacterium. New strategies to develop novel antibiotics based on intrinsic bacterial toxin-antitoxin (TA) systems have been recently reported. Because TA systems are present only in bacteria and not in humans, these distinctive systems are attractive targets for developing antibiotics with new modes of action. S. aureus PemIK is a type II TA system, comprising the toxin protein PemK and the labile antitoxin protein PemI. Here, we determined the crystal structures of both PemK and the PemIK complex, in which PemK is neutralized by PemI. Our biochemical approaches, including fluorescence quenching and polarization assays, identified Glu20, Arg25, Thr48, Thr49, and Arg84 of PemK as being important for RNase function. Our study indicates that the active site and RNA-binding residues of PemK are covered by PemI, leading to unique conformational changes in PemK accompanied by repositioning of the loop between β1 and β2. These changes can interfere with RNA binding by PemK. Overall, PemK adopts particular open and closed forms for precise neutralization by PemI. This structural and functional information on PemIK will contribute to the discovery and development of novel antibiotics in the form of peptides or small molecules inhibiting direct binding between PemI and PemK.

The first member of the pleuromutilin (PLM) class suitable for systemic antibacterial chemotherapy in humans recently entered clinical use, underscoring the need to better understand mechanisms of PLM resistance in disease-causing bacterial genera. Of the proteins reported to mediate PLM resistance in staphylococci, the least-well studied to date is Sal(A), a putative ABC-F NTPase that-by analogy to other proteins of this type-may act to protect the ribosome from PLMs. Here, we establish the importance of Sal proteins as a common source of PLM resistance across multiple species of staphylococci. Sal(A) is revealed as but one member of a larger group of Sal-type ABC-F proteins that vary considerably in their ability to mediate resistance to PLMs and other antibiotics. We find that specific sal genes are intrinsic to particular staphylococcal species, and show that this gene family is likely ancestral to the genus Staphylococcus. Finally, we solve the cryo-EM structure of a representative Sal-type protein (Sal(B)) in complex with the staphylococcal 70S ribosome, revealing that Sal-type proteins bind into the E site to mediate target protection, likely by displacing PLMs and other antibiotics via an allosteric mechanism.

Human genetic studies identified a strong association between loss of function mutations in RBFOX2 and hypoplastic left heart syndrome (HLHS). There are currently no Rbfox2 mouse models that recapitulate HLHS. Therefore, it is still unknown how RBFOX2 as an RNA binding protein contributes to heart development. To address this, we conditionally deleted Rbfox2 in embryonic mouse hearts and found profound defects in cardiac chamber and yolk sac vasculature formation. Importantly, our Rbfox2 conditional knockout mouse model recapitulated several molecular and phenotypic features of HLHS. To determine the molecular drivers of these cardiac defects, we performed RNA-sequencing in Rbfox2 mutant hearts and identified dysregulated alternative splicing (AS) networks that affect cell adhesion to extracellular matrix (ECM) mediated by Rho GTPases. We identified two Rho GTPase cycling genes as targets of RBFOX2. Modulating AS of these two genes using antisense oligos led to cell cycle and cell-ECM adhesion defects. Consistently, Rbfox2 mutant hearts displayed cell cycle defects and inability to undergo endocardial-mesenchymal transition, processes dependent on cell-ECM adhesion and that are seen in HLHS. Overall, our work not only revealed that loss of Rbfox2 leads to heart development defects resembling HLHS, but also identified RBFOX2-regulated AS networks that influence cell-ECM communication vital for heart development.

The DNA damage response (DDR) preserves the genetic integrity of the cell by sensing and repairing damages after a genotoxic stress. Translesion Synthesis (TLS), an error-prone DNA damage tolerance pathway, is controlled by PCNA ubiquitination. In this work, we raise the question whether TLS is controlled locally or globally. Using a recently developed method that allows to follow the bypass of a single lesion inserted into the yeast genome, we show that (i) TLS is controlled locally at each individual lesion by PCNA ubiquitination, (ii) a single lesion is enough to induce PCNA ubiquitination and (iii) PCNA ubiquitination is imperative for TLS to occur. More importantly, we show that the activation of the DDR that follows a genotoxic stress does not increase TLS at individual lesions. We conclude that unlike the SOS response in bacteria, the eukaryotic DDR does not promote TLS and mutagenesis.

Only a select few L1 loci in the human genome are expressed in any given cell line or organ, likely to minimize damage done to the genome. The epigenetic features and requirements of expressed L1 loci are currently unknown. Using human cells and comprehensive epigenetic analysis of individual expressed and unexpressed L1 loci, we determined that endogenous L1 transcription depends on a combination of epigenetic factors, including open chromatin, activating histone modifications, and hypomethylation at the L1 promoter. We demonstrate that the L1 promoter seems to require interaction with enhancer elements for optimal function. We utilize epigenetic context to predict the expression status of L1Hs loci that are poorly mappable with RNA-Seq. Our analysis identified a population of 'transitional' L1 loci that likely have greater potential to be activated during the epigenetic dysregulation seen in tumors and during aging because they are the most responsive to targeted CRISPR-mediated delivery of trans-activating domains. We demonstrate that an engineered increase in endogenous L1 mRNA expression increases Alu mobilization. Overall, our findings present the first global and comprehensive analysis of epigenetic status of individual L1 loci based on their expression status and demonstrate the importance of epigenetic context for L1 expression heterogeneity.

Unlike many other aminoacyl-tRNA synthetases, alanyl-tRNA synthetase (AlaRS) retains a conserved prototype structure throughout biology. While Caenorhabditis elegans cytoplasmic AlaRS (CeAlaRSc) retains the prototype structure, its mitochondrial counterpart (CeAlaRSm) contains only a residual C-terminal domain (C-Ala). We demonstrated herein that the C-Ala domain from CeAlaRSc robustly binds both tRNA and DNA. It bound different tRNAs but preferred tRNAAla. Deletion of this domain from CeAlaRSc sharply reduced its aminoacylation activity, while fusion of this domain to CeAlaRSm selectively and distinctly enhanced its aminoacylation activity toward the elbow-containing (or L-shaped) tRNAAla. Phylogenetic analysis showed that CeAlaRSm once possessed the C-Ala domain but later lost most of it during evolution, perhaps in response to the deletion of the T-arm (part of the elbow) from its cognate tRNA. This study underscores the evolutionary gain of C-Ala for docking AlaRS to the L-shaped tRNAAla.

Conflicts between transcription and replication machinery are a potent source of replication stress and genome instability; however, no technique currently exists to identify endogenous genomic locations prone to transcription-replication interactions. Here, we report a novel method to identify genomic loci prone to transcription-replication interactions termed transcription-replication immunoprecipitation on nascent DNA sequencing, TRIPn-Seq. TRIPn-Seq employs the sequential immunoprecipitation of RNA polymerase 2 phosphorylated at serine 5 (RNAP2s5) followed by enrichment of nascent DNA previously labeled with bromodeoxyuridine. Using TRIPn-Seq, we mapped 1009 unique transcription-replication interactions (TRIs) in mouse primary B cells characterized by a bimodal pattern of RNAP2s5, bidirectional transcription, an enrichment of RNA:DNA hybrids, and a high probability of forming G-quadruplexes. TRIs are highly enriched at transcription start sites and map to early replicating regions. TRIs exhibit enhanced Replication Protein A association and TRI-associated genes exhibit higher replication fork termination than control transcription start sites, two marks of replication stress. TRIs colocalize with double-strand DNA breaks, are enriched for deletions, and accumulate mutations in tumors. We propose that replication stress at TRIs induces mutations potentially contributing to age-related disease, as well as tumor formation and development.

Approximately half of purified mammalian RNA polymerase II (Pol II) is associated with a tightly interacting sub-stoichiometric subunit, Gdown1. Previous studies have established that Gdown1 inhibits transcription initiation through competitive interactions with general transcription factors and blocks the Pol II termination activity of transcription termination factor 2 (TTF2). However, the biological functions of Gdown1 remain poorly understood. Here, we utilized genetic, microscopic, and multi-omics approaches to functionally characterize Gdown1 in three human cell lines. Acute depletion of Gdown1 caused minimal direct effects on transcription. We show that Gdown1 resides predominantly in the cytoplasm of interphase cells, shuttles between the cytoplasm and nucleus, and is regulated by nuclear export. Gdown1 enters the nucleus at the onset of mitosis. Consistently, genetic ablation of Gdown1 is associated with partial de-repression of mitotic transcription, and Gdown1 KO cells present with evidence of aberrant mitoses coupled to p53 pathway activation. Evidence is presented demonstrating that Gdown1 modulates the combined functions of purified productive elongation factors PAF1C, RTF1, SPT6, DSIF and P-TEFb in vitro. Collectively, our findings support a model wherein the Pol II-regulatory function of Gdown1 occurs during mitosis and is required for genome integrity.

Ribosomes have long been thought of as homogeneous macromolecular machines, but recent evidence suggests they are heterogeneous and could be specialised to regulate translation. Here, we have characterised ribosomal protein heterogeneity across 4 tissues of Drosophila melanogaster. We find that testes and ovaries contain the most heterogeneous ribosome populations, which occurs through a combination of paralog-enrichment and paralog-switching. We have solved structures of ribosomes purified from in vivo tissues by cryo-EM, revealing differences in precise ribosomal arrangement for testis and ovary 80S ribosomes. Differences in the amino acid composition of paralog pairs and their localisation on the ribosome exterior indicate paralog-switching could alter the ribosome surface, enabling different proteins to regulate translation. One testis-specific paralog-switching pair is also found in humans, suggesting this is a conserved site of ribosome heterogeneity. Overall, this work allows us to propose that mRNA translation might be regulated in the gonads through ribosome heterogeneity, providing a potential means of ribosome specialisation.

Erratum in Nucleic Acids Res. 2022 Mar 21;50(5):2995.

Tight control of gene expression networks required for adipose tissue formation and plasticity is essential for adaptation to energy needs and environmental cues. However, the mechanisms that orchestrate the global and dramatic transcriptional changes leading to adipocyte differentiation remain to be fully unraveled. We investigated the regulation of nascent transcription by the sumoylation pathway during adipocyte differentiation using SLAMseq and ChIPseq. We discovered that the sumoylation pathway has a dual function in differentiation; it supports the initial downregulation of pre-adipocyte-specific genes, while it promotes the establishment of the mature adipocyte transcriptional program. By characterizing endogenous sumoylome dynamics in differentiating adipocytes by mass spectrometry, we found that sumoylation of specific transcription factors like PPARγ/RXR and their co-factors are associated with the transcription of adipogenic genes. Finally, using RXR as a model, we found that sumoylation may regulate adipogenic transcription by supporting the chromatin occurrence of transcription factors. Our data demonstrate that the sumoylation pathway supports the rewiring of transcriptional networks required for formation of functional adipocytes. This study also provides the scientists in the field of cellular differentiation and development with an in-depth resource of the dynamics of the SUMO-chromatin landscape, SUMO-regulated transcription and endogenous sumoylation sites during adipocyte differentiation.

Translation Complex Profile Sequencing (TCP-seq), a protocol that was developed and implemented on Saccharomyces cerevisiae, provides the footprints of the small subunit (SSU) of the ribosome (with additional factors) across the entire transcriptome of the analyzed organism. In this study, based on the TCP-seq data, we developed for the first-time a predictive model of the SSU density and analyzed the effect of transcript features on the dynamics of the SSU scan in the 5'UTR. Among others, our model is based on novel tools for detecting complex statistical relations tailored to TCP-seq. We quantitatively estimated the effect of several important features, including the context of the upstream AUG, the upstream ORF length and the mRNA folding strength. Specifically, we suggest that around 50% of the variance related to the read counts (RC) distribution near a start codon can be attributed to the AUG context score. We provide the first large scale direct quantitative evidence that shows that indeed AUG context affects the small sub-unit movement. In addition, we suggest that strong folding may cause the detachment of the SSU from the mRNA. We also identified a number of novel sequence motifs that can affect the SSU scan; some of these motifs affect transcription factors and RNA binding proteins. The results presented in this study provide a better understanding of the biophysical aspects related to the SSU scan along the 5'UTR and of translation initiation in S. cerevisiae, a fundamental step toward a comprehensive modeling of initiation.

In the last years, many studies were able to identify associations between common genetic variants and complex diseases. However, the mechanistic biological links explaining these associations are still mostly unknown. Common variants are usually associated with a relatively small effect size, suggesting that interactions among multiple variants might be a major genetic component of complex diseases. Hence, elucidating the presence of functional relations among variants may be fundamental to identify putative variants' interactions. To this aim, we developed Polympact, a web-based resource that allows to explore functional relations among human common variants by exploiting variants' functional element landscape, their impact on transcription factor binding motifs, and their effect on transcript levels of protein-coding genes. Polympact characterizes over 18 million common variants and allows to explore putative relations by combining clustering analysis and innovative similarity and interaction network models. The properties of the network models were studied and the utility of Polympact was demonstrated by analysing the rich sets of Breast Cancer and Alzheimer's GWAS variants. We identified relations among multiple variants, suggesting putative interactions. Polympact is freely available at bcglab.cibio.unitn.it/polympact.

Chromosome replication depends on efficient removal of nucleosomes by accessory factors to ensure rapid access to genomic information. Here, we show this process requires recruitment of the nucleosome reorganization activity of the histone chaperone FACT. Using single-molecule FRET, we demonstrate that reorganization of nucleosomal DNA by FACT requires coordinated engagement by the middle and C-terminal domains of Spt16 and Pob3 but does not require the N-terminus of Spt16. Using structure-guided pulldowns, we demonstrate instead that the N-terminal region is critical for recruitment by the fork protection complex subunit Tof1. Using in vitro chromatin replication assays, we confirm the importance of these interactions for robust replication. Our findings support a mechanism in which nucleosomes are removed through the coordinated engagement of multiple FACT domains positioned at the replication fork by the fork protection complex.

CRISPR-Cas systems provide prokaryotic organisms with an adaptive defense mechanism that acquires immunological memories of infections. This is accomplished by integration of short fragments from the genome of invaders such as phages and plasmids, called 'spacers', into the CRISPR locus of the host. Depending on their genetic composition, CRISPR-Cas systems can be classified into six types, I-VI, however spacer acquisition has been extensively studied only in type I and II systems. Here, we used an inducible spacer acquisition assay to study this process in the type III-A CRISPR-Cas system of Staphylococcus epidermidis, in the absence of phage selection. Similarly to type I and II spacer acquisition, this type III system uses Cas1 and Cas2 to preferentially integrate spacers from the chromosomal terminus and free dsDNA ends produced after DNA breaks, in a manner that is enhanced by the AddAB DNA repair complex. Surprisingly, a different mode of spacer acquisition from rRNA and tRNA loci, which spans only the transcribed sequences of these genes and is not enhanced by AddAB, was also detected. Therefore, our findings reveal both common mechanistic principles that may be conserved in all CRISPR-Cas systems, as well as unique and intriguing features of type III spacer acquisition.

The SARS-CoV-2 coronavirus is the causal agent of the current global pandemic. SARS-CoV-2 belongs to an order, Nidovirales, with very large RNA genomes. It is proposed that the fidelity of coronavirus (CoV) genome replication is aided by an RNA nuclease complex, comprising the non-structural proteins 14 and 10 (nsp14-nsp10), an attractive target for antiviral inhibition. Our results validate reports that the SARS-CoV-2 nsp14-nsp10 complex has RNase activity. Detailed functional characterization reveals nsp14-nsp10 is a versatile nuclease capable of digesting a wide variety of RNA structures, including those with a blocked 3'-terminus. Consistent with a role in maintaining viral genome integrity during replication, we find that nsp14-nsp10 activity is enhanced by the viral RNA-dependent RNA polymerase complex (RdRp) consisting of nsp12-nsp7-nsp8 (nsp12-7-8) and demonstrate that this stimulation is mediated by nsp8. We propose that the role of nsp14-nsp10 in maintaining replication fidelity goes beyond classical proofreading by purging the nascent replicating RNA strand of a range of potentially replication-terminating aberrations. Using our developed assays, we identify drug and drug-like molecules that inhibit nsp14-nsp10, including the known SARS-CoV-2 major protease (Mpro) inhibitor ebselen and the HIV integrase inhibitor raltegravir, revealing the potential for multifunctional inhibitors in COVID-19 treatment.

Pervasive transcription is a universal phenomenon leading to the production of a plethora of non-coding RNAs. If left uncontrolled, pervasive transcription can be harmful for genome expression and stability. However, non-coding transcription can also play important regulatory roles, for instance by promoting the repression of specific genes by a mechanism of transcriptional interference. The efficiency of transcription termination can strongly influence the regulatory capacity of non-coding transcription events, yet very little is known about the mechanisms modulating the termination of non-coding transcription in response to environmental cues. Here, we address this question by investigating the mechanisms that regulate the activity of the main actor in termination of non-coding transcription in budding yeast, the helicase Sen1. We identify a phosphorylation at a conserved threonine of the catalytic domain of Sen1 and we provide evidence that phosphorylation at this site reduces the efficiency of Sen1-mediated termination. Interestingly, we find that this phosphorylation impairs termination at an unannotated non-coding gene, thus repressing the expression of a downstream gene encoding the master regulator of Zn homeostasis, Zap1. Consequently, many additional genes exhibit an expression pattern mimicking conditions of Zn excess, where ZAP1 is naturally repressed. Our findings provide a novel paradigm of gene regulatory mechanism relying on the direct modulation of non-coding transcription termination.

Type III CRISPR-Cas systems have a unique mode of interference, involving crRNA-guided recognition of nascent RNA and leading to DNA and RNA degradation. How type III systems acquire new CRISPR spacers is currently not well understood. Here, we characterize CRISPR spacer uptake by a type III-A system within its native host, Streptococcus thermophilus. Adaptation by the type II-A system in the same host provided a basis for comparison. Cas1 and Cas2 proteins were critical for type III adaptation but deletion of genes responsible for crRNA biogenesis or interference did not detectably change spacer uptake patterns, except those related to host counter-selection. Unlike the type II-A system, type III spacers are acquired in a PAM- and orientation-independent manner. Interestingly, certain regions of plasmids and the host genome were particularly well-sampled during type III-A, but not type II-A, spacer uptake. These regions included the single-stranded origins of rolling-circle replicating plasmids, rRNA and tRNA encoding gene clusters, promoter regions of expressed genes and 5' UTR regions involved in transcription attenuation. These features share the potential to form DNA secondary structures, suggesting a preferred substrate for type III adaptation. Lastly, the type III-A system adapted to and protected host cells from lytic phage infection.

Recombineering assisted multiplex genome editing generally uses single-stranded oligonucleotides for site directed mutational changes. It has proven highly efficient for functional screens and to optimize microbial cell factories. However, this approach is limited to relatively small mutational changes. Here, we addressed the challenges involved in the use of double-stranded DNA substrates for multiplex genome engineering. Recombineering is mediated by phage single-strand annealing proteins annealing ssDNAs into the replication fork. We apply this insight to facilitate the generation of ssDNA from the dsDNA substrate and to alter the speed of replication by elevating the available deoxynucleoside triphosphate (dNTP) levels. Intracellular dNTP concentration was elevated by ribonucleotide reductase overexpression or dNTP addition to establish double-stranded DNA Recombineering-assisted Multiplex Genome Engineering (dReaMGE), which enables rapid and flexible insertional and deletional mutagenesis at multiple sites on kilobase scales in diverse bacteria without the generation of double-strand breaks or disturbance of the mismatch repair system. dReaMGE can achieve combinatorial genome engineering works, for example, alterations to multiple biosynthetic pathways, multiple promoter or gene insertions, variations of transcriptional regulator combinations, within a few days. dReaMGE adds to the repertoire of bacterial genome engineering to facilitate discovery, functional genomics, strain optimization and directed evolution of microbial cell factories.

For many RNA molecules, the secondary structure is essential for the correct function of the RNA. Predicting RNA secondary structure from nucleotide sequences is a long-standing problem in genomics, but the prediction performance has reached a plateau over time. Traditional RNA secondary structure prediction algorithms are primarily based on thermodynamic models through free energy minimization, which imposes strong prior assumptions and is slow to run. Here, we propose a deep learning-based method, called UFold, for RNA secondary structure prediction, trained directly on annotated data and base-pairing rules. UFold proposes a novel image-like representation of RNA sequences, which can be efficiently processed by Fully Convolutional Networks (FCNs). We benchmark the performance of UFold on both within- and cross-family RNA datasets. It significantly outperforms previous methods on within-family datasets, while achieving a similar performance as the traditional methods when trained and tested on distinct RNA families. UFold is also able to predict pseudoknots accurately. Its prediction is fast with an inference time of about 160 ms per sequence up to 1500 bp in length. An online web server running UFold is available at https://ufold.ics.uci.edu. Code is available at https://github.com/uci-cbcl/UFold.

Eukaryotic cells recognize intracellular pathogens through pattern recognition receptors, including sensors of aberrant nucleic acid structures. Sensors of double-stranded RNA (dsRNA) are known to detect replication intermediates of RNA viruses. It has long been suggested that annealing of mRNA from symmetrical transcription of both top and bottom strands of DNA virus genomes can produce dsRNA during infection. Supporting this hypothesis, nearly all DNA viruses encode inhibitors of dsRNA-recognition pathways. However, direct evidence that DNA viruses produce dsRNA is lacking. Contrary to dogma, we show that the nuclear-replicating DNA virus adenovirus (AdV) does not produce detectable levels of dsRNA during infection. In contrast, abundant dsRNA is detected within the nucleus of cells infected with AdV mutants defective for viral RNA processing. In the presence of nuclear dsRNA, the cytoplasmic dsRNA sensor PKR is relocalized and activated within the nucleus. Accumulation of viral dsRNA occurs in the late phase of infection, when unspliced viral transcripts form intron/exon base pairs between top and bottom strand transcripts. We propose that DNA viruses actively limit dsRNA formation by promoting efficient splicing and mRNA processing, thus avoiding detection and restriction by host innate immune sensors of pathogenic nucleic acids.

A critical challenge for the successful development of RNA interference-based therapeutics therapeutics has been the enhancement of their in vivo metabolic stability. In therapeutically relevant, fully chemically modified small interfering RNAs (siRNAs), modification of the two terminal phosphodiester linkages in each strand of the siRNA duplex with phosphorothioate (PS) is generally sufficient to protect against exonuclease degradation in vivo. Since PS linkages are chiral, we systematically studied the properties of siRNAs containing single chiral PS linkages at each strand terminus. We report an efficient and simple method to introduce chiral PS linkages and demonstrate that Rp diastereomers at the 5' end and Sp diastereomers at the 3' end of the antisense siRNA strand improved pharmacokinetic and pharmacodynamic properties in a mouse model. In silico modeling studies provide mechanistic insights into how the Rp isomer at the 5' end and Sp isomer at the 3' end of the antisense siRNA enhance Argonaute 2 (Ago2) loading and metabolic stability of siRNAs in a concerted manner.

G-quadruplex (G4) structures that can form at guanine-rich genomic sites, including telomeres and gene promoters, are actively involved in genome maintenance, replication, and transcription, through finely tuned interactions with protein networks. In the present study, we identified the intermediate filament protein Vimentin as a binder with nanomolar affinity for those G-rich sequences that give rise to at least two adjacent G4 units, named G4 repeats. This interaction is supported by the N-terminal domains of soluble Vimentin tetramers. The selectivity of Vimentin for G4 repeats versus individual G4s provides an unprecedented result. Based on GO enrichment analysis performed on genes having putative G4 repeats within their core promoters, we suggest that Vimentin recruitment at these sites may contribute to the regulation of gene expression during cell development and migration, possibly by reshaping the local higher-order genome topology, as already reported for lamin B.

Giardiasis is a disease caused by the protist Giardia lamblia. As no human vaccines have been approved so far against it, and resistance to current drugs is spreading, new strategies for combating giardiasis need to be developed. The G. lamblia ribosome may provide a promising therapeutic target due to its distinct sequence differences from ribosomes of most eukaryotes and prokaryotes. Here, we report the cryo-electron microscopy structure of the G. lamblia (WB strain) ribosome determined at 2.75 Å resolution. The ribosomal RNA is the shortest known among eukaryotes, and lacks nearly all the eukaryote-specific ribosomal RNA expansion segments. In contrast, the ribosomal proteins are typically eukaryotic with some species-specific insertions/extensions. Most typical inter-subunit bridges are maintained except for one missing contact site. Unique structural features are located mainly at the ribosome's periphery. These may be exploited as target sites for the design of new compounds that inhibit selectively the parasite's ribosomal activity.

The rational design and realisation of simple-to-use genetic control elements that are modular, orthogonal and robust is essential to the construction of predictable and reliable biological systems of increasing complexity. To this effect, we introduce modular Artificial RNA interference (mARi), a rational, modular and extensible design framework that enables robust, portable and multiplexed post-transcriptional regulation of gene expression in Escherichia coli. The regulatory function of mARi was characterised in a range of relevant genetic contexts, demonstrating its independence from other genetic control elements and the gene of interest, and providing new insight into the design rules of RNA based regulation in E. coli, while a range of cellular contexts also demonstrated it to be independent of growth-phase and strain type. Importantly, the extensibility and orthogonality of mARi enables the simultaneous post-transcriptional regulation of multi-gene systems as both single-gene cassettes and poly-cistronic operons. To facilitate adoption, mARi was designed to be directly integrated into the modular BASIC DNA assembly framework. We anticipate that mARi-based genetic control within an extensible DNA assembly framework will facilitate metabolic engineering, layered genetic control, and advanced genetic circuit applications.

DNA methylation modulates telomere function. In Arabidopsis thaliana, telomeric regions have a bimodal chromatin organization with unmethylated telomeres and methylated subtelomeres. To gain insight into this organization we have generated TAIR10-Tel, a modified version of the Arabidopsis reference genome with additional sequences at most chromosome ends. TAIR10-Tel has allowed us to analyse DNA methylation at nucleotide resolution level in telomeric regions. We have analysed the wild-type strain and mutants that encode inactive versions of all currently known relevant methyltransferases involved in cytosine methylation. These analyses have revealed that subtelomeric DNA methylation extends 1 to 2 kbp from Interstitial Telomeric Sequences (ITSs) that abut or are very near to telomeres. However, DNA methylation drops at the telomeric side of the telomere-subtelomere boundaries and disappears at the inner part of telomeres. We present a comprehensive and integrative model for subtelomeric DNA methylation that should help to decipher the mechanisms that govern the epigenetic regulation of telomeres. This model involves a complex network of interactions between methyltransferases and subtelomeric DNA sequences.

Homologous recombination (HR) is critical for error-free repair of DNA double-strand breaks. Chromatin loading of RAD51, a key protein that mediates the recombination, is a crucial step in the execution of the HR repair. Here, we present evidence that SUMOylation of RAD51 is crucial for the RAD51 recruitment to chromatin and HR repair. We found that topoisomerase 1-binding arginine/serine-rich protein (TOPORS) induces the SUMOylation of RAD51 at lysine residues 57 and 70 in response to DNA damaging agents. The SUMOylation was facilitated by an ATM-induced phosphorylation of TOPORS at threonine 515 upon DNA damage. Knockdown of TOPORS or expression of SUMOylation-deficient RAD51 mutants caused reduction in supporting normal RAD51 functions during the HR repair, suggesting the physiological importance of the modification. We found that the SUMOylation-deficient RAD51 reduces the association with its crucial binding partner BRCA2, explaining its deficiency in supporting the HR repair. These findings altogether demonstrate a crucial role for TOPORS-mediated RAD51 SUMOylation in promoting HR repair and genomic maintenance.

A prominent aspect of most, if not all, central nervous systems (CNSs) is that anterior regions (brain) are larger than posterior ones (spinal cord). Studies in Drosophila and mouse have revealed that Polycomb Repressor Complex 2 (PRC2), a protein complex responsible for applying key repressive histone modifications, acts by several mechanisms to promote anterior CNS expansion. However, it is unclear what the full spectrum of PRC2 action is during embryonic CNS development and how PRC2 intersects with the epigenetic landscape. We removed PRC2 function from the developing mouse CNS, by mutating the key gene Eed, and generated spatio-temporal transcriptomic data. To decode the role of PRC2, we developed a method that incorporates standard statistical analyses with probabilistic deep learning to integrate the transcriptomic response to PRC2 inactivation with epigenetic data. This multi-variate analysis corroborates the central involvement of PRC2 in anterior CNS expansion, and also identifies several unanticipated cohorts of genes, such as proliferation and immune response genes. Furthermore, the analysis reveals specific profiles of regulation via PRC2 upon these gene cohorts. These findings uncover a differential logic for the role of PRC2 upon functionally distinct gene cohorts that drive CNS anterior expansion. To support the analysis of emerging multi-modal datasets, we provide a novel bioinformatics package that integrates transcriptomic and epigenetic datasets to identify regulatory underpinnings of heterogeneous biological processes.

Amplification of short interfering RNA (siRNAs) via RNA-dependent RNA polymerases (RdRPs) is of fundamental importance in RNA silencing. Plant microRNA (miRNA) action generally does not involve engagement of RdRPs, in part thanks to a poorly understood activity of the cytoplasmic exosome adaptor SKI2. Here, we show that inactivation of the exosome subunit RRP45B and SKI2 results in similar patterns of miRNA-induced siRNA production. Furthermore, loss of the nuclear exosome adaptor HEN2 leads to secondary siRNA production from miRNA targets largely distinct from those producing siRNAs in ski2. Importantly, mutation of the Release Factor paralogue PELOTA1 required for subunit dissociation of stalled ribosomes causes siRNA production from miRNA targets overlapping with, but distinct from, those affected in ski2 and rrp45b mutants. We also show that in exosome mutants, miRNA targets can be sorted into producers and non-producers of illicit secondary siRNAs based on trigger miRNA levels and miRNA:target affinity rather than on presence of 5'-cleavage fragments. We propose that stalled RNA-Induced Silencing Complex (RISC) and ribosomes, but not mRNA cleavage fragments released from RISC, trigger siRNA production, and that the exosome limits siRNA amplification by reducing RISC dwell time on miRNA target mRNAs while PELOTA1 does so by reducing ribosome stalling.

p53-binding protein 1 (53BP1) regulates the DNA double-strand break (DSB) repair pathway and maintains genomic integrity. Here we found that 53BP1 functions as a molecular scaffold for the nucleoside diphosphate kinase-mediated phosphorylation of ATP-citrate lyase (ACLY) which enhances the ACLY activity. This functional association is critical for promoting global histone acetylation and subsequent transcriptome-wide alterations in gene expression. Specifically, expression of a replication-dependent histone biogenesis factor, stem-loop binding protein (SLBP), is dependent upon 53BP1-ACLY-controlled acetylation at the SLBP promoter. This chain of regulation events carried out by 53BP1, ACLY, and SLBP is crucial for both quantitative and qualitative histone biogenesis as well as for the preservation of genomic integrity. Collectively, our findings reveal a previously unknown role for 53BP1 in coordinating replication-dependent histone biogenesis and highlight a DNA repair-independent function in the maintenance of genomic stability through a regulatory network that includes ACLY and SLBP.

In many cases, the unprecedented availability of data provided by high-throughput sequencing has shifted the bottleneck from a data availability issue to a data interpretation issue, thus delaying the promised breakthroughs in genetics and precision medicine, for what concerns Human genetics, and phenotype prediction to improve plant adaptation to climate change and resistance to bioagressors, for what concerns plant sciences. In this paper, we propose a novel Genome Interpretation paradigm, which aims at directly modeling the genotype-to-phenotype relationship, and we focus on A. thaliana since it is the best studied model organism in plant genetics. Our model, called Galiana, is the first end-to-end Neural Network (NN) approach following the genomes in/phenotypes out paradigm and it is trained to predict 288 real-valued Arabidopsis thaliana phenotypes from Whole Genome sequencing data. We show that 75 of these phenotypes are predicted with a Pearson correlation ≥0.4, and are mostly related to flowering traits. We show that our end-to-end NN approach achieves better performances and larger phenotype coverage than models predicting single phenotypes from the GWAS-derived known associated genes. Galiana is also fully interpretable, thanks to the Saliency Maps gradient-based approaches. We followed this interpretation approach to identify 36 novel genes that are likely to be associated with flowering traits, finding evidence for 6 of them in the existing literature.

The Cas9 nuclease from Staphylococcus aureus (SaCas9) holds great potential for use in gene therapy, and variants with increased fidelity have been engineered. However, we find that existing variants have not reached the greatest accuracy to discriminate base mismatches and exhibited much reduced activity when their mutations were grafted onto the KKH mutant of SaCas9 for editing an expanded set of DNA targets. We performed structure-guided combinatorial mutagenesis to re-engineer KKH-SaCas9 with enhanced accuracy. We uncover that introducing a Y239H mutation on KKH-SaCas9's REC domain substantially reduces off-target edits while retaining high on-target activity when added to a set of mutations on REC and RuvC domains that lessen its interactions with the target DNA strand. The Y239H mutation is modelled to have removed an interaction from the REC domain with the guide RNA backbone in the guide RNA-DNA heteroduplex structure. We further confirmed the greatly improved genome-wide editing accuracy and single-base mismatch discrimination of our engineered variants, named KKH-SaCas9-SAV1 and SAV2, in human cells. In addition to generating broadly useful KKH-SaCas9 variants with unprecedented accuracy, our findings demonstrate the feasibility for multi-domain combinatorial mutagenesis on SaCas9's DNA- and guide RNA- interacting residues to optimize its editing fidelity.

The life of RNA polymerase II (RNAPII) transcripts is shaped by the dynamic formation of mutually exclusive ribonucleoprotein complexes (RNPs) that direct transcript biogenesis and turnover. A key regulator of RNA metabolism in the nucleus is the scaffold protein ARS2 (arsenic resistance protein 2), bound to the cap binding complex (CBC). We report here that alternative splicing of ARS2's intron 5, generates cytoplasmic isoforms that lack 270 amino acids from the N-terminal of the protein and are functionally distinct from nuclear ARS2. Switching of ARS2 isoforms within the CBC in the cytoplasm has dramatic functional consequences, changing ARS2 from a NMD inhibitor to a NMD promoter that enhances the binding of UPF1 to NCBP1 and ERF1, favouring SURF complex formation, SMG7 recruitment and transcript degradation. ARS2 isoform exchange is also relevant during arsenic stress, where cytoplasmic ARS2 promotes a global response to arsenic in a CBC-independent manner. We propose that ARS2 isoform switching promotes the proper recruitment of RNP complexes during NMD and the cellular response to arsenic stress. The existence of non-redundant ARS2 isoforms is relevant for cell homeostasis, and stress response.

Hfq, a bacterial RNA chaperone, stabilizes small regulatory RNAs (sRNAs) and facilitates sRNA base-pairing with target mRNAs. Hfq has a conserved N-terminal domain and a poorly conserved disordered C-terminal domain (CTD). In a transcriptome-wide examination of the effects of a chromosomal CTD deletion (Hfq1-65), the Escherichia coli mutant was most defective for the accumulation of sRNAs that bind the proximal and distal faces of Hfq (Class II sRNAs), but other sRNAs also were affected. There were only modest effects on the levels of mRNAs, suggesting little disruption of sRNA-dependent regulation. However, cells expressing Hfq lacking the CTD in combination with a weak distal face mutation were defective for the function of the Class II sRNA ChiX and repression of mutS, both dependent upon distal face RNA binding. Loss of the region between amino acids 66-72 was critical for this defect. The CTD region beyond amino acid 72 was not necessary for distal face-dependent regulation, but was needed for functions associated with the Hfq rim, seen most clearly in combination with a rim mutant. Our results suggest that the C-terminus collaborates in various ways with different binding faces of Hfq, leading to distinct outcomes for individual sRNAs.

Operons are a hallmark of the genomic and regulatory architecture of prokaryotes. However, the mechanism by which two genes placed far apart gradually come close and form operons remains to be elucidated. Here, we propose a new model of the origin of operons: Mobile genetic elements called insertion sequences can facilitate the formation of operons by consecutive insertion-deletion-excision reactions. This mechanism barely leaves traces of insertion sequences and thus difficult to detect in nature. In this study, as a proof-of-concept, we reproducibly demonstrated operon formation in the laboratory. The insertion sequence IS3 and the insertion sequence excision enhancer are genes found in a broad range of bacterial species. We introduced these genes into insertion sequence-less Escherichia coli and found that, supporting our hypothesis, the activity of the two genes altered the expression of genes surrounding IS3, closed a 2.7 kb gap between a pair of genes, and formed new operons. This study shows how insertion sequences can facilitate the rapid formation of operons through locally increasing the structural mutation rates and highlights how coevolution with mobile elements may shape the organization of prokaryotic genomes and gene regulation.

Proline tRNA 3'-maturation in Escherichia coli occurs through a one-step RNase E endonucleolytic cleavage immediately after the CCA determinant. This processing pathway is distinct from the 3'-end maturation of the other tRNAs by avoiding the widespread use of 3' → 5' exonucleolytic processing, 3'-polyadenylation and subsequent degradation. Here, we show that the cytosine (C) at the mature 5'-terminus of the proK and proL tRNAs is required for both the RNase E cleavage immediately after the CCA determinant and their functionality. Thus, changing the C nucleotide at the mature 5'-terminus of the proL and proK tRNAs to the more common G nucleotide led to RNase E cleavages 1-4 nucleotides downstream of the CCA determinant. Furthermore, the 5'-modified mutant tRNAs required RNase T and RNase PH for their 3'-maturation and became substrates for polyadenylation and degradation. Strikingly, the aminoacylation of the 5'-modified proline tRNAs was blocked due to the change in the recognition element for prolyl-tRNA-synthetase. An analogous modification of the pheV 5'-mature terminus from G to C nucleotide did not support cell viability. This result provides additional support for the importance of first nucleotide of the mature tRNAs in their processing and functionality.

Turnover of nucleoplasmic transcripts by the mammalian multi-subunit RNA exosome is mediated by two adaptors: the Nuclear EXosome Targeting (NEXT) complex and the Poly(A) tail eXosome Targeting (PAXT) connection. Functional analyses of NEXT and PAXT have largely utilized long-term factor depletion strategies, facilitating the appearance of indirect phenotypes. Here, we rapidly deplete NEXT, PAXT and core exosome components, uncovering the direct consequences of their acute losses. Generally, proteome changes are sparse and largely dominated by co-depletion of other exosome and adaptor subunits, reflecting possible subcomplex compositions. While parallel high-resolution 3' end sequencing of newly synthesized RNA confirms previously established factor specificities, it concomitantly demonstrates an inflation of long-term depletion datasets by secondary effects. Most strikingly, a general intron degradation phenotype, observed in long-term NEXT depletion samples, is undetectable upon short-term depletion, which instead emphasizes NEXT targeting of snoRNA-hosting introns. Further analysis of these introns uncovers an unusual mode of core exosome-independent RNA decay. Our study highlights the accumulation of RNAs as an indirect result of long-term decay factor depletion, which we speculate is, at least partly, due to the exhaustion of alternative RNA decay pathways.

Plasmids are mobile genetic elements that play a key role in microbial ecology and evolution by mediating horizontal transfer of important genes, such as antimicrobial resistance genes. Many microbial genomes have been sequenced by short read sequencers and have resulted in a mix of contigs that derive from plasmids or chromosomes. New tools that accurately identify plasmids are needed to elucidate new plasmid-borne genes of high biological importance. We have developed Deeplasmid, a deep learning tool for distinguishing plasmids from bacterial chromosomes based on the DNA sequence and its encoded biological data. It requires as input only assembled sequences generated by any sequencing platform and assembly algorithm and its runtime scales linearly with the number of assembled sequences. Deeplasmid achieves an AUC-ROC of over 89%, and it was more accurate than five other plasmid classification methods. Finally, as a proof of concept, we used Deeplasmid to predict new plasmids in the fish pathogen Yersinia ruckeri ATCC 29473 that has no annotated plasmids. Deeplasmid predicted with high reliability that a long assembled contig is part of a plasmid. Using long read sequencing we indeed validated the existence of a 102 kb long plasmid, demonstrating Deeplasmid's ability to detect novel plasmids.

Single-stranded genomic DNA can fold into G-quadruplex (G4) structures or form DNA:RNA hybrids (R loops). Recent evidence suggests that such non-canonical DNA structures affect gene expression, DNA methylation, replication fork progression and genome stability. When and how G4 structures form and are resolved remains unclear. Here we report the use of Cleavage Under Targets and Tagmentation (CUT&Tag) for mapping native G4 in mammalian cell lines at high resolution and low background. Mild native conditions used for the procedure retain more G4 structures and provide a higher signal-to-noise ratio than ChIP-based methods. We determine the G4 landscape of mouse embryonic stem cells (ESC), observing widespread G4 formation at active promoters, active and poised enhancers. We discover that the presence of G4 motifs and G4 structures distinguishes active and primed enhancers in mouse ESCs. Upon differentiation to neural progenitor cells (NPC), enhancer G4s are lost. Further, performing R-loop CUT&Tag, we demonstrate the genome-wide co-occurrence of single-stranded DNA, G4s and R loops at promoters and enhancers. We confirm that G4 structures exist independent of ongoing transcription, suggesting an intricate relationship between transcription and non-canonical DNA structures.

DNA nanostructures are a promising tool to deliver molecular payloads to cells. DNA origami structures, where long single-stranded DNA is folded into a compact nanostructure, present an attractive approach to package genes; however, effective delivery of genetic material into cell nuclei has remained a critical challenge. Here, we describe the use of DNA nanostructures encoding an intact human gene and a fluorescent protein encoding gene as compact templates for gene integration by CRISPR-mediated homology-directed repair (HDR). Our design includes CRISPR-Cas9 ribonucleoprotein binding sites on DNA nanostructures to increase shuttling into the nucleus. We demonstrate efficient shuttling and genomic integration of DNA nanostructures using transfection and electroporation. These nanostructured templates display lower toxicity and higher insertion efficiency compared to unstructured double-stranded DNA templates in human primary cells. Furthermore, our study validates virus-like particles as an efficient method of DNA nanostructure delivery, opening the possibility of delivering nanostructures in vivo to specific cell types. Together, these results provide new approaches to gene delivery with DNA nanostructures and establish their use as HDR templates, exploiting both their design features and their ability to encode genetic information. This work also opens a door to translate other DNA nanodevice functions, such as biosensing, into cell nuclei.

Canonical eukaryotic mRNA translation requires 5'cap recognition by initiation factor 4E (eIF4E). In contrast, many positive-strand RNA virus genomes lack a 5'cap and promote translation by non-canonical mechanisms. Among plant viruses, PTEs are a major class of cap-independent translation enhancers located in/near the 3'UTR that recruit eIF4E to greatly enhance viral translation. Previous work proposed a single form of PTE characterized by a Y-shaped secondary structure with two terminal stem-loops (SL1 and SL2) atop a supporting stem containing a large, G-rich asymmetric loop that forms an essential pseudoknot (PK) involving C/U residues located between SL1 and SL2. We found that PTEs with less than three consecutive cytidylates available for PK formation have an upstream stem-loop that forms a kissing loop interaction with the apical loop of SL2, important for formation/stabilization of PK. PKs found in both subclasses of PTE assume a specific conformation with a hyperreactive guanylate (G*) in SHAPE structure probing, previously found critical for binding eIF4E. While PTE PKs were proposed to be formed by Watson-Crick base-pairing, alternative chemical probing and 3D modeling indicate that the Watson-Crick faces of G* and an adjacent guanylate have high solvent accessibilities. Thus, PTE PKs are likely composed primarily of non-canonical interactions.

CRISPR (Clustered Regularly Interspaced Short Palindromic Repeat) technology is a powerful tool in biology and medicine. However, the safety and application of this technology is hampered by excessive activity of CRISPR machinery. It is particularly important to develop methods for switching off CRISPR activity in human cells. The current study demonstrates the concept of supramolecular CRISPR-OFF switches by employing host-guest chemistry. We demonstrate that the CRISPR systems show considerable tolerance to adamantoylation on guide RNAs (gRNAs), whereas supramolecular complexation tremendously affects the function of adamantoyl gRNAs. Host-guest chemistry is demonstrated to be novel and effective tools to reduce unwanted excessive activities of CRISPR complexes in human cells. This work indicates considerable potential of supramolecular strategy for controlling and enhancing CRISPR systems.

In mammalian cells, small non-coding RNAs (sncRNAs) negatively regulate gene expression in a pathway known as RNA interference (RNAi). RNAi can be categorized into post-transcriptional gene silencing (PTGS), which involves the cleavage of target messenger RNA (mRNA) or inhibition of translation in the cytoplasm, and transcriptional gene silencing (TGS), which is mediated by the establishment of repressive epigenetic marks at target loci. Transfer RNAs (tRNAs), which are essential for translation, can be processed into small ncRNAs, termed tRNA-derived small RNAs (tsRNAs). The biogenesis of tsRNAs and their role in gene expression regulation has not yet been fully understood. Here, we show that Dicer dependent tsRNAs promote gene silencing through a mechanism distinct from PTGS and TGS. tsRNAs can lead to downregulation of target genes by targeting introns via nascent RNA silencing (NRS) in nuclei. Furthermore, we show that Ago2 slicer activity is required for this mechanism. Synthetic tsRNAs can significantly reduce expression of a target gene at both RNA and protein levels. Target genes regulated by NRS are associated with various diseases, which further underpins its biological significance. Finally, we show that NRS is evolutionarily conserved and has the potential to be explored as a novel synthetic sRNA based therapeutic.

During the course of the COVID-19 pandemic, large-scale genome sequencing of SARS-CoV-2 has been useful in tracking its spread and in identifying variants of concern (VOC). Viral and host factors could contribute to variability within a host that can be captured in next-generation sequencing reads as intra-host single nucleotide variations (iSNVs). Analysing 1347 samples collected till June 2020, we recorded 16 410 iSNV sites throughout the SARS-CoV-2 genome. We found ∼42% of the iSNV sites to be reported as SNVs by 30 September 2020 in consensus sequences submitted to GISAID, which increased to ∼80% by 30th June 2021. Following this, analysis of another set of 1774 samples sequenced in India between November 2020 and May 2021 revealed that majority of the Delta (B.1.617.2) and Kappa (B.1.617.1) lineage-defining variations appeared as iSNVs before getting fixed in the population. Besides, mutations in RdRp as well as RNA-editing by APOBEC and ADAR deaminases seem to contribute to the differential prevalence of iSNVs in hosts. We also observe hyper-variability at functionally critical residues in Spike protein that could alter the antigenicity and may contribute to immune escape. Thus, tracking and functional annotation of iSNVs in ongoing genome surveillance programs could be important for early identification of potential variants of concern and actionable interventions.

Expression of the E3 ligase TRIM21 is increased in a broad spectrum of cancers; however, the functionally relevant molecular pathway targeted by TRIM21 overexpression remains largely unknown. Here, we show that TRIM21 directly interacts with and ubiquitinates CLASPIN, a mediator for ATR-dependent CHK1 activation. TRIM21-mediated K63-linked ubiquitination of CLASPIN counteracts the K6-linked ubiquitination of CLASPIN which is essential for its interaction with TIPIN and subsequent chromatin loading. We further show that overexpression of TRIM21, but not a TRIM21 catalytically inactive mutant, compromises CHK1 activation, leading to replication fork instability and tumorigenesis. Our findings demonstrate that TRIM21 suppresses CHK1 activation by preferentially targeting CLASPIN for K63-linked ubiquitination, providing a potential target for cancer therapy.

Small noncoding RNAs such as piRNAs are guides for Argonaute proteins, enabling sequence-specific, post-transcriptional regulation of gene expression. The piRNAs of Caenorhabditis elegans have been observed to bind targets with high mismatch tolerance and appear to lack specific transposon targets, unlike piRNAs in Drosophila melanogaster and other organisms. These observations support a model in which C. elegans piRNAs provide a broad, indiscriminate net of silencing, competing with siRNAs associated with the CSR-1 Argonaute that specifically protect self-genes from silencing. However, the breadth of piRNA targeting has not been subject to in-depth quantitative analysis, nor has it been explained how piRNAs are distributed across sequence space to achieve complete coverage. Through a bioinformatic analysis of piRNA sequences, incorporating an original data-based metric of piRNA-target distance, we demonstrate that C. elegans piRNAs are functionally random, in that their coverage of sequence space is comparable to that of random sequences. By possessing a sufficient number of distinct, essentially random piRNAs, C. elegans is able to target arbitrary nonself sequences with high probability. We extend this approach to a selection of other nematodes, finding results which elucidate the mechanism by which nonself mRNAs are silenced, and have implications for piRNA evolution and biogenesis.

Transcriptional regulation in response to thyroid hormone (3,5,3'-triiodo-l-thyronine, T3) is a dynamic and cell-type specific process that maintains cellular homeostasis and identity in all tissues. However, our understanding of the mechanisms of thyroid hormone receptor (TR) actions at the molecular level are actively being refined. We used an integrated genomics approach to profile and characterize the cistrome of TRβ, map changes in chromatin accessibility, and capture the transcriptomic changes in response to T3 in normal human thyroid cells. There are significant shifts in TRβ genomic occupancy in response to T3, which are associated with differential chromatin accessibility, and differential recruitment of SWI/SNF chromatin remodelers. We further demonstrate selective recruitment of BAF and PBAF SWI/SNF complexes to TRβ binding sites, revealing novel differential functions in regulating chromatin accessibility and gene expression. Our findings highlight three distinct modes of TRβ interaction with chromatin and coordination of coregulator activity.

Toxin-antitoxin (TA) systems are proposed to play crucial roles in bacterial growth under stress conditions such as phage infection. The type III TA systems consist of a protein toxin whose activity is inhibited by a noncoding RNA antitoxin. The toxin is an endoribonuclease, while the antitoxin consists of multiple repeats of RNA. The toxin assembles with the individual antitoxin repeats into a cyclic complex in which the antitoxin forms a pseudoknot structure. While structure and functions of some type III TA systems are characterized, the complex assembly process is not well understood. Using bioinformatics analysis, we have identified type III TA systems belonging to the ToxIN family across different Escherichia coli strains and found them to be clustered into at least five distinct clusters. Furthermore, we report a 2.097 Å resolution crystal structure of the first E. coli ToxIN complex that revealed the overall assembly of the protein-RNA complex. Isothermal titration calorimetry experiments showed that toxin forms a high-affinity complex with antitoxin RNA resulting from two independent (5' and 3' sides of RNA) RNA binding sites on the protein. These results further our understanding of the assembly of type III TA complexes in bacteria.

The HIV-1 capsid core participates in several replication processes. The mature capsid core is a lattice composed of capsid (CA) monomers thought to assemble first into CA dimers, then into ∼250 CA hexamers and 12 CA pentamers. CA assembly requires conformational flexibility of each unit, resulting in the presence of unique, solvent-accessible surfaces. Significant advances have improved our understanding of the roles of the capsid core in replication; however, the contributions of individual CA assembly forms remain unclear and there are limited tools available to evaluate these forms in vivo. Here, we have selected aptamers that bind CA lattice tubes. We describe aptamer CA15-2, which selectively binds CA lattice, but not CA monomer or CA hexamer, suggesting that it targets an interface present and accessible only on CA lattice. CA15-2 does not compete with PF74 for binding, indicating that it likely binds a non-overlapping site. Furthermore, CA15-2 inhibits HIV-1 replication when expressed in virus producer cells, but not target cells, suggesting that it binds a biologically-relevant site during virus production that is either not accessible during post-entry replication steps or is accessible but unaltered by aptamer binding. Importantly, CA15-2 represents the first aptamer that specifically recognizes the HIV-1 CA lattice.

Multiple RNA polymerases (RNAPs) transcribing a gene have been known to exhibit collective group behavior, causing the transcription elongation rate to increase with the rate of transcription initiation. Such behavior has long been believed to be driven by a physical interaction or 'push' between closely spaced RNAPs. However, recent studies have posited that RNAPs separated by longer distances may cooperate by modifying the DNA segment under transcription. Here, we present a theoretical model incorporating the mechanical coupling between RNAP translocation and the DNA torsional response. Using stochastic simulations, we demonstrate DNA supercoiling-mediated long-range cooperation between co-transcribing RNAPs. We find that inhibiting transcription initiation can slow down the already recruited RNAPs, in agreement with recent experimental observations, and predict that the average transcription elongation rate varies non-monotonically with the rate of transcription initiation. We further show that while RNAPs transcribing neighboring genes oriented in tandem can cooperate, those transcribing genes in divergent or convergent orientations can act antagonistically, and that such behavior holds over a large range of intergenic separations. Our model makes testable predictions, revealing how the mechanical interplay between RNAPs and the DNA they transcribe can govern transcriptional dynamics.

Information about the cellular concentrations of deoxyribonucleoside triphosphates (dNTPs) is instrumental for mechanistic studies of DNA replication and for understanding diseases caused by defects in dNTP metabolism. The dNTPs are measured by methods based on either HPLC or DNA polymerization. An advantage with the HPLC-based techniques is that the parallel analysis of ribonucleoside triphosphates (rNTPs) can serve as an internal quality control of nucleotide integrity and extraction efficiency. We have developed a Freon-free trichloroacetic acid-based method to extract cellular nucleotides and an isocratic reverse phase HPLC-based technique that is able to separate dNTPs, rNTPs and ADP in a single run. The ability to measure the ADP levels improves the control of nucleotide integrity, and the use of an isocratic elution overcomes the shifting baseline problems in previously developed gradient-based reversed phase protocols for simultaneously measuring dNTPs and rNTPs. An optional DNA-polymerase-dependent step is used for confirmation that the dNTP peaks do not overlap with other components of the extracts, further increasing the reliability of the analysis. The method is compatible with a wide range of biological samples and has a sensitivity better than other UV-based HPLC protocols, closely matching that of mass spectrometry-based detection.

DNA chemical modifications, including methylation, are widespread and play important roles in prokaryotes and viruses. However, current knowledge of these modification systems is severely biased towards a limited number of culturable prokaryotes, despite the fact that a vast majority of microorganisms have not yet been cultured. Here, using single-molecule real-time sequencing, we conducted culture-independent 'metaepigenomic' analyses (an integrated analysis of metagenomics and epigenomics) of marine microbial communities. A total of 233 and 163 metagenomic-assembled genomes (MAGs) were constructed from diverse prokaryotes and viruses, respectively, and 220 modified motifs and 276 DNA methyltransferases (MTases) were identified. Most of the MTase genes were not genetically linked with the endonuclease genes predicted to be involved in defense mechanisms against extracellular DNA. The MTase-motif correspondence found in the MAGs revealed 10 novel pairs, 5 of which showed novel specificities and experimentally confirmed the catalytic specificities of the MTases. We revealed novel alternative specificities in MTases that are highly conserved in Alphaproteobacteria, which may enhance our understanding of the co-evolutionary history of the methylation systems and the genomes. Our findings highlight diverse unexplored DNA modifications that potentially affect the ecology and evolution of prokaryotes and viruses in nature.