Proliferating cell nuclear antigen (PCNA) is a sliding clamp that acts as a central co-ordinator for mismatch repair (MMR) as well as DNA replication. Loss of Elg1, the major subunit of the PCNA unloader complex, causes over-accumulation of PCNA on DNA and also increases mutation rate, but it has been unclear if the two effects are linked. Here we show that timely removal of PCNA from DNA by the Elg1 complex is important to prevent mutations. Although premature unloading of PCNA generally increases mutation rate, the mutator phenotype of elg1Δ is attenuated by PCNA mutants PCNA-R14E and PCNA-D150E that spontaneously fall off DNA. In contrast, the elg1Δ mutator phenotype is exacerbated by PCNA mutants that accumulate on DNA due to enhanced electrostatic PCNA-DNA interactions. Epistasis analysis suggests that PCNA over-accumulation on DNA interferes with both MMR and MMR-independent process(es). In elg1Δ, over-retained PCNA hyper-recruits the Msh2-Msh6 mismatch recognition complex through its PCNA-interacting peptide motif, causing accumulation of MMR intermediates. Our results suggest that PCNA retention controlled by the Elg1 complex is critical for efficient MMR: PCNA needs to be on DNA long enough to enable MMR, but if it is retained too long it interferes with downstream repair steps.

Apurinic/apyrimidinic (AP) sites are a class of highly mutagenic and toxic DNA lesions arising in the genome from a number of exogenous and endogenous sources. Repair of AP lesions takes place predominantly by the base excision pathway (BER). However, among chemically heterogeneous AP lesions formed in DNA, some are resistant to the endonuclease APE1 and thus refractory to BER. Here, we employed two types of reporter constructs accommodating synthetic APE1-resistant AP lesions to investigate the auxiliary repair mechanisms in human cells. By combined analyses of recovery of the transcription rate and suppression of transcriptional mutagenesis at specifically positioned AP lesions, we demonstrate that nucleotide excision repair pathway (NER) efficiently removes BER-resistant AP lesions and significantly enhances the repair of APE1-sensitive ones. Our results further indicate that core NER components XPA and XPF are equally required and that both global genome (GG-NER) and transcription coupled (TC-NER) subpathways contribute to the repair.

RNA molecules are able to bind proteins, DNA and other small or long RNAs using information at primary, secondary or tertiary structure level. Recent techniques that use cross-linking and immunoprecipitation of RNAs can detect these interactions and, if followed by high-throughput sequencing, molecules can be analysed to find recurrent elements shared by interactors, such as sequence and/or structure motifs. Many tools are able to find sequence motifs from lists of target RNAs, while others focus on structure using different approaches to find specific interaction elements. In this work, we make a systematic analysis of RBP-RNA and RNA-RNA datasets to better characterize the interaction landscape with information about multi-motifs on the same RNAs. To achieve this goal, we updated our BEAM algorithm to combine both sequence and structure information to create pairs of patterns that model motifs of interaction. This algorithm was applied to several RNA binding proteins and ncRNAs interactors, confirming already known motifs and discovering new ones. This landscape analysis on interaction variability reflects the diversity of target recognition and underlines that often both primary and secondary structure are involved in molecular recognition.

