Global In-Silico Scenario of tRNA Genes and Their Organization in Virus Genomes

Virophages Found in Viromes from Lake Baikal

In this study, a previously little-studied group of viruses-virophages-was searched for and identified in the viromes of the ancient oligotrophic Lake Baikal. Virophages are small dsDNA viruses that parasitize giant viruses (e.g., Mimiviridae), which in turn affect unicellular eukaryotes. We analyzed eight viromes obtained from the deep-water areas of three basins of Lake Baikal and the shallow-water strait Maloye More in different seasons. The sequences of virophages were revealed in all viromes and were dominant after bacteriophages and algal viruses. Sixteen putative complete genomes of virophages were assembled, all of which contained four conserved genes encoding major capsid protein (MCP), minor capsid protein (mCP), maturation cysteine protease (PRO), and FtsK-HerA family DNA-packaging ATPase (ATPase). The MCP-based cluster analysis showed a sequence separation according to seasons, and a dependence on the geographical localization was not detected.

Simultaneous entry as an adaptation to virulence in a novel satellite-helper system infecting Streptomyces species

Satellites are mobile genetic elements that are dependent upon the replication machinery of their helper viruses. Bacteriophages have provided many examples of satellite nucleic acids that utilize their helper morphogenic genes for propagation. Here we describe two novel satellite-helper phage systems, Mulch and Flayer, that infect Streptomyces species. The satellites in these systems encode for encapsidation machinery but have an absence of key replication genes, thus providing the first example of bacteriophage satellite viruses. We also show that codon usage of the satellites matches the tRNA gene content of the helpers. The satellite in one of these systems, Flayer, does not appear to integrate into the host genome, which represents the first example of a virulent satellite phage. The Flayer satellite has a unique tail adaptation that allows it to attach to its helper for simultaneous co-infection. These findings demonstrate an ever-increasing array of satellite strategies for genetic dependence on their helpers in the evolutionary arms race between satellite and helper phages.

Genomic Analysis and Taxonomic Characterization of Seven Bacteriophage Genomes Metagenomic-Assembled from the Dishui Lake

Viruses in aquatic ecosystems exhibit remarkable abundance and diversity. However, scattered studies have been conducted to mine uncultured viruses and identify them taxonomically in lake water. Here, whole genomes (29-173 kbp) of seven uncultured dsDNA bacteriophages were discovered in Dishui Lake, the largest artificial lake in Shanghai. We analyzed their genomic signatures and found a series of viral auxiliary metabolic genes closely associated with protein synthesis and host metabolism. Dishui Lake phages shared more genes with uncultivated environmental viruses than with reference viruses based on the gene-sharing network classification. Phylogeny of proteomes and comparative genomics delineated three new genera within two known viral families of Kyanoviridae and Autographiviridae, and four new families in Caudoviricetes for these seven novel phages. Their potential hosts appeared to be from the dominant bacterial phyla in Dishui Lake. Altogether, our study provides initial insights into the composition and diversity of bacteriophage communities in Dishui Lake, contributing valuable knowledge to the ongoing research on the roles played by viruses in freshwater ecosystems.

New Bacteriophages Members of the Ackermannviridae Family Specific for Klebsiella pneumoniae ST258

Background

Carbapenem-resistant Klebsiella pneumoniae, particularly isolates classified as sequence-type 258 (ST258), are multidrug-resistant strains that are strongly associated with poor-prognosis nosocomial infections, as current therapeutic options are limited and ineffective. In recent years, phage therapy has emerged as a promising treatment option for these scenarios.

Methodology and Results

We report the isolation and characterization of three new phages against Klebsiella pneumoniae ST258 strains recovered from Machángara river wastewater. These new members of the Ackermannviridae family showed stability over a wide temperature and pH range and burst sizes ranging from 6 to 44 plaque-forming units per bacteria. Their genomes were about 157 kilobases, with an average guanine-cytosine content of 46.4% and showed presence of several transfer RNAs, which also allowed us to predict in silico a lytic replicative cycle due to the presence of endolysins and lysozymes.

Conclusion

Three lytic phages of Ackermannviridae family were recovered against Klebsiella pneumoniae ST258 strains from sewage; however, further characterization is needed for future consideration as therapeutic alternatives.

A swapped genetic code prevents viral infections and gene transfer

