Structure-functional characterization of Lactococcus AbiA phage defense system.

Bacterial reverse transcriptases (RTs) are a large and diverse enzyme family. AbiA, AbiK and Abi-P2 are abortive infection system (Abi) RTs that mediate defense against bacteriophages. What sets Abi RTs apart from other RT enzymes is their ability to synthesize long DNA products of random sequences in a template- and primer-independent manner. Structures of AbiK and Abi-P2 representatives have recently been determined, but there are no structural data available for AbiA. Here, we report the crystal structure of Lactococcus AbiA polymerase in complex with a single-stranded polymerization product. AbiA comprises three domains: an RT-like domain, a helical domain that is typical for Abi polymerases, and a higher eukaryotes and prokaryotes nucleotide-binding (HEPN) domain that is common for many antiviral proteins. AbiA forms a dimer that distinguishes it from AbiK and Abi-P2, which form trimers/hexamers. We show the DNA polymerase activity of AbiA in an in vitro assay and demonstrate that it requires the presence of the HEPN domain which is enzymatically inactive. We validate our biochemical and structural results in vivo through bacteriophage infection assays. Finally, our in vivo results suggest that AbiA-mediated phage defense may not rely on AbiA-mediated cell death.

The structural basis of hyperpromiscuity in a core combinatorial network of type II toxin-antitoxin and related phage defense systems.

Toxin-antitoxin (TA) systems are a large group of small genetic modules found in prokaryotes and their mobile genetic elements. Type II TAs are encoded as bicistronic (two-gene) operons that encode two proteins: a toxin and a neutralizing antitoxin. Using our tool NetFlax (standing for Network-FlaGs for toxins and antitoxins), we have performed a large-scale bioinformatic analysis of proteinaceous TAs, revealing interconnected clusters constituting a core network of TA-like gene pairs. To understand the structural basis of toxin neutralization by antitoxins, we have predicted the structures of 3,419 complexes with AlphaFold2. Together with mutagenesis and functional assays, our structural predictions provide insights into the neutralizing mechanism of the hyperpromiscuous Panacea antitoxin domain. In antitoxins composed of standalone Panacea, the domain mediates direct toxin neutralization, while in multidomain antitoxins the neutralization is mediated by other domains, such as PAD1, Phd-C, and ZFD. We hypothesize that Panacea acts as a sensor that regulates TA activation. We have experimentally validated 16 NetFlax TA systems and used domain annotations and metabolic labeling assays to predict their potential mechanisms of toxicity (such as membrane disruption, and inhibition of cell division or protein synthesis) as well as biological functions (such as antiphage defense). We have validated the antiphage activity of a RosmerTA system encoded by Gordonia phage Kita, and used fluorescence microscopy to confirm its predicted membrane-depolarizing activity. The interactive version of the NetFlax TA network that includes structural predictions can be accessed at http://netflax.webflags.se/.

Genome-encoded ABCF factors implicated in intrinsic antibiotic resistance in Gram-positive bacteria: VmlR2, Ard1 and CplR.

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

uORF4u: a tool for annotation of conserved upstream open reading frames.

SUMMARY: Upstream open reading frames (uORFs, often encoding so-called leader peptides) can regulate translation and transcription of downstream main ORFs (mORFs) in prokaryotes and eukaryotes. However, annotation of novel functional uORFs is challenging due to their short size of usually <100 codons. While transcription- and translation-level next-generation sequencing methods can be used for genome-wide functional uORF identification, this data are not available for the vast majority of species with sequenced genomes. At the same time, the exponentially increasing amount of genome assemblies gives us the opportunity to take advantage of evolutionary conservation in our predictions of functional ORFs.Here, we present a tool for conserved uORF annotation in 5' upstream sequences of a user-defined protein of interest or a set of protein homologs. It can also be used to find small conserved ORFs within a set of nucleotide sequences. The output includes publication-quality figures with multiple sequence alignments, sequence logos, and locus annotation of the predicted conserved uORFs in graphical vector format. AVAILABILITY AND IMPLEMENTATION: uORF4u is written in Python3 and runs on Linux and MacOS. The command-line interface covers most practical use cases, while the provided Python API allows usage within a Python program and additional customization. Source code is available from the GitHub page: github.com/GCA-VH-lab/uorf4u. Detailed documentation that includes an example-driven guide available at the software home page: gca-vh-lab.github.io/uorf4u. A web version of uORF4u is available at server.atkinson-lab.com/uorf4u.

A standard knockout procedure alters expression of adjacent loci at the translational level.

