Tyrosine-type site-specific recombinases (Y-SSRs) are versatile tools for genome engineering due to their ability to mediate excision, integration, inversion and exchange of genomic DNA with single nucleotide precision. The ever-increasing need for sophisticated genome engineering is driving efforts to identify novel SSR systems with intrinsic properties more suitable for particular applications. In this work, we develop a systematic computational workflow for annotation of putative Y-SSR systems and apply this pipeline to identify and characterize eight new naturally occurring Cre-type SSR systems. We test their activity in bacterial and mammalian cells and establish selectivity profiles for the new and already established Cre-type SSRs with regard to their ability to mutually recombine their target sites. These data form the basis for sophisticated genome engineering experiments using combinations of Y-SSRs in research fields including advanced genomics and synthetic biology. Finally, we identify putative pseudo-sites and potential off-targets for Y-SSRs in the human and mouse genome. Together with established methods for altering the DNA-binding specificity of this class of enzymes, this work should facilitate the use of Y-SSRs for future genome surgery applications.

The N-(2-deoxy-d-erythro-pentofuranosyl)-urea DNA lesion forms following hydrolytic fragmentation of cis-5R,6S- and trans-5R,6R-dihydroxy-5,6-dihydrothymidine (thymine glycol, Tg) or from oxidation of 7,8-dihydro-8-oxo-deoxyguanosine (8-oxodG) and subsequent hydrolysis. It interconverts between α and β deoxyribose anomers. Synthetic oligodeoxynucleotides containing this adduct are efficiently incised by unedited (K242) and edited (R242) forms of the hNEIL1 glycosylase. The structure of a complex between the active site unedited mutant CΔ100 P2G hNEIL1 (K242) glycosylase and double-stranded (ds) DNA containing a urea lesion reveals a pre-cleavage intermediate, in which the Gly2 N-terminal amine forms a conjugate with the deoxyribose C1' of the lesion, with the urea moiety remaining intact. This structure supports a proposed catalytic mechanism in which Glu3-mediated protonation of O4' facilitates attack at deoxyribose C1'. The deoxyribose is in the ring-opened configuration with the O4' oxygen protonated. The electron density of Lys242 suggests the 'residue 242-in conformation' associated with catalysis. This complex likely arises because the proton transfer steps involving Glu6 and Lys242 are hindered due to Glu6-mediated H-bonding with the Gly2 and the urea lesion. Consistent with crystallographic data, biochemical analyses show that the CΔ100 P2G hNEIL1 (K242) glycosylase exhibits a residual activity against urea-containing dsDNA.

While global transcription factors (TFs) have been studied extensively in Escherichia coli model strains, conservation and diversity in TF regulation between strains is still unknown. Here we use a combination of ChIP-exo-to define ferric uptake regulator (Fur) binding sites-and differential gene expression-to define the Fur regulon in nine E. coli strains. We then define a pan-regulon consisting of 469 target genes that includes all Fur target genes in all nine strains. The pan-regulon is then divided into the core regulon (target genes found in all the strains, n = 36), the accessory regulon (target found in two to eight strains, n = 158) and the unique regulon (target genes found in one strain, n = 275). Thus, there is a small set of Fur regulated genes common to all nine strains, but a large number of regulatory targets unique to a particular strain. Many of the unique regulatory targets are genes unique to that strain. This first-established pan-regulon reveals a common core of conserved regulatory targets and significant diversity in transcriptional regulation amongst E. coli strains, reflecting diverse niche specification and strain history.

The structure of mRNA molecules plays an important role in its interactions with trans-acting factors, notably RNA binding proteins (RBPs), thus contributing to the functional consequences of this interplay. However, current transcriptome-wide experimental methods to chart these interactions are limited by their poor sensitivity. Here we extend the hiCLIP atlas of duplexes bound by Staufen1 (STAU1) ∼10-fold, through careful consideration of experimental assumptions, and the development of bespoke computational methods which we apply to existing data. We present Tosca, a Nextflow computational pipeline for the processing, analysis and visualisation of proximity ligation sequencing data generally. We use our extended duplex atlas to discover insights into the RNA selectivity of STAU1, revealing the importance of structural symmetry and duplex-span-dependent nucleotide composition. Furthermore, we identify heterogeneity in the relationship between transcripts with STAU1-bound 3' UTR duplexes and metabolism of the associated RNAs that we relate to RNA structure: transcripts with short-range proximal 3' UTR duplexes have high degradation rates, but those with long-range duplexes have low rates. Overall, our work enables the integrative analysis of proximity ligation data delivering insights into specific features and effects of RBP-RNA structure interactions.

DNA modifications are critical in fine-tuning the biological processes in model organisms. However, the presence of cytosine methylation (5mC) and the function of the putative DNA methyltransferase, PfDNMT2, in the human malaria pathogen, Plasmodium falciparum, remain controversial. Here, we revisited the 5mC in the parasite genome and the function of PfDNMT2. Low levels of genomic 5mC (0.1-0.2%) during asexual development were identified using a sensitive mass spectrometry procedure. Native PfDNMT2 displayed substantial DNA methylation activities, and disruption or overexpression of PfDNMT2 resulted in reduced or elevated genomic 5mC levels, respectively. PfDNMT2 disruption led to an increased proliferation phenotype, with the parasites having an extended schizont stage and producing a higher number of progenies. Consistent with PfDNMT2's interaction with an AP2 domain-containing transcription factor, transcriptomic analyses revealed that PfDNMT2 disruption led to a drastic alteration in the expression of many genes, some of which provided the molecular basis of enhanced proliferation after PfDNMT2 disruption. Furthermore, levels of tRNAAsp and its methylation rate at position C38, and the translation of a reporter containing an aspartate repeat were significantly reduced after PfDNMT2 disruption, while the levels of tRNAAsp and its C38 methylation were restored after complementation of PfDNMT2. Our study sheds new light on the dual function of PfDNMT2 during P. falciparum asexual development.

Chondrogenesis is a multistep process, in which cartilage progenitor cells generate a tissue with distinct structural and functional properties. Although several approaches to cartilage regeneration rely on the differentiation of implanted progenitor cells, the temporal transcriptomic landscape of in vitro chondrogenesis in different models has not been reported. Using RNA sequencing, we examined differences in gene expression patterns during cartilage formation in micromass cultures of embryonic limb bud-derived progenitors. Principal component and trajectory analyses revealed a progressively different and distinct transcriptome during chondrogenesis. Differentially expressed genes (DEGs), based on pairwise comparisons of samples from consecutive days were classified into clusters and analysed. We confirmed the involvement of the top DEGs in chondrogenic differentiation using pathway analysis and identified several chondrogenesis-associated transcription factors and collagen subtypes that were not previously linked to cartilage formation. Transient gene silencing of ATOH8 or EBF1 on day 0 attenuated chondrogenesis by deregulating the expression of key osteochondrogenic marker genes in micromass cultures. These results provide detailed insight into the molecular mechanism of chondrogenesis in primary micromass cultures and present a comprehensive dataset of the temporal transcriptomic landscape of chondrogenesis, which may serve as a platform for new molecular approaches in cartilage tissue engineering.

The chemical modification of ribosomal RNA and proteins is critical for ribosome assembly, for protein synthesis and may drive ribosome specialisation in development and disease. However, the inability to accurately visualise these modifications has limited mechanistic understanding of the role of these modifications in ribosome function. Here we report the 2.15 Å resolution cryo-EM reconstruction of the human 40S ribosomal subunit. We directly visualise post-transcriptional modifications within the 18S rRNA and four post-translational modifications of ribosomal proteins. Additionally, we interpret the solvation shells in the core regions of the 40S ribosomal subunit and reveal how potassium and magnesium ions establish both universally conserved and eukaryote-specific coordination to promote the stabilisation and folding of key ribosomal elements. This work provides unprecedented structural details for the human 40S ribosomal subunit that will serve as an important reference for unravelling the functional role of ribosomal RNA modifications.

O-Phosphoseryl-tRNASec selenium transferase (SepSecS) catalyzes the terminal step of selenocysteine (Sec) synthesis in archaea and eukaryotes. How the Sec synthetic machinery recognizes and discriminates tRNASec from the tRNA pool is essential to the integrity of the selenoproteome. Previously, we suggested that SepSecS adopts a competent conformation that is pre-ordered for catalysis. Herein, using high-resolution X-ray crystallography, we visualized tRNA-dependent conformational changes in human SepSecS that may be a prerequisite for achieving catalytic competency. We show that tRNASec binding organizes the active sites of the catalytic protomer, while stabilizing the N- and C-termini of the non-catalytic protomer. Binding of large anions to the catalytic groove may further optimize the catalytic site for substrate binding and catalysis. Our biochemical and mutational analyses demonstrate that productive SepSecS•tRNASec complex formation is enthalpically driven and primarily governed by electrostatic interactions between the acceptor-, TΨC-, and variable arms of tRNASec and helices α1 and α14 of SepSecS. The detailed visualization of the tRNA-dependent activation of SepSecS provides a structural basis for a revised model of the terminal reaction of Sec formation in archaea and eukaryotes.

Natural transformation is one of the major mechanisms of horizontal gene transfer. Although it is usually studied using purified DNA in the laboratory, recent studies showed that many naturally competent bacteria acquired exogenous DNA from neighboring donor cells. Our previous work indicates that cell-to-cell natural transformation (CTCNT) using two different Bacillus subtilis strains is a highly efficient process; however, the mechanism is unclear. In this study, we further characterized CTCNT and mapped the transferred DNA in the recombinants using whole genome sequencing. We found that a recombinant strain generated by CTCNT received up to 66 transferred DNA segments; the average length of acquired continuous DNA stretches was approximately 27 kb with a maximum length of 347 kb. Moreover, up to 1.54 Mb genomic DNA (37% of the chromosome) was transferred from the donors into one recipient cell. These results suggest that B. subtilis CTCNT facilitates horizontal gene transfer by increasing the transfer of DNA segments and fostering the exchange of large continuous genomic regions. This indicates that the potency of bacterial natural transformation is underestimated using traditional approaches and reveals that DNA donor cells may play an important role in the transformation process in natural environments.

