Plasmids can be transferred between cells by conjugation, thereby driving bacterial evolution by horizontal gene transfer. Yet, we ignore the molecular mechanisms of transfer for many plasmids because they lack all protein-coding genes required for conjugation. We solved this conundrum by identifying hundreds of plasmids and chromosomes with conjugative origins of transfer in Escherichia coli and Staphylococcus aureus. These plasmids (pOriT) hijack the relaxases of conjugative or mobilizable elements, but not both. The functional dependencies between pOriT and other plasmids explain their co-occurrence: pOriT are abundant in cells with many plasmids, whereas conjugative plasmids are the most common in the others. We systematically characterized plasmid mobility in relation to conjugation and alternative mechanisms of transfer and can now propose a putative mechanism of transfer for ∼90% of them. In most cases, plasmid mobility seems to involve conjugation. Interestingly, the mechanisms of mobility are important determinants of plasmid-encoded accessory traits, since pOriTs have the highest densities of antimicrobial resistance genes, whereas plasmids lacking putative mechanisms of transfer have the lowest. We illuminate the evolutionary relationships between plasmids and suggest that many pOriT may have arisen by gene deletions in other types of plasmids. These results suggest that most plasmids can be transferred by conjugation.

The reversible attachment of ubiquitin governs the interaction, activity and degradation of proteins whereby the type and target of this conjugation determine the biological response. The investigation of this complex and multi-faceted protein ubiquitination mostly relies on painstaking biochemical analyses. Here, we employ recombinant binding domains to probe the ubiquitination of proteins in living cells. We immobilize GFP-fused proteins of interest at a distinct cellular structure and detect their ubiquitination state with red fluorescent ubiquitin binders. With this ubiquitin fluorescent three-hybrid (ubiF3H) assay we identified HP1β as a novel ubiquitination target of UHRF1. The use of linkage specific ubiquitin binding domains enabled the discrimination of K48 and K63 linked protein ubiquitination. To enhance signal-to-noise ratio, we implemented fluorescence complementation (ubiF3Hc) with split YFP. Using in addition a cell cycle marker we could show that HP1β is mostly ubiquitinated by UHRF1 during S phase and deubiquitinated by the protease USP7. With this complementation assay we could also directly detect the ubiquitination of the tumor suppressor p53 and monitor its inhibition by the anti-cancer drug Nutlin-3. Altogether, we demonstrate the utility of the ubiF3H assay to probe the ubiquitination of specific proteins and to screen for ligases, proteases and small molecules controlling this posttranslational modification.

The use of synthetic biological circuits to deal with numerous biological challenges has been proposed in several studies, but its implementation is still remote. A major problem encountered is the complexity of the cellular engineering needed to achieve complex biological circuits and the lack of general-purpose biological systems. The generation of re-programmable circuits can increase circuit flexibility and the scalability of complex cell-based computing devices. Here we present a new architecture to produce reprogrammable biological circuits that allow the development of a variety of different functions with minimal cell engineering. We demonstrate the feasibility of creating several circuits using only a small set of engineered cells, which can be externally reprogrammed to implement simple logics in response to specific inputs. In this regard, depending on the computation needs, a device composed of a number of defined cells can generate a variety of circuits without the need of further cell engineering or rearrangements. In addition, the inclusion of a memory module in the circuits strongly improved the digital response of the devices. The reprogrammability of biological circuits is an intrinsic capacity that is not provided in electronics and it may be used as a tool to solve complex biological problems.

RNA polymerase II (RNAPII) transcribes small nuclear RNA (snRNA) genes in close proximity to Cajal bodies, subnuclear compartments that depend on the SUMO isopeptidase USPL1 for their assembly. We show here that overexpression of USPL1 as well as of another nuclear SUMO isopeptidase, SENP6, alters snRNA 3'-end cleavage, a process carried out by the Integrator complex. Beyond its role in snRNA biogenesis, this complex is responsible for regulating the expression of different RNAPII transcripts. While several subunits of the complex are SUMO conjugation substrates, we found that the SUMOylation of the INTS11 subunit is regulated by USPL1 and SENP6. We defined Lys381, Lys462 and Lys475 as bona fide SUMO attachment sites on INTS11 and observed that SUMOylation of this protein modulates its subcellular localization and is required for Integrator activity. Moreover, while an INTS11 SUMOylation-deficient mutant is still capable of interacting with INTS4 and INTS9, its interaction with other subunits of the complex is affected. These findings point to a regulatory role for SUMO conjugation on Integrator activity and suggest the involvement of INTS11 SUMOylation in the assembly of the complex. Furthermore, this work adds Integrator-dependent RNA processing to the growing list of cellular processes regulated by SUMO conjugation.

The genome consists of non-B-DNA structures such as G-quadruplexes (G4) that are involved in the regulation of genome stability and transcription. Telomeric-repeat containing RNA (TERRA) is capable of folding into G-quadruplex and interacting with chromatin remodeler ATRX. Here we show that TERRA modulates ATRX occupancy on repetitive sequences and over genes, and maintains DNA G-quadruplex structures at TERRA target and non-target sites in mouse embryonic stem cells. TERRA prevents ATRX from binding to subtelomeric regions and represses H3K9me3 formation. G4 ChIP-seq reveals that G4 abundance decreases at accessible chromatin regions, particularly at transcription start sites (TSS) after TERRA depletion; such G4 reduction at TSS is associated with elevated ATRX occupancy and differentially expressed genes. Loss of ATRX alleviates the effect of gene repression caused by TERRA depletion. Immunostaining analyses demonstrate that knockdown of TERRA diminishes DNA G4 signals, whereas silencing ATRX elevates G4 formation. Our results uncover an epigenetic regulation by TERRA that sequesters ATRX and preserves DNA G4 structures.

Despite being a member of the chromodomain helicase DNA-binding protein family, little is known about the exact role of CHD6 in chromatin remodeling or cancer disease. Here we show that CHD6 binds to chromatin to promote broad nucleosome eviction for transcriptional activation of many cancer pathways. By integrating multiple patient cohorts for bioinformatics analysis of over a thousand prostate cancer datasets, we found CHD6 expression elevated in prostate cancer and associated with poor prognosis. Further comprehensive experiments demonstrated that CHD6 regulates oncogenicity of prostate cancer cells and tumor development in a murine xenograft model. ChIP-Seq for CHD6, along with MNase-Seq and RNA-Seq, revealed that CHD6 binds on chromatin to evict nucleosomes from promoters and gene bodies for transcriptional activation of oncogenic pathways. These results demonstrated a key function of CHD6 in evicting nucleosomes from chromatin for transcriptional activation of prostate cancer pathways.

Human gene research generates new biology insights with translational potential, yet few studies have considered the health of the human gene literature. The accessibility of human genes for targeted research, combined with unreasonable publication pressures and recent developments in scholarly publishing, may have created a market for low-quality or fraudulent human gene research articles, including articles produced by contract cheating organizations known as paper mills. This review summarises the evidence that paper mills contribute to the human gene research literature at scale and outlines why targeted gene research may be particularly vulnerable to systematic research fraud. To raise awareness of targeted gene research from paper mills, we highlight features of problematic manuscripts and publications that can be detected by gene researchers and/or journal staff. As improved awareness and detection could drive the further evolution of paper mill-supported publications, we also propose changes to academic publishing to more effectively deter and correct problematic publications at scale. In summary, the threat of paper mill-supported gene research highlights the need for all researchers to approach the literature with a more critical mindset, and demand publications that are underpinned by plausible research justifications, rigorous experiments and fully transparent reporting.

Several basic leucine zipper (bZIP) transcription factors have accessory motifs in their DNA-binding domains, such as the CNC motif of CNC family or the EHR motif of small Maf (sMaf) proteins. CNC family proteins heterodimerize with sMaf proteins to recognize CNC-sMaf binding DNA elements (CsMBEs) in competition with sMaf homodimers, but the functional role of the CNC motif remains elusive. In this study, we report the crystal structures of Nrf2/NFE2L2, a CNC family protein regulating anti-stress transcriptional responses, in a complex with MafG and CsMBE. The CNC motif restricts the conformations of crucial Arg residues in the basic region, which form extensive contact with the DNA backbone phosphates. Accordingly, the Nrf2-MafG heterodimer has approximately a 200-fold stronger affinity for CsMBE than canonical bZIP proteins, such as AP-1 proteins. The high DNA affinity of the CNC-sMaf heterodimer may allow it to compete with the sMaf homodimer on target genes without being perturbed by other low-affinity bZIP proteins with similar sequence specificity.

R-loops are three-stranded nucleic acid structures composed of an RNA:DNA hybrid and displaced DNA strand. These structures can halt DNA replication when formed co-transcriptionally in the opposite orientation to replication fork progression. A recent study has shown that replication forks stalled by co-transcriptional R-loops can be restarted by a mechanism involving fork cleavage by MUS81 endonuclease, followed by ELL-dependent reactivation of transcription, and fork religation by the DNA ligase IV (LIG4)/XRCC4 complex. However, how R-loops are eliminated to allow the sequential restart of transcription and replication in this pathway remains elusive. Here, we identified the human DDX17 helicase as a factor that associates with R-loops and counteracts R-loop-mediated replication stress to preserve genome stability. We show that DDX17 unwinds R-loops in vitro and promotes MUS81-dependent restart of R-loop-stalled forks in human cells in a manner dependent on its helicase activity. Loss of DDX17 helicase induces accumulation of R-loops and the formation of R-loop-dependent anaphase bridges and micronuclei. These findings establish DDX17 as a component of the MUS81-LIG4-ELL pathway for resolution of R-loop-mediated transcription-replication conflicts, which may be involved in R-loop unwinding.

Ubiquitin-like with PHD and RING finger domain-containing protein 1 (UHRF1)-dependent DNA methylation is essential for maintaining cell fate during cell proliferation. Developmental pluripotency-associated 3 (DPPA3) is an intrinsically disordered protein that specifically interacts with UHRF1 and promotes passive DNA demethylation by inhibiting UHRF1 chromatin localization. However, the molecular basis of how DPPA3 interacts with and inhibits UHRF1 remains unclear. We aimed to determine the structure of the mouse UHRF1 plant homeodomain (PHD) complexed with DPPA3 using nuclear magnetic resonance. Induced α-helices in DPPA3 upon binding of UHRF1 PHD contribute to stable complex formation with multifaceted interactions, unlike canonical ligand proteins of the PHD domain. Mutations in the binding interface and unfolding of the DPPA3 helical structure inhibited binding to UHRF1 and its chromatin localization. Our results provide structural insights into the mechanism and specificity underlying the inhibition of UHRF1 by DPPA3.

Neural differentiation of embryonic stem cells (ESCs) requires precisely orchestrated gene regulation, a process governed in part by changes in 3D chromatin structure. How these changes regulate gene expression in this context remains unclear. In this study, we observed enrichment of the transcription factor KLF4 at some poised or closed enhancers at TSS-linked regions of genes associated with neural differentiation. Combination analysis of ChIP, HiChIP and RNA-seq data indicated that KLF4 loss in ESCs induced changes in 3D chromatin structure, including increased chromatin interaction loops between neural differentiation-associated genes and active enhancers, leading to upregulated expression of neural differentiation-associated genes and therefore early neural differentiation. This study suggests KLF4 inhibits early neural differentiation by regulation of 3D chromatin structure, which is a new mechanism of early neural differentiation.

