Mining bacterial NGS data vastly expands the complete genomes of temperate phages

Gut phageome in Mexican Americans: a population at high risk for metabolic dysfunction-associated steatotic liver disease and diabetes

Mexican Americans are disproportionally affected by metabolic dysfunction-associated steatotic liver disease (MASLD), which often co-occurs with diabetes. Despite extensive evidence on the causative role of the gut microbiome in MASLD, studies determining the involvement of the gut phageome are scarce. In this cross-sectional study, we characterized the gut phageome in Mexican Americans of South Texas by stool shotgun metagenomic sequencing of 340 subjects, concurrently screened for liver steatosis by transient elastography. Inter-individual variations in the phageome were associated with gender, country of birth, diabetes, and liver steatosis. The phage signatures for diabetes and liver steatosis were subsequently determined. Enrichment of Inoviridae was associated with both diabetes and liver steatosis. Diabetes was further associated with the enrichment of predominantly temperate Escherichia phages, some of which possessed virulence factors. Liver steatosis was associated with the depletion of Lactococcus phages r1t and BK5-T, and enrichment of the globally prevalent Crassvirales phages, including members of genus cluster IX (Burzaovirus coli, Burzaovirus faecalis) and VI (Kahnovirus oralis). The Lactococcus phages showed strong correlations and co-occurrence with Lactococcus lactis, while the Crassvirales phages, B. coli, B. faecalis, and UAG-readthrough crAss clade correlated and co-occurred with Prevotella copri. In conclusion, we identified the gut phageome signatures for two closely linked metabolic diseases with significant global burden. These phage signatures may have utility in risk modeling and disease prevention in this high-risk population, and identification of potential bacterial targets for phage therapy.IMPORTANCEPhages influence human health and disease by shaping the gut bacterial community. Using stool samples from a high-risk Mexican American population, we provide insights into the gut phageome changes associated with diabetes and liver steatosis, two closely linked metabolic diseases with significant global burden. Common to both diseases was an enrichment of Inoviridae, a group of phages that infect bacterial hosts chronically without lysis, allowing them to significantly influence bacterial growth, virulence, motility, biofilm formation, and horizontal gene transfer. Diabetes was additionally associated with the enrichment of Escherichia coli-infecting phages, some of which contained virulence factors. Liver steatosis was additionally associated with the depletion of Lactococcus lactis-infecting phages, and enrichment of Crassvirales phages, a group of virulent phages with high global prevalence and persistence across generations. These phageome signatures may have utility in risk modeling, as well as identify potential bacterial targets for phage therapy.

The trade-off of Vibrio parahaemolyticus between bacteriophage resistance and growth competitiveness

Vibrio parahaemolyticus is a food-borne pathogen, which is often isolated from various seafood products. In this study, two kinds of bacteriophages was isolated from the offshore sediments samples. The anti-phage mutant strain were obtained after seventeen rounds of co-culture of Vibrio parahaemolyticus and mixed bacteriophage, multigroup sequencing was carried out on spontaneous the anti-phage mutant strain and the wild-type strain. We used the Sanger sequencing to verify the accuracy of the mutation sites. Biolog GEN III MicroPlates were used to evaluate the metabolic capacity of wild-type strains and the anti-phage mutant strain. In this study, we found that with flaG gene (slight homology to N terminus of multiple flagellins) mutated, making the bacteriophage unable to absorb to the cell surface of the host. And, the growth competitiveness of the anti-phage mutant strain is lower than the wild-type strain. These results indicated that the fitness cost, including loss of the growth competitiveness, constitutes a barrier to the prevalence of these defense mechanisms. And the selection pressure on different anti-phage strategies depends on the trade-off between mortality imposed by bacteriophages and fitness cost of the defense strategy under the given environmental conditions. In conclusion, this study provides valuable insights into the phage-host interaction and phage resistance in Vibrio parahaemolyticus. Our study provided knowledge for the evolutionary adaption of bacteria against the bacteriophage, which could add more information to understand the phage resistance mechanism before applying in the industry.

ProkBERT family: genomic language models for microbiome applications