Engineering the genetic code of an organism has been proposed to provide a firewall from natural ecosystems by preventing viral infections and gene transfer1-6. However, numerous viruses and mobile genetic elements encode parts of the translational apparatus7-9, potentially rendering a genetic-code-based firewall ineffective. Here we show that such mobile transfer RNAs (tRNAs) enable gene transfer and allow viral replication in Escherichia coli despite the genome-wide removal of 3 of the 64 codons and the previously essential cognate tRNA and release factor genes. We then establish a genetic firewall by discovering viral tRNAs that provide exceptionally efficient codon reassignment allowing us to develop cells bearing an amino acid-swapped genetic code that reassigns two of the six serine codons to leucine during translation. This amino acid-swapped genetic code renders cells resistant to viral infections by mistranslating viral proteomes and prevents the escape of synthetic genetic information by engineered reliance on serine codons to produce leucine-requiring proteins. As these cells may have a selective advantage over wild organisms due to virus resistance, we also repurpose a third codon to biocontain this virus-resistant host through dependence on an amino acid not found in nature10. Our results may provide the basis for a general strategy to make any organism safely resistant to all natural viruses and prevent genetic information flow into and out of genetically modified organisms.

Isolation, characterization, and comparative genomic analysis of vB_BviS-A10Y, a novel bacteriophage from mangrove sediments

Mangrove is among the most carbon-rich biomes on earth, and viruses are believed to play a significant role in modulating local and global carbon cycling. However, few viruses have been isolated from mangrove sediments to date. Here, we report the isolation of a novel Bacillus phage (named phage vB_BviS-A10Y) from mangrove sediments. Phage vB_BviS-A10Y has a hexameric head with a diameter of ~ 79.22 nm and a tail with a length of ~ 548.56 nm, which are typical features of siphophages. vB_BviS-A10Y initiated host lysis at 3.5 h postinfection with a burst size of 25 plaque-forming units (PFU)/cell. The genome of phage vB_BviS-A10Y is 162,435 bp long with 225 predicted genes, and the GC content is 34.03%. A comparison of the whole genome sequence of phage vB_BviS-A10Y with those of other phages from the NCBI viral genome database showed that phage vB_BviS-A10Y has the highest similarity (73.7% identity with 33% coverage) to Bacillus phage PBC2. Interestingly, abundant auxiliary metabolic genes (AMGs) were identified in the vB_BviS-A10Y genome. The presence of a β-1,3-glucosyltransferase gene in the phage genome supported our previous hypothesis that mangrove viruses may manipulate carbon cycling directly through their encoded carbohydrate-active enzyme (CAZyme) genes. Therefore, our study will contribute to a better understanding of the diversity and potential roles of viruses in mangrove ecosystems.

Characterization and genome analysis of Escherichia phage fBC-Eco01, isolated from wastewater in Tunisia

The rise of antibiotic resistance in bacterial strains has led to vigorous exploration for alternative treatments. To this end, phage therapy has been revisited, and it is gaining increasing attention, as it may represent an efficient alternative for treating multiresistant pathogenic bacteria. Phage therapy is considered safe, and phages do not infect eukaryotic cells. There have been many studies investigating phage-host bacteria interactions and the ability of phages to target specific hosts. Escherichia coli is the causative agent of a multitude of infections, ranging from urinary tract infections to sepsis, with growing antibiotic resistance. In this study, we characterized the Escherichia phage fBC-Eco01, which was isolated from a water sample collected at Oued, Tunis. Electron microscopy showed that fBC-Eco01 phage particles have siphovirus morphology, with an icosahedral head of 61 ± 3 nm in diameter and a non-contractile tail of 94 ± 2 nm in length and 12 ± 0.9 nm in width. The genome of fBC-Eco01 is a linear double-stranded DNA of 43.466 bp with a GC content of 50.4%. Comparison to databases allowed annotation of the functions to 39 of the 78 predicted gene products. A single-step growth curve revealed that fBC-Eco01 has a latent period of 30 minutes and a burst size of 175 plaque-forming units (PFU) per infected cell. Genomic analysis indicated that fBC-Eco01 is a member of the subfamily Guernseyvirinae. It is most closely related to a group of phages of the genus Kagunavirus that infect Enterobacter, Raoultella, and Escherichia strains.

Isolation, Characterization, and Comparative Genomic Analysis of vB_Pd_C23, a Novel Bacteriophage of Pantoea dispersa

Pantoea bacteria species cause human animal infections, and contribute to soil and aquatic environmental pollution. A novel bacteriophage, vB_Pd_C23 was isolated from a Tunisian wastewater system and represents the first new phage infecting P. dispersa. Lysis kinetics, electron microscopy, and genomic analyses revealed that the vB_Pd_C23 phage has a head diameter of 50 nm and contractile tail dimensions of 100 nm by 23 nm; vB_Pd_C23 has a linear double-stranded DNA genome consisting of 44,714-bp and 49.66% GC-content. Predicted functions were assigned to 75 open reading frames (ORFs) encoding proteins and one tRNA, the annotation revealed that 21 ORFs encode for unique proteins of yet unknown function with no reliable homologies. This indicates that the new species vB_Pd_C23 exhibits novel viral genes. Phylogenetic analysis along with comparative analyses generating nucleotide identity and similarity of vB_Pd_C23 whole genome suggests that the phage is a candidate for a new genus within the Caudoviricetes Class. The characteristics of this phage could not be attributed to any previous genera recognized by the International Committee on Taxonomy of Viruses (ICTV).

