Ultraconserved bacteriophage genome sequence identified in 1300-year-old human palaeofaeces.

Bacteriophages are widely recognised as rapidly evolving biological entities. However, knowledge about ancient bacteriophages is limited. Here, we analyse DNA sequence datasets previously generated from ancient palaeofaeces and human gut-content samples, and identify an ancient phage genome nearly identical to present-day Mushuvirus mushu, a virus that infects gut commensal bacteria. The DNA damage patterns of the genome are consistent with its ancient origin and, despite 1300 years of evolution, the ancient Mushuvirus genome shares 97.7% nucleotide identity with its modern counterpart, indicating a long-term relationship between the prophage and its host. In addition, we reconstruct and authenticate 297 other phage genomes from the last 5300 years, including those belonging to unknown families. Our findings demonstrate the feasibility of reconstructing ancient phage genome sequences, thus expanding the known virosphere and offering insights into phage-bacteria interactions spanning several millennia.

A ligand-receptor interactome atlas of the zebrafish.

Studies in zebrafish can unravel the functions of cellular communication and thus identify novel bench-to-bedside drugs targeting cellular communication signaling molecules. Due to the incomplete annotation of zebrafish proteome, the knowledge of zebrafish receptors, ligands, and tools to explore their interactome is limited. To address this gap, we de novo predicted the cellular localization of zebrafish reference proteome using deep learning algorithm. We combined the predicted and existing annotations on cellular localization of zebrafish proteins and created repositories of zebrafish ligands, membrane receptome, and interactome as well as associated diseases and targeting drugs. Unlike other tools, our interactome atlas is based on both the physical interaction data of zebrafish proteome and existing human ligand-receptor pair databases. The resources are available as R and Python scripts. DanioTalk provides a novel resource for researchers interested in targeting cellular communication in zebrafish, as we demonstrate in applications studying synapse and axo-glial interactome. DanioTalk methodology can be applied to build and explore the ligand-receptor atlas of other non-mammalian model organisms.

TRGdb: a universal resource for the exploration of taxonomically restricted genes in bacteria.

The TRGdb database is a resource dedicated to taxonomically restricted genes (TRGs) in bacteria. It provides a comprehensive collection of genes that are specific to different genera and species, according to the latest release of bacterial taxonomy. The user interface allows for easy browsing and searching as well as sequence similarity exploration. The website also provides information on each TRG protein sequence, including its level of disorder, complexity and tendency to aggregate. TRGdb is a valuable resource for gaining a deeper understanding of the TRG-associated, unique features, and characteristics of bacterial organisms. Database URL www.combio.pl/trgdb.

The effects of nature-inspired amino acid substitutions on structural and biochemical properties of the E. coli L-asparaginase EcAIII.

The Escherichia coli enzyme EcAIII catalyzes the hydrolysis of L-Asn to L-Asp and ammonia. Using a nature-inspired mutagenesis approach, we designed and produced five new EcAIII variants (M200I, M200L, M200K, M200T, M200W). The modified proteins were characterized by spectroscopic and crystallographic methods. All new variants were enzymatically active, confirming that the applied mutagenesis procedure has been successful. The determined crystal structures revealed new conformational states of the EcAIII molecule carrying the M200W mutation and allowed a high-resolution observation of an acyl-enzyme intermediate with the M200L mutant. In addition, we performed structure prediction, substrate docking, and molecular dynamics simulations for 25 selected bacterial orthologs of EcAIII, to gain insights into how mutations at the M200 residue affect the active site and substrate binding mode. This comprehensive strategy, including both experimental and computational methods, can be used to guide further enzyme engineering and can be applied to the study of other proteins of medicinal or biotechnological importance.

Taxonomically Restricted Genes in Bacillus may Form Clusters of Homologs and Can be Traced to a Large Reservoir of Noncoding Sequences.

Taxonomically restricted genes (TRGs) are unique for a defined group of organisms and may act as potential genetic determinants of lineage-specific, biological properties. Here, we explore the TRGs of highly diverse and economically important Bacillus bacteria by examining commonly used TRG identification parameters and data sources. We show the significant effects of sequence similarity thresholds, composition, and the size of the reference database in the identification process. Subsequently, we applied stringent TRG search parameters and expanded the identification procedure by incorporating an analysis of noncoding and non-syntenic regions of non-Bacillus genomes. A multiplex annotation procedure minimized the number of false-positive TRG predictions and showed nearly one-third of the alleged TRGs could be mapped to genes missed in genome annotations. We traced the putative origin of TRGs by identifying homologous, noncoding genomic regions in non-Bacillus species and detected sequence changes that could transform these regions into protein-coding genes. In addition, our analysis indicated that Bacillus TRGs represent a specific group of genes mostly showing intermediate sequence properties between genes that are conserved across multiple taxa and nonannotated peptides encoded by open reading frames.