The hybridization kinetic of an oligonucleotide to its template is a fundamental step in many biological processes such as replication arrest, CRISPR recognition, DNA sequencing, DNA origami, etc. Although single kinetic descriptions exist for special cases of this problem, there are no simple general prediction schemes. In this work, we have measured experimentally, with no fluorescent labelling, the displacement of an oligonucleotide from its substrate in two situations: one corresponding to oligonucleotide binding/unbinding on ssDNA and one in which the oligonucleotide is displaced by the refolding of a dsDNA fork. In this second situation, the fork is expelling the oligonucleotide thus significantly reducing its residence time. To account for our data in these two situations, we have constructed a mathematical model, based on the known nearest neighbour dinucleotide free energies, and provided a good estimate of the residence times of different oligonucleotides (DNA, RNA, LNA) of various lengths in different experimental conditions (force, temperature, buffer conditions, presence of mismatches, etc.). This study provides a foundation for the dynamics of oligonucleotide displacement, a process of importance in numerous biological and bioengineering contexts.

In mammalian cells, chromosomal replication starts at thousands of origins at which replisomes are assembled. Replicative stress triggers additional initiation events from 'dormant' origins whose genomic distribution and regulation are not well understood. In this study, we have analyzed origin activity in mouse embryonic stem cells in the absence or presence of mild replicative stress induced by aphidicolin, a DNA polymerase inhibitor, or by deregulation of origin licensing factor CDC6. In both cases, we observe that the majority of stress-responsive origins are also active in a small fraction of the cell population in a normal S phase, and stress increases their frequency of activation. In a search for the molecular determinants of origin efficiency, we compared the genetic and epigenetic features of origins displaying different levels of activation, and integrated their genomic positions in three-dimensional chromatin interaction networks derived from high-depth Hi-C and promoter-capture Hi-C data. We report that origin efficiency is directly proportional to the proximity to transcriptional start sites and to the number of contacts established between origin-containing chromatin fragments, supporting the organization of origins in higher-level DNA replication factories.

Nucleoid-associated proteins (NAPs) play a central role in chromosome organization and environment-responsive transcription regulation. The Bacillus subtilis-encoded NAP Rok binds preferentially AT-rich regions of the genome, which often contain genes of foreign origin that are silenced by Rok binding. Additionally, Rok plays a role in chromosome architecture by binding in genomic clusters and promoting chromosomal loop formation. Based on this, Rok was proposed to be a functional homolog of E. coli H-NS. However, it is largely unclear how Rok binds DNA, how it represses transcription and whether Rok mediates environment-responsive gene regulation. Here, we investigated Rok's DNA binding properties and the effects of physico-chemical conditions thereon. We demonstrate that Rok is a DNA bridging protein similar to prototypical H-NS-like proteins. However, unlike these proteins, the DNA bridging ability of Rok is not affected by changes in physico-chemical conditions. The DNA binding properties of the Rok interaction partner sRok are affected by salt concentration. This suggests that in a minority of Bacillus strains Rok activity can be modulated by sRok, and thus respond indirectly to environmental stimuli. Despite several functional similarities, the absence of a direct response to physico-chemical changes establishes Rok as disparate member of the H-NS family.

Human pre-mRNA processing relies on multi-subunit macromolecular complexes, which recognize specific RNA sequence elements essential for assembly and activity. Canonical pre-mRNA processing proceeds via the recognition of a polyadenylation signal (PAS) and a downstream sequence element (DSE), and produces polyadenylated mature mRNAs, while replication-dependent (RD) histone pre-mRNA processing requires association with a stem-loop (SL) motif and a histone downstream element (HDE), and produces cleaved but non-polyadenylated mature mRNAs. H2AC18 mRNA, a specific H2A RD histone pre-mRNA, can be processed to give either a non-polyadenylated mRNA, ending at the histone SL, or a polyadenylated mRNA. Here, we reveal how H2AC18 captures the two human pre-mRNA processing complexes in a mutually exclusive mode by overlapping a canonical PAS (AAUAAA) sequence element with a HDE. Disruption of the PAS sequence on H2AC18 pre-mRNA prevents recruitment of the canonical complex in vitro, without affecting the histone machinery. This shows how the relative position of cis-acting elements in histone pre-mRNAs allows the selective recruitment of distinct human pre-mRNA complexes, thereby expanding the capability to regulate 3' processing and polyadenylation.

5-Methyl-cytosine (5mC) is one of the most important DNA modifications and plays versatile biological roles. It is well known that 5mC stabilizes DNA duplexes. However, it remains unclear how 5mC affects the kinetics of DNA melting and hybridization. Here, we studied the kinetics of unzipping and rezipping using a 502-bp DNA hairpin by single-molecule magnetic tweezers. Under constant loading rates, 5mC increases the unzipping force but counterintuitively decreases the rezipping force at various salt and temperature conditions. Under constant forces, the non-methylated DNA hops between metastable states during unzipping and rezipping, which implies low energy barriers. Surprisingly, the 5mC DNA can't rezip after fully unzipping unless much lower forces are applied, where it rezips stochastically in a one-step manner, which implies 5mC kinetically hinders DNA hybridization and high energy barriers in DNA hybridization. All-atom molecular dynamics simulations reveal that the 5mC kinetically hinders DNA hybridization due to steric effects rather than electrostatic effects caused by the additional methyl groups of cytosines. Considering the possible high speed of DNA unzipping and zipping during replication and transcription, our findings provide new insights into the biological roles of 5mC.

Bacterial RNases process RNAs until only short oligomers (2-5 nucleotides) remain, which are then processed by one or more specialized enzymes until only nucleoside monophosphates remain. Oligoribonuclease (Orn) is an essential enzyme that acts in this capacity. However, many bacteria do not encode for Orn and instead encode for NanoRNase A (NrnA). Yet, the catalytic mechanism, cellular roles and physiologically relevant substrates have not been fully resolved for NrnA proteins. We herein utilized a common set of reaction assays to directly compare substrate preferences exhibited by NrnA-like proteins from Bacillus subtilis, Enterococcus faecalis, Streptococcus pyogenes and Mycobacterium tuberculosis. While the M. tuberculosis protein specifically cleaved cyclic di-adenosine monophosphate, the B. subtilis, E. faecalis and S. pyogenes NrnA-like proteins uniformly exhibited striking preference for short RNAs between 2-4 nucleotides in length, all of which were processed from their 5' terminus. Correspondingly, deletion of B. subtilis nrnA led to accumulation of RNAs between 2 and 4 nucleotides in length in cellular extracts. Together, these data suggest that many Firmicutes NrnA-like proteins are likely to resemble B. subtilis NrnA to act as a housekeeping enzyme for processing of RNAs between 2 and 4 nucleotides in length.

DNA polymerase α (Polα) is essential for DNA replication initiation and makes a notable contribution to genome mutagenesis. The activity and fidelity of Polα during the early steps of DNA replication have not been well studied. Here we show that at the beginning of DNA synthesis, when extending the RNA primer received from primase, Polα is more mutagenic than during the later DNA elongation steps. Kinetic and binding studies revealed substantially higher activity and affinity to the template:primer when Polα interacts with ribonucleotides of a chimeric RNA-DNA primer. Polα activity greatly varies during first six steps of DNA synthesis, and the bias in the rates of correct and incorrect dNTP incorporation leads to impaired fidelity, especially upon the second step of RNA primer extension. Furthermore, increased activity and stability of Polα/template:primer complexes containing RNA-DNA primers result in higher efficiency of mismatch extension.

RNA helicases are involved in multiple steps of RNA metabolism to direct their roles in gene expression, yet their functions in pluripotency control remain largely unexplored. Starting from an RNA interference (RNAi) screen of RNA helicases, we identified that eIF4A3, a DEAD-box (Ddx) helicase component of the exon junction complex (EJC), is essential for the maintenance of embryonic stem cells (ESCs). Mechanistically, we show that eIF4A3 post-transcriptionally controls the pluripotency-related cell cycle regulators and that its depletion causes the loss of pluripotency via cell cycle dysregulation. Specifically, eIF4A3 is required for the efficient nuclear export of Ccnb1 mRNA, which encodes Cyclin B1, a key component of the pluripotency-promoting pathway during the cell cycle progression of ESCs. Our results reveal a previously unappreciated role for eIF4A3 and its associated EJC in maintaining stem cell pluripotency through post-transcriptional control of the cell cycle.

Most cancer alterations occur in the noncoding portion of the human genome, where regulatory regions control gene expression. The discovery of noncoding mutations altering the cells' regulatory programs has been limited to few examples with high recurrence or high functional impact. Here, we show that transcription factor binding sites (TFBSs) have similar mutation loads to those in protein-coding exons. By combining cancer somatic mutations in TFBSs and expression data for protein-coding and miRNA genes, we evaluate the combined effects of transcriptional and post-transcriptional alterations on the regulatory programs in cancers. The analysis of seven TCGA cohorts culminates with the identification of protein-coding and miRNA genes linked to mutations at TFBSs that are associated with a cascading trans-effect deregulation on the cells' regulatory programs. Our analyses of cis-regulatory mutations associated with miRNAs recurrently predict 12 mature miRNAs (derived from 7 precursors) associated with the deregulation of their target gene networks. The predictions are enriched for cancer-associated protein-coding and miRNA genes and highlight cis-regulatory mutations associated with the dysregulation of key pathways associated with carcinogenesis. By combining transcriptional and post-transcriptional regulation of gene expression, our method predicts cis-regulatory mutations related to the dysregulation of key gene regulatory networks in cancer patients.

