tRNAs are key partners in ribosome-dependent protein synthesis. This process is highly dependent on the fidelity of tRNA aminoacylation by aminoacyl-tRNA synthetases and relies primarily on sets of identities within tRNA molecules composed of determinants and antideterminants preventing mischarging by non-cognate synthetases. Such identity sets were discovered in the tRNAs of a few model organisms, and their properties were generalized as universal identity rules. Since then, the panel of identity elements governing the accuracy of tRNA aminoacylation has expanded considerably, but the increasing number of reported functional idiosyncrasies has led to some confusion. In parallel, the description of other processes involving tRNAs, often well beyond aminoacylation, has progressed considerably, greatly expanding their interactome and uncovering multiple novel identities on the same tRNA molecule. This review highlights key findings on the mechanistics and evolution of tRNA and tRNA-like identities. In addition, new methods and their results for searching sets of multiple identities on a single tRNA are discussed. Taken together, this knowledge shows that a comprehensive understanding of the functional role of individual and collective nucleotide identity sets in tRNA molecules is needed for medical, biotechnological and other applications.

The intestinal epithelial regeneration is driven by intestinal stem cells under homeostatic conditions. Differentiated intestinal epithelial cells, such as Paneth cells, are capable of acquiring multipotency and contributing to regeneration upon the loss of intestinal stem cells. Paneth cells also support intestinal stem cell survival and regeneration. We report here that depletion of an RNA-binding protein named polypyrimidine tract binding protein 1 (PTBP1) in mouse intestinal epithelial cells causes intestinal stem cell death and epithelial regeneration failure. Mechanistically, we show that PTBP1 inhibits neuronal-like splicing programs in intestinal crypt cells, which is critical for maintaining intestinal stem cell stemness. This function is achieved at least in part through promoting the non-productive splicing of its paralog PTBP2. Moreover, PTBP1 inhibits the expression of an AKT inhibitor PHLDA3 in Paneth cells and permits AKT activation, which presumably maintains Paneth cell plasticity and function in supporting intestinal stem cell niche. We show that PTBP1 directly binds to a CU-rich region in the 3' UTR of Phlda3, which we demonstrate to be critical for downregulating the mRNA and protein levels of Phlda3. Our results thus reveal the multifaceted in vivo regulation of intestinal epithelial regeneration by PTBP1 at the post-transcriptional level.

LINE-1 retrotransposons are sequences capable of copying themselves to new genomic loci via an RNA intermediate. New studies implicate LINE-1 in a range of diseases, especially in the context of aging, but without an accurate understanding of where and when LINE-1 is expressed, a full accounting of its role in health and disease is not possible. We therefore developed a method-5' scL1seq-that makes use of a widely available library preparation method (10x Genomics 5' single cell RNA-seq) to measure LINE-1 expression in tens of thousands of single cells. We recapitulated the known pattern of LINE-1 expression in tumors-present in cancer cells, absent from immune cells-and identified hitherto undescribed LINE-1 expression in human epithelial cells and mouse hippocampal neurons. In both cases, we saw a modest increase with age, supporting recent research connecting LINE-1 to age related diseases.

Anti-CRISPR proteins are encoded by phages to inhibit the CRISPR-Cas systems of the hosts. AcrIIC5 inhibits several naturally high-fidelity type II-C Cas9 enzymes, including orthologs from Neisseria meningitidis (Nme1Cas9) and Simonsiella muelleri (SmuCas9). Here, we solve the structure of AcrIIC5 in complex with Nme1Cas9 and sgRNA. We show that AcrIIC5 adopts a novel fold to mimic the size and charge distribution of double-stranded DNA, and uses its negatively charged grooves to bind and occlude the protospacer adjacent motif (PAM) binding site in the target DNA cleft of Cas9. AcrIIC5 is positioned into the crevice between the WED and PI domains of Cas9, and one end of the anti-CRISPR interacts with the phosphate lock loop and a linker between the RuvC and BH domains. We employ biochemical and mutational analyses to build a model for AcrIIC5's mechanism of action, and identify residues on both the anti-CRISPR and Cas9 that are important for their interaction and inhibition. Together, the structure and mechanism of AcrIIC5 reveal convergent evolution among disparate anti-CRISPR proteins that use a DNA-mimic strategy to inhibit diverse CRISPR-Cas surveillance complexes, and provide new insights into a tool for potent inhibition of type II-C Cas9 orthologs.

Pseudomonas aeruginosa is capable of causing acute and chronic infections in various host tissues, which depends on its abilities to effectively utilize host-derived nutrients and produce protein virulence factors and toxic compounds. However, the regulatory mechanisms that direct metabolic intermediates towards production of toxic compounds are poorly understood. We previously identified a regulatory protein PvrA that controls genes involved in fatty acid catabolism by binding to palmitoyl-coenzyme A (CoA). In this study, transcriptomic analyses revealed that PvrA activates the Pseudomonas quinolone signal (PQS) synthesis genes, while suppressing genes for production of polyhydroxyalkanoates (PHAs). When palmitic acid was the sole carbon source, mutation of pvrA reduced production of pyocyanin and rhamnolipids due to defective PQS synthesis, but increased PHA production. We further solved the co-crystal structure of PvrA with palmitoyl-CoA and identified palmitoyl-CoA-binding residues. By using pvrA mutants, we verified the roles of the key palmitoyl-CoA-binding residues in gene regulation in response to palmitic acid. Since the PQS signal molecules, rhamnolipids and PHA synthesis pathways are interconnected by common metabolic intermediates, our results revealed a regulatory mechanism that directs carbon flux from carbon/energy storage to virulence factor production, which might be crucial for the pathogenesis.