Massive annotation of bacterial L-asparaginases reveals their puzzling distribution and frequent gene transfer events.

L-Asparaginases, which convert L-asparagine to L-aspartate and ammonia, come in five types, AI-AV. Some bacterial type AII enzymes are a key element in the treatment of acute lymphoblastic leukemia in children, but new L-asparaginases with better therapeutic properties are urgently needed. Here, we search publicly available bacterial genomes to annotate L-asparaginase proteins belonging to the five known types. We characterize taxonomic, phylogenetic, and genomic patterns of L-asparaginase occurrences pointing to frequent horizontal gene transfer (HGT) events, also occurring multiple times in the same recipient species. We show that the reference AV gene, encoding a protein originally found and structurally studied in Rhizobium etli, was acquired via HGT from Burkholderia. We also describe the sequence variability of the five L-asparaginase types and map the conservation levels on the experimental or predicted structures of the reference enzymes, finding the most conserved residues in the protein core near the active site, and the most variable ones on the protein surface. Additionally, we highlight the most common sequence features of bacterial AII proteins that may aid in selecting therapeutic L-asparaginases. Finally, we point to taxonomic units of bacteria that do not contain recognizable sequences of any of the known L-asparaginase types, implying that those microorganisms most likely contain new, as yet unknown types of L-asparaginases. Such novel enzymes, when properly identified and characterized, could hold promise as antileukemic drugs.

Daily Reports on Phage-Host Interactions.

Understanding phage-host relationships is crucial for the study of virus biology and the application of phages in biotechnology and medicine. However, information concerning the range of hosts for bacterial and archaeal viruses is scattered across numerous databases and is difficult to obtain. Therefore, here we present PHD (Phage & Host Daily), a web application that offers a comprehensive, up-to-date catalog of known phage-host associations that allows users to select viruses targeting specific bacterial and archaeal taxa of interest. Our service combines the latest information on virus-host interactions from seven source databases with current taxonomic classification retrieved directly from the groups and institutions responsible for its maintenance. The web application also provides summary statistics on host and virus diversity, their pairwise interactions, and the host range of deposited phages. PHD is updated daily and available at http://phdaily.info or http://combio.pl/phdaily.

Biogenesis, conservation, and function of miRNA in liverworts.

MicroRNAs (miRNAs) are small non-coding endogenous RNA molecules, 18-24 nucleotides long, that control multiple gene regulatory pathways via post-transcriptional gene silencing in eukaryotes. To develop a comprehensive picture of the evolutionary history of miRNA biogenesis and action in land plants, studies on bryophyte representatives are needed. Here, we review current understanding of liverwort MIR gene structure, miRNA biogenesis, and function, focusing on the simple thalloid Pellia endiviifolia and the complex thalloid Marchantia polymorpha. We review what is known about conserved and non-conserved miRNAs, their targets, and the functional implications of miRNA action in M. polymorpha and P. endiviifolia. We note that most M. polymorpha miRNAs are encoded within protein-coding genes and provide data for 23 MIR gene structures recognized as independent transcriptional units. We identify M. polymorpha genes involved in miRNA biogenesis that are homologous to those identified in higher plants, including those encoding core microprocessor components and other auxiliary and regulatory proteins that influence the stability, folding, and processing of pri-miRNAs. We analyzed miRNA biogenesis proteins and found similar domain architecture in most cases. Our data support the hypothesis that almost all miRNA biogenesis factors in higher plants are also present in liverworts, suggesting that they emerged early during land plant evolution.

PHIST: fast and accurate prediction of prokaryotic hosts from metagenomic viral sequences.

SUMMARY: Phage-Host Interaction Search Tool (PHIST) predicts prokaryotic hosts of viruses based on exact matches between viral and host genomes. It improves host prediction accuracy at species level over current alignment-based tools (on average by 3 percentage points) as well as alignment-free and CRISPR-based tools (by 14-20 percentage points). PHIST is also two orders of magnitude faster than alignment-based tools making it suitable for metagenomics studies. AVAILABILITY AND IMPLEMENTATION: GNU-licensed C++ code wrapped in Python API available at: https://github.com/refresh-bio/phist. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Taxonomy-aware, sequence similarity ranking reliably predicts phage-host relationships.