Animal cell lines often undergo extreme genome restructuring events, including polyploidy and segmental aneuploidy that can impede de novo whole-genome assembly (WGA). In some species like Drosophila, cell lines also exhibit massive proliferation of transposable elements (TEs). To better understand the role of transposition during animal cell culture, we sequenced the genome of the tetraploid Drosophila S2R+ cell line using long-read and linked-read technologies. WGAs for S2R+ were highly fragmented and generated variable estimates of TE content across sequencing and assembly technologies. We therefore developed a novel WGA-independent bioinformatics method called TELR that identifies, locally assembles, and estimates allele frequency of TEs from long-read sequence data (https://github.com/bergmanlab/telr). Application of TELR to a ∼130x PacBio dataset for S2R+ revealed many haplotype-specific TE insertions that arose by transposition after initial cell line establishment and subsequent tetraploidization. Local assemblies from TELR also allowed phylogenetic analysis of paralogous TEs, which revealed that proliferation of TE families in vitro can be driven by single or multiple source lineages. Our work provides a model for the analysis of TEs in complex heterozygous or polyploid genomes that are recalcitrant to WGA and yields new insights into the mechanisms of genome evolution in animal cell culture.

Although long-read sequencing can often enable chromosome-level reconstruction of genomes, it is still unclear how one can routinely obtain gapless assemblies. In the model plant Arabidopsis thaliana, other than the reference accession Col-0, all other accessions de novo assembled with long-reads until now have used PacBio continuous long reads (CLR). Although these assemblies sometimes achieved chromosome-arm level contigs, they inevitably broke near the centromeres, excluding megabases of DNA from analysis in pan-genome projects. Since PacBio high-fidelity (HiFi) reads circumvent the high error rate of CLR technologies, albeit at the expense of read length, we compared a CLR assembly of accession Eyach15-2 to HiFi assemblies of the same sample. The use of five different assemblers starting from subsampled data allowed us to evaluate the impact of coverage and read length. We found that centromeres and rDNA clusters are responsible for 71% of contig breaks in the CLR scaffolds, while relatively short stretches of GA/TC repeats are at the core of >85% of the unfilled gaps in our best HiFi assemblies. Since the HiFi technology consistently enabled us to reconstruct gapless centromeres and 5S rDNA clusters, we demonstrate the value of the approach by comparing these previously inaccessible regions of the genome between the Eyach15-2 accession and the reference accession Col-0.

RNA is modified by hundreds of chemical reactions and folds into innumerable shapes. However, the regulatory role of RNA sequence and structure and how dysregulation leads to diseases remain largely unknown. Here, we uncovered a mechanism where RNA abasic sites in R-loops regulate transcription by pausing RNA polymerase II. We found an enhancer RNA, AANCR, that regulates the transcription and expression of apolipoprotein E (APOE). In some human cells such as fibroblasts, AANCR is folded into an R-loop and modified by N-glycosidic cleavage; in this form, AANCR is a partially transcribed nonfunctional enhancer and APOE is not expressed. In contrast, in other cell types including hepatocytes and under stress, AANCR does not form a stable R-loop as its sequence is not modified, so it is transcribed into a full-length enhancer that promotes APOE expression. DNA sequence variants in AANCR are associated significantly with APOE expression and Alzheimer's Disease, thus AANCR is a modifier of Alzheimer's Disease. Besides AANCR, thousands of noncoding RNAs are regulated by abasic sites in R-loops. Together our data reveal the essentiality of the folding and modification of RNA in cellular regulation and demonstrate that dysregulation underlies common complex diseases such as Alzheimer's disease.

Although single-cell sequencing has provided a powerful tool to deconvolute cellular heterogeneity of diseases like cancer, extrapolating clinical significance or identifying clinically-relevant cells remains challenging. Here, we propose a novel computational method scAB, which integrates single-cell genomics data with clinically annotated bulk sequencing data via a knowledge- and graph-guided matrix factorization model. Once combined, scAB provides a coarse- and fine-grain multiresolution perspective of phenotype-associated cell states and prognostic signatures previously not visible by single-cell genomics. We use scAB to enhance live cancer single-cell RNA-seq data, identifying clinically-relevant previously unrecognized cancer and stromal cell subsets whose signatures show a stronger poor-survival association. The identified fine-grain cell subsets are associated with distinct cancer hallmarks and prognosis power. Furthermore, scAB demonstrates its utility as a biomarker identification tool, with the ability to predict immunotherapy, drug responses and survival when applied to melanoma single-cell RNA-seq datasets and glioma single-cell ATAC-seq datasets. Across multiple single-cell and bulk datasets from different cancer types, we also demonstrate the superior performance of scAB in generating prognosis signatures and survival predictions over existing models. Overall, scAB provides an efficient tool for prioritizing clinically-relevant cell subsets and predictive signatures, utilizing large publicly available databases to improve prognosis and treatments.

The PglZ family of proteins belongs to the alkaline phosphatase superfamily, which consists of metallohydrolases with limited sequence identity but similar metal-coordination architectures in otherwise divergent active sites. Proteins with a well-defined PglZ domain are ubiquitous among prokaryotes as essential components of BREX phage defence systems and two-component systems (TCSs). Whereas other members of the alkaline phosphatase superfamily are well characterized, the activity, structure and biological function of PglZ family proteins remain unclear. We therefore investigated the structure and function of PorX, an orphan response regulator of the Porphyromonas gingivalis TCS containing a putative PglZ effector domain. The crystal structure of PorX revealed a canonical receiver domain, a helical bundle, and an unprecedented PglZ domain, similar to the general organization of the phylogenetically related BREX-PglZ proteins. The PglZ domain of PorX features an active site cleft suitable for large substrates. An extensive search for substrates revealed that PorX is a phosphodiesterase that acts on cyclic and linear oligonucleotides, including signalling molecules such as cyclic oligoadenylates. These results, combined with mutagenesis, biophysical and enzymatic analysis, suggest that PorX coordinates oligonucleotide signalling pathways and indirectly regulates gene expression to control the secretion of virulence factors.

In Escherichia coli, the heat shock protein 15 (Hsp15) is part of the cellular response to elevated temperature. Hsp15 interacts with peptidyl-tRNA-50S complexes that arise upon dissociation of translating 70S ribosomes, and is proposed to facilitate their rescue and recycling. A previous structure of E. coli Hsp15 in complex with peptidyl-tRNA-50S complex reported a binding site located at the central protuberance of the 50S subunit. By contrast, recent structures of RqcP, the Hsp15 homolog in Bacillus subtilis, in complex with peptidyl-tRNA-50S complexes have revealed a distinct site positioned between the anticodon-stem-loop (ASL) of the P-site tRNA and H69 of the 23S rRNA. Here we demonstrate that exposure of E. coli cells to heat shock leads to a decrease in 70S ribosomes and accumulation of 50S subunits, thus identifying a natural substrate for Hsp15 binding. Additionally, we have determined a cryo-EM reconstruction of the Hsp15-50S-peptidyl-tRNA complex isolated from heat shocked E. coli cells, revealing that Hsp15 binds to the 50S-peptidyl-tRNA complex analogously to its B. subtilis homolog RqcP. Collectively, our findings support a model where Hsp15 stabilizes the peptidyl-tRNA in the P-site and thereby promotes access to the A-site for putative rescue factors to release the aberrant nascent polypeptide chain.

Bacterial chromosome topology is controlled by topoisomerases and nucleoid-associated proteins (NAPs). While topoisomerases regulate DNA supercoiling, NAPs introduce bends or coat DNA upon its binding, affecting DNA loop formation. Streptomyces, hyphal, multigenomic bacteria known for producing numerous clinically important compounds, use the highly processive topoisomerase I (TopA) to remove excessive negative DNA supercoils. Elongated vegetative Streptomyces cells contain multiple copies of their linear chromosome, which remain relaxed and relatively evenly distributed. Here, we explored how TopA cooperates with HupA, an HU homologue that is the most abundant Streptomyces NAP. We verified that HupA has an increased affinity for supercoiled DNA in vivo and in vitro. Analysis of mutant strains demonstrated that HupA elimination is detrimental under high DNA supercoiling conditions. The absence of HupA, combined with decreased TopA levels, disrupted chromosome distribution in hyphal cells, eventually inhibiting hyphal growth. We concluded that increased HupA binding to DNA under elevated chromosome supercoiling conditions is critical for the preservation of chromosome organisation.

The action of Type II restriction-modification (RM) systems depends on restriction endonuclease (REase), which cleaves foreign DNA at specific sites, and methyltransferase (MTase), which protects host genome from restriction by methylating the same sites. We here show that protection from phage infection increases as the copy number of plasmids carrying the Type II RM Esp1396I system is increased. However, since increased plasmid copy number leads to both increased absolute intracellular RM enzyme levels and to a decreased MTase/REase ratio, it is impossible to determine which factor determines resistance/susceptibility to infection. By controlled expression of individual Esp1396I MTase or REase genes in cells carrying the Esp1396I system, we show that a shift in the MTase to REase ratio caused by overproduction of MTase or REase leads, respectively, to decreased or increased protection from infection. Consistently, due to stochastic variation of MTase and REase amount in individual cells, bacterial cells that are productively infected by bacteriophage have significantly higher MTase to REase ratios than cells that ward off the infection. Our results suggest that cells with transiently increased MTase to REase ratio at the time of infection serve as entry points for unmodified phage DNA into protected bacterial populations.

The nucleic acid polymerase-catalyzed nucleotidyl transfer reaction associated with polymerase active site closure is a key step in the nucleotide addition cycle (NAC). Two proton transfer events can occur in such a nucleotidyl transfer: deprotonation of the priming nucleotide 3'-hydroxyl nucleophile and protonation of the pyrophosphate (PPi) leaving group. In viral RNA-dependent RNA polymerases (RdRPs), whether and how active site residues participate in this two-proton transfer reaction remained to be clarified. Here we report a 2.5 Å resolution crystal structure of enterovirus 71 (EV71) RdRP in a catalytically closed pre-chemistry conformation, with a proposed proton donor candidate K360 in close contact with the NTP γ-phosphate. Enzymology data reveal that K360 mutations not only reduce RdRP catalytic efficiency but also alter pH dependency profiles in both elongation and pre-elongation synthesis modes. Interestingly, mutations at R174, an RdRP-invariant residue in motif F, had similar effects with additional impact on the Michaelis constant of NTP (KM,NTP). However, direct participation in protonation was not evident for K360 or R174. Our data suggest that both K360 and R174 participate in nucleotidyl transfer, while their possible roles in acid-base or positional catalysis are discussed in comparison with other classes of nucleic acid polymerases.

The human APOBEC family of eleven cytosine deaminases use RNA and single-stranded DNA (ssDNA) as substrates to deaminate cytosine to uracil. This deamination event has roles in lipid metabolism by altering mRNA coding, adaptive immunity by causing evolution of antibody genes, and innate immunity through inactivation of viral genomes. These benefits come at a cost where some family members, primarily from the APOBEC3 subfamily (APOBEC3A-H, excluding E), can cause off-target deaminations of cytosine to form uracil on transiently single-stranded genomic DNA, which induces mutations that are associated with cancer evolution. Since uracil is only promutagenic, the mutations observed in cancer genomes originate only when uracil is not removed by uracil DNA glycosylase (UNG) or when the UNG-induced abasic site is erroneously repaired. However, when ssDNA is present, replication protein A (RPA) binds and protects the DNA from nucleases or recruits DNA repair proteins, such as UNG. Thus, APOBEC enzymes must compete with RPA to access their substrate. Certain APOBEC enzymes can displace RPA, bind and scan ssDNA efficiently to search for cytosines, and can become highly overexpressed in tumor cells. Depending on the DNA replication conditions and DNA structure, RPA can either be in excess or deficient. Here we discuss the interplay between these factors and how despite RPA, multiple cancer genomes have a mutation bias at cytosines indicative of APOBEC activity.

Recognition of single-stranded RNA (ssRNA) by RNA recognition motif (RRM) domains is an important class of protein-RNA interactions. Many such complexes were characterized using nuclear magnetic resonance (NMR) and/or X-ray crystallography techniques, revealing ensemble-averaged pictures of the bound states. However, it is becoming widely accepted that better understanding of protein-RNA interactions would be obtained from ensemble descriptions. Indeed, earlier molecular dynamics simulations of bound states indicated visible dynamics at the RNA-RRM interfaces. Here, we report the first atomistic simulation study of spontaneous binding of short RNA sequences to RRM domains of HuR and SRSF1 proteins. Using a millisecond-scale aggregate ensemble of unbiased simulations, we were able to observe a few dozen binding events. HuR RRM3 utilizes a pre-binding state to navigate the RNA sequence to its partially disordered bound state and then to dynamically scan its different binding registers. SRSF1 RRM2 binding is more straightforward but still multiple-pathway. The present study necessitated development of a goal-specific force field modification, scaling down the intramolecular van der Waals interactions of the RNA which also improves description of the RNA-RRM bound state. Our study opens up a new avenue for large-scale atomistic investigations of binding landscapes of protein-RNA complexes, and future perspectives of such research are discussed.

The formation of triplex DNA is a site-specific recognition method that directly targets duplex DNA. However, triplex DNA formation is generally formed for the GC and AT base pairs of duplex DNA, and there are no natural nucleotides that recognize the CG and TA base pairs, or even the 5-methyl-CG (5mCG) base pair. Moreover, duplex DNA, including 5mCG base pairs, epigenetically regulates gene expression in vivo, and thus targeting strategies are of biological importance. Therefore, the development of triplex-forming oligonucleotides (TFOs) with artificial nucleosides that selectively recognize these base pairs with high affinity is needed. We recently reported that 2'-deoxy-2-aminonebularine derivatives exhibited the ability to recognize 5mCG and CG base pairs in triplex formation; however, this ability was dependent on sequences. Therefore, we designed and synthesized new nucleoside derivatives based on the 2'-deoxy-nebularine (dN) skeleton to shorten the linker length connecting to the hydrogen-bonding unit in formation of the antiparallel motif triplex. We successfully demonstrated that TFOs with 2-guanidinoethyl-2'-deoxynebularine (guanidino-dN) recognized 5mCG and CG base pairs with very high affinity in all four DNA sequences with different adjacent nucleobases of guanidino-dN as well as in the promoter sequences of human genes containing 5mCG base pairs with a high DNA methylation frequency.

RNA folds cotranscriptionally to traverse out-of-equilibrium intermediate structures that are important for RNA function in the context of gene regulation. To investigate this process, here we study the structure and function of the Bacillus subtilis yxjA purine riboswitch, a transcriptional riboswitch that downregulates a nucleoside transporter in response to binding guanine. Although the aptamer and expression platform domain sequences of the yxjA riboswitch do not completely overlap, we hypothesized that a strand exchange process triggers its structural switching in response to ligand binding. In vivo fluorescence assays, structural chemical probing data and experimentally informed secondary structure modeling suggest the presence of a nascent intermediate central helix. The formation of this central helix in the absence of ligand appears to compete with both the aptamer's P1 helix and the expression platform's transcriptional terminator. All-atom molecular dynamics simulations support the hypothesis that ligand binding stabilizes the aptamer P1 helix against central helix strand invasion, thus allowing the terminator to form. These results present a potential model mechanism to explain how ligand binding can induce downstream conformational changes by influencing local strand displacement processes of intermediate folds that could be at play in multiple riboswitch classes.

G-quadruplexes (G4s) are four-stranded nucleic acid structures formed by the stacking of G-tetrads. Here we investigated their formation and function during HIV-1 infection. Using bioinformatics and biophysics analyses we first searched for evolutionary conserved G4-forming sequences in HIV-1 genome. We identified 10 G4s with conservation rates higher than those of HIV-1 regulatory sequences such as RRE and TAR. We then used porphyrin-based G4-binders to probe the formation of the G4s during infection of human cells by native HIV-1. The G4-binders efficiently inhibited HIV-1 infectivity, which is attributed to the formation of G4 structures during HIV-1 replication. Using a qRT-PCR approach, we showed that the formation of viral G4s occurs during the first 2 h post-infection and their stabilization by the G4-binders prevents initiation of reverse transcription. We also used a G4-RNA pull-down approach, based on a G4-specific biotinylated probe, to allow the direct detection and identification of viral G4-RNA in infected cells. Most of the detected G4-RNAs contain crucial regulatory elements such as the PPT and cPPT sequences as well as the U3 region. Hence, these G4s would function in the early stages of infection when the viral RNA genome is being processed for the reverse transcription step.

Exome sequencing is widely used in genetic studies of human diseases and clinical genetic diagnosis. Accurate detection of copy number variants (CNVs) is important to fully utilize exome sequencing data. However, exome data are noisy. None of the existing methods alone can achieve both high precision and recall rate. A common practice is to perform heuristic filtration followed by manual inspection of read depth of putative CNVs. This approach does not scale in large studies. To address this issue, we developed a transfer learning method, CNV-espresso, for in silico confirming rare CNVs from exome sequencing data. CNV-espresso encodes candidate CNVs from exome data as images and uses pretrained convolutional neural network models to classify copy number states. We trained CNV-espresso using an offspring-parents trio exome sequencing dataset, with inherited CNVs as positives and CNVs with Mendelian errors as negatives. We evaluated the performance using additional samples that have both exome and whole-genome sequencing (WGS) data. Assuming the CNVs detected from WGS data as a proxy of ground truth, CNV-espresso significantly improves precision while keeping recall almost intact, especially for CNVs that span a small number of exons. CNV-espresso can effectively replace manual inspection of CNVs in large-scale exome sequencing studies.

In-depth analysis of the transcriptomes of several model organisms has revealed that genomes are pervasively transcribed, giving rise to an abundance of non-canonical and mainly antisense RNA polymerase II-derived transcripts that are produced from almost any genomic context. Pervasive RNAs are degraded by surveillance mechanisms, but the repertoire of proteins that control the fate of these non-productive transcripts is still incomplete. Trypanosomes are single-celled eukaryotes that show constitutive RNA polymerase II transcription and in which initiation and termination of transcription occur at a limited number of sites per chromosome. It is not known whether pervasive transcription exists in organisms with unregulated RNA polymerase II activity, and which factors could be involved in the process. We show here that depletion of RBP33 results in overexpression of ∼40% of all annotated genes in the genome, with a marked accumulation of sense and antisense transcripts derived from silenced regions. RBP33 loss does not result in a significant increase in chromatin accessibility. Finally, we have found that transcripts that increase in abundance upon RBP33 knockdown are significantly more stable in RBP33-depleted trypanosomes, and that the exosome complex is responsible for their degradation. Our results provide strong evidence that RBP33 dampens non-productive transcription in trypanosomes.

Tree- and linear-shaped cell differentiation trajectories have been widely observed in developmental biologies and can be also inferred through computational methods from single-cell RNA-sequencing datasets. However, trajectories with complicated topologies such as loops, disparate lineages and bifurcating hierarchy remain difficult to infer accurately. Here, we introduce a density-based trajectory inference method capable of constructing diverse shapes of topological patterns including the most intriguing bifurcations. The novelty of our method is a step to exploit overlapping probability distributions to identify transition states of cells for determining connectability between cell clusters, and another step to infer a stable trajectory through a base-topology guided iterative fitting. Our method precisely re-constructed various benchmark reference trajectories. As a case study to demonstrate practical usefulness, our method was tested on single-cell RNA sequencing profiles of blood cells of SARS-CoV-2-infected patients. We not only re-discovered the linear trajectory bridging the transition from IgM plasmablast cells to developing neutrophils, and also found a previously-undiscovered lineage which can be rigorously supported by differentially expressed gene analysis.

Meiotic chromosome segregation relies on programmed DNA double-strand break induction. These are in turn repaired by homologous recombination, generating physical attachments between the parental chromosomes called crossovers. A subset of breaks yields recombinant outcomes, while crossover-independent mechanisms repair the majority of lesions. The balance between different repair pathways is crucial to ensure genome integrity. We show that Caenorhabditis elegans BRC-1/BRCA1-BRD-1/BARD1 and PARG-1/PARG form a complex in vivo, essential for accurate DNA repair in the germline. Simultaneous depletion of BRC-1 and PARG-1 causes synthetic lethality due to reduced crossover formation and impaired break repair, evidenced by hindered RPA-1 removal and presence of aberrant chromatin bodies in diakinesis nuclei, whose formation depends on spo-11 function. These factors undergo a similar yet independent loading in developing oocytes, consistent with operating in different pathways. Abrogation of KU- or Theta-mediated end joining elicits opposite effects in brc-1; parg-1 doubles, suggesting a profound impact in influencing DNA repair pathway choice by BRC-1-PARG-1. Importantly, lack of PARG-1 catalytic activity suppresses untimely accumulation of RAD-51 foci in brc-1 mutants but is only partially required for fertility. Our data show that BRC-1/BRD-1-PARG-1 joint function is essential for genome integrity in meiotic cells by regulating multiple DNA repair pathways.

Single-cell ATAC-seq (scATAC-seq) has proven to be a state-of-art approach to investigating gene regulation at the single-cell level. However, existing methods cannot precisely uncover cell-type-specific binding of transcription regulators (TRs) and construct gene regulation networks (GRNs) in single-cell. ChIP-seq has been widely used to profile TR binding sites in the past decades. Here, we developed SCRIP, an integrative method to infer single-cell TR activity and targets based on the integration of scATAC-seq and a large-scale TR ChIP-seq reference. Our method showed improved performance in evaluating TR binding activity compared to the existing motif-based methods and reached a higher consistency with matched TR expressions. Besides, our method enables identifying TR target genes as well as building GRNs at the single-cell resolution based on a regulatory potential model. We demonstrate SCRIP's utility in accurate cell-type clustering, lineage tracing, and inferring cell-type-specific GRNs in multiple biological systems. SCRIP is freely available at https://github.com/wanglabtongji/SCRIP.

Multimodal single-cell sequencing technologies provide unprecedented information on cellular heterogeneity from multiple layers of genomic readouts. However, joint analysis of two modalities without properly handling the noise often leads to overfitting of one modality by the other and worse clustering results than vanilla single-modality analysis. How to efficiently utilize the extra information from single cell multi-omics to delineate cell states and identify meaningful signal remains as a significant computational challenge. In this work, we propose a deep learning framework, named SAILERX, for efficient, robust, and flexible analysis of multi-modal single-cell data. SAILERX consists of a variational autoencoder with invariant representation learning to correct technical noises from sequencing process, and a multimodal data alignment mechanism to integrate information from different modalities. Instead of performing hard alignment by projecting both modalities to a shared latent space, SAILERX encourages the local structures of two modalities measured by pairwise similarities to be similar. This strategy is more robust against overfitting of noises, which facilitates various downstream analysis such as clustering, imputation, and marker gene detection. Furthermore, the invariant representation learning part enables SAILERX to perform integrative analysis on both multi- and single-modal datasets, making it an applicable and scalable tool for more general scenarios.

Erratum in Nucleic Acids Res. 2022 Sep 12;:

Although extended pluripotent stem cells (EPSCs) have the potential to form both embryonic and extraembryonic lineages, how their transcriptional regulatory mechanism differs from that of embryonic stem cells (ESCs) remains unclear. Here, we discovered that YY1 binds to specific open chromatin regions in EPSCs. Yy1 depletion in EPSCs leads to a gene expression pattern more similar to that of ESCs than control EPSCs. Moreover, Yy1 depletion triggers a series of epigenetic crosstalk activities, including changes in DNA methylation, histone modifications and high-order chromatin structures. Yy1 depletion in EPSCs disrupts the enhancer-promoter (EP) interactions of EPSC-specific genes, including Dnmt3l. Yy1 loss results in DNA hypomethylation and dramatically reduces the enrichment of H3K4me3 and H3K27ac on the promoters of EPSC-specific genes by upregulating the expression of Kdm5c and Hdac6 through facilitating the formation of CCCTC-binding factor (CTCF)-mediated EP interactions surrounding their loci. Furthermore, single-cell RNA sequencing (scRNA-seq) experiments revealed that YY1 is required for the derivation of extraembryonic endoderm (XEN)-like cells from EPSCs in vitro. Together, this study reveals that YY1 functions as a key regulator of multidimensional epigenetic crosstalk associated with extended pluripotency.

Identification of protein-coding and non-coding transcripts is paramount for understanding their biological roles. Computational approaches have been addressing this task for over a decade; however, generalized and high-performance models are still unreliable. This benchmark study assessed the performance of 24 tools producing >55 models on the datasets covering a wide range of species. We have collected 135 small and large transcriptomic datasets from existing studies for comparison and identified the potential bottlenecks hampering the performance of current tools. The key insights of this study include lack of standardized training sets, reliance on homogeneous training data, gradual changes in annotated data, lack of augmentation with homology searches, the presence of false positives and negatives in datasets and the lower performance of end-to-end deep learning models. We also derived a new dataset, RNAChallenge, from the benchmark considering hard instances that may include potential false alarms. The best and least well performing models under- and overfit the dataset, respectively, thereby serving a dual purpose. For computational approaches, it will be valuable to develop accurate and unbiased models. The identification of false alarms will be of interest for genome annotators, and experimental study of hard RNAs will help to untangle the complexity of the RNA world.

Bidirectional DNA replication from a chromosome origin requires the asymmetric loading of two helicases, one for each replisome. Our understanding of the molecular mechanisms underpinning helicase loading at bacterial chromosome origins is incomplete. Here we report both positive and negative mechanisms for directing helicase recruitment in the model organism Bacillus subtilis. Systematic characterization of the essential initiation protein DnaD revealed distinct protein interfaces required for homo-oligomerization, interaction with the master initiator protein DnaA, and interaction with the helicase co-loader protein DnaB. Informed by these properties of DnaD, we went on to find that the developmentally expressed repressor of DNA replication initiation, SirA, blocks the interaction between DnaD and DnaA, thereby restricting helicase recruitment from the origin during sporulation to inhibit further initiation events. These results advance our understanding of the mechanisms underpinning DNA replication initiation in B. subtilis, as well as guiding the search for essential cellular activities to target for antimicrobial drug design.

Telomeric repeat-containing RNA (TERRA) is a type of long non-coding RNA transcribed from telomeres, and it forms R-loops by invasion into telomeric DNA. Since either an excessive or inadequate number of R-loops leads to telomere instability, the TERRA levels need to be delicately modulated. In this study, we found that m6A modification presents on the subtelomeric regions of TERRA and stabilizes it, and the loss of METTL3 impacts telomere stability. Mechanically, the m6A modification on TERRA is catalyzed by METTL3, recognized and stabilized by the m6A reader YTHDC1. Knockdown of either METTL3 or YTHDC1 enhances TERRA degradation. The m6A-modified TERRA forms R-loops and promotes homologous recombination which is essential for the alternative lengthening of telomeres (ALT) pathway in cancer cells. METTL3 depletion leads to R-loop reduction, telomere shortening and instability. Altogether, these findings reveal that METTL3 protects telomeres by catalyzing m6A modification on TERRA, indicating that inhibition or deletion of METTL3 is potentially a new avenue for ALT cancer therapy.

Xenobiotic nucleic acids (XNAs) offer tremendous potential for synthetic biology, biotechnology, and molecular medicine but their ability to mimic nucleic acids still needs to be explored. Here, to study the ability of XNA oligonucleotides to mimic tRNA, we synthesized three L-Ala-tXNAs analogs. These molecules were used in a non-ribosomal peptide synthesis involving a bacterial Fem transferase. We compared the ability of this enzyme to use amino-acyl tXNAs containing 1',5'-anhydrohexitol (HNA), 2'-fluoro ribose (2'F-RNA) and 2'-fluoro arabinose. L-Ala-tXNA containing HNA or 2'F-RNA were substrates of the Fem enzyme. The synthesis of peptidyl-XNA and the resolution of their structures in complex with the enzyme show the impact of the XNA on protein binding. For the first time we describe functional tXNA in an in vitro assay. These results invite to test tXNA also as substitute for tRNA in translation.

