CRISPR-Cas12a is an RNA-guided, programmable genome editing enzyme found within bacterial adaptive immune pathways. Unlike CRISPR-Cas9, Cas12a uses only a single catalytic site to both cleave target double-stranded DNA (dsDNA) (cis-activity) and indiscriminately degrade single-stranded DNA (ssDNA) (trans-activity). To investigate how the relative potency of cis- versus trans-DNase activity affects Cas12a-mediated genome editing, we first used structure-guided engineering to generate variants of Lachnospiraceae bacterium Cas12a that selectively disrupt trans-activity. The resulting engineered mutant with the biggest differential between cis- and trans-DNase activity in vitro showed minimal genome editing activity in human cells, motivating a second set of experiments using directed evolution to generate additional mutants with robust genome editing activity. Notably, these engineered and evolved mutants had enhanced ability to induce homology-directed repair (HDR) editing by 2-18-fold compared to wild-type Cas12a when using HDR donors containing mismatches with crRNA at the PAM-distal region. Finally, a site-specific reversion mutation produced improved Cas12a (iCas12a) variants with superior genome editing efficiency at genomic sites that are difficult to edit using wild-type Cas12a. This strategy establishes a pipeline for creating improved genome editing tools by combining structural insights with randomization and selection. The available structures of other CRISPR-Cas enzymes will enable this strategy to be applied to improve the efficacy of other genome-editing proteins.

Telomeres terminate with a 50-300 bases long single-stranded G-rich overhang, which can be misrecognized as a DNA damage repair site. Shelterin plays critical roles in maintaining and protecting telomere ends by regulating access of various physiological agents to telomeric DNA, but the underlying mechanism is not well understood. Here, we measure how shelterin affects the accessibility of long telomeric overhangs by monitoring transient binding events of a short complementary peptide nucleic acid (PNA) probe using FRET-PAINT in vitro. We observed that the POT1 subunit of shelterin reduces the accessibility of the PNA probe by ∼2.5-fold, indicating that POT1 effectively binds to and protects otherwise exposed telomeric sequences. In comparison, a four-component shelterin stabilizes POT1 binding to the overhang by tethering POT1 to the double-stranded telomeric DNA and reduces the accessibility of telomeric overhangs by ∼5-fold. This enhanced protection suggests shelterin restructures the junction between single and double-stranded telomere, which is otherwise the most accessible part of the telomeric overhang.

Genetic interventions on microbiomes, for clinical or biotechnological purposes, remain challenging. Conjugation-based delivery of genetic cargo is still unspecific and limited by low conjugation rates. Here we report an approach to overcome these problems, based on a synthetic bacterial adhesion system. Mating assemblers consist on a synthetic adhesion formed by the expression on the surface of donor and target cells of specific nanobodies (Nb) and their cognate antigen (Ag). The Nb-Ag bridge increased 1-3 logs transfer of a variety of plasmids, especially in liquid media, confirming that cell-cell docking is a main determinant limiting mating efficiency. Synthetic cell-to-cell adhesion allows efficient conjugation to targeted recipients, enhancing delivery of desired genes to a predefined subset of prey species, or even specific pathogenic strains such as enterohemorrhagic Escherichia coli (EHEC), within a bacterial community. The synthetic conjugation enhancer presented here optimizes plasmid delivery by selecting the target hosts with high selectivity.

Pif1 family 5' → 3' DNA helicases are important for replication fork progression and genome stability. The budding yeast Saccharomyces cerevisiae encodes two Pif1 family helicases, Rrm3 and Pif1, both of which are multi-functional. Here we describe novel functions for Rrm3 in promoting mutation avoidance during DNA replication. We show that loss of RRM3 results in elevated spontaneous mutations made by DNA polymerases Pols ϵ and δ, which are subject to DNA mismatch repair. The absence of RRM3 also causes higher mutagenesis by the fourth B-family DNA polymerase Pol ζ. By genome-wide analysis, we show that the mutational consequences due to loss of RRM3 vary depending on the genomic locus. Rrm3 promotes the accuracy of DNA replication by Pols ϵ and δ across the genome, and it is particularly important for preventing Pol ζ-dependent mutagenesis at tRNA genes. In addition, mutation avoidance by Rrm3 depends on its helicase activity, and Pif1 serves as a backup for Rrm3 in suppressing mutagenesis. We present evidence that the sole human Pif1 family helicase in human cells likely also promotes replication fidelity, suggesting that a role for Pif1 family helicases in mutation avoidance may be evolutionarily conserved, a possible underlying mechanism for its potential tumor-suppressor function.

The type III-E CRISPR-Cas systems are newly identified adaptive immune systems in prokaryotes that use a single Cas7-11 protein to specifically cleave target RNA. Cas7-11 could associate with Csx29, a putative caspase-like protein encoded by the gene frequently found in the type III-E loci, suggesting a functional linkage between the RNase and protease activities in type III-E systems. Here, we demonstrated that target RNA recognition would stimulate the proteolytic activity of Csx29, and protein Csx30 is the endogenous substrate. More interestingly, while the cognate target RNA recognition would activate Csx29, non-cognate target RNA with the complementary 3' anti-tag sequence inhibits the enzymatic activity. Csx30 could bind to the sigma factor RpoE, which may initiate the stress response after proteolytic cleavage. Combined with biochemical and structural studies, we have elucidated the mechanisms underlying the target RNA-guided proteolytic activity of Csx29. Our work will guide further developments leveraging this simple RNA targeting system for RNA and protein-related applications.

Splicing factors play an essential role in regulation of alternative pre-mRNA splicing. While much progress has been made in delineating the mechanisms of the splicing machinery, the identity of signal transduction pathways and upstream factors that regulate splicing factor activity is largely unknown. A major challenge in the discovery of upstream regulatory factors of pre-mRNA splicing is the scarcity of functional genomics screening methods to monitor splicing outcomes of endogenous genes. Here, we have developed HiFENS (high throughput FISH detection of endogenous splicing isoforms), a high-throughput imaging assay based on hybridization chain reaction (HCR) and used HiFENS to screen for cellular factors that regulate alternative splicing of endogenous genes. We demonstrate optimized detection with high specificity of endogenous splicing isoforms and multiplexing of probes for accurate detection of splicing outcomes with single cell resolution. As proof-of-principle, we perform an RNAi screen of 702 human kinases and identify potential candidate upstream splicing regulators of the FGFR2 gene. HiFENS should be a useful tool for the unbiased delineation of cellular pathways involved in alternative splicing regulation.

Recent advances in spatial transcriptomics (ST) have brought unprecedented opportunities to understand tissue organization and function in spatial context. However, it is still challenging to precisely dissect spatial domains with similar gene expression and histology in situ. Here, we present DeepST, an accurate and universal deep learning framework to identify spatial domains, which performs better than the existing state-of-the-art methods on benchmarking datasets of the human dorsolateral prefrontal cortex. Further testing on a breast cancer ST dataset, we showed that DeepST can dissect spatial domains in cancer tissue at a finer scale. Moreover, DeepST can achieve not only effective batch integration of ST data generated from multiple batches or different technologies, but also expandable capabilities for processing other spatial omics data. Together, our results demonstrate that DeepST has the exceptional capacity for identifying spatial domains, making it a desirable tool to gain novel insights from ST studies.

