Internal ribosome entry sites (IRESs) are RNA elements capable of initiating translation on an internal portion of a messenger RNA. The intergenic region (IGR) IRES of the Dicistroviridae virus family folds into a triple pseudoknot tertiary structure, allowing it to recruit the ribosome and initiate translation in a structure dependent manner. This IRES has also been reported to drive translation in Escherichia coli and to date is the only described translation initiation signal that functions across domains of life. Here we show that unlike in the eukaryotic context the tertiary structure of the IGR IRES is not required for prokaryotic ribosome recruitment. In E. coli IGR IRES translation efficiency is dependent on ribosomal protein S1 in conjunction with an AG-rich Shine-Dalgarno-like element, supporting a model where the translational activity of the IGR IRESs is due to S1-mediated canonical prokaryotic translation.

Reactive oxygen species attack the structure of DNA, thus altering its base-pairing properties. Consequently, oxidative stress-associated DNA lesions are a major source of the mutation load that gives rise to cancer and other diseases. Base excision repair (BER) is the pathway primarily tasked with repairing DNA base damage, with apurinic/apyrimidinic endonuclease (APE1) having both AP-endonuclease and 3' to 5' exonuclease (exo) DNA cleavage functions. The lesion 8-oxo-7,8-dihydroguanine (8-oxoG) can enter the genome as either a product of direct damage to the DNA, or through polymerase insertion at the 3'-end of a DNA strand during replication or repair. Importantly, 3'-8-oxoG impairs the ligation step of BER and therefore must be removed by the exo activity of a surrogate enzyme to prevent double stranded breaks and cell death. In the present study, we use X-ray crystallography to characterize the exo activity of APE1 on 3'-8-oxoG substrates. These structures support a unified APE1 exo mechanism that differs from its more canonical AP-endonuclease activity. In addition, through complementation of the structural data with enzyme kinetics and binding studies employing both wild-type and rationally designed APE1 mutants, we were able to identify and characterize unique protein: DNA contacts that specifically mediate 8-oxoG removal by APE1.

In eukaryotes, many stable and heritable phenotypes arise from the same DNA sequence, owing to epigenetic regulatory mechanisms relying on the molecular cooperativity of 'reader-writer' enzymes. In this work, we focus on the fundamental, generic mechanisms behind the epigenome memory encoded by post-translational modifications of histone tails. Based on experimental knowledge, we introduce a unified modeling framework, the painter model, describing the mechanistic interplay between sequence-specific recruitment of chromatin regulators, chromatin-state-specific reader-writer processes and long-range spreading mechanisms. A systematic analysis of the model building blocks highlights the crucial impact of tridimensional chromatin organization and state-specific recruitment of enzymes on the stability of epigenomic domains and on gene expression. In particular, we show that enhanced 3D compaction of the genome and enzyme limitation facilitate the formation of ultra-stable, confined chromatin domains. The model also captures how chromatin state dynamics impact the intrinsic transcriptional properties of the region, slower kinetics leading to noisier expression. We finally apply our framework to analyze experimental data, from the propagation of γH2AX around DNA breaks in human cells to the maintenance of heterochromatin in fission yeast, illustrating how the painter model can be used to extract quantitative information on epigenomic molecular processes.

Plasmids of the ColE1 family are among the most frequently used in molecular biology. They were adopted early for many biotechnology applications, and as models to study plasmid biology. Their mechanism of replication is well understood, involving specific interactions between a plasmid encoded sense-antisense gene pair (RNAI and RNAII). Due to such mechanism, two plasmids with the same origin cannot be stably maintained in cells-a process known as incompatibility. While mutations in RNAI and RNAII can make colE1 more compatible, there has been no systematic effort to engineer new compatible colE1 origins, which could bypass technical design constraints for multi-plasmid applications. Here, we show that by diversifying loop regions in RNAI (and RNAII), it is possible to select new viable colE1 origins compatible with the wild-type one. We demonstrate that sequence divergence is not sufficient to enable compatibility and pairwise interactions are not an accurate guide for higher order interactions. We identify potential principles to engineer plasmid copy number independently from other regulatory strategies and we propose plasmid compatibility as a tractable model to study biological orthogonality.

The Smc5/6 complex plays an essential role in the resolution of recombination intermediates formed during mitosis or meiosis, or as a result of the cellular response to replication stress. It also functions as a restriction factor preventing viral replication. Here, we report the cryogenic EM (cryo-EM) structure of the six-subunit budding yeast Smc5/6 holo-complex, reconstituted from recombinant proteins expressed in insect cells - providing both an architectural overview of the entire complex and an understanding of how the Nse1/3/4 subcomplex binds to the hetero-dimeric SMC protein core. In addition, we demonstrate that a region within the head domain of Smc5, equivalent to the 'W-loop' of Smc4 or 'F-loop' of Smc1, mediates an important interaction with Nse1. Notably, mutations that alter the surface-charge profile of the region of Nse1 which accepts the Smc5-loop, lead to a slow-growth phenotype and a global reduction in the chromatin-associated fraction of the Smc5/6 complex, as judged by single molecule localisation microscopy experiments in live yeast. Moreover, when taken together, our data indicates functional equivalence between the structurally unrelated KITE and HAWK accessory subunits associated with SMC complexes.

Metagenomics and total RNA sequencing (total RNA-Seq) have the potential to improve the taxonomic identification of diverse microbial communities, which could allow for the incorporation of microbes into routine ecological assessments. However, these target-PCR-free techniques require more testing and optimization. In this study, we processed metagenomics and total RNA-Seq data from a commercially available microbial mock community using 672 data-processing workflows, identified the most accurate data-processing tools, and compared their microbial identification accuracy at equal and increasing sequencing depths. The accuracy of data-processing tools substantially varied among replicates. Total RNA-Seq was more accurate than metagenomics at equal sequencing depths and even at sequencing depths almost one order of magnitude lower than those of metagenomics. We show that while data-processing tools require further exploration, total RNA-Seq might be a favorable alternative to metagenomics for target-PCR-free taxonomic identifications of microbial communities and might enable a substantial reduction in sequencing costs while maintaining accuracy. This could be particularly an advantage for routine ecological assessments, which require cost-effective yet accurate methods, and might allow for the incorporation of microbes into ecological assessments.

Chemical modifications are important for RNA function and metabolism. N4-acetylcytidine (ac4C) is critical for the translation and stability of mRNA. Although ac4C is found in RNA viruses, the detailed mechanisms through which ac4C affects viral replication are unclear. Here, we reported that the 5' untranslated region of the enterovirus 71 (EV71) genome was ac4C modified by the host acetyltransferase NAT10. Inhibition of NAT10 and mutation of the ac4C sites within the internal ribosomal entry site (IRES) suppressed EV71 replication. ac4C enhanced viral RNA translation via selective recruitment of PCBP2 to the IRES and boosted RNA stability. Additionally, ac4C increased the binding of RNA-dependent RNA polymerase (3D) to viral RNA. Notably, ac4C-deficient mutant EV71 showed reduced pathogenicity in vivo. Our findings highlighted the essential role of ac4C in EV71 infection and provided insights into potential antiviral treatments.

DNA supercoiling is an essential mechanism of bacterial chromosome compaction, whose level is mainly regulated by topoisomerase I and DNA gyrase. Inhibiting either of these enzymes with antibiotics leads to global supercoiling modifications and subsequent changes in global gene expression. In previous studies, genes responding to DNA relaxation induced by DNA gyrase inhibition were categorised as 'supercoiling-sensitive'. Here, we studied the opposite variation of DNA supercoiling in the phytopathogen Dickeya dadantii using the non-marketed antibiotic seconeolitsine. We showed that the drug is active against topoisomerase I from this species, and analysed the first transcriptomic response of a Gram-negative bacterium to topoisomerase I inhibition. We find that the responding genes essentially differ from those observed after DNA relaxation, and further depend on the growth phase. We characterised these genes at the functional level, and also detected distinct patterns in terms of expression level, spatial and orientational organisation along the chromosome. Altogether, these results highlight that the supercoiling-sensitivity is a complex feature, which depends on the action of specific topoisomerases, on the physiological conditions, and on their genomic context. Based on previous in vitro expression data of several promoters, we propose a qualitative model of SC-dependent regulation that accounts for many of the contrasting transcriptomic features observed after DNA gyrase or topoisomerase I inhibition.

Normalization of single cell RNA-seq data remains a challenging task. The performance of different methods can vary greatly between datasets when unwanted factors and biology are associated. Most normalization methods also only remove the effects of unwanted variation for the cell embedding but not from gene-level data typically used for differential expression (DE) analysis to identify marker genes. We propose RUV-III-NB, a method that can be used to remove unwanted variation from both the cell embedding and gene-level counts. Using pseudo-replicates, RUV-III-NB explicitly takes into account potential association with biology when removing unwanted variation. The method can be used for both UMI or read counts and returns adjusted counts that can be used for downstream analyses such as clustering, DE and pseudotime analyses. Using published datasets with different technological platforms, kinds of biology and levels of association between biology and unwanted variation, we show that RUV-III-NB manages to remove library size and batch effects, strengthen biological signals, improve DE analyses, and lead to results exhibiting greater concordance with independent datasets of the same kind. The performance of RUV-III-NB is consistent and is not sensitive to the number of factors assumed to contribute to the unwanted variation.

DDX5 and DDX17 are DEAD-box RNA helicase paralogs which regulate several aspects of gene expression, especially transcription and splicing, through incompletely understood mechanisms. A transcriptome analysis of DDX5/DDX17-depleted human cells confirmed the large impact of these RNA helicases on splicing and revealed a widespread deregulation of 3' end processing. In silico analyses and experiments in cultured cells showed the binding and functional contribution of the genome organizing factor CTCF to chromatin sites at or near a subset of DDX5/DDX17-dependent exons that are characterized by a high GC content and a high density of RNA Polymerase II. We propose the existence of an RNA helicase-dependent relationship between CTCF and the dynamics of transcription across DNA and/or RNA structured regions, that contributes to the processing of internal and terminal exons. Moreover, local DDX5/DDX17-dependent chromatin loops spatially connect RNA helicase-regulated exons with their cognate promoter, and we provide the first direct evidence that de novo gene looping modifies alternative splicing and polyadenylation. Overall our findings uncover the impact of DDX5/DDX17-dependent chromatin folding on pre-messenger RNA processing.

The AP1 transcription factor ΔFOSB, a splice variant of FOSB, accumulates in the brain in response to chronic insults such as exposure to drugs of abuse, depression, Alzheimer's disease and tardive dyskinesias, and mediates subsequent long-term neuroadaptations. ΔFOSB forms heterodimers with other AP1 transcription factors, e.g. JUND, that bind DNA under control of a putative cysteine-based redox switch. Here, we reveal the structural basis of the redox switch by determining a key missing crystal structure in a trio, the ΔFOSB/JUND bZIP domains in the reduced, DNA-free form. Screening a cysteine-focused library containing 3200 thiol-reactive compounds, we identify specific compounds that target the redox switch, validate their activity biochemically and in cell-based assays, and show that they are well tolerated in different cell lines despite their general potential to bind to cysteines covalently. A crystal structure of the ΔFOSB/JUND bZIP domains in complex with a redox-switch-targeting compound reveals a deep compound-binding pocket near the DNA-binding site. We demonstrate that ΔFOSB, and potentially other, related AP1 transcription factors, can be targeted specifically and discriminately by exploiting unique structural features such as the redox switch and the binding partner to modulate biological function despite these proteins previously being thought to be undruggable.

Nestor-Guillermo progeria syndrome (NGPS) is caused by a homozygous alanine-to-threonine mutation at position 12 (A12T) in barrier-to-autointegration factor (BAF). It is characterized by accelerated aging with severe skeletal abnormalities. BAF is an essential protein binding to DNA and nuclear envelope (NE) proteins, involved in NE rupture repair. Here, we assessed the impact of BAF A12T on NE integrity using NGPS-derived patient fibroblasts. We observed a strong defect in lamin A/C accumulation to NE ruptures in NGPS cells, restored upon homozygous reversion of the pathogenic BAF A12T mutation with CRISPR/Cas9. By combining in vitro and cellular assays, we demonstrated that while the A12T mutation does not affect BAF 3D structure and phosphorylation by VRK1, it specifically decreases the interaction between BAF and lamin A/C. Finally, we revealed that the disrupted interaction does not prevent repair of NE ruptures but instead generates weak points in the NE that lead to a higher frequency of NE re-rupturing in NGPS cells. We propose that this NE fragility could directly contribute to the premature aging phenotype in patients.