CRISPR (clustered regularly interspaced short palindromic repeats) technology holds tremendous promise for gene regulation and editing. However, precise control of CRISPR editing is essential to overcome its uncontrollable reaction process and excessive activity that leads to off-target editing. To overcome this problem, we engineered a photoswitch on G-quadruplex gRNA (GqRNA) for precisely controlled gene editing and expression by embedding dicationic azobenzene derivatives (AZD++). Our results demonstrated that rational design of the G-quadruplex onto crRNA conferred higher stability and sequence recognition specificity than unmodified single guide (sgRNA). Light-induced isomerization of AZD++ quickly transformed the on state of GqRNA, which facilitated rapid activation of ribonucleoprotein activity for genome editing of on-target sites in cells with excellent editing efficiency. In turn, AZD++-GqRNA promptly refolded to an off state to inhibit genomic cleavage, and limited the generation of off-target effects and by-products. Therefore, the proposed strategy of a photo-reversible modality presents a new opportunity for CRISPR-Cas9 modulation to improve its safety and applicability.

Two types of glycyl-tRNA synthetase (GlyRS) are known, the α2 and the α2β2 GlyRSs. Both types of synthetase employ a class II catalytic domain to aminoacylate tRNAGly. In plastids and some bacteria, the α and β subunits are fused and are designated as (αβ)2 GlyRSs. While the tRNA recognition and aminoacylation mechanisms are well understood for α2 GlyRSs, little is known about the mechanisms for α2β2/(αβ)2 GlyRSs. Here we describe structures of the (αβ)2 GlyRS from Oryza sativa chloroplast by itself and in complex with cognate tRNAGly. The set of structures reveals that the U-shaped β half of the synthetase selects the tRNA in a two-step manner. In the first step, the synthetase engages the elbow and the anticodon base C35 of the tRNA. In the second step, the tRNA has rotated ∼9° toward the catalytic centre. The synthetase probes the tRNA for the presence of anticodon base C36 and discriminator base C73. This intricate mechanism enables the tRNA to access the active site of the synthetase from a direction opposite to that of most other class II synthetases.

Magnesium, the most abundant divalent cation in cells, catalyzes RNA cleavage but also promotes RNA folding. Because folding can protect RNA from cleavage, we predicted a 'Goldilocks landscape', with local maximum in RNA lifetime at Mg2+ concentrations required for folding. Here, we use simulation and experiment to discover an innate and sophisticated mechanism of control of RNA lifetime. By simulation we characterized RNA Goldilocks landscapes and their dependence on cleavage and folding parameters. Experiments with yeast tRNAPhe and the Tetrahymena ribozyme P4-P6 domain show that structured RNAs can inhabit Goldilocks peaks. The Goldilocks peaks are tunable by differences in folded and unfolded cleavage rate constants, Mg2+ binding cooperativity, and Mg2+ affinity. Different folding and cleavage parameters produce Goldilocks landscapes with a variety of features. Goldilocks behavior allows ultrafine control of RNA chemical lifetime, whereas non-folding RNAs do not display Goldilocks peaks of protection. In sum, the effects of Mg2+ on RNA persistence are expected to be pleomorphic, both protecting and degrading RNA. In evolutionary context, Goldilocks behavior may have been a selectable trait of RNA in an early Earth environment containing Mg2+ and other metals.

Posttranscriptional regulation of the maternal nanos mRNA is essential for the development of the anterior - posterior axis of the Drosophila embryo. The nanos RNA is regulated by the protein Smaug, which binds to Smaug recognition elements (SREs) in the nanos 3'-UTR and nucleates the assembly of a larger repressor complex including the eIF4E-T paralog Cup and five additional proteins. The Smaug-dependent complex represses translation of nanos and induces its deadenylation by the CCR4-NOT deadenylase. Here we report an in vitro reconstitution of the Drosophila CCR4-NOT complex and Smaug-dependent deadenylation. We find that Smaug by itself is sufficient to cause deadenylation by the Drosophila or human CCR4-NOT complexes in an SRE-dependent manner. CCR4-NOT subunits NOT10 and NOT11 are dispensable, but the NOT module, consisting of NOT2, NOT3 and the C-terminal part of NOT1, is required. Smaug interacts with the C-terminal domain of NOT3. Both catalytic subunits of CCR4-NOT contribute to Smaug-dependent deadenylation. Whereas the CCR4-NOT complex itself acts distributively, Smaug induces a processive behavior. The cytoplasmic poly(A) binding protein (PABPC) has a minor inhibitory effect on Smaug-dependent deadenylation. Among the additional constituents of the Smaug-dependent repressor complex, Cup also facilitates CCR4-NOT-dependent deadenylation, both independently and in cooperation with Smaug.

The modified nucleosides 2'-deoxy-7-cyano- and 2'-deoxy-7-amido-7-deazaguanosine (dPreQ0 and dADG, respectively) recently discovered in DNA are the products of the bacterial queuosine tRNA modification pathway and the dpd gene cluster, the latter of which encodes proteins that comprise the elaborate Dpd restriction-modification system present in diverse bacteria. Recent genetic studies implicated the dpdA, dpdB and dpdC genes as encoding proteins necessary for DNA modification, with dpdD-dpdK contributing to the restriction phenotype. Here we report the in vitro reconstitution of the Dpd modification machinery from Salmonella enterica serovar Montevideo, the elucidation of the roles of each protein and the X-ray crystal structure of DpdA supported by small-angle X-ray scattering analysis of DpdA and DpdB, the former bound to DNA. While the homology of DpdA with the tRNA-dependent tRNA-guanine transglycosylase enzymes (TGT) in the queuosine pathway suggested a similar transglycosylase activity responsible for the exchange of a guanine base in the DNA for 7-cyano-7-deazaguanine (preQ0), we demonstrate an unexpected ATPase activity in DpdB necessary for insertion of preQ0 into DNA, and identify several catalytically essential active site residues in DpdA involved in the transglycosylation reaction. Further, we identify a modification site for DpdA activity and demonstrate that DpdC functions independently of DpdA/B in converting preQ0-modified DNA to ADG-modified DNA.

Multimodal single-cell omics technologies enable multiple molecular programs to be simultaneously profiled at a global scale in individual cells, creating opportunities to study biological systems at a resolution that was previously inaccessible. However, the analysis of multimodal single-cell omics data is challenging due to the lack of methods that can integrate across multiple data modalities generated from such technologies. Here, we present Matilda, a multi-task learning method for integrative analysis of multimodal single-cell omics data. By leveraging the interrelationship among tasks, Matilda learns to perform data simulation, dimension reduction, cell type classification, and feature selection in a single unified framework. We compare Matilda with other state-of-the-art methods on datasets generated from some of the most popular multimodal single-cell omics technologies. Our results demonstrate the utility of Matilda for addressing multiple key tasks on integrative multimodal single-cell omics data analysis. Matilda is implemented in Pytorch and is freely available from https://github.com/PYangLab/Matilda.

In Escherichia coli, the single-stranded DNA-binding protein (SSB) acts as a genome maintenance organizational hub by interacting with multiple DNA metabolism proteins. Many SSB-interacting proteins (SIPs) form complexes with SSB by docking onto its carboxy-terminal tip (SSB-Ct). An alternative interaction mode in which SIPs bind to PxxP motifs within an intrinsically-disordered linker (IDL) in SSB has been proposed for the RecG DNA helicase and other SIPs. Here, RecG binding to SSB and SSB peptides was measured in vitro and the RecG/SSB interface was identified. The results show that RecG binds directly and specifically to the SSB-Ct, and not the IDL, through an evolutionarily conserved binding site in the RecG helicase domain. Mutations that block RecG binding to SSB sensitize E. coli to DNA damaging agents and induce the SOS DNA-damage response, indicating formation of the RecG/SSB complex is important in vivo. The broader role of the SSB IDL is also investigated. E. coli ssb mutant strains encoding SSB IDL deletion variants lacking all PxxP motifs retain wildtype growth and DNA repair properties, demonstrating that the SSB PxxP motifs are not major contributors to SSB cellular functions.

The RNA exosome is an essential 3' to 5' exoribonuclease complex that mediates degradation, processing and quality control of virtually all eukaryotic RNAs. The nucleolar RNA exosome, consisting of a nine-subunit core and a distributive 3' to 5' exonuclease EXOSC10, plays a critical role in processing and degrading nucleolar RNAs, including pre-rRNA. However, how the RNA exosome is regulated in the nucleolus is poorly understood. Here, we report that the nucleolar ubiquitin-specific protease USP36 is a novel regulator of the nucleolar RNA exosome. USP36 binds to the RNA exosome through direct interaction with EXOSC10 in the nucleolus. Interestingly, USP36 does not significantly regulate the levels of EXOSC10 and other tested exosome subunits. Instead, it mediates EXOSC10 SUMOylation at lysine (K) 583. Mutating K583 impaired the binding of EXOSC10 to pre-rRNAs, and the K583R mutant failed to rescue the defects in rRNA processing and cell growth inhibition caused by knockdown of endogenous EXOSC10. Furthermore, EXOSC10 SUMOylation is markedly reduced in cells in response to perturbation of ribosomal biogenesis. Together, these results suggest that USP36 acts as a SUMO ligase to promote EXOSC10 SUMOylation critical for the RNA exosome function in ribosome biogenesis.