Many replicative DNA polymerases couple DNA replication and unwinding activities to perform strand displacement DNA synthesis, a critical ability for DNA metabolism. Strand displacement is tightly regulated by partner proteins, such as single-stranded DNA (ssDNA) binding proteins (SSBs) by a poorly understood mechanism. Here, we use single-molecule optical tweezers and biochemical assays to elucidate the molecular mechanism of strand displacement DNA synthesis by the human mitochondrial DNA polymerase, Polγ, and its modulation by cognate and noncognate SSBs. We show that Polγ exhibits a robust DNA unwinding mechanism, which entails lowering the energy barrier for unwinding of the first base pair of the DNA fork junction, by ∼55%. However, the polymerase cannot prevent the reannealing of the parental strands efficiently, which limits by ∼30-fold its strand displacement activity. We demonstrate that SSBs stimulate the Polγ strand displacement activity through several mechanisms. SSB binding energy to ssDNA additionally increases the destabilization energy at the DNA junction, by ∼25%. Furthermore, SSB interactions with the displaced ssDNA reduce the DNA fork reannealing pressure on Polγ, in turn promoting the productive polymerization state by ∼3-fold. These stimulatory effects are enhanced by species-specific functional interactions and have significant implications in the replication of the human mitochondrial DNA.

The potential of synthetic mRNA as a genetic carrier has increased its application in scientific fields. Because the 5' cap regulates the stability and translational activity of mRNAs, there are concerted efforts to search for and synthesize chemically-modified 5' caps that improve the functionality of mRNA. Here, we report an easy and efficient method to synthesize functional mRNAs by modifying multiple 5' cap analogs using a vaccinia virus-capping enzyme. We show that this enzyme can introduce a variety of GTP analogs to the 5' end of RNA to generate 5' cap-modified mRNAs that exhibit different translation levels. Notably, some of these modified mRNAs improve translation efficiency and can be conjugated to chemical structures, further increasing their functionality. Our versatile method to generate 5' cap-modified mRNAs will provide useful tools for RNA therapeutics and biological research.

Structural, regulatory and enzymatic proteins interact with DNA to maintain a healthy and functional genome. Yet, our structural understanding of how proteins interact with DNA is limited. We present MELD-DNA, a novel computational approach to predict the structures of protein-DNA complexes. The method combines molecular dynamics simulations with general knowledge or experimental information through Bayesian inference. The physical model is sensitive to sequence-dependent properties and conformational changes required for binding, while information accelerates sampling of bound conformations. MELD-DNA can: (i) sample multiple binding modes; (ii) identify the preferred binding mode from the ensembles; and (iii) provide qualitative binding preferences between DNA sequences. We first assess performance on a dataset of 15 protein-DNA complexes and compare it with state-of-the-art methodologies. Furthermore, for three selected complexes, we show sequence dependence effects of binding in MELD predictions. We expect that the results presented herein, together with the freely available software, will impact structural biology (by complementing DNA structural databases) and molecular recognition (by bringing new insights into aspects governing protein-DNA interactions).

X chromosome inactivation (XCI) is an essential process, yet it initiates with remarkable diversity in various mammalian species. XIST, the main trigger of XCI, is controlled in the mouse by an interplay of lncRNA genes (LRGs), some of which evolved concomitantly to XIST and have orthologues across all placental mammals. Here, we addressed the functional conservation of human orthologues of two such LRGs, FTX and JPX. By combining analysis of single-cell RNA-seq data from early human embryogenesis with various functional assays in matched human and mouse pluripotent stem- or differentiated post-XCI cells, we demonstrate major functional differences for these orthologues between species, independently of primary sequence conservation. While the function of FTX is not conserved in humans, JPX stands as a major regulator of XIST expression in both species. However, we show that different entities of JPX control the production of XIST at various steps depending on the species. Altogether, our study highlights the functional versatility of LRGs across evolution, and reveals that functional conservation of orthologous LRGs may involve diversified mechanisms of action. These findings represent a striking example of how the evolvability of LRGs can provide adaptative flexibility to constrained gene regulatory networks.

The clustered regularly interspaced short palindromic repeats (CRISPR) Cas system is a powerful tool that has the potential to become a therapeutic gene editor in the near future. Cas9 is the best studied CRISPR system and has been shown to have problems that restrict its use in therapeutic applications. Chromatin structure is a known impactor of Cas9 targeting and there is a gap in knowledge on Cas9's efficacy when targeting such locations. To quantify at a single base pair resolution how chromatin inhibits on-target gene editing relative to off-target editing of exposed mismatching targets, we developed the gene editor mismatch nucleosome inhibition assay (GEMiNI-seq). GEMiNI-seq utilizes a library of nucleosome sequences to examine all target locations throughout nucleosomes in a single assay. The results from GEMiNI-seq revealed that the location of the protospacer-adjacent motif (PAM) sequence on the nucleosome edge drives the ability for Cas9 to access its target sequence. In addition, Cas9 had a higher affinity for exposed mismatched targets than on-target sequences within a nucleosome. Overall, our results show how chromatin structure impacts the fidelity of Cas9 to potential targets and highlight how targeting sequences with exposed PAMs could limit off-target gene editing, with such considerations improving Cas9 efficacy and resolving current limitations.