CRISPR-Cas systems generally provide adaptive immunity in prokaryotes through RNA-guided degradation of foreign genetic elements like bacteriophages and plasmids. Recently, however, transposon-encoded and nuclease-deficient CRISPR-Cas systems were characterized and shown to be co-opted by Tn7-like transposons for CRISPR RNA-guided DNA transposition. As a genome engineering tool, these CRISPR-Cas systems and their associated transposon proteins can be deployed for programmable, site-specific integration of sizable cargo DNA, circumventing the need for DNA cleavage and homology-directed repair involving endogenous repair machinery. Here, we selected a diverse set of type I-F3 CRISPR-associated transposon systems derived from Gammaproteobacteria, predicted all components essential for transposition activity, and deployed them for functionality testing within Escherichia coli. Our results demonstrate that these systems possess a significant range of integration efficiencies with regards to temperature, transposon size, and flexible PAM requirements. Additionally, our findings support the categorization of these systems into functional compatibility groups for efficient and orthogonal RNA-guided DNA integration. This work expands the CRISPR-based toolbox with new CRISPR RNA-guided DNA integrases that can be applied to complex and extensive genome engineering efforts.

Telomeres in Drosophila are composed of sequential non-LTR retrotransposons HeT-A, TART and TAHRE. Although they are repressed by the PIWI-piRNA pathway or heterochromatin in the germline, the regulation of these retrotransposons in somatic cells is poorly understood. In this study, we demonstrated that specific splice variants of Mod(mdg4) repress HeT-A by blocking subtelomeric enhancers in ovarian somatic cells. Among the variants, we found that the Mod(mdg4)-N variant represses HeT-A expression the most efficiently. Subtelomeric sequences bound by Mod(mdg4)-N block enhancer activity within subtelomeric TAS-R repeats. This enhancer-blocking activity is increased by the tandem association of Mod(mdg4)-N to repetitive subtelomeric sequences. In addition, the association of Mod(mdg4)-N couples with the recruitment of RNA polymerase II to the subtelomeres, which reinforces its enhancer-blocking function. Our findings provide novel insights into how telomeric retrotransposons are regulated by the specific variants of insulator proteins associated with subtelomeric sequences.