Identification and analysis of putative tRNA genes in baculovirus genomes

Transfer RNAs (tRNAs) genes are both coded for and arranged along some viral genomes representing the entire virosphere and seem to play different biological functions during infection, other than transferring the correct amino acid to a growing peptide chain. Baculovirus genome description and annotation has focused mostly on protein-coding genes, microRNA, and homologous regions. Here we carried out a large-scale in silico search for putative tRNA genes in baculovirus genomes. Ninety-six of 257 baculovirus genomes analyzed was found to contain at least one putative tRNA gene. We found great diversity in primary and secondary structure, in location within the genome, in intron presence and size, and in anti-codon identity. In some cases, genes of tRNA-containing genomes were found to have a bias for the codons specified by the tRNAs present in such genomes. Moreover, analysis revealed that most of the putative tRNA genes possessed conserved motifs for tRNA type 2 promoters, including the A-box and B-box motifs with few mismatches from the eukaryotic canonical motifs. From publicly available small RNA deep sequencing datasets of baculovirus-infected insect cells, we found evidence that a putative Autographa californica multiple nucleopolyhedrovirus Gln-tRNA gene was transcribed and modified with the addition of the non-templated 3'-CCA tail found at the end of all tRNAs. Further research is needed to determine the expression and functionality of these viral tRNAs.

A Cocktail of Three Virulent Phages Controls Multidrug-Resistant Salmonella Enteritidis Infection in Poultry

Salmonella enterica is not only the most common pathogen of poultry and poultry-derived products but is also a significant foodborne pathogen. In recent years, many S. enterica isolates have exhibited multi-drug resistance, which places huge pressure on global economy and health. Since phages are an attractive alternative to biocontrol pathogens, we isolated a total of 15 Salmonella phages from sewage effluent, sediment, and chicken manure. The GRNsp1, GRNsp3, GRNsp6, GRNsp21, GRNsp27, GRNsp30, GRNsp50, and GRNsp51 phages exhibited a wide host range against S. enterica serovars Enteritidis and Typhimurium in vitro. In particular, GRNsp51 exerted highly efficient lytic effects against a large proportion of S. Enteritidis and S. Typhimurium strains isolated from different regions of China. Meanwhile, GRNsp8 expanded the host range of GRNsp6 and GRNsp51. Based on their host ranges and lytic capacities, GRNsp6, GRNssp8, and GRNsp51 were selected for further investigation. Morphology, one-step growth curves, and stability assays revealed that GRNsp6, GRNsp8, and GRNsp51 all belong to the Caudovirales order and display relatively short latency periods with broad pH and thermal stability. Genomic analysis indicated that the genomes of these three phages contained no genes related to virulence, antibiotic resistance, or lysogeny. In addition, we tested the effectiveness of a cocktail composed of these three phages against S. Enteritidis in a chicken model. Treatment with the oral phage cocktail 24 h before or alongside Salmonella challenge significantly reduced colonization of the intestinal tract and decreased the mRNA expression of IL-6, IFN-γ, and IL-1β in the duodenum. Together, these findings indicate that a cocktail of the GRNsp6, GRNsp8, and GRNsp51 phages could serve as an effective antimicrobial therapeutic agent against multidrug-resistant Salmonella in animal production to mitigate infections by multiple zoonotic Salmonella species.

Novel Salmonella Phage, vB_Sen_STGO-35-1, Characterization and Evaluation in Chicken Meat

Salmonellosis is one of the most frequently reported zoonotic foodborne diseases worldwide, and poultry is the most important reservoir of Salmonella enterica serovar Enteritidis. The use of lytic bacteriophages (phages) to reduce foodborne pathogens has emerged as a promising biocontrol intervention for Salmonella spp. Here, we describe and evaluate the newly isolated Salmonella phage STGO-35-1, including: (i) genomic and phenotypic characterization, (ii) an analysis of the reduction of Salmonella in chicken meat, and (iii) genome plasticity testing. Phage STGO-35-1 represents an unclassified siphovirus, with a length of 47,483 bp, a G + C content of 46.5%, a headful strategy of packaging, and a virulent lifestyle. Phage STGO-35-1 reduced S. Enteritidis counts in chicken meat by 2.5 orders of magnitude at 4 °C. We identified two receptor-binding proteins with affinity to LPS, and their encoding genes showed plasticity during an exposure assay. Phenotypic, proteomic, and genomic characteristics of STGO-35-1, as well as the Salmonella reduction in chicken meat, support the potential use of STGO-35-1 as a targeted biocontrol agent against S. Enteritidis in chicken meat. Additionally, computational analysis and a short exposure time assay allowed us to predict the plasticity of genes encoding putative receptor-binding proteins.

