The development of the CRISPR-Cas9 technology has provided a simple yet powerful system for genome editing. Current gRNA design tools serve as an important platform for the efficient application of the CRISPR systems. However, most of the existing tools are black-box models that suffer from limitations, such as variable performance and unclear mechanism of decision making. Here, we introduce CRISPRedict, an interpretable gRNA efficiency prediction model for CRISPR-Cas9 gene editing. Its strength lies in the fact that it can accurately predict efficient guide RNAs-with equivalent performance to state-of-the-art tools-while being a simple linear model. Implemented as a user-friendly web server, CRISPRedict offers (i) quick and accurate predictions across various experimental conditions (e.g. U6/T7 transcription); (ii) regression and classification models for scoring gRNAs and (iii) multiple visualizations to explain the obtained results. Given its performance, interpretability, and versatility, we expect that it will assist researchers in the gRNA design process and facilitate genome editing research. CRISPRedict is available for use at http://www.crispredict.org/.

Residue coevolution within and between proteins is used as a marker of physical interaction and/or residue functional cooperation. Pairs or groups of coevolving residues are extracted from multiple sequence alignments based on a variety of computational approaches. However, coevolution signals emerging in subsets of sequences might be lost if the full alignment is considered. iBIS2Analyzer is a web server dedicated to a phylogeny-driven coevolution analysis of protein families with different evolutionary pressure. It is based on the iterative version, iBIS2, of the coevolution analysis method BIS, Blocks in Sequences. iBIS2 is designed to iteratively select and analyse subtrees in phylogenetic trees, possibly large and comprising thousands of sequences. With iBIS2Analyzer, openly accessible at http://ibis2analyzer.lcqb.upmc.fr/, the user visualizes, compares and inspects clusters of coevolving residues by mapping them onto sequences, alignments or structures of choice, greatly simplifying downstream analysis steps. A rich and interactive graphic interface facilitates the biological interpretation of the results.

Recent high-throughput omics techniques have produced a large amount of biological data. Visualization of big omics data is essential to answer a wide range of biological problems. As a concise but comprehensive strategy, a heatmap can analyze and visualize high-dimensional and heterogeneous biomolecular expression data in an attractive artwork. In 2014, we developed a stand-alone software package, Heat map Illustrator (HemI 1.0), which implemented three clustering methods and seven distance metrics for heatmap illustration. Here, we significantly improved 1.0 and released the online service of HemI 2.0, in which 7 clustering methods and 22 types of distance metrics were implemented. In HemI 2.0, the clustering results and publication-quality heatmaps can be exported directly. For an in-depth analysis of the data, we further added an option of enrichment analysis for 12 model organisms, with 15 types of functional annotations. The enrichment results can be visualized in five idioms, including bubble chart, bar graph, coxcomb chart, pie chart and word cloud. We anticipate that HemI 2.0 can be a helpful web server for visualization of biomolecular expression data, as well as the additional enrichment analysis. HemI 2.0 is freely available for all users at: https://hemi.biocuckoo.org/.