Combination cancer chemotherapy is one of the most useful treatment methods to achieve a synergistic effect and reduce the toxicity of dosing with a single drug. Here, we use a combination of two well-established anticancer DNA intercalators, actinomycin D (ActD) and echinomycin (Echi), to screen their binding capabilities with DNA duplexes containing different mismatches embedded within Watson-Crick base-pairs. We have found that combining ActD and Echi preferentially stabilised thymine-related T:T mismatches. The enhanced stability of the DNA duplex-drug complexes is mainly due to the cooperative binding of the two drugs to the mismatch duplex, with many stacking interactions between the two different drug molecules. Since the repair of thymine-related mismatches is less efficient in mismatch repair (MMR)-deficient cancer cells, we have also demonstrated that the combination of ActD and Echi exhibits enhanced synergistic effects against MMR-deficient HCT116 cells and synergy is maintained in a MMR-related MLH1 gene knockdown in SW620 cells. We further accessed the clinical potential of the two-drug combination approach with a xenograft mouse model of a colorectal MMR-deficient cancer, which has resulted in a significant synergistic anti-tumour effect. The current study provides a novel approach for the development of combination chemotherapy for the treatment of cancers related to DNA-mismatches.

The RNA-guided Cas9 endonuclease from Staphylococcus aureus (SauCas9) can catalyze multiple-turnover reactions whereas Cas9 from Streptococcus pyogenes (SpyCas9) is a single-turnover enzyme. Here we dissect the mechanism of multiple-turnover catalysis by SauCas9 and elucidate its molecular basis. We show that the multiple-turnover catalysis does not require more than stoichiometric RNA guides to Cas9 nuclease. Rather, the RNA-guide loaded ribonucleoprotein (RNP) is the reactive unity that is slowly released from product and recycled in the subsequent reaction. The mechanism that RNP is recycled for multiple-turnover reaction entails the unwinding of the RNA:DNA duplex in the R-loop. We argue that DNA rehybridization is required for RNP release by supplementing the energy cost in the process. Indeed, turnover is arrested when DNA rehybridization is suppressed. Further, under higher salt conditions, both SauCas9 and SpyCas9 showed increased turnover, and engineered SpyCas9 nucleases that form fewer direct or hydrogen bonding interactions with target DNA became multiple-turnover enzymes. Thus, these results indicate that for both SpyCas9 and SauCas9, turnover is determined by the energetic balance of the post-chemistry RNP-DNA interaction. Due to the conserved protein core folds, the mechanism underpinning turnover we establish here is likely operant in all Cas9 nucleases.

To perform double-stranded DNA passage, type II topoisomerases generate a covalent enzyme-cleaved DNA complex (i.e. cleavage complex). Although this complex is a requisite enzyme intermediate, it is also intrinsically dangerous to genomic stability. Consequently, cleavage complexes are the targets for several clinically relevant anticancer and antibacterial drugs. Human topoisomerase IIα and IIβ and bacterial gyrase maintain higher levels of cleavage complexes with negatively supercoiled over positively supercoiled DNA substrates. Conversely, bacterial topoisomerase IV is less able to distinguish DNA supercoil handedness. Despite the importance of supercoil geometry to the activities of type II topoisomerases, the basis for supercoil handedness recognition during DNA cleavage has not been characterized. Based on the results of benchtop and rapid-quench flow kinetics experiments, the forward rate of cleavage is the determining factor of how topoisomerase IIα/IIβ, gyrase and topoisomerase IV distinguish supercoil handedness in the absence or presence of anticancer/antibacterial drugs. In the presence of drugs, this ability can be enhanced by the formation of more stable cleavage complexes with negatively supercoiled DNA. Finally, rates of enzyme-mediated DNA ligation do not contribute to the recognition of DNA supercoil geometry during cleavage. Our results provide greater insight into how type II topoisomerases recognize their DNA substrates.

More than 170 posttranscriptional RNA modifications are so far known on both coding and noncoding RNA species. Within this group, pseudouridine (Ψ) and queuosine (Q) represent conserved RNA modifications with fundamental functional roles in regulating translation. Current detection methods of these modifications, which both are reverse transcription (RT)-silent, are mostly based on the chemical treatment of RNA prior to analysis. To overcome the drawbacks associated with indirect detection strategies, we have engineered an RT-active DNA polymerase variant called RT-KTq I614Y that produces error RT signatures specific for Ψ or Q without prior chemical treatment of the RNA samples. Combining this polymerase with next-generation sequencing techniques allows the direct identification of Ψ and Q sites of untreated RNA samples using a single enzymatic tool.

The placenta is an essential organ that supports the growth and development of the fetus during pregnancy. However, cell type-specific enhancers and transcription factors (TFs), and the mechanisms underlying the maintenance and differentiation of trophoblast stem cell (TSC) populations in the human placenta remain elusive. Here, using human TSCs as a model system, we identify 31,362 enhancers that are enriched with the motifs of previously reported TSC-pivotal TFs, including TEAD4, GATA2/3 and TFAP2C. Subsequently, we identify 580 super-enhancers (SEs) and 549 SE-associated genes. These genes are robustly expressed in the human placenta and include numerous TFs, implying that SE-associated TFs (SE-TFs) may play crucial roles in placental development. Additionally, we identify the global binding sites of five TSC-pivotal SE-TFs (FOS, GATA2, MAFK, TEAD4 and TFAP2C), revealing that they preferentially co-occupy enhancers, regulate each other and form a trophoblast-active gene regulatory network. Loss-of-function studies unveil that the five TFs promote self-renewal of TSCs by activating proliferation-associated genes while repressing developmental genes. We further reveal that the five TFs exert conserved and unique functions on placental development between humans and mice. Our study provides important insights into the roles of human TSC-pivotal TFs in regulating placenta-specific gene expression programs.

High fidelity tRNA aminoacylation by aminoacyl-tRNA synthetases is essential for cell viability. ProXp-ala is a trans-editing protein that is present in all three domains of life and is responsible for hydrolyzing mischarged Ala-tRNAPro and preventing mistranslation of proline codons. Previous studies have shown that, like bacterial prolyl-tRNA synthetase, Caulobacter crescentus ProXp-ala recognizes the unique C1:G72 terminal base pair of the tRNAPro acceptor stem, helping to ensure deacylation of Ala-tRNAPro but not Ala-tRNAAla. The structural basis for C1:G72 recognition by ProXp-ala is still unknown and was investigated here. NMR spectroscopy, binding, and activity assays revealed two conserved residues, K50 and R80, that likely interact with the first base pair, stabilizing the initial protein-RNA encounter complex. Modeling studies are consistent with direct interaction between R80 and the major groove of G72. A third key contact between A76 of tRNAPro and K45 of ProXp-ala was essential for binding and accommodating the CCA-3' end in the active site. We also demonstrated the essential role that the 2'OH of A76 plays in catalysis. Eukaryotic ProXp-ala proteins recognize the same acceptor stem positions as their bacterial counterparts, albeit with different nucleotide base identities. ProXp-ala is encoded in some human pathogens; thus, these results have the potential to inform new antibiotic drug design.

Thioflavin T (ThT) is a classical fluorescent dye gaining prominence in current research regarding nucleic acid conformations (NACs). However, most NACs with the ability to excite ThT fluorescent are unique or form in demanding conditions, limiting the extensiveness and depth of ThT application in sensing and imaging. Therefore, this study proposed CGG-AAA mismatched cavity hairpin ThT-light nucleic acid switches (CHTLNAS) with excellent fluorescence excitation over 500-fold higher than spontaneous, 17∼20-fold higher than ssDNA and 2.5∼5-fold higher than complementary duplex. Based on the excellent fluorescence excitation, convenient conformation formation, good sequence programmability, and flexible allosteric ability (known as the Worm-crack pod mechanism mediated by the target), it achieved the label- and enzyme-free detection of tetracycline (TET) and berberine (BB) at the pM level within 10 min. Moreover, it was found enable to realize the sensitive tracking of intracellular carriers at the nM level of ThT entry concentration, and prolongated its cell nuclear-entry time of ThT over 8 h, overcoming the non-specific high background signal interference of ThT in the nuclear region, and expanding the diversified application of ThT in cell biology research. Therefore, CHTLNAS is a more universal, practical tool than G-quadruplex or other kinds of NACs for ThT development and utilization in sensing and imaging platforms.

Meiotic recombinases RAD51 and DMC1 mediate strand exchange in the repair of DNA double-strand breaks (DSBs) by homologous recombination. This is a landmark event of meiosis that ensures genetic diversity in sexually reproducing organisms. However, the regulatory mechanism of DMC1/RAD51-ssDNA nucleoprotein filaments during homologous recombination in mammals has remained largely elusive. Here, we show that SPIDR (scaffold protein involved in DNA repair) regulates the assembly or stability of RAD51/DMC1 on ssDNA. Knockout of Spidr in male mice causes complete meiotic arrest, accompanied by defects in synapsis and crossover formation, which leads to male infertility. In females, loss of Spidr leads to subfertility; some Spidr-/- oocytes are able to complete meiosis. Notably, fertility is rescued partially by ablation of the DNA damage checkpoint kinase CHK2 in Spidr-/- females but not in males. Thus, our study identifies SPIDR as an essential meiotic recombination factor in homologous recombination in mammals.

The histone acetyltransferase (HAT) subunit of coactivator complex SAGA, Gcn5, stimulates eviction of promoter nucleosomes at certain highly expressed yeast genes, including those activated by transcription factor Gcn4 in amino acid-deprived cells; however, the importance of other HAT complexes in this process was poorly understood. Analyzing mutations that disrupt the integrity or activity of HAT complexes NuA4 or NuA3, or HAT Rtt109, revealed that only NuA4 acts on par with Gcn5, and functions additively, in evicting and repositioning promoter nucleosomes and stimulating transcription of starvation-induced genes. NuA4 is generally more important than Gcn5, however, in promoter nucleosome eviction, TBP recruitment, and transcription at most other genes expressed constitutively. NuA4 also predominates over Gcn5 in stimulating TBP recruitment and transcription of genes categorized as principally dependent on the cofactor TFIID versus SAGA, except for the most highly expressed subset including ribosomal protein genes, where Gcn5 contributes strongly to PIC assembly and transcription. Both SAGA and NuA4 are recruited to promoter regions of starvation-induced genes in a manner that might be feedback controlled by their HAT activities. Our findings reveal an intricate interplay between these two HATs in nucleosome eviction, PIC assembly, and transcription that differs between the starvation-induced and basal transcriptomes.