It is widely accepted that pooled library CRISPR knockout screens offer greater sensitivity and specificity than prior technologies in detecting genes whose disruption leads to fitness defects, a critical step in identifying candidate cancer targets. However, the assumption that CRISPR screens are saturating has been largely untested. Through integrated analysis of screen data in cancer cell lines generated by the Cancer Dependency Map, we show that a typical CRISPR screen has a ∼20% false negative rate, in addition to library-specific false negatives. Replicability falls sharply as gene expression decreases, while cancer subtype-specific genes within a tissue show distinct profiles compared to false negatives. Cumulative analyses across tissues improves our understanding of core essential genes and suggest only a small number of lineage-specific essential genes, enriched for transcription factors that define pathways of tissue differentiation. To recover false negatives, we introduce a method, Joint Log Odds of Essentiality (JLOE), which builds on our prior work with BAGEL to selectively rescue the false negatives without an increased false discovery rate.

Discovery of oxidative catalysis with G-quadruplex•hemin constructs prompted a range of exciting developments in the field of biosensor design. Thus, G-quadruplex based DNAzymes with peroxidase activity found a niche as signal transduction modules in a wide range of analytical applications. The ability of nucleic acid scaffolds to recognise a variety of practically meaningful markers and to translate the recognition events into conformational changes powers numerous sensor design possibilities. In this work, we establish a catalase activity of G-quadruplex•hemin scaffolds. Catalase activated hydrogen peroxide decomposition generates molecular oxygen that forms bubbles. Observation of bubbles is a truly equipment free signal readout platform that is highly desirable in limited resources or do-it-yourself environments. We take a preliminary insight into a G-quadruplex structure-folding topology-catalase activity correlation and establish efficient operating conditions. Further, we demonstrate the platform's potential as a signal transduction modality for reporting on biomolecular recognition using an oligonucleotide as a proof-of-concept target. Ultimately, activatable catalases based on G-quadruplex•hemin scaffolds promise to become valuable contributors towards accessible molecular diagnostics applications.

During oocyte development in mice, transcripts accumulate in the growth phase and are subsequently degraded during maturation. At the transition point between growth and maturation, oocytes have an intact nucleus or germinal vesicle (GV), and terminal uridylation labels RNA for degradation in meiosis I. By profiling the transcriptome using single-oocyte long-read PacBio RNA sequencing, we document that a small cohort of mRNAs are polyadenylated after terminal uridylation in GV oocytes [designated uridylated-poly(A) RNA]. Because DIS3L2 ribonuclease is known to degrade uridylated transcripts, we established oocyte-specific Dis3l2 knockout mice (Dis3l2cKO). Upon DIS3L2 depletion, uridylated-poly(A) RNAs remain intact which increases their abundance, and they predominate in the transcriptome of Dis3l2cKO oocytes. The abundance of uridylated-poly(A) RNA in Dis3l2cKO oocytes arises not only from insufficient degradation, but also from the stabilizing effect of subsequent polyadenylation. Uridylated-poly(A) RNAs have shorter poly(A) tails and their translation activity decreases in Dis3l2cKO oocytes. Almost all Dis3l2cKO oocytes arrest at the GV stage, and female mice are infertile. Our study demonstrates multiple fates for RNA after terminal uridylation and highlights the role of DIS3L2 ribonuclease in safeguarding the transcriptome and ensuring female fertility.

Ribosomes of Bacteroidia (formerly Bacteroidetes) fail to recognize Shine-Dalgarno (SD) sequences even though they harbor the anti-SD (ASD) of 16S rRNA. Inhibition of SD-ASD pairing is due to sequestration of the 3' tail of 16S rRNA in a pocket formed by bS21, bS18, and bS6 on the 30S platform. Interestingly, in many Flavobacteriales, the gene encoding bS21, rpsU, contains an extended SD sequence. In this work, we present genetic and biochemical evidence that bS21 synthesis in Flavobacterium johnsoniae is autoregulated via a subpopulation of ribosomes that specifically lack bS21. Mutation or depletion of bS21 in the cell increases translation of reporters with strong SD sequences, such as rpsU'-gfp, but has no effect on other reporters. Purified ribosomes lacking bS21 (or its C-terminal region) exhibit higher rates of initiation on rpsU mRNA and lower rates of initiation on other (SD-less) mRNAs than control ribosomes. The mechanism of autoregulation depends on extensive pairing between mRNA and 16S rRNA, and exceptionally strong SD sequences, with predicted pairing free energies of < -13 kcal/mol, are characteristic of rpsU across the Bacteroidota. This work uncovers a clear example of specialized ribosomes in bacteria.