Circadian rhythms are a foundational aspect of biology. These rhythms are found at the molecular level in every cell of every living organism and they play a fundamental role in homeostasis and a variety of physiological processes. As a result, biomedical research of circadian rhythms continues to expand at a rapid pace. To support this research, CircadiOmics (http://circadiomics.igb.uci.edu/) is the largest annotated repository and analytic web server for high-throughput omic (e.g. transcriptomic, metabolomic, proteomic) circadian time series experimental data. CircadiOmics contains over 290 experiments and over 100 million individual measurements, across >20 unique tissues/organs, and 11 different species. Users are able to visualize and mine these datasets by deriving and comparing periodicity statistics for oscillating molecular species including: period, amplitude, phase, P-value and q-value. These statistics are obtained from BIO_CYCLE and JTK_CYCLE and are intuitively aggregated and displayed for comparison. CircadiOmics is the most up-to-date and cutting-edge web portal for searching and analyzing circadian omic data and is used by researchers around the world.

RNA structures play critical roles in regulating gene expression across all domains of life and viruses. Chemical probing methods coupled with massively parallel sequencing have revolutionized the RNA structure field by enabling the assessment of many structures in their native, physiological context. Previously, we developed Dimethyl-Sulfate-based Mutational Profiling and Sequencing (DMS-MaPseq), which uses DMS to label the Watson-Crick face of open and accessible adenine and cytosine bases in the RNA. We used this approach to determine the genome-wide structures of HIV-1 and SARS-CoV-2 in infected cells, which permitted uncovering new biology and identifying therapeutic targets. Due to the simplicity and ease of the experimental procedure, DMS-MaPseq has been adopted by labs worldwide. However, bioinformatic analysis remains a substantial hurdle for labs that often lack the necessary infrastructure and computational expertise. Here we present a modern web-based interface that automates the analysis of chemical probing profiles from raw sequencing files (http://rnadreem.org). The availability of a simple web-based platform for DMS-MaPseq analysis will dramatically expand studies of RNA structure and aid in the design of structure-based therapeutics.

In eukaryotes and in archaea late steps of translation initiation involve the two initiation factors e/aIF5B and e/aIF1A. In eukaryotes, the role of eIF5B in ribosomal subunit joining is established and structural data showing eIF5B bound to the full ribosome were obtained. To achieve its function, eIF5B collaborates with eIF1A. However, structural data illustrating how these two factors interact on the small ribosomal subunit have long been awaited. The role of the archaeal counterparts, aIF5B and aIF1A, remains to be extensively addressed. Here, we study the late steps of Pyrococcus abyssi translation initiation. Using in vitro reconstituted initiation complexes and light scattering, we show that aIF5B bound to GTP accelerates subunit joining without the need for GTP hydrolysis. We report the crystallographic structures of aIF5B bound to GDP and GTP and analyze domain movements associated to these two nucleotide states. Finally, we present the cryo-EM structure of an initiation complex containing 30S bound to mRNA, Met-tRNAiMet, aIF5B and aIF1A at 2.7 Å resolution. Structural data shows how archaeal 5B and 1A factors cooperate to induce a conformation of the initiator tRNA favorable to subunit joining. Archaeal and eukaryotic features of late steps of translation initiation are discussed.

Schlafen-5 (SLFN5) is an interferon-induced protein of the Schlafen family, which are involved in immune responses and oncogenesis. To date, little is known regarding its anti-HIV-1 function. Here, the authors report that overexpression of SLFN5 inhibits HIV-1 replication and reduces viral mRNA levels, whereas depletion of endogenous SLFN5 promotes HIV-1 replication. Moreover, they show that SLFN5 markedly decreases the transcriptional activity of HIV-1 long terminal repeat (LTR) via binding to two sequences in the U5-R region, which consequently represses the recruitment of RNA polymerase II to the transcription initiation site. Mutagenesis studies show the importance of nuclear localization and the N-terminal 1-570 amino acids fragment in the inhibition of HIV-1. Further mechanistic studies demonstrate that SLFN5 interacts with components of the PRC2 complex, G9a and Histone H3, thereby promoting H3K27me2 and H3K27me3 modification leading to silencing HIV-1 transcription. In concert with this, they find that SLFN5 blocks the activation of latent HIV-1. Altogether, their findings demonstrate that SLFN5 is a transcriptional repressor of HIV-1 through epigenetic modulation and a potential determinant of HIV-1 latency.

We have investigated the function of human topoisomerase 1 (TOP1) in regulation of G-quadruplex (G4) formation in the Pu27 region of the MYC P1 promoter. Pu27 is among the best characterized G4 forming sequences in the human genome and it is well known that promoter activity is inhibited upon G4 formation in this region. We found that TOP1 downregulation stimulated transcription from a promoter with wildtype Pu27 but not if the G4 motif in Pu27 was interrupted by mutation(s). The effect was not specific to the MYC promoter and similar results were obtained for the G4 forming promoter element WT21. The other major DNA topoisomerases with relaxation activity, topoisomerases 2α and β, on the other hand, did not affect G4 dependent promoter activity. The cellular studies were supported by in vitro investigations demonstrating a high affinity of TOP1 for wildtype Pu27 but not for mutant sequences unable to form G4. Moreover, TOP1 was able to induce G4 formation in Pu27 inserted in double stranded plasmid DNA in vitro. This is the first time TOP1 has been demonstrated capable of inducing G4 formation in double stranded DNA and of influencing G4 formation in cells.

Human Positive Coactivator 4 (PC4) is a multifaceted chromatin protein involved in diverse cellular processes including genome organization, transcription regulation, replication, DNA repair and autophagy. PC4 exists as a phospho-protein in cells which impinges on its acetylation by p300 and thereby affects its transcriptional co-activator functions via double-stranded DNA binding. Despite the inhibitory effects, the abundance of phosphorylated PC4 in cells intrigued us to investigate its role in chromatin functions in a basal state of the cell. We found that casein kinase-II (CKII)-mediated phosphorylation of PC4 is critical for its interaction with linker histone H1. By employing analytical ultracentrifugation and electron microscopy imaging of in vitro reconstituted nucleosomal array, we observed that phospho-mimic (PM) PC4 displays a superior chromatin condensation potential in conjunction with linker histone H1. ATAC-sequencing further unveiled the role of PC4 phosphorylation to be critical in inducing chromatin compaction of a wide array of coding and non-coding genes in vivo. Concordantly, phospho-PC4 mediated changes in chromatin accessibility led to gene repression and affected global histone modifications. We propose that the abundance of PC4 in its phosphorylated state contributes to genome compaction contrary to its co-activator function in driving several cellular processes like gene transcription and autophagy.

Simultaneous targeting multiple genes is a big advantage of CRISPR (clustered regularly interspaced short palindromic repeats) genome editing but challenging to achieve in CRISPR screening. The crosstalk among genes or gene products is a common and fundamental mechanism to ensure cellular stability and functional diversity. However, the screening approach to map high-order gene combinations to the interesting phenotype is still lacking. Here, we developed a universal in-library ligation strategy and applied it to generate multiplexed CRISPR library, which could perturb four pre-designed targets in a cell. We conducted in vivo CRISPR screening for potential guide RNA (gRNA) combinations inducing anti-tumor immune responses. Simultaneously disturbing a combination of three checkpoints in CD8+ T cells was demonstrated to be more effective than disturbing Pdcd1 only for T cell activation in the tumor environment. This study developed a novel in-library ligation strategy to facilitate the multiplexed CRISPR screening, which could extend our ability to explore the combinatorial outcomes from coordinated gene behaviors.

Dynamic regulation is an effective strategy for control of gene expression in microbial cell factories. In some pathway contexts, several metabolic modules must be controlled in a time dependent or ordered manner to maximize production, while the creation of genetic circuits with ordered regulation capacity still remains a great challenge. In this work, we develop a pathway independent and programmable system that enables multi-modular ordered control of metabolism in Bacillus subtilis. First, a series of thermosensors were created and engineered to expand their thresholds. Then we designed single-input-multi-output circuits for ordered control based on the use of thermosensors with different transition points. Meanwhile, a repression circuit was constructed by combining CRISPRi-based NOT gates. As a proof-of-concept, these genetic circuits were applied for multi-modular ordered control of 2'-fucosyllactose (2'-FL) biosynthesis, resulting in a production of 1839.7 mg/l in shake flask, which is 5.16-times that of the parental strain. In a 5-l bioreactor, the 2'-FL titer reached 28.2 g/l with down-regulation of autolysis. Taken together, this work provides programmable and versatile thermosensitive genetic toolkits for dynamic regulation in B. subtilis and a multi-modular ordered control framework that can be used to improve metabolic modules in other chassis cells and for other compounds.

The integrity and proper expression of genomes are safeguarded by DNA and RNA surveillance pathways. While many RNA surveillance factors have additional functions in the nucleus, little is known about the incidence and physiological impact of converging RNA and DNA signals. Here, using genetic screens and genome-wide analyses, we identified unforeseen SMG-1-dependent crosstalk between RNA surveillance and DNA repair in living animals. Defects in RNA processing, due to viable THO complex or PNN-1 mutations, induce a shift in DNA repair in dividing and non-dividing tissues. Loss of SMG-1, an ATM/ATR-like kinase central to RNA surveillance by nonsense-mediated decay (NMD), restores DNA repair and radio-resistance in THO-deficient animals. Mechanistically, we find SMG-1 and its downstream target SMG-2/UPF1, but not NMD per se, to suppress DNA repair by non-homologous end-joining in favour of single strand annealing. We postulate that moonlighting proteins create short-circuits in vivo, allowing aberrant RNA to redirect DNA repair.

The X family polymerases (PolXs) are specialized DNA polymerases that are found in all domains of life. While the main representatives of eukaryotic PolXs, which have dedicated functions in DNA repair, were studied in much detail, the functions and diversity of prokaryotic PolXs have remained largely unexplored. Here, by combining a comprehensive bioinformatic analysis of prokaryotic PolXs and biochemical experiments involving selected recombinant enzymes, we reveal a previously unrecognized group of PolXs that seem to be lacking DNA polymerase activity. The noncanonical PolXs contain substitutions of the key catalytic residues and deletions in their polymerase and dNTP binding sites in the palm and fingers domains, but contain functional nuclease domains, similar to canonical PolXs. We demonstrate that representative noncanonical PolXs from the Deinococcus genus are indeed inactive as DNA polymerases but are highly efficient as 3'-5' exonucleases. We show that both canonical and noncanonical PolXs are often encoded together with the components of the non-homologous end joining pathway and may therefore participate in double-strand break repair, suggesting an evolutionary conservation of this PolX function. This is a remarkable example of polymerases that have lost their main polymerase activity, but retain accessory functions in DNA processing and repair.

The chromosome dimer resolution machinery of bacteria is generally composed of two tyrosine recombinases, XerC and XerD. They resolve chromosome dimers by adding a crossover between sister copies of a specific site, dif. The reaction depends on a cell division protein, FtsK, which activates XerD by protein-protein interactions. The toxin-linked cryptic satellite phage (TLCΦ) of Vibrio cholerae, which participates in the emergence of cholera epidemic strains, carries a dif-like attachment site (attP). TLCΦ exploits the Xer machinery to integrate into the dif site of its host chromosomes. The TLCΦ integration reaction escapes the control of FtsK because TLCΦ encodes for its own XerD-activation factor, XafT. Additionally, TLCΦ attP is a poor substrate for XerD binding, in apparent contradiction with the high integration efficiency of the phage. Here, we present a sequencing-based methodology to analyse the integration and excision efficiency of thousands of synthetic mini-TLCΦ plasmids with differing attP sites in vivo. This methodology is applicable to the fine-grained analyses of DNA transactions on a wider scale. In addition, we compared the efficiency with which XafT and the XerD-activation domain of FtsK drive recombination reactions in vitro. Our results suggest that XafT not only activates XerD-catalysis but also helps form and/or stabilize synaptic complexes between imperfect Xer recombination sites.

NAT10 is an essential enzyme that catalyzes N4-acetylcytidine (ac4C) in eukaryotic transfer RNA and 18S ribosomal RNA. Recent studies suggested that rRNA acetylation is dependent on SNORD13, a box C/D small nucleolar RNA predicted to base-pair with 18S rRNA via two antisense elements. However, the selectivity of SNORD13-dependent cytidine acetylation and its relationship to NAT10's essential function remain to be defined. Here, we demonstrate that SNORD13 is required for acetylation of a single cytidine of human and zebrafish 18S rRNA. In-depth characterization revealed that SNORD13-dependent ac4C is dispensable for human cell growth, ribosome biogenesis, translation and development. This loss of function analysis inspired a cross-evolutionary survey of the eukaryotic rRNA acetylation 'machinery' that led to the characterization of many novel metazoan SNORD13 genes. This includes an atypical SNORD13-like RNA in Drosophila melanogaster which guides ac4C to 18S rRNA helix 45 despite lacking one of the two rRNA antisense elements. Finally, we discover that Caenorhabditis elegans 18S rRNA is not acetylated despite the presence of an essential NAT10 homolog. Our findings shed light on the molecular mechanisms underlying SNORD13-mediated rRNA acetylation across eukaryotic evolution and raise new questions regarding the biological and evolutionary relevance of this highly conserved rRNA modification.

CRISPR/Cas13 effectors have garnered increasing attention as easily customizable tools for detecting and depleting RNAs of interest. Near perfect complementarity between a target RNA and the Cas13-associated guide RNA is required for activation of Cas13 ribonuclease activity. Nonetheless, the specificity of Cas13 effectors in eukaryotic cells has been debated as the Cas13 nuclease domains can be exposed on the enzyme surface, providing the potential for promiscuous cleavage of nearby RNAs (so-called collateral damage). Here, using co-transfection assays in Drosophila and human cells, we found that the off-target effects of RxCas13d, a commonly used Cas13 effector, can be as strong as the level of on-target RNA knockdown. The extent of off-target effects is positively correlated with target RNA expression levels, and collateral damage can be observed even after reducing RxCas13d/guide RNA levels. The PspCas13b effector showed improved specificity and, unlike RxCas13d, can be used to deplete a Drosophila circular RNA without affecting the expression of the associated linear RNA. PspCas13b nonetheless still can have off-target effects and we notably found that the extent of off-target effects for Cas13 effectors differs depending on the cell type and target RNA examined. In total, these results highlight the need for caution when designing and interpreting Cas13-based knockdown experiments.

CRISPR/Cas12a is a single effector nuclease that, like CRISPR/Cas9, has been harnessed for genome editing based on its ability to generate targeted DNA double strand breaks (DSBs). Unlike the blunt-ended DSB generated by Cas9, Cas12a generates sticky-ended DSB that could potentially aid precise genome editing, but this unique feature has thus far been underutilized. In the current study, we found that a short double-stranded DNA (dsDNA) repair template containing a sticky end that matched one of the Cas12a-generated DSB ends and a homologous arm sharing homology with the genomic region adjacent to the other end of the DSB enabled precise repair of the DSB and introduced a desired nucleotide substitution. We termed this strategy 'Ligation-Assisted Homologous Recombination' (LAHR). Compared to the single-stranded oligo deoxyribonucleotide (ssODN)-mediated homology directed repair (HDR), LAHR yields relatively high editing efficiency as demonstrated for both a reporter gene and endogenous genes. We found that both HDR and microhomology-mediated end joining (MMEJ) mechanisms are involved in the LAHR process. Our LAHR genome editing strategy, extends the repertoire of genome editing technologies and provides a broader understanding of the type and role of DNA repair mechanisms involved in genome editing.

Erratum in Nucleic Acids Res. 2022 Sep 9;50(16):9597.

Translational readthrough (TR) occurs when the ribosome decodes a stop codon as a sense codon, resulting in two protein isoforms synthesized from the same mRNA. TR has been identified in several eukaryotic organisms; however, its biological significance and mechanism remain unclear. Here, we quantify TR of several candidate genes in Drosophila melanogaster and characterize the regulation of TR in the large Maf transcription factor Traffic jam (Tj). Using CRISPR/Cas9-generated mutant flies, we show that the TR-generated Tj isoform is expressed in a subset of neural cells of the central nervous system and is excluded from the somatic cells of gonads. Control of TR in Tj is critical for preservation of neuronal integrity and maintenance of reproductive health. The tissue-specific distribution of a release factor splice variant, eRF1H, plays a critical role in modulating differential TR of leaky stop codon contexts. Fine-tuning of gene regulatory functions of transcription factors by TR provides a potential mechanism for cell-specific regulation of gene expression.

Reprogramming of transcription is critical for the survival under cellular stress. Heat shock has provided an excellent model to investigate nascent transcription in stressed cells, but the molecular mechanisms orchestrating RNA synthesis during other types of stress are unknown. We utilized PRO-seq and ChIP-seq to study how Heat Shock Factors, HSF1 and HSF2, coordinate transcription at genes and enhancers upon oxidative stress and heat shock. We show that pause-release of RNA polymerase II (Pol II) is a universal mechanism regulating gene transcription in stressed cells, while enhancers are activated at the level of Pol II recruitment. Moreover, besides functioning as conventional promoter-binding transcription factors, HSF1 and HSF2 bind to stress-induced enhancers to trigger Pol II pause-release from poised gene promoters. Importantly, HSFs act at distinct genes and enhancers in a stress type-specific manner. HSF1 binds to many chaperone genes upon oxidative and heat stress but activates them only in heat-shocked cells. Under oxidative stress, HSF1 localizes to a unique set of promoters and enhancers to trans-activate oxidative stress-specific genes. Taken together, we show that HSFs function as multi-stress-responsive factors that activate distinct genes and enhancers when encountering changes in temperature and redox state.

Genome-scale metabolic models have been recognised as useful tools for better understanding living organisms' metabolism. merlin (https://www.merlin-sysbio.org/) is an open-source and user-friendly resource that hastens the models' reconstruction process, conjugating manual and automatic procedures, while leveraging the user's expertise with a curation-oriented graphical interface. An updated and redesigned version of merlin is herein presented. Since 2015, several features have been implemented in merlin, along with deep changes in the software architecture, operational flow, and graphical interface. The current version (4.0) includes the implementation of novel algorithms and third-party tools for genome functional annotation, draft assembly, model refinement, and curation. Such updates increased the user base, resulting in multiple published works, including genome metabolic (re-)annotations and model reconstructions of multiple (lower and higher) eukaryotes and prokaryotes. merlin version 4.0 is the only tool able to perform template based and de novo draft reconstructions, while achieving competitive performance compared to state-of-the art tools both for well and less-studied organisms.

Nucleobase deamination, such as A-to-I editing, represents an important posttranscriptional modification of RNA. When deamination affects guanosines, a xanthosine (X) containing RNA is generated. However, the biological significance and chemical consequences on RNA are poorly understood. We present a comprehensive study on the preparation and biophysical properties of X-modified RNA. Thermodynamic analyses revealed that base pairing strength is reduced to a level similar to that observed for a G•U replacement. Applying NMR spectroscopy and X-ray crystallography, we demonstrate that X can form distinct wobble geometries with uridine depending on the sequence context. In contrast, X pairing with cytidine occurs either through wobble geometry involving protonated C or in Watson-Crick-like arrangement. This indicates that the different pairing modes are of comparable stability separated by low energetic barriers for switching. Furthermore, we demonstrate that the flexible pairing properties directly affect the recognition of X-modified RNA by reverse transcription enzymes. Primer extension assays and PCR-based sequencing analysis reveal that X is preferentially read as G or A and that the ratio depends on the type of reverse transcriptase. Taken together, our results elucidate important properties of X-modified RNA paving the way for future studies on its biological significance.

Heterogenous nuclear ribonucleoproteins (hnRNPs) are abundant proteins implicated in various steps of RNA processing that assemble on nuclear RNA into larger complexes termed 40S hnRNP particles. Despite their initial discovery 55 years ago, our understanding of these intriguing macromolecular assemblies remains limited. Here, we report the biochemical purification of native 40S hnRNP particles and the determination of their complete protein composition by label-free quantitative mass spectrometry, identifying A-group and C-group hnRNPs as the major protein constituents. Isolated 40S hnRNP particles dissociate upon RNA digestion and can be reconstituted in vitro on defined RNAs in the presence of the individual protein components, demonstrating a scaffolding role for RNA in nucleating particle formation. Finally, we revealed their nanometer scale, condensate-like nature, promoted by intrinsically disordered regions of A-group hnRNPs. Collectively, we identify nuclear 40S hnRNP particles as novel dynamic biomolecular condensates.

CRISPR-Cas12a (Cpf1) is a bacterial RNA-guided nuclease that cuts double-stranded DNA (dsDNA) at sites specified by a CRISPR RNA (crRNA) guide. Additional activities have been ascribed to this enzyme in vitro: site-specific (cis) single-stranded DNA (ssDNA) cleavage and indiscriminate (trans) degradation of ssDNA, RNA, and dsDNA after activation by a complementary target. The ability of Cas12a to cleave nucleic acids indiscriminately has been harnessed for many applications, including diagnostics, but it remains unknown if it contributes to bacterial immunity. Here, we provide evidence that cleavage of ssDNA in cis or in trans by Cas12a is insufficient to impact immunity. Using LbCas12a expressed in either Pseudomonas aeruginosa or Escherichia coli, we observed that cleavage of dsDNA targets did not elicit cell death or dormancy, suggesting insignificant levels of collateral damage against host RNA or DNA. Canonical immunity against invasive dsDNA also had no impact on the replicative fitness of co-infecting dsDNA phage, ssDNA phage or plasmid in trans. Lastly, crRNAs complementary to invasive ssDNA did not provide protection, suggesting that ssDNA cleavage does not occur in vivo or is insignificant. Overall, these results suggest that CRISPR-Cas12a immunity predominantly occurs via canonical targeting of dsDNA, and that the other activities do not significantly impact infection outcomes.

In every domain of life, NusG-like proteins bind to the elongating RNA polymerase (RNAP) to support processive RNA synthesis and to couple transcription to ongoing cellular processes. Structures of factor-bound transcription elongation complexes (TECs) reveal similar contacts to RNAP, consistent with a shared mechanism of action. However, NusG homologs differ in their regulatory roles, modes of recruitment, and effects on RNA synthesis. Some of these differences could be due to conformational changes in RNAP and NusG-like proteins, which cannot be captured in static structures. Here, we employed hydrogen-deuterium exchange mass spectrometry to investigate changes in local and non-local structural dynamics of Escherichia coli NusG and its paralog RfaH, which have opposite effects on expression of xenogenes, upon binding to TEC. We found that NusG and RfaH regions that bind RNAP became solvent-protected in factor-bound TECs, whereas RNAP regions that interact with both factors showed opposite deuterium uptake changes when bound to NusG or RfaH. Additional changes far from the factor-binding site were observed only with RfaH. Our results provide insights into differences in structural dynamics exerted by NusG and RfaH during binding to TEC, which may explain their different functional outcomes and allosteric regulation of transcriptional pausing by RfaH.

DNA transcription is regulated by a range of diverse mechanisms and primarily by transcription factors that recruit the RNA polymerase complex to the promoter region on the DNA. Protein binding to DNA at nearby or distant sites can synergistically affect this process in a variety of ways, but mainly through direct interactions between DNA-binding proteins. Here we show that a Transcription Activator-Like Effector (TALE), which lacks an activation domain, can enhance transcription in mammalian cells when it binds in the vicinity of and without direct interaction with several different dimeric or monomeric transcription factors. This effect was observed for several TALEs regardless of the recognition sequences and their DNA-bound orientation. TALEs can exert an effect over the distance of tens of nucleotides and it also potentiated KRAB-mediated repression. The augmentation of transcriptional regulation of another transcription factor is characteristic of TALEs, as it was not observed for dCas9/gRNA, zinc finger, or Gal4 DNA-binding domains. We propose that this mechanism involves an allosteric effect exerted on DNA structure or dynamics. This mechanism could be used to modulate transcription but may also play a role in the natural context of TALEs.

Poly(ADP-ribose) polymerase-1 (PARP-1) is a DNA damage sensor and contributes to both DNA repair and cell death processes. However, how PARP-1 signaling is regulated to switch its function from DNA repair to cell death remains largely unknown. Here, we found that PARP-1 plays a central role in alkylating agent-induced PARthanatic cancer cell death. Lysine demethylase 6B (KDM6B) was identified as a key regulator of PARthanatos. Loss of KDM6B protein or its demethylase activity conferred cancer cell resistance to PARthanatic cell death in response to alkylating agents. Mechanistically, KDM6B knockout suppressed methylation at the promoter of O6-methylguanine-DNA methyltransferase (MGMT) to enhance MGMT expression and its direct DNA repair function, thereby inhibiting DNA damage-evoked PARP-1 hyperactivation and subsequent cell death. Moreover, KDM6B knockout triggered sustained Chk1 phosphorylation and activated a second XRCC1-dependent repair machinery to fix DNA damage evading from MGMT repair. Inhibition of MGMT or checkpoint response re-sensitized KDM6B deficient cells to PARthanatos induced by alkylating agents. These findings provide new molecular insights into epigenetic regulation of PARP-1 signaling mediating DNA repair or cell death and identify KDM6B as a biomarker for prediction of cancer cell vulnerability to alkylating agent treatment.

Reverse transcriptases (RTs) are enzymes capable of synthesizing DNA using RNA as a template. Within the last few years, a burst of research has led to the discovery of novel prokaryotic RTs with diverse antiviral properties, such as DRTs (Defense-associated RTs), which belong to the so-called group of unknown RTs (UG) and are closely related to the Abortive Infection system (Abi) RTs. In this work, we performed a systematic analysis of UG and Abi RTs, increasing the number of UG/Abi members up to 42 highly diverse groups, most of which are predicted to be functionally associated with other gene(s) or domain(s). Based on this information, we classified these systems into three major classes. In addition, we reveal that most of these groups are associated with defense functions and/or mobile genetic elements, and demonstrate the antiphage role of four novel groups. Besides, we highlight the presence of one of these systems in novel families of human gut viruses infecting members of the Bacteroidetes and Firmicutes phyla. This work lays the foundation for a comprehensive and unified understanding of these highly diverse RTs with enormous biotechnological potential.

During translation, nascent polypeptide chains travel from the peptidyl transferase center through the nascent polypeptide exit tunnel (NPET) to emerge from 60S subunits. The NPET includes portions of five of the six 25S/5.8S rRNA domains and ribosomal proteins uL4, uL22, and eL39. Internal loops of uL4 and uL22 form the constriction sites of the NPET and are important for both assembly and function of ribosomes. Here, we investigated the roles of eL39 in tunnel construction, 60S biogenesis, and protein synthesis. We show that eL39 is important for proper protein folding during translation. Consistent with a delay in processing of 27S and 7S pre-rRNAs, eL39 functions in pre-60S assembly during middle nucleolar stages. Our biochemical assays suggest the presence of eL39 in particles at these stages, although it is not visualized in them by cryo-electron microscopy. This indicates that eL39 takes part in assembly even when it is not fully accommodated into the body of pre-60S particles. eL39 is also important for later steps of assembly, rotation of the 5S ribonucleoprotein complex, likely through long range rRNA interactions. Finally, our data strongly suggest the presence of alternative pathways of ribosome assembly, previously observed in the biogenesis of bacterial ribosomal subunits.

Most genomes harbor a large number of transposons, and they play an important role in evolution and gene regulation. They are also of interest to clinicians as they are involved in several diseases, including cancer and neurodegeneration. Although several methods for transposon identification are available, they are often highly specialised towards specific tasks or classes of transposons, and they lack common standards such as a unified taxonomy scheme and output file format. We present TransposonUltimate, a powerful bundle of three modules for transposon classification, annotation, and detection of transposition events. TransposonUltimate comes as a Conda package under the GPL-3.0 licence, is well documented and it is easy to install through https://github.com/DerKevinRiehl/TransposonUltimate. We benchmark the classification module on the large TransposonDB covering 891,051 sequences to demonstrate that it outperforms the currently best existing solutions. The annotation and detection modules combine sixteen existing softwares, and we illustrate its use by annotating Caenorhabditis elegans, Rhizophagus irregularis and Oryza sativa subs. japonica genomes. Finally, we use the detection module to discover 29 554 transposition events in the genomes of 20 wild type strains of C. elegans. Databases, assemblies, annotations and further findings can be downloaded from (https://doi.org/10.5281/zenodo.5518085).

In eukaryotes, fine-scale maps of meiotic recombination events have greatly advanced our understanding of the factors that affect genomic variation patterns and evolution of traits. However, in bacteria that lack natural systems for sexual reproduction, unbiased characterization of recombination landscapes has remained challenging due to variable rates of genetic exchange and influence of natural selection. Here, to overcome these limitations and to gain a genome-wide view on recombination, we crossed Bacillus strains with different genetic distances using protoplast fusion. The offspring displayed complex inheritance patterns with one of the parents consistently contributing the major part of the chromosome backbone and multiple unselected fragments originating from the second parent. Our results demonstrate that this bias was in part due to the action of restriction-modification systems, whereas genome features like GC content and local nucleotide identity did not affect distribution of recombination events around the chromosome. Furthermore, we found that recombination occurred uniformly across the genome without concentration into hotspots. Notably, our results show that species-level genetic distance did not affect genome-wide recombination. This study provides a new insight into the dynamics of recombination in bacteria and a platform for studying recombination patterns in diverse bacterial species.

Prime editors consisting of Cas9-nickase and reverse transcriptase enable targeted precise editing of small DNA pieces, including all 12 kinds of base substitutions, insertions and deletions, while without requiring double-strand breaks or donor templates. Current optimized prime editing strategy (PE3) uses two guide RNAs to guide the performance of prime editor. One guide RNA carrying both spacer and templating sequences (pegRNA) guides prime editor to produce ssDNA break and subsequent extension, and the other one produces a nick in the complementary strand. Here, we demonstrated that positioning the nick sgRNA nearby the templating sequences of the pegRNA facilitated targeted large fragment deletion and that engineering both guide RNAs to be pegRNAs to achieve bi-direction prime editing (Bi-PE) further increase the efficiency by up to 16 times and improved the accuracy of editing products by 60 times. In addition, we showed that Bi-PE strategy also increased the efficiency of simultaneous conversion of multiple bases but not single base conversion over PE3. In conclusion, Bi-PE strategy expanded the editing scope and improved the efficiency and the accuracy of prime editing system, which might have a wide range of potential applications.

ARMC5 is implicated in several pathological conditions, but its function remains unknown. We have previously identified CUL3 and RPB1 (the largest subunit of RNA polymerase II (Pol II) as potential ARMC5-interacting proteins. Here, we show that ARMC5, CUL3 and RBX1 form an active E3 ligase complex specific for RPB1. ARMC5, CUL3, and RBX1 formed an active E3 specific for RPB1. Armc5 deletion caused a significant reduction in RPB1 ubiquitination and an increase in an accumulation of RPB1, and hence an enlarged Pol II pool in normal tissues and organs. The compromised RPB1 degradation did not cause generalized Pol II stalling nor depressed transcription in the adrenal glands but did result in dysregulation of a subset of genes, with most upregulated. We found RPB1 to be highly expressed in the adrenal nodules from patients with primary bilateral macronodular adrenal hyperplasia (PBMAH) harboring germline ARMC5 mutations. Mutant ARMC5 had altered binding with RPB1. In summary, we discovered that wildtype ARMC5 was part of a novel RPB1-specific E3. ARMC5 mutations resulted in an enlarged Pol II pool, which dysregulated a subset of effector genes. Such an enlarged Pol II pool and gene dysregulation was correlated to adrenal hyperplasia in humans and KO mice.

At the time of writing, although siRNA therapeutics are approved for human use, no official regulatory guidance specific to this modality is available. In the absence of guidance, preclinical development for siRNA followed a hybrid of the small molecule and biologics guidance documents. However, siRNA differs significantly from small molecules and protein-based biologics in its physicochemical, absorption, distribution, metabolism and excretion properties, and its mechanism of action. Consequently, certain reports typically included in filing packages for small molecule or biologics may benefit from adaption, or even omission, from an siRNA filing. In this white paper, members of the 'siRNA working group' in the IQ Consortium compile a list of reports included in approved siRNA filing packages and discuss the relevance of two in vitro reports-the plasma protein binding evaluation and the drug-drug interaction risk assessment-to support siRNA regulatory filings. Publicly available siRNA approval packages and the literature were systematically reviewed to examine the role of siRNA plasma protein binding and drug-drug interactions in understanding pharmacokinetic/pharmacodynamic relationships, safety and translation. The findings are summarized into two decision trees to help guide industry decide when in vitro siRNA plasma protein binding and drug-drug interaction studies are warranted.

DNA mismatch repair removes mis-incorporated bases after DNA replication and reduces the error rate a 100-1000-fold. After recognition of a mismatch, a large section of up to a thousand nucleotides is removed from the daughter strand followed by re-synthesis. How these opposite activities are coordinated is poorly understood. Here we show that the Escherichia coli MutL protein binds to the 3' end of the resected strand and blocks access of Pol I and Pol III. The cryo-EM structure of an 85-kDa MutL-DNA complex, determined to 3.7 Å resolution, reveals a unique DNA binding mode that positions MutL at the 3' end of a primer-template, but not at a 5' resected DNA end or a blunt DNA end. Hence, our work reveals a novel role for MutL in the final stages of mismatch repair by preventing premature DNA synthesis during removal of the mismatched strand.

Box C/D small nucleolar RNAs (snoRNAs) are a conserved class of RNA known for their role in guiding ribosomal RNA 2'-O-ribose methylation. Recently, C/D snoRNAs were also implicated in regulating the expression of non-ribosomal genes through different modes of binding. Large scale RNA-RNA interaction datasets detect many snoRNAs binding messenger RNA, but are limited by specific experimental conditions. To enable a more comprehensive study of C/D snoRNA interactions, we created snoGloBe, a human C/D snoRNA interaction predictor based on a gradient boosting classifier. SnoGloBe considers the target type, position and sequence of the interactions, enabling it to outperform existing predictors. Interestingly, for specific snoRNAs, snoGloBe identifies strong enrichment of interactions near gene expression regulatory elements including splice sites. Abundance and splicing of predicted targets were altered upon the knockdown of their associated snoRNA. Strikingly, the predicted snoRNA interactions often overlap with the binding sites of functionally related RNA binding proteins, reinforcing their role in gene expression regulation. SnoGloBe is also an excellent tool for discovering viral RNA targets, as shown by its capacity to identify snoRNAs targeting the heavily methylated SARS-CoV-2 RNA. Overall, snoGloBe is capable of identifying experimentally validated binding sites and predicting novel sites with shared regulatory function.

In the presence of physiological monovalent cations, thousands of RNA G-rich sequences can form parallel G-quadruplexes (G4s) unless RNA-binding proteins inhibit, destabilize, or resolve the formation of such secondary RNA structures. Here, we have used a disease-relevant model system to investigate the biophysical properties of the RNA-binding protein HNRNPH1's interaction with G-rich sequences. We demonstrate the importance of two EWSR1-exon 8 G-rich regions in mediating the exclusion of this exon from the oncogenic EWS-FLI1 transcripts expressed in a subset of Ewing sarcomas, using complementary analysis of tumor data, long-read sequencing, and minigene studies. We determined that HNRNPH1 binds the EWSR1-exon 8 G-rich sequences with low nM affinities irrespective of whether in a non-G4 or G4 state but exhibits different kinetics depending on RNA structure. Specifically, HNRNPH1 associates and dissociates from G4-folded RNA faster than the identical sequences in a non-G4 state. Importantly, we demonstrate using gel shift and spectroscopic assays that HNRNPH1, particularly the qRRM1-qRRM2 domains, destabilizes the G4s formed by the EWSR1-exon 8 G-rich sequences in a non-catalytic fashion. Our results indicate that HNRNPH1's binding of G-rich sequences favors the accumulation of RNA in a non-G4 state and that this contributes to its regulation of RNA processing.

Since antibiotic resistance is often associated with a fitness cost, bacteria employ multi-layered regulatory mechanisms to ensure that expression of resistance factors is restricted to times of antibiotic challenge. In Bacillus subtilis, the chromosomally-encoded ABCF ATPase VmlR confers resistance to pleuromutilin, lincosamide and type A streptogramin translation inhibitors. Here we show that vmlR expression is regulated by translation attenuation and transcription attenuation mechanisms. Antibiotic-induced ribosome stalling during translation of an upstream open reading frame in the vmlR leader region prevents formation of an anti-antiterminator structure, leading to the formation of an antiterminator structure that prevents intrinsic termination. Thus, transcription in the presence of antibiotic induces vmlR expression. We also show that NusG-dependent RNA polymerase pausing in the vmlR leader prevents leaky expression in the absence of antibiotic. Furthermore, we demonstrate that induction of VmlR expression by compromised protein synthesis does not require the ability of VmlR to rescue the translational defect, as exemplified by constitutive induction of VmlR by ribosome assembly defects. Rather, the specificity of induction is determined by the antibiotic's ability to stall the ribosome on the regulatory open reading frame located within the vmlR leader. Finally, we demonstrate the involvement of (p)ppGpp-mediated signalling in antibiotic-induced VmlR expression.

The human primosome, a four-subunit complex of primase and DNA polymerase alpha (Polα), synthesizes chimeric RNA-DNA primers of a limited length for DNA polymerases delta and epsilon to initiate DNA replication on both chromosome strands. Despite recent structural insights into the action of its two catalytic centers, the mechanism of DNA synthesis termination is still unclear. Here we report results of functional and structural studies revealing how the human primosome counts RNA-DNA primer length and timely terminates DNA elongation. Using a single-turnover primer extension assay, we defined two factors that determine a mature primer length (∼35-mer): (i) a tight interaction of the C-terminal domain of the DNA primase large subunit (p58C) with the primer 5'-end, and (ii) flexible tethering of p58C and the DNA polymerase alpha catalytic core domain (p180core) to the primosome platform domain by extended linkers. The obtained data allow us to conclude that p58C is a key regulator of all steps of RNA-DNA primer synthesis. The above-described findings provide a notable insight into the mechanism of DNA synthesis termination by a eukaryotic primosome, an important process for ensuring successful primer handover to replication DNA polymerases and for maintaining genome integrity.

Homologous recombination (HR) serves multiple roles in DNA repair that are essential for maintaining genomic stability, including double-strand DNA break (DSB) repair. The central HR protein, RAD51, is frequently overexpressed in human malignancies, thereby elevating HR proficiency and promoting resistance to DNA-damaging therapies. Here, we find that the non-canonical NF-κB factors p100/52, but not RelB, control the expression of RAD51 in various human cancer subtypes. While p100/p52 depletion inhibits HR function in human tumor cells, it does not significantly influence the proficiency of non-homologous end joining, the other key mechanism of DSB repair. Clonogenic survival assays were performed using a pair DLD-1 cell lines that differ only in their expression of the key HR protein BRCA2. Targeted silencing of p100/p52 sensitizes the HR-competent cells to camptothecin, while sensitization is absent in HR-deficient control cells. These results suggest that p100/p52-dependent signaling specifically controls HR activity in cancer cells. Since non-canonical NF-κB signaling is known to be activated after various forms of genomic crisis, compensatory HR upregulation may represent a natural consequence of DNA damage. We propose that p100/p52-dependent signaling represents a promising oncologic target in combination with DNA-damaging treatments.

Translation initiation in eukaryotes starts with the recognition of the mRNA 5'-cap by eIF4F, a hetero-trimeric complex of eIF4E, the cap-binding protein, eIF4A, a DEAD-box helicase, and eIF4G, a scaffold protein. eIF4G comprises eIF4E- and eIF4A-binding domains (4E-BD, 4A-BD) and three RNA-binding regions (RNA1-RNA3), and interacts with eIF4A, eIF4E, and with the mRNA. Within the eIF4F complex, the helicase activity of eIF4A is increased. We showed previously that RNA3 of eIF4G is important for the stimulation of the eIF4A conformational cycle and its ATPase and helicase activities. Here, we dissect the interplay between the eIF4G domains and the role of the eIF4E/cap interaction in eIF4A activation. We show that RNA2 leads to an increase in the fraction of eIF4A in the closed state, an increased RNA affinity, and faster RNA unwinding. This stimulatory effect is partially reduced when the 4E-BD is present. eIF4E binding to the 4E-BD then further inhibits the helicase activity and closing of eIF4A, but does not affect the RNA-stimulated ATPase activity of eIF4A. The 5'-cap renders the functional interaction of mRNA with eIF4A less efficient. Overall, the activity of eIF4A at the 5'-cap is thus fine-tuned by a delicately balanced network of stimulatory and inhibitory interactions.

Chromosome 17q23 amplification occurs in 20% of primary breast tumours and is associated with poor outcome. The TBX2 gene is located on 17q23 and is often over-expressed in this breast tumour subset. TBX2 is an anti-senescence gene, promoting cell growth and survival through repression of Tumour Suppressor Genes (TSGs), such as NDRG1 and CST6. Previously we found that TBX2 cooperates with the PRC2 complex to repress several TSGs, and that PRC2 inhibition restored NDRG1 expression to impede cellular proliferation. Here, we now identify CoREST proteins, LSD1 and ZNF217, as novel interactors of TBX2. Genetic or pharmacological targeting of CoREST emulated TBX2 loss, inducing NDRG1 expression and abolishing breast cancer growth in vitro and in vivo. Furthermore, we uncover that TBX2/CoREST targeting of NDRG1 is achieved by recruitment of TBX2 to the NDRG1 promoter by Sp1, the abolishment of which resulted in NDRG1 upregulation and diminished cancer cell proliferation. Through ChIP-seq we reveal that 30% of TBX2-bound promoters are shared with ZNF217 and identify novel targets repressed by TBX2/CoREST; of these targets a lncRNA, LINC00111, behaves as a negative regulator of cell proliferation. Overall, these data indicate that inhibition of CoREST proteins represents a promising therapeutic intervention for TBX2-addicted breast tumours.

Poaceae plants can locally accumulate iron to suppress pathogen infection. It remains unknown how pathogens overcome host-derived iron stress during their successful infections. Here, we report that Fusarium graminearum (Fg), a destructive fungal pathogen of cereal crops, is challenged by host-derived high-iron stress. Fg infection induces host alkalinization, and the pH-dependent transcription factor FgPacC undergoes a proteolytic cleavage into the functional isoform named FgPacC30 under alkaline host environment. Subsequently FgPacC30 binds to a GCCAR(R = A/G)G element at the promoters of the genes involved in iron uptake and inhibits their expression, leading to adaption of Fg to high-iron stress. Mechanistically, FgPacC30 binds to FgGcn5 protein, a catalytic subunit of Spt-Ada-Gcn5 Acetyltransferase (SAGA) complex, leading to deregulation of histone acetylation at H3K18 and H2BK11, and repression of iron uptake genes. Moreover, we identified a protein kinase FgHal4, which is highly induced by extracellular high-iron stress and protects FgPacC30 against 26S proteasome-dependent degradation by promoting FgPacC30 phosphorylation at Ser2. Collectively, this study uncovers a novel inhibitory mechanism of the SAGA complex by a transcription factor that enables a fungal pathogen to adapt to dynamic microenvironments during infection.

Antisense peptide nucleic acids (PNAs) that target mRNAs of essential bacterial genes exhibit specific bactericidal effects in several microbial species, but our mechanistic understanding of PNA activity and their target gene spectrum is limited. Here, we present a systematic analysis of PNAs targeting 11 essential genes with varying expression levels in uropathogenic Escherichia coli (UPEC). We demonstrate that UPEC is susceptible to killing by peptide-conjugated PNAs, especially when targeting the widely-used essential gene acpP. Our evaluation yields three additional promising target mRNAs for effective growth inhibition, i.e.dnaB, ftsZ and rpsH. The analysis also shows that transcript abundance does not predict target vulnerability and that PNA-mediated growth inhibition is not universally associated with target mRNA depletion. Global transcriptomic analyses further reveal PNA sequence-dependent but also -independent responses, including the induction of envelope stress response pathways. Importantly, we show that 9mer PNAs are generally as effective in inhibiting bacterial growth as their 10mer counterparts. Overall, our systematic comparison of a range of PNAs targeting mRNAs of different essential genes in UPEC suggests important features for PNA design, reveals a general bacterial response to PNA conjugates and establishes the feasibility of using PNA antibacterials to combat UPEC.

Human telomere overhang composed of tandem repeats of TTAGGG folds into G-quadruplex (G4). Unlike in an experimental setting in the test tube in which the entire length is allowed to fold at once, inside the cell, the overhang is expected to fold as it is synthesized directionally (5' to 3') and released segmentally by a specialized enzyme, the telomerase. To mimic such vectorial G4 folding process, we employed a superhelicase, Rep-X which can unwind DNA to release the TTAGGG repeats in 5' to 3' direction. We demonstrate that the folded conformation achieved by the refolding of full sequence is significantly different from that of the vectorial folding for two to eight TTAGGG repeats. Strikingly, the vectorially folded state leads to a remarkably higher accessibility to complementary C-rich strand and the telomere binding protein POT1, reflecting a less stably folded state resulting from the vectorial folding. Importantly, our study points to an inherent difference between the co-polymerizing and post-polymerized folding of telomere overhang that can impact telomere architecture and downstream processes.

Methylation of cytosine to 5-methylcytosine (mC) at CpG sites is a prevalent reversible epigenetic mark in vertebrates established by DNA methyltransferases (MTases); the attached methyl groups can alter local structure of DNA and chromatin as well as binding of dedicated proteins. Nucleosome assembly on methylated DNA has been studied extensively, however little is known how the chromatin structure is affected by larger chemical variations in the major groove of DNA. Here, we studied the nucleosome formation in vitro on DNA containing an extended 5mC analog, 5-(6-azidohex-2-ynyl)cytosine (ahyC) installed at biological relevant CpG sites. We found that multiple ahyC residues on 80-Widom and Hsp70 promoter DNA fragments proved compatible with nucleosome assembly. Moreover, unlike mC, ahyC increases the affinity of histones to the DNA, partially altering nucleosome positioning, stability, and the action of chromatin remodelers. Based on molecular dynamics calculations, we suggest that these new features are due to increased DNA flexibility at ahyC-modified sites. Our findings provide new insights into the biophysical behavior of modified DNA and open new ways for directed design of synthetic nucleosomes.

HSUR1 and HSUR2, two noncoding RNAs expressed by the oncogenic Herpesvirus saimiri, bind host microRNAs miR-142-3p, miR-16, and miR-27 with different purposes. While binding of miR-27 to HSUR1 triggers the degradation of the microRNA, miR-16 is tethered by HSUR2 to target host mRNAs to repress their expression. Here we show that the interaction with miR-142-3p is required for the activity of both HSURs. Coimmunoprecipitation experiments revealed that miR-142-3p allosterically regulates the binding of miR-27 and miR-16 to HSUR1 and HSUR2, respectively. The binding of two different miRNAs to each HSUR is not cooperative. HSURs can be engineered to be regulated by other miRNAs, indicating that the identity of the binding miRNA is not important for HSUR regulation. Our results uncover a mechanism for allosteric regulation of noncoding RNA function and a previously unappreciated way in which microRNAs can regulate gene expression.

Alternative polyadenylation increases transcript diversities at the 3' end, regulating biological processes including cell differentiation, embryonic development and cancer progression. Here, we present a Bayesian method SCAPE, which enables de novo identification and quantification of polyadenylation (pA) sites at single-cell level by utilizing insert size information. We demonstrated its accuracy and robustness and identified 31 558 sites from 36 mouse organs, 43.8% (13 807) of which were novel. We illustrated that APA isoforms were associated with miRNAs binding and regulated in tissue-, cell type-and tumor-specific manners where no difference was found at gene expression level, providing an extra layer of information for cell clustering. Furthermore, we found genome-wide dynamic changes of APA usage during erythropoiesis and induced pluripotent stem cell (iPSC) differentiation, suggesting APA contributes to the functional flexibility and diversity of single cells. We expect SCAPE to aid the analyses of cellular dynamics and diversities in health and disease.

Alternative lengthening of telomeres (ALT) occurs in ∼10% of cancer entities. However, little is known about the heterogeneity of ALT activity since robust ALT detection assays with high-throughput in situ readouts are lacking. Here, we introduce ALT-FISH, a method to quantitate ALT activity in single cells from the accumulation of single-stranded telomeric DNA and RNA. It involves a one-step fluorescent in situ hybridization approach followed by fluorescence microscopy imaging. Our method reliably identified ALT in cancer cell lines from different tumor entities and was validated in three established models of ALT induction and suppression. Furthermore, we successfully applied ALT-FISH to spatially resolve ALT activity in primary tissue sections from leiomyosarcoma and neuroblastoma tumors. Thus, our assay provides insights into the heterogeneity of ALT tumors and is suited for high-throughput applications, which will facilitate screening for ALT-specific drugs.

Transcriptional regulators select their targets from a large pool of similar genomic sites. The binding of the Drosophila dosage compensation complex (DCC) exclusively to the male X chromosome provides insight into binding site selectivity rules. Previous studies showed that the male-specific organizer of the complex, MSL2, and ubiquitous DNA-binding protein CLAMP directly interact and play an important role in the specificity of X chromosome binding. Here, we studied the highly specific interaction between the intrinsically disordered region of MSL2 and the N-terminal zinc-finger C2H2-type (C2H2) domain of CLAMP. We obtained the NMR structure of the CLAMP N-terminal C2H2 zinc finger, which has a classic C2H2 zinc-finger fold with a rather unusual distribution of residues typically used in DNA recognition. Substitutions of residues in this C2H2 domain had the same effect on the viability of males and females, suggesting that it plays a general role in CLAMP activity. The N-terminal C2H2 domain of CLAMP is highly conserved in insects. However, the MSL2 region involved in the interaction is conserved only within the Drosophila genus, suggesting that this interaction emerged during the evolution of a mechanism for the specific recruitment of the DCC on the male X chromosome in Drosophilidae.

During eukaryotic ribosome biogenesis, pre-ribosomes travel from the nucleolus, where assembly is initiated, to the nucleoplasm and then are exported to the cytoplasm, where assembly concludes. Although nuclear export of pre-ribosomes has been extensively investigated, the release of pre-ribosomes from the nucleolus is an understudied phenomenon. Initial data indicate that unfolded rRNA interacts in trans with nucleolar components and that, when rRNA folds due to ribosomal protein (RP) binding, the number of trans interactions drops below the threshold necessary for nucleolar retention. To validate and expand on this idea, we performed a bioinformatic analysis of the protein components of the Saccharomyces cerevisiae ribosome assembly pathway. We found that ribosome biogenesis factors (RiBi factors) contain significantly more predicted trans interacting regions than RPs. We also analyzed cryo-EM structures of ribosome assembly intermediates to determine how nucleolar pre-ribosomes differ from post-nucleolar pre-ribosomes, specifically the capacity of RPs, RiBi factors, and rRNA components to interact in trans. We observed a significant decrease in the theoretical trans-interacting capability of pre-ribosomes between nucleolar and post-nucleolar stages of assembly due to the release of RiBi factors from particles and the folding of rRNA. Here, we provide a mechanism for the release of pre-ribosomes from the nucleolus.

Next-generation sequencing (NGS) enables the identification of functional nucleic acid sequences from in vitro selection/evolution experiments and illuminates the evolutionary process at single-nucleotide resolution. However, analyzing the vast output from NGS can be daunting, especially with limited programming skills. We developed REVERSE (Rapid EValuation of Experimental RNA Selection/Evolution) (https://www.reverseserver.org/), a web server that implements an integrated computational pipeline through a graphical user interface, which performs both pre-processing and detailed sequence level analyses within minutes. Raw FASTQ files are quality-filtered, dereplicated, and trimmed before being analyzed by either of two pipelines. The first pipeline counts, sorts, and tracks enrichment of unique sequences and user-defined sequence motifs. It also identifies mutational intermediates present in the sequence data that connect two input sequences. The second pipeline sorts similar sequences into clusters and tracks enrichment of peak sequences. It also performs nucleotide conservation analysis on the cluster of choice and generates a consensus sequence. Both pipelines generate downloadable spreadsheets and high-resolution figures. Collectively, REVERSE is a one-stop solution for the rapid analysis of NGS data obtained from in vitro selection/evolution experiments that obviates the need for computational expertise.

Elongation by RNA polymerase is dynamically modulated by accessory factors. The transcription-repair coupling factor (TRCF) recognizes paused/stalled RNAPs and either rescues transcription or initiates transcription termination. Precisely how TRCFs choose to execute either outcome remains unclear. With Escherichia coli as a model, we used single-molecule assays to study dynamic modulation of elongation by Mfd, the bacterial TRCF. We found that nucleotide-bound Mfd converts the elongation complex (EC) into a catalytically poised state, presenting the EC with an opportunity to restart transcription. After long-lived residence in this catalytically poised state, ATP hydrolysis by Mfd remodels the EC through an irreversible process leading to loss of the RNA transcript. Further, biophysical studies revealed that the motor domain of Mfd binds and partially melts DNA containing a template strand overhang. The results explain pathway choice determining the fate of the EC and provide a molecular mechanism for transcription modulation by TRCF.

The exon junction complex (EJC) is an RNA-binding multi-protein complex with critical functions in post-transcriptional gene regulation. It is deposited on the mRNA during splicing and regulates diverse processes including pre-mRNA splicing and nonsense-mediated mRNA decay (NMD) via various interacting proteins. The peripheral EJC-binding protein RNPS1 was reported to serve two insufficiently characterized functions: suppressing mis-splicing of cryptic splice sites and activating NMD in the cytoplasm. The analysis of transcriptome-wide effects of EJC and RNPS1 knockdowns in different human cell lines supports the conclusion that RNPS1 can moderately influence NMD activity, but is not a globally essential NMD factor. However, numerous aberrant splicing events strongly suggest that the main function of RNPS1 is splicing regulation. Rescue analyses revealed that the RRM and C-terminal domain of RNPS1 both contribute partially to regulate RNPS1-dependent splicing events. We defined the RNPS1 core interactome using complementary immunoprecipitations and proximity labeling, which identified interactions with splicing-regulatory factors that are dependent on the C-terminus or the RRM domain of RNPS1. Thus, RNPS1 emerges as a multifunctional splicing regulator that promotes correct and efficient splicing of different vulnerable splicing events via the formation of diverse splicing-promoting complexes.

Maternal-to-zygotic transition (MZT) is the first and key step in the control of animal development and intimately related to changes in chromatin structure and histone modifications. H2AK119ub1, an important epigenetic modification in regulating chromatin configuration and function, is primarily catalyzed by PRC1 and contributes to resistance to transcriptional reprogramming in mouse embryos. In this study, the genome-wide dynamic distribution of H2AK119ub1 during MZT in mice was investigated using chromosome immunoprecipitation and sequencing. The results indicated that H2AK119ub1 accumulated in fully grown oocytes and was enriched at the TSSs of maternal genes, but was promptly declined after meiotic resumption at genome-wide including the TSSs of early zygotic genes, by a previously unidentified mechanism. Genetic evidences indicated that ubiquitin-specific peptidase 16 (USP16) is the major deubiquitinase for H2AK119ub1 in mouse oocytes. Conditional knockout of Usp16 in oocytes did not impair their survival, growth, or meiotic maturation. However, oocytes lacking USP16 have defects when undergoing zygotic genome activation or gaining developmental competence after fertilization, potentially associated with high levels of maternal H2AK119ub1 deposition on the zygotic genomes. Taken together, H2AK119ub1 level is declined during oocyte maturation by an USP16-dependent mechanism, which ensures zygotic genome reprogramming and transcriptional activation of essential early zygotic genes.

Homologous recombination is the predominant DNA repair pathway used in the gonad. Of the excess DNA double-strand breaks formed in meiosis, only a subset matures into crossovers, with the remainder repaired as non-crossovers. The conserved BTR complex (comprising Bloom helicase, topoisomerase 3 and RMI1/2 scaffold proteins) acts at multiple steps during recombination to dismantle joint DNA molecules, thereby mediating the non-crossover outcome and chromosome integrity. Furthermore, the complex displays a role at the crossover site that is less well understood. Besides catalytic and TOPRIM domains, topoisomerase 3 enzymes contain a variable number of carboxy terminal zinc finger (ZnF) domains. Here, we studied the Caenorhabditis elegans mutant, in which the single ZnF domain is deleted. In contrast to the gene disruption allele, the top-3-ZnF mutant is viable, with no replication defects; the allele appears to be a hypomorph. The TOP-3-ZnF protein is recruited into foci but the mutant has increased numbers of crossovers along its chromosomes, with minor defects in repressing heterologous recombination, and a marked delay in the maturation/processing of recombination intermediates after loading of the RAD-51 recombinase. The ZnF domain cooperates with the RMI1 homolog RMH-2 to stabilize association of the BTR complex with recombination intermediates and to prevent recombination between heterologous DNA sequences.

The spread of drug-resistant bacteria represents one of the most significant medical problems of our time. Bacterial fitness loss associated with drug resistance can be counteracted by acquisition of secondary mutations, thereby enhancing the virulence of such bacteria. Antibiotic rifampicin (Rif) targets cellular RNA polymerase (RNAP). It is potent broad spectrum drug used for treatment of bacterial infections. We have investigated the compensatory mechanism of the secondary mutations alleviating Rif resistance (Rifr) on biochemical, structural and fitness indices. We find that substitutions in RNAP genes compensating for the growth defect caused by βQ513P and βT563P Rifr mutations significantly enhanced bacterial relative growth rate. By assaying RNAP purified from these strains, we show that compensatory mutations directly stimulated basal transcriptional machinery (2-9-fold) significantly improving promoter clearance step of the transcription pathway as well as elongation rate. Molecular modeling suggests that compensatory mutations affect transcript retention, substrate loading, and nucleotidyl transfer catalysis. Strikingly, one of the identified compensatory substitutions represents mutation conferring rifampicin resistance on its own. This finding reveals an evolutionary process that creates more virulent species by simultaneously improving the fitness and augmenting bacterial drug resistance.

Human Long Intergenic Noncoding RNA-p21 (LincRNA-p21) is a regulatory noncoding RNA that plays an important role in promoting apoptosis. LincRNA-p21 is also critical in down-regulating many p53 target genes through its interaction with a p53 repressive complex. The interaction between LincRNA-p21 and the repressive complex is likely dependent on the RNA tertiary structure. Previous studies have determined the two-dimensional secondary structures of the sense and antisense human LincRNA-p21 AluSx1 IRs using SHAPE. However, there were no insights into its three-dimensional structure. Therefore, we in vitro transcribed the sense and antisense regions of LincRNA-p21 AluSx1 Inverted Repeats (IRs) and performed analytical ultracentrifugation, size exclusion chromatography, light scattering, and small angle X-ray scattering (SAXS) studies. Based on these studies, we determined low-resolution, three-dimensional structures of sense and antisense LincRNA-p21. By adapting previously known two-dimensional information, we calculated their sense and antisense high-resolution models and determined that they agree with the low-resolution structures determined using SAXS. Thus, our integrated approach provides insights into the structure of LincRNA-p21 Alu IRs. Our study also offers a viable pipeline for combining the secondary structure information with biophysical and computational studies to obtain high-resolution atomistic models for long noncoding RNAs.

Deep sequencing has become one of the most popular tools for transcriptome profiling in biomedical studies. While an abundance of computational methods exists for 'normalizing' sequencing data to remove unwanted between-sample variations due to experimental handling, there is no consensus on which normalization is the most suitable for a given data set. To address this problem, we developed 'DANA'-an approach for assessing the performance of normalization methods for microRNA sequencing data based on biology-motivated and data-driven metrics. Our approach takes advantage of well-known biological features of microRNAs for their expression pattern and chromosomal clustering to simultaneously assess (i) how effectively normalization removes handling artifacts and (ii) how aptly normalization preserves biological signals. With DANA, we confirm that the performance of eight commonly used normalization methods vary widely across different data sets and provide guidance for selecting a suitable method for the data at hand. Hence, it should be adopted as a routine preprocessing step (preceding normalization) for microRNA sequencing data analysis. DANA is implemented in R and publicly available at https://github.com/LXQin/DANA.

Non-homologous end joining (NHEJ) is the major pathway that mediates the repair of DNA double-strand breaks (DSBs) generated by ionizing radiation (IR). Previously, the DNA helicase RECQL4 was implicated in promoting NHEJ, but its role in the pathway remains unresolved. In this study, we report that RECQL4 stabilizes the NHEJ machinery at DSBs to promote repair. Specifically, we find that RECQL4 interacts with the NHEJ core factor DNA-PKcs and the interaction is increased following IR. RECQL4 promotes DNA end bridging mediated by DNA-PKcs and Ku70/80 in vitro and the accumulation/retention of NHEJ factors at DSBs in vivo. Moreover, interaction between DNA-PKcs and the other core NHEJ proteins following IR treatment is attenuated in the absence of RECQL4. These data indicate that RECQL4 promotes the stabilization of the NHEJ factors at DSBs to support formation of the NHEJ long-range synaptic complex. In addition, we observed that the kinase activity of DNA-PKcs is required for accumulation of RECQL4 to DSBs and that DNA-PKcs phosphorylates RECQL4 at six serine/threonine residues. Blocking phosphorylation at these sites reduced the recruitment of RECQL4 to DSBs, attenuated the interaction between RECQL4 and NHEJ factors, destabilized interactions between the NHEJ machinery, and resulted in decreased NHEJ. Collectively, these data illustrate reciprocal regulation between RECQL4 and DNA-PKcs in NHEJ.

UPF3 is a key nonsense-mediated mRNA decay (NMD) factor required for mRNA surveillance and eukaryotic gene expression regulation. UPF3 exists as two paralogs (A and B) which are differentially expressed depending on cell type and developmental stage and believed to regulate NMD activity based on cellular requirements. UPF3B mutations cause intellectual disability. The underlying molecular mechanisms remain elusive, as many of the mutations lie in the poorly characterized middle-domain of UPF3B. Here, we show that UPF3A and UPF3B share structural and functional homology to paraspeckle proteins comprising an RNA-recognition motif-like domain (RRM-L), a NONA/paraspeckle-like domain (NOPS-L), and extended α-helical domain. These domains are essential for RNA/ribosome-binding, RNA-induced oligomerization and UPF2 interaction. Structures of UPF2's third middle-domain of eukaryotic initiation factor 4G (MIF4GIII) in complex with either UPF3B or UPF3A reveal unexpectedly intimate binding interfaces. UPF3B's disease-causing mutation Y160D in the NOPS-L domain displaces Y160 from a hydrophobic cleft in UPF2 reducing the binding affinity ∼40-fold compared to wildtype. UPF3A, which is upregulated in patients with the UPF3B-Y160D mutation, binds UPF2 with ∼10-fold higher affinity than UPF3B reliant mainly on NOPS-L residues. Our characterization of RNA- and UPF2-binding by UPF3's middle-domain elucidates its essential role in NMD.

The structure and properties of DNA depend on the environment, in particular the ion atmosphere. Here, we investigate how DNA twist -one of the central properties of DNA- changes with concentration and identity of the surrounding ions. To resolve how cations influence the twist, we combine single-molecule magnetic tweezer experiments and extensive all-atom molecular dynamics simulations. Two interconnected trends are observed for monovalent alkali and divalent alkaline earth cations. First, DNA twist increases monotonously with increasing concentration for all ions investigated. Second, for a given salt concentration, DNA twist strongly depends on cation identity. At 100 mM concentration, DNA twist increases as Na+ < K+ < Rb+ < Ba2+ < Li+ ≈ Cs+ < Sr2+ < Mg2+ < Ca2+. Our molecular dynamics simulations reveal that preferential binding of the cations to the DNA backbone or the nucleobases has opposing effects on DNA twist and provides the microscopic explanation of the observed ion specificity. However, the simulations also reveal shortcomings of existing force field parameters for Cs+ and Sr2+. The comprehensive view gained from our combined approach provides a foundation for understanding and predicting cation-induced structural changes both in nature and in DNA nanotechnology.

MicroRNA-dependent mRNA decay plays an important role in gene silencing by facilitating posttranscriptional and translational repression. Inspired by this intrinsic nature of microRNA-mediated mRNA cleavage, here, we describe a microRNA-targeting mRNA as a switch platform called mRNA bridge mimetics to regulate the translocation of proteins. We applied the mRNA bridge mimetics platform to Cas9 protein to confer it the ability to translocate into the nucleus via cleavage of the nuclear export signal. This system performed programmed gene editing in vitro and in vivo. Combinatorial treatment with cisplatin and miR-21-EZH2 axis-targeting CRISPR Self Check-In improved sensitivity to chemotherapeutic drugs in vivo. Using the endogenous microRNA-mediated mRNA decay mechanism, our platform is able to remodel a cell's natural biology to allow the entry of precise drugs into the nucleus, devoid of non-specific translocation. The mRNA bridge mimetics strategy is promising for applications in which the reaction must be controlled via intracellular stimuli and modulates Cas9 proteins to ensure safe genome modification in diseased conditions.

In the last decade, multiple studies demonstrated that cells maintain a balance of mRNA production and degradation, but the mechanisms by which cells implement this balance remain unknown. Here, we monitored cells' total and recently-transcribed mRNA profiles immediately following an acute depletion of Xrn1-the main 5'-3' mRNA exonuclease-which was previously implicated in balancing mRNA levels. We captured the detailed dynamics of the adaptation to rapid degradation of Xrn1 and observed a significant accumulation of mRNA, followed by a delayed global reduction in transcription and a gradual return to baseline mRNA levels. We found that this transcriptional response is not unique to Xrn1 depletion; rather, it is induced earlier when upstream factors in the 5'-3' degradation pathway are perturbed. Our data suggest that the mRNA feedback mechanism monitors the accumulation of inputs to the 5'-3' exonucleolytic pathway rather than its outputs.

Emerging single-cell technologies provide high-resolution measurements of distinct cellular modalities opening new avenues for generating detailed cellular atlases of many and diverse tissues. The high dimensionality, sparsity, and inaccuracy of single cell sequencing measurements, however, can obscure discriminatory information, mask cellular subtype variations and complicate downstream analyses which can limit our understanding of cell function and tissue heterogeneity. Here, we present a novel pre-processing method (scPSD) inspired by power spectral density analysis that enhances the accuracy for cell subtype separation from large-scale single-cell omics data. We comprehensively benchmarked our method on a wide range of single-cell RNA-sequencing datasets and showed that scPSD pre-processing, while being fast and scalable, significantly reduces data complexity, enhances cell-type separation, and enables rare cell identification. Additionally, we applied scPSD to transcriptomics and chromatin accessibility cell atlases and demonstrated its capacity to discriminate over 100 cell types across the whole organism and across different modalities of single-cell omics data.

DDX58 encodes RIG-I, a cytosolic RNA sensor that ensures immune surveillance of nonself RNAs. Individuals with RIG-IE510V and RIG-IQ517H mutations have increased susceptibility to Singleton-Merten syndrome (SMS) defects, resulting in tissue-specific (mild) and classic (severe) phenotypes. The coupling between RNA recognition and conformational changes is central to RIG-I RNA proofreading, but the molecular determinants leading to dissociated disease phenotypes remain unknown. Herein, we employed hydrogen/deuterium exchange mass spectrometry (HDX-MS) and single molecule magnetic tweezers (MT) to precisely examine how subtle conformational changes in the helicase insertion domain (HEL2i) promote impaired ATPase and erroneous RNA proofreading activities. We showed that the mutations cause a loosened latch-gate engagement in apo RIG-I, which in turn gradually dampens its self RNA (Cap2 moiety:m7G cap and N1-2-2'-O-methylation RNA) proofreading ability, leading to increased immunopathy. These results reveal HEL2i as a unique checkpoint directing two specialized functions, i.e. stabilizing the CARD2-HEL2i interface and gating the helicase from incoming self RNAs; thus, these findings add new insights into the role of HEL2i in the control of antiviral innate immunity and autoimmunity diseases.

Comment in 10.1093/nar/gkac018.

The cell-wall recycling process is important for bacterial survival in nutrient-limited conditions and, in certain cases, is directly involved in antibiotic resistance. In the sophisticated cell-wall recycling process in Escherichia coli, the transcriptional repressor MurR controls the expression of murP and murQ, which are involved in transporting and metabolizing N-acetylmuramic acid (MurNAc), generating N-acetylmuramic acid-6-phosphate (MurNAc-6-P) and N-acetylglucosamine-6-phosphate (GlcNAc-6-P). Here, we report that both MurNAc-6-P and GlcNAc-6-P can bind to MurR and weaken the DNA binding ability of MurR. Structural characterizations of MurR in complex with MurNAc-6-P or GlcNAc-6-P as well as in the apo form revealed the detailed ligand recognition chemistries. Further studies showed that only MurNAc-6-P, but not GlcNAc-6-P, is capable of derepressing the expression of murQP controlled by MurR in cells and clarified the substrate specificity through the identification of key residues responsible for ligand binding in the complex structures. In summary, this study deciphered the molecular mechanism of the cell wall recycling process regulated by MurR in E. coli.

Triplex-forming oligonucleotides (TFOs) are short, single-stranded oligomers that hybridise to a specific sequence of duplex DNA. TFOs can block transcription and thereby inhibit protein production, making them highly appealing in the field of antigene therapeutics. In this work, a primer extension protocol was developed to enzymatically prepare chemical nuclease TFO hybrid constructs, with gene-silencing applications. Click chemistry was employed to generate novel artificial metallo-nuclease (AMN)-dNTPs, which were selectively incorporated into the TFO strand by a DNA polymerase. This purely enzymatic protocol was then extended to facilitate the construction of 5-methylcytosine (5mC) modified TFOs that displayed increased thermal stability. The utility of the enzymatically synthesised di-(2-picolyl)amine (DPA)-TFOs was assessed and compared to a specifically prepared solid-phase synthesis counterpart through gel electrophoresis, quantitative PCR, and Sanger sequencing, which revealed similar recognition and damage properties to target genes. The specificity was then enhanced through coordinated designer intercalators-DPQ and DPPZ-and high-precision DNA cleavage was achieved. To our knowledge, this is the first example of the enzymatic production of an AMN-TFO hybrid and is the largest base modification incorporated using this method. These results indicate how chemical nuclease-TFOs may overcome limitations associated with non-molecularly targeted metallodrugs and open new avenues for artificial gene-editing technology.

The segmented negative-sense RNA genome of influenza A virus is assembled into ribonucleoprotein complexes (RNP) with viral RNA-dependent RNA polymerase and nucleoprotein (NP). It is in the context of these RNPs that the polymerase transcribes and replicates viral RNA (vRNA). Host acidic nuclear phosphoprotein 32 (ANP32) family proteins play an essential role in vRNA replication by mediating the dimerization of the viral polymerase via their N-terminal leucine-rich repeat (LRR) domain. However, whether the C-terminal low-complexity acidic region (LCAR) plays a role in RNA synthesis remains unknown. Here, we report that the LCAR is required for viral genome replication during infection. Specifically, we show that the LCAR directly interacts with NP and this interaction is mutually exclusive with RNA. Furthermore, we show that the replication of a short vRNA-like template that can be replicated in the absence of NP is less sensitive to LCAR truncations compared with the replication of full-length vRNA segments which is NP-dependent. We propose a model in which the LCAR interacts with NP to promote NP recruitment to nascent RNA during influenza virus replication, ensuring the co-replicative assembly of RNA into RNPs.

Efficient splicing requires a balance between high-fidelity splice-site (SS) selection and speed. In Saccharomyces cerevisiae, Pre-mRNA processing factor 8 (Prp8) helps to balance precise SS selection and rapid, efficient intron excision and exon joining. argonaute1-52 (ago1-52) and incurvata13 (icu13) are hypomorphic alleles of the Arabidopsis thaliana genes ARGONAUTE1 (AGO1) and AUXIN RESISTANT6 (AXR6) that harbor point mutations creating a novel 3'SS and 5'SS, respectively. The spliceosome recognizes these novel SSs, as well as the intact genuine SSs, producing a mixture of wild-type and aberrant mature mRNAs. Here, we characterized five novel mutant alleles of PRP8 (one of the two Arabidopsis co-orthologs of yeast Prp8), naming these alleles morphology of ago1-52 suppressed5 (mas5). In the mas5-1 background, the spliceosome preferentially recognizes the intact genuine 3'SS of ago1-52 and 5'SS of icu13. Since point mutations that damage genuine SSs make the spliceosome prone to recognizing cryptic SSs, we also tested alleles of four genes carrying damaged genuine SSs, finding that mas5-1 did not suppress their missplicing. The mas5-1 and mas5-3 mutations represent a novel class of missplicing suppressors that increase splicing fidelity by hampering the use of novel SSs, but do not alter general pre-mRNA splicing.

T-box riboswitches (T-boxes) are essential RNA regulatory elements with a remarkable structural diversity, especially among bacterial pathogens. In staphylococci, all glyS T-boxes synchronize glycine supply during synthesis of nascent polypeptides and cell wall formation and are characterized by a conserved and unique insertion in their antiterminator/terminator domain, termed stem Sa. Interestingly, in Staphylococcus aureus the stem Sa can accommodate binding of specific antibiotics, which in turn induce robust and diverse effects on T-box-mediated transcription. In the present study, domain swap mutagenesis and probing analysis were performed to decipher the role of stem Sa. Deletion of stem Sa significantly reduces both the S. aureus glyS T-box-mediated transcription readthrough levels and the ability to discriminate among tRNAGly isoacceptors, both in vitro and in vivo. Moreover, the deletion inverted the previously reported stimulatory effects of specific antibiotics. Interestingly, stem Sa insertion in the terminator/antiterminator domain of Geobacillus kaustophilus glyS T-box, which lacks this domain, resulted in elevated transcription in the presence of tigecycline and facilitated discrimination among proteinogenic and nonproteinogenic tRNAGly isoacceptors. Overall, stem Sa represents a lineage-specific structural feature required for efficient staphylococcal glyS T-box-mediated transcription and it could serve as a species-selective druggable target through its ability to modulate antibiotic binding.

The assembly of high molecular mass ribonucleoprotein complexes typically relies on the binary interaction of defined RNA sequences or precisely folded RNA motifs with dedicated RNA-binding domains on the protein side. Here we describe a new molecular recognition principle of RNA molecules by a high molecular mass protein complex. By chemically probing the solvent accessibility of mitochondrial pre-mRNAs when bound to the Trypanosoma brucei editosome, we identified multiple similar but non-identical RNA motifs as editosome contact sites. However, by treating the different motifs as mathematical graph objects we demonstrate that they fit a consensus 2D-graph consisting of 4 vertices (V) and 3 edges (E) with a Laplacian eigenvalue of 0.5477 (λ2). We establish that synthetic 4V(3E)-RNAs are sufficient to compete for the editosomal pre-mRNA binding site and that they inhibit RNA editing in vitro. Furthermore, we demonstrate that only two topological indices are necessary to predict the binding of any RNA motif to the editosome with a high level of confidence. Our analysis corroborates that the editosome has adapted to the structural multiplicity of the mitochondrial mRNA folding space by recognizing a fuzzy continuum of RNA folds that fit a consensus graph descriptor.

Chemical synthesis of RNA conjugates has opened new strategies to study enzymatic mechanisms in RNA biology. To gain insights into poorly understood RNA nucleotide methylation processes, we developed a new method to synthesize RNA-conjugates for the study of RNA recognition and methyl-transfer mechanisms of SAM-dependent m6A RNA methyltransferases. These RNA conjugates contain a SAM cofactor analogue connected at the N6-atom of an adenosine within dinucleotides, a trinucleotide or a 13mer RNA. Our chemical route is chemo- and regio-selective and allows flexible modification of the RNA length and sequence. These compounds were used in crystallization assays with RlmJ, a bacterial m6A rRNA methyltransferase. Two crystal structures of RlmJ in complex with RNA-SAM conjugates were solved and revealed the RNA-specific recognition elements used by RlmJ to clamp the RNA substrate in its active site. From these structures, a model of a trinucleotide bound in the RlmJ active site could be built and validated by methyltransferase assays on RlmJ mutants. The methyl transfer by RlmJ could also be deduced. This study therefore shows that RNA-cofactor conjugates are potent molecular tools to explore the active site of RNA modification enzymes.

Generation of haploid gametes depends on a modified version of homologous recombination in meiosis. Meiotic recombination is initiated by single-stranded DNA (ssDNA) ends originating from programmed DNA double-stranded breaks (DSBs) that are generated by the topoisomerase-related SPO11 enzyme. Meiotic recombination involves chromosomal synapsis, which enhances recombination-mediated DSB repair, and thus, crucially contributes to genome maintenance in meiocytes. Synapsis defects induce oocyte apoptosis ostensibly due to unrepaired DSBs that persist in asynaptic chromosomes. In mice, SPO11-deficient oocytes feature asynapsis, apoptosis and, surprisingly, numerous foci of the ssDNA-binding recombinase RAD51, indicative of DSBs of unknown origin. Hence, asynapsis is suggested to trigger apoptosis due to inefficient DSB repair even in mutants that lack programmed DSBs. By directly detecting ssDNAs, we discovered that RAD51 is an unreliable marker for DSBs in oocytes. Further, SPO11-deficient oocytes have fewer persistent ssDNAs than wild-type oocytes. These observations suggest that oocyte quality is safeguarded in mammals by a synapsis surveillance mechanism that can operate without persistent ssDNAs.

Axonally synthesized proteins support nerve regeneration through retrograde signaling and local growth mechanisms. RNA binding proteins (RBP) are needed for this and other aspects of post-transcriptional regulation of neuronal mRNAs, but only a limited number of axonal RBPs are known. We used targeted proteomics to profile RBPs in peripheral nerve axons. We detected 76 proteins with reported RNA binding activity in axoplasm, and levels of several change with axon injury and regeneration. RBPs with altered levels include KHSRP that decreases neurite outgrowth in developing CNS neurons. Axonal KHSRP levels rapidly increase after injury remaining elevated up to 28 days post axotomy. Khsrp mRNA localizes into axons and the rapid increase in axonal KHSRP is through local translation of Khsrp mRNA in axons. KHSRP can bind to mRNAs with 3'UTR AU-rich elements and targets those transcripts to the cytoplasmic exosome for degradation. KHSRP knockout mice show increased axonal levels of KHSRP target mRNAs, Gap43, Snap25, and Fubp1, following sciatic nerve injury and these mice show accelerated nerve regeneration in vivo. Together, our data indicate that axonal translation of the RNA binding protein Khsrp mRNA following nerve injury serves to promote decay of other axonal mRNAs and slow axon regeneration.

Biomolecular associations forged by specific interaction among structural scaffolds are fundamental to the control and regulation of cell processes. One such structural architecture, characterized by HEAT repeats, is involved in a multitude of cellular processes, including intracellular transport, signaling, and protein synthesis. Here, we review the multitude and versatility of HEAT domains in the regulation of mRNA translation initiation. Structural and cellular biology approaches, as well as several biophysical studies, have revealed that a number of HEAT domain-mediated interactions with a host of protein factors and RNAs coordinate translation initiation. We describe the basic structural architecture of HEAT domains and briefly introduce examples of the cellular processes they dictate, including nuclear transport by importin and RNA degradation. We then focus on proteins in the translation initiation system featuring HEAT domains, specifically the HEAT domains of eIF4G, DAP5, eIF5, and eIF2Bϵ. Comparative analysis of their remarkably versatile interactions, including protein-protein and protein-RNA recognition, reveal the functional importance of flexible regions within these HEAT domains. Here we outline how HEAT domains orchestrate fundamental aspects of translation initiation and highlight open mechanistic questions in the area.

Advances in transcriptomic and translatomic techniques enable in-depth studies of RNA activity profiles and RNA-based regulatory mechanisms. Ribosomal RNA (rRNA) sequences are highly abundant among cellular RNA, but if the target sequences do not include polyadenylation, these cannot be easily removed in library preparation, requiring their post-hoc removal with computational techniques to accelerate and improve downstream analyses. Here, we describe RiboDetector, a novel software based on a Bi-directional Long Short-Term Memory (BiLSTM) neural network, which rapidly and accurately identifies rRNA reads from transcriptomic, metagenomic, metatranscriptomic, noncoding RNA, and ribosome profiling sequence data. Compared with state-of-the-art approaches, RiboDetector produced at least six times fewer misclassifications on the benchmark datasets. Importantly, the few false positives of RiboDetector were not enriched in certain Gene Ontology (GO) terms, suggesting a low bias for downstream functional profiling. RiboDetector also demonstrated a remarkable generalizability for detecting novel rRNA sequences that are divergent from the training data with sequence identities of <90%. On a personal computer, RiboDetector processed 40M reads in less than 6 min, which was ∼50 times faster in GPU mode and ∼15 times in CPU mode than other methods. RiboDetector is available under a GPL v3.0 license at https://github.com/hzi-bifo/RiboDetector.

Directed evolution is a powerful method for engineering biology in the absence of detailed sequence-function relationships. To enable directed evolution of complex phenotypes encoded by multigene pathways, we require large library sizes for DNA sequences >5-10 kb in length, elimination of genomic hitchhiker mutations, and decoupling of diversification and screening steps. To meet these challenges, we developed Inducible Directed Evolution (IDE), which uses a temperate bacteriophage to package large plasmids and transfer them to naive cells after intracellular mutagenesis. To demonstrate IDE, we evolved a 5-gene pathway from Bacillus licheniformis that accelerates tagatose catabolism in Escherichia coli, resulting in clones with 65% shorter lag times during growth on tagatose after only two rounds of evolution. Next, we evolved a 15.4 kb, 10-gene pathway from Bifidobacterium breve UC2003 that aids E. coli's utilization of melezitose. After three rounds of IDE, we isolated evolved pathways that both reduced lag time by more than 2-fold and enabled 150% higher final optical density. Taken together, this work enhances the capacity and utility of a whole pathway directed evolution approach in E. coli.

Comment in 10.1093/nar/gkac037.

Author information: (1)Aix-Marseille Université, Université de Toulon, IRD, CNRS, Mediterranean Institute of Oceanography (MIO) UM 110, Marseille, France. (2)Research Federation for the study of Global Ocean systems ecology and evolution, FR2022/Tara Oceans-GOSEE, Paris, France. (3)SIP, OSU PYTHEAS, Marseille, France. (4)Department of Marine Biology and Oceanography, Institute of Marine Sciences (ICM), CSIC, Barcelona, Spain. (5)Department of Biology, Institute of Microbiology and Swiss Institute of Bioinformatics, ETH Zurich, Zurich, Switzerland. (6)Génomique Métabolique, Genoscope, Institut de Biologie François-Jacob, CEA, CNRS, Univ Evry, Univ Paris-Saclay, 91057 Evry, France.

Saccharomyces cerevisiae is an important model eukaryotic microorganism and widely applied in fundamental research and the production of various chemicals. Its ability to efficiently and precisely control the expression of multiple genes is valuable for metabolic engineering. The clustered regularly interspaced short palindromic repeats (CRISPR)-mediated regulation enables complex gene expression programming; however, the regulation efficiency is often limited by the efficiency of pertinent regulators. Here, we developed CRISPR-mediated protein-tagging signal amplification system for simultaneous multiplexed gene activation and repression in S. cerevisiae. By introducing protein scaffolds (SPY and SunTag systems) to recruit multiple copies of regulators to different nuclease-deficient CRISPR proteins and design optimization, our system amplified gene regulation efficiency significantly. The gene activation and repression efficiencies reached as high as 34.9-fold and 95%, respectively, being 3.8- and 8.6-fold higher than those observed on the direct fusion of regulators with nuclease-deficient CRISPR proteins, respectively. We then applied the orthogonal bifunctional CRISPR-mediated transcriptional regulation system to regulate the expression of genes associated with 3-hydroxypropanoic acid production to deduce that CRISPR-associated regulator recruiting systems represent a robust method for simultaneously regulating multiple genes and rewiring metabolic pathways.

Rob, which serves as a paradigm of the large AraC/XylS family transcription activators, regulates diverse subsets of genes involved in multidrug resistance and stress response. However, the underlying mechanism of how it engages bacterial RNA polymerase and promoter DNA to finely respond to environmental stimuli is still elusive. Here, we present two cryo-EM structures of Rob-dependent transcription activation complex (Rob-TAC) comprising of Escherichia coli RNA polymerase (RNAP), Rob-regulated promoter and Rob in alternative conformations. The structures show that a single Rob engages RNAP by interacting with RNAP αCTD and σ70R4, revealing their generally important regulatory roles. Notably, by occluding σ70R4 from binding to -35 element, Rob specifically binds to the conserved Rob binding box through its consensus HTH motifs, and retains DNA bending by aid of the accessory acidic loop. More strikingly, our ligand docking and biochemical analysis demonstrate that the large Rob C-terminal domain (Rob CTD) shares great structural similarity with the global Gyrl-like domains in effector binding and allosteric regulation, and coordinately promotes formation of competent Rob-TAC. Altogether, our structural and biochemical data highlight the detailed molecular mechanism of Rob-dependent transcription activation, and provide favorable evidences for understanding the physiological roles of the other AraC/XylS-family transcription factors.

Replication fork reversal occurs via a two-step process that entails reversal initiation and reversal extension. DNA topoisomerase IIalpha (TOP2A) facilitates extensive fork reversal, on one hand through resolving the topological stress generated by the initial reversal, on the other hand via its role in recruiting the SUMO-targeted DNA translocase PICH to stalled forks in a manner that is dependent on its SUMOylation by the SUMO E3 ligase ZATT. However, how TOP2A activities at stalled forks are precisely regulated remains poorly understood. Here we show that, upon replication stress, the SUMO-targeted ubiquitin E3 ligase RNF4 accumulates at stalled forks and targets SUMOylated TOP2A for ubiquitination and degradation. Downregulation of RNF4 resulted in aberrant activation of the ZATT-TOP2A-PICH complex at stalled forks, which in turn led to excessive reversal and elevated frequencies of fork collapse. These results uncover a previously unidentified regulatory mechanism that regulates TOP2A activities at stalled forks and thus the extent of fork reversal.

Transcription elongation factor Spt6 associates with RNA polymerase II (Pol II) and acts as a histone chaperone, which promotes the reassembly of nucleosomes following the passage of Pol II. The precise mechanism of nucleosome reassembly mediated by Spt6 remains unclear. In this study, we used a hybrid approach combining cryo-electron microscopy and small-angle X-ray scattering to visualize the architecture of Spt6 from Saccharomyces cerevisiae. The reconstructed overall architecture of Spt6 reveals not only the core of Spt6, but also its flexible N- and C-termini, which are critical for Spt6's function. We found that the acidic N-terminal region of Spt6 prevents the binding of Spt6 not only to the Pol II CTD and Pol II CTD-linker, but also to pre-formed intact nucleosomes and nucleosomal DNA. The N-terminal region of Spt6 self-associates with the tSH2 domain and the core of Spt6 and thus controls binding to Pol II and nucleosomes. Furthermore, we found that Spt6 promotes the assembly of nucleosomes in vitro. These data indicate that the cooperation between the intrinsically disordered and structured regions of Spt6 regulates nucleosome and Pol II CTD binding, and also nucleosome assembly.

A major pharmacological strategy toward HIV cure aims to reverse latency in infected cells as a first step leading to their elimination. While the unbiased identification of molecular targets physically associated with the latent HIV-1 provirus would be highly valuable to unravel the molecular determinants of HIV-1 transcriptional repression and latency reversal, due to technical limitations, this has been challenging. Here we use a dCas9 targeted chromatin and histone enrichment strategy coupled to mass spectrometry (Catchet-MS) to probe the differential protein composition of the latent and activated HIV-1 5'LTR. Catchet-MS identified known and novel latent 5'LTR-associated host factors. Among these, IKZF1 is a novel HIV-1 transcriptional repressor, required for Polycomb Repressive Complex 2 recruitment to the LTR. We find the clinically advanced thalidomide analogue iberdomide, and the FDA approved analogues lenalidomide and pomalidomide, to be novel LRAs. We demonstrate that, by targeting IKZF1 for degradation, these compounds reverse HIV-1 latency in CD4+ T-cells isolated from virally suppressed people living with HIV-1 and that they are able to synergize with other known LRAs.

Heterochromatic silencing is thought to occur through a combination of transcriptional silencing and RNA degradation, but the relative contribution of each pathway is not known. In this study, we analyzed RNA Polymerase II (RNA Pol II) occupancy and levels of nascent and steady-state RNA in different mutants of Schizosaccharomyces pombe, in order to quantify the contribution of each pathway to heterochromatic silencing. We found that transcriptional silencing consists of two components, reduced RNA Pol II accessibility and, unexpectedly, reduced transcriptional efficiency. Heterochromatic loci showed lower transcriptional output compared to euchromatic loci, even when comparable amounts of RNA Pol II were present in both types of regions. We determined that the Ccr4-Not complex and H3K9 methylation are required for reduced transcriptional efficiency in heterochromatin and that a subset of heterochromatic RNA is degraded more rapidly than euchromatic RNA. Finally, we quantified the contribution of different chromatin modifiers, RNAi and RNA degradation to each silencing pathway. Our data show that several pathways contribute to heterochromatic silencing in a locus-specific manner and reveal transcriptional efficiency as a new mechanism of silencing.

Synechococcus elongatus, formerly known as Anacystis nidulans, is a representative species of cyanobacteria. It is also a model organism for the study of photoreactivation, which can be fully photoreactivated even after receiving high UV doses. However, for a long time, only one photolyase was found in S. elongatus that is only able to photorepair UV induced cyclobutane pyrimidine dimers (CPDs) in DNA. Here, we characterize another photolyase in S. elongatus, which belongs to iron-sulfur bacterial cryptochromes and photolyases (FeS-BCP), a subtype of prokaryotic 6-4 photolyases. This photolyase was named SePhrB that could efficiently photorepair 6-4 photoproducts in DNA. Chemical analyses revealed that SePhrB contains a catalytic FAD cofactor and an iron-sulfur cluster. All of previously reported FeS-BCPs contain 6,7-dimethyl-8-ribityllumazine (DMRL) as their antenna chromophores. Here, we first demonstrated that SePhrB possesses 7,8-didemethyl-8-hydroxy-5-deazariboflavin (8-HDF) as an antenna chromophore. Nevertheless, SePhrB could be photoreduced without external electron donors. After being photoreduced, the reduced FAD cofactor in SePhrB was extremely stable against air oxidation. These results suggest that FeS-BCPs are more diverse than expected which deserve further investigation.

The aryl hydrocarbon receptor (AHR) regulates the expression of numerous genes in response to activation by agonists including xenobiotics. Although it is well appreciated that environmental signals and cell intrinsic features may modulate this transcriptional response, how it is mechanistically achieved remains poorly understood. We show that hexokinase 2 (HK2) a metabolic enzyme fuelling cancer cell growth, is a transcriptional target of AHR as well as a modulator of its activity. Expression of HK2 is positively regulated by AHR upon exposure to agonists both in human cells and in mice lung tissues. Conversely, over-expression of HK2 regulates the abundance of many proteins involved in the regulation of AHR signalling and these changes are linked with altered AHR expression levels and transcriptional activity. HK2 expression also shows a negative correlation with AHR promoter methylation in tumours, and these tumours with high HK2 expression and low AHR methylation are associated with a worse overall survival in patients. In sum, our study provides novel insights into how AHR signalling is regulated which may help our understanding of the context-specific effects of this pathway and may have implications in cancer.

During fasting, hepatocytes produce glucose in response to hormonal signals. Glucagon and glucocorticoids are principal fasting hormones that cooperate in regulating glucose production via gluconeogenesis. However, how these hormone signals are integrated and interpreted to a biological output is unknown. Here, we use genome-wide profiling of gene expression, enhancer dynamics and transcription factor (TF) binding in primary mouse hepatocytes to uncover the mode of cooperation between glucagon and glucocorticoids. We found that compared to a single treatment with each hormone, a dual treatment directs hepatocytes to a pro-gluconeogenic gene program by synergistically inducing gluconeogenic genes. The cooperative mechanism driving synergistic gene expression is based on 'assisted loading' whereby a glucagon-activated TF (cAMP responsive element binding protein; CREB) leads to enhancer activation which facilitates binding of the glucocorticoid receptor (GR) upon glucocorticoid stimulation. Glucagon does not only activate single enhancers but also activates enhancer clusters, thereby assisting the loading of GR also across enhancer units within the cluster. In summary, we show that cells integrate extracellular signals by an enhancer-specific mechanism: one hormone-activated TF activates enhancers, thereby assisting the loading of a TF stimulated by a second hormone, leading to synergistic gene induction and a tailored transcriptional response to fasting.

Auxilliary splicing sequences in exons, known as enhancers (ESEs) and silencers (ESSs), have been subject to strong selection pressures at the RNA and protein level. The protein component of this splicing code is substantial, recently estimated at ∼50% of the total information within ESEs, but remains poorly understood. The ESE/ESS profiles were previously associated with the Irving-Williams (I-W) stability series for divalent metals, suggesting that the ESE/ESS evolution was shaped by metal binding sites. Here, we have examined splicing activities of exonic sequences that encode protein binding sites for Ca2+, a weak binder in the I-W affinity order. We found that predicted exon inclusion levels for the EF-hand motifs and for Ca2+-binding residues in nonEF-hand proteins were higher than for average exons. For canonical EF-hands, the increase was centred on the EF-hand chelation loop and, in particular, on Ca2+-coordinating residues, with a 1>12>3∼5>9 hierarchy in the 12-codon loop consensus and usage bias at codons 1 and 12. The same hierarchy but a lower increase was observed for noncanonical EF-hands, except for S100 proteins. EF-hand loops preferentially accumulated exon splits in two clusters, one located in their N-terminal halves and the other around codon 12. Using splicing assays and published crosslinking and immunoprecipitation data, we identify candidate trans-acting factors that preferentially bind conserved GA-rich motifs encoding negatively charged amino acids in the loops. Together, these data provide evidence for the high capacity of codons for Ca2+-coordinating residues to be retained in mature transcripts, facilitating their exon-level expansion during eukaryotic evolution.

Despite the rapid and broad implementation of CRISPR-Cas9-based technologies, convenient tools to modulate dose, timing, and precision remain limited. Building on methods using synthetic peptide nucleic acids (PNAs) to bind RNA with unusually high affinity, we describe guide RNA (gRNA) spacer-targeted, or 'antispacer', PNAs as a tool to modulate Cas9 binding and activity in cells in a sequence-specific manner. We demonstrate that PNAs rapidly and efficiently target complexed gRNA spacer sequences at low doses and without design restriction for sequence-selective Cas9 inhibition. We further show that short PAM-proximal antispacer PNAs achieve potent cleavage inhibition (over 2000-fold reduction) and that PAM-distal PNAs modify gRNA affinity to promote on-target specificity. Finally, we apply antispacer PNAs for temporal regulation of two dCas9-fusion systems. These results present a novel rational approach to nucleoprotein engineering and describe a rapidly implementable antisense platform for CRISPR-Cas9 modulation to improve spatiotemporal versatility and safety across applications.

Deciphering the cellular composition in genome-wide spatially resolved transcriptomic data is a critical task to clarify the spatial context of cells in a tissue. In this study, we developed a method, CellDART, which estimates the spatial distribution of cells defined by single-cell level data using domain adaptation of neural networks and applied it to the spatial mapping of human lung tissue. The neural network that predicts the cell proportion in a pseudospot, a virtual mixture of cells from single-cell data, is translated to decompose the cell types in each spatial barcoded region. First, CellDART was applied to a mouse brain and a human dorsolateral prefrontal cortex tissue to identify cell types with a layer-specific spatial distribution. Overall, the proposed approach showed more stable and higher accuracy with short execution time compared to other computational methods to predict the spatial location of excitatory neurons. CellDART was capable of decomposing cellular proportion in mouse hippocampus Slide-seq data. Furthermore, CellDART elucidated the cell type predominance defined by the human lung cell atlas across the lung tissue compartments and it corresponded to the known prevalent cell types. CellDART is expected to help to elucidate the spatial heterogeneity of cells and their close interactions in various tissues.

Spatiotemporal gene expression patterns are governed to a large extent by the activity of enhancer elements, which engage in physical contacts with their target genes. Identification of enhancer-promoter (EP) links that are functional only in a specific subset of cell types is a key challenge in understanding gene regulation. We introduce CT-FOCS (cell type FOCS), a statistical inference method that uses linear mixed effect models to infer EP links that show marked activity only in a single or a small subset of cell types out of a large panel of probed cell types. Analyzing 808 samples from FANTOM5, covering 472 cell lines, primary cells and tissues, CT-FOCS inferred such EP links more accurately than recent state-of-the-art methods. Furthermore, we show that strictly cell type-specific EP links are very uncommon in the human genome.

Type II toxin-antitoxin (TA) systems are two-gene modules widely distributed among prokaryotes. GNAT toxins associated with the DUF1778 antitoxins represent a large family of type II TAs. GNAT toxins inhibit cell growth by disrupting translation via acetylation of aminoacyl-tRNAs. In this work, we explored the evolutionary trajectory of GNAT toxins. Using LC/MS detection of acetylated aminoacyl-tRNAs combined with ribosome profiling, we systematically investigated the in vivo substrate specificity of an array of diverse GNAT toxins. Our functional data show that the majority of GNAT toxins are specific to Gly-tRNA isoacceptors. However, the phylogenetic analysis shows that the ancestor of GNAT toxins was likely a relaxed specificity enzyme capable of acetylating multiple elongator tRNAs. Together, our data provide a remarkable snapshot of the evolution of substrate specificity.

Recent advances in machine learning and natural language processing have made it possible to profoundly advance our ability to accurately predict protein structures and their functions. While such improvements are significantly impacting the fields of biology and biotechnology at large, such methods have the downside of high demands in terms of computing power and runtime, hampering their applicability to large datasets. Here, we present NetSurfP-3.0, a tool for predicting solvent accessibility, secondary structure, structural disorder and backbone dihedral angles for each residue of an amino acid sequence. This NetSurfP update exploits recent advances in pre-trained protein language models to drastically improve the runtime of its predecessor by two orders of magnitude, while displaying similar prediction performance. We assessed the accuracy of NetSurfP-3.0 on several independent test datasets and found it to consistently produce state-of-the-art predictions for each of its output features, with a runtime that is up to to 600 times faster than the most commonly available methods performing the same tasks. The tool is freely available as a web server with a user-friendly interface to navigate the results, as well as a standalone downloadable package.