Streptococcus pneumoniae: a Plethora of Temperate Bacteriophages With a Role in Host Genome Rearrangement

Bacteriophages (phages) are viruses that infect bacteria. They are the most abundant biological entity on Earth (current estimates suggest there to be perhaps 10³¹ particles) and are found nearly everywhere. Temperate phages can integrate into the chromosome of their host, and prophages have been found in abundance in sequenced bacterial genomes. Prophages may modulate the virulence of their host in different ways, e.g., by the secretion of phage-encoded toxins or by mediating bacterial infectivity. Some 70% of Streptococcus pneumoniae (the pneumococcus)—a frequent cause of otitis media, pneumonia, bacteremia and meningitis—isolates harbor one or more prophages. In the present study, over 4000 S. pneumoniae genomes were examined for the presence of prophages, and nearly 90% were found to contain at least one prophage, either defective (47%) or present in full (43%). More than 7000 complete putative integrases, either of the tyrosine (6243) or serine (957) families, and 1210 full-sized endolysins (among them 1180 enzymes corresponding to 318 amino acid-long N-acetylmuramoyl-L-alanine amidases [LytA_PPH]) were found. Based on their integration site, 26 different pneumococcal prophage groups were documented. Prophages coding for tRNAs, putative virulence factors and different methyltransferases were also detected. The members of one group of diverse prophages (PPH090) were found to integrate into the 3’ end of the host lytA_Spn gene encoding the major S. pneumoniae autolysin without disrupting it. The great similarity of the lytA_Spnand lytA_PPH genes (85–92% identity) allowed them to recombine, via an apparent integrase-independent mechanism, to produce different DNA rearrangements within the pneumococcal chromosome. This study provides a complete dataset that can be used to further analyze pneumococcal prophages, their evolutionary relationships, and their role in the pathogenesis of pneumococcal disease.

Characterization of Erwinia tracheiphila Bacteriophage FBB1 Isolated from Spotted Cucumber Beetles that Vector E. tracheiphila

Erwinia tracheiphila, the causal pathogen of bacterial wilt of cucurbit crops, is disseminated by cucumber beetles. A bacteriophage, designated FBB1 (Fu-Beattie-Beetle-1), was isolated from spotted cucumber beetles (Diabrotica undecimpunctata) that were collected from a field in which E. tracheiphila is endemic. FBB1 was classified into the Myoviridae family based on its morphology, which includes an elongated icosahedral head (106 × 82 nm) and a putatively contractile tail (120 nm). FBB1 infected all 62 E. tracheiphila strains examined and three Pantoea spp. strains. FBB1 virions were stable at 55°C for 1 h and tolerated a pH range from 3 to 12. FBB1 has a genome of 175,994 bp with 316 predicted coding sequences and a GC content of 36.5%. The genome contains genes for a major bacterial outer-membrane protein, a putative exopolysaccharide depolymerase, and 22 predicted transfer RNAs. The morphology and genome indicate that FBB1 is a T4-like virus and thus in the Tevenvirinae subfamily. FBB1 is the first virulent phage of E. tracheiphila to be reported and, to date, is one of only two bacteriophages to be isolated from insect vectors of phytopathogens. Collectively, the results support FBB1 as a promising candidate for biocontrol of E. tracheiphila based on its virulent (lytic) rather than lysogenic lifestyle, its infection of all E. tracheiphila strains examined to date, and its infection of a few nonpathogenic bacteria that could be used to support phage populations when pathogen numbers are low.

Characterization and Genome Study of Novel Lytic Bacteriophages against Prevailing Saprophytic Bacterial Microflora of Minimally Processed Plant-Based Food Products