Although CRISPR-Cas12a [clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein 12a] combining pre-amplification technology has the advantage of high sensitivity in biosensing, its generality and specificity are insufficient, which greatly restrains its application range. Here, we discovered a new targeting substrate for LbaCas12a (Lachnospiraceae bacterium Cas12a), namely double-stranded DNA (dsDNA) with a sticky-end region (PAM-SE+ dsDNA). We discovered that CRISPR-Cas12a had special enzymatic properties for this substrate DNA, including the ability to recognize and cleave it without needing a protospacer adjacent motif (PAM) sequence and a high sensitivity to single-base mismatches in that substrate. Further mechanism studies revealed that guide RNA (gRNA) formed a triple-stranded flap structure with the substrate dsDNA. We also discovered the property of low-temperature activation of CRISPR-Cas12a and, by coupling with the unique DNA hybridization kinetics at low temperature, we constructed a complete workflow for low-abundance point mutation detection in real samples, which was fast, convenient and free of single-stranded DNA (ssDNA) transformation. The detection limits were 0.005-0.01% for synthesized strands and 0.01-0.05% for plasmid genomic DNA, and the mutation abundances provided by our system for 28 clinical samples were in accordance with next-generation sequencing results. We believe that our work not only reveals novel information about the target recognition mechanism of the CRISPR-Cas12a system, but also greatly broadens its application scenarios.

UV-DDB is a DNA damage recognition protein recently discovered to participate in the removal of 8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxoG) by stimulating multiple steps of base excision repair (BER). In this study, we examined whether UV-DDB has a wider role in BER besides oxidized bases and found it has specificity for two known DNA substrates of alkyladenine glycosylase (AAG)/N-methylpurine DNA glycosylase (MPG): 1, N6-ethenoadenine (ϵA) and hypoxanthine. Gel mobility shift assays show that UV-DDB recognizes these two lesions 4-5 times better than non-damaged DNA. Biochemical studies indicated that UV-DDB stimulated AAG activity on both substrates by 4- to 5-fold. Native gels indicated UV-DDB forms a transient complex with AAG to help facilitate release of AAG from the abasic site product. Single molecule experiments confirmed the interaction and showed that UV-DDB can act to displace AAG from abasic sites. Cells when treated with methyl methanesulfonate resulted in foci containing AAG and UV-DDB that developed over the course of several hours after treatment. While colocalization did not reach 100%, foci containing AAG and UV-DDB reached a maximum at three hours post treatment. Together these data indicate that UV-DDB plays an important role in facilitating the repair of AAG substrates.

The Hippo-TEAD pathway regulates cellular proliferation and function. The existing paradigm is that TEAD co-activators, YAP and TAZ, and co-repressor, VGLL4, bind to the pocket region of TEAD1 to enable transcriptional activation or repressive function. Here we demonstrate a pocket-independent transcription repression mechanism whereby TEAD1 controls cell proliferation in both non-malignant mature differentiated cells and in malignant cell models. TEAD1 overexpression can repress tumor cell proliferation in distinct cancer cell lines. In pancreatic β cells, conditional knockout of TEAD1 led to a cell-autonomous increase in proliferation. Genome-wide analysis of TEAD1 functional targets via transcriptomic profiling and cistromic analysis revealed distinct modes of target genes, with one class of targets directly repressed by TEAD1. We further demonstrate that TEAD1 controls target gene transcription in a motif-dependent and orientation-independent manner. Mechanistically, we show that TEAD1 has a pocket region-independent, direct repressive function via interfering with RNA polymerase II (POLII) binding to target promoters. Our study reveals that TEAD1 target genes constitute a mutually restricted regulatory loop to control cell proliferation and uncovers a novel direct repression mechanism involved in its transcriptional control that could be leveraged in future studies to modulate cell proliferation in tumors and potentially enhance the proliferation of normal mature cells.

The model plant Arabidopsis thaliana encodes as many as ten Argonaute proteins (AGO1-10) with different functions. Each AGO selectively loads a set of small RNAs by recognizing their length and 5' nucleotide identity to properly regulate target genes. Previous studies showed that AGO4 and AGO6, key factors in DNA methylation, incorporate 24-nt small-interfering RNAs with 5' adenine (24A siRNAs). However, it has been unclear how these AGOs specifically load 24A siRNAs. Here, we biochemically investigated the siRNA preference of AGO4, AGO6 and their chimeric mutants. We found that AGO4 and AGO6 use distinct mechanisms to preferentially load 24A siRNAs. Moreover, we showed that the 5' A specificity of AGO4 and AGO6 is not determined by the previously known nucleotide specificity loop in the MID domain but rather by the coordination of the MID and PIWI domains. These findings advance our mechanistic understanding of how small RNAs are accurately sorted into different AGO proteins in plants.

Rearrangement hot spot (Rhs) proteins are members of the broad family of polymorphic toxins. Polymorphic toxins are modular proteins composed of an N-terminal region that specifies their mode of secretion into the medium or into the target cell, a central delivery module, and a C-terminal domain that has toxic activity. Here, we structurally and functionally characterize the C-terminal toxic domain of the antibacterial Rhsmain protein, TreTu, which is delivered by the type VI secretion system of Salmonella enterica Typhimurium. We show that this domain adopts an ADP-ribosyltransferase fold and inhibits protein synthesis by transferring an ADP-ribose group from NAD+ to the elongation factor Tu (EF-Tu). This modification is specifically placed on the side chain of the conserved D21 residue located on the P-loop of the EF-Tu G-domain. Finally, we demonstrate that the TriTu immunity protein neutralizes TreTu activity by acting like a lid that closes the catalytic site and traps the NAD+.

Antibiotics have been widely used for plasmid-mediated cell engineering. However, continued use of antibiotics increases the metabolic burden, horizontal gene transfer risks, and biomanufacturing costs. There are limited approaches to maintaining multiple plasmids without antibiotics. Herein, we developed an inverter cascade using CRISPRi by building a plasmid containing a single guide RNA (sgRNA) landing pad (pSLiP); this inhibited host cell growth by repressing an essential cellular gene. Anti-sgRNAs on separate plasmids restored cell growth by blocking the expression of growth-inhibitory sgRNAs in pSLiP. We maintained three plasmids in Escherichia coli with a single antibiotic selective marker. To completely avoid antibiotic use and maintain the CRISPRi-based logic inverter cascade, we created a novel d-glutamate auxotrophic E. coli. This enabled the stable maintenance of the plasmid without antibiotics, enhanced the production of the terpenoid, (-)-α-bisabolol, and generation of an antibiotic-resistance gene-free plasmid. CRISPRi is therefore widely applicable in genetic circuits and may allow for antibiotic-free biomanufacturing.

