Poly(ADP-ribose) polymerases (PARPs) facilitate the repair of DNA single-strand breaks (SSBs). When PARPs are inhibited, unrepaired SSBs colliding with replication forks give rise to cytotoxic double-strand breaks. These are normally rescued by homologous recombination (HR), but, in cells with suboptimal HR, PARP inhibition leads to genomic instability and cell death, a phenomenon currently exploited in the therapy of ovarian cancers in BRCA1/2 mutation carriers. In spite of their promise, resistance to PARP inhibitors (PARPis) has already emerged. In order to identify the possible underlying causes of the resistance, we set out to identify the endogenous source of DNA damage that activates PARPs. We argued that if the toxicity of PARPis is indeed caused by unrepaired SSBs, these breaks must arise spontaneously, because PARPis are used as single agents. We now show that a significant contributor to PARPi toxicity is oxygen metabolism. While BRCA1-depleted or -mutated cells were hypersensitive to the clinically approved PARPi olaparib, its toxicity was significantly attenuated by depletion of OGG1 or MYH DNA glycosylases, as well as by treatment with reactive oxygen species scavengers, growth under hypoxic conditions or chemical OGG1 inhibition. Thus, clinical resistance to PARPi therapy may emerge simply through reduced efficiency of oxidative damage repair.

Ribonucleotides represent a threat to DNA genome stability and transmission. Two types of Ribonuclease H (RNase H) excise ribonucleotides when they form part of the DNA strand, or hydrolyse RNA when it base-pairs with DNA in structures termed R-loops. Loss of either RNase H is lethal in mammals, whereas yeast survives the absence of both enzymes. RNase H1 loss is tolerated by the parasite Trypanosoma brucei but no work has examined the function of RNase H2. Here we show that loss of T. brucei RNase H2 (TbRH2A) leads to growth and cell cycle arrest that is concomitant with accumulation of nuclear damage at sites of RNA polymerase (Pol) II transcription initiation, revealing a novel and critical role for RNase H2. Differential gene expression analysis reveals limited overall changes in RNA levels for RNA Pol II genes after TbRH2A loss, but increased perturbation of nucleotide metabolic genes. Finally, we show that TbRH2A loss causes R-loop and DNA damage accumulation in telomeric RNA Pol I transcription sites, also leading to altered gene expression. Thus, we demonstrate separation of function between two nuclear T. brucei RNase H enzymes during RNA Pol II transcription, but overlap in function during RNA Pol I-mediated gene expression during host immune evasion.