The food industry is still searching for novel solutions to effectively ensure the microbiological safety of food, especially fresh and minimally processed food products. Nowadays, the use of bacteriophages as potential biological control agents in microbiological food safety and preservation is a promising strategy. The aim of the study was the isolation and comprehensive characterization of novel bacteriophages with lytic activity against saprophytic bacterial microflora of minimally processed plant-based food products, such as mixed leaf salads. From 43 phages isolated from municipal sewage, four phages, namely Enterobacter phage KKP 3263, Citrobacter phage KKP 3664, Enterobacter phage KKP 3262, and Serratia phage KKP 3264 have lytic activity against Enterobacter ludwigii KKP 3083, Citrobacter freundii KKP 3655, Enterobacter cloacae KKP 3082, and Serratia fonticola KKP 3084 bacterial strains, respectively. Transmission electron microscopy (TEM) and whole-genome sequencing (WGS) identified Enterobacter phage KKP 3263 as an Autographiviridae, and Citrobacter phage KKP 3664, Enterobacter phage KKP 3262, and Serratia phage KKP 3264 as members of the Myoviridae family. Genome sequencing revealed that these phages have linear double-stranded DNA (dsDNA) with sizes of 39,418 bp (KKP 3263), 61,608 bp (KKP 3664), 84,075 bp (KKP 3262), and 148,182 bp (KKP 3264). No antibiotic resistance genes, virulence factors, integrase, recombinase, or repressors, which are the main markers of lysogenic viruses, were annotated in phage genomes. Serratia phage KKP 3264 showed the greatest growth inhibition of Serratia fonticola KKP 3084 strain. The use of MOI 1.0 caused an almost 5-fold decrease in the value of the specific growth rate coefficient. The phages retained their lytic activity in a wide range of temperatures (from -20 °C to 50 °C) and active acidity values (pH from 4 to 11). All phages retained at least 70% of lytic activity at 60 °C. At 80 °C, no lytic activity against tested bacterial strains was observed. Serratia phage KKP 3264 was the most resistant to chemical factors, by maintaining high lytic activity across a broader range of pH from 3 to 11. The results indicated that these phages could be a potential biological control agent against saprophytic bacterial microflora of minimally processed plant-based food products.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01146-6.

Isolation, Genomic Analysis, and Preliminary Application of a Bovine Klebsiella pneumoniae Bacteriophage vB_Kpn_B01

Klebsiella pneumoniae is an important pathogen that can infect both humans and cattle. The widespread K. pneumoniae and its high drug resistance make it difficult to treat Klebsiella infections/diseases. In this study, a lytic K. pneumoniae bacteriophage vB_Kpn_B01 was isolated from a dairy farm trough in Sichuan Province, and its biological properties were studied, and the entire genome of vB_Kpn_B01 was sequenced. The therapeutic effects of the phage on disease-causing mice were preliminarily tested. Phages found in this study are double-stranded DNA bacterial viruses belonging to the family Siphoviridae, Sugarlandvirus. The results suggest that vB_Kpn_B01 has strong specificity and low adaptability to different adverse conditions. Meanwhile, the predicted gene products of phage vB_Kpn_B01 comprised 149 coding sequences (CDS) and 25 tRNAs, of which 34 CDS had known functions. Of course, vB_Kpn_B01 did not contain any known antibiotic-resistant or virulent genes. The pathological sections of the liver and lungs of mice showed that the inflammatory scores of the treatment group were lower than in the bacterial group. Phage vB_Kpn_B01 alleviated the inflammatory response in the organs of the infected mice, and the organ tissue bacterial load of the treatment group was significantly lower than that of the bacterial group. Therefore, vB_Kpn_B01 can inhibit the proliferation of K. pneumoniae 18 in vivo and can alleviate the inflammation of target organs caused by infectious bacteria, which preliminarily indicates that vB_Kpn_B01 has a certain therapeutic effect on laboratory-infected mice.

Role of a cryptic tRNA gene operon in survival under translational stress

As compared to eukaryotes, bacteria have a reduced tRNA gene set encoding between 30 and 220 tRNAs. Although in most bacterial phyla tRNA genes are dispersed in the genome, many species from distinct phyla also show genes forming arrays. Here, we show that two types of arrays with distinct evolutionary origins exist. This work focuses on long tRNA gene arrays (L-arrays) that encompass up to 43 genes, which disseminate by horizontal gene transfer and contribute supernumerary tRNA genes to the host. Although in the few cases previously studied these arrays were reported to be poorly transcribed, here we show that the L-array of the model cyanobacterium Anabaena sp. PCC 7120, encoding 23 functional tRNAs, is largely induced upon impairment of the translation machinery. The cellular response to this challenge involves a global reprogramming of the transcriptome in two phases. tRNAs encoded in the array are induced in the second phase of the response, directly contributing to cell survival. Results presented here show that in some bacteria the tRNA gene set may be partitioned between a housekeeping subset, which constantly sustains translation, and an inducible subset that is generally silent but can provide functionality under particular conditions.

Characterization of Two New Shiga Toxin-Producing Escherichia coli O103-Infecting Phages Isolated from an Organic Farm