Programmable RNA editing enables rewriting gene expression without changing genome sequences. Current tools for specific RNA editing dependent on the assembly of guide RNA into an RNA/protein complex, causing delivery barrier and low editing efficiency. We report a new gRNA-free system, RNA editing with individual RNA-binding enzyme (REWIRE), to perform precise base editing with a single engineered protein. This artificial enzyme contains a human-originated programmable PUF domain to specifically recognize RNAs and different deaminase domains to achieve efficient A-to-I or C-to-U editing, which achieved 60-80% editing rate in human cells, with a few non-specific editing sites in the targeted region and a low level off-target effect globally. The RNA-binding domain in REWIREs was further optimized to improve editing efficiency and minimize off-target effects. We applied the REWIREs to correct disease-associated mutations and achieve both types of base editing in mice. As a single-component system originated from human proteins, REWIRE presents a precise and efficient RNA editing platform with broad applicability.

Recent discovery of ectopic repeats (outside CRISPR arrays) provided unprecedented insights into the nondefense roles of CRISPR-Cas. A striking example is the addiction module CreTA (CRISPR-regulated toxin-antitoxins), where one or two (in most cases) ectopic repeats produce CRISPR-resembling antitoxic (CreA) RNAs that direct the CRISPR effector Cascade to transcriptionally repress a toxic RNA (CreT). Here, we demonstrated that CreTA repeats are extensively degenerated in sequence, with the first repeat (ψR1) being more diverged than the second one (ψR2). As a result, such addiction modules become highly specific to their physically-linked CRISPR-Cas loci, and in most cases, CreA could not harness a heterologous CRISPR-Cas to suppress its cognate toxin. We further disclosed that this specificity primarily derives from the degeneration of ψR1, and could generally be altered by modifying this repeat element. We also showed that the degenerated repeats of CreTA were insusceptible to recombination and thus more stable compared to a typical CRISPR array, which could be exploited to develop highly stable CRISPR-based tools. These data illustrated that repeat degeneration (a common feature of ectopic repeats) improves the stability and specificity of CreTA in protecting CRISPR-Cas, which could have contributed to the widespread occurrence and deep diversification of CRISPR systems.

Pseudouridine (Ψ) at position 55 in tRNAs plays an important role in their structure and function. This modification is catalyzed by TruB/Pus4/Cbf5 family of pseudouridine synthases in bacteria and yeast. However, the mechanism of TRUB family underlying the formation of Ψ55 in the mammalian tRNAs is largely unknown. In this report, the CMC/reverse transcription assays demonstrated the presence of Ψ55 in the human mitochondrial tRNAAsn, tRNAGln, tRNAGlu, tRNAPro, tRNAMet, tRNALeu(UUR) and tRNASer(UCN). TRUB1 knockout (KO) cell lines generated by CRISPR/Cas9 technology exhibited the loss of Ψ55 modification in mitochondrial tRNAAsn, tRNAGln, tRNAGlu and tRNAPro but did not affect other 18 mitochondrial tRNAs. An in vitro assay revealed that recombinant TRUB1 protein can catalyze the efficient formation of Ψ55 in tRNAAsn and tRNAGln, but not in tRNAMet and tRNAArg. Notably, the overexpression of TRUB1 cDNA reversed the deficient Ψ55 modifications in these tRNAs in TRUB1KO HeLa cells. TRUB1 deficiency affected the base-pairing (18A/G-Ψ55), conformation and stability but not aminoacylation capacity of these tRNAs. Furthermore, TRUB1 deficiency impacted mitochondrial translation and biogenesis of oxidative phosphorylation system. Our findings demonstrated that human TRUB1 is a highly conserved mitochondrial pseudouridine synthase responsible for the Ψ55 modification in the mitochondrial tRNAAsn, tRNAGln, tRNAGlu and tRNAPro.

Aptamers' vast conformation ensemble consisting of interconverting substates severely impairs their performance and applications in biomedicine. Therefore, developing new chemistries stabilizing aptamer conformation and exploring the conformation-performance relationship are highly desired. Herein, we developed an 8-methoxypsoralen-based photochemically covalent lock to stabilize aptamer conformation via crosslinking the inter-stranded thymine nucleotides at TpA sites. Systematical studies and molecular dynamics simulations were performed to explore the conformation-performance relationship of aptamers, revealing that conformation-stabilized aptamers displayed better ability to bind targets, adapt to physiological environment, resist macrophage uptake, prolong circulation half-life, accumulate in and penetrate into tumor than their counterparts. As expected, conformation-stabilized aptamers efficiently improved the therapeutic efficacy of aptamer-drug conjugation on tumor-bearing mice. Collectively, our study has developed a general, simple and economic strategy to stabilize aptamer conformation and shed light on the conformation-performance relationship of aptamers, laying a basis for promoting their basic researches and applications in biomedicine.

While most research suggests mitochondrial DNA (mtDNA) harbors low or no methylation, a few studies claim to report evidence of high-level methylation in the mtDNA. The reasons behind these contradictory results are likely to be methodological but remain largely unexplored. Here, we critically reanalyzed a recent study by Patil et al. (2019) reporting extensive methylation in human mtDNA in a non-CpG context. Our analyses refute the original findings and show that these do not reflect the biology of the tested samples, but rather stem from a combination of methodological and technical pitfalls. The authors employ an oversimplified model that defines as methylated all reference positions with methylation proportions above an arbitrary cutoff of 9%. This substantially exacerbates the overestimation of methylated cytosines due to the selective degradation of unmethylated cytosine-rich regions. Additional limitations are the small sample sizes and lack of sample-specific controls for bisulfite conversion efficiency. In conclusion, using the same dataset employed in the original study by Patil et al., we find no evidence supporting the existence of extensive non-CpG methylation in the human mtDNA.

Cells are constantly challenged by genotoxic stresses that can lead to genome instability. The integrity of the nuclear genome is preserved by the DNA damage response (DDR) and repair. Additionally, these stresses can induce mitochondria to transiently hyperfuse; however, it remains unclear whether canonical DDR is linked to these mitochondrial morphological changes. Here, we report that the abolition of mitochondrial fusion causes a substantial defect in the ATM-mediated DDR signaling. This deficiency is overcome by the restoration of mitochondria fusion. In cells with fragmented mitochondria, genotoxic stress-induced activation of JNK and its translocation to DNA lesion are lost. Importantly, the mitochondrial fusion machinery of MFN1/MFN2 associates with Sab (SH3BP5) and JNK, and these interactions are indispensable for the Sab-mediated activation of JNK and the ATM-mediated DDR signaling. Accordingly, the formation of BRCA1 and 53BP1 foci, as well as homology and end-joining repair are impaired in cells with fragmented mitochondria. Together, these data show that mitochondrial fusion-dependent JNK signaling is essential for the DDR, providing vital insight into the integration of nuclear and cytoplasmic stress signals.

The effects of rapid acute depletion of components of RNA polymerase II (Pol II) general transcription factors (GTFs) that are thought to be critical for formation of preinitiation complexes (PICs) and initiation in vitro were quantified in HAP1 cells using precision nuclear run-on sequencing (PRO-Seq). The average dependencies for each factor across >70 000 promoters varied widely even though levels of depletions were similar. Some of the effects could be attributed to the presence or absence of core promoter elements such as the upstream TBP-specificity motif or downstream G-rich sequences, but some dependencies anti-correlated with such sequences. While depletion of TBP had a large effect on most Pol III promoters only a small fraction of Pol II promoters were similarly affected. TFIIB depletion had the largest general effect on Pol II and also correlated with apparent termination defects downstream of genes. Our results demonstrate that promoter activity is combinatorially influenced by recruitment of TFIID and sequence-specific transcription factors. They also suggest that interaction of the preinitiation complex (PIC) with nucleosomes can affect activity and that recruitment of TFIID containing TBP only plays a positive role at a subset of promoters.

Recent studies have revealed multiple mechanisms that can lead to heterogeneity in ribosomal composition. This heterogeneity can lead to preferential translation of specific panels of mRNAs, and is defined in large part by the ribosomal protein (RP) content, amongst other things. However, it is currently unknown to what extent ribosomal composition is heterogeneous across tissues, which is compounded by a lack of tools available to study it. Here we present dripARF, a method for detecting differential RP incorporation into the ribosome using Ribosome Profiling (Ribo-seq) data. We combine the 'waste' rRNA fragment data generated in Ribo-seq with the known 3D structure of the human ribosome to predict differences in the composition of ribosomes in the material being studied. We have validated this approach using publicly available data, and have revealed a potential role for eS25/RPS25 in development. Our results indicate that ribosome heterogeneity can be detected in Ribo-seq data, providing a new method to study this phenomenon. Furthermore, with dripARF, previously published Ribo-seq data provides a wealth of new information, allowing the identification of RPs of interest in many disease and normal contexts. dripARF is available as part of the ARF R package and can be accessed through https://github.com/fallerlab/ARF.

While amber suppression is the most common approach to introduce noncanonical amino acids into proteins in live cells, quadruplet codon decoding has potential to enable a greatly expanded genetic code with up to 256 new codons for protein biosynthesis. Since triplet codons are the predominant form of genetic code in nature, quadruplet codon decoding often displays limited efficiency. In this work, we exploited a new approach to significantly improve quadruplet UAGN and AGGN (N = A, U, G, C) codon decoding efficiency by using recoding signals imbedded in mRNA. With representative recoding signals, the expression level of mutant proteins containing UAGN and AGGN codons reached 48% and 98% of that of the wild-type protein, respectively. Furthermore, this strategy mitigates a common concern of reading-through endogenous stop codons with amber suppression-based system. Since synthetic recoding signals are rarely found near the endogenous UAGN and AGGN sequences, a low level of undesirable suppression is expected. Our strategy will greatly enhance the utility of noncanonical amino acid mutagenesis in live-cell studies.

A proportion of previously defined benign variants or variants of uncertain significance in humans, which are challenging to identify, may induce an abnormal splicing process. An increasing number of methods have been developed to predict splicing variants, but their performance has not been completely evaluated using independent benchmarks. Here, we manually sourced ∼50 000 positive/negative splicing variants from > 8000 studies and selected the independent splicing variants to evaluate the performance of prediction methods. These methods showed different performances in recognizing splicing variants in donor and acceptor regions, reminiscent of different weight coefficient applications to predict novel splicing variants. Of these methods, 66.67% exhibited higher specificities than sensitivities, suggesting that more moderate cut-off values are necessary to distinguish splicing variants. Moreover, the high correlation and consistent prediction ratio validated the feasibility of integration of the splicing prediction method in identifying splicing variants. We developed a splicing analytics platform called SPCards, which curates splicing variants from publications and predicts splicing scores of variants in genomes. SPCards also offers variant-level and gene-level annotation information, including allele frequency, non-synonymous prediction and comprehensive functional information. SPCards is suitable for high-throughput genetic identification of splicing variants, particularly those located in non-canonical splicing regions.

In this report, we investigated the molecular mechanism underlying a deafness-associated m.5783C > T mutation that affects the canonical C50-G63 base-pairing of TΨC stem of tRNACys and immediately adjacent to 5' end of light-strand origin of mitochondrial DNA (mtDNA) replication (OriL). Two dimensional agarose gel electrophoresis revealed marked decreases in the replication intermediates including ascending arm of Y-fork arcs spanning OriL in the mutant cybrids bearing m.5783C > T mutation. mtDNA replication alterations were further evidenced by decreased levels of PolγA, Twinkle and SSBP1, newly synthesized mtDNA and mtDNA contents in the mutant cybrids. The m.5783C > T mutation altered tRNACys structure and function, including decreased melting temperature, conformational changes, instability and deficient aminoacylation of mutated tRNACys. The m.5783C > T mutation impaired the 5' end processing efficiency of tRNACys precursors and reduced the levels of tRNACys and downstream tRNATyr. The aberrant tRNA metabolism impaired mitochondrial translation, which was especially pronounced effects in the polypeptides harboring higher numbers of cysteine and tyrosine codons. These alterations led to deficient oxidative phosphorylation including instability and reduced activities of the respiratory chain enzyme complexes I, III, IV and intact supercomplexes overall. Our findings highlight the impact of mitochondrial dysfunction on deafness arising from defects in mitochondrial DNA replication and tRNA metabolism.

In mammals, m7G-adjacent nucleotides undergo extensive modifications. Ribose of the first or first and second transcribed nucleotides can be subjected to 2'-O-methylation to form cap1 or cap2, respectively. When the first transcribed nucleotide is 2'-O-methylated adenosine, it can be additionally modified to N6,2'-O-dimethyladenosine (m6Am). Recently, the crucial role of cap1 in distinguishing between 'self' and 'non-self' in mammalian cells during viral infection was revealed. Here, we attempted to understand the impact of cap methylations on RNA-related processes. Therefore, we synthesized tetranucleotide cap analogues and used them for RNA capping during in vitro transcription. Using this tool, we found that 2'-O-methylation of the second transcribed nucleotide within the mRNA 5' cap influences protein production levels in a cell-specific manner. This modification can strongly hamper protein biosynthesis or have no influence on protein production levels, depending on the cell line. Interestingly, 2'-O-methylation of the second transcribed nucleotide and the presence of m6Am as the first transcribed nucleotide serve as determinants that define transcripts as 'self' and contribute to transcript escape from the host innate immune response. Additionally, cap methylation status does not influence transcript affinity towards translation initiation factor eIF4E or in vitro susceptibility to decapping by DCP2; however, we observe the resistance of cap2-RNA to DXO (decapping exoribonuclease)-mediated decapping and degradation.

