Host life-history traits influence the distribution of prophages and the genes they carry

Mobile genetic elements define the non-random structure of the Salmonella enterica serovar Typhi pangenome

Bacterial relatedness measured using select chromosomal loci forms the basis of public health genomic surveillance. While approximating vertical evolution through this approach has proven exceptionally valuable for understanding pathogen dynamics, it excludes a fundamental dimension of bacterial evolution-horizontal gene transfer. Incorporating the accessory genome is the logical remediation and has recently shown promise in expanding epidemiological resolution for enteric pathogens. Employing k-mer-based Jaccard index analysis, and a novel genome length distance metric, we computed pangenome (i.e., core and accessory) relatedness for the globally important pathogen Salmonella enterica serotype Typhi (Typhi), and graphically express both vertical (homology-by-descent) and horizontal (homology-by-admixture) evolutionary relationships in a reticulate network of over 2,200 U.S. Typhi genomes. This analysis revealed non-random structure in the Typhi pangenome that is driven predominantly by the gain and loss of mobile genetic elements, confirming and expanding upon known epidemiological patterns, revealing novel plasmid dynamics, and identifying avenues for further genomic epidemiological exploration. With an eye to public health application, this work adds important biological context to the rapidly improving ways of analyzing bacterial genetic data and demonstrates the value of the accessory genome to infer pathogen epidemiology and evolution.IMPORTANCEGiven bacterial evolution occurs in both vertical and horizontal dimensions, inclusion of both core and accessory genetic material (i.e., the pangenome) is a logical step toward a more thorough understanding of pathogen dynamics. With an eye to public, and indeed, global health relevance, we couple contemporary tools for genomic analysis with decades of research on mobile genetic elements to demonstrate the value of the pangenome, known and unknown, annotated, and hypothetical, for stratification of Salmonella enterica serovar Typhi (Typhi) populations. We confirm and expand upon what is known about Typhi epidemiology, plasmids, and antimicrobial resistance dynamics, and offer new avenues of exploration to further deduce Typhi ecology and evolution, and ultimately to reduce the incidence of human disease.

Exploring the diversity and evolutionary strategies of prophages in Hyphomicrobiales, comparing animal-associated with non-animal-associated bacteria

The Hyphomicrobiales bacterial order (previously Rhizobiales) exhibits a wide range of lifestyle characteristics, including free-living, plant-association, nitrogen-fixing, and association with animals (Bartonella and Brucella). This study explores the diversity and evolutionary strategies of bacteriophages within the Hyphomicrobiales order, comparing animal-associated (AAB) with non-animal-associated bacteria (NAAB). We curated 560 high-quality complete genomes of 58 genera from this order and used the PHASTER server for prophage annotation and classification. For 19 genera with representative genomes, we curated 96 genomes and used the Defense-Finder server to summarize the type of anti-phage systems (APS) found in this order. We analyzed the genetic repertoire and length distributions of prophages, estimating evolutionary rates and comparing intact, questionable, and incomplete prophages in both groups. Analyses of best-fit parameters and bootstrap sensitivity were used to understand the evolutionary processes driving prophage gene content. A total of 1860 prophages distributed in Hyphomicrobiales were found, 695 in AAB and 1165 in the NAAB genera. The results revealed a similar number of prophages per genome in AAB and NAAB and a similar length distribution, suggesting shared mechanisms of genetic acquisition of prophage genes. Changes in the frequency of specific gene classes were observed between incomplete and intact prophages, indicating preferential loss or enrichment in both groups. The analysis of best-fit parameters and bootstrap sensitivity tests indicated a higher selection coefficient, induction rate, and turnover in NAAB genomes. We found 68 types of APS in Hyphomicrobiales; restriction modification (RM) and abortive infection (Abi) were the most frequent APS found for all Hyphomicrobiales, and within the AAB group. This classification of APS showed that NAAB genomes have a greater diversity of defense systems compared to AAB, which could be related to the higher rates of prophage induction and turnover in the latter group. Our study provides insights into the distributions of both prophages and APS in Hyphomicrobiales genomes, demonstrating that NAAB carry more defense systems against phages, while AAB show increased prophage stability and an increased number of incomplete prophages. These results suggest a greater role for domesticated prophages within animal-associated bacteria in Hyphomicrobiales.

Diverse Prophage Elements of Salmonella enterica Serovars Show Potential Roles in Bacterial Pathogenicity

Nontyphoidal salmonellosis is an important foodborne and zoonotic infection that causes significant global public health concern. Diverse serovars are multidrug-resistant and encode several virulence indicators; however, little is known on the role prophages play in driving these traits. Here, we extracted prophages from seventy-five Salmonella genomes which represent the fifteen important serovars in the United Kingdom. We analyzed the intact prophages for the presence of virulence genes and established their genomic relationships. We identified 615 prophages from the Salmonella strains, from which 195 prophages are intact, 332 are incomplete, while 88 are questionable. The average prophage carriage was found to be 'extreme' in S. Heidelberg, S. Inverness, and S. Newport (10.2-11.6 prophages/strain), 'high' in S. Infantis, S. Stanley, S. Typhimurium, and S. Virchow (8.2-9.0 prophages/strain), 'moderate' in S. Agona, S. Braenderup, S. Bovismorbificans, S. Choleraesuis, S. Dublin, and S. Java (6.0-7.8 prophages/strain), and 'low' in S. Javiana and S. Enteritidis (5.8 prophages/strain). Cumulatively, 61 virulence genes (1500 gene copies) were detected from representative intact prophages and linked to Salmonella delivery/secretion system (42.62%), adherence (32.7%), magnesium uptake (3.88%), regulation (5%), stress/survival (1.6%), toxins (10%), and antivirulence (1.6%). Diverse clusters were formed among the intact prophages and with bacteriophages of other enterobacteria, suggesting different lineages and associations. Our work provides a strong body of data to support the contributions diverse prophages make to the pathogenicity of Salmonella, including thirteen previously unexplored serovars.