Shiga toxin-producing Escherichia coli (STEC) O103 strains have been recently attributed to various foodborne outbreaks in the United States. Due to the emergence of antibiotic-resistant strains, lytic phages are considered as alternative biocontrol agents. This study was to biologically and genomically characterize two STEC O103-infecting bacteriophages, vB_EcoP-Ro103C3lw (or Ro103C3lw) and vB_EcoM-Pr103Blw (or Pr103Blw), isolated from an organic farm. Based on genomic and morphological analyses, phages Ro103C3lw and Pr103Blw belonged to Autographiviridae and Myoviridae families, respectively. Ro103C3lw contained a 39,389-bp double-stranded DNA and encoded a unique tail fiber with depolymerase activity, resulting in huge plaques. Pr103Blw had an 88,421-bp double-stranded DNA with 26 predicted tRNAs associated with the enhancement of the phage fitness. Within each phage genome, no virulence, antibiotic-resistant, and lysogenic genes were detected. Additionally, Ro103C3lw had a short latent period (2 min) and a narrow host range, infecting only STEC O103 strains. By contrast, Pr103Blw had a large burst size (152 PFU/CFU) and a broad host range against STEC O103, O26, O111, O157:H7, and Salmonella Javiana strains. Furthermore, both phages showed strong antimicrobial activities against STEC O103:H2 strains. The findings provide valuable insight into these two phages’ genomic features with the potential antimicrobial activities against STEC O103.

Comparative Analysis of Felixounavirus Genomes Including Two New Members of the Genus That Infect Salmonella Infantis

Salmonella spp. is one of the most common foodborne pathogens worldwide; therefore, its control is highly relevant for the food industry. Phages of the Felixounavirus genus have the characteristic that one phage can infect a large number of different Salmonella serovars and, thus, are proposed as an alternative to antimicrobials in food production. Here, we describe two new members of the Felixounavirus genus named vB_Si_35FD and vB_Si_DR94, which can infect Salmonella Infantis. These new members were isolated and sequenced, and a subsequent comparative genomic analysis was conducted including 23 publicly available genomes of Felixounaviruses that infect Salmonella. The genomes of vB_Si_35FD and vB_Si_DR94 are 85,818 and 85,730 bp large and contain 129 and 125 coding sequences, respectively. The genomes did not show genes associated with virulence or antimicrobial resistance, which could be useful for candidates to use as biocontrol agents. Comparative genomics revealed that closely related Felixounavirus are found in distinct geographical locations and that this genus has a conserved genomic structure despite its worldwide distribution. Our study revealed a highly conserved structure of the phage genomes, and the two newly described phages could represent promising biocontrol candidates against Salmonella spp. from a genomic viewpoint.

Translational adaptation of human viruses to the tissues they infect

Viruses need to hijack the translational machinery of the host cell for a productive infection to happen. However, given the dynamic landscape of tRNA pools among tissues, it is unclear whether different viruses infecting different tissues have adapted their codon usage toward their tropism. Here, we collect the coding sequences of 502 human-infecting viruses and determine that tropism explains changes in codon usage. Using the tRNA abundances across 23 human tissues from The Cancer Genome Atlas (TCGA), we build an in silico model of translational efficiency that validates the correspondence of the viral codon usage with the translational machinery of their tropism. For instance, we detect that severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is specifically adapted to the upper respiratory tract and alveoli. Furthermore, this correspondence is specifically defined in early viral proteins. The observed tissue-specific translational efficiency could be useful for the development of antiviral therapies and vaccines.

On the Track of the Missing tRNA Genes: A Source of Non-Canonical Functions?

Cellular tRNAs appear today as a diverse population of informative macromolecules with conserved general elements ensuring essential common functions and different and distinctive features securing specific interactions and activities. Their differential expression and the variety of post-transcriptional modifications they are subject to, lead to the existence of complex repertoires of tRNA populations adjusted to defined cellular states. Despite the tRNA-coding genes redundancy in prokaryote and eukaryote genomes, it is surprising to note the absence of genes coding specific translational-active isoacceptors throughout the phylogeny. Through the analysis of different releases of tRNA databases, this review aims to provide a general summary about those “missing tRNA genes.” This absence refers to both tRNAs that are not encoded in the genome, as well as others that show critical sequence variations that would prevent their activity as canonical translation adaptor molecules. Notably, while a group of genes are universally missing, others are absent in particular kingdoms. Functional information available allows to hypothesize that the exclusion of isodecoding molecules would be linked to: 1) reduce ambiguities of signals that define the specificity of the interactions in which the tRNAs are involved; 2) ensure the adaptation of the translational apparatus to the cellular state; 3) divert particular tRNA variants from ribosomal protein synthesis to other cellular functions. This leads to consider the “missing tRNA genes” as a source of putative non-canonical tRNA functions and to broaden the concept of adapter molecules in ribosomal-dependent protein synthesis.