Mutations in DKC1 (encoding dyskerin) cause telomere diseases including dyskeratosis congenita (DC) by decreasing steady-state levels of TERC, the non-coding RNA component of telomerase. How DKC1 mutations variably impact numerous other snoRNAs remains unclear, which is a barrier to understanding disease mechanisms in DC beyond impaired telomere maintenance. Here, using DC patient iPSCs, we show that mutations in the dyskerin N-terminal extension domain (NTE) dysregulate scaRNA13. In iPSCs carrying the del37L NTE mutation or engineered to carry NTE mutations via CRISPR/Cas9, but not in those with C-terminal mutations, we found scaRNA13 transcripts with aberrant 3' extensions, as seen when the exoribonuclease PARN is mutated in DC. Biogenesis of scaRNA13 was rescued by repair of the del37L DKC1 mutation by genome-editing, or genetic or pharmacological inactivation of the polymerase PAPD5, which counteracts PARN. Inspection of the human telomerase cryo-EM structure revealed that in addition to mediating intermolecular dyskerin interactions, the NTE interacts with terminal residues of the associated snoRNA, indicating a role for this domain in 3' end definition. Our results provide mechanistic insights into the interplay of dyskerin and the PARN/PAPD5 axis in the biogenesis and accumulation of snoRNAs beyond TERC, broadening our understanding of ncRNA dysregulation in human diseases.

Mitochondrial tRNAs are indispensable for the intra-mitochondrial translation of genes related to respiratory subunits, and mutations in mitochondrial tRNA genes have been identified in various disease patients. However, the molecular mechanism underlying pathogenesis remains unclear due to the lack of animal models. Here, we established a mouse model, designated 'mito-mice tRNALeu(UUR)2748', that carries a pathogenic A2748G mutation in the tRNALeu(UUR) gene of mitochondrial DNA (mtDNA). The A2748G mutation is orthologous to the human A3302G mutation found in patients with mitochondrial diseases and diabetes. A2748G mtDNA was maternally inherited, equally distributed among tissues in individual mice, and its abundance did not change with age. At the molecular level, A2748G mutation is associated with aberrant processing of precursor mRNA containing tRNALeu(UUR) and mt-ND1, leading to a marked decrease in the steady-levels of ND1 protein and Complex I activity in tissues. Mito-mice tRNALeu(UUR)2748 with ≥50% A2748G mtDNA exhibited age-dependent metabolic defects including hyperglycemia, insulin insensitivity, and hepatic steatosis, resembling symptoms of patients carrying the A3302G mutation. This work demonstrates a valuable mouse model with an inheritable pathological A2748G mutation in mt-tRNALeu(UUR) that shows metabolic syndrome-like phenotypes at high heteroplasmy level. Furthermore, our findings provide molecular basis for understanding A3302G mutation-mediated mitochondrial disorders.

The genome of influenza A virus (IAV) consists of eight unique viral RNA segments. This genome organization allows genetic reassortment between co-infecting IAV strains, whereby new IAVs with altered genome segment compositions emerge. While it is known that reassortment events can create pandemic IAVs, it remains impossible to anticipate reassortment outcomes with pandemic prospects. Recent research indicates that reassortment is promoted by a viral genome packaging mechanism that delivers the eight genome segments as a supramolecular complex into the virus particle. This finding holds promise of predicting pandemic IAVs by understanding the intermolecular interactions governing this genome packaging mechanism. Here, we critically review the prevailing mechanistic model postulating that IAV genome packaging is orchestrated by a network of intersegmental RNA-RNA interactions. Although we find supporting evidence, including segment-specific packaging signals and experimentally proposed RNA-RNA interaction networks, this mechanistic model remains debatable due to a current shortage of functionally validated intersegmental RNA-RNA interactions. We speculate that identifying such functional intersegmental RNA-RNA contacts might be hampered by limitations of the utilized probing techniques and the inherent complexity of the genome packaging mechanism. Nevertheless, we anticipate that improved probing strategies combined with a mutagenesis-based validation could facilitate their discovery.

Precise maintenance of PTEN dosage is crucial for tumor suppression across a wide variety of cancers. Post-transcriptional regulation of Pten heavily relies on regulatory elements encoded by its 3'UTR. We previously reported the important diversity of 3'UTR isoforms of Pten mRNAs produced through alternative polyadenylation (APA). Here, we reveal the direct regulation of Pten APA by the mammalian cleavage factor I (CFIm) complex, which in turn contributes to PTEN protein dosage. CFIm consists of the UGUA-binding CFIm25 and APA regulatory subunits CFIm59 or CFIm68. Deep sequencing analyses of perturbed (KO and KD) cell lines uncovered the differential regulation of Pten APA by CFIm59 and CFIm68 and further revealed that their divergent functions have widespread impact for APA in transcriptomes. Differentially regulated genes include numerous factors within the phosphoinositide 3-kinase (PI3K)/protein kinase B (Akt) signalling pathway that PTEN counter-regulates. We further reveal a stratification of APA dysregulation among a subset of PTEN-driven cancers, with recurrent alterations among PI3K/Akt pathway genes regulated by CFIm. Our results refine the transcriptome selectivity of the CFIm complex in APA regulation, and the breadth of its impact in PTEN-driven cancers.

The static and dynamic structures of DNA duplexes affected by 5S-Tg (Tg, Thymine glycol) epimers were studied using MD simulations and Markov State Models (MSMs) analysis. The results show that the 5S,6S-Tg base caused little perturbation to the helix, and the base-flipping barrier was determined to be 4.4 kcal mol-1 through the use of enhanced sampling meta-eABF calculations, comparable to 5.4 kcal mol-1 of the corresponding thymine flipping. Two conformations with the different hydrogen bond structures between 5S,6R-Tg and A19 were identified in several independent MD trajectories. The 5S,6R-Tg:O6HO6•••N1:A19 hydrogen bond is present in the high-energy conformation displaying a clear helical distortion, and near barrier-free Tg base flipping. The low-energy conformation always maintains Watson-Crick base pairing between 5S,6R-Tg and A19, and 5S-Tg base flipping is accompanied by a small barrier of ca. 2.0 KBT (T = 298 K). The same conformations are observed in the MSMs analysis. Moreover, the transition path and metastable structures of the damaged base flipping are for the first time verified through MSMs analysis. The data clearly show that the epimers have completely different influence on the stability of the DNA duplex, thus implying different enzymatic mechanisms for DNA repair.

Protein synthesis in eukaryotic cell is spatially and structurally compartmentalized that ensures high efficiency of this process. One of the distinctive features of higher eukaryotes is the existence of stable multi-protein complexes of aminoacyl-tRNA synthetases and translation elongation factors. Here, we report a quaternary organization of the human guanine-nucleotide exchange factor (GEF) complex, eEF1B, comprising α, β and γ subunits that specifically associate into a heterotrimeric form eEF1B(αβγ)3. As both the eEF1Bα and eEF1Bβ proteins have structurally conserved GEF domains, their total number within the complex is equal to six. Such, so far, unique structural assembly of the guanine-nucleotide exchange factors within a stable complex may be considered as a 'GEF hub' that ensures efficient maintenance of the translationally active GTP-bound conformation of eEF1A in higher eukaryotes.

Improvements in cryo-electron tomography sample preparation, electron-microscopy instrumentations, and image processing algorithms have advanced the structural analysis of macromolecules in situ. Beyond such analyses of individual macromolecules, the study of their interactions with functionally related neighbors in crowded cellular habitats, i.e. 'molecular sociology', is of fundamental importance in biology. Here we present a NEighboring Molecule TOpology Clustering (NEMO-TOC) algorithm. We optimized this algorithm for the detection and profiling of polyribosomes, which play both constitutive and regulatory roles in gene expression. Our results suggest a model where polysomes are formed by connecting multiple nonstochastic blocks, in which translation is likely synchronized.

Algorithms have been widely used to predict G-quadruplexes (G4s)-prone sequences. However, an experimental validation of these predictions is generally required. We previously reported a high-throughput technique to evidence G4 formation in vitro called FRET-MC. This method, while convenient and reproducible, has one known weakness: its inability to pin point G4 motifs of low thermal stability. As such quadruplexes may still be biologically relevant if formed at physiological temperature, we wanted to develop an independent assay to overcome this limitation. To this aim, we introduced an isothermal version of the competition assay, called iso-FRET, based on a duplex-quadruplex competition and a well-characterized bis-quinolinium G4 ligand, PhenDC3. G4-forming competitors act as decoys for PhenDC3, lowering its ability to stabilize the G4-forming motif reporter oligonucleotide conjugated to a fluorescence quencher (37Q). The decrease in available G4 ligand concentration restores the ability of 37Q to hybridize to its FAM-labeled short complementary C-rich strand (F22), leading to a decrease in fluorescence signal. In contrast, when no G4-forming competitor is present, PhenDC3 remains available to stabilize the 37Q quadruplex, preventing the formation of the F22 + 37Q complex. Iso-FRET was first applied to a reference panel of 70 sequences, and then used to investigate 23 different viral sequences.

ATRX is a chromatin remodeler, which is mutated in ATRX syndrome, a neurodevelopmental disorder. ATRX mutations that alter histone binding or chromatin remodeling activities cluster in the PHD finger or the helicase domain respectively. Using engineered mouse embryonic stem cells that exclusively express ATRX protein with mutations in the PHD finger (PHDmut) or helicase domains (K1584R), we examine how specific ATRX mutations affect neurodifferentiation. ATRX PHDmut and K1584R proteins interact with the DAXX histone chaperone but show reduced localization to pericentromeres. Neurodifferentiation is both delayed and compromised in PHDmut and K1584R, and manifest differently from complete ATRX loss. We observe reduced enrichment of PHDmut protein to ATRX targets, while K1584R accumulates at these sites. Interestingly, ATRX mutations have distinct effects on the genome-wide localization of the polycomb repressive complex 2 (PRC2), with PHDmut and ATRX knockout showing reduced PRC2 binding at polycomb targets and K1584R showing loss at some sites and gains at others. Notably, each mutation associated with unique gene signatures, suggesting distinct pathways leading to impaired neurodifferentiation. Our results indicate that the histone binding and chromatin remodeling functions of ATRX play non-redundant roles in neurodevelopment, and when mutated lead to ATRX syndrome through separate regulatory pathways.

Post-translational modifications (PTMs) of core histones are important epigenetic determinants that correlate with functional chromatin states. However, despite multiple linker histone H1s PTMs have been identified, little is known about their genomic distribution and contribution to the epigenetic regulation of chromatin. Here, we address this question in Drosophila that encodes a single somatic linker histone, dH1. We previously reported that dH1 is dimethylated at K27 (dH1K27me2). Here, we show that dH1K27me2 is a major PTM of Drosophila heterochromatin. At mitosis, dH1K27me2 accumulates at pericentromeric heterochromatin, while, in interphase, it is also detected at intercalary heterochromatin. ChIPseq experiments show that >98% of dH1K27me2 enriched regions map to heterochromatic repetitive DNA elements, including transposable elements, simple DNA repeats and satellite DNAs. Moreover, expression of a mutated dH1K27A form, which impairs dH1K27me2, alters heterochromatin organization, upregulates expression of heterochromatic transposable elements and results in the accumulation of RNA:DNA hybrids (R-loops) in heterochromatin, without affecting H3K9 methylation and HP1a binding. The pattern of dH1K27me2 is H3K9 methylation independent, as it is equally detected in flies carrying a H3K9R mutation, and is not affected by depletion of Su(var)3-9, HP1a or Su(var)4-20. Altogether these results suggest that dH1K27me2 contributes to heterochromatin organization independently of H3K9 methylation.