The momentum of scRNA-seq datasets prompts for simple and powerful tools exploring their meaningful signatures. Here we present Single-Cell_Signature_Explorer (https://sites.google.com/site/fredsoftwares/products/single-cell-signature-explorer), the first method for qualitative and high-throughput scoring of any gene set-based signature at the single cell level and its visualization using t-SNE or UMAP. By scanning datasets for single or combined signatures, it rapidly maps any multi-gene feature, exemplified here with signatures of cell lineages, biological hallmarks and metabolic pathways in large scRNAseq datasets of human PBMC, melanoma, lung cancer and adult testis.

DNA nanotechnology allows for the design of programmable DNA-built nanodevices which controllably interact with biological membranes and even mimic the function of natural membrane proteins. Hydrophobic modifications, covalently linked to the DNA, are essential for targeted interfacing of DNA nanostructures with lipid membranes. However, these hydrophobic tags typically induce undesired aggregation eliminating structural control, the primary advantage of DNA nanotechnology. Here, we study the aggregation of cholesterol-modified DNA nanostructures using a combined approach of non-denaturing polyacrylamide gel electrophoresis, dynamic light scattering, confocal microscopy and atomistic molecular dynamics simulations. We show that the aggregation of cholesterol-tagged ssDNA is sequence-dependent, while for assembled DNA constructs, the number and position of the cholesterol tags are the dominating factors. Molecular dynamics simulations of cholesterol-modified ssDNA reveal that the nucleotides wrap around the hydrophobic moiety, shielding it from the environment. Utilizing this behavior, we demonstrate experimentally that the aggregation of cholesterol-modified DNA nanostructures can be controlled by the length of ssDNA overhangs positioned adjacent to the cholesterol. Our easy-to-implement method for tuning cholesterol-mediated aggregation allows for increased control and a closer structure-function relationship of membrane-interfacing DNA constructs - a fundamental prerequisite for employing DNA nanodevices in research and biomedicine.

Microrchidia family CW-type zinc finger 2 (MORC2) is a newly identified chromatin remodeling enzyme with an emerging role in DNA damage response (DDR), but the underlying mechanism remains largely unknown. Here, we show that poly(ADP-ribose) polymerase 1 (PARP1), a key chromatin-associated enzyme responsible for the synthesis of poly(ADP-ribose) (PAR) polymers in mammalian cells, interacts with and PARylates MORC2 at two residues within its conserved CW-type zinc finger domain. Following DNA damage, PARP1 recruits MORC2 to DNA damage sites and catalyzes MORC2 PARylation, which stimulates its ATPase and chromatin remodeling activities. Mutation of PARylation residues in MORC2 results in reduced cell survival after DNA damage. MORC2, in turn, stabilizes PARP1 through enhancing acetyltransferase NAT10-mediated acetylation of PARP1 at lysine 949, which blocks its ubiquitination at the same residue and subsequent degradation by E3 ubiquitin ligase CHFR. Consequently, depletion of MORC2 or expression of an acetylation-defective PARP1 mutant impairs DNA damage-induced PAR production and PAR-dependent recruitment of DNA repair proteins to DNA lesions, leading to enhanced sensitivity to genotoxic stress. Collectively, these findings uncover a previously unrecognized mechanistic link between MORC2 and PARP1 in the regulation of cellular response to DNA damage.

RNA interference represents a potent intervention for cancer treatment but requires a robust delivery agent for transporting gene-modulating molecules, such as small interfering RNAs (siRNAs). Although numerous molecular approaches for siRNA delivery are adequate in vitro, delivery to therapeutic targets in vivo is limited by payload integrity, cell targeting, efficient cell uptake, and membrane penetration. We constructed nonviral biomaterials to transport small nucleic acids to cell targets, including tumor cells, on the basis of the self-assembling and cell-penetrating activities of the adenovirus capsid penton base. Our recombinant penton base chimera contains polypeptide domains designed for noncovalent assembly with anionic molecules and tumor homing. Here, structural modeling, molecular dynamics simulations, and functional assays suggest that it forms pentameric units resembling viral capsomeres that assemble into larger capsid-like structures when combined with siRNA cargo. Pentamerization forms a barrel lined with charged residues mediating pH-responsive dissociation and exposing masked domains, providing insight on the endosomolytic mechanism. The therapeutic impact was examined on tumors expressing high levels of HER3/ErbB3 that are resistant to clinical inhibitors. Our findings suggest that our construct may utilize ligand mimicry to avoid host attack and target the siRNA to HER3+ tumors by forming multivalent capsid-like structures.

Staphylococcus aureus ArlRS is a key two-component regulatory system necessary for adhesion, biofilm formation, and virulence. The response regulator ArlR consists of a C-terminal DNA-binding effector domain and an N-terminal receiver domain that is phosphorylated by ArlS, the cognate transmembrane sensor histidine kinase. We demonstrate that the receiver domain of ArlR adopts the canonical α5β5 response regulator assembly, which dimerizes upon activation, using beryllium trifluoride as an aspartate phosphorylation mimic. Activated ArlR recognizes a 20-bp imperfect inverted repeat sequence in the ica operon, which is involved in intercellular adhesion polysaccharide production. Crystal structures of the inactive and activated forms reveal that activation induces a significant conformational change in the β4-α4 and β5-α5-connecting loops, in which the α4 and α5 helices constitute the homodimerization interface. Crystal structures of the DNA-binding ArlR effector domain indicate that it is able to dimerize via a non-canonical β1-β2 hairpin domain swapping, raising the possibility of a new mechanism for signal transduction from the receiver domain to effector domain. Taken together, the current study provides structural insights into the activation of ArlR and its recognition, adding to the diversity of response regulation mechanisms that may inspire novel antimicrobial strategies specifically targeting Staphylococcus.

CancerGeneNet (https://signor.uniroma2.it/CancerGeneNet/) is a resource that links genes that are frequently mutated in cancers to cancer phenotypes. The resource takes advantage of a curation effort aimed at embedding a large fraction of the gene products that are found altered in cancer cells into a network of causal protein relationships. Graph algorithms, in turn, allow to infer likely paths of causal interactions linking cancer associated genes to cancer phenotypes thus offering a rational framework for the design of strategies to revert disease phenotypes. CancerGeneNet bridges two interaction layers by connecting proteins whose activities are affected by cancer drivers to proteins that impact on the 'hallmarks of cancer'. In addition, CancerGeneNet annotates curated pathways that are relevant to rationalize the pathological consequences of cancer driver mutations in selected common cancers and 'MiniPathways' illustrating regulatory circuits that are frequently altered in different cancers.

The leaf senescence database (LSD) is a comprehensive resource of senescence-associated genes (SAGs) and their corresponding mutants. Through manual curation and extensive annotation, we updated the LSD to a new version LSD 3.0, which contains 5853 genes and 617 mutants from 68 species. To provide sustainable and reliable services for the plant research community, LSD 3.0 (https://bigd.big.ac.cn/lsd/) has been moved to and maintained by the National Genomics Data Center at Beijing Institute of Genomics, Chinese Academy of Sciences. In the current release, we added some new features: (i) Transcriptome data of leaf senescence in poplar were integrated; (ii) Leaf senescence-associated transcriptome data information in Arabidopsis, rice and soybean were included; (iii) Senescence-differentially expressed small RNAs (Sen-smRNA) in Arabidopsis were identified; (iv) Interaction pairs between Sen-smRNAs and senescence-associated transcription factors (Sen-TF) were established; (v) Senescence phenotypes of 90 natural accessions (ecotypes) and 42 images of ecotypes in Arabidopsis were incorporated; (vi) Mutant seed information of SAGs in rice obtained from Kitbase was integrated; (vii) New options of search engines for ecotypes and transcriptome data were implemented. Together, the updated database bears great utility to continue to provide users with useful resources for studies of leaf senescence.

We present a multi-laboratory effort to describe the structural and dynamical properties of duplex B-DNA under physiological conditions. By processing a large amount of atomistic molecular dynamics simulations, we determine the sequence-dependent structural properties of DNA as expressed in the equilibrium distribution of its stochastic dynamics. Our analysis includes a study of first and second moments of the equilibrium distribution, which can be accurately captured by a harmonic model, but with nonlocal sequence-dependence. We characterize the sequence-dependent choreography of backbone and base movements modulating the non-Gaussian or anharmonic effects manifested in the higher moments of the dynamics of the duplex when sampling the equilibrium distribution. Contrary to prior assumptions, such anharmonic deformations are not rare in DNA and can play a significant role in determining DNA conformation within complexes. Polymorphisms in helical geometries are particularly prevalent for certain tetranucleotide sequence contexts and are always coupled to a complex network of coordinated changes in the backbone. The analysis of our simulations, which contain instances of all tetranucleotide sequences, allow us to extend Calladine-Dickerson rules used for decades to interpret the average geometry of DNA, leading to a set of rules with quantitative predictive power that encompass nonlocal sequence-dependence and anharmonic fluctuations.

The binding of p53-binding protein 1 (53BP1) to damaged chromatin is a critical event in non-homologous DNA end joining (NHEJ)-mediated DNA damage repair. Although several molecular pathways explaining how 53BP1 binds damaged chromatin have been described, the precise underlying mechanisms are still unclear. Here we report that a newly identified H4K16 monomethylation (H4K16me1) mark is involved in 53BP1 binding activity in the DNA damage response (DDR). During the DDR, H4K16me1 rapidly increases as a result of catalyzation by the histone methyltransferase G9a-like protein (GLP). H4K16me1 shows an increased interaction level with 53BP1, which is important for the timely recruitment of 53BP1 to DNA double-strand breaks. Differing from H4K16 acetylation, H4K16me1 enhances the 53BP1-H4K20me2 interaction at damaged chromatin. Consistently, GLP knockdown markedly attenuates 53BP1 foci formation, leading to impaired NHEJ-mediated repair and decreased cell survival. Together, these data support a novel axis of the DNA damage repair pathway based on H4K16me1 catalysis by GLP, which promotes 53BP1 recruitment to permit NHEJ-mediated DNA damage repair.

The GWAS Central resource provides a toolkit for integrative access and visualization of a uniquely extensive collection of genome-wide association study data, while ensuring safe open access to prevent research participant identification. GWAS Central is the world's most comprehensive openly accessible repository of summary-level GWAS association information, providing over 70 million P-values for over 3800 studies investigating over 1400 unique phenotypes. The database content comprises direct submissions received from GWAS authors and consortia, in addition to actively gathered data sets from various public sources. GWAS data are discoverable from the perspective of genetic markers, genes, genome regions or phenotypes, via graphical visualizations and detailed downloadable data reports. Tested genetic markers and relevant genomic features can be visually interrogated across up to sixteen multiple association data sets in a single view using the integrated genome browser. The semantic standardization of phenotype descriptions with Medical Subject Headings and the Human Phenotype Ontology allows the precise identification of genetic variants associated with diseases, phenotypes and traits of interest. Harmonization of the phenotype descriptions used across several GWAS-related resources has extended the phenotype search capabilities to enable cross-database study discovery using a range of ontologies. GWAS Central is updated regularly and available at https://www.gwascentral.org.

The nucleolus is a nuclear sub-domain containing the most highly transcribed genes in the genome. Hundreds of human ribosomal RNA (rRNA) genes, located in the nucleolus, rely on constant maintenance. DNA double-strand breaks (DSBs) in rRNA genes activate the ATM kinase, repress rRNA transcription and induce nucleolar cap formation. Yet how ribosomal-DNA (rDNA) lesions are detected and processed remains elusive. Here, we use CRISPR/Cas9-mediated induction of DSBs and report a chromatin response unique to rDNA depending on ATM-phosphorylation of the nucleolar protein TCOF1 and recruitment of the MRE11-RAD50-NBS1 (MRN) complex via the NBS1-subunit. NBS1- and MRE11-depleted cells fail to suppress rRNA transcription and to translocate rDNA into nucleolar caps. Furthermore, the DNA damage response (DDR) kinase ATR operates downstream of the ATM-TCOF1-MRN interplay and is required to fully suppress rRNA transcription and complete DSB-induced nucleolar restructuring. Unexpectedly, we find that DSBs in rDNA neither activate checkpoint kinases CHK1/CHK2 nor halt cell-cycle progression, yet the nucleolar-DDR protects against genomic aberrations and cell death. Our data highlight the concept of a specialized nucleolar DNA damage response (n-DDR) with a distinct protein composition, spatial organization and checkpoint communication. The n-DDR maintains integrity of ribosomal RNA genes, with implications for cell physiology and disease.

Ribosome biogenesis is an essential process in all living cells, which entails countless highly sequential and dynamic structural reorganization events. These include formation of dozens RNA helices through Watson-Crick base-pairing within ribosomal RNAs (rRNAs) and between rRNAs and small nucleolar RNAs (snoRNAs), transient association of hundreds of proteinaceous assembly factors to nascent precursor (pre-)ribosomes, and stable assembly of ribosomal proteins. Unsurprisingly, the largest group of ribosome assembly factors are energy-consuming proteins (NTPases) including 25 RNA helicases in budding yeast. Among these, the DEAH-box Dhr1 is essential to displace the box C/D snoRNA U3 from the pre-rRNAs where it is bound in order to prevent premature formation of the central pseudoknot, a dramatic irreversible long-range interaction essential to the overall folding of the small ribosomal subunit. Here, we report the crystal structure of the Dhr1 helicase module, revealing the presence of a remarkable carboxyl-terminal domain essential for Dhr1 function in ribosome biogenesis in vivo and important for its interaction with its coactivator Utp14 in vitro. Furthermore, we report the functional consequences on ribosome biogenesis of DHX37 (human Dhr1) mutations found in patients suffering from microcephaly and other neurological diseases.

Argonaute (Ago) proteins are key players in RNA interference in eukaryotes, where they function as RNA-guided RNA endonucleases. Prokaryotic Argonautes (pAgos) are much more diverse than their eukaryotic counterparts but their cellular functions and mechanisms of action remain largely unknown. Some pAgos were shown to use small DNA guides for endonucleolytic cleavage of complementary DNA in vitro. However, previously studied pAgos from thermophilic prokaryotes function at elevated temperatures, which limits their potential use as a tool in genomic applications. Here, we describe two pAgos from mesophilic bacteria, Clostridium butyricum (CbAgo) and Limnothrix rosea (LrAgo), that act as DNA-guided DNA nucleases at physiological temperatures. In comparison with previously studied pAgos, CbAgo and LrAgo do not show strong preferences for the 5'-nucleotide in guide DNA and can use not only 5'-phosphorylated but also 5'-hydroxyl DNA guides. Both CbAgo and LrAgo can tolerate guide/target mismatches in the seed region, but are sensitive to mismatches in the 3'-guide region. Both pAgos can perform programmable endonucleolytic cleavage of double-stranded DNA substrates, showing enhanced activity at AT-rich regions and at elevated temperatures. The biochemical characterization of mesophilic pAgo proteins paves the way for their use for DNA manipulations both in vitro and in vivo.

The frequency by which genes are studied correlates with the prior knowledge accumulated about them. This leads to an imbalance in research attention where some genes are highly investigated while others are ignored. Geneshot is a search engine developed to illuminate this gap and to promote attention to the under-studied genome. Through a simple web interface, Geneshot enables researchers to enter arbitrary search terms, to receive ranked lists of genes relevant to the search terms. Returned ranked gene lists contain genes that were previously published in association with the search terms, as well as genes predicted to be associated with the terms based on data integration from multiple sources. The search results are presented with interactive visualizations. To predict gene function, Geneshot utilizes gene-gene similarity matrices from processed RNA-seq data, or from gene-gene co-occurrence data obtained from multiple sources. In addition, Geneshot can be used to analyze the novelty of gene sets and augment gene sets with additional relevant genes. The Geneshot web-server and API are freely and openly available from https://amp.pharm.mssm.edu/geneshot.

Identifying active prophages is critical for studying coevolution of phage and bacteria, investigating phage physiology and biochemistry, and engineering designer phages for diverse applications. We present Prophage Hunter, a tool aimed at hunting for active prophages from whole genome assembly of bacteria. Combining sequence similarity-based matching and genetic features-based machine learning classification, we developed a novel scoring system that exhibits higher accuracy than current tools in predicting active prophages on the validation datasets. The option of skipping similarity matching is also available so that there's higher chance for novel phages to be discovered. Prophage Hunter provides a one-stop web service to extract prophage genomes from bacterial genomes, evaluate the activity of the prophages, identify phylogenetically related phages, and annotate the function of phage proteins. Prophage Hunter is freely available at https://pro-hunter.bgi.com/.

The structure, dynamics and, ultimately, biological function of proteins and nucleic acids are determined by the physicochemical properties of their primary sequences. Such properties are frequently captured via one-dimensional profile plots depicting a given physicochemical variable as a function of sequence position. Hydrophobicity, charge or structural disorder in proteins or nucleobase-density in nucleic acids are routinely visualized in this manner to analyze sequences at a glance. Such visualizations, however, are typically created case-by-case in a purely static manner, employ fixed visualization parameters only and do not enable a quantitative comparison between different sequences. Here, we present VOLPES (volpes.univie.ac.at), a user-friendly web server and the corresponding JavaScript library that enable a fully interactive, multifunctional visualization, analysis and comparison of the physicochemical properties of protein and nucleic-acid sequences, allowing unprecedented insight into biological sequence data and creating a starting point for further in-depth exploration.

Transfer RNAs (tRNAs) are divided into two types, type I with a short variable loop and type II with a long variable loop. Aminoacylation of type I or type II tRNALeu is catalyzed by their cognate leucyl-tRNA synthetases (LeuRSs). However, in Streptomyces coelicolor, there are two types of tRNALeu and only one LeuRS (ScoLeuRS). We found that the enzyme could leucylate both types of ScotRNALeu, and had a higher catalytic efficiency for type II ScotRNALeu(UAA) than for type I ScotRNALeu(CAA). The results from tRNA and enzyme mutagenesis showed that ScoLeuRS did not interact with the canonical discriminator A73. The number of nucleotides, rather than the type of base of the variable loop in the two types of ScotRNALeus, was determined as important for aminoacylation. In vitro and in vivo assays showed that the tertiary structure formed by the D-loop and TψC-loop is more important for ScotRNALeu(UAA). We showed that the leucine-specific domain (LSD) of ScoLeuRS could help LeuRS, which originally only leucylates type II tRNALeu, to aminoacylate type I ScotRNALeu(CAA) and identified the crucial amino acid residues at the C-terminus of the LSD to recognize type I ScotRNALeu(CAA). Overall, our findings identified a rare recognition mechanism of LeuRS to tRNALeu.

WebGestalt is a popular tool for the interpretation of gene lists derived from large scale -omics studies. In the 2019 update, WebGestalt supports 12 organisms, 342 gene identifiers and 155 175 functional categories, as well as user-uploaded functional databases. To address the growing and unique need for phosphoproteomics data interpretation, we have implemented phosphosite set analysis to identify important kinases from phosphoproteomics data. We have completely redesigned result visualizations and user interfaces to improve user-friendliness and to provide multiple types of interactive and publication-ready figures. To facilitate comprehension of the enrichment results, we have implemented two methods to reduce redundancy between enriched gene sets. We introduced a web API for other applications to get data programmatically from the WebGestalt server or pass data to WebGestalt for analysis. We also wrapped the core computation into an R package called WebGestaltR for users to perform analysis locally or in third party workflows. WebGestalt can be freely accessed at http://www.webgestalt.org.