DNA transposon systems are widely used in mammalian cells for genetic modification experiments, but their regulation remains poorly understood. We used biochemical and cell-based assays together with AlphaFold modeling and rational protein redesign to evaluate aspects of piggyBac transposition including the previously unexplained role of the transposase N-terminus and the need for asymmetric transposon ends for cellular activity. We found that phosphorylation at predicted casein kinase II sites in the transposase N-terminus inhibits transposition, most likely by preventing transposase-DNA interactions. Deletion of the region containing these sites releases inhibition thereby enhancing activity. We also found that the N-terminal domain promotes transposase dimerization in the absence of transposon DNA. When the N-terminus is deleted, the transposase gains the ability to carry out transposition using symmetric transposon left ends. This novel activity is also conferred by appending a second C-terminal domain. When combined, these modifications together result in a transposase that is highly active when symmetric transposon ends are used. Our results demonstrate that transposase N-terminal phosphorylation and the requirement for asymmetric transposon ends both negatively regulate piggyBac transposition in mammalian cells. These novel insights into the mechanism and structure of the piggyBac transposase expand its potential use for genomic applications.

The glucocorticoid receptor (GR) is a ubiquitously expressed transcription factor that controls metabolic and homeostatic processes essential for life. Although numerous crystal structures of the GR ligand-binding domain (GR-LBD) have been reported, the functional oligomeric state of the full-length receptor, which is essential for its transcriptional activity, remains disputed. Here we present five new crystal structures of agonist-bound GR-LBD, along with a thorough analysis of previous structural work. We identify four distinct homodimerization interfaces on the GR-LBD surface, which can associate into 20 topologically different homodimers. Biologically relevant homodimers were identified by studying a battery of GR point mutants including crosslinking assays in solution, quantitative fluorescence microscopy in living cells, and transcriptomic analyses. Our results highlight the relevance of non-canonical dimerization modes for GR, especially of contacts made by loop L1-3 residues such as Tyr545. Our work illustrates the unique flexibility of GR's LBD and suggests different dimeric conformations within cells. In addition, we unveil pathophysiologically relevant quaternary assemblies of the receptor with important implications for glucocorticoid action and drug design.

Single-stranded DNA-binding proteins (SSBs) interact with single-stranded DNA (ssDNA) to form filamentous structures with various degrees of cooperativity, as a result of intermolecular interactions between neighboring SSB subunits on ssDNA. However, it is still challenging to perform structural studies on SSB-ssDNA filaments at high resolution using the most studied SSB models, largely due to the intrinsic flexibility of these nucleoprotein complexes. In this study, HaLEF-3, an SSB protein from Helicoverpa armigera nucleopolyhedrovirus, was used for in vitro assembly of SSB-ssDNA filaments, which were structurally studied at atomic resolution using cryo-electron microscopy. Combined with the crystal structure of ssDNA-free HaLEF-3 octamers, our results revealed that the three-dimensional rearrangement of HaLEF-3 induced by an internal hinge-bending movement is essential for the formation of helical SSB-ssDNA complexes, while the contacting interface between adjacent HaLEF-3 subunits remains basically intact. We proposed a local cooperative SSB-ssDNA binding model, in which, triggered by exposure to oligonucleotides, HaLEF-3 molecules undergo ring-to-helix transition to initiate continuous SSB-SSB interactions along ssDNA. Unique structural features revealed by the assembly of HaLEF-3 on ssDNA suggest that HaLEF-3 may represent a new class of SSB.

Genotoxic agents, that are used in cancer therapy, elicit the reprogramming of the transcriptome of cancer cells. These changes reflect the cellular response to stress and underlie some of the mechanisms leading to drug resistance. Here, we profiled genome-wide changes in pre-mRNA splicing induced by cisplatin in breast cancer cells. Among the set of cisplatin-induced alternative splicing events we focused on COASY, a gene encoding a mitochondrial enzyme involved in coenzyme A biosynthesis. Treatment with cisplatin induces the production of a short isoform of COASY lacking exons 4 and 5, whose depletion impedes mitochondrial function and decreases sensitivity to cisplatin. We identified RBM39 as a major effector of the cisplatin-induced effect on COASY splicing. RBM39 also controls a genome-wide set of alternative splicing events partially overlapping with the cisplatin-mediated ones. Unexpectedly, inactivation of RBM39 in response to cisplatin involves its interaction with the AP-1 family transcription factor c-Jun that prevents RBM39 binding to pre-mRNA. Our findings therefore uncover a novel cisplatin-induced interaction between a splicing regulator and a transcription factor that has a global impact on alternative splicing and contributes to drug resistance.

Ribosome profiling (Ribo-seq) is a powerful method for the transcriptome-wide assessment of protein synthesis rates and the study of translational control mechanisms. Yet, Ribo-seq also has limitations. These include difficulties with the analysis of translation-modulating molecules such as antibiotics, which are often toxic or challenging to deliver into living cells. Here, we have developed in vitro Ribo-seq (INRI-seq), a cell-free method to analyze the translational landscape of a fully customizable synthetic transcriptome. Using Escherichia coli as an example, we show how INRI-seq can be used to analyze the translation initiation sites of a transcriptome of interest. We also study the global impact of direct translation inhibition by antisense peptide nucleic acid (PNA) to analyze PNA off-target effects. Overall, INRI-seq presents a scalable, sensitive method to study translation initiation in a transcriptome-wide manner without the potentially confounding effects of extracting ribosomes from living cells.

Understanding how modifications to the ribosome affect function has implications for studying ribosome biogenesis, building minimal cells, and repurposing ribosomes for synthetic biology. However, efforts to design sequence-modified ribosomes have been limited because point mutations in the ribosomal RNA (rRNA), especially in the catalytic active site (peptidyl transferase center; PTC), are often functionally detrimental. Moreover, methods for directed evolution of rRNA are constrained by practical considerations (e.g. library size). Here, to address these limitations, we developed a computational rRNA design approach for screening guided libraries of mutant ribosomes. Our method includes in silico library design and selection using a Rosetta stepwise Monte Carlo method (SWM), library construction and in vitro testing of combined ribosomal assembly and translation activity, and functional characterization in vivo. As a model, we apply our method to making modified ribosomes with mutant PTCs. We engineer ribosomes with as many as 30 mutations in their PTCs, highlighting previously unidentified epistatic interactions, and show that SWM helps identify sequences with beneficial phenotypes as compared to random library sequences. We further demonstrate that some variants improve cell growth in vivo, relative to wild type ribosomes. We anticipate that SWM design and selection may serve as a powerful tool for rRNA engineering.

Disruptive mutations in the chromodomain helicase DNA-binding protein 8 gene (CHD8) have been recurrently associated with autism spectrum disorders (ASDs). Here we investigated how chromatin reacts to CHD8 suppression by analyzing a panel of histone modifications in induced pluripotent stem cell-derived neural progenitors. CHD8 suppression led to significant reduction (47.82%) in histone H3K36me3 peaks at gene bodies, particularly impacting on transcriptional elongation chromatin states. H3K36me3 reduction specifically affects highly expressed, CHD8-bound genes and correlates with altered alternative splicing patterns of 462 genes implicated in 'regulation of RNA splicing' and 'mRNA catabolic process'. Mass spectrometry analysis uncovered a novel interaction between CHD8 and the splicing regulator heterogeneous nuclear ribonucleoprotein L (hnRNPL), providing the first mechanistic insights to explain the CHD8 suppression-derived splicing phenotype, partly implicating SETD2, a H3K36me3 methyltransferase. In summary, our results point toward broad molecular consequences of CHD8 suppression, entailing altered histone deposition/maintenance and RNA processing regulation as important regulatory processes in ASD.