Positive transcription elongation factor b (P-TEFb) is the crucial player in RNA polymerase II (Pol II) pause release that has emerged as a promising target in cancer. Because single-agent therapy may fail to deliver durable clinical response, targeting of P-TEFb shall benefit when deployed as a combination therapy. We screened a comprehensive oncology library and identified clinically relevant antimetabolites and Mouse double minute 2 homolog (MDM2) inhibitors as top compounds eliciting p53-dependent death of colorectal cancer cells in synergy with selective inhibitors of P-TEFb. While the targeting of P-TEFb augments apoptosis by anti-metabolite 5-fluorouracil, it switches the fate of cancer cells by the non-genotoxic MDM2 inhibitor Nutlin-3a from cell-cycle arrest to apoptosis. Mechanistically, the fate switching is enabled by the induction of p53-dependent pro-apoptotic genes and repression of P-TEFb-dependent pro-survival genes of the PI3K-AKT signaling cascade, which stimulates caspase 9 and intrinsic apoptosis pathway in BAX/BAK-dependent manner. Finally, combination treatments trigger apoptosis of cancer cell spheroids. Together, co-targeting of P-TEFb and suppressors of intrinsic apoptosis could become a viable strategy to eliminate cancer cells.

Promoter-proximal Polymerase II (Pol II) pausing is a key rate-limiting step for gene expression. DNA and RNA-binding trans-acting factors regulating the extent of pausing have been identified. However, we lack a quantitative model of how interactions of these factors determine pausing, therefore the relative importance of implicated factors is unknown. Moreover, previously unknown regulators might exist. Here we address this gap with a machine learning model that accurately predicts the extent of promoter-proximal Pol II pausing from large-scale genome and transcriptome binding maps and gene annotation and sequence composition features. We demonstrate high accuracy and generalizability of the model by validation on an independent cell line which reveals the model's cell line agnostic character. Model interpretation in light of prior knowledge about molecular functions of regulatory factors confirms the interconnection of pausing with other RNA processing steps. Harnessing underlying feature contributions, we assess the relative importance of each factor, quantify their predictive effects and systematically identify previously unknown regulators of pausing. We additionally identify 16 previously unknown 7SK ncRNA interacting RNA-binding proteins predictive of pausing. Our work provides a framework to further our understanding of the regulation of the critical early steps in transcriptional elongation.

Overexpression of androgen receptor (AR) is the primary cause of castration-resistant prostate cancer, although mechanisms upregulating AR transcription in this context are not well understood. Our RNA-seq studies revealed that SMAD3 knockdown decreased levels of AR and AR target genes, whereas SMAD4 or SMAD2 knockdown had little or no effect. ChIP-seq analysis showed that SMAD3 knockdown decreased global binding of AR to chromatin. Mechanistically, we show that SMAD3 binds to intron 3 of the AR gene to promote AR expression. Targeting these binding sites by CRISPRi reduced transcript levels of AR and AR targets. In addition, ∼50% of AR and SMAD3 ChIP-seq peaks overlapped, and SMAD3 may also cooperate with or co-activate AR for AR target expression. Functionally, AR re-expression in SMAD3-knockdown cells partially rescued AR target expression and cell growth defects. The SMAD3 peak in AR intron 3 overlapped with H3K27ac ChIP-seq and ATAC-seq peaks in datasets of prostate cancer. AR and SMAD3 mRNAs were upregulated in datasets of metastatic prostate cancer and CRPC compared with primary prostate cancer. A SMAD3 PROTAC inhibitor reduced levels of AR, AR-V7 and AR targets in prostate cancer cells. This study suggests that SMAD3 could be targeted to inhibit AR in prostate cancer.

Cell cycle regulation is crucial for all living organisms and is often targeted by viruses to facilitate their own propagation, yet cell cycle progression control is largely underexplored in archaea. In this work, we reveal a cell cycle regulator (aCcr1) carrying a ribbon-helix-helix (RHH) domain and ubiquitous in the Thermoproteota of the order Sulfolobales and their viruses. Overexpression of several aCcr1 members including gp21 of rudivirus SIRV2 and its host homolog SiL_0190 of Saccharolobus islandicus LAL14/1 results in impairment of cell division, evidenced by growth retardation, cell enlargement and an increase in cellular DNA content. Additionally, both gp21 and SiL_0190 can bind to the motif AGTATTA conserved in the promoter of several genes involved in cell division, DNA replication and cellular metabolism thereby repressing or inducing their transcription. Our results suggest that aCcr1 silences cell division and drives progression to the S-phase in Sulfolobales, a function exploited by viruses to facilitate viral propagation.

Altered eIF4A1 activity promotes translation of highly structured, eIF4A1-dependent oncogene mRNAs at root of oncogenic translational programmes. It remains unclear how these mRNAs recruit and activate eIF4A1 unwinding specifically to facilitate their preferential translation. Here, we show that single-stranded RNA sequence motifs specifically activate eIF4A1 unwinding allowing local RNA structural rearrangement and translation of eIF4A1-dependent mRNAs in cells. Our data demonstrate that eIF4A1-dependent mRNAs contain AG-rich motifs within their 5'UTR which specifically activate eIF4A1 unwinding of local RNA structure to facilitate translation. This mode of eIF4A1 regulation is used by mRNAs encoding components of mTORC-signalling and cell cycle progression, and renders these mRNAs particularly sensitive to eIF4A1-inhibition. Mechanistically, we show that binding of eIF4A1 to AG-rich sequences leads to multimerization of eIF4A1 with eIF4A1 subunits performing distinct enzymatic activities. Our structural data suggest that RNA-binding of multimeric eIF4A1 induces conformational changes in the RNA resulting in an optimal positioning of eIF4A1 proximal to the RNA duplex enabling efficient unwinding. Our data proposes a model in which AG-motifs in the 5'UTR of eIF4A1-dependent mRNAs specifically activate eIF4A1, enabling assembly of the helicase-competent multimeric eIF4A1 complex, and positioning these complexes proximal to stable localised RNA structure allowing ribosomal subunit scanning.