Topoisomerase IIA (TOP2a) has traditionally been known as an important nuclear enzyme that resolves entanglements and relieves torsional stress of DNA double strands. However, its function in genomic transcriptional regulation remains largely unknown, especially during adult neurogenesis. Here, we show that TOP2a is preferentially expressed in neurogenic niches in the brain of adult mice, such as the subventricular zone (SVZ). Conditional knockout of Top2a in adult neural stem cells (NSCs) of the SVZ significantly inhibits their self-renewal and proliferation, and ultimately reduces neurogenesis. To gain insight into the molecular mechanisms by which TOP2a regulates adult NSCs, we perform RNA-sequencing (RNA-Seq) plus chromatin immunoprecipitation sequencing (ChIP-Seq) and identify ubiquitin-specific protease 37 (Usp37) as a direct TOP2a target gene. Importantly, overexpression of Usp37 is sufficient to rescue the impaired self-renewal ability of adult NSCs caused by Top2a knockdown. Taken together, this proof-of-principle study illustrates a TOP2a/Usp37-mediated novel molecular mechanism in adult neurogenesis, which will significantly expand our understanding of the function of topoisomerase in the adult brain.

We present a physics-based machine learning approach to predict in vitro transcription factor binding affinities from structural and mechanical DNA properties directly derived from atomistic molecular dynamics simulations. The method is able to predict affinities obtained with techniques as different as uPBM, gcPBM and HT-SELEX with an excellent performance, much better than existing algorithms. Due to its nature, the method can be extended to epigenetic variants, mismatches, mutations, or any non-coding nucleobases. When complemented with chromatin structure information, our in vitro trained method provides also good estimates of in vivo binding sites in yeast.

Enhancers establish proximity with distant target genes to regulate temporospatial gene expression and specify cell identity. Lim domain binding protein 1 (LDB1) is a conserved and widely expressed protein that functions as an enhancer looping factor. Previous studies in erythroid cells and neuronal cells showed that LDB1 forms protein complexes with different transcription factors to regulate cell-specific gene expression. Here, we show that LDB1 regulates expression of liver genes by occupying enhancer elements and cooperating with hepatic transcription factors HNF4A, FOXA1, TCF7 and GATA4. Using the glucose transporter SLC2A2 gene, encoding GLUT2, as an example, we find that LDB1 regulates gene expression by mediating enhancer-promoter interactions. In vivo, we find that LDB1 deficiency in primary mouse hepatocytes dysregulates metabolic gene expression and changes the enhancer landscape. Conditional deletion of LDB1 in adult mouse liver induces glucose intolerance. However, Ldb1 knockout hepatocytes show improved liver pathology under high-fat diet conditions associated with increased expression of genes related to liver fatty acid metabolic processes. Thus, LDB1 is linked to liver metabolic functions under normal and obesogenic conditions.

RNA function relies heavily on posttranscriptional modifications. Recently, it was shown that certain PARPs and TRPT1 can ADP-ribosylate RNA in vitro. Traditionally, intracellular ADP-ribosylation has been considered mainly as a protein posttranslational modification. To date, it is not clear whether RNA ADP-ribosylation occurs in cells. Here we present evidence that different RNA species are ADP-ribosylated in human cells. The modification of cellular RNA is mediated by several transferases such as TRPT1, PARP10, PARP11, PARP12 and PARP15 and is counteracted by different hydrolases including TARG1, PARG and ARH3. In addition, diverse cellular stressors can modulate the content of ADP-ribosylated RNA in cells. We next investigated potential consequences of ADP-ribosylation for RNA and found that ADPr-capped mRNA is protected against XRN1 mediated degradation but is not translated. T4 RNA ligase 1 can ligate ADPr-RNA in absence of ATP, resulting in the incorporation of an abasic site. We thus provide the first evidence of RNA ADP-ribosylation in mammalian cells and postulate potential functions of this novel RNA modification.

Dynamic regulation of RNA folding and structure is critical for the biogenesis and function of RNAs and ribonucleoprotein (RNP) complexes. Through their nucleotide triphosphate-dependent remodelling functions, RNA helicases are key modulators of RNA/RNP structure. While some RNA helicases are dedicated to a specific target RNA, others are multifunctional and engage numerous substrate RNAs in different aspects of RNA metabolism. The discovery of such multitasking RNA helicases raises the intriguing question of how these enzymes can act on diverse RNAs but also maintain specificity for their particular targets within the RNA-dense cellular environment. Furthermore, the identification of RNA helicases that sit at the nexus between different aspects of RNA metabolism raises the possibility that they mediate cross-regulation of different cellular processes. Prominent and extensively characterized multifunctional DEAH/RHA-box RNA helicases are DHX15 and its Saccharomyces cerevisiae (yeast) homologue Prp43. Due to their central roles in key cellular processes, these enzymes have also served as prototypes for mechanistic studies elucidating the mode of action of this type of enzyme. Here, we summarize the current knowledge on the structure, regulation and cellular functions of Prp43/DHX15, and discuss the general concept and implications of RNA helicase multifunctionality.

La-related proteins (LARPs) comprise a family of RNA-binding proteins involved in a wide range of posttranscriptional regulatory activities. LARPs share a unique tandem of two RNA-binding domains, La motif (LaM) and RNA recognition motif (RRM), together referred to as a La-module, but vary in member-specific regions. Prior structural studies of La-modules reveal they are pliable platforms for RNA recognition in diverse contexts. Here, we characterize the La-module of LARP1, which plays an important role in regulating synthesis of ribosomal proteins in response to mTOR signaling and mRNA stabilization. LARP1 has been well characterized functionally but no structural information exists for its La-module. We show that unlike other LARPs, the La-module in LARP1 does not contain an RRM domain. The LaM alone is sufficient for binding poly(A) RNA with submicromolar affinity and specificity. Multiple high-resolution crystal structures of the LARP1 LaM domain in complex with poly(A) show that it is highly specific for the RNA 3'-end, and identify LaM residues Q333, Y336 and F348 as the most critical for binding. Use of a quantitative mRNA stabilization assay and poly(A) tail-sequencing demonstrate functional relevance of LARP1 RNA binding in cells and provide novel insight into its poly(A) 3' protection activity.

Failure to prevent accumulation of the non-canonical nucleotide inosine triphosphate (ITP) by inosine triphosphate pyrophosphatase (ITPase) during nucleotide synthesis results in misincorporation of inosine into RNA and can cause severe and fatal developmental anomalies in humans. While the biochemical activity of ITPase is well understood, the pathogenic basis of ITPase deficiency and the molecular and cellular consequences of ITP misincorporation into RNA remain cryptic. Here, we demonstrate that excess ITP in the nucleotide pool during in vitro transcription results in T7 polymerase-mediated inosine misincorporation in luciferase RNA. In vitro translation of inosine-containing luciferase RNA reduces resulting luciferase activity, which is only partly explained by reduced abundance of the luciferase protein produced. Using Oxford Nanopore Direct RNA sequencing, we reveal inosine misincorporation to be stochastic but biased largely towards misincorporation in place of guanosine, with evidence for misincorporation also in place of cytidine, adenosine and uridine. Inosine misincorporation into RNA is also detected in Itpa-null mouse embryonic heart tissue as an increase in relative variants compared with the wild type using Illumina RNA sequencing. By generating CRISPR/Cas9 rat H9c2 Itpa-null cardiomyoblast cells, we validate a translation defect in cells that accumulate inosine within endogenous RNA. Furthermore, we observe hindered cellular translation of transfected luciferase RNA containing misincorporated inosine in both wild-type and Itpa-null cells. We therefore conclude that inosine misincorporation into RNA perturbs translation, thus providing mechanistic insight linking ITPase deficiency, inosine accumulation and pathogenesis.

Stress-responsive genes are lowly transcribed under normal conditions and robustly induced in response to stress. The significant difference between basal and induced transcription indicates that the general transcriptional machinery requires a mechanism to distinguish each transcription state. However, what factors specifically function in basal transcription remains poorly understood. Using a classic model stress-responsive gene (Drosophila MtnA), we found that knockdown of the DEAD-box helicase Hlc resulted in a significant transcription attenuation of MtnA under normal, but not stressed, conditions. Mechanistically, Hlc directly binds to the MtnA locus to maintain the accessibility of chromatin near the transcriptional start site, which allows the recruitment of RNA polymerase II and subsequent MtnA transcription. Using RNA-seq, we then identified plenty of additional stress-responsive genes whose basal transcription was reduced upon knockdown of Hlc. Taken together, these data suggest that Hlc-mediated basal transcription regulation is an essential and widespread mechanism for precise control of stress-responsive genes.

The HAV nonstructural protein 2C is essential for virus replication; however, its precise function remains elusive. Although HAV 2C shares 24-27% sequence identity with other 2Cs, key motifs are conserved. Here, we demonstrate that HAV 2C is an ATPase but lacking helicase activity. We identified an ATPase-independent nuclease activity of HAV 2C with a preference for polyuridylic single-stranded RNAs. We determined the crystal structure of an HAV 2C fragment to 2.2 Å resolution, containing an ATPase domain, a region equivalent to enterovirus 2C zinc-finger (ZFER) and a C-terminal amphipathic helix (PBD). The PBD of HAV 2C occupies a hydrophobic pocket (Pocket) in the adjacent 2C, and we show the PBD-Pocket interaction is vital for 2C functions. We identified acidic residues that are essential for the ribonuclease activity and demonstrated mutations at these sites abrogate virus replication. We built a hexameric-ring model of HAV 2C, revealing the ribonuclease-essential residues clustering around the central pore of the ring, whereas the ATPase active sites line up at the gaps between adjacent 2Cs. Finally, we show the ribonuclease activity is shared by other picornavirus 2Cs. Our findings identified a previously unfound activity of picornavirus 2C, providing novel insights into the mechanisms of virus replication.

The tripartite ParABS system mediates chromosome segregation in a wide range of bacteria. Dimeric ParB was proposed to nucleate on parS sites and spread to neighboring DNA. However, how properly distributed ParB dimers further compact chromosomal DNA into a higher-order nucleoprotein complex for partitioning remains poorly understood. Here, using a single-molecule approach, we show that tens of Bacillus subtilis ParB (Spo0J) proteins can stochastically multimerize on and stably bind to nonspecific DNA. The introduction of CTP promotes the formation and diffusion of the multimeric ParB along DNA, offering an opportunity for ParB proteins to further forgather and cluster. Intriguingly, ParB multimers can recognize parS motifs and are more inclined to remain immobile on them. Importantly, the ParB multimer features distinct capabilities of not only bridging two independent DNA molecules but also mediating their transportation, both of which are enhanced by the presence of either CTP or parS in the DNA. These findings shed new light on ParB dynamics in self-multimerization and DNA organization and help to better comprehend the assembly of the ParB-DNA partition complex.

DNA methylation, specifically, methylation of cytosine (C) nucleotides at the 5-carbon position (5-mC), is the most studied and significant epigenetic modification. Here we developed a chemoenzymatic procedure to fluorescently label non-methylated cytosines in CpG context, allowing epigenetic profiling of single DNA molecules spanning hundreds of thousands of base pairs. We used a CpG methyltransferase with a synthetic S-adenosyl-l-methionine cofactor analog to transfer an azide to cytosines instead of the natural methyl group. A fluorophore was then clicked onto the DNA, reporting on the amount and position of non-methylated CpGs. We found that labeling efficiency was increased up to 2-fold by the addition of a nucleosidase, presumably by degrading the inactive by-product of the cofactor after labeling, preventing its inhibitory effect. We used the method to determine the decline in global DNA methylation in a chronic lymphocytic leukemia patient and then performed whole-genome methylation mapping of the model plant Arabidopsis thaliana. Our genome maps show high concordance with published bisulfite sequencing methylation maps. Although mapping resolution is limited by optical detection to 500-1000 bp, the labeled DNA molecules produced by this approach are hundreds of thousands of base pairs long, allowing access to long repetitive and structurally variable genomic regions.

Analyzing single-cell transcriptomes promises to decipher the plasticity, heterogeneity, and rapid switches in developmental cellular state transitions. Such analyses require the identification of gene markers for semi-stable transition states. However, there are nontrivial challenges such as unexplainable stochasticity, variable population sizes, and alternative trajectory constructions. By advancing current tipping-point theory-based models with feature selection, network decomposition, accurate estimation of correlations, and optimization, we developed BioTIP to overcome these challenges. BioTIP identifies a small group of genes, called critical transition signal (CTS), to characterize regulated stochasticity during semi-stable transitions. Although methods rooted in different theories converged at the same transition events in two benchmark datasets, BioTIP is unique in inferring lineage-determining transcription factors governing critical transition. Applying BioTIP to mouse gastrulation data, we identify multiple CTSs from one dataset and validated their significance in another independent dataset. We detect the established regulator Etv2 whose expression change drives the haemato-endothelial bifurcation, and its targets together in CTS across three datasets. After comparing to three current methods using six datasets, we show that BioTIP is accurate, user-friendly, independent of pseudo-temporal trajectory, and captures significantly interconnected and reproducible CTSs. We expect BioTIP to provide great insight into dynamic regulations of lineage-determining factors.

Erratum in Nucleic Acids Res. 2023 Jan 10;:

Xeno-nucleic acids (XNAs) are synthetic genetic polymers with backbone structures composed of non-ribose or non-deoxyribose sugars. Phosphonomethylthreosyl nucleic acid (pTNA), a type of XNA that does not base pair with DNA or RNA, has been suggested as a possible genetic material for storing synthetic biology information in cells. A critical step in this process is the synthesis of XNA episomes using laboratory-evolved polymerases to copy DNA information into XNA. Here, we investigate the polymerase recognition of pTNA nucleotides using X-ray crystallography to capture the post-catalytic complex of engineered polymerases following the sequential addition of two pTNA nucleotides onto the 3'-end of a DNA primer. High-resolution crystal structures reveal that the polymerase mediates Watson-Crick base pairing between the extended pTNA adducts and the DNA template. Comparative analysis studies demonstrate that the sugar conformation and backbone position of pTNA are structurally more similar to threose nucleic acid than DNA even though pTNA and DNA share the same six-atom backbone repeat length. Collectively, these findings provide new insight into the structural determinants that guide the enzymatic synthesis of an orthogonal genetic polymer, and may lead to the discovery of new variants that function with enhanced activity.

Protein aggregates and abnormal proteins are toxic and associated with neurodegenerative diseases. There are several mechanisms to help cells get rid of aggregates but little is known on how cells prevent aggregate-prone proteins from being synthesised. The EBNA1 of the Epstein-Barr virus (EBV) evades the immune system by suppressing its own mRNA translation initiation in order to minimize the production of antigenic peptides for the major histocompatibility (MHC) class I pathway. Here we show that the emerging peptide of the disordered glycine-alanine repeat (GAr) within EBNA1 dislodges the nascent polypeptide-associated complex (NAC) from the ribosome. This results in the recruitment of nucleolin to the GAr-encoding mRNA and suppression of mRNA translation initiation in cis. Suppressing NAC alpha (NACA) expression prevents nucleolin from binding to the GAr mRNA and overcomes GAr-mediated translation inhibition. Taken together, these observations suggest that EBNA1 exploits a nascent protein quality control pathway to regulate its own rate of synthesis that is based on sensing the nascent GAr peptide by NAC followed by the recruitment of nucleolin to the GAr-encoding RNA sequence.

High mobility group (HMG) proteins are chromatin regulators with essential functions in development, cell differentiation and cell proliferation. The protein HMG20A is predicted by the AlphaFold2 software to contain three distinct structural elements, which we have functionally characterized: i) an amino-terminal, intrinsically disordered domain with transactivation activity; ii) an HMG box with higher binding affinity for double-stranded, four-way-junction DNA than for linear DNA; and iii) a long coiled-coil domain. Our proteomic study followed by a deletion analysis and structural modeling demonstrates that HMG20A forms a complex with the histone reader PHF14, via the establishment of a two-stranded alpha-helical coiled-coil structure. siRNA-mediated knockdown of either PHF14 or HMG20A in MDA-MB-231 cells causes similar defects in cell migration, invasion and homotypic cell-cell adhesion ability, but neither affects proliferation. Transcriptomic analyses demonstrate that PHF14 and HMG20A share a large subset of targets. We show that the PHF14-HMG20A complex modulates the Hippo pathway through a direct interaction with the TEAD1 transcription factor. PHF14 or HMG20A deficiency increases epithelial markers, including E-cadherin and the epithelial master regulator TP63 and impaired normal TGFβ-trigged epithelial-to-mesenchymal transition. Taken together, these data indicate that PHF14 and HMG20A cooperate in regulating several pathways involved in epithelial-mesenchymal plasticity.

Nucleases derived from the prokaryotic defense system CRISPR-Cas are frequently re-purposed for gene editing and molecular diagnostics. Hence, an in-depth understanding of the molecular mechanisms of these enzymes is of crucial importance. We focused on Cas12a from Francisella novicida (FnCas12a) and investigated the functional role of helix 1, a structural element that together with the bridge helix (BH) connects the recognition and the nuclease lobes of FnCas12a. Helix 1 is structurally connected to the lid domain that opens upon DNA target loading thereby activating the active site of FnCas12a. We probed the structural states of FnCas12a variants altered in helix 1 and/or the bridge helix using single-molecule FRET measurements and assayed the pre-crRNA processing, cis- and trans-DNA cleavage activity. We show that helix 1 and not the bridge helix is the predominant structural element that confers conformational stability of FnCas12a. Even small perturbations in helix 1 lead to a decrease in DNA cleavage activity while the structural integrity is not affected. Our data, therefore, implicate that the concerted remodeling of helix 1 and the bridge helix upon DNA binding is structurally linked to the opening of the lid and therefore involved in the allosteric activation of the active site.

Mutations in the lamin A/C gene (LMNA) cause laminopathies such as the premature aging Hutchinson Gilford progeria syndrome (HGPS) and altered lamin A/C levels are found in diverse malignancies. The underlying lamin-associated mechanisms remain poorly understood. Here we report that lamin A/C-null mouse embryo fibroblasts (Lmna-/- MEFs) and human progerin-expressing HGPS fibroblasts both display reduced NAD+ levels, unstable mitochondrial DNA and attenuated bioenergetics. This mitochondrial dysfunction is associated with reduced chromatin recruitment (Lmna-/- MEFs) or low levels (HGPS) of PGC1α, the key transcription factor for mitochondrial homeostasis. Lmna-/- MEFs showed reduced expression of the NAD+-biosynthesis enzyme NAMPT and attenuated activity of the NAD+-dependent deacetylase SIRT1. We find high PARylation in lamin A/C-aberrant cells, further decreasing the NAD+ pool and consistent with impaired DNA base excision repair in both cell models, a condition that fuels DNA damage-induced PARylation under oxidative stress. Further, ATAC-sequencing revealed a substantially altered chromatin landscape in Lmna-/- MEFs, including aberrantly reduced accessibility at the Nampt gene promoter. Thus, we identified a new role of lamin A/C as a key modulator of mitochondrial function through impairments of PGC1α and the NAMPT-NAD+ pathway, with broader implications for the aging process.

Trinucleotide repeat (TNR) diseases are caused by the aberrant expansion of CXG (X = C, A, G and T) sequences in genomes. We have reported two small molecules binding to TNR, NCD, and NA, which strongly bind to CGG repeat (responsible sequence of fragile X syndrome) and CAG repeat (Huntington's disease). The NMR structure of NA binding to the CAG/CAG triad has been clarified, but the structure of NCD bound to the CGG/CGG triad remained to be addressed. We here report the structural determination of the NCD-CGG/CGG complex by NMR spectroscopy and the comparison with the NA-CAG/CAG complex. While the NCD-CGG/CGG structure shares the binding characteristics with that of the NA-CAG/CAG complex, a significant difference was found in the overall structure caused by the structural fluctuation at the ligand-bound site. The NCD-CGG/CGG complex was suggested in the equilibrium between stacked and kinked structures, although NA-CAG/CAG complex has only the stacked structures. The dynamic fluctuation of the NCD-CGG/CGG structure at the NCD-binding site suggested room for optimization in the linker structure of NCD to gain improved affinity to the CGG/CGG triad.

Chromatin complexes control a vast number of epigenetic developmental processes. Filamentous fungi present an important clade of microbes with poor understanding of underlying epigenetic mechanisms. Here, we describe a chromatin binding complex in the fungus Aspergillus nidulans composing of a H3K4 histone demethylase KdmB, a cohesin acetyltransferase (EcoA), a histone deacetylase (RpdA) and a histone reader/E3 ligase protein (SntB). In vitro and in vivo evidence demonstrate that this KERS complex is assembled from the EcoA-KdmB and SntB-RpdA heterodimers. KdmB and SntB play opposing roles in regulating the cellular levels and stability of EcoA, as KdmB prevents SntB-mediated degradation of EcoA. The KERS complex is recruited to transcription initiation start sites at active core promoters exerting promoter-specific transcriptional effects. Interestingly, deletion of any one of the KERS subunits results in a common negative effect on morphogenesis and production of secondary metabolites, molecules important for niche securement in filamentous fungi. Consequently, the entire mycotoxin sterigmatocystin gene cluster is downregulated and asexual development is reduced in the four KERS mutants. The elucidation of the recruitment of epigenetic regulators to chromatin via the KERS complex provides the first mechanistic, chromatin-based understanding of how development is connected with small molecule synthesis in fungi.

Increasing numbers of small, regulatory RNAs (sRNAs) corresponding to 3' untranslated regions (UTR) are being discovered in bacteria. One such sRNA, denoted GlnZ, corresponds to the 3' UTR of the Escherichia coli glnA mRNA encoding glutamine synthetase. Several forms of GlnZ, processed from the glnA mRNA, are detected in cells growing with limiting ammonium. GlnZ levels are regulated transcriptionally by the NtrC transcription factor and post-transcriptionally by RNase III. Consistent with the expression, E. coli cells lacking glnZ show delayed outgrowth from nitrogen starvation compared to wild type cells. Transcriptome-wide RNA-RNA interactome datasets indicated that GlnZ binds to multiple target RNAs. Immunoblots and assays of fusions confirmed GlnZ-mediated repression of glnP and sucA, encoding proteins that contribute to glutamine transport and the citric acid cycle, respectively. Although the overall sequences of GlnZ from E. coli K-12, Enterohemorrhagic E. coli and Salmonella enterica have significant differences due to various sequence insertions, all forms of the sRNA were able to regulate the two targets characterized. Together our data show that GlnZ impacts growth of E. coli under low nitrogen conditions by modulating genes that affect carbon and nitrogen flux.

Due to genome segmentation, rotaviruses must co-package eleven distinct genomic RNAs. The packaging is mediated by virus-encoded RNA chaperones, such as the rotavirus NSP2 protein. While the activities of distinct RNA chaperones are well studied on smaller RNAs, little is known about their global effect on the entire viral transcriptome. Here, we used Selective 2'-hydroxyl Acylation Analyzed by Primer Extension and Mutational Profiling (SHAPE-MaP) to examine the secondary structure of the rotavirus transcriptome in the presence of increasing amounts of NSP2. SHAPE-MaP data reveals that despite the well-documented helix-unwinding activity of NSP2 in vitro, its incubation with cognate rotavirus transcripts does not induce a significant change in the SHAPE reactivities. However, a quantitative analysis of mutation rates measured by mutational profiling reveals a global 5-fold rate increase in the presence of NSP2. We demonstrate that the normalization procedure used in deriving SHAPE reactivities from mutation rates can mask an important global effect of an RNA chaperone. Analysis of the mutation rates reveals a larger effect on stems rather than loops. Together, these data provide the first experimentally derived secondary structure model of the rotavirus transcriptome and reveal that NSP2 acts by globally increasing RNA backbone flexibility in a concentration-dependent manner.

Prostate cancer (PC) relies on androgen receptor (AR) signaling. While hormonal therapy (HT) is efficacious, most patients evolve to an incurable castration-resistant stage (CRPC). To date, most proposed mechanisms of acquired resistance to HT have focused on AR transcriptional activity. Herein, we uncover a new role for the AR in alternative cleavage and polyadenylation (APA). Inhibition of the AR by Enzalutamide globally regulates APA in PC cells, with specific enrichment in genes related to transcription and DNA topology, suggesting their involvement in transcriptome reprogramming. AR inhibition selects promoter-distal polyadenylation sites (pAs) enriched in cis-elements recognized by the cleavage and polyadenylation specificity factor (CPSF) complex. Conversely, promoter-proximal intronic pAs relying on the cleavage stimulation factor (CSTF) complex are repressed. Mechanistically, Enzalutamide induces rearrangement of APA subcomplexes and impairs the interaction between CPSF and CSTF. AR inhibition also induces co-transcriptional CPSF recruitment to gene promoters, predisposing the selection of pAs depending on this complex. Importantly, the scaffold CPSF160 protein is up-regulated in CRPC cells and its depletion represses HT-induced APA patterns. These findings uncover an unexpected role for the AR in APA regulation and suggest that APA-mediated transcriptome reprogramming represents an adaptive response of PC cells to HT.

Eukaryotic ribosome synthesis involves more than 200 assembly factors, which promote ribosomal RNA (rRNA) processing, modification and folding, and assembly of ribosomal proteins. The formation and maturation of the earliest pre-60S particles requires structural remodeling by the Npa1 complex, but is otherwise still poorly understood. Here, we introduce Rbp95 (Ycr016w), a constituent of early pre-60S particles, as a novel ribosome assembly factor. We show that Rbp95 is both genetically and physically linked to most Npa1 complex members and to ribosomal protein Rpl3. We demonstrate that Rbp95 is an RNA-binding protein containing two independent RNA-interacting domains. In vivo, Rbp95 associates with helix H95 in the 3' region of the 25S rRNA, in close proximity to the binding sites of Npa1 and Rpl3. Additionally, Rbp95 interacts with several snoRNAs. The absence of Rbp95 results in alterations in the protein composition of early pre-60S particles. Moreover, combined mutation of Rbp95 and Npa1 complex members leads to a delay in the maturation of early pre-60S particles. We propose that Rbp95 acts together with the Npa1 complex during early pre-60S maturation, potentially by promoting pre-rRNA folding events within pre-60S particles.

Small RNAs (sRNAs), including microRNAs (miRNAs) and small interfering RNAs (siRNAs), are essential gene regulators for plant and animal development. The loading of sRNA duplexes into the proper ARGONAUTE (AGO) protein is a key step to forming a functional silencing complex. In Arabidopsis thaliana, the specific loading of miR166/165 into AGO10 (AtAGO10) is critical for the maintenance of the shoot apical meristem, the source of all shoot organs, but the mechanism by which AtAGO10 distinguishes miR166/165 from other cellular miRNAs is not known. Here, we show purified AtAGO10 alone lacks loading selectivity towards miR166/165 duplexes. However, phosphate and HSP chaperone systems reshape the selectivity of AtAGO10 to its physiological substrates. A loop in the AtAGO10 central cleft is essential for recognizing specific mismatches opposite the guide strand 3' region in miR166/165 duplexes. Replacing this loop with the equivalent loop from Homo sapiens AGO2 (HsAGO2) changes AtAGO10 miRNA loading behavior such that 3' region mismatches are ignored and mismatches opposite the guide 5' end instead drive loading, as in HsAGO2. Thus, this study uncovers the molecular mechanism underlying the miR166/165 selectivity of AtAGO10, essential for plant development, and provides new insights into how miRNA duplex structures are recognized for sRNA sorting.

Alternative splicing (AS) is necessary for viral proliferation in host cells and a critical regulatory component of viral gene expression. Conventional RNA-seq approaches provide incomplete coverage of AS due to their short read lengths and are susceptible to biases and artifacts introduced in prevailing library preparation methodologies. Moreover, viral splicing studies are often conducted separately from host cell transcriptome analysis, precluding an assessment of the viral manipulation of host splicing machinery. To address current limitations, we developed a quantitative full-length direct cDNA sequencing strategy to simultaneously profile viral and host cell transcripts. This nanopore-based approach couples processive reverse transcriptases with a novel one-step chemical ablation of 3' RNA ends (termed CASPR), which decreases ribosomal RNA reads and enriches polyadenylated coding sequences. We extensively validate our approach using synthetic reference transcripts and show that CASPR doubles the breadth of coverage per transcript and increases detection of long transcripts (>4 kb), while being functionally equivalent to PolyA+ selection for transcript quantification. We used our approach to interrogate host cell and HIV-1 transcript dynamics during viral reactivation and identified novel putative HIV-1 host factors containing exon skipping or novel intron retentions and delineated the HIV-1 transcriptional state associated with these differentially regulated host factors.

SERPINA1 mRNAs encode the protease inhibitor α-1-antitrypsin and are regulated through post-transcriptional mechanisms. α-1-antitrypsin deficiency leads to chronic obstructive pulmonary disease (COPD) and liver cirrhosis, and specific variants in the 5'-untranslated region (5'-UTR) are associated with COPD. The NM_000295.4 transcript is well expressed and translated in lung and blood and features an extended 5'-UTR that does not contain a competing upstream open reading frame (uORF). We show that the 5'-UTR of NM_000295.4 folds into a well-defined multi-helix structural domain. We systematically destabilized mRNA structure across the NM_000295.4 5'-UTR, and measured changes in (SHAPE quantified) RNA structure and cap-dependent translation relative to a native-sequence reporter. Surprisingly, despite destabilizing local RNA structure, most mutations either had no effect on or decreased translation. Most structure-destabilizing mutations retained native, global 5'-UTR structure. However, those mutations that disrupted the helix that anchors the 5'-UTR domain yielded three groups of non-native structures. Two of these non-native structure groups refolded to create a stable helix near the translation initiation site that decreases translation. Thus, in contrast to the conventional model that RNA structure in 5'-UTRs primarily inhibits translation, complex folding of the NM_000295.4 5'-UTR creates a translation-optimized message by promoting accessibility at the translation initiation site.