Cancer cells achieve immortality by employing either homology-directed repair (HDR) or the telomerase enzyme to maintain telomeres. ALT (alternative lengthening of telomeres) refers to the subset of cancer cells that employ HDR. Many ALT features are conserved from yeast to human cells, with the yeast equivalent being referred to as survivors. The non-coding RNA TERRA, and its ability to form RNA-DNA hybrids, has been implicated in ALT/survivor maintenance by promoting HDR. It is not understood which telomeres in ALT/survivors engage in HDR, nor is it clear which telomeres upregulate TERRA. Using yeast survivors as a model for ALT, we demonstrate that HDR only occurs at telomeres when they become critically short. Moreover, TERRA levels steadily increase as telomeres shorten and decrease again following HDR-mediated recombination. We observe that survivors undergo cycles of senescence, in a similar manner to non-survivors following telomerase loss, which we refer to as survivor associated senescence (SAS). Similar to 'normal' senescence, we report that RNA-DNA hybrids slow the rate of SAS, likely through the elongation of critically short telomeres, however decreasing the rate of telomere shortening may contribute to this effect. In summary, TERRA RNA-DNA hybrids regulate telomere dysfunction-induced senescence before and after survivor formation.

Triplex nanostructures can be formed in vitro in the promoter region of DNA templates, and it is commonly accepted that these assemblies inhibit the transcription of the downstream genes. Herein, a proof of concept highlighting the possibility of the up- or downregulation of RNA transcription is presented. Hybrid DNA-RNA triplex nanostructures were rationally designed to produce bacterial transcription units with switchable promoters. The rate of RNA production was measured using the signal of a transcribed fluorescent RNA aptamer (i.e. Broccoli). Indeed, several designed bacterial promoters showed the ability of induced transcriptional inhibition, while other properly tailored sequences demonstrated switchable enhancement of transcriptional activity, representing an unprecedented feature to date. The use of RNA-regulated transcription units and fluorescent RNA aptamers as readouts will allow the realization of biocomputation circuits characterized by a strongly reduced set of components. Triplex forming RNA oligonucleotides are proposed as smart tools for transcriptional modulation and represent an alternative to current methods for producing logic gates using protein-based components.

Single-stranded DNA binding proteins (SSBs) avidly bind ssDNA and yet enzymes that need to act during DNA replication and repair are not generally impeded by SSB, and are often stimulated by SSB. Here, the effects of Escherichia coli SSB on the activities of the DNA polymerase processivity clamp loader were investigated. SSB enhances binding of the clamp loader to DNA by increasing the lifetime on DNA. Clamp loading was measured on DNA substrates that differed in length of ssDNA overhangs to permit SSB binding in different binding modes. Even though SSB binds DNA adjacent to single-stranded/double-stranded DNA junctions where clamps are loaded, the rate of clamp loading on DNA was not affected by SSB on any of the DNA substrates. Direct measurements of the relative timing of DNA-SSB remodeling and enzyme-DNA binding showed that the clamp loader rapidly remodels SSB on DNA such that SSB has little effect on DNA binding rates. However, when SSB was mutated to reduce protein-protein interactions with the clamp loader, clamp loading was inhibited by impeding binding of the clamp loader to DNA. Thus, protein-protein interactions between the clamp loader and SSB facilitate rapid DNA-SSB remodeling to allow rapid clamp loader-DNA binding and clamp loading.

Aggregation of the microtubule-associated protein tau characterizes tauopathies, including Alzheimer's disease and frontotemporal lobar degeneration (FTLD-Tau). Gene expression regulation of tau is complex and incompletely understood. Here we report that the human tau gene (MAPT) generates two circular RNAs (circRNAs) through backsplicing of exon 12 to either exon 7 (12→7 circRNA) or exon 10 (12→10 circRNA). Both circRNAs lack stop codons. The 12→7 circRNA contains one start codon and is translated in a rolling circle, generating a protein consisting of multimers of the microtubule-binding repeats R1-R4. For the 12→10 circRNA, a start codon can be introduced by two FTLD-Tau mutations, generating a protein consisting of multimers of the microtubule-binding repeats R2-R4, suggesting that mutations causing FTLD may act in part through tau circRNAs. Adenosine to inosine RNA editing dramatically increases translation of circRNAs and, in the 12→10 circRNA, RNA editing generates a translational start codon by changing AUA to AUI. Circular tau proteins self-aggregate and promote aggregation of linear tau proteins. Our data indicate that adenosine to inosine RNA editing initiates translation of human circular tau RNAs, which may contribute to tauopathies.

Sam68, also known as KHDRBS1, is a member of the STAR family of proteins that directly link signal transduction with post-transcriptional gene regulation. Sam68 controls the alternative splicing of many oncogenic proteins and its role is modulated by post-translational modifications, including serine/threonine phosphorylation, that differ at various stages of the cell cycle. However, the molecular basis and mechanisms of these modulations remain largely unknown. Here, we combined mass spectrometry, nuclear magnetic resonance spectroscopy and cell biology techniques to provide a comprehensive post-translational modification mapping of Sam68 at different stages of the cell cycle in HEK293 and HCT116 cells. We established that Sam68 is specifically phosphorylated at T33 and T317 by Cdk1, and demonstrated that these phosphorylation events reduce the binding of Sam68 to RNA, control its cellular localization and reduce its alternative splicing activity, leading to a reduction in the induction of apoptosis and an increase in the proliferation of HCT116 cells.

Friedreich's ataxia is an incurable disease caused by frataxin (FXN) protein deficiency, which is mostly induced by GAA repeat expansion in intron 1 of the FXN gene. Here, we identified antisense oligonucleotides (ASOs), complementary to two regions within the first intron of FXN pre-mRNA, which could increase FXN mRNA by ∼2-fold in patient fibroblasts. The increase in FXN mRNA was confirmed by the identification of multiple overlapping FXN-activating ASOs at each region, two independent RNA quantification assays, and normalization by multiple housekeeping genes. Experiments on cells with the ASO-binding sites deleted indicate that the ASO-induced FXN activation was driven by indirect effects. RNA sequencing analyses showed that the two ASOs induced similar transcriptome-wide changes, which did not resemble the transcriptome of wild-type cells. This RNA-seq analysis did not identify directly base-paired off-target genes shared across ASOs. Mismatch studies identified two guanosine-rich motifs (CCGG and G4) within the ASOs that were required for FXN activation. The phosphorodiamidate morpholino oligomer analogs of our ASOs did not activate FXN, pointing to a PS-backbone-mediated effect. Our study demonstrates the importance of multiple, detailed control experiments and target validation in oligonucleotide studies employing novel mechanisms such as gene activation.