Comprehensive in silico survey of the Mycolicibacterium mobilome reveals an as yet underexplored diversity

The mobilome plays a crucial role in bacterial adaptation and is therefore a starting point to understand and establish the gene flow occurring in the process of bacterial evolution. This is even more so if we consider that the mobilome of environmental bacteria can be the reservoir of genes that may later appear in the clinic. Recently, new genera have been proposed in the family Mycobacteriaceae, including the genus Mycolicibacterium, which encompasses dozens of species of agricultural, biotechnological, clinical and ecological importance, being ubiquitous in several environments. The current scenario in the Mycobacteriaceae mobilome has some bias because most of the characterized mycobacteriophages were isolated using a single host strain, and the few plasmids reported mainly relate to the genus Mycobacterium. To fill in the gaps in these issues, we performed a systematic in silico study of these mobile elements based on 242 available genomes of the genus Mycolicibacterium. The analyses identified 156 putative plasmids (19 conjugative, 45 mobilizable and 92 non-mobilizable) and 566 prophages in 86 and 229 genomes, respectively. Moreover, a contig was characterized by resembling an actinomycete integrative and conjugative element (AICE). Within this diversity of mobile genetic elements, there is a pool of genes associated with several canonical functions, in addition to adaptive traits, such as virulence and resistance to antibiotics and metals (mercury and arsenic). The type-VII secretion system was a common feature in the predicted plasmids, being associated with genes encoding virulent proteins (EsxA, EsxB, PE and PPE). In addition to the characterization of plasmids and prophages of the family Mycobacteriaceae, this study showed an abundance of these genetic elements in a dozen species of the genus Mycolicibacterium.

Genomic Characterization of Two Shiga Toxin-Converting Bacteriophages Induced From Environmental Shiga Toxin-Producing Escherichia coli

Shiga toxin (Stx), encoded by stx genes located in prophage sequences, is the major agent responsible for the pathogenicity of Shiga toxin-producing Escherichia coli (STEC) and is closely associated with the development of hemolytic uremic syndrome (HUS). Although numerous Stx prophage sequences have been reported as part of STEC bacterial genomes, the information about the genomic characterization of Stx-converting bacteriophages induced from STEC strains is relatively scarce. The objectives of this study were to genomically characterize two Stx-converting phages induced from environmental STEC strains and to evaluate their correlations with published Stx-converting phages and STEC strains of different origins. The Stx1-converting phage Lys8385Vzw and the Stx2-converting phage Lys19259Vzw were induced from E. coli O103:H11 (RM8385) and E. coli O157:H7 (RM19259), respectively. Whole-genome sequencing of these phages was conducted on a MiSeq sequencer for genomic characterization. Phylogenetic analysis and comparative genomics were performed to determine the correlations between these two Stx-converting phages, 13 reference Stx-converting phages, and 10 reference STEC genomes carrying closely related Stx prophages. Both Stx-converting phages Lys8385Vzw and Lys19259Vzw had double-stranded DNA, with genome sizes of 50,953 and 61,072 bp, respectively. Approximately 40% of the annotated coding DNA sequences with the predicted functions were likely associated with the fitness for both phages and their bacterial hosts. The whole-genome–based phylogenetic analysis of these two Stx-converting phages and 13 reference Stx-converting phages revealed that the 15 Stx-converting phages were divided into three distinct clusters, and those from E. coli O157:H7, in particular, were distributed in each cluster, demonstrating the high genomic diversity of these Stx-converting phages. The genomes of Stx-converting phage Lys8385Vzw and Lys19259Vzw shared a high-nucleotide similarity with the prophage sequences of the selected STEC isolates from the clinical and environmental origin. The findings demonstrate the genomic diversity of Stx-converting phages induced from different STEC strains and provide valuable insights into the dissemination of stx genes among E. coli population via the lysogenization of Stx-converting phages.

Diversity of tRNA Clusters in the Chloroviruses