Phage-plasmids promote recombination and emergence of phages and plasmids

Phages and plasmids are regarded as distinct types of mobile genetic elements that drive bacterial evolution by horizontal gene transfer. However, the distinction between both types is blurred by the existence of elements known as prophage-plasmids or phage-plasmids, which transfer horizontally between cells as viruses and vertically within cellular lineages as plasmids. Here, we study gene flow between the three types of elements. We show that the gene repertoire of phage-plasmids overlaps with those of phages and plasmids. By tracking recent recombination events, we find that phage-plasmids exchange genes more frequently with plasmids than with phages, and that direct gene exchange between plasmids and phages is less frequent in comparison. The results suggest that phage-plasmids can mediate gene flow between plasmids and phages, including exchange of mobile element core functions, defense systems, and antibiotic resistance. Moreover, a combination of gene transfer and gene inactivation may result in the conversion of elements. For example, gene loss turns P1-like phage-plasmids into integrative prophages or into plasmids (that are no longer phages). Remarkably, some of the latter have acquired conjugation-related functions to became mobilisable by conjugation. Thus, our work indicates that phage-plasmids can play a key role in the transfer of genes across mobile elements within their hosts, and can act as intermediates in the conversion of one type of element into another.

Novel lineages of single-stranded DNA phages that coevolved with the symbiotic bacteria Rhizobium

This study describes novel single-stranded DNA phages isolated from common bean agriculture soils by infection of the nitrogen-fixing symbiotic bacteria Rhizobium etli and R. phaseoli. A total of 29 phages analyzed have 4.3–6 kb genomes in size and GC 59–60%. They belong to different clades unrelated to other Microviridae subfamilies. Three-dimensional models of the major capsid protein (MCP) showed a conserved β-barrel structural “jelly-roll” fold. A variable-length loop in the MCPs distinguished three Rhizobium microvirus groups. Microviridae subfamilies were consistent with viral clusters determined by the protein-sharing network. All viral clusters, except for Bullavirinae, included mostly microviruses identified in metagenomes from distinct ecosystems. Two Rhizobium microvirus clusters, chaparroviruses, and chicoviruses, were included within large viral unknown clusters with microvirus genomes identified in diverse metagenomes. A third Rhizobium microvirus cluster belonged to the subfamily Amoyvirinae. Phylogenetic analysis of the MCP confirms the divergence of the Rhizobium microviruses into separate clades. The phylogeny of the bacterial hosts matches the microvirus MCP phylogeny, suggesting a coevolutionary history between the phages and their bacterial host. This study provided essential biological information on cultivated microvirus for understanding the evolution and ecological diversification of the Microviridae family in diverse microbial ecosystems.

Genetic Characteristics and Enzymatic Activities of Bacillus velezensis KS04AU as a Stable Biocontrol Agent against Phytopathogens

Bacillus velezensis has a broad application in the agricultural and industrial sectors for its biocontrol properties and its potential active secondary metabolites. The defined phenotypic characteristics of a strain vary according to its ecosystem. We report the complete genomic analysis of B. velezensis KS04AU compared to four strains of B. velezensis (SRCM102752, ONU-553, FZB42, and JS25R) and two closely related Bacillus amyloliquefaciens (LL3 and IT-45). A total of 4771 protein coding genes comprises the KS04AU genome, in comparison with 3334 genes core genes found in the six other strains and the remaining 1437 shell genes. Average nucleotide identity of the target strain to the six other strains showed 99.65% to B. velezensis ONU-553, sharing 60 orthologous genes. Secondary metabolite gene cluster analysis of all strains showed that KS04AU has a mersacidin cluster gene, which is absent in the genome of the other strains. PHASTER analysis also showed KS04AU harboring two phages (Aeribacllus AP45 NC_048651 and Paenibacillus_Tripp NC_028930), which were also unique in comparison with the other strains. Analysis on anti-microbial resistance genes showed no difference in the genome of KS04AU to any of the other genomes, with the exception of B. amyloliquefaciens IT-45 which had one unique small multidrug-resistance antibiotic efflux-pump gene (qacJ). The CRISPR-Cas systems in the strains were also compared showing one CRISPR gene found only in KS04AU. Hydrolytic activity, antagonistic activity against phytopathogens (Fusarium oxysporum, Fusarium graminearum, Alternaria alternata and Pseudomonas syringae) and biocontrol against tomato foot and root rot experiments were carried out. B. velezensis KS04AU inhibits the growth of all phytopathogens tested, produces hydrolytic activity, and reduces Fusarium oxysporum f.sp. radicis-lycopersici (Forl) ZUM2407 lesions up to 46.02 ± 0.12%. The obtained results confirm B. velezensis KS04AU as a potential biocontrol strain for plant protection.