Chromosome axis-associated HORMA domain proteins (HORMADs), e.g. ASY1 in Arabidopsis, are crucial for meiotic recombination. ASY1, as other HORMADs, is assembled on the axis at early meiosis and depleted when homologous chromosomes synapse. Puzzlingly, both processes are catalyzed by AAA+ ATPase PCH2 together with its cofactor COMET. Here, we show that the ASY1 remodeling complex is temporally and spatially differently assembled. While PCH2 and COMET appear to directly interact in the cytoplasm in early meiosis, PCH2 is recruited by the transverse filament protein ZYP1 and brought to the ASY1-bound COMET assuring the timely removal of ASY1 during chromosome synapsis. Since we found that the PCH2 homolog TRIP13 also binds to the ZYP1 homolog SYCP1 in mouse, we postulate that this mechanism is conserved among eukaryotes. Deleting the PCH2 binding site of ZYP1 led to a failure of ASY1 removal. Interestingly, the placement of one obligatory crossover per homologous chromosome pair, compromised by ZYP1 depletion, is largely restored in this separation-of-function zyp1 allele suggesting that crossover assurance is promoted by synapsis. In contrast, this zyp1 allele, similar to the zyp1 null mutant, showed elevated type I crossover numbers indicating that PCH2-mediated eviction of ASY1 from the axis restricts crossover formation.

Control of mRNA translation is key for stress responses. Translation initiation is usually rate-limiting and, in eukaryotes, involves mRNA scanning by the small ribosomal subunit. Despite its importance, many aspects of translation in vivo have not been explored fully, especially at the transcriptome-wide level. A recent method termed translation-complex profiling (TCP-seq) allows transcriptome-wide views of scanning ribosomal subunits. We applied TCP-seq to nutritional stress in the fission yeast Schizosaccharomyces pombe. At initiation sites, we observed multiple complexes resembling those of mammals, and consistent with queuing of scanning subunits. In 5' UTRs, small subunit accumulations were common and may reflect impediments to scanning. A key mediator of stress responses in S. pombe is the Fil1 transcription factor, which is regulated translationally by a poorly-understood mechanism involving upstream Open Reading Frames (uORFs). TCP-seq data of fil1 shows that stress allows scanning subunits to by-pass specific uORFs and reach the fil1 coding sequence. The integration of these observations with reporter assays revealed that fil1 translational control is mediated by a combination of scanning reinitiation-repressive and permissive uORFs, and establishes fil1 as a model for uORF-mediated translational control. Altogether, our transcriptome-wide study reveals general and gene-specific features of translation in a model eukaryote.

Mitochondrial RNA metabolism is suggested to occur in identified compartmentalized foci, i.e. mitochondrial RNA granules (MRGs). Mitochondrial aminoacyl-tRNA synthetases (mito aaRSs) catalyze tRNA charging and are key components in mitochondrial gene expression. Mutations of mito aaRSs are associated with various human disorders. However, the suborganelle distribution, interaction network and regulatory mechanism of mito aaRSs remain largely unknown. Here, we found that all mito aaRSs partly colocalize with MRG, and this colocalization is likely facilitated by tRNA-binding capacity. A fraction of human mitochondrial AlaRS (hmtAlaRS) and hmtSerRS formed a direct complex via interaction between catalytic domains in vivo. Aminoacylation activities of both hmtAlaRS and hmtSerRS were fine-tuned upon complex formation in vitro. We further established a full spectrum of interaction networks via immunoprecipitation and mass spectrometry for all mito aaRSs and discovered interactions between hmtSerRS and hmtAsnRS, between hmtSerRS and hmtTyrRS and between hmtThrRS and hmtArgRS. The activity of hmtTyrRS was also influenced by the presence of hmtSerRS. Notably, hmtSerRS utilized the same catalytic domain in mediating several interactions. Altogether, our results systematically analyzed the suborganelle localization and interaction network of mito aaRSs and discovered several mito aaRS-containing complexes, deepening our understanding of the functional and regulatory mechanisms of mito aaRSs.

Rpc31 is a subunit in the TFIIE-related Rpc82/34/31 heterotrimeric subcomplex of Saccharomyces cerevisiae RNA polymerase III (pol III). Structural analyses of pol III have indicated that the N-terminal region of Rpc31 anchors on Rpc82 and further interacts with the polymerase core and stalk subcomplex. However, structural and functional information for the C-terminal region of Rpc31 is sparse. We conducted a mutational analysis on Rpc31, which uncovered a functional peptide adjacent to the highly conserved Asp-Glu-rich acidic C-terminus. This C-terminal peptide region, termed 'pre-acidic', is important for optimal cell growth, tRNA synthesis, and stable association of Rpc31 in the pre-initiation complex (PIC). Our site-directed photo-cross-linking to map protein interactions within the PIC reveal that this pre-acidic region specifically targets Rpc34 during transcription initiation, but also interacts with the DNA entry surface in free pol III. Thus, we have uncovered a switchable Rpc31 C-terminal region that functions in an initiation-specific protein interaction for pol III transcription.

In cyanobacteria DNA supercoiling varies over the diurnal cycle and is integrated with temporal programs of transcription and replication. We manipulated DNA supercoiling in Synechocystis sp. PCC 6803 by CRISPRi-based knockdown of gyrase subunits and overexpression of topoisomerase I (TopoI). Cell division was blocked but cell growth continued in all strains. The small endogenous plasmids were only transiently relaxed, then became strongly supercoiled in the TopoI overexpression strain. Transcript abundances showed a pronounced 5'/3' gradient along transcription units, incl. the rRNA genes, in the gyrase knockdown strains. These observations are consistent with the basic tenets of the homeostasis and twin-domain models of supercoiling in bacteria. TopoI induction initially led to downregulation of G+C-rich and upregulation of A+T-rich genes. The transcriptional response quickly bifurcated into six groups which overlap with diurnally co-expressed gene groups. Each group shows distinct deviations from a common core promoter structure, where helically phased A-tracts are in phase with the transcription start site. Together, our data show that major co-expression groups (regulons) in Synechocystis all respond differentially to DNA supercoiling, and suggest to re-evaluate the long-standing question of the role of A-tracts in bacterial promoters.

The replicative DNA helicase translocates on single-stranded DNA to drive replication forks during chromosome replication. In most bacteria the ubiquitous replicative helicase, DnaB, co-evolved with the accessory subunit DciA, but how they function remains incompletely understood. Here, using the model bacterium Caulobacter crescentus, we demonstrate that DciA plays a prominent role in DNA replication fork maintenance. Cell cycle analyses using a synchronized Caulobacter cell population showed that cells devoid of DciA exhibit a severe delay in fork progression. Biochemical characterization revealed that the DnaB helicase in its default state forms a hexamer that inhibits self-loading onto single-stranded DNA. We found that upon binding to DciA, the DnaB hexamer undergoes conformational changes required for encircling single-stranded DNA, thereby establishing the replication fork. Further investigation of the functional structure of DciA revealed that the C-terminus of DciA includes conserved leucine residues responsible for DnaB binding and is essential for DciA in vivo functions. We propose that DciA stimulates loading of DnaB onto single strands through topological isomerization of the DnaB structure, thereby ensuring fork progression. Given that the DnaB-DciA modules are widespread among eubacterial species, our findings suggest that a common mechanism underlies chromosome replication.

