We previously reported that human Rev1 (hRev1) bound to a parallel-stranded G-quadruplex (G4) from the c-MYC promoter with high affinity. We have extended those results to include other G4 motifs, finding that hRev1 exhibited stronger affinity for parallel-stranded G4 than either anti-parallel or hybrid folds. Amino acids in the αE helix of insert-2 were identified as being important for G4 binding. Mutating E466 and Y470 to alanine selectively perturbed G4 binding affinity. The E466K mutant restored wild-type G4 binding properties. Using a forward mutagenesis assay, we discovered that loss of hRev1 increased G4 mutation frequency >200-fold compared to the control sequence. Base substitutions and deletions occurred around and within the G4 motif. Pyridostatin (PDS) exacerbated this effect, as the mutation frequency increased >700-fold over control and deletions upstream of the G4 site more than doubled. Mutagenic replication of G4 DNA (±PDS) was partially rescued by wild-type and E466K hRev1. The E466A or Y470A mutants failed to suppress the PDS-induced increase in G4 mutation frequency. These findings have implications for the role of insert-2, a motif conserved in vertebrates but not yeast or plants, in Rev1-mediated suppression of mutagenesis during G4 replication.

Antibiotic resistance is a serious problem which may be caused by bacterial dormancy. It has been suggested that bacterial toxin-antitoxin systems induce dormancy. We analyzed the genome-wide role of Staphylococcus aureus endoribonuclease toxin MazF using RNA-Seq, Ribo-Seq and quantitative proteomics. We characterized changes in transcriptome, translatome and proteome caused by MazF, and proposed that MazF decreases translation directly by cleaving mRNAs, and indirectly, by decreasing translation factors and by promoting ribosome hibernation. Important pathways affected during the early stage of MazF induction were identified: MazF increases cell wall thickness and decreases cell division; MazF activates SsrA-system which rescues stalled ribosomes, appearing as a result of MazF mRNA cleavage. These pathways may be promising targets for new antibacterial drugs that prevent bacteria dormancy. Finally, we described the overall impact of MazF on S. aureus cell physiology, and propose one of the mechanisms by which MazF might regulate cellular changes leading to dormancy.

A set of >300 nonredundant high-resolution RNA-protein complexes were rigorously searched for π-contacts between an amino acid side chain (W, H, F, Y, R, E and D) and an RNA nucleobase (denoted π-π interaction) or ribose moiety (denoted sugar-π). The resulting dataset of >1500 RNA-protein π-contacts were visually inspected and classified based on the interaction type, and amino acids and RNA components involved. More than 80% of structures searched contained at least one RNA-protein π-interaction, with π-π contacts making up 59% of the identified interactions. RNA-protein π-π and sugar-π contacts exhibit a range in the RNA and protein components involved, relative monomer orientations and quantum mechanically predicted binding energies. Interestingly, π-π and sugar-π interactions occur more frequently with RNA (4.8 contacts/structure) than DNA (2.6). Moreover, the maximum stability is greater for RNA-protein contacts than DNA-protein interactions. In addition to highlighting distinct differences between RNA and DNA-protein binding, this work has generated the largest dataset of RNA-protein π-interactions to date, thereby underscoring that RNA-protein π-contacts are ubiquitous in nature, and key to the stability and function of RNA-protein complexes.

We recently showed that site-specific incorporation of 2'-modifications or neutral linkages in the oligo-deoxynucleotide gap region of toxic phosphorothioate (PS) gapmer ASOs can enhance therapeutic index and safety. In this manuscript, we determined if introducing substitution at the 5'-position of deoxynucleotide monomers in the gap can also enhance therapeutic index. Introducing R- or S-configured 5'-Me DNA at positions 3 and 4 in the oligodeoxynucleotide gap enhanced the therapeutic profile of the modified ASOs suggesting a different positional preference as compared to the 2'-OMe gap modification strategy. The generality of these observations was demonstrated by evaluating R-5'-Me and R-5'-Ethyl DNA modifications in multiple ASOs targeting HDAC2, FXI and Dynamin2 mRNA in the liver. The current work adds to a growing body of evidence that small structural changes can modulate the therapeutic properties of PS ASOs and ushers a new era of chemical optimization with a focus on enhancing the therapeutic profile as opposed to nuclease stability, RNA-affinity and pharmacokinetic properties. The 5'-methyl DNA modified ASOs exhibited excellent safety and antisense activity in mice highlighting the therapeutic potential of this class of nucleic acid analogs for next generation ASO designs.

Glucocorticoid receptor (GR) is an essential transcription factor (TF), controlling metabolism, development and immune responses. SUMOylation regulates chromatin occupancy and target gene expression of GR in a locus-selective manner, but the mechanism of regulation has remained elusive. Here, we identify the protein network around chromatin-bound GR by using selective isolation of chromatin-associated proteins and show that the network is affected by receptor SUMOylation, with several nuclear receptor coregulators and chromatin modifiers preferring interaction with SUMOylation-deficient GR and proteins implicated in transcriptional repression preferring interaction with SUMOylation-competent GR. This difference is reflected in our chromatin binding, chromatin accessibility and gene expression data, showing that the SUMOylation-deficient GR is more potent in binding and opening chromatin at glucocorticoid-regulated enhancers and inducing expression of target loci. Blockage of SUMOylation by a SUMO-activating enzyme inhibitor (ML-792) phenocopied to a large extent the consequences of GR SUMOylation deficiency on chromatin binding and target gene expression. Our results thus show that SUMOylation modulates the specificity of GR by regulating its chromatin protein network and accessibility at GR-bound enhancers. We speculate that many other SUMOylated TFs utilize a similar regulatory mechanism.

GNRA tetraloop-binding receptor interactions are key components in the macromolecular assembly of a variety of functional RNAs. In nature, there is an apparent bias for GAAA/11nt receptor and GYRA/helix interactions, with the former interaction being thermodynamically more stable than the latter. While past in vitro selections allowed isolation of novel GGAA and GUGA receptors, we report herein an in vitro selection that revealed several novel classes of specific GUAA receptors with binding affinities comparable to those from natural GAAA/11nt interactions. These GUAA receptors have structural homology with double-locked bulge RNA modules naturally occurring in ribosomal RNAs. They display mutational robustness that enables exploration of the sequence/phenotypic space associated to GNRA/receptor interactions through epistasis. Their thermodynamic self-assembly fitness landscape is characterized by a rugged neutral network with possible evolutionary trajectories toward natural GNRA/receptor interactions. High throughput sequencing analysis revealed synergetic mutations located away from the tertiary interactions that positively contribute to assembly fitness. Our study suggests that the repertoire of GNRA/receptor interactions is much larger than initially thought from the analysis of natural stable RNA molecules and also provides clues for their evolution towards natural GNRA/receptors.

Short non-coding RNA molecules (sRNAs) play a fundamental role in gene regulation and development in higher organisms. They act as molecular postcodes and guide AGO proteins to target nucleic acids. In plants, sRNA-targeted mRNAs are degraded, reducing gene expression. In contrast, sRNA-targeted DNA sequences undergo cytosine methylation referred to as RNA-directed DNA methylation (RdDM). Cytosine methylation can suppress transcription, thus sRNAs are potent regulators of gene expression. sRNA-mediated RdDM is involved in genome stability through transposon silencing, mobile signalling for epigenetic gene control and hybrid vigour. Since cytosine methylation can be passed on to subsequent generations, RdDM contributes to transgenerational inheritance of the epigenome. Using a novel approach, which can differentiate between primary (inducer) and secondary (amplified) sRNAs, we show that initiation of heritable RdDM does not require complete sequence complementarity between the sRNAs and their nuclear target sequences. sRNAs with up to four regularly interspaced mismatches are potent inducers of RdDM, however, the number and disruptive nature of nucleotide polymorphisms negatively correlate with their efficacy. Our findings contribute to understanding how sRNA can directly shape the epigenome and may be used in designing the next generation of RNA silencing constructs.

PARP-1 is a key early responder to DNA damage in eukaryotic cells. An allosteric mechanism links initial sensing of DNA single-strand breaks by PARP-1's F1 and F2 domains via a process of further domain assembly to activation of the catalytic domain (CAT); synthesis and attachment of poly(ADP-ribose) (PAR) chains to protein sidechains then signals for assembly of DNA repair components. A key component in transmission of the allosteric signal is the HD subdomain of CAT, which alone bridges between the assembled DNA-binding domains and the active site in the ART subdomain of CAT. Here we present a study of isolated CAT domain from human PARP-1, using NMR-based dynamics experiments to analyse WT apo-protein as well as a set of inhibitor complexes (with veliparib, olaparib, talazoparib and EB-47) and point mutants (L713F, L765A and L765F), together with new crystal structures of the free CAT domain and inhibitor complexes. Variations in both dynamics and structures amongst these species point to a model for full-length PARP-1 activation where first DNA binding and then substrate interaction successively destabilise the folded structure of the HD subdomain to the point where its steric blockade of the active site is released and PAR synthesis can proceed.

Regulatory protein access to the DNA duplex 'interior' depends on local DNA 'breathing' fluctuations, and the most fundamental of these are thermally-driven base stacking-unstacking interactions. The smallest DNA unit that can undergo such transitions is the dinucleotide, whose structural and dynamic properties are dominated by stacking, while the ion condensation, cooperative stacking and inter-base hydrogen-bonding present in duplex DNA are not involved. We use dApdA to study stacking-unstacking at the dinucleotide level because the fluctuations observed are likely to resemble those of larger DNA molecules, but in the absence of constraints introduced by cooperativity are likely to be more pronounced, and thus more accessible to measurement. We study these fluctuations with a combination of Molecular Dynamics simulations on the microsecond timescale and Markov State Model analyses, and validate our results by calculations of circular dichroism (CD) spectra, with results that agree well with the experimental spectra. Our analyses show that the CD spectrum of dApdA is defined by two distinct chiral conformations that correspond, respectively, to a Watson-Crick form and a hybrid form with one base in a Hoogsteen configuration. We find also that ionic structure and water orientation around dApdA play important roles in controlling its breathing fluctuations.

In Staphylococcus aureus, de novo methionine biosynthesis is regulated by a unique hierarchical pathway involving stringent-response controlled CodY repression in combination with a T-box riboswitch and RNA decay. The T-box riboswitch residing in the 5' untranslated region (met leader RNA) of the S. aureus metICFE-mdh operon controls downstream gene transcription upon interaction with uncharged methionyl-tRNA. met leader and metICFE-mdh (m)RNAs undergo RNase-mediated degradation in a process whose molecular details are poorly understood. Here we determined the secondary structure of the met leader RNA and found the element to harbor, beyond other conserved T-box riboswitch structural features, a terminator helix which is target for RNase III endoribonucleolytic cleavage. As the terminator is a thermodynamically highly stable structure, it also forms posttranscriptionally in met leader/ metICFE-mdh read-through transcripts. Cleavage by RNase III releases the met leader from metICFE-mdh mRNA and initiates RNase J-mediated degradation of the mRNA from the 5'-end. Of note, metICFE-mdh mRNA stability varies over the length of the transcript with a longer lifespan towards the 3'-end. The obtained data suggest that coordinated RNA decay represents another checkpoint in a complex regulatory network that adjusts costly methionine biosynthesis to current metabolic requirements.