DNA lesions in S phase threaten genome stability. The DNA damage tolerance (DDT) pathways overcome these obstacles and allow completion of DNA synthesis by the use of specialised translesion (TLS) DNA polymerases or through recombination-related processes. However, how these mechanisms coordinate with each other and with bulk replication remains elusive. To address these issues, we monitored the variation of replication intermediate architecture in response to ultraviolet irradiation using transmission electron microscopy. We show that the TLS polymerase η, able to accurately bypass the major UV lesion and mutated in the skin cancer-prone xeroderma pigmentosum variant (XPV) syndrome, acts at the replication fork to resolve uncoupling and prevent post-replicative gap accumulation. Repriming occurs as a compensatory mechanism when this on-the-fly mechanism cannot operate, and is therefore predominant in XPV cells. Interestingly, our data support a recombination-independent function of RAD51 at the replication fork to sustain repriming. Finally, we provide evidence for the post-replicative commitment of recombination in gap repair and for pioneering observations of in vivo recombination intermediates. Altogether, we propose a chronology of UV damage tolerance in human cells that highlights the key role of polη in shaping this response and ensuring the continuity of DNA synthesis.

tRNAs that are transcribed in the nucleus are exported to the cytoplasm to perform their iterative essential function in translation. However, the complex set of tRNA post-transcriptional processing and subcellular trafficking steps are not completely understood. In particular, proteins involved in tRNA nuclear export remain unknown since the canonical tRNA nuclear exportin, Los1/Exportin-t, is unessential in all tested organisms. We previously reported that budding yeast Mex67-Mtr2, a mRNA nuclear exporter, co-functions with Los1 in tRNA nuclear export. Here we employed in vivo co-purification of tRNAs with endogenously expressed nuclear exporters to document that Crm1 also is a bona fide tRNA nuclear exporter. We document that Los1, Mex67-Mtr2 and Crm1 possess individual tRNA preferences for forming nuclear export complexes with members of the 10 families of intron-containing pre-tRNAs. Remarkably, Mex67-Mtr2, but not Los1 or Crm1, is error-prone, delivering tRNAs to the cytoplasm prior to 5' leader removal. tRNA retrograde nuclear import functions to monitor the aberrant leader-containing spliced tRNAs, returning them to the nucleus where they are degraded by 3' to 5' exonucleases. Overall, our work identifies a new tRNA nuclear exporter, uncovers exporter preferences for specific tRNA families, and documents contribution of tRNA nuclear import to tRNA quality control.

Development of floral organs exhibits complex molecular mechanisms involving the co-regulation of many genes specialized and precisely functioning in various tissues and developing stages. Advance in spatial transcriptome technologies allows for quantitative measurement of spatially localized gene abundance making it possible to bridge complex scenario of flower organogenesis with genome-wide molecular phenotypes. Here, we apply the 10× Visium technology in the study of the formation of floral organs through development in an orchid plant, Phalaenopsis Big Chili. Cell-types of early floral development including inflorescence meristems, primordia of floral organs and identity determined tissues, are recognized based on spatial expression distribution of thousands of genes in high resolution. In addition, meristematic cells on the basal position of floral organs are found to continuously function in multiple developmental stages after organ initiation. Particularly, the development of anther, which primordium starts from a single spot to multiple differentiated cell-types in later stages including pollinium and other vegetative tissues, is revealed by well-known MADS-box genes and many other downstream regulators. The spatial transcriptome analyses provide comprehensive information of gene activity for understanding the molecular architecture of flower organogenesis and for future genomic and genetic studies of specific cell-types.

Trypanosomes diverged from the main eukaryotic lineage about 600 million years ago, and display some unusual genomic and epigenetic properties that provide valuable insight into the early processes employed by eukaryotic ancestors to regulate chromatin-mediated functions. We analysed Trypanosoma brucei core histones by high mass accuracy middle-down mass spectrometry to map core histone post-translational modifications (PTMs) and elucidate cis-histone combinatorial PTMs (cPTMs). T. brucei histones are heavily modified and display intricate cPTMs patterns, with numerous hypermodified cPTMs that could contribute to the formation of non-repressive euchromatic states. The Trypanosoma brucei H2A C-terminal tail is hyperacetylated, containing up to five acetylated lysine residues. MNase-ChIP-seq revealed a striking enrichment of hyperacetylated H2A at Pol II transcription start regions, and showed that H2A histones that are hyperacetylated in different combinations localised to different genomic regions, suggesting distinct epigenetic functions. Our genomics and proteomics data provide insight into the complex epigenetic mechanisms used by this parasite to regulate a genome that lacks the transcriptional control mechanisms found in later-branched eukaryotes. The findings further demonstrate the complexity of epigenetic mechanisms that were probably shared with the last eukaryotic common ancestor.

Each of the three similar RNA Editing Catalytic Complexes (RECCs) that perform gRNA-directed uridine insertion and deletion during Trypanosoma brucei mitochondrial (mt) mRNA editing has a distinct endonuclease activity that requires two related RNase III proteins, with only one competent for catalysis. We identified multiple loss-of-function mutations in the RNase III and other motifs of the non-catalytic KREPB6, KREPB7, and KREPB8 components by random mutagenesis and screening. These mutations had various effects on growth, editing, and both the abundances and RECC associations of these RNase III protein pairs in bloodstream form (BF) and procyclic form (PF) cells. Protein structure modelling predicted that the Zinc Finger (ZnF) of each paired RNase III protein contacts RNA positioned at the heterodimeric active site which is flanked by helices of a novel RNase III-Associated Motif (RAM). The results indicate that the protein domains of the non-catalytic subunits function together in RECC integrity, substrate binding, and editing site recognition during the multistep RNA editing process. Additionally, several mutants display distinct functional consequences in different life cycle stages. These results highlight the complementary roles of protein pairs and three RECCs within the complicated T. brucei mRNA editing machinery that matures mt mRNAs differentially between developmental stages.

Adenine base editors (ABEs) catalyze A-to-G conversions, offering therapeutic options to treat the major class of human pathogenic single nucleotide polymorphisms (SNPs). However, robust and precise editing at diverse genome loci remains challenging. Here, using high-throughput chemical screening, we identified and validated SB505124, a selective ALK5 inhibitor, as an ABE activator. Treating cells with SB505124 enhanced on-target editing at multiple genome loci, including epigenetically refractory regions, and showed little effect on off-target conversion on the genome. Furthermore, SB505124 facilitated the editing of disease-associated genes in vitro and in vivo. Intriguingly, SB505124 served as a specific activator by selectively promoting ABE activity. Mechanistically, SB505124 promotes ABE editing, at least in part, by enhancing ABE expression and modulating DNA repair-associated genes. Our findings reveal the role of the canonical transforming growth factor-β pathway in gene editing and equip ABEs with precise chemical control.

Interactions between epigenetic readers and histone modifications play a pivotal role in gene expression regulation and aberrations can enact etiopathogenic roles in both developmental and acquired disorders like cancer. Typically, epigenetic interactions are studied by mass spectrometry or chromatin immunoprecipitation sequencing. However, in these methods, spatial information is completely lost. Here, we devise an expansion microscopy based method, termed Expansion Microscopy for Epigenetics or ExEpi, to preserve spatial information and improve resolution. We calculated relative co-localization ratios for two epigenetic readers, lens epithelium derived growth factor (LEDGF) and bromodomain containing protein 4 (BRD4), with marks for heterochromatin (H3K9me3 and H3K27me3) and euchromatin (H3K36me2, H3K36me3 and H3K9/14ac). ExEpi confirmed their preferred epigenetic interactions, showing co-localization for LEDGF with H3K36me3/me2 and for BRD4 with H3K9/14ac. Moreover addition of JQ1, a known BET-inhibitor, abolished BRD4 interaction with H3K9/14ac with an IC50 of 137 nM, indicating ExEpi could serve as a platform for epigenetic drug discovery. Since ExEpi retains spatial information, the nuclear localization of marks and readers was determined, which is one of the main advantages of ExEpi. The heterochromatin mark, H3K9me3, is located in the nuclear rim whereas LEDGF co-localization with H3K36me3 and BRD4 co-localization with H3K9/14ac occur further inside the nucleus.