The Notch pathway transmits signals between neighboring cells to elicit downstream transcriptional programs. Notch is a major regulator of cell fate specification, proliferation, and apoptosis, such that aberrant signaling leads to a pleiotropy of human diseases, including developmental disorders and cancers. The pathway signals through the transcription factor CSL (RBPJ in mammals), which forms an activation complex with the intracellular domain of the Notch receptor and the coactivator Mastermind. CSL can also function as a transcriptional repressor by forming complexes with one of several different corepressor proteins, such as FHL1 or SHARP in mammals and Hairless in Drosophila. Recently, we identified L3MBTL3 as a bona fide RBPJ-binding corepressor that recruits the repressive lysine demethylase LSD1/KDM1A to Notch target genes. Here, we define the RBPJ-interacting domain of L3MBTL3 and report the 2.06 Å crystal structure of the RBPJ-L3MBTL3-DNA complex. The structure reveals that L3MBTL3 interacts with RBPJ via an unusual binding motif compared to other RBPJ binding partners, which we comprehensively analyze with a series of structure-based mutants. We also show that these disruptive mutations affect RBPJ and L3MBTL3 function in cells, providing further insights into Notch mediated transcriptional regulation.

Bacteria can adapt in response to numerous stress conditions. One such stress condition is zinc depletion. The zinc-sensing transcription factor Zur regulates the way numerous bacterial species respond to severe changes in zinc availability. Under zinc sufficient conditions, Zn-loaded Zur (Zn2-Zur) is well-known to repress transcription of genes encoding zinc uptake transporters and paralogues of a few ribosomal proteins. Here, we report the discovery and mechanistic basis for the ability of Zur to up-regulate expression of the ribosomal protein L31 in response to zinc in E. coli. Through genetic mutations and reporter gene assays, we find that Zur achieves the up-regulation of L31 through a double repression cascade by which Zur first represses the transcription of L31p, a zinc-lacking paralogue of L31, which in turn represses the translation of L31. Mutational analyses show that translational repression by L31p requires an RNA hairpin structure within the l31 mRNA and involves the N-terminus of the L31p protein. This work uncovers a new genetic network that allows bacteria to respond to host-induced nutrient limiting conditions through a sophisticated ribosomal protein switching mechanism.

The development of advanced genetic tools is boosting microbial engineering which can potentially tackle wide-ranging challenges currently faced by our society. Here we present SURE editing, a multi-recombinase engineering rationale combining oligonucleotide recombineering with the selective capacity of antibiotic resistance via transient insertion of selector plasmids. We test this method in Mycoplasma pneumoniae, a bacterium with a very inefficient native recombination machinery. Using SURE editing, we can seamlessly generate, in a single step, a wide variety of genome modifications at high efficiencies, including the largest possible deletion of this genome (30 Kb) and the targeted complementation of essential genes in the deletion of a region of interest. Additional steps can be taken to remove the selector plasmid from the edited area, to obtain markerless or even scarless edits. Of note, SURE editing is compatible with different site-specific recombinases for mediating transient plasmid integration. This battery of selector plasmids can be used to select different edits, regardless of the target sequence, which significantly reduces the cloning load associated to genome engineering projects. Given the proven functionality in several microorganisms of the machinery behind the SURE editing logic, this method is likely to represent a valuable advance for the synthetic biology field.

Chromatin remodelers have been thought to be crucial in creating an accessible chromatin environment before transcription activation. However, it is still unclear how chromatin remodelers recognize and bind to the active regions. In this study, we found that chromatin remodelers SPLAYED (SYD) and BRAHMA (BRM) interact and co-occupy with Suppressor of Ty6-like (SPT6L), a core subunit of the transcription machinery, at thousands of the transcription start sites (TSS). The association of SYD and BRM to chromatin is dramatically reduced in spt6l and can be restored mainly by SPT6LΔtSH2, which binds to TSS in a RNA polymerase II (Pol II)-independent manner. Furthermore, SPT6L and SYD/BRM are involved in regulating the nucleosome and Pol II occupancy around TSS. The presence of SPT6L is sufficient to restore the association of the chromatin remodeler SYD to chromatin and maintain normal nucleosome occupancy. Our findings suggest that the two chromatin remodelers can form protein complexes with the core subunit of the transcription machinery and regulate nucleosome occupancy in the early transcription stage.

The development of single-cell RNA-seq has broadened the spectrum for biological research by providing a high-resolution analysis of cellular heterogeneity. However, the requirement for sophisticated devices for the compartmentalization of cells has limited its widespread applicability. Here, we develop Onepot-Seq, a device-free method, that harnesses the transient localization of mRNA after lysis to capture single-cell transcriptomes simultaneously in a continuous fluid medium. In mixed-species experiments, we obtained high-quality single-cell profiles. Further, cell type-specific poly(A)-conjugated antibodies allow Onepot-Seq to effectively capture target cells in complex populations. Chemical perturbations to cells can be profiled by Onepot-Seq at single-cell resolution. Onepot-Seq should allow routine transcriptional profiling at single-cell resolution, accelerating clinical and scientific discoveries in many fields of science.

Analysis of the protein coding transcriptome by the RNA sequencing requires either enrichment of the desired fraction of coding transcripts or depletion of the abundant non-coding fraction consisting mainly of rRNA. We propose an alternative mRNA enrichment strategy based on the RNA-binding properties of the human IFIT1, an antiviral protein recognizing cap 0 RNA. Here, we compare for Saccharomyces cerevisiae an IFIT1-based mRNA pull-down with yeast targeted rRNA depletion by the RiboMinus method. IFIT1-based RNA capture depletes rRNA more effectively, producing high quality RNA-seq data with an excellent coverage of the protein coding transcriptome, while depleting cap-less transcripts such as mitochondrial or some non-coding RNAs. We propose IFIT1 as a cost effective and versatile tool to prepare mRNA libraries for a variety of organisms with cap 0 mRNA ends, including diverse plants, fungi and eukaryotic microbes.

Heterochromatin protein 1α (HP1α) is a crucial element of chromatin organization. It has been proposed that HP1α functions through liquid-liquid phase separation (LLPS), which allows it to compact chromatin into transcriptionally repressed heterochromatin regions. In vitro, HP1α can undergo phase separation upon phosphorylation of its N-terminus extension (NTE) and/or through interactions with DNA and chromatin. Here, we combine computational and experimental approaches to elucidate the molecular interactions that drive these processes. In phosphorylation-driven LLPS, HP1α can exchange intradimer hinge-NTE interactions with interdimer contacts, which also leads to a structural change from a compacted to an extended HP1α dimer conformation. This process can be enhanced by the presence of positively charged HP1α peptide ligands and disrupted by the addition of negatively charged or neutral peptides. In DNA-driven LLPS, both positively and negatively charged peptide ligands can perturb phase separation. Our findings demonstrate the importance of electrostatic interactions in HP1α LLPS where binding partners can modulate the overall charge of the droplets and screen or enhance hinge region interactions through specific and non-specific effects. Our study illuminates the complex molecular framework that can fine-tune the properties of HP1α and that can contribute to heterochromatin regulation and function.