The Saccharomyces cerevisiae gene deletion collection is widely used for functional gene annotation and genetic interaction analyses. However, the standard G418-resistance cassette used to produce knockout mutants delivers strong regulatory elements into the target genetic loci. To date, its side effects on the expression of neighboring genes have never been systematically assessed. Here, using ribosome profiling data, RT-qPCR, and reporter expression, we investigated perturbations induced by the KanMX module. Our analysis revealed significant alterations in the transcription efficiency of neighboring genes and, more importantly, severe impairment of their mRNA translation, leading to changes in protein abundance. In the 'head-to-head' orientation of the deleted and neighboring genes, knockout often led to a shift of the transcription start site of the latter, introducing new uAUG codon(s) into the expanded 5' untranslated region (5' UTR). In the 'tail-to-tail' arrangement, knockout led to activation of alternative polyadenylation signals in the neighboring gene, thus altering its 3' UTR. These events may explain the so-called neighboring gene effect (NGE), i.e. false genetic interactions of the deleted genes. We estimate that in as much as ∼1/5 of knockout strains the expression of neighboring genes may be substantially (>2-fold) deregulated at the level of translation.

Human Tissues Exhibit Diverse Composition of Translation Machinery.

While protein synthesis is vital for the majority of cell types of the human body, diversely differentiated cells require specific translation regulation. This suggests the specialization of translation machinery across tissues and organs. Using transcriptomic data from GTEx, FANTOM, and Gene Atlas, we systematically explored the abundance of transcripts encoding translation factors and aminoacyl-tRNA synthetases (ARSases) in human tissues. We revised a few known and identified several novel translation-related genes exhibiting strict tissue-specific expression. The proteins they encode include eEF1A1, eEF1A2, PABPC1L, PABPC3, eIF1B, eIF4E1B, eIF4ENIF1, and eIF5AL1. Furthermore, our analysis revealed a pervasive tissue-specific relative abundance of translation machinery components (e.g., PABP and eRF3 paralogs, eIF2B and eIF3 subunits, eIF5MPs, and some ARSases), suggesting presumptive variance in the composition of translation initiation, elongation, and termination complexes. These conclusions were largely confirmed by the analysis of proteomic data. Finally, we paid attention to sexual dimorphism in the repertoire of translation factors encoded in sex chromosomes (eIF1A, eIF2γ, and DDX3), and identified the testis and brain as organs with the most diverged expression of translation-associated genes.

svist4get: a simple visualization tool for genomic tracks from sequencing experiments.

BACKGROUND: High-throughput sequencing often provides a foundation for experimental analyses in the life sciences. For many such methods, an intermediate layer of bioinformatics data analysis is the genomic signal track constructed by short read mapping to a particular genome assembly. There are many software tools to visualize genomic tracks in a web browser or with a stand-alone graphical user interface. However, there are only few command-line applications suitable for automated usage or production of publication-ready visualizations. RESULTS: Here we present svist4get, a command-line tool for customizable generation of publication-quality figures based on data from genomic signal tracks. Similarly to generic genome browser software, svist4get visualizes signal tracks at a given genomic location and is able to aggregate data from several tracks on a single plot along with the transcriptome annotation. The resulting plots can be saved as the vector or high-resolution bitmap images. We demonstrate practical use cases of svist4get for Ribo-Seq and RNA-Seq data. CONCLUSIONS: svist4get is implemented in Python 3 and runs on Linux. The command-line interface of svist4get allows for easy integration into bioinformatics pipelines in a console environment. Extra customization is possible through configuration files and Python API. For convenience, svist4get is provided as pypi package. The source code is available at https://bitbucket.org/artegorov/svist4get/.

The ABCF ATPase New1 resolves translation termination defects associated with specific tRNAArg and tRNALys isoacceptors in the P site.

The efficiency of translation termination is determined by the nature of the stop codon as well as its context. In eukaryotes, recognition of the A-site stop codon and release of the polypeptide are mediated by release factors eRF1 and eRF3, respectively. Translation termination is modulated by other factors which either directly interact with release factors or bind to the E-site and modulate the activity of the peptidyl transferase center. Previous studies suggested that the Saccharomyces cerevisiae ABCF ATPase New1 is involved in translation termination and/or ribosome recycling, however, the exact function remained unclear. Here, we have applied 5PSeq, single-particle cryo-EM and readthrough reporter assays to provide insight into the biological function of New1. We show that the lack of New1 results in ribosomal stalling at stop codons preceded by a lysine or arginine codon and that the stalling is not defined by the nature of the C-terminal amino acid but rather by the identity of the tRNA isoacceptor in the P-site. Collectively, our results suggest that translation termination is inefficient when ribosomes have specific tRNA isoacceptors in the P-site and that the recruitment of New1 rescues ribosomes at these problematic termination contexts.

Mechanism of phage sensing and restriction by toxin-antitoxin-chaperone systems.