BACKGROUND: Characterizing phage-host interactions is critical to understanding the ecological role of both partners and effective isolation of phage therapeuticals. Unfortunately, experimental methods for studying these interactions are markedly slow, low-throughput, and unsuitable for phages or hosts difficult to maintain in laboratory conditions. Therefore, a number of in silico methods emerged to predict prokaryotic hosts based on viral sequences. One of the leading approaches is the application of the BLAST tool that searches for local similarities between viral and microbial genomes. However, this prediction method has three major limitations: (i) top-scoring sequences do not always point to the actual host; (ii) mosaic virus genomes may match to many, typically related, bacteria; and (iii) viral and host sequences may diverge beyond the point where their relationship can be detected by a BLAST alignment. RESULTS: We created an extension to BLAST, named Phirbo, that improves host prediction quality beyond what is obtainable from standard BLAST searches. The tool harnesses information concerning sequence similarity and bacteria relatedness to predict phage-host interactions. Phirbo was evaluated on three benchmark sets of known virus-host pairs, and it improved precision and recall by 11-40 percentage points over currently available, state-of-the-art, alignment-based, alignment-free, and machine-learning host prediction tools. Moreover, the discriminatory power of Phirbo for the recognition of virus-host relationships surpassed the results of other tools by at least 10 percentage points (area under the curve = 0.95), yielding a mean host prediction accuracy of 57% and 68% at the genus and family levels, respectively, and drops by 12 percentage points when using only a fraction of viral genome sequences (3 kb). Finally, we provide insights into a repertoire of protein and ncRNA genes that are shared between phages and hosts and may be prone to horizontal transfer during infection. CONCLUSIONS: Our results suggest that Phirbo is a simple and effective tool for predicting phage-host relationships.

RaFAH: Host prediction for viruses of Bacteria and Archaea based on protein content.

Culture-independent approaches have recently shed light on the genomic diversity of viruses of prokaryotes. One fundamental question when trying to understand their ecological roles is: which host do they infect? To tackle this issue we developed a machine-learning approach named Random Forest Assignment of Hosts (RaFAH), that uses scores to 43,644 protein clusters to assign hosts to complete or fragmented genomes of viruses of Archaea and Bacteria. RaFAH displayed performance comparable with that of other methods for virus-host prediction in three different benchmarks encompassing viruses from RefSeq, single amplified genomes, and metagenomes. RaFAH was applied to assembled metagenomic datasets of uncultured viruses from eight different biomes of medical, biotechnological, and environmental relevance. Our analyses led to the identification of 537 sequences of archaeal viruses representing unknown lineages, whose genomes encode novel auxiliary metabolic genes, shedding light on how these viruses interfere with the host molecular machinery. RaFAH is available at https://sourceforge.net/projects/rafah/.

Quantitative Analysis of Plant miRNA Primary Transcripts.

MicroRNAs control plant development and are key regulators of plant responses to biotic and abiotic stresses. Thus, their expression must be carefully controlled since both excess and deficiency of a given microRNA may be deleterious to plant cell. MicroRNA expression regulation can occur at several stages of their biogenesis pathway. One of the most important of these regulatory checkpoints is transcription efficiency. mirEX database is a tool for exploration and visualization of plant pri-miRNA expression profiles. It includes results obtained using high-throughput RT-qPCR platform designed to monitor pri-miRNA expression in different miRNA biogenesis mutants and developmental stages of Arabidopsis, barley, and Pellia plants. A step-by-step instruction for browsing the database and detailed protocol for high-throughput RT-qPCR experiments, including list of primers designed for the amplification of pri-miRNAs, are presented.

TOXiTAXi: a web resource for toxicity of Bacillus thuringiensis protein compositions towards species of various taxonomic groups.