Viruses rely on their host’s translation machinery for the synthesis of their own proteins. Problems belie viral translation when the host has a codon usage bias (CUB) that is different from an infecting virus due to differences in the GC content between the host and virus genomes. Here, we examine the hypothesis that chloroviruses adapted to host CUB by acquisition and selection of tRNAs that at least partially favor their own CUB. The genomes of 41 chloroviruses comprising three clades, each infecting a different algal host, have been sequenced, assembled and annotated. All 41 viruses not only encode tRNAs, but their tRNA genes are located in clusters. While differences were observed between clades and even within clades, seven tRNA genes were common to all three clades of chloroviruses, including the tRNA^Arg gene, which was found in all 41 chloroviruses. By comparing the codon usage of one chlorovirus algal host, in which the genome has been sequenced and annotated (67% GC content), to that of two of its viruses (40% GC content), we found that the viruses were able to at least partially overcome the host’s CUB by encoding tRNAs that recognize AU-rich codons. Evidence presented herein supports the hypothesis that a chlorovirus tRNA cluster was present in the most recent common ancestor (MRCA) prior to divergence into three clades. In addition, the MRCA encoded a putative isoleucine lysidine synthase (TilS) that remains in 39/41 chloroviruses examined herein, suggesting a strong evolutionary pressure to retain the gene. TilS alters the anticodon of tRNA^Met that normally recognizes AUG to then recognize AUA, a codon for isoleucine. This is advantageous to the chloroviruses because the AUA codon is 12–13 times more common in the chloroviruses than their host, further helping the chloroviruses to overcome CUB. Among large DNA viruses infecting eukaryotes, the presence of tRNA genes and tRNA clusters appear to be most common in the Phycodnaviridae and, to a lesser extent, in the Mimiviridae.

Complete genome sequence of a novel Bacillus phage, P59, that infects Bacillus oceanisediminis

P59, a virulent phage of Bacillus oceanisediminis, was isolated from the sediment of Weiming Lake at Peking University (Beijing, China). P59 showed the typical morphology of myovirids. The complete genome sequence of P59 is 159,363 bp in length with a G+C content of 42.34%. The genome sequence has very low similarity to the other phage genome sequences in the GenBank database, suggesting that P59 is a new phage. A total of 261 open reading frames and 15 tRNA genes were predicted. Based on its morphological and genetic traits, we propose phage P59 to be a new member of the family Herelleviridae.

Pantoeaagglomerans-Infecting Bacteriophage vB_PagS_AAS21: A Cold-Adapted Virus Representing a Novel Genus within the Family Siphoviridae

A novel cold-adapted siphovirus, vB_PagS_AAS21 (AAS21), was isolated in Lithuania using Pantoea agglomerans as the host for phage propagation. AAS21 has an isometric head (~85 nm in diameter) and a non-contractile flexible tail (~174 × 10 nm). With a genome size of 116,649 bp, bacteriophage AAS21 is the largest Pantoea-infecting siphovirus sequenced to date. The genome of AAS21 has a G+C content of 39.0% and contains 213 putative protein-encoding genes and 29 genes for tRNAs. A comparative sequence analysis revealed that 89 AAS21 open reading frames (ORFs) code for unique proteins that have no reliable identity to database entries. In total, 63 AAS21 ORFs were functionally annotated, including those coding for the proteins responsible for virion morphogenesis, phage-host interactions, and DNA metabolism. Proteomic analysis led to the experimental identification of 19 virion proteins, including 11 that were predicted by bioinformatics approaches. Based on comparative phylogenetic analysis, AAS21 cannot be assigned to any genus currently recognized by ICTV and may represents a new branch of viruses within the family Siphoviridae.

The rare lncRNA GOLLD is widespread and structurally conserved among Mycobacterium tRNA arrays

Noncoding RNA (ncRNA) genes produce transcripts involved in a wide range of functions, including catalytic and regulatory functions. Besides, some transcripts have highly complex structures that may impact their activities. Among the largest bacterial ncRNAs, there is the rare GOLLD RNA, which is associated with tRNA genes and supposed to be chromosome- and phage-encoded in specialized groups of bacteria, including those from Lactobacillales and Actinomycetales orders. The only GOLLD structure was inferred from a variety of sequences, including many marine metagenomes. To explore GOLLD RNA in bacterial genomes, we mined the GOLLD gene in thousands of Mycobacterium and virus genomes using Infernal software. We identified this gene in 350 mycobacteria, including megaplasmids, and 39 bacteriophages, mainly in the genomic context of tRNA arrays. Mycobacterium GOLLD genes presented a high diversity and were distributed in three phylogenetic groups: (i) Mycobacterium exclusive; (ii) Mycobacterium and mycobacteriophages; and (iii) mycobacteriophage exclusive. We also determined the GOLLD secondary structure of each group using R2 R software based on GOLLD alignments generated by Infernal software. All GOLLD groups displayed a 3' half conserved structure, including utter E-loops pseudoknots substructures, also shared by non-Mycobacterium GOLLD while the 5' half motif was different among the groups. Here, we showed that the lncRNA GOLLD is widespread in Mycobacterium within tRNA arrays and corroborated the previously predicted GOLLD secondary structure.