Background

In the evolving landscape of microbiology and microbiome analysis, the integration of machine learning is crucial for understanding complex microbial interactions, and predicting and recognizing novel functionalities within extensive datasets. However, the effectiveness of these methods in microbiology faces challenges due to the complex and heterogeneous nature of microbial data, further complicated by low signal-to-noise ratios, context-dependency, and a significant shortage of appropriately labeled datasets. This study introduces the ProkBERT model family, a collection of large language models, designed for genomic tasks. It provides a generalizable sequence representation for nucleotide sequences, learned from unlabeled genome data. This approach helps overcome the above-mentioned limitations in the field, thereby improving our understanding of microbial ecosystems and their impact on health and disease.

Methods

ProkBERT models are based on transfer learning and self-supervised methodologies, enabling them to use the abundant yet complex microbial data effectively. The introduction of the novel Local Context-Aware (LCA) tokenization technique marks a significant advancement, allowing ProkBERT to overcome the contextual limitations of traditional transformer models. This methodology not only retains rich local context but also demonstrates remarkable adaptability across various bioinformatics tasks.

Results

In practical applications such as promoter prediction and phage identification, the ProkBERT models show superior performance. For promoter prediction tasks, the top-performing model achieved a Matthews Correlation Coefficient (MCC) of 0.74 for E. coli and 0.62 in mixed-species contexts. In phage identification, ProkBERT models consistently outperformed established tools like VirSorter2 and DeepVirFinder, achieving an MCC of 0.85. These results underscore the models' exceptional accuracy and generalizability in both supervised and unsupervised tasks.

Conclusions