Following antigenic challenge, activated B cells rapidly expand and undergo somatic hypermutation, yielding groups of clonally related B cells with diversified immunoglobulin receptors. Inference of clonal relationships based on the receptor sequence is an essential step in many adaptive immune receptor repertoire sequencing studies. These relationships are typically identified by a multi-step process that involves: (i) grouping sequences based on shared V and J gene assignments, and junction lengths and (ii) clustering these sequences using a junction-based distance. However, this approach is sensitive to the initial gene assignments, which are error-prone, and fails to identify clonal relatives whose junction length has changed through accumulation of indels. Through defining a translation-invariant feature space in which we cluster the sequences, we develop an alignment free clonal identification method that does not require gene assignments and is not restricted to a fixed junction length. This alignment free approach has higher sensitivity compared to a typical junction-based distance method without loss of specificity and PPV. While the alignment free procedure identifies clones that are broadly consistent with the junction-based distance method, it also identifies clones with characteristics (multiple V or J gene assignments or junction lengths) that are not detectable with the junction-based distance method.

In arthropods, zinc finger-associated domains (ZADs) are found at the N-termini of many DNA-binding proteins with tandem arrays of Cys2-His2 zinc fingers (ZAD-C2H2 proteins). ZAD-C2H2 proteins undergo fast evolutionary lineage-specific expansion and functional diversification. Here, we show that all ZADs from Drosophila melanogaster form homodimers, but only certain ZADs with high homology can also heterodimerize. CG2712, for example, is unable to heterodimerize with its paralog, the previously characterized insulator protein Zw5, with which it shares 46% homology. We obtained a crystal structure of CG2712 protein's ZAD domain that, in spite of a low sequence homology, has similar spatial organization with the only known ZAD structure (from Grauzone protein). Steric clashes prevented the formation of heterodimers between Grauzone and CG2712 ZADs. Using detailed structural analysis, site-directed mutagenesis, and molecular dynamics simulations, we demonstrated that rapid evolutionary acquisition of interaction specificity was mediated by the more energy-favorable formation of homodimers in comparison to heterodimers, and that this specificity was achieved by multiple amino acid substitutions resulting in the formation or breaking of stabilizing interactions. We speculate that specific homodimerization of ZAD-C2H2 proteins is important for their architectural role in genome organization.

RcsB is a transcriptional regulator that controls expression of numerous genes in enteric bacteria. RcsB accomplishes this role alone or in combination with auxiliary transcriptional factors independently or dependently of phosphorylation. To understand the mechanisms by which RcsB regulates such large number of genes, we performed structural studies as well as in vitro and in vivo functional studies with different RcsB variants. Our structural data reveal that RcsB binds promoters of target genes such as rprA and flhDC in a dimeric active conformation. In this state, the RcsB homodimer docks the DNA-binding domains into the major groove of the DNA, facilitating an initial weak read-out of the target sequence. Interestingly, comparative structural analyses also show that DNA binding may stabilize an active conformation in unphosphorylated RcsB. Furthermore, RNAseq performed in strains expressing wild-type or several RcsB variants provided new insights into the contribution of phosphorylation to gene regulation and assign a potential role of RcsB in controlling iron metabolism. Finally, we delimited the RcsB box for homodimeric active binding to DNA as the sequence TN(G/A)GAN4TC(T/C)NA. This RcsB box was found in promoter, intergenic and intragenic regions, facilitating both increased or decreased gene transcription.

Long noncoding RNAs (lncRNAs) play a wide range of roles in the epigenetic regulation of crucial biological processes, but the functions of lncRNAs in cortical development are poorly understood. Using human embryonic stem cell (hESC)-based 2D neural differentiation approach and 3D cerebral organoid system, we identified that the lncRNA PAUPAR, which is adjacent to PAX6, plays essential roles in cortical differentiation by interacting with PAX6 to regulate the expression of a large number of neural genes. Mechanistic studies showed that PAUPAR confers PAX6 proper binding sites on the target neural genes by directly binding the genomic regions of these genes. Moreover, PAX6 recruits the histone methyltransferase NSD1 through its C-terminal PST enrichment domain, then regulate H3K36 methylation and the expression of target genes. Collectively, our data reveal that the PAUPAR/PAX6/NSD1 complex plays a critical role in the epigenetic regulation of hESC cortical differentiation and highlight the importance of PAUPAR as an intrinsic regulator of cortical differentiation.

Bacteria deploy multiple defenses to prevent mobile genetic element (MGEs) invasion. CRISPR-Cas immune systems use RNA-guided nucleases to target MGEs, which counter with anti-CRISPR (Acr) proteins. Our understanding of the biology and co-evolutionary dynamics of the common Type I-C CRISPR-Cas subtype has lagged because it lacks an in vivo phage-host model system. Here, we show the anti-phage function of a Pseudomonas aeruginosa Type I-C CRISPR-Cas system encoded on a conjugative pKLC102 island, and its Acr-mediated inhibition by distinct MGEs. Seven genes with anti-Type I-C function (acrIC genes) were identified, many with highly acidic amino acid content, including previously described DNA mimic AcrIF2. Four of the acr genes were broad spectrum, also inhibiting I-E or I-F P. aeruginosa CRISPR-Cas subtypes. Dual inhibition comes at a cost, however, as simultaneous expression of Type I-C and I-F systems renders phages expressing the dual inhibitor AcrIF2 more sensitive to targeting. Mutagenesis of numerous acidic residues in AcrIF2 did not impair anti-I-C or anti-I-F function per se but did exacerbate inhibition defects during competition, suggesting that excess negative charge may buffer DNA mimics against competition. Like AcrIF2, five of the Acr proteins block Cascade from binding DNA, while two function downstream, likely preventing Cas3 recruitment or activity. One such inhibitor, AcrIC3, is found in an 'anti-Cas3' cluster within conjugative elements, encoded alongside bona fide Cas3 inhibitors AcrIF3 and AcrIE1. Our findings demonstrate an active battle between an MGE-encoded CRISPR-Cas system and its diverse MGE targets.

The integrity and regulation of the nuclear lamina is essential for nuclear organization and chromatin stability, with its dysregulation being linked to laminopathy diseases and cancer. Although numerous posttranslational modifications have been identified on lamins, few have been ascribed a regulatory function. Here, we establish that lamin B1 (LMNB1) acetylation at K134 is a molecular toggle that controls nuclear periphery stability, cell cycle progression, and DNA repair. LMNB1 acetylation prevents lamina disruption during herpesvirus type 1 (HSV-1) infection, thereby inhibiting virus production. We also demonstrate the broad impact of this site on laminar processes in uninfected cells. LMNB1 acetylation negatively regulates canonical nonhomologous end joining by impairing the recruitment of 53BP1 to damaged DNA. This defect causes a delay in DNA damage resolution and a persistent activation of the G1/S checkpoint. Altogether, we reveal LMNB1 acetylation as a mechanism for controlling DNA repair pathway choice and stabilizing the nuclear periphery.

Vast G-quadruplexes (GQs) are primarily folded by one, two, or four G-rich oligomers, rarely with an exception. Here, we present the first NMR solution structure of a trimolecular GQ (tri-GQ) that is solely assembled by the self-trimerization of d(GTTAGG), preferentially in Na+ solution tolerant to an equal amount of K+ cation. Eight guanines from three asymmetrically folded strands of d(GTTAGG) are organized into a two-tetrad core, which features a broken G-column and two width-irregular grooves. Fast strand exchanges on a timescale of second at 17°C spontaneously occur between folded tri-GQ and unfolded single-strand of d(GTTAGG) that both species coexist in dynamic equilibrium. Thus, this tri-GQ is not just simply a static assembly but rather a dynamic assembly. Moreover, another minor tetra-GQ that has putatively tetrameric (2+2) antiparallel topology becomes noticeable only at an extremely high strand concentration above 18 mM. The major tri-GQ and minor tetra-GQ are considered to be mutually related, and their reversible interconversion pathways are proposed accordingly. The sequence d(GTTAGG) could be regarded as either a reading frame shifted single repeat of human telomeric DNA or a 1.5 repeat of Bombyx mori telomeric DNA. Overall, our findings provide new insight into GQs and expect more functional applications.

Animal models are crucial for advancing our knowledge about the molecular pathways involved in human diseases. However, it remains unclear to what extent tissue expression of pathways in healthy individuals is conserved between species. In addition, organism-specific information on pathways in animal models is often lacking. Within these limitations, we explore the possibilities that arise from publicly available data for the animal models mouse, rat, and pig. We approximate the animal pathways activity by integrating the human counterparts of curated pathways with tissue expression data from the models. Specifically, we compare whether the animal orthologs of the human genes are expressed in the same tissue. This is complicated by the lower coverage and worse quality of data in rat and pig as compared to mouse. Despite that, from 203 human KEGG pathways and the seven tissues with best experimental coverage, we identify 95 distinct pathways, for which the tissue expression in one animal model agrees better with human than the others. Our systematic pathway-tissue comparison between human and three animal modes points to specific similarities with human and to distinct differences among the animal models, thereby suggesting the most suitable organism for modeling a human pathway or tissue.

During homologous recombination, Dbl2 protein is required for localisation of Fbh1, an F-box helicase that efficiently dismantles Rad51-DNA filaments. RNA-seq analysis of dbl2Δ transcriptome showed that the dbl2 deletion results in upregulation of more than 500 loci in Schizosaccharomyces pombe. Compared with the loci with no change in expression, the misregulated loci in dbl2Δ are closer to long terminal and long tandem repeats. Furthermore, the misregulated loci overlap with antisense transcripts, retrotransposons, meiotic genes and genes located in subtelomeric regions. A comparison of the expression profiles revealed that Dbl2 represses the same type of genes as the HIRA histone chaperone complex. Although dbl2 deletion does not alleviate centromeric or telomeric silencing, it suppresses the silencing defect at the outer centromere caused by deletion of hip1 and slm9 genes encoding subunits of the HIRA complex. Moreover, our analyses revealed that cells lacking dbl2 show a slight increase of nucleosomes at transcription start sites and increased levels of methylated histone H3 (H3K9me2) at centromeres, subtelomeres, rDNA regions and long terminal repeats. Finally, we show that other proteins involved in homologous recombination, such as Fbh1, Rad51, Mus81 and Rad54, participate in the same gene repression pathway.

Erratum in Nucleic Acids Res. 2021 Apr 6;49(6):3599. Gerdes, Kenn [added].