Inefficient endosomal escape remains the primary barrier to the broad application of oligonucleotide therapeutics. Liver uptake after systemic administration is sufficiently robust that a therapeutic effect can be achieved but targeting extrahepatic tissues remains challenging. Prior attempts to improve oligonucleotide activity using small molecules that increase the leakiness of endosomes have failed due to unacceptable toxicity. Here, we show that the well-tolerated and orally bioavailable synthetic sphingolipid analog, SH-BC-893, increases the activity of antisense oligonucleotides (ASOs) and small interfering RNAs (siRNAs) up to 200-fold in vitro without permeabilizing endosomes. SH-BC-893 treatment trapped endocytosed oligonucleotides within extra-lysosomal compartments thought to be more permeable due to frequent membrane fission and fusion events. Simultaneous disruption of ARF6-dependent endocytic recycling and PIKfyve-dependent lysosomal fusion was necessary and sufficient for SH-BC-893 to increase non-lysosomal oligonucleotide levels and enhance their activity. In mice, oral administration of SH-BC-893 increased ASO potency in the liver by 15-fold without toxicity. More importantly, SH-BC-893 enabled target RNA knockdown in the CNS and lungs of mice treated subcutaneously with cholesterol-functionalized duplexed oligonucleotides or unmodified ASOs, respectively. Together, these results establish the feasibility of using a small molecule that disrupts endolysosomal trafficking to improve the activity of oligonucleotides in extrahepatic tissues.

Direct RNA sequencing with a commercial nanopore platform was used to sequence RNA containing uridine (U), pseudouridine (Ψ) or N1-methylpseudouridine (m1Ψ) in >100 different 5-nucleotide contexts. The base calling data for Ψ or m1Ψ were similar but different from U allowing their detection. Understanding the nanopore signatures for Ψ and m1Ψ enabled a running start T7 RNA polymerase assay to study the selection of UTP versus ΨTP or m1ΨTP competing mixtures in all possible adjacent sequence contexts. A significant sequence context dependency was observed for T7 RNA polymerase with insertion yields for ΨTP versus UTP spanning a range of 20-65%, and m1ΨTP versus UTP producing variable yields that differ by 15-70%. Experiments with SP6 RNA polymerase, as well as chemically-modified triphosphates and DNA templates provide insight to explain the observations. The SP6 polymerase introduced m1ΨTP when competed with UTP with a smaller window of yields (15-30%) across all sequence contexts studied. These results may aid in future efforts that employ RNA polymerases to make therapeutic mRNAs with sub-stoichiometric amounts of m1Ψ.

Conjugation of DNA relies on multicomponent protein complexes bridging two bacterial cytoplasmic compartments. Whereas plasmid conjugation systems have been well documented, those of integrative and conjugative elements (ICEs) have remained poorly studied. We characterize here the conjugation system of the ICEclc element in Pseudomonas putida UWC1 that is a model for a widely distributed family of ICEs. By in frame deletion and complementation, we show the importance on ICE transfer of 22 genes in a 20-kb conserved ICE region. Protein comparisons recognized seven homologs to plasmid type IV secretion system components, another six homologs to frequent accessory proteins, and the rest without detectable counterparts. Stationary phase imaging of P. putida ICEclc with in-frame fluorescent protein fusions to predicted type IV components showed transfer-competent cell subpopulations with multiple fluorescent foci, largely overlapping in dual-labeled subcomponents, which is suggestive for multiple conjugation complexes per cell. Cross-dependencies between subcomponents in ICE-type IV secretion system assembly were revealed by quantitative foci image analysis in a variety of ICEclc mutant backgrounds. In conclusion, the ICEclc family presents an evolutionary distinct type IV conjugative system with transfer competent cells specialized in efficient transfer.

Translation is a key determinant of gene expression and an important biotechnological engineering target. In bacteria, 5'-untranslated region (5'-UTR) and coding sequence (CDS) are well-known mRNA parts controlling translation and thus cellular protein levels. However, the complex interaction of 5'-UTR and CDS has so far only been studied for few sequences leading to non-generalisable and partly contradictory conclusions. Herein, we systematically assess the dynamic translation from over 1.2 million 5'-UTR-CDS pairs in Escherichia coli to investigate their collective effect using a new method for ultradeep sequence-function mapping. This allows us to disentangle and precisely quantify effects of various sequence determinants of translation. We find that 5'-UTR and CDS individually account for 53% and 20% of variance in translation, respectively, and show conclusively that, contrary to a common hypothesis, tRNA abundance does not explain expression changes between CDSs with different synonymous codons. Moreover, the obtained large-scale data provide clear experimental evidence for a base-pairing interaction between initiator tRNA and mRNA beyond the anticodon-codon interaction, an effect that is often masked for individual sequences and therefore inaccessible to low-throughput approaches. Our study highlights the indispensability of ultradeep sequence-function mapping to accurately determine the contribution of parts and phenomena involved in gene regulation.