Epigenetic alterations, including 5-methylcytosine (5mC), 5-hydroxymethylcytosine (5hmC) and nucleosome positioning (NP), in cell-free DNA (cfDNA) have been widely observed in human diseases, and many available cfDNA-based epigenome-wide profiles exhibit high sensitivity and specificity in disease detection and classification. However, due to the lack of efficient collection, standardized quality control, and analysis procedures, efficiently integrating and reusing these data remain considerable challenges. Here, we introduce CFEA (http://www.bio-data.cn/CFEA), a cell-free epigenome database dedicated to three types of widely adopted epigenetic modifications (5mC, 5hmC and NP) involved in 27 human diseases. We developed bioinformatic pipelines for quality control and standard data processing and an easy-to-use web interface to facilitate the query, visualization and download of these cell-free epigenome data. We also manually curated related biological and clinical information for each profile, allowing users to better browse and compare cfDNA epigenomes at a specific stage (such as early- or metastasis-stage) of cancer development. CFEA provides a comprehensive and timely resource to the scientific community and supports the development of liquid biopsy-based biomarkers for various human diseases.

Cellular response to oxidative stress is a crucial mechanism that promotes the survival of Pseudomonas aeruginosa during infection. However, the translational regulation of oxidative stress response remains largely unknown. Here, we reveal a tRNA modification-mediated translational response to H2O2 in P. aeruginosa. We demonstrated that the P. aeruginosa trmB gene encodes a tRNA guanine (46)-N7-methyltransferase that catalyzes the formation of m7G46 in the tRNA variable loop. Twenty-three tRNA substrates of TrmB with a guanosine residue at position 46 were identified, including 11 novel tRNA substrates. We showed that loss of trmB had a strong negative effect on the translation of Phe- and Asp-enriched mRNAs. The trmB-mediated m7G modification modulated the expression of the catalase genes katA and katB, which are enriched with Phe/Asp codons at the translational level. In response to H2O2 exposure, the level of m7G modification increased, consistent with the increased translation efficiency of Phe- and Asp-enriched mRNAs. Inactivation of trmB led to decreased KatA and KatB protein abundance and decreased catalase activity, resulting in H2O2-sensitive phenotype. Taken together, our observations reveal a novel role of m7G46 tRNA modification in oxidative stress response through translational regulation of Phe- and Asp-enriched genes, such as katA and katB.

Conventional high-throughput genomic technologies for mapping regulatory element activities in bulk samples such as ChIP-seq, DNase-seq and FAIRE-seq cannot analyze samples with small numbers of cells. The recently developed low-input and single-cell regulome mapping technologies such as ATAC-seq and single-cell ATAC-seq (scATAC-seq) allow analyses of small-cell-number and single-cell samples, but their signals remain highly discrete or noisy. Compared to these regulome mapping technologies, transcriptome profiling by RNA-seq is more widely used. Transcriptome data in single-cell and small-cell-number samples are more continuous and often less noisy. Here, we show that one can globally predict chromatin accessibility and infer regulatory element activities using RNA-seq. Genome-wide chromatin accessibility predicted by RNA-seq from 30 cells can offer better accuracy than ATAC-seq from 500 cells. Predictions based on single-cell RNA-seq (scRNA-seq) can more accurately reconstruct bulk chromatin accessibility than using scATAC-seq. Integrating ATAC-seq with predictions from RNA-seq increases the power and value of both methods. Thus, transcriptome-based prediction provides a new tool for decoding gene regulatory circuitry in samples with limited cell numbers.

In response to a viral infection, the plant's RNA silencing machinery processes viral RNAs into a huge number of small interfering RNAs (siRNAs). However, a very few of these siRNAs actually interfere with viral replication. A reliable approach to identify these immunologically effective siRNAs (esiRNAs) and to define the characteristics underlying their activity has not been available so far. Here, we develop a novel screening approach that enables a rapid functional identification of antiviral esiRNAs. Tests on the efficacy of such identified esiRNAs of a model virus achieved a virtual full protection of plants against a massive subsequent infection in transient applications. We find that the functionality of esiRNAs depends crucially on two properties: the binding affinity to Argonaute proteins and the ability to access the target RNA. The ability to rapidly identify functional esiRNAs could be of great benefit for all RNA silencing-based plant protection measures against viruses and other pathogens.

HMGB1 is a ubiquitous non-histone protein, which biological effects depend on its expression and subcellular location. Inside the nucleus, HMGB1 is engaged in many DNA events such as DNA repair, transcription and telomere maintenance. HMGB1 has been reported to bind preferentially to bent DNA as well as to noncanonical DNA structures like 4-way junctions and, more recently, to G-quadruplexes. These are four-stranded conformations of nucleic acids involved in important cellular processes, including telomere maintenance. In this frame, G-quadruplex recognition by specific proteins represents a key event to modulate physiological or pathological pathways. Herein, to get insights into the telomeric G-quadruplex DNA recognition by HMGB1, we performed detailed biophysical studies complemented with biological analyses. The obtained results provided information about the molecular determinants for the interaction and showed that the structural variability of human telomeric G-quadruplex DNA may have significant implications in HMGB1 recognition. The biological data identified HMGB1 as a telomere-associated protein in both telomerase-positive and -negative tumor cells and showed that HMGB1 gene silencing in such cells induces telomere DNA damage foci. Altogether, these findings provide a deeper understanding of telomeric G-quadruplex recognition by HMGB1 and suggest that this protein could actually represent a new target for cancer therapy.

HMCES and yedK were recently identified as sensors of abasic sites in ssDNA. In this study, we present multiple crystal structures captured in the apo-, nonspecific-substrate-binding, specific-substrate-binding, and product-binding states of yedK. In combination with biochemical data, we unveil the molecular basis of AP site sensing in ssDNA by yedK. Our results indicate that yedK has a strong preference for AP site-containing ssDNA over native ssDNA and that the conserved Glu105 residue is important for identifying AP sites in ssDNA. Moreover, our results reveal that a thiazolidine linkage is formed between yedK and AP sites in ssDNA, with the residues that stabilize the thiazolidine linkage important for the formation of DNA-protein crosslinks between yedK and the AP sites. We propose that our findings offer a unique platform to develop yedK and other SRAP domain-containing proteins as tools for detecting abasic sites in vitro and in vivo.

Chemical similarity searching is a basic research tool that can be used to find small molecules which are similar in shape to known active molecules. Despite its popularity, the retrieval of local molecular features that are critical to functional activity related to target binding often fails. To overcome this limitation, we developed a novel machine learning-based chemical binding similarity score by using various evolutionary relationships of binding targets. The chemical similarity was defined by the probability of chemical compounds binding to identical targets. Comprehensive and heterogeneous multiple target-binding chemical data were integrated into a paired data format and processed using multiple classification similarity-learning models with various levels of target evolutionary information. Encoding evolutionary information to chemical compounds through their binding targets substantially expanded available chemical-target interaction data and significantly improved model performance. The output probability of our integrated model, referred to as ensemble evolutionary chemical binding similarity (ensECBS), was effective for finding hidden chemical relationships. The developed method can serve as a novel chemical similarity tool that uses evolutionarily conserved target binding information.

The association of reverse transcriptases (RTs) with CRISPR-Cas system has recently attracted interest because the RT activity appears to facilitate the RT-dependent acquisition of spacers from RNA molecules. However, our understanding of this spacer acquisition process remains limited. We characterized the in vivo acquisition of spacers mediated by an RT-Cas1 fusion protein linked to a type III-D system from Vibrio vulnificus strain YJ016, and showed that the adaptation module, consisting of the RT-Cas1 fusion, two different Cas2 proteins (A and B) and one of the two CRISPR arrays, was completely functional in a heterologous host. We found that mutations of the active site of the RT domain significantly decreased the acquisition of new spacers and showed that this RT-Cas1-associated adaptation module was able to incorporate spacers from RNA molecules into the CRISPR array. We demonstrated that the two Cas2 proteins of the adaptation module were required for spacer acquisition. Furthermore, we found that several sequence-specific features were required for the acquisition and integration of spacers derived from any region of the genome, with no bias along the 5'and 3'ends of coding sequences. This study provides new insight into the RT-Cas1 fusion protein-mediated acquisition of spacers from RNA molecules.

U2 auxiliary factor 1 (U2AF1) functions in 3'-splice site selection during pre-mRNA processing. Alternative usage of duplicated tandem exons in U2AF1 produces two isoforms, U2AF1a and U2AF1b, but their functional differences are unappreciated due to their homology. Through integrative approaches of genome editing, customized-transcriptome profiling and crosslinking-mediated interactome analyses, we discovered that the expression of U2AF1 isoforms is controlled by mTOR and they exhibit a distinctive molecular profile for the splice site and protein interactomes. Mechanistic dissection of mutually exclusive alternative splicing events revealed that U2AF1 isoforms' inherent differential preferences of nucleotide sequences and their stoichiometry determine the 3'-splice site. Importantly, U2AF1a-driven transcriptomes feature alternative splicing events in the 5'-untranslated region (5'-UTR) that are favorable for translation. These findings unveil distinct roles of duplicated tandem exon-derived U2AF1 isoforms in the regulation of the transcriptome and suggest U2AF1a-driven 5'-UTR alternative splicing as a molecular mechanism of mTOR-regulated translational control.

RNA-modifying enzymes targeting mRNA poly(A) tails are universal regulators of post-transcriptional gene expression programs. Current data suggest that an RNA-binding protein (RBP) directed tug-of-war between tail shortening and re-elongating enzymes operates in the cytoplasm to repress or activate specific mRNA targets. While this concept is widely accepted, it was primarily described in the final meiotic stages of frog oogenesis and relies molecularly on a single class of RBPs, i.e. CPEBs, the deadenylase PARN and cytoplasmic poly(A) polymerase GLD-2. Using the spatial and temporal resolution of female gametogenesis in the nematode C. elegans, we determined the distinct roles of known deadenylases throughout germ cell development and discovered that the Ccr4-Not complex is the main antagonist to GLD-2-mediated mRNA regulation. We find that the Ccr4-Not/GLD-2 balance is critical for essentially all steps of oocyte production and reiteratively employed by various classes of RBPs. Interestingly, its two deadenylase subunits appear to affect mRNAs stage specifically: while a Caf1/GLD-2 antagonism regulates mRNA abundance during all stages of oocyte production, a Ccr4/GLD-2 antagonism regulates oogenesis in an mRNA abundance independent manner. Our combined data suggests that the Ccr4-Not complex represents the evolutionarily conserved molecular opponent to GLD-2 providing an antagonistic framework of gene-specific poly(A)-tail regulation.

As the use of next-generation sequencing (NGS) for the Mendelian diseases diagnosis is expanding, the performance of this method has to be improved in order to achieve higher quality. Typically, performance measures are considered to be designed in the context of each application and, therefore, account for a spectrum of clinically relevant variants. We present EphaGen, a new computational methodology for bioinformatics quality control (QC). Given a single NGS dataset in BAM format and a pre-compiled VCF-file of targeted clinically relevant variants it associates this dataset with a single arbiter parameter. Intrinsically, EphaGen estimates the probability to miss any variant from the defined spectrum within a particular NGS dataset. Such performance measure virtually resembles the diagnostic sensitivity of given NGS dataset. Here we present case studies of the use of EphaGen in context of BRCA1/2 and CFTR sequencing in a series of 14 runs across 43 blood samples and 504 publically available NGS datasets. EphaGen is superior to conventional bioinformatics metrics such as coverage depth and coverage uniformity. We recommend using this software as a QC step in NGS studies in the clinical context. Availability: https://github.com/m4merg/EphaGen or https://hub.docker.com/r/m4merg/ephagen.

Alternative polyadenylation (APA) is emerging as an important layer of gene regulation because the majority of mammalian protein-coding genes contain multiple polyadenylation (pA) sites in their 3' UTR. By alteration of 3' UTR length, APA can considerably affect post-transcriptional gene regulation. Yet, our understanding of APA remains rudimentary. Novel single-cell RNA sequencing (scRNA-seq) techniques allow molecular characterization of different cell types to an unprecedented degree. Notably, the most popular scRNA-seq protocols specifically sequence the 3' end of transcripts. Building on this property, we implemented a method for analysing patterns of APA regulation from such data. Analyzing multiple datasets from diverse tissues, we identified widespread modulation of APA in different cell types resulting in global 3' UTR shortening/lengthening and enhanced cleavage at intronic pA sites. Our results provide a proof-of-concept demonstration that the huge volume of scRNA-seq data that accumulates in the public domain offers a unique resource for the exploration of APA based on a very broad collection of cell types and biological conditions.

Delivery of double-stranded RNA (dsRNA) into animals can silence genes of matching sequence in diverse cell types through mechanisms that have been collectively called RNA interference. In the nematode Caenorhabditis elegans, dsRNA from multiple sources can trigger the amplification of silencing signals. Amplification occurs through the production of small RNAs by two RNA-dependent RNA polymerases (RdRPs) that are thought to be tissue-specific - EGO-1 in the germline and RRF-1 in somatic cells. Here we demonstrate that EGO-1 can compensate for the lack of RRF-1 when dsRNA from neurons is used to silence genes in intestinal cells. However, the lineal origins of cells that can use EGO-1 varies. This variability could be because random sets of cells can either receive different amounts of dsRNA from the same source or use different RdRPs to perform the same function. Variability is masked in wild-type animals, which show extensive silencing by neuronal dsRNA. As a result, cells appear similarly functional despite underlying differences that vary from animal to animal. This functional mosaicism cautions against inferring uniformity of mechanism based on uniformity of outcome. We speculate that functional mosaicism could contribute to escape from targeted therapies and could allow developmental systems to drift over evolutionary time.

Circular RNAs (circRNA) are a unique class of transcripts that can only be identified from sequence alignments spanning discordant junctions, commonly referred to as backsplice junctions (BSJ). Canonical splicing is also linked with circRNA biogenesis either from the parental transcript or internal to the circRNA, and is not fully utilized in circRNA software. Here we present Ularcirc, a software tool that integrates the visualization of both BSJ and forward splicing junctions and provides downstream analysis of selected circRNA candidates. Ularcirc utilizes the output of CIRI, circExplorer, or raw chimeric output of the STAR aligner and assembles BSJ count table to allow multi-sample analysis. We used Ularcirc to identify and characterize circRNA from public and in-house generated data sets and demonstrate how it can be used to (i) discover novel splicing patterns of parental transcripts, (ii) detect internal splicing patterns of circRNA, and (iii) reveal the complexity of BSJ formation. Furthermore, we identify circRNA that have potential open reading frames longer than their linear sequence. Finally, we detected and validated the presence of a novel class of circRNA generated from ApoA4 transcripts whose BSJ derive from multiple non-canonical splicing sites within coding exons. Ularcirc is accessed via https://github.com/VCCRI/Ularcirc.

One molecule of ADP-ribose 1'',2''-cyclic phosphate (Appr>p) is formed during each of the approximately 500 000 tRNA splicing events per Saccharomyces cerevisiae generation. The metabolism of Appr>p remains poorly defined. A cyclic phosphodiesterase (Cpd1p) has been shown to convert Appr>p to ADP-ribose-1''-phosphate (Appr1p). We used a biochemical genomics approach to identify two yeast phosphatases that can convert Appr1p to ADP-ribose: the product of ORF YBR022w (now Poa1p), which is completely unrelated to other known phosphatases; and Hal2p, a known 3'-phosphatase of 5',3'-pAp. Poa1p is highly specific for Appr1p, and thus likely acts on this molecule in vivo. Poa1 has a relatively low K(M) for Appr1p (2.8 microM) and a modest kcat (1.7 min(-1)), but no detectable activity on several other substrates. Furthermore, Poa1p is strongly inhibited by ADP-ribose (K(I), 17 microM), modestly inhibited by other nucleotides containing an ADP-ribose moiety and not inhibited at all by other tested molecules. In contrast, Hal2p is much more active on pAp than on Appr1p, and several other tested molecules were Hal2p substrates or inhibitors. poa1-Delta mutants have no obvious growth defect at different temperatures in rich media, and analysis of yeast extracts suggests that approximately 90% of Appr1p processing activity originates from Poa1p.

HIV-1 replication can be efficiently inhibited by intracellular expression of an siRNA targeting the viral RNA. However, HIV-1 escape variants emerged after prolonged culturing. These RNAi-resistant viruses contain nucleotide substitutions or deletions in or near the targeted sequence. We observed an inverse correlation between the level of resistance and the stability of the siRNA/target-RNA duplex. However, two escape variants showed a higher level of resistance than expected based on the duplex stability. We demonstrate that these mutations induce alternative folding of the RNA such that the target sequence is occluded from binding to the siRNA, resulting in reduced RNAi efficiency. HIV-1 can thus escape from RNAi-mediated inhibition not only through nucleotide substitutions or deletions in the siRNA target sequence, but also through mutations that alter the local RNA secondary structure. The results highlight the enormous genetic flexibility of HIV-1 and provide detailed molecular insight into the sequence specificity of RNAi and the impact of target RNA secondary structure.

We describe a physical mRNA mapping strategy employing fluorescent self-quenching reporter molecules (SQRMs) that facilitates the identification of mRNA sequence accessible for hybridization with antisense nucleic acids in vitro and in vivo, real time. SQRMs are 20-30 base oligodeoxynucleotides with 5-6 bp complementary ends to which a 5' fluorophore and 3' quenching group are attached. Alone, the SQRM complementary ends form a stem that holds the fluorophore and quencher in contact. When the SQRM forms base pairs with its target, the structure separates the fluorophore from the quencher. This event can be reported by fluorescence emission when the fluorophore is excited. The stem-loop of the SQRM suggests that SQRM be made to target natural stem-loop structures formed during mRNA synthesis. The general utility of this method is demonstrated by SQRM identification of targetable sequence within c-myb and bcl-6 mRNA. Corresponding antisense oligonucleotides reduce these gene products in cells.