Comment in 10.1093/nar/gkac910.

Quadruplex-duplex (Q-D) junctions constitute unique structural motifs in genomic sequences. Through comprehensive calorimetric as well as high-resolution NMR structural studies, Q-D junctions with a hairpin-type snapback loop coaxially stacked onto an outer G-tetrad were identified to be most effective binding sites for various polycyclic quadruplex ligands. The Q-D interface is readily recognized by intercalation of the ligand aromatic core structure between G-tetrad and the neighboring base pair. Based on the thermodynamic and structural data, guidelines for the design of ligands with enhanced selectivity towards a Q-D interface emerge. Whereas intercalation at Q-D junctions mostly outcompete stacking at the quadruplex free outer tetrad or intercalation between duplex base pairs to varying degrees, ligand side chains considerably contribute to the selectivity for a Q-D target over other binding sites. In contrast to common perceptions, an appended side chain that additionally interacts within the duplex minor groove may confer only poor selectivity. Rather, the Q-D selectivity is suggested to benefit from an extension of the side chain towards the exposed part of the G-tetrad at the junction. The presented results will support the design of selective high-affinity binding ligands for targeting Q-D interfaces in medicinal but also technological applications.

Massively parallel reporter assay (MPRA) is a high-throughput method that enables the study of the regulatory activities of tens of thousands of DNA oligonucleotides in a single experiment. While MPRA experiments have grown in popularity, their small sample sizes compared to the scale of the human genome limits our understanding of the regulatory effects they detect. To address this, we develop a deep learning model, MpraNet, to distinguish potential MPRA targets from the background genome. This model achieves high discriminative performance (AUROC = 0.85) at differentiating MPRA positives from a set of control variants that mimic the background genome when applied to the lymphoblastoid cell line. We observe that existing functional scores represent very distinct functional effects, and most of them fail to characterize the regulatory effect that MPRA detects. Using MpraNet, we predict potential MPRA functional variants across the genome and identify the distributions of MPRA effect relative to other characteristics of genetic variation, including allele frequency, alternative functional annotations specified by FAVOR, and phenome-wide associations. We also observed that the predicted MPRA positives are not uniformly distributed across the genome; instead, they are clumped together in active regions comprising 9.95% of the genome and inactive regions comprising 89.07% of the genome. Furthermore, we propose our model as a screen to filter MPRA experiment candidates at genome-wide scale, enabling future experiments to be more cost-efficient by increasing precision relative to that observed from previous MPRAs.

PARP1 mediates poly-ADP-ribosylation of proteins on chromatin in response to different types of DNA lesions. PARP inhibitors are used for the treatment of BRCA1/2-deficient breast, ovarian, and prostate cancer. Loss of DNA replication fork protection is proposed as one mechanism that contributes to the vulnerability of BRCA1/2-deficient cells to PARP inhibitors. However, the mechanisms that regulate PARP1 activity at stressed replication forks remain poorly understood. Here, we performed proximity proteomics of PARP1 and isolation of proteins on stressed replication forks to map putative PARP1 regulators. We identified TPX2 as a direct PARP1-binding protein that regulates the auto-ADP-ribosylation activity of PARP1. TPX2 interacts with DNA damage response proteins and promotes homology-directed repair of DNA double-strand breaks. Moreover, TPX2 mRNA levels are increased in BRCA1/2-mutated breast and prostate cancers, and high TPX2 expression levels correlate with the sensitivity of cancer cells to PARP-trapping inhibitors. We propose that TPX2 confers a mitosis-independent function in the cellular response to replication stress by interacting with PARP1.

Breast cancers are known to be driven by the transcription factor estrogen receptor and its ligand estrogen. While the receptor's cis-binding elements are known to vary between tumors, heterogeneity of hormone signaling at a single-cell level is unknown. In this study, we systematically tracked estrogen response across time at a single-cell level in multiple cell line and organoid models. To accurately model these changes, we developed a computational tool (TITAN) that quantifies signaling gradients in single-cell datasets. Using this approach, we found that gene expression response to estrogen is non-uniform, with distinct cell groups expressing divergent transcriptional networks. Pathway analysis suggested the two most distinct signatures are driven separately by ER and FOXM1. We observed that FOXM1 was indeed activated by phosphorylation upon estrogen stimulation and silencing of FOXM1 attenuated the relevant gene signature. Analysis of scRNA-seq data from patient samples confirmed the existence of these divergent cell groups, with the FOXM1 signature predominantly found in ER negative cells. Further, multi-omic single-cell experiments indicated that the different cell groups have distinct chromatin accessibility states. Our results provide a comprehensive insight into ER biology at the single-cell level and potential therapeutic strategies to mitigate resistance to therapy.

Telomere elongation is coupled with genome replication, raising the question of the repair of short telomeres in post-mitotic cells. We investigated the fate of a telomere-repeat capped end that mimics a single short telomere in quiescent fission yeast cells. We show that telomerase is able to elongate this single short telomere during quiescence despite the binding of Ku to the proto-telomere. While Taz1 and Rap1 repress telomerase in vegetative cells, both shelterin proteins are required for efficient telomere extension in quiescent cells, underscoring a distinct mode of telomerase control. We further show that Rad3ATR and Tel1ATM are redundantly required for telomere elongation in quiescence through the phosphorylation of Ccq1 and that Rif1 and its associated-PP1 phosphatases negatively regulate telomerase activity by opposing Ccq1 phosphorylation. The distinct mode of telomerase regulation in quiescent fission yeast cells may be relevant to that in human stem and progenitor cells.

Understanding the function of non-coding genomic sequence variants represents a challenge for biomedicine. Many diseases are products of gene-by-environment interactions with complex mechanisms. This study addresses these themes by mechanistic characterization of non-coding variants that influence gene expression only after drug or hormone exposure. Using glucocorticoid signaling as a model system, we integrated genomic, transcriptomic, and epigenomic approaches to unravel mechanisms by which variant function could be revealed by hormones or drugs. Specifically, we identified cis-regulatory elements and 3D interactions underlying ligand-dependent associations between variants and gene expression. One-quarter of the glucocorticoid-modulated variants that we identified had already been associated with clinical phenotypes. However, their affected genes were 'unmasked' only after glucocorticoid exposure and often with function relevant to the disease phenotypes. These diseases involved glucocorticoids as risk factors or therapeutic agents and included autoimmunity, metabolic and mood disorders, osteoporosis and cancer. For example, we identified a novel breast cancer risk gene, MAST4, with expression that was repressed by glucocorticoids in cells carrying the risk genotype, repression that correlated with MAST4 expression in breast cancer and treatment outcomes. These observations provide a mechanistic framework for understanding non-coding genetic variant-chemical environment interactions and their role in disease risk and drug response.