Bacteriophage exclusion ('BREX') systems are multi-protein complexes encoded by a variety of bacteria and archaea that restrict phage by an unknown mechanism. One BREX factor, termed BrxL, has been noted to display sequence similarity to various AAA+ protein factors including Lon protease. In this study we describe multiple CryoEM structures of BrxL that demonstrate it to be a chambered, ATP-dependent DNA binding protein. The largest BrxL assemblage corresponds to a dimer of heptamers in the absence of bound DNA, versus a dimer of hexamers when DNA is bound in its central pore. The protein displays DNA-dependent ATPase activity, and ATP binding promotes assembly of the complex on DNA. Point mutations within several regions of the protein-DNA complex alter one or more in vitro behaviors and activities, including ATPase activity and ATP-dependent association with DNA. However, only the disruption of the ATPase active site fully eliminates phage restriction, indicating that other mutations can still complement BrxL function within the context of an otherwise intact BREX system. BrxL displays significant structural homology to MCM subunits (the replicative helicase in archaea and eukaryotes), implying that it and other BREX factors may collaborate to disrupt initiation of phage DNA replication.

A haploid of the fission yeast Schizosaccharomyces pombe expresses either the P or M mating-type, determined by the active, euchromatic, mat1 cassette. Mating-type is switched by Rad51-driven gene conversion of mat1 using a heterochromatic donor cassette, mat2-P or mat3-M. The Swi2-Swi5 complex, a mating-type switching factor, is central to this process by designating a preferred donor in a cell-type-specific manner. Swi2-Swi5 selectively enables one of two cis-acting recombination enhancers, SRE2 adjacent to mat2-P or SRE3 adjacent to mat3-M. Here, we identified two functionally important motifs in Swi2, a Swi6 (HP1 homolog)-binding site and two DNA-binding AT-hooks. Genetic analysis demonstrated that the AT-hooks were required for Swi2 localization at SRE3 to select the mat3-M donor in P cells, while the Swi6-binding site was required for Swi2 localization at SRE2 to select mat2-P in M cells. In addition, the Swi2-Swi5 complex promoted Rad51-driven strand exchange in vitro. Taken together, our results show how the Swi2-Swi5 complex would localize to recombination enhancers through a cell-type specific binding mechanism and stimulate Rad51-driven gene conversion at the localization site.

During every cell cycle, both the genome and the associated chromatin must be accurately replicated. Chromatin Assembly Factor-1 (CAF-1) is a key regulator of chromatin replication, but how CAF-1 functions in relation to the DNA replication machinery is unknown. Here, we reveal that this crosstalk differs between the leading and lagging strand at replication forks. Using biochemical reconstitutions, we show that DNA and histones promote CAF-1 recruitment to its binding partner PCNA and reveal that two CAF-1 complexes are required for efficient nucleosome assembly under these conditions. Remarkably, in the context of the replisome, CAF-1 competes with the leading strand DNA polymerase epsilon (Polϵ) for PCNA binding. However, CAF-1 does not affect the activity of the lagging strand DNA polymerase Delta (Polδ). Yet, in cells, CAF-1 deposits newly synthesized histones equally on both daughter strands. Thus, on the leading strand, chromatin assembly by CAF-1 cannot occur simultaneously to DNA synthesis, while on the lagging strand these processes may be coupled. We propose that these differences may facilitate distinct parental histone recycling mechanisms and accommodate the inherent asymmetry of DNA replication.