The ProkBERT model family is a compact yet powerful tool in the field of microbiology and bioinformatics. Its capacity for rapid, accurate analyses and its adaptability across a spectrum of tasks marks a significant advancement in machine learning applications in microbiology. The models are available on GitHub (https://github.com/nbrg-ppcu/prokbert) and HuggingFace (https://huggingface.co/nerualbioinfo) providing an accessible tool for the community.

PhageScope: a well-annotated bacteriophage database with automatic analyses and visualizations

Bacteriophages are viruses that infect bacteria or archaea. Understanding the diverse and intricate genomic architectures of phages is essential to study microbial ecosystems and develop phage therapy strategies. However, the existing phage databases are short of meticulous annotations. To this end, we propose PhageScope (https://phagescope.deepomics.org), an online phage database with comprehensive annotations. PhageScope harbors a collection of 873 718 phage sequences from various sources. Applying fifteen state-of-the-art tools to perform systematic annotations and analyses, PhageScope provides annotations on genome completeness, host range, lifestyle information, taxonomy classification, nine types of structural and functional genetic elements, and three types of comparative genomic studies for curated phages. Additionally, PhageScope incorporates automatic analyses and visualizations for curated and customized phages, serving as an efficient platform for phage study.

Identification of hidden N4-like viruses and their interactions with hosts

IMPORTANCE

The findings of this study are significant, as N4-like viruses represent a unique viral lineage with a distinct replication mechanism and a conserved core genome. This work has resulted in a comprehensive global map of the entire N4-like viral lineage, including information on their distribution in different biomes, evolutionary divergence, genomic diversity, and the potential for viral-mediated host metabolic reprogramming. As such, this work significantly contributes to our understanding of the ecological function and viral-host interactions of bacteriophages.

Genetic characteristics and integration specificity of Salmonella enterica temperate phages

Introduction

Temperate phages can engage in the horizontal transfer of functional genes to their bacterial hosts. Thus, their genetic material becomes an intimate part of bacterial genomes and plays essential roles in bacterial mutation and evolution. Specifically, temperate phages can naturally transmit genes by integrating their genomes into the bacterial host genomes via integrases. Our previous study showed that Salmonella enterica contains the largest number of temperate phages among all publicly available bacterial species. S. enterica is an important pathogen that can cause serious systemic infections and even fatalities.

Methods

Initially, we extracted all S. enterica temperate phages from the extensively developed temperate phage database established in our previous study. Subsequently, we conducted an in-depth analysis of the genetic characteristics and integration specificity exhibited by these S. enterica temperate phages.

Results

Here we identified 8,777 S. enterica temperate phages, all of which have integrases in their genomes. We found 491 non-redundant S. enterica temperate phage integrases (integrase entries). S. enterica temperate phage integrases were classified into three types: intA, intS, and phiRv2. Correlation analysis showed that the sequence lengths of S. enterica integrase and core regions of attB and attP were strongly correlated. Further phylogenetic analysis and taxonomic classification indicated that both the S. enterica temperate phage genomes and the integrase gene sequences were of high diversities.

Discussion

Our work provides insight into the essential integration specificity and genetic diversity of S. enterica temperate phages. This study paves the way for a better understanding of the interactions between phages and S. enterica. By analyzing a large number of S. enterica temperate phages and their integrases, we provide valuable insights into the genetic diversity and prevalence of these elements. This knowledge has important implications for developing targeted therapeutic interventions, such as phage therapy, to combat S. enterica infections. By harnessing the lytic capabilities of temperate phages, they can be engineered or utilized in phage cocktails to specifically target and eradicate S. enterica strains, offering an alternative or complementary approach to traditional antibiotic treatments. Our study has implications for public health and holds potential significance in combating clinical infections caused by S. enterica.

Occurrence and genetic diversity of prophage sequences identified in the genomes of L. casei group bacteria

It is widely believed that microorganisms belonging to L. casei group can have positive effects on the human body. Therefore, these bacteria are used in many industrial processes, including the production of dietary supplements and probiotic preparations. When using live microorganisms in technological processes, it is important to use those without phage sequences within their genomes that can ultimately lead to lysis of the bacteria. It has been shown that many prophages have a benign nature, meaning that they don't directly lead to lysis or inhibit microbial growth. Moreover, the presence of phage sequences in the genomes of these bacteria increases their genetic diversity, which may contribute to easier colonization of new ecological niches. In the 439 analyzed genomes of the L. casei group, 1509 sequences of prophage origin were detected. The average length of intact prophage sequences analyzed was just under 36 kb. GC content of tested sequences was similar for all analyzed species (44.6 ± 0.9%). Analyzing the protein coding sequences collectively, it was found that there was an average of 44 putative ORFs per genome, while the ORF density of all phage genomes varied from 0.5 to 2.1. The average nucleotide identity calculated on sequence alignments for analyzed sequences was 32.7%. Of the 56 L. casei strains used in the next part of the study, 32 did not show culture growth above the OD600 value of 0.5, even at a mitomycin C concentration of 0.25 μg/ml. Primers used for this study allowed for the detection of prophage sequences for over 90% of tested bacterial strains. Finally, prophages of selected strains were induced using mitomycin C, phage particles were isolated and then genomes of viruses obtained were sequenced and analyzed.

Characterization of the novel temperate Staphylococcus haemolyticus phage IME1365_01

The presence of a novel functional prophage, IME1365_01, was predicted from bacterial high-throughput sequencing data and then successfully induced from Staphylococcus haemolyticus by mitomycin C treatment. Transmission electron microscopy showed that phage IME1365_01 has an icosahedral head (43 nm in diameter) and a long tail (172 nm long). This phage possesses a double-stranded DNA genome of 44,875 bp with a G+C content of 35.35%. A total of 63 putative open reading frames (ORFs) were identified in its genome. BLASTn analysis revealed that IME1365_01 is similar to Staphylococcus phage vB_SepS_E72, but with a genome homology coverage of only 26%. The phage genome does not have fixed termini. In ORF24 of phage IME1365_01, a conserved Toll-interleukin-1 receptor domain of the TIR_2 superfamily (accession no. c123749) is located at its N-terminus, and this might serve as a component of an anti-bacterial system. In conclusion, we developed a platform to obtain active temperate phage from prediction, identification, and induction from its bacterial host. After mass screening using this platform, numerous temperate phages and their innate anti-bacterial elements can provide extensive opportunities for therapy against bacterial (especially drug-resistant bacterial) infections.