De novo mutations (DNMs) are an important cause of genetic disorders. The accurate identification of DNMs from sequencing data is therefore fundamental to rare disease research and diagnostics. Unfortunately, identifying reliable DNMs remains a major challenge due to sequence errors, uneven coverage, and mapping artifacts. Here, we developed a deep convolutional neural network (CNN) DNM caller (DeNovoCNN), that encodes the alignment of sequence reads for a trio as 160$ \times$164 resolution images. DeNovoCNN was trained on DNMs of 5616 whole exome sequencing (WES) trios achieving total 96.74% recall and 96.55% precision on the test dataset. We find that DeNovoCNN has increased recall/sensitivity and precision compared to existing DNM calling approaches (GATK, DeNovoGear, DeepTrio, Samtools) based on the Genome in a Bottle reference dataset and independent WES and WGS trios. Validations of DNMs based on Sanger and PacBio HiFi sequencing confirm that DeNovoCNN outperforms existing methods. Most importantly, our results suggest that DeNovoCNN is likely robust against different exome sequencing and analyses approaches, thereby allowing the application on other datasets. DeNovoCNN is freely available as a Docker container and can be run on existing alignment (BAM/CRAM) and variant calling (VCF) files from WES and WGS without a need for variant recalling.

RNA molecules harbor diverse modifications that play important regulatory roles in a variety of biological processes. Over 150 modifications have been identified in RNA molecules. N6-methyladenosine (m6A) and 1-methyladenosine (m1A) are prevalent modifications occurring in various RNA species of mammals. Apart from the single methylation of adenosine (m6A and m1A), dual methylation modification occurring in the nucleobase of adenosine, such as N6,N6-dimethyladenosine (m6,6A), also has been reported to be present in RNA of mammals. Whether there are other forms of dual methylation modification occurring in the nucleobase of adenosine other than m6,6A remains elusive. Here, we reported the existence of a novel adenosine dual methylation modification, i.e. 1,N6-dimethyladenosine (m1,6A), in tRNAs of living organisms. We confirmed that m1,6A is located at position 58 of tRNAs and is prevalent in mammalian cells and tissues. The measured level of m1,6A ranged from 0.0049% to 0.047% in tRNAs. Furthermore, we demonstrated that TRMT6/61A could catalyze the formation of m1,6A in tRNAs and m1,6A could be demethylated by ALKBH3. Collectively, the discovery of m1,6A expands the diversity of RNA modifications and may elicit a new tRNA modification-mediated gene regulation pathway.

Autogenous interactions between mRNAs and the proteins they encode are implicated in cellular feedback-loop regulation, but their extent and mechanistic foundation are unclear. It was recently hypothesized that such interactions may be common, reflecting the role of intrinsic nucleobase-amino acid affinities in shaping the genetic code's structure. Here we analyze a comprehensive set of CLIP-seq experiments involving multiple protocols and report on widespread autogenous interactions across different organisms. Specifically, 230 of 341 (67%) studied RNA-binding proteins (RBPs) interact with their own mRNAs, with a heavy enrichment among high-confidence hits and a preference for coding sequence binding. We account for different confounding variables, including physical (overexpression and proximity during translation), methodological (difference in CLIP protocols, peak callers and cell types) and statistical (treatment of null backgrounds). In particular, we demonstrate a high statistical significance of autogenous interactions by sampling null distributions of fixed-margin interaction matrices. Furthermore, we study the dependence of autogenous binding on the presence of RNA-binding motifs and structured domains in RBPs. Finally, we show that intrinsic nucleobase-amino acid affinities favor co-aligned binding between mRNA coding regions and the proteins they encode. Our results suggest a central role for autogenous interactions in RBP regulation and support the possibility of a fundamental connection between coding and binding.

RNA editing processes are strikingly different in animals and plants. Up to thousands of specific cytidines are converted into uridines in plant chloroplasts and mitochondria whereas up to millions of adenosines are converted into inosines in animal nucleo-cytosolic RNAs. It is unknown whether these two different RNA editing machineries are mutually incompatible. RNA-binding pentatricopeptide repeat (PPR) proteins are the key factors of plant organelle cytidine-to-uridine RNA editing. The complete absence of PPR mediated editing of cytosolic RNAs might be due to a yet unknown barrier that prevents its activity in the cytosol. Here, we transferred two plant mitochondrial PPR-type editing factors into human cell lines to explore whether they could operate in the nucleo-cytosolic environment. PPR56 and PPR65 not only faithfully edited their native, co-transcribed targets but also different sets of off-targets in the human background transcriptome. More than 900 of such off-targets with editing efficiencies up to 91%, largely explained by known PPR-RNA binding properties, were identified for PPR56. Engineering two crucial amino acid positions in its PPR array led to predictable shifts in target recognition. We conclude that plant PPR editing factors can operate in the entirely different genetic environment of the human nucleo-cytosol and can be intentionally re-engineered towards new targets.

Hybridization and strand displacement kinetics determine the evolution of the base paired configurations of mixtures of oligonucleotides over time. Although much attention has been focused on the thermodynamics of DNA and RNA base pairing in the scientific literature, much less work has been done on the time dependence of interactions involving multiple strands, especially in RNA. Here we provide a study of oligoribonucleotide interaction kinetics and show that it is possible to calculate the association, dissociation and strand displacement rates displayed by short oligonucleotides (5nt-12nt) that exhibit no expected secondary structure as simple functions of oligonucleotide length, CG content, ΔG of hybridization and ΔG of toehold binding. We then show that the resultant calculated kinetic parameters are consistent with the experimentally observed time dependent changes in concentrations of the different species present in mixtures of multiple competing RNA strands. We show that by changing the mixture composition, it is possible to create and tune kinetic traps that extend by orders of magnitude the typical sub-second hybridization timescale of two complementary oligonucleotides. We suggest that the slow equilibration of complex oligonucleotide mixtures may have facilitated the nonenzymatic replication of RNA during the origin of life.

Natural methylotrophs are attractive methanol utilization hosts, but lack flexible expression tools. In this study, we developed yeast transcriptional device libraries for precise synthesis of value-added chemicals from methanol. We synthesized transcriptional devices by fusing bacterial DNA-binding proteins (DBPs) with yeast transactivation domains, and linking bacterial binding sequences (BSs) with the yeast core promoter. Three DBP-BS pairs showed good activity when working with transactivation domains and the core promoter of PAOX1 in the methylotrophic yeast, Pichia pastoris. Fine-tuning of the tandem BSs, spacers and differentiated input promoters further enabled a constitutive transcriptional device library (cTRDL) composed of 126 transcriptional devices with an expression strength of 16-520% and an inducible TRDL (iTRDL) composed of 162 methanol-inducible transcriptional devices with an expression strength of 30-500%, compared with PAOX1. Selected devices from iTRDL were adapted to the dihydromonacolin L biosynthetic pathway by orthogonal experimental design, reaching 5.5-fold the production from the PAOX1-driven pathway. The full factorial design of the selected devices from the cTRDL was adapted to the downstream pathway of dihydromonacolin L to monacolin J. Monacolin J production from methanol reached 3.0-fold the production from the PAOX1-driven pathway. Our engineered toolsets ensured multilevel pathway control of chemical synthesis in methylotrophic yeasts.

Intron diversity facilitates regulated gene expression and alternative splicing. Spliceosomes excise introns after recognizing their splicing signals: the 5'-splice site (5'ss), branchpoint (BP) and 3'-splice site (3'ss). The latter two signals are recognized by U2 small nuclear ribonucleoprotein (snRNP) and its accessory factors (U2AFs), but longer spacings between them result in weaker splicing. Here, we show that excision of introns with a BP-distant 3'ss (e.g. rap1 intron 2) requires the ubiquitin-fold-activated splicing regulator Sde2 in Schizosaccharomyces pombe. By monitoring splicing-specific ura4 reporters in a collection of S. pombe mutants, Cay1 and Tls1 were identified as additional regulators of this process. The role of Sde2, Cay1 and Tls1 was further confirmed by increasing BP-3'ss spacings in a canonical tho5 intron. We also examined BP-distant exons spliced independently of these factors and observed that RNA secondary structures possibly bridged the gap between the two signals. These proteins may guide the 3'ss towards the spliceosome's catalytic centre by folding the RNA between the BP and 3'ss. Orthologues of Sde2, Cay1 and Tls1, although missing in the intron-poor Saccharomyces cerevisiae, are present in intron-rich eukaryotes, including humans. This type of intron-specific pre-mRNA splicing appears to have evolved for regulated gene expression and alternative splicing of key heterochromatin factors.

Nucleus-mitochondria crosstalk is essential for cellular and organismal homeostasis. Although anterograde (nucleus-to-mitochondria) pathways have been well characterized, retrograde (mitochondria-to-nucleus) pathways remain to be clarified. Here, we found that mitochondrial dysfunction triggered a retrograde signaling via unique transcriptional and chromatin factors in hepatic cells. Our transcriptomic analysis revealed that the loss of mitochondrial transcription factor A led to mitochondrial dysfunction and dramatically induced expression of amphiregulin (AREG) and other secretory protein genes. AREG expression was also induced by various mitochondria stressors and was upregulated in murine liver injury models, suggesting that AREG expression is a hallmark of mitochondrial damage. Using epigenomic and informatic approaches, we identified that mitochondrial dysfunction-responsive enhancers of AREG gene were activated by c-JUN/YAP1/TEAD axis and were repressed by chromatin remodeler BRG1. Furthermore, while mitochondrial dysfunction-activated enhancers were enriched with JUN and TEAD binding motifs, the repressed enhancers possessed the binding motifs for hepatocyte nuclear factor 4α, suggesting that both stress responsible and cell type-specific enhancers were reprogrammed. Our study revealed that c-JUN and YAP1-mediated enhancer activation shapes the mitochondrial stress-responsive phenotype, which may shift from metabolism to stress adaptation including protein secretion under such stressed conditions.

The reversible post-translational modification (PTM) of proteins plays an important role in many cellular processes. Lysine crotonylation (Kcr) is a newly identified PTM, but its functional significance remains unclear. Here, we found that Kcr is involved in the replication stress response. We show that crotonylation of histone H2A at lysine 119 (H2AK119) and ubiquitination of H2AK119 are reversibly regulated by replication stress. Decrotonylation of H2AK119 by SIRT1 is a prerequisite for subsequent ubiquitination of H2AK119 by BMI1. Accumulation of ubiquitinated H2AK119 at reversed replication forks leads to the release of RNA Polymerase II and transcription repression in the vicinity of stalled replication forks. These effects attenuate transcription-replication conflicts (TRCs) and TRC-associated R-loop formation and DNA double-strand breaks. These findings suggest that decrotonylation and ubiquitination of H2A at lysine 119 act together to resolve replication stress-induced TRCs and protect genome stability.

Neurospora crassa protein QDE-1, a member of the two-barrel polymerase superfamily, possesses both DNA- and RNA-dependent RNA polymerase (DdRP and RdRP) activities. The dual activities are essential for the production of double-stranded RNAs (dsRNAs), the precursors of small interfering RNAs (siRNAs) in N. crassa. Here, we report five complex structures of N-terminal truncated QDE-1 (QDE-1ΔN), representing four different reaction states: DNA/RNA-templated elongation, the de novo initiation of RNA synthesis, the first step of nucleotide condensation during de novo initiation and initial NTP loading. The template strand is aligned by a bridge-helix and double-psi beta-barrels 2 (DPBB2), the RNA product is held by DPBB1 and the slab domain. The DNA template unpairs with the RNA product at position -7, but the RNA template remains paired. The NTP analog coordinates with cations and is precisely positioned at the addition site by a rigid trigger loop and a proline-containing loop in the active center. The unique C-terminal tail from the QDE-1 dimer partner inserts into the substrate-binding cleft and plays regulatory roles in RNA synthesis. Collectively, this work elucidates the conserved mechanisms for DNA/RNA-dependent dual activities by QDE-1 and other two-barrel polymerase superfamily members.

Retrograde bone morphogenetic protein (BMP) signaling at the Drosophila neuromuscular junction (NMJ) has served as a paradigm to study TGF-β-dependent synaptic function and maturation. Yet, how retrograde BMP signaling transcriptionally regulates these functions remains unresolved. Here, we uncover a gene network, enriched for neurotransmission-related genes, that is controlled by retrograde BMP signaling in motor neurons through two Smad-binding cis-regulatory motifs, the BMP-activating (BMP-AE) and silencer (BMP-SE) elements. Unpredictably, both motifs mediate direct gene activation, with no involvement of the BMP derepression pathway regulators Schnurri and Brinker. Genome editing of candidate BMP-SE and BMP-AE within the locus of the active zone gene bruchpilot, and a novel Ly6 gene witty, demonstrated the role of these motifs in upregulating genes required for the maturation of pre- and post-synaptic NMJ compartments. Our findings uncover how Smad-dependent transcriptional mechanisms specific to motor neurons directly orchestrate a gene network required for synaptic maturation by retrograde BMP signaling.