In situ capturing technologies add tissue context to gene expression data, with the potential of providing a greater understanding of complex biological systems. However, splicing variants and full-length sequence heterogeneity cannot be characterized at spatial resolution with current transcriptome profiling methods. To that end, we introduce spatial isoform transcriptomics (SiT), an explorative method for characterizing spatial isoform variation and sequence heterogeneity using long-read sequencing. We show in mouse brain how SiT can be used to profile isoform expression and sequence heterogeneity in different areas of the tissue. SiT reveals regional isoform switching of Plp1 gene between different layers of the olfactory bulb, and the use of external single-cell data allows the nomination of cell types expressing each isoform. Furthermore, SiT identifies differential isoform usage for several major genes implicated in brain function (Snap25, Bin1, Gnas) that are independently validated by in situ sequencing. SiT also provides for the first time an in-depth A-to-I RNA editing map of the adult mouse brain. Data exploration can be performed through an online resource (https://www.isomics.eu), where isoform expression and RNA editing can be visualized in a spatial context.

Biobank projects are generating genomic data for many thousands of individuals. Computational methods are needed to handle these massive data sets, including genetic ancestry (GA) inference tools. Current methods for GA inference do not scale to biobank-size genomic datasets. We present Rye-a new algorithm for GA inference at biobank scale. We compared the accuracy and runtime performance of Rye to the widely used RFMix, ADMIXTURE and iAdmix programs and applied it to a dataset of 488221 genome-wide variant samples from the UK Biobank. Rye infers GA based on principal component analysis of genomic variant samples from ancestral reference populations and query individuals. The algorithm's accuracy is powered by Metropolis-Hastings optimization and its speed is provided by non-negative least squares regression. Rye produces highly accurate GA estimates for three-way admixed populations-African, European and Native American-compared to RFMix and ADMIXTURE (${R}^2 = \ 0.998 - 1.00$), and shows 50× runtime improvement compared to ADMIXTURE on the UK Biobank dataset. Rye analysis of UK Biobank samples demonstrates how it can be used to infer GA at both continental and subcontinental levels. We discuss user consideration and options for the use of Rye; the program and its documentation are distributed on the GitHub repository: https://github.com/healthdisparities/rye.

DNA inverted repeats (IRs) are widespread across many eukaryotic genomes. Their ability to form stable hairpin/cruciform secondary structures is causative in triggering chromosome instability leading to several human diseases. Distance and sequence divergence between IRs are inversely correlated with their ability to induce gross chromosomal rearrangements (GCRs) because of a lesser probability of secondary structure formation and chromosomal breakage. In this study, we demonstrate that structural parameters that normally constrain the instability of IRs are overcome when the repeats interact in single-stranded DNA (ssDNA). We established a system in budding yeast whereby >73 kb of ssDNA can be formed in cdc13-707fs mutants. We found that in ssDNA, 12 bp or 30 kb spaced Alu-IRs show similarly high levels of GCRs, while heterology only beyond 25% suppresses IR-induced instability. Mechanistically, rearrangements arise after cis-interaction of IRs leading to a DNA fold-back and the formation of a dicentric chromosome, which requires Rad52/Rad59 for IR annealing as well as Rad1-Rad10, Slx4, Msh2/Msh3 and Saw1 proteins for nonhomologous tail removal. Importantly, using structural characteristics rendering IRs permissive to DNA fold-back in yeast, we found that ssDNA regions mapped in cancer genomes contain a substantial number of potentially interacting and unstable IRs.

MINORg is an offline gRNA design tool that generates the smallest possible combination of gRNA capable of covering all desired targets in multiple non-reference genomes. As interest in pangenomic research grows, so does the workload required for large screens in multiple individuals. MINORg aims to lessen this workload by capitalising on sequence homology to favour multi-target gRNA while simultaneously screening multiple genetic backgrounds in order to generate reusable gRNA panels. We demonstrated the practical application of MINORg by knocking out 11 homologous genes tandemly arrayed in a multi-gene cluster in two Arabidopsis thaliana lineages using three gRNA output by MINORg. We also described a new PCR-free modular cloning system for multiplexing gRNA, and used it to knockout three tandemly arrayed genes in another multi-gene cluster with gRNA designed by MINORg. Source code is freely available at https://github.com/rlrq/MINORg.

16S rRNA gene sequence clustering is an important tool in characterizing the diversity of microbial communities. As 16S rRNA gene data sets are growing in size, existing sequence clustering algorithms increasingly become an analytical bottleneck. Part of this bottleneck is due to the substantial computational cost expended on small clusters and singleton sequences. We propose an iterative sampling-based 16S rRNA gene sequence clustering approach that targets the largest clusters in the data set, allowing users to stop the clustering process when sufficient clusters are available for the specific analysis being targeted. We describe a probabilistic analysis of the iterative clustering process that supports the intuition that the clustering process identifies the larger clusters in the data set first. Using real data sets of 16S rRNA gene sequences, we show that the iterative algorithm, coupled with an adaptive sampling process and a mode-shifting strategy for identifying cluster representatives, substantially speeds up the clustering process while being effective at capturing the large clusters in the data set. The experiments also show that SCRAPT (Sample, Cluster, Recruit, AdaPt and iTerate) is able to produce operational taxonomic units that are less fragmented than popular tools: UCLUST, CD-HIT and DNACLUST. The algorithm is implemented in the open-source package SCRAPT. The source code used to generate the results presented in this paper is available at https://github.com/hsmurali/SCRAPT.

PBRM1 is a subunit of the PBAF chromatin remodeling complex, which is mutated in 40-50% of clear cell renal cell carcinoma patients. It is thought to largely function as a chromatin binding subunit of the PBAF complex, but the molecular mechanism underlying this activity is not fully known. PBRM1 contains six tandem bromodomains which are known to cooperate in binding of nucleosomes acetylated at histone H3 lysine 14 (H3K14ac). Here, we demonstrate that the second and fourth bromodomains from PBRM1 also bind nucleic acids, selectively associating with double stranded RNA elements. Disruption of the RNA binding pocket is found to compromise PBRM1 chromatin binding and inhibit PBRM1-mediated cellular growth effects.

The study of cellular networks mediated by ligand-receptor interactions has attracted much attention recently owing to single-cell omics. However, rich collections of bulk data accompanied with clinical information exists and continue to be generated with no equivalent in single-cell so far. In parallel, spatial transcriptomic (ST) analyses represent a revolutionary tool in biology. A large number of ST projects rely on multicellular resolution, for instance the Visium™ platform, where several cells are analyzed at each location, thus producing localized bulk data. Here, we describe BulkSignalR, a R package to infer ligand-receptor networks from bulk data. BulkSignalR integrates ligand-receptor interactions with downstream pathways to estimate statistical significance. A range of visualization methods complement the statistics, including functions dedicated to spatial data. We demonstrate BulkSignalR relevance using different datasets, including new Visium liver metastasis ST data, with experimental validation of protein colocalization. A comparison with other ST packages shows the significantly higher quality of BulkSignalR inferences. BulkSignalR can be applied to any species thanks to its built-in generic ortholog mapping functionality.

CRISPR-associated transposases (CASTs) direct DNA integration downstream of target sites using the RNA-guided DNA binding activity of nuclease-deficient CRISPR-Cas systems. Transposition relies on several key protein-protein and protein-DNA interactions, but little is known about the explicit sequence requirements governing efficient transposon DNA integration activity. Here, we exploit pooled library screening and high-throughput sequencing to reveal novel sequence determinants during transposition by the Type I-F Vibrio cholerae CAST system (VchCAST). On the donor DNA, large transposon end libraries revealed binding site nucleotide preferences for the TnsB transposase, as well as an additional conserved region that encoded a consensus binding site for integration host factor (IHF). Remarkably, we found that VchCAST requires IHF for efficient transposition, thus revealing a novel cellular factor involved in CRISPR-associated transpososome assembly. On the target DNA, we uncovered preferred sequence motifs at the integration site that explained previously observed heterogeneity with single-base pair resolution. Finally, we exploited our library data to design modified transposon variants that enable in-frame protein tagging. Collectively, our results provide new clues about the assembly and architecture of the paired-end complex formed between TnsB and the transposon DNA, and inform the design of custom payload sequences for genome engineering applications with CAST systems.

Adenosine deaminase acting on RNA ADAR1 promotes A-to-I conversion in double-stranded and structured RNAs. ADAR1 has two isoforms transcribed from different promoters: cytoplasmic ADAR1p150 is interferon-inducible while ADAR1p110 is constitutively expressed and primarily localized in the nucleus. Mutations in ADAR1 cause Aicardi - Goutières syndrome (AGS), a severe autoinflammatory disease associated with aberrant IFN production. In mice, deletion of ADAR1 or the p150 isoform leads to embryonic lethality driven by overexpression of interferon-stimulated genes. This phenotype is rescued by deletion of the cytoplasmic dsRNA-sensor MDA5 indicating that the p150 isoform is indispensable and cannot be rescued by ADAR1p110. Nevertheless, editing sites uniquely targeted by ADAR1p150 remain elusive. Here, by transfection of ADAR1 isoforms into ADAR-less mouse cells we detect isoform-specific editing patterns. Using mutated ADAR variants, we test how intracellular localization and the presence of a Z-DNA binding domain-α affect editing preferences. These data show that ZBDα only minimally contributes to p150 editing-specificity while isoform-specific editing is primarily directed by the intracellular localization of ADAR1 isoforms. Our study is complemented by RIP-seq on human cells ectopically expressing tagged-ADAR1 isoforms. Both datasets reveal enrichment of intronic editing and binding by ADAR1p110 while ADAR1p150 preferentially binds and edits 3'UTRs.

Cadmium (Cd) is one of the most toxic heavy metals. Exposure to Cd can impair the functions of the kidney, respiratory system, reproductive system and skeletal system. Cd2+-binding aptamers have been extensively utilized in the development of Cd2+-detecting devices; however, the underlying mechanisms remain elusive. This study reports four Cd2+-bound DNA aptamer structures, representing the only Cd2+-specific aptamer structures available to date. In all the structures, the Cd2+-binding loop (CBL-loop) adopts a compact, double-twisted conformation and the Cd2+ ion is mainly coordinated with the G9, C12 and G16 nucleotides. Moreover, T11 and A15 within the CBL-loop form one regular Watson-Crick pair and stabilize the conformation of G9. The conformation of G16 is stabilized by the G8-C18 pair of the stem. By folding and/or stabilizing the CBL-loop, the other four nucleotides of the CBL-loop also play important roles in Cd2+ binding. Similarly to the native sequence, crystal structures, circular dichroism spectrum and isothermal titration calorimetry analysis confirm that several variants of the aptamer can recognize Cd2+. This study not only reveals the underlying basis for the binding of Cd2+ ions with the aptamer, but also extends the sequence for the construction of novel metal-DNA complex.

The human gut microbiome has been linked to health and disease. Investigation of the human microbiome has largely employed 16S amplicon sequencing, with limited ability to distinguish microbes at the species level. Herein, we describe the development of Reference-based Exact Mapping (RExMap) of microbial amplicon variants that enables mapping of microbial species from standard 16S sequencing data. RExMap analysis of 16S data captures ∼75% of microbial species identified by whole-genome shotgun sequencing, despite hundreds-fold less sequencing depth. RExMap re-analysis of existing 16S data from 29,349 individuals across 16 regions from around the world reveals a detailed landscape of gut microbial species across populations and geography. Moreover, RExMap identifies a core set of fifteen gut microbes shared by humans. Core microbes are established soon after birth and closely associate with BMI across multiple independent studies. RExMap and the human microbiome dataset are presented as resources with which to explore the role of the human microbiome.

The functions of non-coding RNAs usually depend on their 3D structures. Therefore, comparing RNA 3D structures is critical in analyzing their functions. We noticed an interesting phenomenon that two non-coding RNAs may share similar substructures when rotating their sequence order. To the best of our knowledge, no existing RNA 3D structural alignment tools can detect this type of matching. In this article, we defined the RNA 3D structure circular matching problem and developed a software tool named CircularSTAR3D to solve this problem. CircularSTAR3D first uses the conserved stacks (consecutive base pairs with similar 3D structures) in the input RNAs to identify the circular matched internal loops and multiloops. Then it performs a local extension iteratively to obtain the whole circular matched substructures. The computational experiments conducted on a non-redundant RNA structure dataset show that circular matching is ubiquitous. Furthermore, we demonstrated the utility of CircularSTAR3D by detecting the conserved substructures missed by regular alignment tools, including structural motifs and conserved structures between riboswitches and ribozymes from different classes. We anticipate CircularSTAR3D to be a valuable supplement to the existing RNA 3D structural analysis techniques.

In response to different stimuli many transcription factors (TFs) display different activation dynamics that trigger the expression of specific sets of target genes, suggesting that promoters have a way to decode dynamics. Here, we use optogenetics to directly manipulate the nuclear localization of a synthetic TF in mammalian cells without affecting other processes. We generate pulsatile or sustained TF dynamics and employ live cell microscopy and mathematical modelling to analyse the behaviour of a library of reporter constructs. We find decoding of TF dynamics occurs only when the coupling between TF binding and transcription pre-initiation complex formation is inefficient and that the ability of a promoter to decode TF dynamics gets amplified by inefficient translation initiation. Using the knowledge acquired, we build a synthetic circuit that allows obtaining two gene expression programs depending solely on TF dynamics. Finally, we show that some of the promoter features identified in our study can be used to distinguish natural promoters that have previously been experimentally characterized as responsive to either sustained or pulsatile p53 and NF-κB signals. These results help elucidate how gene expression is regulated in mammalian cells and open up the possibility to build complex synthetic circuits steered by TF dynamics.

Precise control of gene expression is essential for flux redistribution in metabolic pathways. Although the CRISPR interference (CRISPRi) system can effectively repress gene expression at the transcriptional level, it has still been difficult to precisely control the level without loss of specificity or an increase in cell toxicity. In this study, we developed a tunable CRISPRi system that performs transcriptional regulation at various levels. We constructed a single-guide RNA (sgRNA) library targeting repeat, tetraloop, and anti-repeat regions to modulate the binding affinity against dCas9. Each screened sgRNA could regulate the gene expression at a certain level between fully-repressing and non-repressing states (>45-fold). These sgRNAs also enabled modular regulation with various target DNA sequences. We applied this system to redistribute the metabolic flux to produce violacein derivatives in a predictable ratio and optimize lycopene production. This system would help accelerate the flux optimization processes in metabolic engineering and synthetic biology.

The long non-coding telomeric RNA transcript TERRA, in the form of an RNA-DNA duplex, regulates telomere recombination. In a screen for nucleases that affects telomere recombination, mutations in DNA2, EXO1, MRE11 and SAE2 cause severe delay in type II survivor formation, indicating that type II telomere recombination is mediated through a mechanism similar to repairing double-strand breaks. On the other hand, mutation in RAD27 results in early formation of type II recombination, suggesting that RAD27 acts as a negative regulator in telomere recombination. RAD27 encodes a flap endonuclease that plays a role in DNA metabolism, including replication, repair and recombination. We demonstrate that Rad27 suppresses the accumulation of the TERRA-associated R-loop and selectively cleaves TERRA of R-loop and double-flapped structures in vitro. Moreover, we show that Rad27 negatively regulates single-stranded C-rich telomeric DNA circles (C-circles) in telomerase-deficient cells, revealing a close correlation between R-loop and C-circles during telomere recombination. These results demonstrate that Rad27 participates in telomere recombination by cleaving TERRA in the context of an R-loop or flapped RNA-DNA duplex, providing mechanistic insight into how Rad27 maintains chromosome stability by restricting the accumulation of the R-loop structure within the genome.

A methyltransferase ribozyme (MTR1) was selected in vitro to catalyze alkyl transfer from exogenous O6-methylguanine (O6mG) to a target adenine N1, and recently, high-resolution crystal structures have become available. We use a combination of classical molecular dynamics, ab initio quantum mechanical/molecular mechanical (QM/MM) and alchemical free energy (AFE) simulations to elucidate the atomic-level solution mechanism of MTR1. Simulations identify an active reactant state involving protonation of C10 that hydrogen bonds with O6mG:N1. The deduced mechanism involves a stepwise mechanism with two transition states corresponding to proton transfer from C10:N3 to O6mG:N1 and rate-controlling methyl transfer (19.4  kcal·mol-1 barrier). AFE simulations predict the pKa for C10 to be 6.3, close to the experimental apparent pKa of 6.2, further implicating it as a critical general acid. The intrinsic rate derived from QM/MM simulations, together with pKa calculations, enables us to predict an activity-pH profile that agrees well with experiment. The insights gained provide further support for a putative RNA world and establish new design principles for RNA-based biochemical tools.

DNA sequence composition determines the topology and stability of G-quadruplexes (G4s). Bulged G-quadruplex structures (G4-Bs) are a subset of G4s characterized by 3D conformations with bulges. Current search algorithms fail to capture stable G4-B, making their genome-wide study infeasible. Here, we introduced a large family of computationally defined and experimentally verified potential G4-B forming sequences (pG4-BS). We found 478 263 pG4-BS regions that do not overlap 'canonical' G4-forming sequences in the human genome and are preferentially localized in transcription regulatory regions including R-loops and open chromatin. Over 90% of protein-coding genes contain pG4-BS in their promoter or gene body. We observed generally higher pG4-BS content in R-loops and their flanks, longer genes that are associated with brain tissue, immune and developmental processes. Also, the presence of pG4-BS on both template and non-template strands in promoters is associated with oncogenesis, cardiovascular disease and stemness. Our G4-BS models predicted G4-forming ability in vitro with 91.5% accuracy. Analysis of G4-seq and CUT&Tag data strongly supports the existence of G4-BS conformations genome-wide. We reconstructed a novel G4-B 3D structure located in the E2F8 promoter. This study defines a large family of G4-like sequences, offering new insights into the essential biological functions and potential future therapeutic uses of G4-B.

The long interspersed element 1 (LINE-1 or L1) integration is affected by many cellular factors through various mechanisms. Some of these factors are required for L1 amplification, while others either suppress or enhance specific steps during L1 propagation. Previously, TRIM28 has been identified to suppress transposable elements, including L1 expression via its canonical role in chromatin remodeling. Here, we report that TRIM28 through its B box domain increases L1 retrotransposition and facilitates shorter cDNA and L1 insert generation in cultured cells. Consistent with the latter, we observe that tumor specific L1 inserts are shorter in endometrial, ovarian, and prostate tumors with higher TRIM28 mRNA expression than in those with lower TRIM28 expression. We determine that three amino acids in the B box domain that are involved in TRIM28 multimerization are critical for its effect on both L1 retrotransposition and cDNA synthesis. We provide evidence that B boxes from the other two members in the Class VI TRIM proteins, TRIM24 and TRIM33, also increase L1 retrotransposition. Our findings could lead to a better understanding of the host/L1 evolutionary arms race in the germline and their interplay during tumorigenesis.

Precise genome editing requires the resolution of nuclease-induced DNA double strand breaks (DSBs) via the homology-directed repair (HDR) pathway. In mammals, this is typically outcompeted by non-homologous end-joining (NHEJ) that can generate potentially genotoxic insertion/deletion mutations at DSB sites. Because of higher efficacy, clinical genome editing has been restricted to imperfect but efficient NHEJ-based approaches. Hence, strategies that promote DSB resolution via HDR are essential to facilitate clinical transition of HDR-based editing strategies and increase safety. Here we describe a novel platform that consists of a Cas9 fused to DNA repair factors to synergistically inhibit NHEJ and favor HDR for precise repairing of Cas-induced DSBs. Compared to canonical CRISPR/Cas9, the increase in error-free editing ranges from 1.5-fold to 7-fold in multiple cell lines and in primary human cells. This novel CRISPR/Cas9 platform accepts clinically relevant repair templates, such as oligodeoxynucleotides (ODNs) and adeno-associated virus (AAV)-based vectors, and has a lower propensity to induce chromosomal translocations as compared to benchmark CRISPR/Cas9. The observed reduced mutational burden, resulting from diminished indel formation at on- and off-target sites, provides a remarkable gain in safety and advocates this novel CRISPR system as an attractive tool for therapeutic applications depending on precision genome editing.

Interferon regulatory factor 1 (IRF1) is a critical component of cell-intrinsic innate immunity that regulates both constitutive and induced antiviral defenses. Due to its short half-life, IRF1 function is generally considered to be regulated by its synthesis. However, how IRF1 activity is controlled post-translationally has remained poorly characterized. Here, we employed a proteomics approach to identify proteins interacting with IRF1, and found that CSNK2B, a regulatory subunit of casein kinase 2, interacts directly with IRF1 and constitutively modulates its transcriptional activity. Genome-wide CUT&RUN analysis of IRF1 binding loci revealed that CSNK2B acts generally to enhance the binding of IRF1 to chromatin, thereby enhancing transcription of key antiviral genes, such as PLAAT4 (also known as RARRES3/RIG1/TIG3). On the other hand, depleting CSNK2B triggered abnormal accumulation of IRF1 at AFAP1 loci, thereby down-regulating transcription of AFAP1, revealing contrary effects of CSNK2B on IRF1 binding at different loci. AFAP1 encodes an actin crosslinking factor that mediates Src activation. Importantly, CSNK2B was also found to mediate phosphorylation-dependent activation of AFAP1-Src signaling and exert suppressive effects against flaviviruses, including dengue virus. These findings reveal a previously unappreciated mode of IRF1 regulation and identify important effector genes mediating multiple cellular functions governed by CSNK2B and IRF1.

Tryptophanyl-tRNA synthetase (TrpRS) links tryptophan to tRNATrp, thereby playing an indispensable role in protein translation. Unlike most class I aminoacyl-tRNA synthetases (AARSs), TrpRS functions as a homodimer. Herein, we captured an 'open-closed' asymmetric structure of Escherichia coli TrpRS (EcTrpRS) with one active site occupied by a copurified intermediate product and the other remaining empty, providing structural evidence for the long-discussed half-of-the-sites reactivity of bacterial TrpRS. In contrast to its human counterpart, bacterial TrpRS may rely on this asymmetric conformation to functionally bind with substrate tRNA. As this asymmetric conformation is probably a dominant form of TrpRS purified from bacterial cells, we performed fragment screening against asymmetric EcTrpRS to support antibacterial discovery. Nineteen fragment hits were identified, and 8 of them were successfully cocrystallized with EcTrpRS. While a fragment named niraparib bound to the L-Trp binding site of the 'open' subunit, the other 7 fragments all bound to an unprecedented pocket at the interface between two TrpRS subunits. Binding of these fragments relies on residues specific to bacterial TrpRS, avoiding undesired interactions with human TrpRS. These findings improve our understanding of the catalytic mechanism of this important enzyme and will also facilitate the discovery of bacterial TrpRS inhibitors with therapeutic potential.

Mobile genetic elements (MGEs) are key promoters of microbial evolution. These elements can be located extrachromosomally or integrated into the chromosome. Well-known examples of chromosomally integrated MGEs (ciMGEs) are integrative and conjugative/mobilizable elements (ICEs and IMEs), and most studies to date have focused on the biological mechanisms that shape their lifestyle. It is crucial to profile the diversity and understand their distribution across the microbial community, as the number of genome sequences increases exponentially. Herein, I scanned a collection of >20 000 bacterial and archaeal non-redundant genomes and found over 13 000 ciMGEs across multiple phyla, representing a massive increase in the number of ciMGEs available in public databases (<1000). Although ICEs are the most important ciMGEs for the accretion of defense systems, virulence, and antimicrobial resistance (AMR) genes, IMEs outnumbered ICEs. Moreover, defense systems, AMR, and virulence genes were negatively correlated in both ICEs and IMEs. Multiple ciMGEs form heterogeneous communities and challenge inter-phylum barriers. Finally, I observed that the functional landscape of ICEs was populated by uncharacterized proteins. Altogether, this study provides a comprehensive catalog of nucleotide sequences and associated metadata for ciMGEs from 34 phyla across the bacterial and archaeal domains.

Herein, we report the systematic investigation of stereopure phosphorothioate (PS) and phosphoryl guanidine (PN) linkages on siRNA-mediated silencing. The incorporation of appropriately positioned and configured stereopure PS and PN linkages to N-acetylgalactosamine (GalNAc)-conjugated siRNAs based on multiple targets (Ttr and HSD17B13) increased potency and durability of mRNA silencing in mouse hepatocytes in vivo compared with reference molecules based on clinically proven formats. The observation that the same modification pattern had beneficial effects on unrelated transcripts suggests that it may be generalizable. The effect of stereopure PN modification on silencing is modulated by 2'-ribose modifications in the vicinity, particularly on the nucleoside 3' to the linkage. These benefits corresponded with both an increase in thermal instability at the 5'-end of the antisense strand and improved Argonaute 2 (Ago2) loading. Application of one of our most effective designs to generate a GalNAc-siRNA targeting human HSD17B13 led to ∼80% silencing that persisted for at least 14 weeks after administration of a single 3 mg/kg subcutaneous dose in transgenic mice. The judicious use of stereopure PN linkages improved the silencing profile of GalNAc-siRNAs without disrupting endogenous RNA interference pathways and without elevating serum biomarkers for liver dysfunction, suggesting they may be suitable for therapeutic application.

Type IIS restriction endonucleases contain separate DNA recognition and catalytic domains and cleave their substrates at well-defined distances outside their target sequences. They are employed in biotechnology for a variety of purposes, including the creation of gene-targeting zinc finger and TAL effector nucleases and DNA synthesis applications such as Golden Gate assembly. The most thoroughly studied Type IIS enzyme, FokI, has been shown to require multimerization and engagement with multiple DNA targets for optimal cleavage activity; however, details of how it or similar enzymes forms a DNA-bound reaction complex have not been described at atomic resolution. Here we describe biochemical analyses of DNA cleavage by the Type IIS PaqCI restriction endonuclease and a series of molecular structures in the presence and absence of multiple bound DNA targets. The enzyme displays a similar tetrameric organization of target recognition domains in the absence or presence of bound substrate, with a significant repositioning of endonuclease domains in a trapped DNA-bound complex that is poised to deliver the first of a series of double-strand breaks. PaqCI and FokI share similar structural mechanisms of DNA cleavage, but considerable differences in their domain organization and quaternary architecture, facilitating comparisons between distinct Type IIS enzymes.

RNA-binding proteins (RPBs) are deeply involved in fundamental cellular processes in bacteria and are vital for their survival. Despite this, few studies have so far been dedicated to direct and global identification of bacterial RBPs. We have adapted the RNA interactome capture (RIC) technique, originally developed for eukaryotic systems, to globally identify RBPs in bacteria. RIC takes advantage of the base pairing potential of poly(A) tails to pull-down RNA-protein complexes. Overexpressing poly(A) polymerase I in Escherichia coli drastically increased transcriptome-wide RNA polyadenylation, enabling pull-down of crosslinked RNA-protein complexes using immobilized oligo(dT) as bait. With this approach, we identified 169 putative RBPs, roughly half of which are already annotated as RNA-binding. We experimentally verified the RNA-binding ability of a number of uncharacterized RBPs, including YhgF, which is exceptionally well conserved not only in bacteria, but also in archaea and eukaryotes. We identified YhgF RNA targets in vivo using CLIP-seq, verified specific binding in vitro, and reveal a putative role for YhgF in regulation of gene expression. Our findings present a simple and robust strategy for RBP identification in bacteria, provide a resource of new bacterial RBPs, and lay the foundation for further studies of the highly conserved RBP YhgF.

The importance of non-canonical DNA structures such as G-quadruplexes (G4) and intercalating-motifs (iMs) in the fine regulation of a variety of cellular processes has been recently demonstrated. As the crucial roles of these structures are being unravelled, it is becoming more and more important to develop tools that allow targeting these structures with the highest possible specificity. While targeting methodologies have been reported for G4s, this is not the case for iMs, as evidenced by the limited number of specific ligands able to bind the latter and the total absence of selective alkylating agents for their covalent targeting. Furthermore, strategies for the sequence-specific covalent targeting of G4s and iMs have not been reported thus far. Herein, we describe a simple methodology to achieve sequence-specific covalent targeting of G4 and iM DNA structures based on the combination of (i) a peptide nucleic acid (PNA) recognizing a specific sequence of interest, (ii) a pro-reactive moiety enabling a controlled alkylation reaction, and (iii) a G4 or iM ligand orienting the alkylating warhead to the reactive residues. This multi-component system allows for the targeting of specific G4 or iM sequences of interest in the presence of competing DNA sequences and under biologically relevant conditions.

Crossing over between homologs is critical for the stable segregation of chromosomes during the first meiotic division. Saccharomyces cerevisiae Mer3 (HFM1 in mammals) is a SF2 helicase and member of the ZMM group of proteins, that facilitates the formation of the majority of crossovers during meiosis. Here, we describe the structural organisation of Mer3 and using AlphaFold modelling and XL-MS we further characterise the previously described interaction with Mlh1-Mlh2. We find that Mer3 also forms a previously undescribed complex with the recombination regulating factors Top3 and Rmi1 and that this interaction is competitive with Sgs1BLM helicase. Using in vitro reconstituted D-loop assays we show that Mer3 inhibits the anti-recombination activity of Sgs1 helicase, but only in the presence of Dmc1. Thus we provide a mechanism whereby Mer3 interacts with a network of proteins to protect Dmc1 derived D-loops from dissolution.

BRCA1 mutations are associated with increased breast and ovarian cancer risk. BRCA1-mutant tumors are high-grade, recurrent, and often become resistant to standard therapies. Herein, we performed a targeted CRISPR-Cas9 screen and identified MEPCE, a methylphosphate capping enzyme, as a synthetic lethal interactor of BRCA1. Mechanistically, we demonstrate that depletion of MEPCE in a BRCA1-deficient setting led to dysregulated RNA polymerase II (RNAPII) promoter-proximal pausing, R-loop accumulation, and replication stress, contributing to transcription-replication collisions. These collisions compromise genomic integrity resulting in loss of viability of BRCA1-deficient cells. We also extend these findings to another RNAPII-regulating factor, PAF1. This study identifies a new class of synthetic lethal partners of BRCA1 that exploit the RNAPII pausing regulation and highlight the untapped potential of transcription-replication collision-inducing factors as unique potential therapeutic targets for treating cancers associated with BRCA1 mutations.

The transcription factor BTB and CNC homology 1(BACH1) has been linked to coronary artery disease risk by human genome-wide association studies, but little is known about the role of BACH1 in vascular smooth muscle cell (VSMC) phenotype switching and neointima formation following vascular injury. Therefore, this study aims to explore the role of BACH1 in vascular remodeling and its underlying mechanisms. BACH1 was highly expressed in human atherosclerotic plaques and has high transcriptional factor activity in VSMCs of human atherosclerotic arteries. VSMC-specific loss of Bach1 in mice inhibited the transformation of VSMC from contractile to synthetic phenotype and VSMC proliferation and attenuated the neointimal hyperplasia induced by wire injury. Mechanistically, BACH1 suppressed chromatin accessibility at the promoters of VSMC marker genes via recruiting histone methyltransferase G9a and cofactor YAP and maintaining the H3K9me2 state, thereby repressing VSMC marker genes expression in human aortic smooth muscle cells (HASMCs). BACH1-induced repression of VSMC marker genes was abolished by the silencing of G9a or YAP. Thus, these findings demonstrate a crucial regulatory role of BACH1 in VSMC phenotypic transition and vascular homeostasis and shed light on potential future protective vascular disease intervention via manipulation of BACH1.

YAP, the key protein effector of the Hippo pathway, is a transcriptional co-activator that controls the expression of cell cycle genes, promotes cell growth and proliferation and regulates organ size. YAP modulates gene transcription by binding to distal enhancers, but the mechanisms of gene regulation by YAP-bound enhancers remain poorly understood. Here we show that constitutive active YAP5SA leads to widespread changes in chromatin accessibility in untransformed MCF10A cells. Newly accessible regions include YAP-bound enhancers that mediate activation of cycle genes regulated by the Myb-MuvB (MMB) complex. By CRISPR-interference we identify a role for YAP-bound enhancers in phosphorylation of Pol II at Ser5 at MMB-regulated promoters, extending previously published studies that suggested YAP primarily regulates the pause-release step and transcriptional elongation. YAP5SA also leads to less accessible 'closed' chromatin regions, which are not directly YAP-bound but which contain binding motifs for the p53 family of transcription factors. Diminished accessibility at these regions is, at least in part, a consequence of reduced expression and chromatin-binding of the p53 family member ΔNp63 resulting in downregulation of ΔNp63-target genes and promoting YAP-mediated cell migration. In summary, our studies uncover changes in chromatin accessibility and activity that contribute to the oncogenic activities of YAP.

Human papillomavirus (HPV) integration is a critical step in cervical cancer development; however, the oncogenic mechanism at the genome-wide transcriptional level is still poorly understood. In this study, we employed integrative analysis on multi-omics data of six HPV-positive and three HPV-negative cell lines. Through HPV integration detection, super-enhancer (SE) identification, SE-associated gene expression and extrachromosomal DNA (ecDNA) investigation, we aimed to explore the genome-wide transcriptional influence of HPV integration. We identified seven high-ranking cellular SEs generated by HPV integration in total (the HPV breakpoint-induced cellular SEs, BP-cSEs), leading to intra-chromosomal and inter-chromosomal regulation of chromosomal genes. The pathway analysis revealed that the dysregulated chromosomal genes were correlated to cancer-related pathways. Importantly, we demonstrated that BP-cSEs existed in the HPV-human hybrid ecDNAs, explaining the above transcriptional alterations. Our results suggest that HPV integration generates cellular SEs that function as ecDNA to regulate unconstrained transcription, expanding the tumorigenic mechanism of HPV integration and providing insights for developing new diagnostic and therapeutic strategies.

Cryogenic electron microscopy (cryo-EM) is a promising method for characterizing the structure of larger RNA structures and complexes. However, the structure of individual aptamers is difficult to solve by cryo-EM due to their low molecular weight and a high signal-to-noise ratio. By placing RNA aptamers on larger RNA scaffolds, the contrast for cryo-EM can be increased to allow the determination of the tertiary structure of the aptamer. Here we use the RNA origami method to scaffold two fluorescent aptamers (Broccoli and Pepper) in close proximity and show that their cognate fluorophores serve as donor and acceptor for FRET. Next, we use cryo-EM to characterize the structure of the RNA origami with the two aptamers to a resolution of 4.4 Å. By characterizing the aptamers with and without ligand, we identify two distinct modes of ligand binding, which are further supported by selective chemical probing. 3D variability analysis of the cryo-EM data show that the relative position between the two bound fluorophores on the origami fluctuate by only 3.5 Å. Our results demonstrate a general approach for using RNA origami scaffolds for characterizing small RNA motifs by cryo-EM and for positioning functional RNA motifs with high spatial precision.

A panel of unnatural base pairs is developed to expand genetic alphabets. One or more unnatural base pairs (UBPs) can be inserted to enlarge the capacity, diversity, and functionality of canonical DNA, so monitoring the multiple-UBPs-containing DNA by simple and convenient approaches is essential. Herein, we report a bridge-base approach to repurpose the capability of determining TPT3-NaM UBPs. The success of this approach depends on the design of isoTAT that can simultaneously pair with NaM and G as a bridge base, as well as the discovering of the transformation of NaM to A in absence of its complementary base. TPT3-NaM can be transferred to C-G or A-T by simple PCR assays with high read-through ratios and low sequence-dependent properties, permitting for the first time to dually locate the multiple sites of TPT3-NaM pairs. Then we show the unprecedented capacity of this approach to trace accurate changes and retention ratios of multiple TPT3-NaM UPBs during in vivo replications. In addition, the method can also be applied to identify multiple-site DNA lesions, transferring TPT3-NaM makers to different natural bases. Taken together, our work presents the first general and convenient approach capable of locating, tracing, and sequencing site- and number-unlimited TPT3-NaM pairs.

Genome-encoded antibiotic resistance (ARE) ATP-binding cassette (ABC) proteins of the F subfamily (ARE-ABCFs) mediate intrinsic resistance in diverse Gram-positive bacteria. The diversity of chromosomally-encoded ARE-ABCFs is far from being fully experimentally explored. Here we characterise phylogenetically diverse genome-encoded ABCFs from Actinomycetia (Ard1 from Streptomyces capreolus, producer of the nucleoside antibiotic A201A), Bacilli (VmlR2 from soil bacterium Neobacillus vireti) and Clostridia (CplR from Clostridium perfringens, Clostridium sporogenes and Clostridioides difficile). We demonstrate that Ard1 is a narrow spectrum ARE-ABCF that specifically mediates self-resistance against nucleoside antibiotics. The single-particle cryo-EM structure of a VmlR2-ribosome complex allows us to rationalise the resistance spectrum of this ARE-ABCF that is equipped with an unusually long antibiotic resistance determinant (ARD) subdomain. We show that CplR contributes to intrinsic pleuromutilin, lincosamide and streptogramin A resistance in Clostridioides, and demonstrate that C. difficile CplR (CDIF630_02847) synergises with the transposon-encoded 23S ribosomal RNA methyltransferase Erm to grant high levels of antibiotic resistance to the C. difficile 630 clinical isolate. Finally, assisted by uORF4u, our novel tool for detection of upstream open reading frames, we dissect the translational attenuation mechanism that controls the induction of cplR expression upon an antibiotic challenge.

Genetic screening based on the clustered regularly interspaced palindromic repeat (CRISPR) system has been indicated to be a powerful tool for identifying regulatory genes or cis-elements. However, when applying CRISPR screens to pinpoint functional elements at particular loci, a large number of guide RNA (gRNA) spacers may be required to achieve saturated coverage. Here, we present a controlled template-dependent elongation (CTDE) method relying on reversible terminators to synthesize gRNA libraries with genomic regions of interest. By applying this approach to H3K4me3 chromatin immunoprecipitation (ChIP)-derived DNA of mammalian cells, mega-sized gRNA libraries were synthesized in a tissue-specific manner, with which we conducted screening experiments to annotate essential sites for cell proliferation. Additionally, we confirmed that an essential site within the intron of LINC00339 regulates its own mRNA and that LINC00339 is a novel regulator of the cell cycle that maintains HepG2 proliferation. The CTDE method has the potential to be automated with high efficiency at low cost, and will be widely used to identify functional elements in mammalian genomes.

The pandemic caused by SARS-CoV-2 has called for concerted efforts to generate new insights into the biology of betacoronaviruses to inform drug screening and development. Here, we establish a workflow to determine the RNA recognition and druggability of the nucleocapsid N-protein of SARS-CoV-2, a highly abundant protein crucial for the viral life cycle. We use a synergistic method that combines NMR spectroscopy and protein-RNA cross-linking coupled to mass spectrometry to quickly determine the RNA binding of two RNA recognition domains of the N-protein. Finally, we explore the druggability of these domains by performing an NMR fragment screening. This workflow identified small molecule chemotypes that bind to RNA binding interfaces and that have promising properties for further fragment expansion and drug development.

Understanding the pathological impact of non-coding genetic variation is a major challenge in medical genetics. Accumulating evidences indicate that a significant fraction of genetic alterations, including structural variants (SVs), can cause human disease by altering the function of non-coding regulatory elements, such as enhancers. In the case of SVs, described pathomechanisms include changes in enhancer dosage and long-range enhancer-gene communication. However, there is still a clear gap between the need to predict and interpret the medical impact of non-coding variants, and the existence of tools to properly perform these tasks. To reduce this gap, we have developed POSTRE (Prediction Of STRuctural variant Effects), a computational tool to predict the pathogenicity of SVs implicated in a broad range of human congenital disorders. By considering disease-relevant cellular contexts, POSTRE identifies SVs with either coding or long-range pathological consequences with high specificity and sensitivity. Furthermore, POSTRE not only identifies pathogenic SVs, but also predicts the disease-causative genes and the underlying pathological mechanism (e.g, gene deletion, enhancer disconnection, enhancer adoption, etc.). POSTRE is available at https://github.com/vicsanga/Postre.

Numerous viruses utilize essential long-range RNA-RNA genome interactions, specifically flaviviruses. Using Japanese encephalitis virus (JEV) as a model system, we computationally predicted and then biophysically validated and characterized its long-range RNA-RNA genomic interaction. Using multiple RNA computation assessment programs, we determine the primary RNA-RNA interacting site among JEV isolates and numerous related viruses. Following in vitro transcription of RNA, we provide, for the first time, characterization of an RNA-RNA interaction using size-exclusion chromatography coupled with multi-angle light scattering and analytical ultracentrifugation. Next, we demonstrate that the 5' and 3' terminal regions of JEV interact with nM affinity using microscale thermophoresis, and this affinity is significantly reduced when the conserved cyclization sequence is not present. Furthermore, we perform computational kinetic analyses validating the cyclization sequence as the primary driver of this RNA-RNA interaction. Finally, we examined the 3D structure of the interaction using small-angle X-ray scattering, revealing a flexible yet stable interaction. This pathway can be adapted and utilized to study various viral and human long-non-coding RNA-RNA interactions and determine their binding affinities, a critical pharmacological property of designing potential therapeutics.

Kinetoplastid parasites, such as trypanosomes or leishmania, rely on RNA-templated RNA editing to mature mitochondrial cryptic pre-mRNAs into functional protein-coding transcripts. Processive pan-editing of multiple editing blocks within a single transcript is dependent on the 20-subunit RNA editing substrate binding complex (RESC) that serves as a platform to orchestrate the interactions between pre-mRNA, guide RNAs (gRNAs), the catalytic RNA editing complex (RECC), and a set of RNA helicases. Due to the lack of molecular structures and biochemical studies with purified components, neither the spacio-temporal interplay of these factors nor the selection mechanism for the different RNA components is understood. Here we report the cryo-EM structure of Trypanosoma brucei RESC1-RESC2, a central hub module of the RESC complex. The structure reveals that RESC1 and RESC2 form an obligatory domain-swapped dimer. Although the tertiary structures of both subunits closely resemble each other, only RESC2 selectively binds 5'-triphosphate-nucleosides, a defining characteristic of gRNAs. We therefore propose RESC2 as the protective 5'-end binding site for gRNAs within the RESC complex. Overall, our structure provides a starting point for the study of the assembly and function of larger RNA-bound kinetoplast RNA editing modules and might aid in the design of anti-parasite drugs.

N6-methyladenosine (m6A) is the most prevalent RNA modification in eukaryotic mRNAs. Currently available detection methods for locus-specific m6A marks rely on RT-qPCR, radioactive methods, or high-throughput sequencing. Here, we develop a non-qPCR, ultrasensitive, isothermal, and naked-eye visible method for m6A detection based on rolling circle amplification (RCA) and loop-mediated isothermal amplification (LAMP), named m6A-Rol-LAMP, to verify putative m6A sites in transcripts obtained from the high-throughput data. When padlock probes hybridize to the potential m6A sites on targets, they are converted to circular form by DNA ligase in the absence of m6A modification, while m6A modification hinders the sealing of padlock probes. Subsequently, Bst DNA polymerase-mediated RCA and LAMP allow the amplification of the circular padlock probe to achieve the locus-specific detection of m6A. Following optimization and validation, m6A-Rol-LAMP can ultra-sensitively and quantitatively determine the existence of m6A modification on a specific target site as low as 100 amol under isothermal conditions. Detections of m6A can be performed on rRNA, mRNA, lincRNA, lncRNA and pre-miRNA from biological samples with naked-eye observations after dye incubation. Together, we provide a powerful tool for locus-specific detection of m6A, which can simply, quickly, sensitively, specifically, and visually determine putative m6A modification on RNA.

Long noncoding RNAs (lncRNAs) are >200 nt RNA transcripts without protein-coding potential. LncRNAs can be categorized into intergenic, intronic, bidirectional, sense, and antisense lncRNAs based on the genomic localization to nearby protein-coding genes. The current CRISPR-based lncRNA knockout strategy works efficiently for lncRNAs distant from the protein-coding gene, whereas it causes genomic perturbance inevitably due to technical limitations. In this study, we introduce a novel lncRNA knockout strategy, BESST, by deleting the genomic DNA fragment from the branch point to the 3' splicing site in the last intron of the target lncRNA. The BESST knockout exhibited comparable or superior repressive efficiency to RNA silencing or conventional promoter-exon1 deletion. Significantly, the BESST knockout strategy minimized the intervention of adjacent/overlap protein-coding genes by removing an average of ∼130 bp from genomic DNA. Our data also found that the BESST knockout strategy causes lncRNA nuclear retention, resulting in decapping and deadenylation of the lncRNA poly(A) tail. Further study revealed that PABPN1 is essential for the BESST-mediated decay and subsequent poly(A) deadenylation and decapping. Together, the BESST knockout strategy provides a versatile tool for investigating gene function by generating knockout cells or animals with high specificity and efficiency.