The 3D genome is characterized by a complex organization made of genomic and epigenomic layers with profound implications on gene regulation and cell function. However, the understanding of the fundamental mechanisms driving the crosstalk between nuclear architecture and (epi)genomic information is still lacking. The plant Arabidopsis thaliana is a powerful model organism to address these questions owing to its compact genome for which we have a rich collection of microscopy, chromosome conformation capture (Hi-C) and ChIP-seq experiments. Using polymer modelling, we investigate the roles of nucleolus formation and epigenomics-driven interactions in shaping the 3D genome of A. thaliana. By validation of several predictions with published data, we demonstrate that self-attracting nucleolar organizing regions and repulsive constitutive heterochromatin are major mechanisms to regulate the organization of chromosomes. Simulations also suggest that interphase chromosomes maintain a partial structural memory of the V-shapes, typical of (sub)metacentric chromosomes in anaphase. Additionally, self-attraction between facultative heterochromatin regions facilitates the formation of Polycomb bodies hosting H3K27me3-enriched gene-clusters. Since nucleolus and heterochromatin are highly-conserved in eukaryotic cells, our findings pave the way for a comprehensive characterization of the generic principles that are likely to shape and regulate the 3D genome in many species.

Methods for the detection of m6A by RNA-Seq technologies are increasingly sought after. We here present NOseq, a method to detect m6A residues in defined amplicons by virtue of their resistance to chemical deamination, effected by nitrous acid. Partial deamination in NOseq affects all exocyclic amino groups present in nucleobases and thus also changes sequence information. The method uses a mapping algorithm specifically adapted to the sequence degeneration caused by deamination events. Thus, m6A sites with partial modification levels of ∼50% were detected in defined amplicons, and this threshold can be lowered to ∼10% by combination with m6A immunoprecipitation. NOseq faithfully detected known m6A sites in human rRNA, and the long non-coding RNA MALAT1, and positively validated several m6A candidate sites, drawn from miCLIP data with an m6A antibody, in the transcriptome of Drosophila melanogaster. Conceptually related to bisulfite sequencing, NOseq presents a novel amplicon-based sequencing approach for the validation of m6A sites in defined sequences.

New genetic tools and strategies are currently under development to facilitate functional genomics analyses. Here, we describe an active member of the Tc1/mariner transposon superfamily, named ZB, which invaded the zebrafish genome very recently. ZB exhibits high activity in vertebrate cells, in the range of those of the widely used transposons piggyBac (PB), Sleeping Beauty (SB) and Tol2. ZB has a similar structural organization and target site sequence preference to SB, but a different integration profile with respect to genome-wide preference among mammalian functional annotation features. Namely, ZB displays a preference for integration into transcriptional regulatory regions of genes. Accordingly, we demonstrate the utility of ZB for enhancer trapping in zebrafish embryos and in the mouse germline. These results indicate that ZB may be a powerful tool for genetic manipulation in vertebrate model species.

G-quadruplexes are noncanonical nucleic acid structures formed by stacked guanine tetrads. They are capable of a range of functions and thought to play widespread biological roles. This diversity raises an important question: what determines the biochemical specificity of G-quadruplex structures? The answer is particularly important from the perspective of biological regulation because genomes can contain hundreds of thousands of G-quadruplexes with a range of functions. Here we analyze the specificity of each sequence in a 496-member library of variants of a reference G-quadruplex with respect to five functions. Our analysis shows that the sequence requirements of G-quadruplexes with these functions are different from one another, with some mutations altering biochemical specificity by orders of magnitude. Mutations in tetrads have larger effects than mutations in loops, and changes in specificity are correlated with changes in multimeric state. To complement our biochemical data we determined the solution structure of a monomeric G-quadruplex from the library. The stacked and accessible tetrads rationalize why monomers tend to promote a model peroxidase reaction and generate fluorescence. Our experiments support a model in which the sequence requirements of G-quadruplexes with different functions are overlapping but distinct. This has implications for biological regulation, bioinformatics, and drug design.

Replication forks often stall at damaged DNA. To overcome these obstructions and complete the DNA duplication in a timely fashion, replication can be restarted downstream of the DNA lesion. In mammalian cells, this repriming of replication can be achieved through the activities of primase and polymerase PrimPol. PrimPol is stimulated in DNA synthesis through interaction with PolDIP2, however the exact mechanism of this PolDIP2-dependent stimulation is still unclear. Here, we show that PrimPol uses a flexible loop to interact with the C-terminal ApaG-like domain of PolDIP2, and that this contact is essential for PrimPol's enhanced processivity. PolDIP2 increases primer-template and dNTP binding affinities of PrimPol, which concomitantly enhances its nucleotide incorporation efficiency. This stimulation is dependent on a unique arginine cluster in PolDIP2. Since the polymerase activity of PrimPol alone is very limited, this mechanism, where the affinity for dNTPs gets increased by PolDIP2 binding, might be critical for the in vivo function of PrimPol in tolerating DNA lesions at physiological nucleotide concentrations.

Interferon regulatory factor 4 (IRF4) is a key transcription factor (TF) in the regulation of immune cells, including B and T cells. It acts by binding DNA as both a homodimer and, in conjunction with other TFs, as a heterodimer. The choice of homo and heterodimeric/ DNA interactions is a critical aspect in the control of the transcriptional program and cell fate outcome. To characterize the nature of this interaction in the homodimeric complex, we have determined the crystal structure of the IRF4/ISRE homodimeric complex. We show that the complex formation is aided by a substantial DNA deformation with co-operative binding achieved exclusively through protein-DNA contact. This markedly contrasts with the heterodimeric form where DNA bound IRF4 is shown to physically interact with PU.1 TF to engage EICE1. We also show that the hotspot residues (Arg98, Cys99 and Asn102) contact both consensus and non-consensus sequences with the L1 loop exhibiting marked flexibility. Additionally, we identified that IRF4L116R, a mutant associated with chronic lymphocytic leukemia, binds more robustly to DNA thereby providing a rationale for the observed gain of function. Together, we demonstrate key structural differences between IRF4 homo and heterodimeric complexes, thereby providing molecular insights into IRF4-mediated transcriptional regulation.

We recently showed that Saccharomyces cerevisiae telomeric DNA can fold into an unprecedented pseudocircular G-hairpin (PGH) structure. However, the formation of PGHs in the context of extended sequences, which is a prerequisite for their function in vivo and their applications in biotechnology, has not been elucidated. Here, we show that despite its 'circular' nature, PGHs tolerate single-stranded (ss) protrusions. High-resolution NMR structure of a novel member of PGH family reveals the atomistic details on a junction between ssDNA and PGH unit. Identification of new sequences capable of folding into one of the two forms of PGH helped in defining minimal sequence requirements for their formation. Our time-resolved NMR data indicate a possibility that PGHs fold via a complex kinetic partitioning mechanism and suggests the existence of K+ ion-dependent PGH folding intermediates. The data not only provide an explanation of cation-type-dependent formation of PGHs, but also explain the unusually large hysteresis between PGH melting and annealing noted in our previous study. Our findings have important implications for DNA biology and nanotechnology. Overrepresentation of sequences able to form PGHs in the evolutionary-conserved regions of the human genome implies their functionally important biological role(s).

Ribosome stalling at tandem CGA codons or poly(A) sequences activates quality controls for nascent polypeptides including ribosome-associated quality control (RQC) and no-go mRNA decay (NGD). In RQC pathway, Hel2-dependent uS10 ubiquitination and the RQC-trigger (RQT) complex are essential for subunit dissociation, and Ltn1-dependent ubiquitination of peptidyl-tRNA in the 60S subunit requires Rqc2. Here, we report that polytryptophan sequences induce Rqc2-independent RQC. More than 11 consecutive tryptophan residues induced RQC in a manner dependent on Hel2-mediated ribosome ubiquitination and the RQT complex. Polytryptophan sequence-mediated RQC was not coupled with CAT-tailing, and Rqc2 was not required for Ltn1-dependent degradation of the arrest products. Eight consecutive tryptophan residues located at the region proximal to the peptidyl transferase center in the ribosome tunnel inhibited CAT-tailing by tandem CGA codons. Polytryptophan sequences also induced Hel2-mediated canonical RQC-coupled NGD and RQC-uncoupled NGD outside the stalled ribosomes. We propose that poly-tryptophan sequences induce Rqc2-independent RQC, suggesting that CAT-tailing in the 60S subunit could be modulated by the polypeptide in the ribosome exit tunnel.

Maintenance of stem-cell identity requires proper regulation of enhancer activity. Both transcription factors OCT4/SOX2/NANOG and histone methyltransferase complexes MLL/SET1 were shown to regulate enhancer activity, but how they are regulated in embryonic stem cells (ESCs) remains further studies. Here, we report a transcription factor BACH1, which directly interacts with OCT4/SOX2/NANOG (OSN) and MLL/SET1 methyltransferase complexes and maintains pluripotency in mouse ESCs (mESCs). BTB domain and bZIP domain of BACH1 are required for these interactions and pluripotency maintenance. Loss of BACH1 reduced the interaction between NANOG and MLL1/SET1 complexes, and decreased their occupancy on chromatin, and further decreased H3 lysine 4 trimethylation (H3K4me3) level on gene promoters and (super-) enhancers, leading to decreased enhancer activity and transcription activity, especially on stemness-related genes. Moreover, BACH1 recruited NANOG through chromatin looping and regulated remote NANOG binding, fine-tuning enhancer-promoter activity and gene expression. Collectively, these observations suggest that BACH1 maintains pluripotency in ESCs by recruiting NANOG and MLL/SET1 complexes to chromatin and maintaining the trimethylated state of H3K4 and enhancer-promoter activity, especially on stemness-related genes.

Dysfunction of splicing factors often result in abnormal cell differentiation and apoptosis, especially in neural tissues. Mutations in pre-mRNAs processing factor 31 (PRPF31) cause autosomal dominant retinitis pigmentosa, a progressive retinal degeneration disease. The transcriptome-wide splicing events specifically regulated by PRPF31 and their biological roles in the development and maintenance of retina are still unclear. Here, we showed that the differentiation and viability of retinal progenitor cells (RPCs) are severely perturbed in prpf31 knockout zebrafish when compared with other tissues at an early embryonic stage. At the cellular level, significant mitotic arrest and DNA damage were observed. These defects could be rescued by the wild-type human PRPF31 rather than the disease-associated mutants. Further bioinformatic analysis and experimental verification uncovered that Prpf31 deletion predominantly causes the skipping of exons with a weak 5' splicing site. Moreover, genes necessary for DNA repair and mitotic progression are most enriched among the differentially spliced events, which may explain the cellular and tissular defects in prpf31 mutant retinas. This is the first time that Prpf31 is demonstrated to be essential for the survival and differentiation of RPCs during retinal neurogenesis by specifically modulating the alternative splicing of genes involved in DNA repair and mitosis.

Escherichia coli RecO is a recombination mediator protein that functions in the RecF pathway of homologous recombination, in concert with RecR, and interacts with E. coli single stranded (ss) DNA binding (SSB) protein via the last 9 amino acids of the C-terminal tails (SSB-Ct). Structures of the E. coli RecR and RecOR complexes are unavailable; however, crystal structures from other organisms show differences in RecR oligomeric state and RecO stoichiometry. We report analytical ultracentrifugation studies of E. coli RecR assembly and its interaction with RecO for a range of solution conditions using both sedimentation velocity and equilibrium approaches. We find that RecR exists in a pH-dependent dimer-tetramer equilibrium that explains the different assembly states reported in previous studies. RecO binds with positive cooperativity to a RecR tetramer, forming both RecR4O and RecR4O2 complexes. We find no evidence of a stable RecO complex with RecR dimers. However, binding of RecO to SSB-Ct peptides elicits an allosteric effect, eliminating the positive cooperativity and shifting the equilibrium to favor a RecR4O complex. These studies suggest a mechanism for how SSB binding to RecO influences the distribution of RecOR complexes to facilitate loading of RecA onto SSB coated ssDNA to initiate homologous recombination.