The mammalian transcriptome comprises a vast family of long noncoding (lnc)RNAs implicated in physiologic processes such as myogenesis, through which muscle forms during embryonic development and regenerates in the adult. However, the specific molecular mechanisms by which lncRNAs regulate human myogenesis are poorly understood. Here, we identified a novel muscle-specific lncRNA, lncFAM71E1-2:2 (lncFAM), which increased robustly during early human myogenesis. Overexpression of lncFAM promoted differentiation of human myoblasts into myotubes, while silencing lncFAM suppressed this process. As lncFAM resides in the nucleus, chromatin isolation by RNA purification followed by mass spectrometry (ChIRP-MS) analysis was employed to identify the molecular mechanisms whereby it might promote myogenesis. Analysis of lncFAM-interacting proteins revealed that lncFAM recruited the RNA-binding protein HNRNPL to the promoter of MYBPC2, in turn increasing MYBPC2 mRNA transcription and enhancing production of the myogenic protein MYBPC2. These results highlight a mechanism whereby a novel ribonucleoprotein complex, lncFAM-HNRNPL, elevates MYBPC2 expression transcriptionally to promote myogenesis.

Quick growth restart after upon encountering favourable environmental conditions is a major fitness contributor in natural environment. It is widely assumed that the time required to restart growth after nutritional upshift is determined by how long it takes for cells to synthesize enough ribosomes to produce the proteins required to reinitiate growth. Here we show that a reduction in the capacity to synthesize ribosomes by reducing number of ribosomal RNA (rRNA) operons (rrn) causes a longer transition from stationary phase to growth of Escherichia coli primarily due to high mortality rates. Cell death results from DNA replication blockage and massive DNA breakage at the sites of the remaining rrn operons that become overloaded with RNA polymerases (RNAPs). Mortality rates and growth restart duration can be reduced by preventing R-loop formation and improving DNA repair capacity. The same molecular mechanisms determine the duration of the recovery phase after ribosome-damaging stresses, such as antibiotics, exposure to bile salts or high temperature. Our study therefore suggests that a major function of rrn operon multiplicity is to ensure that individual rrn operons are not saturated by RNAPs, which can result in catastrophic chromosome replication failure and cell death during adaptation to environmental fluctuations.

Efforts have been made to reduce the genomes of living cells, but phage genome reduction remains challenging. It is of great interest to investigate whether genome reduction can make phages obtain new infectious properties. We developed a CRISPR/Cas9-based iterative phage genome reduction (CiPGr) approach and applied this to four distinct phages, thereby obtaining heterogeneous genome-reduced mutants. We isolated and sequenced 200 mutants with loss of up to 8-23% (3.3-35 kbp) of the original sequences. This allowed the identification of non-essential genes for phage propagation, although loss of these genes is mostly detrimental to phage fitness to various degrees. Notwithstanding this, mutants with higher infectious efficiency than their parental strains were characterized, indicating a trade-off between genome reduction and infectious fitness for phages. In conclusion, this study provides a foundation for future work to leverage the information generated by CiPGr in phage synthetic biology research.

Sulfuration of uridine 8, in bacterial and archaeal tRNAs, is catalyzed by enzymes formerly known as ThiI, but renamed here TtuI. Two different classes of TtuI proteins, which possess a PP-loop-containing pyrophosphatase domain that includes a conserved cysteine important for catalysis, have been identified. The first class, as exemplified by the prototypic Escherichia coli enzyme, possesses an additional C-terminal rhodanese domain harboring a second cysteine, which serves to form a catalytic persulfide. Among the second class of TtuI proteins that do not possess the rhodanese domain, some archaeal proteins display a conserved CXXC + C motif. We report here spectroscopic and enzymatic studies showing that TtuI from Methanococcus maripaludis and Pyrococcus furiosus can assemble a [4Fe-4S] cluster that is essential for tRNA sulfuration activity. Moreover, structural modeling studies, together with previously reported mutagenesis experiments of M. maripaludis TtuI, indicate that the [4Fe-4S] cluster is coordinated by the three cysteines of the CXXC + C motif. Altogether, our results raise a novel mechanism for U8-tRNA sulfuration, in which the cluster is proposed to catalyze the transfer of sulfur atoms to the activated tRNA substrate.

Drugs are designed to bind their target proteins in physiologically relevant tissues and organs to modulate biological functions and elicit desirable clinical outcomes. Information about target engagement at cellular and subcellular resolution is therefore critical for guiding compound optimization in drug discovery, and for probing resistance mechanisms to targeted therapies in clinical samples. We describe a target engagement-mediated amplification (TEMA) technology, where oligonucleotide-conjugated drugs are used to visualize and measure target engagement in situ, amplified via rolling-circle replication of circularized oligonucleotide probes. We illustrate the TEMA technique using dasatinib and gefitinib, two kinase inhibitors with distinct selectivity profiles. In vitro binding by the dasatinib probe to arrays of displayed proteins accurately reproduced known selectivity profiles, while their differential binding to fixed adherent cells agreed with expectations from expression profiles of the cells. We also introduce a proximity ligation variant of TEMA to selectively investigate binding to specific target proteins of interest. This form of the assay serves to improve resolution of binding to on- and off-target proteins. In conclusion, TEMA has the potential to aid in drug development and clinical routine by conferring valuable insights in drug-target interactions at spatial resolution in protein arrays, cells and in tissues.

Erratum in Nucleic Acids Res. 2022 Dec 09;:

The ability to create stimuli-responsive DNA nanostructures has played a prominent role in dynamic DNA nanotechnology. Primary among these is the process of toehold-based strand displacement, where a nucleic acid molecule can act as a trigger to cause conformational changes in custom-designed DNA nanostructures. Here, we add another layer of control to strand displacement reactions through a 'toehold clipping' process. By designing DNA complexes with a photocleavable linker-containing toehold or an RNA toehold, we show that we can use light (UV) or enzyme (ribonuclease) to eliminate the toehold, thus preventing strand displacement reactions. We use molecular dynamics simulations to analyze the structural effects of incorporating a photocleavable linker in DNA complexes. Beyond simple DNA duplexes, we also demonstrate the toehold clipping process in a model DNA nanostructure, by designing a toehold containing double-bundle DNA tetrahedron that disassembles when an invading strand is added, but stays intact after the toehold clipping process even in the presence of the invading strand. This work is an example of combining multiple physical or molecular stimuli to provide additional remote control over DNA nanostructure reconfiguration, advances that hold potential use in biosensing, drug delivery or molecular computation.

Nucleotide excision repair (NER) is critical for removing bulky DNA base lesions and avoiding diseases. NER couples lesion recognition by XPC to strand separation by XPB and XPD ATPases, followed by lesion excision by XPF and XPG nucleases. Here, we describe key regulatory mechanisms and roles of XPG for and beyond its cleavage activity. Strikingly, by combing single-molecule imaging and bulk cleavage assays, we found that XPG binding to the 7-subunit TFIIH core (coreTFIIH) stimulates coreTFIIH-dependent double-strand (ds)DNA unwinding 10-fold, and XPG-dependent DNA cleavage by up to 700-fold. Simultaneous monitoring of rates for coreTFIIH single-stranded (ss)DNA translocation and dsDNA unwinding showed XPG acts by switching ssDNA translocation to dsDNA unwinding as a likely committed step. Pertinent to the NER pathway regulation, XPG incision activity is suppressed during coreTFIIH translocation on DNA but is licensed when coreTFIIH stalls at the lesion or when ATP hydrolysis is blocked. Moreover, ≥15 nucleotides of 5'-ssDNA is a prerequisite for efficient translocation and incision. Our results unveil a paired coordination mechanism in which key lesion scanning and DNA incision steps are sequentially coordinated, and damaged patch removal is only licensed after generation of ≥15 nucleotides of 5'-ssDNA, ensuring the correct ssDNA bubble size before cleavage.