CRISPR (clustered regularly interspaced short palindromic repeats) technology has achieved great breakthroughs in terms of convenience and sensitivity; it is becoming the most promising molecular tool. However, only two CRISPR activation modes (single and double stranded) are available, and they have specificity and universality bottlenecks that limit the application of CRISPR technology in high-precision molecular recognition. Herein, we proposed a novel CRISPR/Cas12a unrestricted activation mode to greatly improve its performance. The new mode totally eliminates the need for a protospacer adjacent motif and accurately activates Cas12a through toehold-mediated strand displacement and branch migration, which is highly universal and ultra-specific. With this mode, we discriminated all mismatch types and detected the EGFR T790M and L858R mutations in very low abundance. Taken together, our activation mode is deeply incorporated with DNA nanotechnology and extensively broadens the application boundaries of CRISPR technology in biomedical and molecular reaction networks.

Cas12c is the recently characterized dual RNA-guided DNase effector of type V-C CRISPR-Cas (clustered regularly interspaced short palindromic repeats and CRISPR-associated protein) systems. Due to minimal requirements for a protospacer adjacent motif (PAM), Cas12c is an attractive candidate for genome editing. Here we report the crystal structure of Cas12c1 in complex with single guide RNA (sgRNA) and target double-stranded DNA (dsDNA) containing the 5'-TG-3' PAM. Supported by biochemical and mutation assays, this study reveals distinct structural features of Cas12c1 and the associated sgRNA, as well as the molecular basis for PAM recognition, target dsDNA unwinding, heteroduplex formation and recognition, and cleavage of non-target and target DNA strands. Cas12c1 recognizes the PAM through a mechanism that is interdependent on sequence identity and Cas12c1-induced conformational distortion of the PAM region. Another special feature of Cas12c1 is the cleavage of both non-target and target DNA strands at a single, uniform site with indistinguishable cleavage capacity and order. Location of the sgRNA seed region and minimal length of target DNA required for triggering Cas12c1 DNase activity were also determined. Our findings provide valuable information for developing the CRISPR-Cas12c1 system into an efficient, high-fidelity genome editing tool.

Streptomyces are prolific producers of specialized metabolites with applications in medicine and agriculture. These bacteria possess a large linear chromosome genetically compartmentalized: core genes are grouped in the central part, while terminal regions are populated by poorly conserved genes. In exponentially growing cells, chromosome conformation capture unveiled sharp boundaries formed by ribosomal RNA (rrn) operons that segment the chromosome into multiple domains. Here we further explore the link between the genetic distribution of rrn operons and Streptomyces genetic compartmentalization. A large panel of genomes of species representative of the genus diversity revealed that rrn operons and core genes form a central skeleton, the former being identifiable from their core gene environment. We implemented a new nomenclature for Streptomyces genomes and trace their rrn-based evolutionary history. Remarkably, rrn operons are close to pericentric inversions. Moreover, the central compartment delimited by rrn operons has a very dense, nearly invariant core gene content. Finally, this compartment harbors genes with the highest expression levels, regardless of gene persistence and distance to the origin of replication. Our results highlight that rrn operons are structural boundaries of a central functional compartment prone to transcription in Streptomyces.

We previously identified RBPMS as a master regulator of alternative splicing in differentiated smooth muscle cells (SMCs). RBPMS is transcriptionally downregulated during SMC dedifferentiation, but we hypothesized that RBPMS protein activity might be acutely downregulated by post-translational modifications. Publicly available phosphoproteomic datasets reveal that Thr113 and Thr118 immediately adjacent to the RRM domain are commonly both phosphorylated. An RBPMS T113/118 phosphomimetic T/E mutant showed decreased splicing regulatory activity both in transfected cells and in a cell-free in vitro assay, while a non-phosphorylatable T/A mutant retained full activity. Loss of splicing activity was associated with a modest reduction in RNA affinity but significantly reduced RNA binding in nuclear extract. A lower degree of oligomerization of the T/E mutant might cause lower avidity of multivalent RNA binding. However, NMR analysis also revealed that the T113/118E peptide acts as an RNA mimic which can loop back and antagonize RNA-binding by the RRM domain. Finally, we identified ERK2 as the most likely kinase responsible for phosphorylation at Thr113 and Thr118. Collectively, our data identify a potential mechanism for rapid modulation of the SMC splicing program in response to external signals during the vascular injury response and atherogenesis.

Twinkle is a mitochondrial replicative helicase which can self-load onto and unwind mitochondrial DNA. Nearly 60 mutations on Twinkle have been linked to human mitochondrial diseases. Using cryo-electron microscopy (cryo-EM) and high-speed atomic force microscopy (HS-AFM), we obtained the atomic-resolution structure of a vertebrate Twinkle homolog with DNA and captured in real-time how Twinkle is self-loaded onto DNA. Our data highlight the important role of the non-catalytic N-terminal domain of Twinkle. The N-terminal domain directly contacts the C-terminal helicase domain, and the contact interface is a hotspot for disease-related mutations. Mutations at the interface destabilize Twinkle hexamer and reduce helicase activity. With HS-AFM, we observed that a highly dynamic Twinkle domain, which is likely to be the N-terminal domain, can protrude ∼5 nm to transiently capture nearby DNA and initialize Twinkle loading onto DNA. Moreover, structural analysis and subunit doping experiments suggest that Twinkle hydrolyzes ATP stochastically, which is distinct from related helicases from bacteriophages.

In the adaptive immune system, V(D)J recombination initiates the production of a diverse antigen receptor repertoire in developing B and T cells. Recombination activating proteins, RAG1 and RAG2 (RAG1/2), catalyze V(D)J recombination by cleaving adjacent to recombination signal sequences (RSSs) that flank antigen receptor gene segments. Previous studies defined the consensus RSS as containing conserved heptamer and nonamer sequences separated by a less conserved 12 or 23 base-pair spacer sequence. However, many RSSs deviate from the consensus sequence. Here, we developed a cell-based, massively parallel assay to evaluate V(D)J recombination activity on thousands of RSSs where the 12-RSS heptamer and adjoining spacer region contained randomized sequences. While the consensus heptamer sequence (CACAGTG) was marginally preferred, V(D)J recombination was highly active on a wide range of non-consensus sequences. Select purine/pyrimidine motifs that may accommodate heptamer unwinding in the RAG1/2 active site were generally preferred. In addition, while different coding flanks and nonamer sequences affected recombination efficiency, the relative dependency on the purine/pyrimidine motifs in the RSS heptamer remained unchanged. Our results suggest RAG1/2 specificity for RSS heptamers is primarily dictated by DNA structural features dependent on purine/pyrimidine pattern, and to a lesser extent, RAG:RSS base-specific interactions.

Dynamic intron retention (IR) in vertebrate cells is of widespread biological importance. Aberrant IR is associated with numerous human diseases including several cancers. Despite consistent reports demonstrating that intrinsic sequence features can help introns evade splicing, conflicting findings about cell type- or condition-specific IR regulation by trans-regulatory and epigenetic mechanisms demand an unbiased and systematic analysis of IR in a controlled experimental setting. We integrated matched mRNA sequencing (mRNA-Seq), whole-genome bisulfite sequencing (WGBS), nucleosome occupancy methylome sequencing (NOMe-Seq) and chromatin immunoprecipitation sequencing (ChIP-Seq) data from primary human myeloid and lymphoid cells. Using these multi-omics data and machine learning, we trained two complementary models to determine the role of epigenetic factors in the regulation of IR in cells of the innate immune system. We show that increased chromatin accessibility, as revealed by nucleosome-free regions, contributes substantially to the retention of introns in a cell-specific manner. We also confirm that intrinsic characteristics of introns are key for them to evade splicing. This study suggests an important role for chromatin architecture in IR regulation. With an increasing appreciation that pathogenic alterations are linked to RNA processing, our findings may provide useful insights for the development of novel therapeutic approaches that target aberrant splicing.

We present a systematic structural and energetic characterization of phosphate(OP)-nucleobase anion…π stacking interactions in RNAs. We observed OP-nucleobase stacking contacts in a variety of structural motifs other than regular helices and spanning broadly diverse sequence distances. Apart from the stacking between a phosphate and a guanine or a uracil two-residue upstream in specific U-turns, such interactions in RNA have been scarcely characterized to date. Our QM calculations showed an energy minimum at a distance between the OP atom and the nucleobase plane centroid slightly below 3 Å for all the nucleobases. By sliding the OP atom over the nucleobase plane we localized the optimal mutual positioning of the stacked moieties, corresponding to an energy minimum below -6 kcal•mol-1, for all the nucleobases, consistently with the projections of the OP atoms over the different π-rings we observed in experimental occurrences. We also found that the strength of the interaction clearly correlates with its electrostatic component, pointing to it as the most relevant contribution. Finally, as OP-uracil and OP-guanine interactions represent together 86% of the instances we detected, we also proved their stability under dynamic conditions in model systems simulated by state-of-the art DFT-MD calculations.

The conserved RNA helicase UPF1 coordinates nonsense-mediated mRNA decay (NMD) by engaging with mRNAs, RNA decay machinery and the terminating ribosome. UPF1 ATPase activity is implicated in mRNA target discrimination and completion of decay, but the mechanisms through which UPF1 enzymatic activities such as helicase, translocase, RNP remodeling, and ATPase-stimulated dissociation influence NMD remain poorly defined. Using high-throughput biochemical assays to quantify UPF1 enzymatic activities, we show that UPF1 is only moderately processive (<200 nt) in physiological contexts and undergoes ATPase-stimulated dissociation from RNA. We combine an in silico screen with these assays to identify and characterize known and novel UPF1 mutants with altered helicase, ATPase, and RNA binding properties. We find that UPF1 mutants with substantially impaired processivity (E797R, G619K/A546H), faster (G619K) or slower (K547P, E797R, G619K/A546H) unwinding rates, and/or reduced mechanochemical coupling (i.e. the ability to harness ATP hydrolysis for work; K547P, R549S, G619K, G619K/A546H) can still support efficient NMD of well-characterized targets in human cells. These data are consistent with a central role for UPF1 ATPase activity in driving cycles of RNA binding and dissociation to ensure accurate NMD target selection.

Some transcription factors bind DNA motifs containing direct or inverted sequence repeats. Preference for each of these DNA topologies is dictated by structural constraints. Most prokaryotic regulators form symmetric oligomers, which require operators with a dyad structure. Binding to direct repeats requires breaking the internal symmetry, a property restricted to a few regulators, most of them from the AraC family. The KorA family of transcriptional repressors, involved in plasmid propagation and stability, includes members that form symmetric dimers and recognize inverted repeats. Our structural analyses show that ArdK, a member of this family, can form a symmetric dimer similar to that observed for KorA, yet it binds direct sequence repeats as a non-symmetric dimer. This is possible by the 180° rotation of one of the helix-turn-helix domains. We then probed and confirmed that ArdK shows affinity for an inverted repeat, which, surprisingly, is also recognized by a non-symmetrical dimer. Our results indicate that structural flexibility at different positions in the dimerization interface constrains transcription factors to bind DNA sequences with one of these two alternative DNA topologies.

Increasing studies have revealed that a subset of circular RNAs (circRNAs) harbor an open reading frame and can act as protein-coding templates to generate functional proteins that are closely associated with multiple physiological and disease-relevant processes, and thus proper regulation of synthesis of these circRNA-derived proteins is a fundamental cellular process required for homeostasis maintenance. However, how circRNA translation initiation is coordinated by different trans-acting factors remains poorly understood. In particular, the impact of different eukaryotic translation initiation factors (eIFs) on circRNA translation and the physiological relevance of this distinct regulation have not yet been characterized. In this study, we screened all 43 Drosophila eIFs and revealed the conflicting functions of eIF3 subunits in the translational control of the translatable circRNA circSfl: eIF3 is indispensable for circSfl translation, while the eIF3-associated factor eIF3j is the most potent inhibitor. Mechanistically, the binding of eIF3j to circSfl promotes the disassociation of eIF3. The C-terminus of eIF3j and an RNA regulon within the circSfl untranslated region (UTR) are essential for the inhibitory effect of eIF3j. Moreover, we revealed the physiological relevance of eIF3j-mediated circSfl translation repression in response to heat shock. Finally, additional translatable circRNAs were identified to be similarly regulated in an eIF3j-dependent manner. Altogether, our study provides a significant insight into the field of cap-independent translational regulation and undiscovered functions of eIF3.