Next-generation sequencing-based metagenomics has enabled to identify microorganisms in characteristic habitats without the need for lengthy cultivation. Importantly, clinically relevant phenomena such as resistance to medication, virulence or interactions with the environment can vary already within species. Therefore, a major current challenge is to reconstruct individual genomes from the sequencing reads at the level of strains, and not just the level of species. However, strains of one species can differ only by minor amounts of variants, which makes it difficult to distinguish them. Despite considerable recent progress, related approaches have remained fragmentary so far. Here, we present StrainXpress, as a comprehensive solution to the problem of strain aware metagenome assembly from next-generation sequencing reads. In experiments, StrainXpress reconstructs strain-specific genomes from metagenomes that involve up to >1000 strains and proves to successfully deal with poorly covered strains. The amount of reconstructed strain-specific sequence exceeds that of the current state-of-the-art approaches by on average 26.75% across all data sets (first quartile: 18.51%, median: 26.60%, third quartile: 35.05%).

Arrayed libraries of defined mutants have been used to elucidate gene function in the post-genomic era. Yeast haploid gene deletion libraries have pioneered this effort, but are costly to construct, do not reveal phenotypes that may occur with partial gene function and lack essential genes required for growth. We therefore devised an efficient method to construct a library of barcoded insertion mutants with a wider range of phenotypes that can be generalized to other organisms or collections of DNA samples. We developed a novel but simple three-dimensional pooling and multiplexed sequencing approach that leveraged sequence information to reduce the number of required sequencing reactions by orders of magnitude, and were able to identify the barcode sequences and DNA insertion sites of 4391 Schizosaccharomyces pombe insertion mutations with only 40 sequencing preparations. The insertion mutations are in the genes and untranslated regions of nonessential, essential and noncoding RNA genes, and produced a wider range of phenotypes compared to the cognate deletion mutants, including novel phenotypes. This mutant library represents both a proof of principle for an efficient method to produce novel mutant libraries and a valuable resource for the S. pombe research community.

Queuosine (Q) is a conserved tRNA modification at the wobble anticodon position of tRNAs that read the codons of amino acids Tyr, His, Asn, and Asp. Q-modification in tRNA plays important roles in the regulation of translation efficiency and fidelity. Queuosine tRNA modification is synthesized de novo in bacteria, whereas in mammals the substrate for Q-modification in tRNA is queuine, the catabolic product of the Q-base of gut bacteria. This gut microbiome dependent tRNA modification may play pivotal roles in translational regulation in different cellular contexts, but extensive studies of Q-modification biology are hindered by the lack of high throughput sequencing methods for its detection and quantitation. Here, we describe a periodate-treatment method that enables single base resolution profiling of Q-modification in tRNAs by Nextgen sequencing from biological RNA samples. Periodate oxidizes the Q-base, which results in specific deletion signatures in the RNA-seq data. Unexpectedly, we found that periodate-treatment also enables the detection of several 2-thio-modifications including τm5s2U, mcm5s2U, cmnm5s2U, and s2C by sequencing in human and E. coli tRNA. We term this method periodate-dependent analysis of queuosine and sulfur modification sequencing (PAQS-seq). We assess Q- and 2-thio-modifications at the tRNA isodecoder level, and 2-thio modification changes in stress response. PAQS-seq should be widely applicable in the biological studies of Q- and 2-thio-modifications in mammalian and microbial tRNAs.

Eukaryotic chromosomes contain regions of varying accessibility, yet DNA replication factors must access all regions. The first replication step is loading MCM complexes to license replication origins during the G1 cell cycle phase. It is not yet known how mammalian MCM complexes are adequately distributed to both accessible euchromatin regions and less accessible heterochromatin regions. To address this question, we combined time-lapse live-cell imaging with immunofluorescence imaging of single human cells to quantify the relative rates of MCM loading in euchromatin and heterochromatin throughout G1. We report here that MCM loading in euchromatin is faster than that in heterochromatin in early G1, but surprisingly, heterochromatin loading accelerates relative to euchromatin loading in middle and late G1. This differential acceleration allows both chromatin types to begin S phase with similar concentrations of loaded MCM. The different loading dynamics require ORCA-dependent differences in origin recognition complex distribution. A consequence of heterochromatin licensing dynamics is that cells experiencing a truncated G1 phase from premature cyclin E expression enter S phase with underlicensed heterochromatin, and DNA damage accumulates preferentially in heterochromatin in the subsequent S/G2 phase. Thus, G1 length is critical for sufficient MCM loading, particularly in heterochromatin, to ensure complete genome duplication and to maintain genome stability.

Ewing sarcoma is a prototypical fusion transcription factor-associated pediatric cancer that expresses EWS/FLI or a highly related FET/ETS chimera. EWS/FLI dysregulates transcription to induce and maintain sarcomagenesis, but the mechanisms utilized are not fully understood. We therefore sought to define the global effects of EWS/FLI on chromatin conformation and transcription in Ewing sarcoma cells using a well-validated 'knock-down/rescue' model of EWS/FLI function in combination with next generation sequencing assays to evaluate how the chromatin landscape changes with loss, and recovery, of EWS/FLI expression. We found that EWS/FLI (and EWS/ERG) genomic localization is largely conserved across multiple patient-derived Ewing sarcoma cell lines. This EWS/FLI binding signature is associated with establishment of topologically-associated domain (TAD) boundaries, compartment activation, enhancer-promoter looping that involve both intra- and inter-TAD interactions, and gene activation. In addition, EWS/FLI co-localizes with the loop-extrusion factor cohesin to promote chromatin loops and TAD boundaries. Importantly, local chromatin features provide the basis for transcriptional heterogeneity in regulation of direct EWS/FLI target genes across different Ewing sarcoma cell lines. These data demonstrate a key role of EWS/FLI in mediating genome-wide changes in chromatin configuration and support the notion that fusion transcription factors serve as master regulators of three-dimensional reprogramming of chromatin.

Cells respond to double-strand breaks (DSBs) by activating DNA damage response pathways, including cell cycle arrest. We have previously shown that a single double-strand break generated via CRISPR/Cas9 is sufficient to delay cell cycle progression and compromise cell viability. However, we also found that the cellular response to DSBs can vary, independent of the number of lesions. This implies that not all DSBs are equally toxic, and raises the question if the location of a single double-strand break could influence its toxicity. To systematically investigate if DSB-location is a determinant of toxicity we performed a CRISPR/Cas9 screen targeting 6237 single sites in the human genome. Next, we developed a data-driven framework to design CRISPR/Cas9 sgRNA (crRNA) pools targeting specific chromatin features. The chromatin context was defined using ChromHMM states, Lamin-B1 DAM-iD, DNAseI hypersensitivity, and RNA-sequencing data. We computationally designed 6 distinct crRNA pools, each containing 10 crRNAs targeting the same chromatin state. We show that the toxicity of a DSB is highly similar across the different ChromHMM states. Rather, we find that the major determinants of toxicity of a sgRNA are cutting efficiency and off-target effects. Thus, chromatin features have little to no effect on the toxicity of a single CRISPR/Cas9-induced DSB.

During DNA replication, DNA lesions in lagging strand templates are initially encountered by DNA polymerase δ (pol δ) holoenzymes comprised of pol δ and the PCNA processivity sliding clamp. These encounters are thought to stall replication of an afflicted template before the lesion, activating DNA damage tolerance (DDT) pathways that replicate the lesion and adjacent DNA sequence, allowing pol δ to resume. However, qualitative studies observed that human pol δ can replicate various DNA lesions, albeit with unknown proficiencies, which raises issues regarding the role of DDT in replicating DNA lesions. To address these issues, we re-constituted human lagging strand replication to quantitatively characterize initial encounters of pol δ holoenzymes with DNA lesions. The results indicate pol δ holoenzymes support dNTP incorporation opposite and beyond multiple lesions and the extent of these activities depends on the lesion and pol δ proofreading. Furthermore, after encountering a given DNA lesion, subsequent dissociation of pol δ is distributed around the lesion and a portion does not dissociate. The distributions of these events are dependent on the lesion and pol δ proofreading. Collectively, these results reveal complexity and heterogeneity in the replication of lagging strand DNA lesions, significantly advancing our understanding of human DDT.

Erratum in Nucleic Acids Res. 2022 Nov 11;50(20):11992.

Abortive infection (Abi) is a bacterial antiphage defense strategy involving suicide of the infected cell. Some Abi pathways involve polymerases that are related to reverse transcriptases. They are unique in the way they combine the ability to synthesize DNA in a template-independent manner with protein priming. Here, we report crystal and cryo-electron microscopy structures of two Abi polymerases: AbiK and Abi-P2. Both proteins adopt a bilobal structure with an RT-like domain that comprises palm and fingers subdomains and a unique helical domain. AbiK and Abi-P2 adopt a hexameric and trimeric configuration, respectively, which is unprecedented for reverse transcriptases. Biochemical experiments showed that the formation of these oligomers is required for the DNA polymerization activity. The structure of the AbiK-DNA covalent adduct visualized interactions between the 3' end of DNA and the active site and covalent attachment of the 5' end of DNA to a tyrosine residue used for protein priming. Our data reveal a structural basis of the mechanism of highly unusual template-independent protein-priming polymerases.

Pseudomonas aeruginosa is an opportunistic pathogen and major cause of hospital-acquired infections. The virulence of P. aeruginosa is largely determined by its transcriptional regulatory network (TRN). We used 411 transcription profiles of P. aeruginosa from diverse growth conditions to construct a quantitative TRN by identifying independently modulated sets of genes (called iModulons) and their condition-specific activity levels. The current study focused on the use of iModulons to analyze the biofilm production and antibiotic resistance of P. aeruginosa. Our analysis revealed: (i) 116 iModulons, 81 of which show strong association with known regulators; (ii) novel roles of regulators in modulating antibiotics efflux pumps; (iii) substrate-efflux pump associations; (iv) differential iModulon activity in response to beta-lactam antibiotics in bacteriological and physiological media; (v) differential activation of 'Cell Division' iModulon resulting from exposure to different beta-lactam antibiotics and (vi) a role of the PprB iModulon in the stress-induced transition from planktonic to biofilm lifestyle. In light of these results, the construction of an iModulon-based TRN provides a transcriptional regulatory basis for key aspects of P. aeruginosa infection, such as antibiotic stress responses and biofilm formation. Taken together, our results offer a novel mechanistic understanding of P. aeruginosa virulence.

The arrangement of nucleosomes inside chromatin is of extensive interest. While in vitro experiments have revealed the formation of 30 nm fibers, most in vivo studies have failed to confirm their presence in cell nuclei. To reconcile the diverging experimental findings, we characterized chromatin organization using a residue-level coarse-grained model. The computed force-extension curve matches well with measurements from single-molecule experiments. Notably, we found that a dodeca-nucleosome in the two-helix zigzag conformation breaks into structures with nucleosome clutches and a mix of trimers and tetramers under tension. Such unfolded configurations can also be stabilized through trans interactions with other chromatin chains. Our study suggests that unfolding from chromatin fibers could contribute to the irregularity of in vivo chromatin configurations. We further revealed that chromatin segments with fibril or clutch structures engaged in distinct binding modes and discussed the implications of these inter-chain interactions for a potential sol-gel phase transition.

The introduction of phosphorothioate (PS) linkages to the backbone of therapeutic nucleic acids substantially increases their stability and potency. It also affects their interactions with cellular proteins, but the molecular mechanisms that underlie this effect are poorly understood. Here, we report structural and biochemical studies of interactions between annexin A2, a protein that does not possess any known canonical DNA binding domains, and phosphorothioate-modified antisense oligonucleotides. We show that a unique mode of hydrophobic interactions between a sulfur atom of the phosphorothioate group and lysine and arginine residues account for the enhanced affinity of modified nucleic acid for the protein. Our results demonstrate that this mechanism of interaction is observed not only for nucleic acid-binding proteins but can also account for the association of PS oligonucleotides with other proteins. Using the anomalous diffraction of sulfur, we showed that preference for phosphorothioate stereoisomers is determined by the hydrophobic environment around the PS linkage that comes not only from protein but also from additional structural features within the ASO such as 5-Me groups on cytosine nucleobases.

Site-specific DNA recombinases play a variety of biological roles, often related to the dissemination of antibiotic resistance, and are also useful synthetic biology tools. The simplest site-specific recombination systems will recombine any two cognate sites regardless of context. Other systems have evolved elaborate mechanisms, often sensing DNA topology, to ensure that only one of multiple possible recombination products is produced. The closely related resolvases from the Tn3 and γδ transposons have historically served as paradigms for the regulation of recombinase activity by DNA topology. However, despite many proposals, models of the multi-subunit protein-DNA complex (termed the synaptosome) that enforces this regulation have been unsatisfying due to a lack of experimental constraints and incomplete concordance with experimental data. Here, we present new structural and biochemical data that lead to a new, detailed model of the Tn3 synaptosome, and discuss how it harnesses DNA topology to regulate the enzymatic activity of the recombinase.