Toxin-antitoxins (TAs) are prokaryotic two-gene systems composed of a toxin neutralized by an antitoxin. Toxin-antitoxin-chaperone (TAC) systems additionally include a SecB-like chaperone that stabilizes the antitoxin by recognizing its chaperone addiction (ChAD) element. TACs mediate antiphage defense, but the mechanisms of viral sensing and restriction are unexplored. We identify two Escherichia coli antiphage TAC systems containing host inhibition of growth (HigBA) and CmdTA TA modules, HigBAC and CmdTAC. HigBAC is triggered through recognition of the gpV major tail protein of phage λ. Chaperone HigC recognizes gpV and ChAD via analogous aromatic molecular patterns, with gpV outcompeting ChAD to trigger toxicity. For CmdTAC, the CmdT ADP-ribosyltransferase toxin modifies mRNA to halt protein synthesis and limit phage propagation. Finally, we establish the modularity of TACs by creating a hybrid broad-spectrum antiphage system combining the CmdTA TA warhead with a HigC chaperone phage sensor. Collectively, these findings reveal the potential of TAC systems in broad-spectrum antiphage defense.

Assessing Ribosome Distribution Along Transcripts with Polarity Scores and Regression Slope Estimates.

During translation, the rate of ribosome movement along mRNA varies. This leads to a non-uniform ribosome distribution along the transcript, depending on local mRNA sequence, structure, tRNA availability, and translation factor abundance, as well as the relationship between the overall rates of initiation, elongation, and termination. Stress, antibiotics, and genetic perturbations affecting composition and properties of translation machinery can alter the ribosome positional distribution dramatically. Here, we offer a computational protocol for analyzing positional distribution profiles using ribosome profiling (Ribo-Seq) data. The protocol uses papolarity, a new Python toolkit for the analysis of transcript-level short read coverage profiles. For a single sample, for each transcript papolarity allows for computing the classic polarity metric which, in the case of Ribo-Seq, reflects ribosome positional preferences. For comparison versus a control sample, papolarity estimates an improved metric, the relative linear regression slope of coverage along transcript length. This involves de-noising by profile segmentation with a Poisson model and aggregation of Ribo-Seq coverage within segments, thus achieving reliable estimates of the regression slope. The papolarity software and the associated protocol can be conveniently used for Ribo-Seq data analysis in the command-line Linux environment. Papolarity package is available through Python pip package manager. The source code is available at https://github.com/autosome-ru/papolarity .

Ribo-Seq and RNA-Seq of TMA46 ( DFRP1) and GIR2 ( DFRP2) knockout yeast strains.

In eukaryotes, stalled and collided ribosomes are recognized by several conserved multicomponent systems, which either block protein synthesis in situ and resolve the collision locally, or trigger a general stress response. Yeast ribosome-binding GTPases RBG1 (DRG1 in mammals) and RBG2 (DRG2) form two distinct heterodimers with TMA46 (DFRP1) and GIR2 (DFRP2), respectively, both involved in mRNA translation. Accumulated evidence suggests that the dimers play partially redundant roles in elongation processivity and resolution of ribosome stalling and collision events, as well as in the regulation of GCN1-mediated signaling involved in ribosome-associated quality control (RQC). They also genetically interact with SLH1 (ASCC3) helicase, a key component of RQC trigger (RQT) complex disassembling collided ribosomes. Here, we present RNA-Seq and ribosome profiling (Ribo-Seq) data from S. cerevisiae strains with individual deletions of the TMA46 and GIR2 genes. Raw RNA-Seq and Ribo-Seq data as well as gene-level read counts are available in NCBI Gene Expression Omnibus (GEO) repository under GEO accession GSE185458 and GSE185286.

Translatome and transcriptome analysis of TMA20 (MCT-1) and TMA64 (eIF2D) knockout yeast strains.

TMA20 (MCT-1), TMA22 (DENR) and TMA64 (eIF2D) are eukaryotic translation factors involved in ribosome recycling and re-initiation. They operate with P-site bound tRNA in post-termination or (re-)initiation translation complexes, thus participating in the removal of 40S ribosomal subunit from mRNA stop codons after termination and controlling translation re-initiation on mRNAs with upstream open reading frames (uORFs), as well as de novo initiation on some specific mRNAs. Here we report ribosomal profiling data of S.cerevisiae strains with individual deletions of TMA20, TMA64 or both TMA20 and TMA64 genes. We provide RNA-Seq and Ribo-Seq data from yeast strains grown in the rich YPD or minimal SD medium. We illustrate our data by plotting differential distribution of ribosomal-bound mRNA fragments throughout uORFs in 5'-untranslated region (5' UTR) of GCN4 mRNA and on mRNA transcripts encoded in MAT locus in the mutant and wild-type strains, thus providing a basis for investigation of the role of these factors in the stress response, mating and sporulation. We also document a shift of transcription start site of the APC4 gene which occurs when the neighboring TMA64 gene is replaced by the standard G418-resistance cassette used for the creation of the Yeast Deletion Library. This shift results in dramatic deregulation of the APC4 gene expression, as revealed by our Ribo-Seq data, which can be probably used to explain strong genetic interactions of TMA64 with genes involved in the cell cycle and mitotic checkpoints. Raw RNA-Seq and Ribo-Seq data as well as all gene counts are available in NCBI Gene Expression Omnibus (GEO) repository under GEO accession GSE122039 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE122039).