RNA G-quadruplexes (rG4s) are RNA secondary structures, which are formed by guanine-rich sequences and have important cellular functions. Existing computational tools for rG4 prediction rely on specific sequence features and/or were trained on small datasets, without considering rG4 stability information, and are therefore sub-optimal. Here, we developed rG4detector, a convolutional neural network to identify potential rG4s in transcriptomics data. rG4detector outperforms existing methods in both predicting rG4 stability and in detecting rG4-forming sequences. To demonstrate the biological-relevance of rG4detector, we employed it to study RNAs that are bound by the RNA-binding protein G3BP1. G3BP1 is central to the induction of stress granules (SGs), which are cytoplasmic biomolecular condensates that form in response to a variety of cellular stresses. Unexpectedly, rG4detector revealed a dynamic enrichment of rG4s bound by G3BP1 in response to cellular stress. In addition, we experimentally characterized G3BP1 cross-talk with rG4s, demonstrating that G3BP1 is a bona fide rG4-binding protein and that endogenous rG4s are enriched within SGs. Furthermore, we found that reduced rG4 availability impairs SG formation. Hence, we conclude that rG4s play a direct role in SG biology via their interactions with RNA-binding proteins and that rG4detector is a novel useful tool for rG4 transcriptomics data analyses.

We describe a purified biochemical system to produce monoclonal antibodies (Abs) in vitro using activation-induced deoxycytidine deaminase (AID) and DNA polymerase η (Polη) to diversify immunoglobulin variable gene (IgV) libraries within a phage display format. AID and Polη function during B-cell affinity maturation by catalyzing somatic hypermutation (SHM) of immunoglobulin variable genes (IgV) to generate high-affinity Abs. The IgV mutational motif specificities observed in vivo are conserved in vitro. IgV mutations occurred in antibody complementary determining regions (CDRs) and less frequently in framework (FW) regions. A unique feature of our system is the use of AID and Polη to perform repetitive affinity maturation on libraries reconstructed from a preceding selection step. We have obtained scFv Abs against human glucagon-like peptide-1 receptor (GLP-1R), a target in the treatment of type 2 diabetes, and VHH nanobodies targeting Fatty Acid Amide Hydrolase (FAAH), involved in chronic pain, and artemin, a neurotropic factor that regulates cold pain. A round of in vitro affinity maturation typically resulted in a 2- to 4-fold enhancement in Ab-Ag binding, demonstrating the utility of the system. We tested one of the affinity matured nanobodies and found that it reduced injury-induced cold pain in a mouse model.

Profiling gametes of an individual enables the construction of personalised haplotypes and meiotic crossover landscapes, now achievable at larger scale than ever through the availability of high-throughput single-cell sequencing technologies. However, high-throughput single-gamete data commonly have low depth of coverage per gamete, which challenges existing gamete-based haplotype phasing methods. In addition, haplotyping a large number of single gametes from high-throughput single-cell DNA sequencing data and constructing meiotic crossover profiles using existing methods requires intensive processing. Here, we introduce efficient software tools for the essential tasks of generating personalised haplotypes and calling crossovers in gametes from single-gamete DNA sequencing data (sgcocaller), and constructing, visualising, and comparing individualised crossover landscapes from single gametes (comapr). With additional data pre-possessing, the tools can also be applied to bulk-sequenced samples. We demonstrate that sgcocaller is able to generate impeccable phasing results for high-coverage datasets, on which it is more accurate and stable than existing methods, and also performs well on low-coverage single-gamete sequencing datasets for which current methods fail. Our tools achieve highly accurate results with user-friendly installation, comprehensive documentation, efficient computation times and minimal memory usage.

The honey bee genome has the capacity to produce three phenotypically distinct organisms (two diploid female castes: queen and worker, and a haploid male drone). Previous studies have implicated metabolic flux acting via epigenetic regulation in directing nutrition-driven phenotypic plasticity in the honey bee. However, the cis-acting DNA regulatory elements that establish tissue and polyphenism -specific epigenomes and gene expression programmes, remain unclear. Using a high resolution multiomic approach including assay for transposase-accessible chromatin by sequencing (ATAC-seq), RNA-seq and ChIP-seq, we produce the first genome-wide maps of the regulatory landscape across all three adult honey bee phenotypes identifying > 5000 regulatory regions in queen, 7500 in worker and 6500 in drone, with the vast majority of these sites located within intronic regions. These regions are defined by positive enrichment of H3K27ac and depletion of H3K4me3 and show a positive correlation with gene expression. Using ATAC-seq footprinting we determine queen, worker and drone -specific transcription factor occupancy and uncover novel phenotype-specific regulatory networks identifying two key nuclear receptors that have previously been implicated in caste-determination and adult behavioural maturation in honey bees; ecdysone receptor and ultraspiracle. Collectively, this study provides novel insights into key gene regulatory networks that are associated with these distinct polyphenisms in the honey bee.

Bacterial pathogens employ a plethora of virulence factors for host invasion, and their use is tightly regulated to maximize infection efficiency and manage resources in a nutrient-limited environment. Here we show that during Escherichia coli stationary phase the 3' UTR-derived small non-coding RNA FimR2 regulates fimbrial and flagellar biosynthesis at the post-transcriptional level, leading to biofilm formation as the dominant mode of survival under conditions of nutrient depletion. FimR2 interacts with the translational regulator CsrA, antagonizing its functions and firmly tightening control over motility and biofilm formation. Generated through RNase E cleavage, FimR2 regulates stationary phase biology by fine-tuning target mRNA levels independently of the chaperones Hfq and ProQ. The Salmonella enterica orthologue of FimR2 induces effector protein secretion by the type III secretion system and stimulates infection, thus linking the sRNA to virulence. This work reveals the importance of bacterial sRNAs in modulating various aspects of bacterial physiology including stationary phase and virulence.

Nonsense-mediated RNA decay (NMD) is a highly conserved and selective RNA turnover pathway that depends on the endonuclease SMG6. Here, we show that SMG6 is essential for male germ cell differentiation in mice. Germ-cell conditional knockout (cKO) of Smg6 induces extensive transcriptome misregulation, including a failure to eliminate meiotically expressed transcripts in early haploid cells, and accumulation of NMD target mRNAs with long 3' untranslated regions (UTRs). Loss of SMG6 in the male germline results in complete arrest of spermatogenesis at the early haploid cell stage. We find that SMG6 is strikingly enriched in the chromatoid body (CB), a specialized cytoplasmic granule in male germ cells also harboring PIWI-interacting RNAs (piRNAs) and the piRNA-binding protein PIWIL1. This raises the possibility that SMG6 and the piRNA pathway function together, which is supported by several findings, including that Piwil1-KO mice phenocopy Smg6-cKO mice and that SMG6 and PIWIL1 co-regulate many genes in round spermatids. Together, our results demonstrate that SMG6 is an essential regulator of the male germline transcriptome, and highlight the CB as a molecular platform coordinating RNA regulatory pathways to control sperm production and fertility.

Differentiated oligodendrocytes produce myelin and thereby ensure rapid nerve impulse conduction and efficient information processing in the vertebrate central nervous system. The Krüppel-like transcription factor KLF9 enhances oligodendrocyte differentiation in culture, but appears dispensable in vivo. Its mode of action and role within the oligodendroglial gene regulatory network are unclear. Here we show that KLF9 shares its expression in differentiating oligodendrocytes with the closely related KLF13 protein. Both KLF9 and KLF13 bind to regulatory regions of genes that are important for oligodendrocyte differentiation and equally recognized by the central differentiation promoting transcription factors SOX10 and MYRF. KLF9 and KLF13 physically interact and synergistically activate oligodendrocyte-specific regulatory regions with SOX10 and MYRF. Similar to KLF9, KLF13 promotes differentiation and myelination in primary oligodendroglial cultures. Oligodendrocyte differentiation is also altered in KLF13-deficient mice as demonstrated by a transiently reduced myelin gene expression during the first postnatal week. Considering mouse phenotypes, the similarities in expression pattern and genomic binding and the behaviour in functional assays, KLF9 and KLF13 are important and largely redundant components of the gene regulatory network in charge of oligodendrocyte differentiation and myelination.

The oncogenic Epstein-Barr virus (EBV) evades the immune system but has an Achilles heel: its genome maintenance protein EBNA1. Indeed, EBNA1 is essential for viral genome maintenance but is also highly antigenic. Hence, EBV seemingly evolved a system in which the glycine-alanine repeat (GAr) of EBNA1 limits the translation of its own mRNA to the minimal level to ensure its essential function, thereby, at the same time, minimizing immune recognition. Therefore, defining intervention points at which to interfere with GAr-based inhibition of translation is an important step to trigger an immune response against EBV-carrying cancers. The host protein nucleolin (NCL) plays a critical role in this process via a direct interaction with G-quadruplexes (G4) formed in the GAr-encoding sequence of the viral EBNA1 mRNA. Here we show that the C-terminal arginine-glycine-rich (RGG) motif of NCL is crucial for its role in GAr-based inhibition of translation by mediating interaction of NCL with G4 of EBNA1 mRNA. We also show that this interaction depends on the type I arginine methyltransferase family, notably PRMT1 and PRMT3: drugs or small interfering RNA that target these enzymes prevent efficient binding of NCL on G4 of EBNA1 mRNA and relieve GAr-based inhibition of translation and of antigen presentation. Hence, this work defines type I arginine methyltransferases as therapeutic targets to interfere with EBNA1 and EBV immune evasion.

Gram-positive bacteria are a nascent platform for synthetic biology and metabolic engineering that can provide new opportunities for the production of biomolecules. However, the lack of standardized methods and genetic parts is a major obstacle towards attaining the acceptance and widespread use of Gram-positive bacterial chassis for industrial bioproduction. In this study, we have engineered a novel mRNA leader sequence containing more than one ribosomal binding site (RBS) which could initiate translation from multiple sites, vastly enhancing the translation efficiency of the Gram-positive industrial strain Bacillus licheniformis. This is the first report elucidating the impact of more than one RBS to initiate translation and enhance protein output in B. licheniformis. We also explored the application of more than one RBS for both intracellular and extracellular protein production in B. licheniformis to demonstrate its efficiency, consistency and potential for biotechnological applications. Moreover, we applied these concepts for use in other industrially relevant Gram-positive bacteria, such as Bacillus subtilis and Corynebacterium glutamicum. In all, a highly efficient and robust broad-host expression element has been designed to strengthen and fine-tune the protein outputs for the use of bioproduction in microbial cell factories.

Initiation factor 3 (IF3) regulates the fidelity of bacterial translation initiation by debarring the use of non-canonical start codons or non-initiator tRNAs and prevents premature docking of the 50S ribosomal subunit to the 30S pre-initiation complex (PIC). The C-terminal domain (CTD) of IF3 can carry out most of the known functions of IF3 and sustain Escherichia coli growth. However, the roles of the N-terminal domain (NTD) have remained unclear. We hypothesized that the interaction between NTD and initiator tRNAfMet (i-tRNA) is essential to coordinate the movement of the two domains during the initiation pathway to ensure fidelity of the process. Here, using atomistic molecular dynamics (MD) simulation, we show that R25A/Q33A/R66A mutations do not impact NTD structure but disrupt its interaction with i-tRNA. These NTD residues modulate the fidelity of translation initiation and are crucial for bacterial growth. Our observations also implicate the role of these interactions in the subunit dissociation activity of CTD of IF3. Overall, the study shows that the interactions between NTD of IF3 and i-tRNA are crucial for coupling the movements of NTD and CTD of IF3 during the initiation pathway and in imparting growth fitness to E. coli.