Sensitive detection of microsatellite instability (MSI) in tissue or liquid biopsies using next generation sequencing (NGS) has growing prognostic and predictive applications in cancer. However, the complexities of NGS make it cumbersome as compared to established multiplex-PCR detection of MSI. We present a new approach to detect MSI using inter-Alu-PCR followed by targeted NGS, that combines the practical advantages of multiplexed-PCR with the breadth of information provided by NGS. Inter-Alu-PCR employs poly-adenine repeats of variable length present in every Alu element and provides a massively-parallel, rapid approach to capture poly-A-rich genomic fractions within short 80-150bp amplicons generated from adjacent Alu-sequences. A custom-made software analysis tool, MSI-tracer, enables Alu-associated MSI detection from tissue biopsies or MSI-tracing at low-levels in circulating-DNA. MSI-associated indels at somatic-indel frequencies of 0.05-1.5% can be detected depending on the availability of matching normal tissue and the extent of instability. Due to the high Alu copy-number in human genomes, a single inter-Alu-PCR retrieves enough information for identification of MSI-associated-indels from ∼100 pg circulating-DNA, reducing current limits by ∼2-orders of magnitude and equivalent to circulating-DNA obtained from finger-sticks. The combined practical and informational advantages of inter-Alu-PCR make it a powerful tool for identifying tissue-MSI-status or tracing MSI-associated-indels in liquid biopsies.

CRISPR-Cas are adaptive immune systems that degrade foreign genetic elements in archaea and bacteria. In carrying out their immune functions, CRISPR-Cas systems heavily rely on RNA components. These CRISPR (cr) RNAs are repeat-spacer units that are produced by processing of pre-crRNA, the transcript of CRISPR arrays, and guide Cas protein(s) to the cognate invading nucleic acids, enabling their destruction. Several bioinformatics tools have been developed to detect CRISPR arrays based solely on DNA sequences, but all these tools employ the same strategy of looking for repetitive patterns, which might correspond to CRISPR array repeats. The identified patterns are evaluated using a fixed, built-in scoring function, and arrays exceeding a cut-off value are reported. Here, we instead introduce a data-driven approach that uses machine learning to detect and differentiate true CRISPR arrays from false ones based on several features. Our CRISPR detection tool, CRISPRidentify, performs three steps: detection, feature extraction and classification based on manually curated sets of positive and negative examples of CRISPR arrays. The identified CRISPR arrays are then reported to the user accompanied by detailed annotation. We demonstrate that our approach identifies not only previously detected CRISPR arrays, but also CRISPR array candidates not detected by other tools. Compared to other methods, our tool has a drastically reduced false positive rate. In contrast to the existing tools, our approach not only provides the user with the basic statistics on the identified CRISPR arrays but also produces a certainty score as a practical measure of the likelihood that a given genomic region is a CRISPR array.

Misregulation of BCL2 expression has been observed with many diseases and is associated with cellular exposure to reactive oxygen species. A region upstream of the P1 promoter in the human BCL2 gene plays a major role in regulating transcription. This G/C-rich region is highly polymorphic and capable of forming G-quadruplex structures. Herein we report that an oxidative event simulated with an 8-oxo-7,8-dihydroguanine (oxoG) substitution within a long G-tract results in a reduction of structural polymorphism. Surprisingly, oxoG within a 25-nt construct boosts thermal stability of the resulting G-quadruplex. This is achieved by distinct hydrogen bonding properties of oxoG, which facilitate formation of an antiparallel basket-type G-quadruplex with a three G-quartet core and a G·oxoG·C base triad. While oxoG has previously been considered detrimental for G-quadruplex formation, its stabilizing effect within a promoter described in this study suggests a potential novel regulatory role of oxidative stress in general and specifically in BCL2 gene transcription.