Bioinsecticides consisting of different sets of Bacillus thuringiensis (Bt) Cry, Cyt and Vip toxins are broadly used in pest control. Possible interactions (synergistic, additive or antagonistic) between these proteins can not only influence the overall efficacy of certain Bt-based bioinsecticide, but also raise questions regarding environmental safety. Here, we assemble, summarize and analyze the outcomes of experiments published over 30 years, investigating combinatorial effects among Bt Cry, Cyt and Vip toxins. We collected the results on 118 various two-to-five-component combinations that have been bioassayed against 38 invertebrate species. Synergism, additive effect and antagonism was indicated in 54%, 32% and 14% of experiments, respectively. Synergism was noted most frequently for Cry/Cyt combinations, followed by Cyt/Vip and Cry/Cry. In Cry/Vip combinations, antagonism is more frequent and higher in magnitude compared to other categories. Despite a significant number of tested Bt toxin combinations, most of them have been bioassayed only against one pest species. To aid the research on Bt pesticidal protein activity, we present TOXiTAXi ( http://www.combio.pl/toxitaxi/ ), a universal database and a dedicated web tool to conveniently gather and analyze the existing and future bioassay results on biocidal activity of toxins against various taxonomic groups.

Silencing of HvGSK1.1-A GSK3/SHAGGY-Like Kinase-Enhances Barley (Hordeum vulgare L.) Growth in Normal and in Salt Stress Conditions.

Glycogen synthase kinase 3 (GSK3) is a highly conserved kinase present in all eukaryotes and functions as a key regulator of a wide range of physiological and developmental processes. The kinase, known in land plants as GSK3/SHAGGY-like kinase (GSK), is a key player in the brassinosteroid (BR) signaling pathway. The GSK genes, through the BRs, affect diverse developmental processes and modulate responses to environmental factors. In this work, we describe functional analysis of HvGSK1.1, which is one of the GSK3/SHAGGY-like orthologs in barley. The RNAi-mediated silencing of the target HvGSK1.1 gene was associated with modified expression of its paralogs HvGSK1.2, HvGSK2.1, HvGSK3.1, and HvGSK4.1 in plants grown in normal and in salt stress conditions. Low nucleotide similarity between the silencing fragment and barley GSK genes and the presence of BR-dependent transcription factors' binding sites in promoter regions of barley and rice GSK genes imply an innate mechanism responsible for co-regulation of the genes. The results of the leaf inclination assay indicated that silencing of HvGSK1.1 and the changes of GSK paralogs enhanced the BR-dependent signaling in the plants. The strongest phenotype of transgenic lines with downregulated HvGSK1.1 and GSK paralogs had greater biomass of the seedlings grown in normal conditions and salt stress as well as elevated kernel weight of plants grown in normal conditions. Both traits showed a strong negative correlation with the transcript level of the target gene and the paralogs. The characteristics of barley lines with silenced expression of HvGSK1.1 are compatible with the expected phenotypes of plants with enhanced BR signaling. The results show that manipulation of the GSK-encoding genes provides data to explore their biological functions and confirm it as a feasible strategy to generate plants with improved agricultural traits.

TaWAK6 encoding wall-associated kinase is involved in wheat resistance to leaf rust similar to adult plant resistance.

In wheat, adult plant resistance (APR) to leaf rust (Puccinia triticina), is effective in restricting pathogen growth and provides durable resistance against a wide range of virulent forms of P. triticina. Despite the importance, there is limited knowledge on the molecular basis of this type of resistance. We isolated and characterized the wall-associated kinase encoding gene in wheat, and assigned it as TaWAK6. Localization of TaWAK6 homeologs in A and B wheat subgenomes was consistent with the presence of the gene's orthologs in T. urartu (AA) and T. dicoccoides (AABB) and with the absence of its orthologs in Aegilops tauschii (DD). Overexpression of TaWAK6 did not change the wheat phenotype, nor did it affect seedling resistance. However, the adult plants overexpressing TaWAK6 showed that important parameters of APR were significantly elevated. Infection types scored on the first (flag), second and third leaves indicated elevated resistance, which significantly correlated with expression of TaWAK6. Analysis of plant-pathogen interactions showed a lower number of uredinia and higher rates of necrosis at the infection sites and this was associated with smaller size of uredinia and a longer latent period. The results indicated a role of TaWAK6 in quantitative partial resistance similar to APR in wheat. It is proposed that TaWAK6, which is a non-arginine-aspartate (non-RD) kinase, represents a novel class of quantitative immune receptors in monocots.

Benchmarking of alignment-free sequence comparison methods.

BACKGROUND: Alignment-free (AF) sequence comparison is attracting persistent interest driven by data-intensive applications. Hence, many AF procedures have been proposed in recent years, but a lack of a clearly defined benchmarking consensus hampers their performance assessment. RESULTS: Here, we present a community resource (http://afproject.org) to establish standards for comparing alignment-free approaches across different areas of sequence-based research. We characterize 74 AF methods available in 24 software tools for five research applications, namely, protein sequence classification, gene tree inference, regulatory element detection, genome-based phylogenetic inference, and reconstruction of species trees under horizontal gene transfer and recombination events. CONCLUSION: The interactive web service allows researchers to explore the performance of alignment-free tools relevant to their data types and analytical goals. It also allows method developers to assess their own algorithms and compare them with current state-of-the-art tools, accelerating the development of new, more accurate AF solutions.

Annotation and profiling of barley GLYCOGEN SYNTHASE3/Shaggy-like genes indicated shift in organ-preferential expression.

GLYCOGEN SYNTHASE KINASE3/Shaggy-like kinases (GSKs) represent a highly conserved group of proteins found in all eukaryotes. In plants they are encoded by multigene families and integrate signaling of brassinosteroids, auxin and abscisic acid in wide range of physiological and developmental processes with a strong impact on plant responses to environmental and biotic factors. Based on comprehensively studied structures of 10 Arabidopsis thaliana GSK genes and encoded proteins we report identification and phylogenetic reconstruction of 7 transcriptionally active GSK genes in barley. We re-evaluated annotation of the GSK genes in the current barley genome (Hv_IBSC_PGSB_v2) and provided data that a single gene annotated in the previous barley genome ensemble should be retained in the current one. The novel structure of another GSK, predicted in Hv_IBSC_PGSB_v2 to encode both GSK and amine oxidase domains, was proposed and experimentally confirmed based on the syntenic region in Brachypodium distachyon. The genes were assigned to 4 groups based on their encoded amino acid sequences and protein kinase domains. The analysis confirmed high level of conservation of functional protein domains and motifs among plant GSKs and the identified barley orthologs. Each of the seven identified HvGSK genes was expressed indicating semi-constitutive regulation in all tested organs and developmental stages. Regulation patterns of GSKs from the indicated groups showed a shift in organ-preferential expression in A. thaliana and barley illustrating diversification of biological roles of individual HvGSKs in different plant species.

tRex: A Web Portal for Exploration of tRNA-Derived Fragments in Arabidopsis thaliana.

tRNA-derived fragments (tRFs) constitute a new class of short regulatory RNAs that are a product of nascent or mature tRNA processing. tRF sequences have been identified in all domains of life; however, most published research pertains to human, yeast and some bacterial organisms. Despite growing interest in plant tRFs and accumulating evidence of their function in plant development and stress responses, no public, web-based repository dedicated to these molecules is currently available. Here, we introduce tRex (http://combio.pl/trex)-the first comprehensive data-driven online resource specifically dedicated to tRFs in the model plant Arabidopsis thaliana. The portal is based on verified Arabidopsis tRNA annotation and includes in-house-generated and publicly available small RNA sequencing experiments from various tissues, ecotypes, genotypes and stress conditions. The provided web-based tools are designed in a user-friendly manner and allow for seamless exploration of the data that are presented in the form of dynamic tables and cumulative coverage profiles. The tRex database is connected to external genomic and citation resources, which makes it a one-stop solution for Arabidopsis tRF-related research.

Alignment-free sequence comparison: benefits, applications, and tools.

Alignment-free sequence analyses have been applied to problems ranging from whole-genome phylogeny to the classification of protein families, identification of horizontally transferred genes, and detection of recombined sequences. The strength of these methods makes them particularly useful for next-generation sequencing data processing and analysis. However, many researchers are unclear about how these methods work, how they compare to alignment-based methods, and what their potential is for use for their research. We address these questions and provide a guide to the currently available alignment-free sequence analysis tools.

Identification and Analysis of WG/GW ARGONAUTE-Binding Domains.

WG/GW domains recruit ARGONAUTE (AGO) proteins to distinct silencing effector complexes using combinations of just two amino acids: tryptophan (W) and glycine (G), forming a wide arsenal of highly simplified interaction surfaces. These unstructured domains exhibit very low sequence identity and excessive length polymorphism, which makes identification of new AGO-binding proteins a challenging task as they escape detection with standard sequence comparison-based methods (e.g., BLAST, HMMER).In this chapter, we explain the use of tools for prediction of AGO-binding WG/GW domains in protein sequences. We also show how to computationally explore an up-to-date information about AGO-interacting proteins and discover new properties of WG/GW domains. Finally, we encourage readers to explore the game-like web application for in silico designing/modifying AGO-binding sequences as well as modeling mutagenesis experiments and predicting their potential effect on AGO-binding activity.