The enteroviral 2C protein is a therapeutic target, but the absence of a mechanistic framework for this enzyme limits our understanding of inhibitor mechanisms. Here, we use poliovirus 2C and a derivative thereof to elucidate the first biochemical mechanism for this enzyme and confirm the applicability of this mechanism to other members of the enterovirus genus. Our biochemical data are consistent with a dimer forming in solution, binding to RNA, which stimulates ATPase activity by increasing the rate of hydrolysis without impacting affinity for ATP substantially. Both RNA and DNA bind to the same or overlapping site on 2C, driven by the phosphodiester backbone, but only RNA stimulates ATP hydrolysis. We propose that RNA binds to 2C driven by the backbone, with reorientation of the ribose hydroxyls occurring in a second step to form the catalytically competent state. 2C also uses a two-step mechanism for binding to ATP. Initial binding is driven by the α and β phosphates of ATP. In the second step, the adenine base and other substituents of ATP are used to organize the active site for catalysis. These studies provide the first biochemical description of determinants driving specificity and catalytic efficiency of a picornaviral 2C ATPase.

The spliceosome undergoes extensive rearrangements as it assembles onto precursor messenger RNAs. In the earliest assembly step, U1snRNA identifies the 5' splice site. However, U1snRNA leaves the spliceosome relatively early in assembly, and 5' splice site identity is subsequently maintained through interactions with U6snRNA, protein factor PRP8, and other components during the rearrangements that build the catalytic site. Using a forward genetic screen in Caenorhabditis elegans, we have identified suppressors of a locomotion defect caused by a 5'ss mutation. Here we report three new suppressor alleles from this screen, two in PRP8 and one in SNRNP200/BRR2. mRNASeq studies of these suppressor strains indicate that they also affect specific native alternative 5'ss, especially for suppressor PRP8 D1549N. A strong suppressor at the unstructured N-terminus of SNRNP200, N18K, indicates a novel role for this region. By examining distinct changes in the splicing of native genes, examining double mutants between suppressors, comparing these new suppressors to previously identified splicing suppressors from yeast, and mapping conserved suppressor residues onto cryoEM structural models of assembling human spliceosomes, we conclude that there are multiple interactions at multiple stages in spliceosome assembly responsible for maintaining the initial 5'ss identified by U1snRNA for entry into the catalytic core.

Current therapies for Duchenne muscular dystrophy (DMD) use phosphorodiamidate morpholino oligomers (PMO) to induce exon skipping in the dystrophin pre-mRNA, enabling the translation of a shortened but functional dystrophin protein. This strategy has been hampered by insufficient delivery of PMO to cardiac and skeletal muscle. To overcome these limitations, we developed the FORCETM platform consisting of an antigen-binding fragment, which binds the transferrin receptor 1, conjugated to an oligonucleotide. We demonstrate that a single dose of the mouse-specific FORCE-M23D conjugate enhances muscle delivery of exon skipping PMO (M23D) in mdx mice, achieving dose-dependent and robust exon skipping and durable dystrophin restoration. FORCE-M23D-induced dystrophin expression reached peaks of 51%, 72%, 62%, 90% and 77%, of wild-type levels in quadriceps, tibialis anterior, gastrocnemius, diaphragm, and heart, respectively, with a single 30 mg/kg PMO-equivalent dose. The shortened dystrophin localized to the sarcolemma, indicating expression of a functional protein. Conversely, a single 30 mg/kg dose of unconjugated M23D displayed poor muscle delivery resulting in marginal levels of exon skipping and dystrophin expression. Importantly, FORCE-M23D treatment resulted in improved functional outcomes compared with administration of unconjugated M23D. Our results suggest that FORCE conjugates are a potentially effective approach for the treatment of DMD.

Comment in 10.1093/nar/gkac961.

Mitochondrial translation is of high significance for cellular energy homeostasis. Aminoacyl-tRNA synthetases (aaRSs) are crucial translational components. Mitochondrial aaRS variants cause various human diseases. However, the pathogenesis of the vast majority of these diseases remains unknown. Here, we identified two novel SARS2 (encoding mitochondrial seryl-tRNA synthetase) variants that cause a multisystem disorder. c.654-14T > A mutation induced mRNA mis-splicing, generating a peptide insertion in the active site; c.1519dupC swapped a critical tRNA-binding motif in the C-terminus due to stop codon readthrough. Both mutants exhibited severely diminished tRNA binding and aminoacylation capacities. A marked reduction in mitochondrial tRNASer(AGY) was observed due to RNA degradation in patient-derived induced pluripotent stem cells (iPSCs), causing impaired translation and comprehensive mitochondrial function deficiencies. These impairments were efficiently rescued by wild-type SARS2 overexpression. Either mutation caused early embryonic fatality in mice. Heterozygous mice displayed reduced muscle tissue-specific levels of tRNASers. Our findings elucidated the biochemical and cellular consequences of impaired translation mediated by SARS2, suggesting that reduced abundance of tRNASer(AGY) is a key determinant for development of SARS2-related diseases.

The chromatin associated with the nuclear lamina (NL) is referred to as lamina-associated domains (LADs). Here, we present an adaptation of the tyramide-signal amplification sequencing (TSA-seq) protocol, which we call chromatin pull down-based TSA-seq (cTSA-seq), that can be used to map chromatin regions at or near the NL from as little as 50 000 cells. The cTSA-seq mapped regions are composed of previously defined LADs and smaller chromatin regions that fall within the Hi-C defined B-compartment containing nuclear peripheral heterochromatin. We used cTSA-seq to map chromatin at or near the assembling NL in cultured cells progressing through early G1. cTSA-seq revealed that the distal ends of chromosomes are near or at the reassembling NL during early G1, a feature similar to those found in senescent cells. We expand the use of cTSA-seq to the mapping of chromatin at or near the NL from fixed-frozen mouse cerebellar tissue sections. This mapping reveals a general conservation of NL-associated chromatin and identifies global and local changes during cerebellar development. The cTSA-seq method reported here is useful for analyzing chromatin at or near the NL from small numbers of cells derived from both in vitro and in vivo sources.

Synthetic mRNA has recently moved into the focus of therapeutic and vaccination efforts. Incorporation of modified nucleotides during in vitro transcription can improve translation and attenuate immunogenicity, but is limited to triphosphate nucleotides which are accepted by RNA polymerases, and their incorporation is either random or complete. In contrast, site-specific modification, herein termed 'point modification' in analogy to point mutations, holds significant technical challenge. We developed fundamental techniques for isolation of long, translatable and internally point-modified mRNAs. Enabling concepts include three-way-one-pot splint ligations, and isolation of mRNA by real-time elution from agarose gels. The use of blue light permitted visualization of mRNA in pre-stained gels without the photochemical damage associated with the use of hard UV-radiation. This allowed visualization of the mRNA through its migration in the agarose gel, which in turn, was a prerequisite for its recovery by electroelution into precast troughs. Co-eluting agarose particles were quantified and found to not be detrimental to mRNA translation in vitro. Translation of EGFP-coding mRNA into functional protein was quantified by incorporation of 35S-labelled methionine and by in-gel EGFP fluorescence. This enabled the functional analysis of point modifications, specifically of ribose methylations in the middle of a 1371 nt long mRNA.

RNA molecules can form secondary and tertiary structures that can regulate their localization and function. Using enzymatic or chemical probing together with high-throughput sequencing, secondary structure can be mapped across the entire transcriptome. However, a limiting factor is that only population averages can be obtained since each read is an independent measurement. Although long-read sequencing has recently been used to determine RNA structure, these methods still used aggregate signals across the strands to detect structure. Averaging across the population also means that only limited information about structural heterogeneity across molecules or dependencies within each molecule can be obtained. Here, we present Single-Molecule Structure sequencing (SMS-seq) that combines structural probing with native RNA sequencing to provide non-amplified, structural profiles of individual molecules with novel analysis methods. Our new approach using mutual information enabled single molecule structural interrogation. Each RNA is probed at numerous bases enabling the discovery of dependencies and heterogeneity of structural features. We also show that SMS-seq can capture tertiary interactions, dynamics of riboswitch ligand binding, and mRNA structural features.

Paraspeckles are ribonucleoprotein granules assembled by NEAT1_2 lncRNA, an isoform of Nuclear Paraspeckle Assembly Transcript 1 (NEAT1). Dysregulation of NEAT1_2/paraspeckles has been linked to multiple human diseases making them an attractive drug target. However currently NEAT1_2/paraspeckle-focused translational research and drug discovery are hindered by a limited toolkit. To fill this gap, we developed and validated a set of tools for the identification of NEAT1_2 binders and modulators comprised of biochemical and cell-based assays. The NEAT1_2 triple helix stability element was utilized as the target in the biochemical assays, and the cellular assay ('ParaQuant') was based on high-content imaging of NEAT1_2 in fixed cells. As a proof of principle, these assays were used to screen a 1,200-compound FDA-approved drug library and a 170-compound kinase inhibitor library and to confirm the screening hits. The assays are simple to establish, use only commercially-available reagents and are scalable for higher throughput. In particular, ParaQuant is a cost-efficient assay suitable for any cells growing in adherent culture and amenable to multiplexing. Using ParaQuant, we identified dual PI3K/mTOR inhibitors as potent negative modulators of paraspeckles. The tools we describe herein should boost paraspeckle studies and help guide the search, validation and optimization of NEAT1_2/paraspeckle-targeted small molecules.

Tandem repeats of simple sequence motifs, also known as microsatellites, are abundant in the genome. Because their repeat structure makes replication error-prone, variant microsatellite lengths are often generated during germline and other somatic expansions. As such, microsatellite length variations can serve as markers for cancer. However, accurate error-free measurement of microsatellite lengths is difficult with current methods precisely because of this high error rate during amplification. We have solved this problem by using partial mutagenesis to disrupt enough of the repeat structure of initial templates so that their sequence lengths replicate faithfully. In this work, we use bisulfite mutagenesis to convert a C to a U, later read as T. Compared to untreated templates, we achieve three orders of magnitude reduction in the error rate per round of replication. By requiring agreement from two independent first copies of an initial template, we reach error rates below one in a million. We apply this method to a thousand microsatellite loci from the human genome, revealing microsatellite length distributions not observable without mutagenesis.

RNA 2'O-methylation is a 'self' epitranscriptomic modification allowing discrimination between host and pathogen. Indeed, human immunodeficiency virus 1 (HIV-1) induces 2'O-methylation of its genome by recruiting the cellular FTSJ3 methyltransferase, thereby impairing detection by RIG-like receptors. Here, we show that RNA 2'O-methylations interfere with the antiviral activity of interferon-stimulated gene 20-kDa protein (ISG20). Biochemical experiments showed that ISG20-mediated degradation of 2'O-methylated RNA pauses two nucleotides upstream of and at the methylated residue. Structure-function analysis indicated that this inhibition is due to steric clash between ISG20 R53 and D90 residues and the 2'O-methylated nucleotide. We confirmed that hypomethylated HIV-1 genomes produced in FTSJ3-KO cells were more prone to in vitro degradation by ISG20 than those produced in cells expressing FTSJ3. Finally, we found that reverse-transcription of hypomethylated HIV-1 was impaired in T cells by interferon-induced ISG20, demonstrating the direct antagonist effect of 2'O-methylation on ISG20-mediated antiviral activity.