Spatial visualization of single-cell transcripts is limited by signal specificity and multiplexing. Here, we report hierarchical DNA branch assembly-encoded fluorescent nanoladders, which achieve denoised and highly multiplexed signal amplification for single-molecule transcript imaging. This method first offers independent RNA-primed rolling circle amplification without nonspecific amplification based on circular DNAzyme. It then executes programmable DNA branch assembly on these amplicons to encode virtual signals for visualizing numbers of targets by FISH. In theory, more virtual signals can be encoded via the increase of detection spectral channels and repeats of the same sequences on barcode. Our method almost eliminates nonspecific amplification in fixed cells (reducing nonspecific spots of single cells from 16 to nearly zero), and achieves simultaneous quantitation of nine transcripts by using only two detection spectral channels. We demonstrate accurate RNA profiling in different cancer cells, and reveal diverse localization patterns for spatial regulation of transcripts.

Aptamers are nucleic acid bioreceptors that have been used in various applications including medical diagnostics and as therapeutic agents. Identifying the most optimal aptamer for a particular application is very challenging. Here, we for the first time have developed a high-throughput method for accurately quantifying aptamer binding affinity, specificity, and cross-reactivity via the kinetics of aptamer digestion by exonucleases. We demonstrate the utility of this approach by isolating a set of new aptamers for fentanyl and its analogs, and then characterizing the binding properties of 655 aptamer-ligand pairs using our exonuclease digestion assay and validating the results with gold-standard methodologies. These data were used to select optimal aptamers for the development of new sensors that detect fentanyl and its analogs in different analytical contexts. Our approach dramatically accelerates the aptamer characterization process and streamlines sensor development, and if coupled with robotics, could enable high-throughput quantitative analysis of thousands of aptamer-ligand pairs.

The vast majority of disease-associated single nucleotide polymorphisms (SNP) identified from genome-wide association studies (GWAS) are localized in non-coding regions. A significant fraction of these variants impact transcription factors binding to enhancer elements and alter gene expression. To functionally interrogate the activity of such variants we developed snpSTARRseq, a high-throughput experimental method that can interrogate the functional impact of hundreds to thousands of non-coding variants on enhancer activity. snpSTARRseq dramatically improves signal-to-noise by utilizing a novel sequencing and bioinformatic approach that increases both insert size and the number of variants tested per loci. Using this strategy, we interrogated known prostate cancer (PCa) risk-associated loci and demonstrated that 35% of them harbor SNPs that significantly altered enhancer activity. Combining these results with chromosomal looping data we could identify interacting genes and provide a mechanism of action for 20 PCa GWAS risk regions. When benchmarked to orthogonal methods, snpSTARRseq showed a strong correlation with in vivo experimental allelic-imbalance studies whereas there was no correlation with predictive in silico approaches. Overall, snpSTARRseq provides an integrated experimental and computational framework to functionally test non-coding genetic variants.

Proper RNA localisation is essential for physiological gene expression. Various kinds of genome-wide approaches permit to comprehensively profile subcellular transcriptomes. Among them, cell fractionation methods, that couple RNase treatment of isolated organelles to the sequencing of protected transcripts, remain most widely used, mainly because they do not require genetic modification of the studied system and can be easily implemented in any cells or tissues, including in non-model species. However, they suffer from numerous false-positives since incompletely digested contaminant RNAs can still be captured and erroneously identified as resident transcripts. Here we introduce Controlled Level of Contamination coupled to deep sequencing (CoLoC-seq) as a new subcellular transcriptomics approach that efficiently bypasses this caveat. CoLoC-seq leverages classical enzymatic kinetics and tracks the depletion dynamics of transcripts in a gradient of an exogenously added RNase, with or without organellar membranes. By means of straightforward mathematical modelling, CoLoC-seq infers the localisation topology of RNAs and robustly distinguishes between genuinely resident, luminal transcripts and merely abundant surface-attached contaminants. Our generic approach performed well on human mitochondria and is in principle applicable to other membrane-bounded organelles, including plastids, compartments of the vacuolar system, extracellular vesicles, and viral particles.

Transfer RNA (tRNA) utilizes multiple properties of abundance, modification, and aminoacylation in translational regulation. These properties were typically studied one-by-one; however, recent advance in high throughput tRNA sequencing enables their simultaneous assessment in the same sequencing data. How these properties are coordinated at the transcriptome level is an open question. Here, we develop a single-read tRNA analysis pipeline that takes advantage of the pseudo single-molecule nature of tRNA sequencing in NGS libraries. tRNAs are short enough that a single NGS read can represent one tRNA molecule, and can simultaneously report on the status of multiple modifications, aminoacylation, and fragmentation of each molecule. We find correlations among modification-modification, modification-aminoacylation and modification-fragmentation. We identify interdependencies among one of the most common tRNA modifications, m1A58, as coordinators of tissue-specific gene expression. Our method, SingLe-read Analysis of Crosstalks (SLAC), reveals tRNAome-wide networks of modifications, aminoacylation, and fragmentation. We observe changes of these networks under different stresses, and assign a function for tRNA modification in translational regulation and fragment biogenesis. SLAC leverages the richness of the tRNA-seq data and provides new insights on the coordination of tRNA properties.

The increasing quantity of multi-omic data, such as methylomic and transcriptomic profiles collected on the same specimen or even on the same cell, provides a unique opportunity to explore the complex interactions that define cell phenotype and govern cellular responses to perturbations. We propose a network approach based on Gaussian Graphical Models (GGMs) that facilitates the joint analysis of paired omics data. This method, called DRAGON (Determining Regulatory Associations using Graphical models on multi-Omic Networks), calibrates its parameters to achieve an optimal trade-off between the network's complexity and estimation accuracy, while explicitly accounting for the characteristics of each of the assessed omics 'layers.' In simulation studies, we show that DRAGON adapts to edge density and feature size differences between omics layers, improving model inference and edge recovery compared to state-of-the-art methods. We further demonstrate in an analysis of joint transcriptome - methylome data from TCGA breast cancer specimens that DRAGON can identify key molecular mechanisms such as gene regulation via promoter methylation. In particular, we identify Transcription Factor AP-2 Beta (TFAP2B) as a potential multi-omic biomarker for basal-type breast cancer. DRAGON is available as open-source code in Python through the Network Zoo package (netZooPy v0.8; netzoo.github.io).

CRISPR/Cas-based approaches have largely replaced conventional gene targeting strategies. However, homology-directed repair (HDR) in the mouse genome is not very efficient, and precisely inserting longer sequences using HDR remains challenging given that donor constructs preferentially integrate as concatemers. Here, we showed that injecting 5' biotinylated donor DNA into mouse embryos at the two-cell stage led to efficient single-copy HDR (scHDR) allele generation. Our dedicated genotyping strategy showed that these alleles occurred with frequencies of 19%, 20%, and 26% at three independent gene loci, indicating that scHDR was dramatically increased by 5' biotinylation. Thus, we suggest that the combination of a 5' biotinylated donor and diligent analysis of concatemer integration are prerequisites for efficiently and reliably generating conditional alleles or other large fragment knock-ins in the mouse genome.