In eukaryotes, many DNA/RNA-binding proteins possess intrinsically disordered regions (IDRs) with large negative charge, some of which involve a consecutive sequence of aspartate (D) or glutamate (E) residues. We refer to them as D/E repeats. The functional role of D/E repeats is not well understood, though some of them are known to cause autoinhibition through intramolecular electrostatic interaction with functional domains. In this work, we investigated the impacts of D/E repeats on the target DNA search kinetics for the high-mobility group box 1 (HMGB1) protein and the artificial protein constructs of the Antp homeodomain fused with D/E repeats of varied lengths. Our experimental data showed that D/E repeats of particular lengths can accelerate the target association in the overwhelming presence of non-functional high-affinity ligands ('decoys'). Our coarse-grained molecular dynamics (CGMD) simulations showed that the autoinhibited proteins can bind to DNA and transition into the uninhibited complex with DNA through an electrostatically driven induced-fit process. In conjunction with the CGMD simulations, our kinetic model can explain how D/E repeats can accelerate the target association process in the presence of decoys. This study illuminates an unprecedented role of the negatively charged IDRs in the target search process.

Dissemination of antibiotic resistance, a current societal challenge, is often driven by horizontal gene transfer through bacterial conjugation. During conjugative plasmid transfer, single-stranded (ss) DNA is transferred from the donor to the recipient cell. Subsequently, a complete double-stranded (ds) plasmid molecule is generated and plasmid-encoded genes are expressed, allowing successful establishment of the transconjugant cell. Such dynamics of transmission can be modulated by host- or plasmid-encoded factors, either in the donor or in the recipient cell. We applied transposon insertion sequencing to identify host-encoded factors that affect conjugative transfer frequency in Escherichia coli. Disruption of the recipient uvrD gene decreased the acquisition frequency of conjugative plasmids belonging to different incompatibility groups. Results from various UvrD mutants suggested that dsDNA binding activity and interaction with RNA polymerase are dispensable, but ATPase activity is required for successful plasmid establishment of transconjugant cells. Live-cell microscopic imaging showed that the newly transferred ssDNA within a uvrD- recipient often failed to be converted to dsDNA. Our work suggested that in addition to its role in maintaining genome integrity, UvrD is also key for the establishment of horizontally acquired plasmid DNA that drives genome diversity and evolution.

Epithelial mesenchymal plasticity (EMP) is a complex cellular reprogramming event that plays a major role in tissue homeostasis. Recently we observed the unfolded protein response (UPR) triggers EMP through the inositol-requiring protein 1 (IRE1α)-X-box-binding protein 1 spliced (XBP1s) axis, enhancing glucose shunting to protein N glycosylation. To better understand the genomic targets of XBP1s, we identified its genomic targets using Cleavage Under Targets and Release Using Nuclease (CUT&RUN) of a FLAG-epitope tagged XBP1s in RSV infection. CUT&RUN identified 7086 binding sites in chromatin that were enriched in AP-1 motifs and GC-sequences. Of these binding sites, XBP1s peaks mapped to 4827 genes controlling Rho-GTPase signaling, N-linked glycosylation and ER-Golgi transport. Strikingly, XBP1s peaks were within 1 kb of transcription start sites of 2119 promoters. In addition to binding core mesenchymal transcription factors SNAI1 and ZEB1, we observed that hexosamine biosynthetic pathway (HBP) enzymes were induced and contained proximal XBP1s peaks. We demonstrate that IRE1α -XBP1s signaling is necessary and sufficient to activate core enzymes by recruiting elongation-competent phospho-Ser2 CTD modified RNA Pol II. We conclude that the IRE1α-XBP1s pathway coordinately regulates mesenchymal transcription factors and hexosamine biosynthesis in EMP by a mechanism involving recruitment of activated pSer2-Pol II to GC-rich promoters.

The exceptionally high positive charge of the histones, concentrated in the N- and C-terminal tails, is believed to contribute to the stability of the nucleosome by neutralizing the negative charge of the nucleosomal DNA. We find, on the contrary, that the high positive charge contributes to instability, performing an essential function in chromatin remodeling. We show that the tails are required for removal of the histone octamer by the RSC chromatin remodeling complex, and this function is not due to direct RSC-tail interaction. We also show that the tails are required for histone octamer transfer from nucleosomes to DNA, and this activity of the tails is a consequence of their positive charge. Thus, the histone tails, intrinsically disordered protein regions, perform a critical role in chromatin structure and transcription, unrelated to their well-known role in regulation through posttranscriptional modification.

In mammals, many germline genes are epigenetically repressed to prevent their illegitimate expression in somatic cells. To advance our understanding of the mechanisms restricting the expression of germline genes, we analyzed their chromatin signature and performed a CRISPR-Cas9 knock-out screen for genes involved in germline gene repression using a Dazl-GFP reporter system in mouse embryonic stem cells (mESCs). We show that the repression of germline genes mainly depends on the polycomb complex PRC1.6 and DNA methylation, which function additively in mESCs. Furthermore, we validated novel genes involved in the repression of germline genes and characterized three of them: Usp7, Shfm1 (also known as Sem1) and Erh. Inactivation of Usp7, Shfm1 or Erh led to the upregulation of germline genes, as well as retrotransposons for Shfm1, in mESCs. Mechanistically, USP7 interacts with PRC1.6 components, promotes PRC1.6 stability and presence at germline genes, and facilitates DNA methylation deposition at germline gene promoters for long term repression. Our study provides a global view of the mechanisms and novel factors required for silencing germline genes in embryonic stem cells.