G-quadruplex DNA structures have become attractive drug targets, and native mass spectrometry can provide detailed characterization of drug binding stoichiometry and affinity, potentially at high throughput. However, the G-quadruplex DNA polymorphism poses problems for interpreting ligand screening assays. In order to establish standardized MS-based screening assays, we studied 28 sequences with documented NMR structures in (usually ∼100 mM) potassium, and report here their circular dichroism (CD), melting temperature (Tm), NMR spectra and electrospray mass spectra in 1 mM KCl/100 mM trimethylammonium acetate. Based on these results, we make a short-list of sequences that adopt the same structure in the MS assay as reported by NMR, and provide recommendations on using them for MS-based assays. We also built an R-based open-source application to build and consult a database, wherein further sequences can be incorporated in the future. The application handles automatically most of the data processing, and allows generating custom figures and reports. The database is included in the g4dbr package (https://github.com/EricLarG4/g4dbr) and can be explored online (https://ericlarg4.github.io/G4_database.html).

CRISPR-based base editors (BEs) are widely used to induce nucleotide substitutions in living cells and organisms without causing the damaging DNA double-strand breaks and DNA donor templates. Cytosine BEs that induce C:G to T:A conversion and adenine BEs that induce A:T to G:C conversion have been developed. Various attempts have been made to increase the efficiency of both BEs; however, their activities need to be improved for further applications. Here, we describe a fluorescent reporter-based drug screening platform to identify novel chemicals with the goal of improving adenine base editing efficiency. The reporter system revealed that histone deacetylase inhibitors, particularly romidepsin, enhanced base editing efficiencies by up to 4.9-fold by increasing the expression levels of proteins and target accessibility. The results support the use of romidepsin as a viable option to improve base editing efficiency in biomedical research and therapeutic genome engineering.

Acquisition of foreign DNA by Staphylococcus aureus, including vancomycin resistance genes, is thwarted by the ATP-dependent endonuclease SauUSI. Deciphering the mechanism of action of SauUSI could unravel the reason how it singularly plays a major role in preventing horizontal gene transfer (HGT) in S. aureus. Here, we report a detailed biochemical and structural characterization of SauUSI, which reveals that in the presence of ATP, the enzyme can cleave DNA having a single or multiple target site/s. Remarkably, in the case of multiple target sites, the entire region of DNA flanked by two target sites is shred into smaller fragments by SauUSI. Crystal structure of SauUSI reveals a stable dimer held together by the nuclease domains, which are spatially arranged to hydrolyze the phosphodiester bonds of both strands of the duplex. Thus, the architecture of the dimeric SauUSI facilitates cleavage of either single-site or multi-site DNA. The structure also provides insights into the molecular basis of target recognition by SauUSI. We show that target recognition activates ATP hydrolysis by the helicase-like ATPase domain, which powers active directional movement (translocation) of SauUSI along the DNA. We propose that a pile-up of multiple translocating SauUSI molecules against a stationary SauUSI bound to a target site catalyzes random double-stranded breaks causing shredding of the DNA between two target sites. The extensive and irreparable damage of the foreign DNA by shredding makes SauUSI a potent barrier against HGT.

The tRNA modification N6-threonylcarbamoyladenosine (t6A) is universally conserved in all organisms. In bacteria, the biosynthesis of t6A requires four proteins (TsaBCDE) that catalyze the formation of t6A via the unstable intermediate l-threonylcarbamoyl-adenylate (TC-AMP). While the formation and stability of this intermediate has been studied in detail, the mechanism of its transfer to A37 in tRNA is poorly understood. To investigate this step, the structure of the TsaBD heterodimer from Escherichia coli has been solved bound to a stable phosphonate isosteric mimic of TC-AMP. The phosphonate inhibits t6A synthesis in vitro with an IC50 value of 1.3 μM in the presence of millimolar ATP and L-threonine. The inhibitor binds to TsaBD by coordination to the active site Zn atom via an oxygen atom from both the phosphonate and the carboxylate moieties. The bound conformation of the inhibitor suggests that the catalysis exploits a putative oxyanion hole created by a conserved active site loop of TsaD and that the metal essentially serves as a binding scaffold for the intermediate. The phosphonate bound crystal structure should be useful for the rational design of potent, drug-like small molecule inhibitors as mechanistic probes or potentially novel antibiotics.

Biomarker-driven trials hold promise for therapeutic development in chronic diseases, such as muscular dystrophy. Myotonic dystrophy type 1 (DM1) involves RNA toxicity, where transcripts containing expanded CUG-repeats (CUGexp) accumulate in nuclear foci and sequester splicing factors in the Muscleblind-like (Mbnl) family. Oligonucleotide therapies to mitigate RNA toxicity have emerged but reliable measures of target engagement are needed. Here we examined muscle transcriptomes in mouse models of DM1 and found that CUGexp expression or Mbnl gene deletion cause similar dysregulation of alternative splicing. We selected 35 dysregulated exons for further study by targeted RNA sequencing. Across a spectrum of mouse models, the individual splice events and a composite index derived from all events showed a graded response to decrements of Mbnl or increments of CUGexp. Antisense oligonucleotides caused prompt reduction of CUGexp RNA and parallel correction of the splicing index, followed by subsequent elimination of myotonia. These results suggest that targeted splice sequencing may provide a sensitive and reliable way to assess therapeutic impact in DM1.

Methylglyoxal (MG) is a byproduct of glycolysis that functions in diverse mammalian developmental processes and diseases and in plant responses to various stresses, including salt stress. However, it is unknown whether MG-regulated gene expression is associated with an epigenetic modification. Here we report that MG methylglyoxalates H3 including H3K4 and increases chromatin accessibility, consistent with the result that H3 methylglyoxalation positively correlates with gene expression. Salt stress also increases H3 methylglyoxalation at salt stress responsive genes correlated to their higher expression. Following exposure to salt stress, salt stress responsive genes were expressed at higher levels in the Arabidopsis glyI2 mutant than in wild-type plants, but at lower levels in 35S::GLYI2 35S::GLYII4 plants, consistent with the higher and lower MG accumulation and H3 methylglyoxalation of target genes in glyI2 and 35S::GLYI2 35S::GLYII4, respectively. Further, ABI3 and MYC2, regulators of salt stress responsive genes, affect the distribution of H3 methylglyoxalation at salt stress responsive genes. Thus, MG functions as a histone-modifying group associated with gene expression that links glucose metabolism and epigenetic regulation.

Erratum in Nucleic Acids Res. 2021 Sep 7;49(15):9003.

Hemin [Fe(III)-protoporphyrin IX] is known to bind tightly to single-stranded DNA and RNA molecules that fold into G-quadruplexes (GQ). Such complexes are strongly activated for oxidative catalysis. These heme•DNAzymes and ribozymes have found broad utility in bioanalytical and medicinal chemistry and have also been shown to occur within living cells. However, how a GQ is able to activate hemin is poorly understood. Herein, we report fast kinetic measurements (using stopped-flow UV-vis spectrophotometry) to identify the H2O2-generated activated heme species within a heme•DNAzyme that is active for the oxidation of a thioether substrate, dibenzothiophene (DBT). Singular value decomposition and global fitting analysis was used to analyze the kinetic data, with the results being consistent with the heme•DNAzyme's DBT oxidation being catalyzed by the initial Fe(III)heme-H2O2 complex. Such a complex has been predicted computationally to be a powerful oxidant for thioether substrates. In the heme•DNAzyme, the DNA GQ enhances both the kinetics of formation of the active intermediate as well as the oxidation step of DBT by the active intermediate. We show, using both stopped flow spectrophotometry and EPR measurements, that a classic Compound I is not observable during the catalytic cycle for thioether sulfoxidation.

Multiple gene activities control complex biological processes such as cell fate specification during development and cellular reprogramming. Investigating the manifold gene functions in biological systems requires also simultaneous depletion of two or more gene activities. RNA interference-mediated knockdown (RNAi) is commonly used in Caenorhabditis elegans to assess essential genes, which otherwise lead to lethality or developmental arrest upon full knockout. RNAi application is straightforward by feeding worms with RNAi plasmid-containing bacteria. However, the general approach of mixing bacterial RNAi clones to deplete two genes simultaneously often yields poor results. To address this issue, we developed a bacterial conjugation-mediated double RNAi technique 'CONJUDOR'. It allows combining RNAi bacteria for robust double RNAi with high-throughput. To demonstrate the power of CONJUDOR for large scale double RNAi screens we conjugated RNAi against the histone chaperone gene lin-53 with more than 700 other chromatin factor genes. Thereby, we identified the Set1/MLL methyltransferase complex member RBBP-5 as a novel germ cell reprogramming barrier. Our findings demonstrate that CONJUDOR increases efficiency and versatility of RNAi screens to examine interconnected biological processes in C. elegans with high-throughput.

RNA structure formation in vivo happens co-transcriptionally while the transcript is being made. The corresponding co-transcriptional folding pathway typically involves transient RNA structure features that are not part of the final, functional RNA structure. These transient features can play important functional roles of their own and also influence the formation of the final RNA structure in vivo. We here present CoBold, a computational method for identifying different functional classes of transient RNA structure features that can either aid or hinder the formation of a known reference RNA structure. Our method takes as input either a single RNA or a corresponding multiple-sequence alignment as well as a known reference RNA secondary structure and identifies different classes of transient RNA structure features that could aid or prevent the formation of the given RNA structure. We make CoBold available via a web-server which includes dedicated data visualisation.

Operations with nucleic acids are among the main means of studying the mechanisms of gene function and developing novel methods of molecular medicine and gene therapy. These endeavours usually imply the necessity of nucleic acid storage and delivery into eukaryotic cells. In spite of diversity of the existing dedicated techniques, all of them have their limitations. Thus, a recent notion of using ionic liquids in manipulations of nucleic acids has been attracting significant attention lately. Due to their unique physicochemical properties, in particular, their micro-structuring impact and tunability, ionic liquids are currently applied as solvents and stabilizing media in chemical synthesis, electrochemistry, biotechnology, and other areas. Here, we review the current knowledge on interactions between nucleic acids and ionic liquids and discuss potential advantages of applying the latter in delivery of the former into eukaryotic cells.

DNA double-strand breaks (DSBs) at ribosomal gene loci trigger inhibition of ribosomal DNA (rDNA) transcription and extensive nucleolar reorganization, including the formation of nucleolar caps where rDNA DSBs engage with canonical DSB signaling and repair factors. While these nucleolar responses underlie maintenance of rDNA stability, the molecular components that drive each of these events remain to be defined. Here we report that full suppression of rRNA synthesis requires the DYRK1B kinase, a nucleolar DSB response that can be uncoupled from ATM-mediated DSB signaling events at the nucleolar periphery. Indeed, by targeting DSBs onto rDNA arrays, we uncovered that chemical inhibition or genetic inactivation of DYRK1B led to sustained nucleolar transcription. Not only does DYRK1B exhibit robust nucleolar accumulation following laser micro-irradiation across cell nuclei, we further showed that DYRK1B is required for rDNA DSB repair and rDNA copy number maintenance, and that DYRK1B-inactivated cells are hypersensitised to DSBs induced at the rDNA arrays. Together, our findings not only identify DYRK1B as a key signaling intermediate that coordinates DSB repair and rDNA transcriptional activities, but also support the idea of specialised DSB responses that operate within the nucleolus to preserve rDNA integrity.

G-quadruplexes (G4s) are four-stranded, guanine-rich nucleic acid structures that can influence a variety of biological processes such as the transcription and translation of genes and DNA replication. In many cases, a single G4-forming nucleic acid sequence can adopt multiple different folded conformations that interconvert on biologically relevant timescales, entropically stabilizing the folded state. The coexistence of different folded conformations also suggests that there are multiple pathways leading from the unfolded to the folded state ensembles, potentially modulating the folding rate and biological activity. We have developed an experimental method for quantifying the contributions of individual pathways to the folding of conformationally heterogeneous G4s that is based on mutagenesis, thermal hysteresis kinetic experiments and global analysis, and validated our results using photocaged kinetic NMR experiments. We studied the regulatory Pu22 G4 from the c-myc oncogene promoter, which adopts at least four distinct folded isomers. We found that the presence of four parallel pathways leads to a 2.5-fold acceleration in folding; that is, the effective folding rate from the unfolded to folded ensembles is 2.5 times as large as the rate constant for the fastest individual pathway. Since many G4 sequences can adopt many more than four isomers, folding accelerations of more than an order of magnitude are possible via this mechanism.

Thymine DNA glycosylase (TDG), as a repair enzyme, plays essential roles in maintaining the genome integrity by correcting several mismatched/damaged nucleobases. TDG acquires an efficient strategy to search for the lesions among a vast number of cognate base pairs. Currently, atomic-level details of how TDG translocates along DNA as it approaches the lesion site and the molecular mechanisms of the interplay between TDG and DNA are still elusive. Here, by constructing the Markov state model based on hundreds of molecular dynamics simulations with an integrated simulation time of ∼25 μs, we reveal the rotation-coupled sliding dynamics of TDG along a 9 bp DNA segment containing one G·T mispair. We find that TDG translocates along DNA at a relatively faster rate when distant from the lesion site, but slows down as it approaches the target, accompanied by deeply penetrating into the minor-groove, opening up the mismatched base pair and significantly sculpturing the DNA shape. Moreover, the electrostatic interactions between TDG and DNA are found to be critical for mediating the TDG translocation. Notably, several uncharacterized TDG residues are identified to take part in regulating the conformational switches of TDG occurred in the site-transfer process, which warrants further experimental validations.

DNA damage response is a fundamental mechanism to maintain genome stability. The ATR-WEE1 kinase module plays a central role in response to replication stress. Although the ATR-WEE1 pathway has been well studied in yeasts and animals, how ATR-WEE1 functions in plants remains unclear. Through a genetic screen for suppressors of the Arabidopsis atr mutant, we found that loss of function of PRL1, a core subunit of the evolutionarily conserved MAC complex involved in alternative splicing, suppresses the hypersensitivity of atr and wee1 to replication stress. Biochemical studies revealed that WEE1 directly interacts with and phosphorylates PRL1 at Serine 145, which promotes PRL1 ubiquitination and subsequent degradation. In line with the genetic and biochemical data, replication stress induces intron retention of cell cycle genes including CYCD1;1 and CYCD3;1, which is abolished in wee1 but restored in wee1 prl1. Remarkably, co-expressing the coding sequences of CYCD1;1 and CYCD3;1 partially restores the root length and HU response in wee1 prl1. These data suggested that the ATR-WEE1 module inhibits the MAC complex to regulate replication stress responses. Our study discovered PRL1 or the MAC complex as a key downstream regulator of the ATR-WEE1 module and revealed a novel cell cycle control mechanism.

Human DNA ligase I (LIG1) is the main replicative ligase and it also seals DNA breaks to complete DNA repair and recombination pathways. Immune compromised patients harbor hypomorphic LIG1 alleles encoding substitutions of conserved arginine residues, R771W and R641L, that compromise LIG1 activity through poorly defined mechanisms. To understand the molecular basis of LIG1 syndrome mutations, we determined high resolution X-ray structures and performed systematic biochemical characterization of LIG1 mutants using steady-state and pre-steady state kinetic approaches. Our results unveil a cooperative network of plastic DNA-LIG1 interactions that connect DNA substrate engagement with productive binding of Mg2+ cofactors for catalysis. LIG1 syndrome mutations destabilize this network, compromising Mg2+ binding affinity, decreasing ligation efficiency, and leading to elevated abortive ligation that may underlie the disease pathology. These findings provide novel insights into the fundamental mechanism by which DNA ligases engage with a nicked DNA substrate, and they suggest that disease pathology of LIG1 syndrome could be modulated by Mg2+ levels.

Mammalian circRNAs can influence different cellular processes by interacting with proteins and other nucleic acids. Here, we used ribonucleoprotein immunoprecipitation (RIP) analysis to identify systematically the circRNAs associated with the cancer-related protein AUF1. Among the circRNAs interacting with AUF1 in HeLa (human cervical carcinoma) cells, we focused on hsa_circ_0032434 (circPCNX), an abundant target of AUF1. Overexpression of circPCNX specifically interfered with the binding of AUF1 to p21 (CDKN1A) mRNA, thereby promoting p21 mRNA stability and elevating the production of p21, a major inhibitor of cell proliferation. Conversely, silencing circPCNX increased AUF1 binding to p21 mRNA, reducing p21 production and promoting cell division. Importantly, eliminating the AUF1-binding region of circPCNX abrogated the rise in p21 levels and rescued proliferation. Therefore, we propose that the interaction of circPCNX with AUF1 selectively prevents AUF1 binding to p21 mRNA, leading to enhanced p21 mRNA stability and p21 protein production, thereby suppressing cell growth.

Sequence-specific protein-DNA interactions are at the heart of the response of the tumor-suppressor p53 to numerous physiological and stress-related signals. Large variability has been previously reported in p53 binding to and transactivating from p53 response elements (REs) due, at least in part, to changes in direct (base) and indirect (shape) readouts of p53 REs. Here, we dissect p53 REs to decipher the mechanism by which p53 optimizes this highly regulated variable level of interaction with its DNA binding sites. We show that hemi-specific binding is more prevalent in p53 REs than previously envisioned. We reveal that sequences flanking the REs modulate p53 binding and activity and show that these effects extend to 4-5 bp from the REs. Moreover, we show here that the arrangement of p53 half-sites within its REs, relative to transcription direction, has been fine-tuned by selection pressure to optimize and regulate the response levels from p53 REs. This directionality in the REs arrangement is at least partly encoded in the structural properties of the REs. Furthermore, we show here that in the p21-5' RE the orientation of the half-sites is such that the effect of the flanking sequences is minimized and we discuss its advantages.

Escherichia coli SSB (EcSSB) is a model single-stranded DNA (ssDNA) binding protein critical in genome maintenance. EcSSB forms homotetramers that wrap ssDNA in multiple conformations to facilitate DNA replication and repair. Here we measure the binding and wrapping of many EcSSB proteins to a single long ssDNA substrate held at fixed tensions. We show EcSSB binds in a biphasic manner, where initial wrapping events are followed by unwrapping events as ssDNA-bound protein density passes critical saturation and high free protein concentration increases the fraction of EcSSBs in less-wrapped conformations. By destabilizing EcSSB wrapping through increased substrate tension, decreased substrate length, and protein mutation, we also directly observe an unstable bound but unwrapped state in which ∼8 nucleotides of ssDNA are bound by a single domain, which could act as a transition state through which rapid reorganization of the EcSSB-ssDNA complex occurs. When ssDNA is over-saturated, stimulated dissociation rapidly removes excess EcSSB, leaving an array of stably-wrapped complexes. These results provide a mechanism through which otherwise stably bound and wrapped EcSSB tetramers are rapidly removed from ssDNA to allow for DNA maintenance and replication functions, while still fully protecting ssDNA over a wide range of protein concentrations.

DNA gyrase, a type II topoisomerase found predominantly in bacteria, is the target for a variety of 'poisons', namely natural product toxins (e.g. albicidin, microcin B17) and clinically important synthetic molecules (e.g. fluoroquinolones). Resistance to both groups can be mediated by pentapeptide repeat proteins (PRPs). Despite long-term studies, the mechanism of action of these protective PRPs is not known. We show that a PRP, QnrB1 provides specific protection against fluoroquinolones, which strictly requires ATP hydrolysis by gyrase. QnrB1 binds to the GyrB protein and stimulates ATPase activity of the isolated N-terminal ATPase domain of GyrB (GyrB43). We probed the QnrB1 binding site using site-specific incorporation of a photoreactive amino acid and mapped the crosslinks to the GyrB43 protein. We propose a model in which QnrB1 binding allosterically promotes dissociation of the fluoroquinolone molecule from the cleavage complex.

Erratum in Nucleic Acids Res. 2021 Apr 19;49(7):4198.

The construction of microbial cell factories for sustainable production of chemicals and pharmaceuticals requires extensive genome engineering. Using Saccharomyces cerevisiae, this study proposes synthetic neochromosomes as orthogonal expression platforms for rewiring native cellular processes and implementing new functionalities. Capitalizing the powerful homologous recombination capability of S. cerevisiae, modular neochromosomes of 50 and 100 kb were fully assembled de novo from up to 44 transcriptional-unit-sized fragments in a single transformation. These assemblies were remarkably efficient and faithful to their in silico design. Neochromosomes made of non-coding DNA were stably replicated and segregated irrespective of their size without affecting the physiology of their host. These non-coding neochromosomes were successfully used as landing pad and as exclusive expression platform for the essential glycolytic pathway. This work pushes the limit of DNA assembly in S. cerevisiae and paves the way for de novo designer chromosomes as modular genome engineering platforms in S. cerevisiae.

Chromatin immunoprecipitation (IP) followed by sequencing (ChIP-seq) is the gold standard to detect transcription-factor (TF) binding sites in the genome. Its success depends on appropriate controls removing systematic biases. The predominantly used controls, i.e. DNA input, correct for uneven sonication, but not for nonspecific interactions of the IP antibody. Another type of controls, 'mock' IP, corrects for both of the issues, but is not widely used because it is considered susceptible to technical noise. The tradeoff between the two control types has not been investigated systematically. Therefore, we generated comparable DNA input and mock IP experiments. Because mock IPs contain only nonspecific interactions, the sites predicted from them using DNA input indicate the spurious-site abundance. This abundance is highly correlated with the 'genomic activity' (e.g. chromatin openness). In particular, compared to cell lines, complex samples such as whole organisms have more spurious sites-probably because they contain multiple cell types, resulting in more expressed genes and more open chromatin. Consequently, DNA input and mock IP controls performed similarly for cell lines, whereas for complex samples, mock IP substantially reduced the number of spurious sites. However, DNA input is still informative; thus, we developed a simple framework integrating both controls, improving binding site detection.

Single-cell RNA sequencing enables us to characterize the cellular heterogeneity in single cell resolution with the help of cell type identification algorithms. However, the noise inherent in single-cell RNA-sequencing data severely disturbs the accuracy of cell clustering, marker identification and visualization. We propose that clustering based on feature density profiles can distinguish informative features from noise. We named such strategy as 'entropy subspace' separation and designed a cell clustering algorithm called ENtropy subspace separation-based Clustering for nOise REduction (ENCORE) by integrating the 'entropy subspace' separation strategy with a consensus clustering method. We demonstrate that ENCORE performs superiorly on cell clustering and generates high-resolution visualization across 12 standard datasets. More importantly, ENCORE enables identification of group markers with biological significance from a hard-to-separate dataset. With the advantages of effective feature selection, improved clustering, accurate marker identification and high-resolution visualization, we present ENCORE to the community as an important tool for scRNA-seq data analysis to study cellular heterogeneity and discover group markers.

The gene and cell therapy fields are advancing rapidly, with a potential to treat and cure a wide range of diseases, and lentivirus-based gene transfer agents are the vector of choice for many investigators. Early cases of insertional mutagenesis caused by gammaretroviral vectors highlighted that integration site (IS) analysis was a major safety and quality control checkpoint for lentiviral applications. The methods established to detect lentiviral integrations using next-generation sequencing (NGS) are limited by short read length, inadvertent PCR bias, low yield, or lengthy protocols. Here, we describe a new method to sequence IS using Amplification-free Integration Site sequencing (AFIS-Seq). AFIS-Seq is based on amplification-free, Cas9-mediated enrichment of high-molecular-weight chromosomal DNA suitable for long-range Nanopore MinION sequencing. This accessible and low-cost approach generates long reads enabling IS mapping with high certainty within a single day. We demonstrate proof-of-concept by mapping IS of lentiviral vectors in a variety of cell models and report up to 1600-fold enrichment of the signal. This method can be further extended to sequencing of Cas9-mediated integration of genes and to in vivo analysis of IS. AFIS-Seq uses long-read sequencing to facilitate safety evaluation of preclinical lentiviral vector gene therapies by providing IS analysis with improved confidence.

Single mRNA molecules are frequently detected by single molecule fluorescence in situ hybridization (smFISH) using branched DNA technology. While providing strong and background-reduced signals, the method is inefficient in detecting mRNAs within dense structures, in monitoring mRNA compactness and in quantifying abundant mRNAs. To overcome these limitations, we have hybridized slices of high pressure frozen, freeze-substituted and LR White embedded cells (LR White smFISH). mRNA detection is physically restricted to the surface of the resin. This enables single molecule detection of RNAs with accuracy comparable to RNA sequencing, irrespective of their abundance, while at the same time providing spatial information on RNA localization that can be complemented with immunofluorescence and electron microscopy, as well as array tomography. Moreover, LR White embedding restricts the number of available probe pair recognition sites for each mRNA to a small subset. As a consequence, differences in signal intensities between RNA populations reflect differences in RNA structures, and we show that the method can be employed to determine mRNA compactness. We apply the method to answer some outstanding questions related to trans-splicing, RNA granules and mitochondrial RNA editing in single-cellular trypanosomes and we show an example of differential gene expression in the metazoan Caenorhabditis elegans.

The occurrence of accidental mutations or deletions caused by genome editing with CRISPR/Cas9 system remains a critical unsolved problem of the technology. Blocking excess or prolonged Cas9 activity in cells is considered as one means of solving this problem. Here, we report the development of an inhibitory DNA aptamer against Cas9 by means of in vitro selection (systematic evolution of ligands by exponential enrichment) and subsequent screening with an in vitro cleavage assay. The inhibitory aptamer could bind to Cas9 at low nanomolar affinity and partially form a duplex with CRISPR RNA, contributing to its inhibitory activity. We also demonstrated that improving the inhibitory aptamer with locked nucleic acids efficiently suppressed Cas9-directed genome editing in cells and reduced off-target genome editing. The findings presented here might enable the development of safer and controllable genome editing for biomedical research and gene therapy.

Regulatory nascent peptides participate in the regulation of cellular functions by the mechanisms involving regulated translation arrest. A class of them in bacteria, called monitoring substrates, feedback-regulates the expression of a specific component of protein localization machinery. Three monitoring substrates, SecM, MifM and VemP have previously been identified. Here, we attempt at identifying additional arrest peptides in bacteria. Our bioinformatic searches over more than 400 bacterial genomic sequences for proteins that have the common characteristic features shared by the known monitoring substrates and subsequent in vitro and in vivo characterization of the highlighted sequences allowed the identification of three arrest peptides termed ApcA, ApdA and ApdP. ApcA and ApdA homologs are conserved among a subset of actinobacteria, whereas ApdP has homologs in a subset of α-proteobacteria. We demonstrate that these arrest peptides, in their ribosome-tethered nascent states, inhibit peptidyl transfer. The elongation arrest occurs at a specific codon near the 3' end of the coding region, in a manner depending on the amino acid sequence of the nascent chain. Interestingly, the arrest sequences of ApcA, ApdA and ApdP share a sequence R-A-P-G/P that is essential for the elongation arrest.

G-quadruplex (G4) DNA structures with a left-handed backbone progression have unique and conserved structural features. Studies on sequence dependency of the structures revealed the prerequisites and some minimal motifs required for left-handed G4 formation. To extend the boundaries, we explore the adaptability of left-handed G4s towards the existence of bulges. Here we present two X-ray crystal structures and an NMR solution structure of left-handed G4s accommodating one, two and three bulges. Bulges in left-handed G4s show distinct characteristics as compared to those in right-handed G4s. The elucidation of intricate structural details will help in understanding the possible roles and limitations of these unique structures.

As an effective programmable DNA targeting tool, CRISPR-Cas9 system has been adopted in varieties of biotechnological applications. However, the off-target effects, derived from the tolerance towards guide-target mismatches, are regarded as the major problems in engineering CRISPR systems. To understand this, we constructed two sgRNA libraries carrying saturated single- and double-nucleotide mismatches in living bacteria cells, and profiled the comprehensive landscape of in vivo binding affinity of dCas9 toward DNA target guided by each individual sgRNA with particular mismatches. We observed a synergistic effect in seed, where combinatorial double mutations caused more severe activity loss compared with the two corresponding single mutations. Moreover, we found that a particular mismatch type, dDrG (D = A, T, G), only showed moderate impairment on binding. To quantitatively understand the causal relationship between mismatch and binding behaviour of dCas9, we further established a biophysical model, and found that the thermodynamic properties of base-pairing coupled with strand invasion process, to a large extent, can account for the observed mismatch-activity landscape. Finally, we repurposed this model, together with a convolutional neural network constructed based on the same mechanism, as a predictive tool to guide the rational design of sgRNA in bacterial CRISPR interference.

The Ebola virus is a deadly human pathogen responsible for several outbreaks in Africa. Its genome encodes the 'large' L protein, an essential enzyme that has polymerase, capping and methyltransferase activities. The methyltransferase activity leads to RNA co-transcriptional modifications at the N7 position of the cap structure and at the 2'-O position of the first transcribed nucleotide. Unlike other Mononegavirales viruses, the Ebola virus methyltransferase also catalyses 2'-O-methylation of adenosines located within the RNA sequences. Herein, we report the crystal structure at 1.8 Å resolution of the Ebola virus methyltransferase domain bound to a fragment of a camelid single-chain antibody. We identified structural determinants and key amino acids specifically involved in the internal adenosine-2'-O-methylation from cap-related methylations. These results provide the first high resolution structure of an ebolavirus L protein domain, and the framework to investigate the effects of epitranscriptomic modifications and to design possible antiviral drugs against the Filoviridae family.

We describe a synthetic riboswitch element that implements a regulatory principle which directly addresses an essential tRNA maturation step. Constructed using a rational in silico design approach, this riboswitch regulates RNase P-catalyzed tRNA 5'-processing by either sequestering or exposing the single-stranded 5'-leader region of the tRNA precursor in response to a ligand. A single base pair in the 5'-leader defines the regulatory potential of the riboswitch both in vitro and in vivo. Our data provide proof for prior postulates on the importance of the structure of the leader region for tRNA maturation. We demonstrate that computational predictions of ligand-dependent structural rearrangements can address individual maturation steps of stable non-coding RNAs, thus making them amenable as promising target for regulatory devices that can be used as functional building blocks in synthetic biology.

Homologous recombination dominates as the major form of DNA repair in Trypanosoma brucei, and is especially important for recombination of the subtelomeric variant surface glycoprotein during antigenic variation. RAD50, a component of the MRN complex (MRE11, RAD50, NBS1), is central to homologous recombination through facilitating resection and governing the DNA damage response. The function of RAD50 in trypanosomes is untested. Here we report that RAD50 and MRE11 are required for RAD51-dependent homologous recombination and phosphorylation of histone H2A following a DNA double strand break (DSB), but neither MRE11 nor RAD50 substantially influence DSB resection at a chromosome-internal locus. In addition, we reveal intrinsic separation-of-function between T. brucei RAD50 and MRE11, with only RAD50 suppressing DSB repair using donors with short stretches of homology at a subtelomeric locus, and only MRE11 directing DSB resection at the same locus. Finally, we show that loss of either MRE11 or RAD50 causes a greater diversity of expressed VSG variants following DSB repair. We conclude that MRN promotes stringent homologous recombination at subtelomeric loci and restrains antigenic variation.

Replication protein A (RPA) binds to single-stranded DNA (ssDNA) and interacts with over three dozen enzymes and serves as a recruitment hub to coordinate most DNA metabolic processes. RPA binds ssDNA utilizing multiple oligosaccharide/oligonucleotide binding domains and based on their individual DNA binding affinities are classified as high versus low-affinity DNA-binding domains (DBDs). However, recent evidence suggests that the DNA-binding dynamics of DBDs better define their roles. Utilizing hydrogen-deuterium exchange mass spectrometry (HDX-MS), we assessed the ssDNA-driven dynamics of the individual domains of human RPA. As expected, ssDNA binding shows HDX changes in DBDs A, B, C, D and E. However, DBD-A and DBD-B are dynamic and do not show robust DNA-dependent protection. DBD-C displays the most extensive changes in HDX, suggesting a major role in stabilizing RPA on ssDNA. Slower allosteric changes transpire in the protein-protein interaction domains and linker regions, and thus do not directly interact with ssDNA. Within a dynamics-based model for RPA, we propose that DBD-A and -B act as the dynamic half and DBD-C, -D and -E function as the less-dynamic half. Thus, segments of ssDNA buried under the dynamic half are likely more readily accessible to RPA-interacting proteins.

Pre-mRNA splicing catalyzed by the spliceosome represents a critical step in the regulation of gene expression contributing to transcriptome and proteome diversity. The spliceosome consists of five small nuclear ribonucleoprotein particles (snRNPs), the biogenesis of which remains only partially understood. Here we define the evolutionarily conserved protein Ecdysoneless (Ecd) as a critical regulator of U5 snRNP assembly and Prp8 stability. Combining Drosophila genetics with proteomic approaches, we demonstrate the Ecd requirement for the maintenance of adult healthspan and lifespan and identify the Sm ring protein SmD3 as a novel interaction partner of Ecd. We show that the predominant task of Ecd is to deliver Prp8 to the emerging U5 snRNPs in the cytoplasm. Ecd deficiency, on the other hand, leads to reduced Prp8 protein levels and compromised U5 snRNP biogenesis, causing loss of splicing fidelity and transcriptome integrity. Based on our findings, we propose that Ecd chaperones Prp8 to the forming U5 snRNP allowing completion of the cytoplasmic part of the U5 snRNP biogenesis pathway necessary to meet the cellular demand for functional spliceosomes.

Erratum in Nucleic Acids Res. 2021 Apr 6;49(6):3600.

The RNA-dependent RNA polymerases (RdRPs) encoded by RNA viruses represent a unique class of nucleic acid polymerases. RdRPs are essential in virus life cycle due to their central role in viral genome replication/transcription processes. However, their contribution in host adaption has not been well documented. By solving the RdRP crystal structure of the tick-borne encephalitis virus (TBEV), a tick-borne flavivirus, and comparing the structural and sequence features with mosquito-borne flavivirus RdRPs, we found that a region between RdRP catalytic motifs B and C, namely region B-C, clearly bears host-related diversity. Inter-virus substitutions of region B-C sequence were designed in both TBEV and mosquito-borne Japanese encephalitis virus backbones. While region B-C substitutions only had little or moderate effect on RdRP catalytic activities, virus proliferation was not supported by these substitutions in both virus systems. Importantly, a TBEV replicon-derived viral RNA replication was significantly reduced but not abolished by the substitution, suggesting the involvement of region B-C in viral and/or host processes beyond RdRP catalysis. A systematic structural analysis of region B-C in viral RdRPs further emphasizes its high level of structure and length diversity, providing a basis to further refine its relevance in RNA virus-host interactions in a general context.

Classification of high-throughput genomic data is a powerful method to assign samples to subgroups with specific molecular profiles. Consensus partitioning is the most widely applied approach to reveal subgroups by summarizing a consensus classification from a list of individual classifications generated by repeatedly executing clustering on random subsets of the data. It is able to evaluate the stability of the classification. We implemented a new R/Bioconductor package, cola, that provides a general framework for consensus partitioning. With cola, various parameters and methods can be user-defined and easily integrated into different steps of an analysis, e.g., feature selection, sample classification or defining signatures. cola provides a new method named ATC (ability to correlate to other rows) to extract features and recommends spherical k-means clustering (skmeans) for subgroup classification. We show that ATC and skmeans have better performance than other commonly used methods by a comprehensive benchmark on public datasets. We also benchmark key parameters in the consensus partitioning procedure, which helps users to select optimal parameter values. Moreover, cola provides rich functionalities to apply multiple partitioning methods in parallel and directly compare their results, as well as rich visualizations. cola can automate the complete analysis and generates a comprehensive HTML report.

Erratum in Nucleic Acids Res. 2021 Feb 22;49(3):1801.

Erratum in Nucleic Acids Res. 2021 Apr 19;49(7):4197.

Many modification-dependent restriction endonucleases (MDREs) are fusions of a PUA superfamily modification sensor domain and a nuclease catalytic domain. EVE domains belong to the PUA superfamily, and are present in MDREs in combination with HNH nuclease domains. Here, we present a biochemical characterization of the EVE-HNH endonuclease VcaM4I and crystal structures of the protein alone, with EVE domain bound to either 5mC modified dsDNA or to 5mC/5hmC containing ssDNA. The EVE domain is moderately specific for 5mC/5hmC containing DNA according to EMSA experiments. It flips the modified nucleotide, to accommodate it in a hydrophobic pocket of the enzyme, primarily formed by P24, W82 and Y130 residues. In the crystallized conformation, the EVE domain and linker helix between the two domains block DNA binding to the catalytic domain. Removal of the EVE domain and inter-domain linker, but not of the EVE domain alone converts VcaM4I into a non-specific toxic nuclease. The role of the key residues in the EVE and HNH domains of VcaM4I is confirmed by digestion and restriction assays with the enzyme variants that differ from the wild-type by changes to the base binding pocket or to the catalytic residues.

Alkylating drugs are among the most often used chemotherapeutics. While cancer cells frequently develop resistance to alkylation treatments, detailed understanding of mechanisms that lead to the resistance is limited. Here, by using genome-wide CRISPR-Cas9 based screen, we identify transcriptional Mediator complex subunit 13 (MED13) as a novel modulator of alkylation response. The alkylation exposure causes significant MED13 downregulation, while complete loss of MED13 results in reduced apoptosis and resistance to alkylating agents. Transcriptome analysis identified cyclin D1 (CCND1) as one of the highly overexpressed genes in MED13 knock-out (KO) cells, characterized by shorter G1 phase. MED13 is able to bind to CCND1 regulatory elements thus influencing the expression. The resistance of MED13 KO cells is directly dependent on the cyclin D1 overexpression, and its down-regulation is sufficient to re-sensitize the cells to alkylating agents. We further demonstrate the therapeutic potential of MED13-mediated response, by applying combinatory treatment with CDK8/19 inhibitor Senexin A. Importantly, the treatment with Senexin A stabilizes MED13, and in combination with alkylating agents significantly reduces viability of cancer cells. In summary, our findings identify novel alkylation stress response mechanism dependent on MED13 and cyclin D1 that can serve as basis for development of innovative therapeutic strategies.

DNA methylation is important for the epigenetic regulation of gene expression and plays a critical role in mammalian development. However, the dynamic regulation of genome-wide DNA methylation in skeletal muscle development remains largely unknown. Here, we generated the first single-base resolution DNA methylome and transcriptome maps of porcine skeletal muscle across 27 developmental stages. The overall methylation level decreased from the embryo to the adult, which was highly correlated with the downregulated expression of DNMT1 and an increase in partially methylated domains. Notably, we identified over 40 000 developmentally differentially methylated CpGs (dDMCs) that reconstitute the developmental trajectory of skeletal muscle and associate with muscle developmental genes and transcription factors (TFs). The dDMCs were significantly under-represented in promoter regulatory regions but strongly enriched as enhancer histone markers and in chromatin-accessible regions. Integrative analysis revealed the negative regulation of both promoter and gene body methylation in genes associated with muscle contraction and insulin signaling during skeletal muscle development. Mechanistically, DNA methylation affected the expression of muscle-related genes by modulating the accessibly of upstream myogenesis TF binding, indicating the involvement of the DNA methylation/SP1/IGF2BP3 axis in skeletal myogenesis. Our results highlight the function and regulation of dynamic DNA methylation in skeletal muscle development.

Underlying higher order chromatin organization are Structural Maintenance of Chromosomes (SMC) complexes, large protein rings that entrap DNA. The molecular mechanism by which SMC complexes organize chromatin is as yet incompletely understood. Two prominent models posit that SMC complexes actively extrude DNA loops (loop extrusion), or that they sequentially entrap two DNAs that come into proximity by Brownian motion (diffusion capture). To explore the implications of these two mechanisms, we perform biophysical simulations of a 3.76 Mb-long chromatin chain, the size of the long Schizosaccharomyces pombe chromosome I left arm. On it, the SMC complex condensin is modeled to perform loop extrusion or diffusion capture. We then compare computational to experimental observations of mitotic chromosome formation. Both loop extrusion and diffusion capture can result in native-like contact probability distributions. In addition, the diffusion capture model more readily recapitulates mitotic chromosome axis shortening and chromatin compaction. Diffusion capture can also explain why mitotic chromatin shows reduced, as well as more anisotropic, movements, features that lack support from loop extrusion. The condensin distribution within mitotic chromosomes, visualized by stochastic optical reconstruction microscopy (STORM), shows clustering predicted from diffusion capture. Our results inform the evaluation of current models of mitotic chromosome formation.

Sensing of environmental cues is crucial for cell survival. To adapt to changes in their surroundings cells need to tightly control the repertoire of genes expressed at any time. Regulation of translation is key, especially in organisms in which transcription is hardly controlled, like Trypanosoma brucei. In this study, we describe the shortening of the bulk of the cellular tRNAs during stress at the expense of the conserved 3' CCA-tail. This tRNA shortening is specific for nutritional stress and renders tRNAs unsuitable substrates for translation. We uncovered the nuclease LCCR4 (Tb927.4.2430), a homologue of the conserved deadenylase Ccr4, as being responsible for tRNA trimming. Once optimal growth conditions are restored tRNAs are rapidly repaired by the trypanosome tRNA nucleotidyltransferase thus rendering the recycled tRNAs amenable for translation. This mechanism represents a fast and efficient way to repress translation during stress, allowing quick reactivation with a low energy input.

Super-enhancers (SEs) mediate high transcription levels of target genes. Previous studies have shown that SEs recruit transcription complexes and generate enhancer RNAs (eRNAs). We characterized transcription at the human and murine β-globin locus control region (LCR) SE. We found that the human LCR is capable of recruiting transcription complexes independently from linked globin genes in transgenic mice. Furthermore, LCR hypersensitive site 2 (HS2) initiates the formation of bidirectional transcripts in transgenic mice and in the endogenous β-globin gene locus in murine erythroleukemia (MEL) cells. HS2 3'eRNA is relatively unstable and remains in close proximity to the globin gene locus. Reducing the abundance of HS2 3'eRNA leads to a reduction in β-globin gene transcription and compromises RNA polymerase II (Pol II) recruitment at the promoter. The Integrator complex has been shown to terminate eRNA transcription. We demonstrate that Integrator interacts downstream of LCR HS2. Inducible ablation of Integrator function in MEL or differentiating primary human CD34+ cells causes a decrease in expression of the adult β-globin gene and accumulation of Pol II and eRNA at the LCR. The data suggest that transcription complexes are assembled at the LCR and transferred to the globin genes by mechanisms that involve Integrator mediated release of Pol II and eRNA from the LCR.

In most bacteria, efficient use of carbohydrates is primarily mediated by the phosphoenolpyruvate (PEP):carbohydrate phosphotransferase system (PTS), which concomitantly phosphorylates the substrates during import. Therefore, transcription of the PTS-encoding genes is precisely regulated by transcriptional regulators, depending on the availability of the substrate. Fructose is transported mainly through the fructose-specific PTS (PTSFru) and simultaneously converted into fructose 1-phosphate (F1P). In Gammaproteobacteria such as Escherichia coli and Pseudomonas putida, transcription of the fru operon encoding two PTSFru components, FruA and FruB, and the 1-phosphofructokinase FruK is repressed by FruR in the absence of the inducer F1P. Here, we show that, contrary to the case in other Gammaproteobacteria, FruR acts as a transcriptional activator of the fru operon and is indispensable for the growth of Vibrio cholerae on fructose. Several lines of evidence suggest that binding of the FruR-F1P complex to an operator which is located between the -35 and -10 promoter elements changes the DNA structure to facilitate RNA polymerase binding to the promoter. We discuss the mechanism by which the highly conserved FruR regulates the expression of its target operon encoding the highly conserved PTSFru and FruK in a completely opposite direction among closely related families of bacteria.

The endonuclease activity within the influenza virus cap-snatching process is a proven therapeutic target. The anti-influenza drug baloxavir is highly effective, but is associated with resistance mutations that threaten its clinical efficacy. The endonuclease resides within the N-terminal domain of the PA subunit (PAN) of the influenza RNA dependent RNA polymerase, and we report here complexes of PAN with RNA and DNA oligonucleotides to understand its specificity and the structural basis of baloxavir resistance mutations. The RNA and DNA oligonucleotides bind within the substrate binding groove of PAN in a similar fashion, explaining the ability of the enzyme to cleave both substrates. The individual nucleotides occupy adjacent conserved pockets that flank the two-metal active site. However, the 2' OH of the RNA ribose moieties engage in additional interactions that appear to optimize the binding and cleavage efficiency for the natural substrate. The major baloxavir resistance mutation at position 38 is at the core of the substrate binding site, but structural studies and modeling suggest that it maintains the necessary virus fitness via compensating interactions with RNA. These studies will facilitate the development of new influenza therapeutics that spatially match the substrate and are less likely to elicit resistance mutations.

Human telomeres contain the repeat DNA sequence 5'-d(TTAGGG), with duplex regions that are several kilobases long terminating in a 3' single-stranded overhang. The structure of the single-stranded overhang is not known with certainty, with disparate models proposed in the literature. We report here the results of an integrated structural biology approach that combines small-angle X-ray scattering, circular dichroism (CD), analytical ultracentrifugation, size-exclusion column chromatography and molecular dynamics simulations that provide the most detailed characterization to date of the structure of the telomeric overhang. We find that the single-stranded sequences 5'-d(TTAGGG)n, with n = 8, 12 and 16, fold into multimeric structures containing the maximal number (2, 3 and 4, respectively) of contiguous G4 units with no long gaps between units. The G4 units are a mixture of hybrid-1 and hybrid-2 conformers. In the multimeric structures, G4 units interact, at least transiently, at the interfaces between units to produce distinctive CD signatures. Global fitting of our hydrodynamic and scattering data to a worm-like chain (WLC) model indicates that these multimeric G4 structures are semi-flexible, with a persistence length of ∼34 Å. Investigations of its flexibility using MD simulations reveal stacking, unstacking, and coiling movements, which yield unique sites for drug targeting.

Approximately 13% of the human genome can fold into non-canonical (non-B) DNA structures (e.g. G-quadruplexes, Z-DNA, etc.), which have been implicated in vital cellular processes. Non-B DNA also hinders replication, increasing errors and facilitating mutagenesis, yet its contribution to genome-wide variation in mutation rates remains unexplored. Here, we conducted a comprehensive analysis of nucleotide substitution frequencies at non-B DNA loci within noncoding, non-repetitive genome regions, their ±2 kb flanking regions, and 1-Megabase windows, using human-orangutan divergence and human single-nucleotide polymorphisms. Functional data analysis at single-base resolution demonstrated that substitution frequencies are usually elevated at non-B DNA, with patterns specific to each non-B DNA type. Mirror, direct and inverted repeats have higher substitution frequencies in spacers than in repeat arms, whereas G-quadruplexes, particularly stable ones, have higher substitution frequencies in loops than in stems. Several non-B DNA types also affect substitution frequencies in their flanking regions. Finally, non-B DNA explains more variation than any other predictor in multiple regression models for diversity or divergence at 1-Megabase scale. Thus, non-B DNA substantially contributes to variation in substitution frequencies at small and large scales. Our results highlight the role of non-B DNA in germline mutagenesis with implications to evolution and genetic diseases.

The maternal mode of mitochondrial DNA (mtDNA) inheritance is central to human genetics. Recently, evidence for bi-parental inheritance of mtDNA was claimed for individuals of three pedigrees that suffered mitochondrial disorders. We sequenced mtDNA using both direct Sanger and Massively Parallel Sequencing in several tissues of eleven maternally related and other affiliated healthy individuals of a family pedigree and observed mixed mitotypes in eight individuals. Cells without nuclear DNA, i.e. thrombocytes and hair shafts, only showed the mitotype of haplogroup (hg) V. Skin biopsies were prepared to generate ρ° cells void of mtDNA, sequencing of which resulted in a hg U4c1 mitotype. The position of the Mega-NUMT sequence was determined by fluorescence in situ hybridization and two different quantitative PCR assays were used to determine the number of contributing mtDNA copies. Thus, evidence for the presence of repetitive, full mitogenome Mega-NUMTs matching haplogroup U4c1 in various tissues of eight maternally related individuals was provided. Multi-copy Mega-NUMTs mimic mixtures of mtDNA that cannot be experimentally avoided and thus may appear in diverse fields of mtDNA research and diagnostics. We demonstrate that hair shaft mtDNA sequencing provides a simple but reliable approach to exclude NUMTs as source of misleading results.

Ribosomes are intricate molecular machines ensuring proper protein synthesis in every cell. Ribosome biogenesis is a complex process which has been intensively analyzed in bacteria and eukaryotes. In contrast, our understanding of the in vivo archaeal ribosome biogenesis pathway remains less characterized. Here, we have analyzed the in vivo role of the almost universally conserved ribosomal RNA dimethyltransferase KsgA/Dim1 homolog in archaea. Our study reveals that KsgA/Dim1-dependent 16S rRNA dimethylation is dispensable for the cellular growth of phylogenetically distant archaea. However, proteomics and functional analyses suggest that archaeal KsgA/Dim1 and its rRNA modification activity (i) influence the expression of a subset of proteins and (ii) contribute to archaeal cellular fitness and adaptation. In addition, our study reveals an unexpected KsgA/Dim1-dependent variability of rRNA modifications within the archaeal phylum. Combining structure-based functional studies across evolutionary divergent organisms, we provide evidence on how rRNA structure sequence variability (re-)shapes the KsgA/Dim1-dependent rRNA modification status. Finally, our results suggest an uncoupling between the KsgA/Dim1-dependent rRNA modification completion and its release from the nascent small ribosomal subunit. Collectively, our study provides additional understandings into principles of molecular functional adaptation, and further evolutionary and mechanistic insights into an almost universally conserved step of ribosome synthesis.

We have identified chemical probes that simultaneously inhibit cancer cell progression and an immune checkpoint. Using the computational Site Identification by Ligand Competitive Saturation (SILCS) technology, structural biology and cell-based assays, we identify small molecules that directly and selectively bind to the RNA Recognition Motif (RRM) of hnRNP A18, a regulator of protein translation in cancer cells. hnRNP A18 recognizes a specific RNA signature motif in the 3'UTR of transcripts associated with cancer cell progression (Trx, VEGF, RPA) and, as shown here, a tumor immune checkpoint (CTLA-4). Post-transcriptional regulation of immune checkpoints is a potential therapeutic strategy that remains to be exploited. The probes target hnRNP A18 RRM in vitro and in cells as evaluated by cellular target engagement. As single agents, the probes specifically disrupt hnRNP A18-RNA interactions, downregulate Trx and CTLA-4 protein levels and inhibit proliferation of several cancer cell lines without affecting the viability of normal epithelial cells. These first-in-class chemical probes will greatly facilitate the elucidation of the underexplored biological function of RNA Binding Proteins (RBPs) in cancer cells, including their effects on proliferation and immune checkpoint activation.