The Mousepost 1.0 online search tool, launched in 2017, allowed to search for variations in all protein-coding gene sequences of 36 sequenced mouse inbred strains, compared to the reference strain C57BL/6J, which could be linked to strain-specific phenotypes and modifier effects. Because recently these genome sequences have been significantly updated and sequences of 16 extra strains added by the Mouse Genomes Project, a profound update, correction and expansion of the Mousepost 1.0 database has been performed and is reported here. Moreover, we have added a new class of protein disturbing sequence polymorphisms (besides stop codon losses, stop codon gains, small insertions and deletions, and missense mutations), namely start codon mutations. The current version, Mousepost 2.0 (https://mousepost.be), therefore is a significantly updated and invaluable tool available to the community and is described here and foreseen by multiple examples.

Inference of global gene regulatory networks from omics data is a long-term goal of systems biology. Most methods developed for inferring transcription factor (TF)-gene interactions either relied on a small dataset or used snapshot data which is not suitable for inferring a process that is inherently temporal. Here, we developed a new computational method that combines neural networks and multi-task learning to predict RNA velocity rather than gene expression values. This allows our method to overcome many of the problems faced by prior methods leading to more accurate and more comprehensive set of identified regulatory interactions. Application of our method to atlas scale single cell data from 6 HuBMAP tissues led to several validated and novel predictions and greatly improved on prior methods proposed for this task.

The uridine at the 34th position of tRNA, which is able to base pair with the 3'-end codon on mRNA, is usually modified to influence many aspects of decoding properties during translation. Derivatives of 5-methyluridine (xm5U), which include methylaminomethyl (mnm-) or carboxymethylaminomethyl (cmnm-) groups at C5 of uracil base, are widely conserved at the 34th position of many prokaryotic tRNAs. In Gram-negative bacteria such as Escherichia coli, a bifunctional MnmC is involved in the last two reactions of the biosynthesis of mnm5(s2)U, in which the enzyme first converts cmnm5(s2)U to 5-aminomethyl-(2-thio)uridine (nm5(s2)U) and subsequently installs the methyl group to complete the formation of mnm5(s2)U. Although mnm5s2U has been identified in tRNAs of Gram-positive bacteria and plants as well, their genomes do not contain an mnmC ortholog and the gene(s) responsible for this modification is unknown. We discovered that MnmM, previously known as YtqB, is the methyltransferase that converts nm5s2U to mnm5s2U in Bacillus subtilis through comparative genomics, gene complementation experiments, and in vitro assays. Furthermore, we determined X-ray crystal structures of MnmM complexed with anticodon stem loop of tRNAGln. The structures provide the molecular basis underlying the importance of U33-nm5s2U34-U35 as the key determinant for the specificity of MnmM.

In Escherichia coli, replication of both strands of genomic DNA is carried out by a single replicase-DNA polymerase III holoenzyme (pol III HE). However, in certain genetic backgrounds, the low-fidelity TLS polymerase, DNA polymerase V (pol V) gains access to undamaged genomic DNA where it promotes elevated levels of spontaneous mutagenesis preferentially on the lagging strand. We employed active site mutants of pol III (pol IIIα_S759N) and pol V (pol V_Y11A) to analyze ribonucleotide incorporation and removal from the E. coli chromosome on a genome-wide scale under conditions of normal replication, as well as SOS induction. Using a variety of methods tuned to the specific properties of these polymerases (analysis of lacI mutational spectra, lacZ reversion assay, HydEn-seq, alkaline gel electrophoresis), we present evidence that repair of ribonucleotides from both DNA strands in E. coli is unequal. While RNase HII plays a primary role in leading-strand Ribonucleotide Excision Repair (RER), the lagging strand is subject to other repair systems (RNase HI and under conditions of SOS activation also Nucleotide Excision Repair). Importantly, we suggest that RNase HI activity can also influence the repair of single ribonucleotides incorporated by the replicase pol III HE into the lagging strand.

Riboswitches regulate downstream gene expression by binding cellular metabolites. Regulation of translation initiation by riboswitches is posited to occur by metabolite-mediated sequestration of the Shine-Dalgarno sequence (SDS), causing bypass by the ribosome. Recently, we solved a co-crystal structure of a prequeuosine1-sensing riboswitch from Carnobacterium antarcticum that binds two metabolites in a single pocket. The structure revealed that the second nucleotide within the gene-regulatory SDS, G34, engages in a crystal contact, obscuring the molecular basis of gene regulation. Here, we report a co-crystal structure wherein C10 pairs with G34. However, molecular dynamics simulations reveal quick dissolution of the pair, which fails to reform. Functional and chemical probing assays inside live bacterial cells corroborate the dispensability of the C10-G34 pair in gene regulation, leading to the hypothesis that the compact pseudoknot fold is sufficient for translation attenuation. Remarkably, the C. antarcticum aptamer retained significant gene-regulatory activity when uncoupled from the SDS using unstructured spacers up to 10 nucleotides away from the riboswitch-akin to steric-blocking employed by sRNAs. Accordingly, our work reveals that the RNA fold regulates translation without SDS sequestration, expanding known riboswitch-mediated gene-regulatory mechanisms. The results infer that riboswitches exist wherein the SDS is not embedded inside a stable fold.

Transposons are mobile genetic elements prevalent in the genomes of most species. The distribution of transposons within a genome reflects the actions of two opposing processes: initial insertion site selection, and selective pressure from the host. By analyzing whole-genome sequencing data from transposon-activated Drosophila melanogaster, we identified 43 316 de novo and 237 germline insertions from four long-terminal-repeat (LTR) transposons, one LINE transposon (I-element), and one DNA transposon (P-element). We found that all transposon types favored insertion into promoters de novo, but otherwise displayed distinct insertion patterns. De novo and germline P-element insertions preferred replication origins, often landing in a narrow region around transcription start sites and in regions of high chromatin accessibility. De novo LTR transposon insertions preferred regions with high H3K36me3, promoters and exons of active genes; within genes, LTR insertion frequency correlated with gene expression. De novo I-element insertion density increased with distance from the centromere. Germline I-element and LTR transposon insertions were depleted in promoters and exons, suggesting strong selective pressure to remove transposons from functional elements. Transposon movement is associated with genome evolution and disease; therefore, our results can improve our understanding of genome and disease biology.

Epigenetic information defines tissue identity and is largely inherited in development through DNA methylation. While studied mostly for mean differences, methylation also encodes stochastic change, defined as entropy in information theory. Analyzing allele-specific methylation in 49 human tissue sample datasets, we find that methylation entropy is associated with specific DNA binding motifs, regulatory DNA, and CpG density. Then applying information theory to 42 mouse embryo methylation datasets, we find that the contribution of methylation entropy to time- and tissue-specific patterns of development is comparable to the contribution of methylation mean, and methylation entropy is associated with sequence and chromatin features conserved with human. Moreover, methylation entropy is directly related to gene expression variability in development, suggesting a role for epigenetic entropy in developmental plasticity.

Transcriptional pause is essential for all types of termination. In this single-molecule study on bacterial Rho factor-dependent terminators, we confirm that the three Rho-dependent termination routes operate compatibly together in a single terminator, and discover that their termination efficiencies depend on the terminational pauses in unexpected ways. Evidently, the most abundant route is that Rho binds nascent RNA first and catches up with paused RNA polymerase (RNAP) and this catch-up Rho mediates simultaneous releases of transcript RNA and template DNA from RNAP. The fastest route is that the catch-up Rho effects RNA-only release and leads to 1D recycling of RNAP on DNA. The slowest route is that the RNAP-prebound stand-by Rho facilitates only the simultaneous rather than sequential releases. Among the three routes, only the stand-by Rho's termination efficiency positively correlates with pause duration, contrary to a long-standing speculation, invariably in the absence or presence of NusA/NusG factors, competitor RNAs or a crowding agent. Accordingly, the essential terminational pause does not need to be long for the catch-up Rho's terminations, and long pauses benefit only the stand-by Rho's terminations. Furthermore, the Rho-dependent termination of mgtA and ribB riboswitches is controlled mainly by modulation of the stand-by rather than catch-up termination.

DNA double-strand breaks (DSBs) can be repaired by either homologous recombination (HR) or non-homologous end-joining (NHEJ). NHEJ is induced by the binding to DSBs of the Ku70-Ku80 heterodimer, which acts as a hub for the recruitment of downstream NHEJ components. An important issue in DSB repair is the maintenance of the DSB ends in close proximity, a function that in yeast involves the MRX complex and Sae2. Here, we provide evidence that Ku contributes to keep the DNA ends tethered to each other. The ku70-C85Y mutation, which increases Ku affinity for DNA and its persistence very close to the DSB ends, enhances DSB end-tethering and suppresses the end-tethering defect of sae2Δ cells. Impairing histone removal around DSBs either by eliminating Tel1 kinase activity or nucleosome remodelers enhances Ku persistence at DSBs and DSB bridging, suggesting that Tel1 antagonizes the Ku function in supporting end-tethering by promoting nucleosome removal and possibly Ku sliding inwards. As Ku provides a block to DSB resection, this Tel1 function can be important to regulate the mode by which DSBs are repaired.

Integrative and conjugative elements (ICEs) are mobile genetic elements that can transfer by conjugation to recipient cells. Some ICEs integrate into a unique site in the genome of their hosts. We studied quantitatively the process by which an ICE searches for its unique integration site in the Bacillus subtilis chromosome. We followed the motion of both ICEBs1 and the chromosomal integration site in real time within individual cells. ICEBs1 exhibited a wide spectrum of dynamical behaviors, ranging from rapid sub-diffusive displacements crisscrossing the cell, to kinetically trapped states. The chromosomal integration site moved sub-diffusively and exhibited pronounced dynamical asymmetry between longitudinal and transversal motions, highlighting the role of chromosomal structure and the heterogeneity of the bacterial interior in the search. The successful search for and subsequent recombination into the integration site is a key step in the acquisition of integrating mobile genetic elements. Our findings provide new insights into intracellular transport processes involving large DNA molecules.