PanACoTA: a modular tool for massive microbial comparative genomics

The Spread of Antibiotic Resistance Is Driven by Plasmids Among the Fastest Evolving and of Broadest Host Range

Microorganisms endure novel challenges for which other microorganisms in other biomes may have already evolved solutions. This is the case of nosocomial bacteria under antibiotic therapy because antibiotics are of ancient natural origin and resistances to them have previously emerged in environmental bacteria. In such cases, the rate of adaptation crucially depends on the acquisition of genes by horizontal transfer of plasmids from distantly related bacteria in different biomes. We hypothesized that such processes should be driven by plasmids among the most mobile and evolvable. We confirmed these predictions by showing that plasmid species encoding antibiotic resistance are very mobile, have broad host ranges, while showing higher rates of homologous recombination and faster turnover of gene repertoires than the other plasmids. These characteristics remain outstanding when we remove resistance plasmids from our dataset, suggesting that antibiotic resistance genes are preferentially acquired and carried by plasmid species that are intrinsically very mobile and plastic. Evolvability and mobility facilitate the transfer of antibiotic resistance, and presumably of other phenotypes, across distant taxonomic groups and biomes. Hence, plasmid species, and possibly those of other mobile genetic elements, have differentiated and predictable roles in the spread of novel traits.

Adaptive genomic plasticity in large-genome, broad-host-range vibrio phages

The host range of a bacteriophage-the diversity of hosts it can infect-is central to understanding phage ecology and applications. Whereas most well-characterized phages have narrow host ranges, broad-host-range phages represent an intriguing component of marine ecosystems. The genetic and evolutionary mechanisms driving their generalism remain poorly understood. In this study, we analyzed Schizotequatroviruses and their Vibrio crassostreae hosts, collected from an oyster farm. Schizotequatroviruses exhibit broad host ranges, large genomes (~252 kbp) encoding 26 transfer ribonucleic acids, and conserved genomic organization interspersed with recombination hotspots. These recombination events, particularly in regions encoding receptor-binding proteins and antidefense systems, highlight their adaptability to host resistance. Some lineages demonstrated the ability of receptor-switching between OmpK and LamB. Despite their broad host range, Schizotequatroviruses were rare in the environment. Their scarcity could not be attributed to burst size, which was comparable to other phages in vitro, but may result from ecological constraints or fitness trade-offs, such as their preference for targeting generalist vibrios in seawater rather than the patho-phylotypes selected in oyster farms. Our findings clarify the genetic and ecological variables shaping Schizotequatrovirus generalism and provide a foundation for future phage applications in aquaculture and beyond.

Leptospira gorisiae sp. nov, L. cinconiae sp. nov, L. mgodei sp. nov, L. milleri sp. nov and L. iowaensis sp. nov: five new species isolated from water sources in the Midwestern United States

Isolates of Leptospira spp. were cultured from water sources at five different sites in central Iowa in the Midwestern United States and characterized by whole-genome sequencing. Isolates were helix-shaped and motile. Genome sequence analyses determined that the isolates could be clearly distinguished from other species described in the genus Leptospira and included one species that belonged to the pathogen subclade P1, one species that belonged to the pathogen subclade P2 and three species that belonged to the saprophyte subclade S1. The names Leptospira gorisiae sp. nov. (type strain WS92.C1T=NVSL-WS92.C1T=KIT0303T), Leptospira cinconiae sp. nov. (type strain WS58.C1T=NVSL-WS58.C1T=KIT0304T), Leptospira mgodei sp. nov. (type strain WS4.C2T=NVSL.WS4.C2T=KIT0305T), Leptospira iowaensis sp. nov. (type strain WS39.C2T=NVSL-WS39.C2T=KIT0306T) and Leptospira milleri sp. nov. (type strain WS60.C2T=NVSL-WS60.C2T=KIT0307T) are proposed.

Genomic Characterization of Probiotic Purple Nonsulfur Bacteria Cereibacter sphaeroides Strains S3W10 and SS15: Implications for Enhanced Shrimp Aquaculture

Cereibacter sphaeroides strains S3W10 and SS15, isolated from shrimp ponds, exhibit potential probiotic benefits for aquaculture. In this study, the genomic features of S3W10 and SS15 were thoroughly characterized to evaluate their probiotic properties and safety for aquaculture use. The genomes of S3W10 and SS15 consist of 130 and 74 contigs, with sizes of 4.6 Mb and 4.4 Mb and GC contents of 69.2%. Average nucleotide identity (ANI), digital DNA-DNA hybridization (dDDH), and phylogenomic analyses confirmed that these strains belong to C. sphaeroides. Genome annotation predicted 4260 coding sequences (CDS) in S3W10 and 4086 CDS in SS15, including genes associated with stress tolerance, nutrient absorption, and antioxidant activity. Notably, genes related to vitamin B12 synthesis, digestive enzyme production, and carotenoid biosynthesis, which support shrimp health, were identified in both genomes. CAZyme analysis identified 116 and 115 carbohydrate-active enzymes in S3W10 and SS15, respectively, supporting adaptation to gastrointestinal environments and the host immune response. Pan-genome analysis across C. sphaeroides strains revealed 7918 gene clusters, highlighting the open pan-genome structure of this species and its high genetic diversity. Further bioinformatic analyses assessing mobile genetic elements, antibiotic-resistance genes, and virulence factors demonstrated the safety of both strains for aquaculture, as no plasmids or virulence genes were identified. The genomic insights in this study provide a deeper understanding of the strains' adaptability and functional potential, aligning with previous in vitro and in vivo studies and highlighting their potential for use in shrimp cultivation.

Metabolic capabilities are highly conserved among human nasal-associated Corynebacterium species in pangenomic analyses

Corynebacterium species are globally ubiquitous in human nasal microbiota across the lifespan. Moreover, nasal microbiota profiles typified by higher relative abundances of Corynebacterium are often positively associated with health. Among the most common human nasal Corynebacterium species are C. propinquum, C. pseudodiphtheriticum, C. accolens, and C. tuberculostearicum. To gain insight into the functions of these four species, we identified genomic, phylogenomic, and pangenomic properties and estimated the metabolic capabilities of 87 distinct human nasal Corynebacterium strain genomes: 31 from Botswana and 56 from the United States. C. pseudodiphtheriticum had geographically distinct clades consistent with localized strain circulation, whereas some strains from the other species had wide geographic distribution spanning Africa and North America. All species had similar genomic and pangenomic structures. Gene clusters assigned to all COG metabolic categories were overrepresented in the persistent versus accessory genome of each species indicating limited strain-level variability in metabolic capacity. Based on prevalence data, at least two Corynebacterium species likely coexist in the nasal microbiota of 82% of adults. So, it was surprising that core metabolic capabilities were highly conserved among the four species indicating limited species-level metabolic variation. Strikingly, strains in the U.S. clade of C. pseudodiphtheriticum lacked genes for assimilatory sulfate reduction present in most of the strains in the Botswana clade and in the other studied species, indicating a recent, geographically related loss of assimilatory sulfate reduction. Overall, the minimal species and strain variability in metabolic capacity implies coexisting strains might have limited ability to occupy distinct metabolic niches.

IMPORTANCE

Pangenomic analysis with estimation of functional capabilities facilitates our understanding of the full biologic diversity of bacterial species. We performed systematic genomic, phylogenomic, and pangenomic analyses with qualitative estimation of the metabolic capabilities of four common human nasal Corynebacterium species, along with focused experimental validations, generating a foundational resource. The prevalence of each species in human nasal microbiota is consistent with the common coexistence of at least two species. We identified a notably high level of metabolic conservation within and among species indicating limited options for species to occupy distinct metabolic niches, highlighting the importance of investigating interactions among nasal Corynebacterium species. Comparing strains from two continents, C. pseudodiphtheriticum had restricted geographic strain distribution characterized by an evolutionarily recent loss of assimilatory sulfate reduction in U.S. strains. Our findings contribute to understanding the functions of Corynebacterium within human nasal microbiota and to evaluating their potential for future use as biotherapeutics.

Prediction of strain level phage-host interactions across the Escherichia genus using only genomic information

Predicting bacteriophage infection of specific bacterial strains promises advancements in phage therapy and microbial ecology. Whether the dynamics of well-established phage-host model systems generalize to the wide diversity of microbes is currently unknown. Here we show that we could accurately predict the outcomes of phage-bacteria interactions at the strain level in natural isolates from the genus Escherichia using only genomic data (area under the receiver operating characteristic curve (AUROC) of 86%). We experimentally established a dataset of interactions between 403 diverse Escherichia strains and 96 phages. Most interactions are explained by adsorption factors as opposed to antiphage systems which play a marginal role. We trained predictive algorithms and pinpoint poorly predicted interactions to direct future research efforts. Finally, we established a pipeline to recommend tailored phage cocktails, demonstrating efficiency on 100 pathogenic E. coli isolates. This work provides quantitative insights into phage-host specificity and supports the use of predictive algorithms in phage therapy.

Whole-Genome Sequencing of Three Lactiplantibacillus plantarum Strains Reveals Potential Metabolites for Boosting Host Immunity Safely

In response to the growing demand for immune-related products, this study evaluated the safety and immune-modulating potential of three newly discovered Lactiplantibacillus plantarum strains (GKM3, GKK1, and GKD7) through toxicity tests and whole-genome sequencing. Safety evaluations, including the analysis of antimicrobial resistance genes, virulence factors, plasmids, and prophages, classified these strains as safe for human consumption. Acute oral toxicity tests further supported their safety. To evaluate their immune-modulating potential, dendritic cells were exposed to these strains, and the secretion of key cytokines (IFN-β and IL-12) was measured. Among the strains, GKK1 exhibited the highest enhancement of IFN-β and IL-12 production, suggesting its potential as an immune-stimulating probiotic. Bioinformatics analysis revealed potential metabolic pathways and secondary metabolites, including predicted bacteriocins, associated with immune modulation. The presence of a nitrate reductase region in the GKK1 strain indicated its ability to produce nitric oxide, a critical molecule involved in immune regulation and host defense. The presence of glucorhamnan-related gene clusters in GKK1 also suggested immune-enhancing effects. Nitrate reductase expression was confirmed using qPCR, with the highest levels detected in GKK1. Moreover, this study is the first to show an anti-inflammatory effect of plantaricin A, linked to its presence in strain GKM3 and its potential therapeutic applications due to sequence similarity to known anti-inflammatory peptides. Overall, these three L. plantarum strains demonstrated a safe profile and GKK1 showed potential as an immunity-enhancing probiotic. However, additional investigation is required to confirm the involvement of specific metabolic pathways, secondary metabolites, and bacteriocins in immune responses.

C. Rocha E. Intragenomic conflicts with plasmids and chromosomal mobile genetic elements drive the evolution of natural transformation within species

Natural transformation is the only mechanism of genetic exchange controlled by the recipient bacteria. We quantified its rates in 786 clinical strains of the human pathogens Legionella pneumophila (Lp) and 496 clinical and environmental strains of Acinetobacter baumannii (Ab). The analysis of transformation rates in the light of phylogeny revealed they evolve by a mixture of frequent small changes and a few large quick jumps across 6 orders of magnitude. In standard conditions close to half of the strains of Lp and a more than a third in Ab are below the detection limit and thus presumably non-transformable. Ab environmental strains tend to have higher transformation rates than the clinical ones. Transitions to non-transformability were frequent and usually recent, suggesting that they are deleterious and subsequently purged by natural selection. Accordingly, we find that transformation decreases genetic linkage in both species, which might accelerate adaptation. Intragenomic conflicts with chromosomal mobile genetic elements (MGEs) and plasmids could explain these transitions and a GWAS confirmed systematic negative associations between transformation and MGEs: plasmids and other conjugative elements in Lp, prophages in Ab, and transposable elements in both. In accordance with the hypothesis of modulation of transformation rates by genetic conflicts, transformable strains have fewer MGEs in both species and some MGEs inactivate genes implicated in the transformation with heterologous DNA (in Ab). Innate defense systems against MGEs are associated with lower transformation rates, especially restriction-modification systems. In contrast, CRISPR-Cas systems are associated with higher transformation rates suggesting that adaptive defense systems may facilitate cell protection from MGEs while preserving genetic exchanges by natural transformation. Ab and Lp have different lifestyles, gene repertoires, and population structure. Nevertheless, they exhibit similar trends in terms of variation of transformation rates and its determinants, suggesting that genetic conflicts could drive the evolution of natural transformation in many bacteria.

Proposal of Lactobacillus amylovorus subsp. animalis subsp. nov. and an emended description of Lactobacillus amylovorus

Seven novel lactic acid bacterial strains (BF125T, BF186, TKL145, YK3, YK6, YK10 and NSK) were isolated from the fresh faeces of Japanese black beef cattle and weanling piglets, spent mushroom substrates, or steeping water of a corn starch production plant. These strains are rod-shaped, Gram-stain-positive, non-motile, non-spore-forming, catalase-negative, cytochrome oxidase-negative, facultatively anaerobic, and homofermentative. Strain BF125T did not produce any gas from glucose; both d- and l-lactate were produced as end-products of glucose (D/L, 40 : 60). Growth occurred at 30-45 °C (optimum, 37 °C), pH 5.0-8.0 (optimum, pH 6.0), and with NaCl concentration of 1.0-3.0% (w/v). The G+C content of genomic DNA of strain BF125T was 37.8 mol% (whole-genome analysis). The major fatty acids were C16 : 0, C18 : 1 ω9c, C19 cyclopropane 9, 10, and summed feature 10. The 16S rRNA gene in strain BF125T showed high similarity to that of the type strain of Lactobacillus amylovorus (99.93%), and the other isolates were also identified as L. amylovorus based on these similarities. A phylogenetic tree based on the core genomes of L. amylovorus strains (n=54), including the seven isolates, showed that they could be divided into two clusters. Strains YK3, YK6, YK10, and NSK were in the first cluster, along with the type strain DSM 20531T, while the second cluster included isolates BF125T, BF186, TKL145, and other strains isolated from various animal origins. Phenotypic differences in fermentability were observed for lactose, salicin, and gentiobiose between these two groups. The intergroup digital DNA-DNA hybridization values (72.9-78.6%) and intergroup average nucleotide identity values (95.64-96.92%) were comparable to values calculated using datasets of other valid subspecies of the genus (ex-) Lactobacillus. In light of the physiological, genotypic, and phylogenetic evidence, we propose a novel subspecies of L. amylovorus, named Lactobacillus amylovorus subsp. animalis subsp. nov. (type strain BF125T=MAFF 212522T=DSM 115528T). Our findings also led to the automatic creation of Lactobacillus amylovorus subsp. amylovorus subsp. nov. and an emended description of the species L. amylovorus.

Integrated OMICs approach reveals energy metabolism pathway is vital for Salmonella Pullorum survival within the egg white

Eggs, an important part of a healthy daily diet, can protect chicken embryo development due to the shell barrier and various antibacterial components within the egg white. Our previous study demonstrated that Salmonella Pullorum, highly adapted to chickens, can survive in the egg white and, therefore, be passed to newly hatched chicks. However, the survival strategy of Salmonella Pullorum in antibacterial conditions remains unknown. The overall transcripts in the egg white showed a large-scale shift compared to LB broth. The expression of common response genes and pathways, such as those involved in iron uptake, biotin biosynthesis, and virulence, was significantly changed, consistent with the other transovarial transmission serovar Enteritidis. Notably, membrane stress response, amino acid metabolism, and carbohydrate metabolism were specifically affected. Additional upregulated functionally relevant genes (JI728_13095, JI728_13100, JI728_17960, JI728_10085, JI728_15605, and nhaA) as mutants confirmed the susceptible phenotype. Furthermore, fim deletion resulted in an increased survival capacity in the egg white, consistent with the downregulated expression. The second-round RNA-Seq analysis of the Δfim mutant in the egg white revealed significantly upregulated genes compared with the wild type in the egg white responsible for energy metabolism located on the hyc and hyp operons regulated by FhlA, indicating the Δfim mutant cannot receive enough oxygen and switched to fermentative growth due to its inability to attach to the albumen surface. Together, this study provides a first estimate of the global transcriptional response of Salmonella Pullorum under antibacterial egg white and highlights the new potential role of fim deletion in optimizing energy metabolism pathways that may assist vertical transmission.

IMPORTANCE

Pullorum disease, causing serious embryo death and chick mortality, results in substantial economic losses worldwide due to transovarial transmission. Egg-borne outbreaks are frequently reported in many countries. The present study has filled the knowledge gap regarding how the specific chicken-adapted pathogen Salmonella Pullorum behaves within the challenging environment of egg white. The deletion of the fim fimbrial system can increase survival in the albumen, possibly by reprogramming metabolism-related gene products, which reveals a new adaptive strategy of pathogens. Moreover, the comparison, including previous research on Salmonella Enteritidis, capable of vertical transmission, aims to provide diversified data sets in the field and further help to implement reasonable and effective measures to improve both food safety and animal health.

Lacrimispora brassicae sp. nov. isolated from fermented cabbage, and proposal of Clostridium indicum Gundawar et al. 2019 and Clostridium methoxybenzovorans Mechichi et al. 1999 as heterotypic synonyms of Lacrimispora amygdalina (Parshina et al. 2003) Haas and Blanchard 2020 and Lacrimispora indolis (McClung and McCoy 1957) Haas and Blanchard 2020, respectively

A Gram-stain-negative, endospore-forming, rod-shaped, indole-producing bacterial strain, designated YZC6T, was isolated from fermented cabbage. Strain YZC6T grew at 10-37  °C, pH 5.5-8.5, and with up to 2  % (w/v) NaCl. The major cellular fatty acids were C16 : 0 and C18 : 1 cis 11 dimethyl acetal. Phylogenetic analysis of the 16S rRNA gene revealed that strain YZC6T belonged to the genus Lacrimispora and was closely related to Lacrimispora aerotolerans DSM 5434T (98.3  % sequence similarity), Lacrimispora saccharolytica WM1T (98.1  %), and Lacrimispora algidixylanolytica SPL73T (98.1  %). The average nucleotide identity based on blast (below 87.8  %) and digital DNA-DNA hybridization (below 36.1 %) values between the novel isolate and its corresponding relatives showed that strain YZC6T could be readily distinguished from its closely related species. Based on genotypic, phenotypic, and chemotaxonomic data, a novel Lacrimispora species, Lacrimispora brassicae sp. nov., was proposed, with YZC6T as the type strain (=MAFF 212518T=JCM 32810T=DSM 112100T). This study also proposed Clostridium indicum Gundawar et al. 2019 as a later heterotypic synonym of Lacrimispora amygdalina (Parshina et al. 2003) Haas and Blanchard 2020 and Clostridium methoxybenzovorans Mechichi et al. 1999 as a later heterotypic synonym of Lacrimispora indolis (McClung and McCpy 1957) Haas and Blanchard 2020.

Capsules and their traits shape phage susceptibility and plasmid conjugation efficiency

Bacterial evolution is affected by mobile genetic elements like phages and conjugative plasmids, offering new adaptive traits while incurring fitness costs. Their infection is affected by the bacterial capsule. Yet, its importance has been difficult to quantify because of the high diversity of confounding mechanisms in bacterial genomes such as anti-viral systems and surface receptor modifications. Swapping capsule loci between Klebsiella pneumoniae strains allowed us to quantify their impact on plasmid and phage infection independently of genetic background. Capsule swaps systematically invert phage susceptibility, revealing serotypes as key determinants of phage infection. Capsule types also influence conjugation efficiency in both donor and recipient cells, a mechanism shaped by capsule volume and conjugative pilus structure. Comparative genomics confirmed that more permissive serotypes in the lab correspond to the strains acquiring more conjugative plasmids in nature. The least capsule-sensitive pili (F-like) are the most frequent in the species' plasmids, and are the only ones associated with both antibiotic resistance and virulence factors, driving the convergence between virulence and antibiotics resistance in the population. These results show how traits of cellular envelopes define slow and fast lanes of infection by mobile genetic elements, with implications for population dynamics and horizontal gene transfer.

Evolution and emergence of Mycobacterium tuberculosis

Tuberculosis (TB) remains one of the deadliest infectious diseases in human history, prevailing even in the 21st century. The causative agents of TB are represented by a group of closely related bacteria belonging to the Mycobacterium tuberculosis complex (MTBC), which can be subdivided into several lineages of human- and animal-adapted strains, thought to have shared a last common ancestor emerged by clonal expansion from a pool of recombinogenic Mycobacterium canettii-like tubercle bacilli. A better understanding of how MTBC populations evolved from less virulent mycobacteria may allow for discovering improved TB control strategies and future epidemiologic trends. In this review, we highlight new insights into the evolution of mycobacteria at the genus level, describing different milestones in the evolution of mycobacteria, with a focus on the genomic events that have likely enabled the emergence and the dominance of the MTBC. We also review the recent literature describing the various MTBC lineages and highlight their particularities and differences with a focus on host preferences and geographic distribution. Finally, we discuss on putative mechanisms driving the evolution of tubercle bacilli and mycobacteria in general, by taking the mycobacteria-specific distributive conjugal transfer as an example.

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.

Phage-inducible chromosomal minimalist islands (PICMIs), a novel family of small marine satellites of virulent phages

Phage satellites are bacterial genetic elements that co-opt phage machinery for their own dissemination. Here we identify a family of satellites, named Phage-Inducible Chromosomal Minimalist Islands (PICMIs), that are broadly distributed in marine bacteria of the family Vibrionaceae. A typical PICMI is characterized by reduced gene content, does not encode genes for capsid remodelling, and packages its DNA as a concatemer. PICMIs integrate in the bacterial host genome next to the fis regulator, and encode three core proteins necessary for excision and replication. PICMIs are dependent on virulent phage particles to spread to other bacteria, and protect their hosts from other competitive phages without interfering with their helper phage. Thus, our work broadens our understanding of phage satellites and narrows down the minimal number of functions necessary to hijack a tailed phage.

Assembly of a unique membrane complex in type VI secretion systems of Bacteroidota

The type VI secretion system (T6SS) of Gram-negative bacteria inhibits competitor cells through contact-dependent translocation of toxic effector proteins. In Proteobacteria, the T6SS is anchored to the cell envelope through a megadalton-sized membrane complex (MC). However, the genomes of Bacteroidota with T6SSs appear to lack genes encoding homologs of canonical MC components. Here, we identify five genes in Bacteroides fragilis (tssNQOPR) that are essential for T6SS function and encode a Bacteroidota-specific MC. We purify this complex, reveal its dimensions using electron microscopy, and identify a protein-protein interaction network underlying the assembly of the MC including the stoichiometry of the five TssNQOPR components. Protein TssN mediates the connection between the Bacteroidota MC and the conserved baseplate. Although MC gene content and organization varies across the phylum Bacteroidota, no MC homologs are detected outside of T6SS loci, suggesting ancient co-option and functional convergence with the non-homologous MC of Pseudomonadota.

Machine learning and phylogenetic analysis allow for predicting antibiotic resistance in M. tuberculosis

BACKGROUND

Antimicrobial resistance (AMR) poses a significant global health threat, and an accurate prediction of bacterial resistance patterns is critical for effective treatment and control strategies. In recent years, machine learning (ML) approaches have emerged as powerful tools for analyzing large-scale bacterial AMR data. However, ML methods often ignore evolutionary relationships among bacterial strains, which can greatly impact performance of the ML methods, especially if resistance-associated features are attempted to be detected. Genome-wide association studies (GWAS) methods like linear mixed models accounts for the evolutionary relationships in bacteria, but they uncover only highly significant variants which have already been reported in literature.

RESULTS

In this work, we introduce a novel phylogeny-related parallelism score (PRPS), which measures whether a certain feature is correlated with the population structure of a set of samples. We demonstrate that PRPS can be used, in combination with SVM- and random forest-based models, to reduce the number of features in the analysis, while simultaneously increasing models' performance. We applied our pipeline to publicly available AMR data from PATRIC database for Mycobacterium tuberculosis against six common antibiotics.

CONCLUSIONS

Using our pipeline, we re-discovered known resistance-associated mutations as well as new candidate mutations which can be related to resistance and not previously reported in the literature. We demonstrated that taking into account phylogenetic relationships not only improves the model performance, but also yields more biologically relevant predicted most contributing resistance markers.

Comparison of gene clustering criteria reveals intrinsic uncertainty in pangenome analyses

BACKGROUND

A key step for comparative genomics is to group open reading frames into functionally and evolutionarily meaningful gene clusters. Gene clustering is complicated by intraspecific duplications and horizontal gene transfers that are frequent in prokaryotes. In consequence, gene clustering methods must deal with a trade-off between identifying vertically transmitted representatives of multicopy gene families, which are recognizable by synteny conservation, and retrieving complete sets of species-level orthologs. We studied the implications of adopting homology, orthology, or synteny conservation as formal criteria for gene clustering by performing comparative analyses of 125 prokaryotic pangenomes.

RESULTS

Clustering criteria affect pangenome functional characterization, core genome inference, and reconstruction of ancestral gene content to different extents. Species-wise estimates of pangenome and core genome sizes change by the same factor when using different clustering criteria, allowing robust cross-species comparisons regardless of the clustering criterion. However, cross-species comparisons of genome plasticity and functional profiles are substantially affected by inconsistencies among clustering criteria. Such inconsistencies are driven not only by mobile genetic elements, but also by genes involved in defense, secondary metabolism, and other accessory functions. In some pangenome features, the variability attributed to methodological inconsistencies can even exceed the effect sizes of ecological and phylogenetic variables.

CONCLUSIONS

Choosing an appropriate criterion for gene clustering is critical to conduct unbiased pangenome analyses. We provide practical guidelines to choose the right method depending on the research goals and the quality of genome assemblies, and a benchmarking dataset to assess the robustness and reproducibility of future comparative studies.

Genomic characterization of the antiviral arsenal of Actinobacteria

Phages are ubiquitous in nature, and bacteria with very different genomics, metabolisms, and lifestyles are subjected to their predation. Yet, the defence systems that allow bacteria to resist their phages have rarely been explored experimentally outside a very limited number of model organisms. Actinobacteria (Actinomycetota) are a phylum of GC-rich Gram-positive bacteria, which often produce an important diversity of secondary metabolites. Despite being ubiquitous in a wide range of environments, from soil to fresh and sea water but also the gut microbiome, relatively little is known about the anti-phage arsenal of Actinobacteria. In this work, we used DefenseFinder to systematically detect 131 anti-phage defence systems in 22803 fully sequenced prokaryotic genomes, among which are 2253 Actinobacteria of more than 700 species. We show that, like other bacteria, Actinobacteria encode many diverse anti-phage systems that are often encoded on mobile genetic elements. We further demonstrate that most detected defence systems are absent or rarer in Actinobacteria than in other bacteria, while a few rare systems are enriched (notably gp29-gp30 and Wadjet). We characterize the spatial distribution of anti-phage systems on Streptomyces chromosomes and show that some defence systems (e.g. RM systems) tend to be encoded in the core region, while others (e.g. Lamassu and Wadjet) are enriched towards the extremities. Overall, our results suggest that Actinobacteria might be a source of novel anti-phage systems and provide clues to characterize mechanistic aspects of known anti-phage systems.

Identification and characterization of thousands of bacteriophage satellites across bacteria

Bacteriophage-bacteria interactions are affected by phage satellites, elements that exploit phages for transfer between bacteria. Satellites can encode defense systems, antibiotic resistance genes, and virulence factors, but their number and diversity are unknown. We developed SatelliteFinder to identify satellites in bacterial genomes, detecting the four best described families: P4-like, phage inducible chromosomal islands (PICI), capsid-forming PICI, and PICI-like elements (PLE). We vastly expanded the number of described elements to ∼5000, finding bacterial genomes with up to three different families of satellites. Most satellites were found in Proteobacteria and Firmicutes, but some are in novel taxa such as Actinobacteria. We characterized the gene repertoires of satellites, which are variable in size and composition, and their genomic organization, which is very conserved. Phylogenies of core genes in PICI and cfPICI indicate independent evolution of their hijacking modules. There are few other homologous core genes between other families of satellites, and even fewer homologous to phages. Hence, phage satellites are ancient, diverse, and probably evolved multiple times independently. Given the many bacteria infected by phages that still lack known satellites, and the recent proposals for novel families, we speculate that we are at the beginning of the discovery of massive numbers and types of satellites.

PlantTribes2: Tools for comparative gene family analysis in plant genomics

Plant genome-scale resources are being generated at an increasing rate as sequencing technologies continue to improve and raw data costs continue to fall; however, the cost of downstream analyses remains large. This has resulted in a considerable range of genome assembly and annotation qualities across plant genomes due to their varying sizes, complexity, and the technology used for the assembly and annotation. To effectively work across genomes, researchers increasingly rely on comparative genomic approaches that integrate across plant community resources and data types. Such efforts have aided the genome annotation process and yielded novel insights into the evolutionary history of genomes and gene families, including complex non-model organisms. The essential tools to achieve these insights rely on gene family analysis at a genome-scale, but they are not well integrated for rapid analysis of new data, and the learning curve can be steep. Here we present PlantTribes2, a scalable, easily accessible, highly customizable, and broadly applicable gene family analysis framework with multiple entry points including user provided data. It uses objective classifications of annotated protein sequences from existing, high-quality plant genomes for comparative and evolutionary studies. PlantTribes2 can improve transcript models and then sort them, either genome-scale annotations or individual gene coding sequences, into pre-computed orthologous gene family clusters with rich functional annotation information. Then, for gene families of interest, PlantTribes2 performs downstream analyses and customizable visualizations including, (1) multiple sequence alignment, (2) gene family phylogeny, (3) estimation of synonymous and non-synonymous substitution rates among homologous sequences, and (4) inference of large-scale duplication events. We give examples of PlantTribes2 applications in functional genomic studies of economically important plant families, namely transcriptomics in the weedy Orobanchaceae and a core orthogroup analysis (CROG) in Rosaceae. PlantTribes2 is freely available for use within the main public Galaxy instance and can be downloaded from GitHub or Bioconda. Importantly, PlantTribes2 can be readily adapted for use with genomic and transcriptomic data from any kind of organism.

Role of Klebsiella pneumoniae Type VI secretion system (T6SS) in long-term gastrointestinal colonization

Type VI secretion systems (T6SS), recently described in hypervirulent K. pneumoniae (hvKp) strains, are involved in bacterial warfare but their role in classical clinical strains (cKp) has been little investigated. In silico analysis indicated the presence of T6SS clusters (from zero to four), irrespective of the strains origin or virulence, with a high prevalence in the K. pneumoniae species (98%). In the strain CH1157, two T6SS-apparented pathogenicity islands were detected, T6SS-1 and -2, harboring a phospholipase-encoding gene (tle1) and a potential new effector-encoding gene named tke (Type VI Klebsiella effector). Tle1 expression in Escherichia coli periplasm affected cell membrane permeability. T6SS-1 isogenic mutants colonized the highest gastrointestinal tract of mice less efficiently than their parental strain, at long term. Comparative analysis of faecal 16S sequences indicated that T6SS-1 impaired the microbiota richness and its resilience capacity. Oscillospiraceae family members could be specific competitors for the long-term gut establishment of K. pneumoniae.

Adaptation to novel spatially-structured environments is driven by the capsule and alters virulence-associated traits

The extracellular capsule is a major virulence factor, but its ubiquity in free-living bacteria with large environmental breadths suggests that it shapes adaptation to novel niches. Yet, how it does so, remains unexplored. Here, we evolve three Klebsiella strains and their capsule mutants in parallel. Their comparison reveals different phenotypic and genotypic evolutionary changes that alter virulence-associated traits. Non-capsulated populations accumulate mutations that reduce exopolysaccharide production and increase biofilm formation and yield, whereas most capsulated populations become hypermucoviscous, a signature of hypervirulence. Hence, adaptation to novel environments primarily occurs by fine-tuning expression of the capsular locus. The same evolutionary conditions selecting for mutations in the capsular gene wzc leading to hypermucoviscosity also result in increased susceptibility to antibiotics by mutations in the ramA regulon. This implies that general adaptive processes outside the host can affect capsule evolution and its role in virulence and infection outcomes may be a by-product of such adaptation.

Phenotypic and genotypic evolution in worrisome Klebsiella spp. is influenced by the capsule. Here the authors show that adaptation outside the host can impact virulence-associated traits, including de novo emergence of hypermucoviscosity.

Intestinal microbial communities and Holdemanella isolated from HIV+/- men who have sex with men increase frequencies of lamina propria CCR5+ CD4+ T cells

Men who have sex with men (MSM), regardless of HIV infection status, have an intestinal microbiome that is compositionally distinct from men who have sex with women (MSW) and women. We recently showed HIV-negative MSM have elevated levels of intestinal CD4⁺ T cells expressing CCR5, a critical co-receptor for HIV. Whether elevated expression of CCR5 is driven by the altered gut microbiome composition in MSM has not been explored. Here we used in vitro stimulation of gut Lamina Propria Mononuclear Cells (LPMCs) with whole intact microbial cells isolated from stool to demonstrate that fecal bacterial communities (FBCs) from HIV-positive/negative MSM induced higher frequencies of CCR5⁺ CD4⁺ T cells compared to FBCs from HIV-negative MSW and women. To identify potential microbial drivers, we related the frequency of CCR5⁺ CD4⁺ T cells to the abundance of individual microbial taxa in rectal biopsy of HIV-positive/negative MSM and controls, and Holdemanella biformis was strongly associated with increased frequency of CCR5⁺ CD4⁺ T cells. We used in vitro stimulation of gut LPMCs with the type strain of H. biformis, a second strain of H.biformis and an isolate of the closely related Holdemanella porci , cultured from either a HIV-positive or a HIV-negative MSM stool. H. porci elevated the frequency of both CCR5⁺ CD4⁺ T cells and the ratio of TNF-α/IL-10 Genomic comparisons of the 3 Holdemanella isolates revealed unique cell wall and capsular components, which may be responsible for their differences in immunogenicity. These findings describe a novel mechanism potentially linking intestinal dysbiosis in MSM to HIV transmission and mucosal pathogenesis.

To catch a hijacker: abundance, evolution and genetic diversity of P4-like bacteriophage satellites

Bacteriophages (phages) are bacterial parasites that can themselves be parasitized by phage satellites. The molecular mechanisms used by satellites to hijack phages are sometimes understood in great detail, but the origins, abundance, distribution and composition of these elements are poorly known. Here, we show that P4-like elements are present in more than 30% of the genomes of Enterobacterales, and in almost half of those of Escherichia coli, sometimes in multiple distinct copies. We identified over 1000 P4-like elements with very conserved genetic organization of the core genome and a few hotspots with highly variable genes. These elements are never found in plasmids and have very little homology to known phages, suggesting an independent evolutionary origin. Instead, they are scattered across chromosomes, possibly because their integrases are often exchanged with other elements. The rooted phylogenies of hijacking functions are correlated and suggest longstanding coevolution. They also reveal broad host ranges in P4-like elements, as almost identical elements can be found in distinct bacterial genera. Our results show that P4-like phage satellites constitute a very distinct, widespread and ancient family of mobile genetic elements. They pave the way for studying the molecular evolution of antagonistic interactions between phages and their satellites. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

Genomic Stability and Genetic Defense Systems in Dolosigranulum pigrum, a Candidate Beneficial Bacterium from the Human Microbiome

Dolosigranulum pigrum is positively associated with indicators of health in multiple epidemiological studies of human nasal microbiota. Knowledge of the basic biology of D. pigrum is a prerequisite for evaluating its potential for future therapeutic use; however, such data are very limited. To gain insight into D. pigrum's chromosomal structure, pangenome, and genomic stability, we compared the genomes of 28 D. pigrum strains that were collected across 20 years. Phylogenomic analysis showed closely related strains circulating over this period and closure of 19 genomes revealed highly conserved chromosomal synteny. Gene clusters involved in the mobilome and in defense against mobile genetic elements (MGEs) were enriched in the accessory genome versus the core genome. A systematic analysis for MGEs identified the first candidate D. pigrum prophage and insertion sequence. A systematic analysis for genetic elements that limit the spread of MGEs, including restriction modification (RM), CRISPR-Cas, and deity-named defense systems, revealed strain-level diversity in host defense systems that localized to specific genomic sites, including one RM system hot spot. Analysis of CRISPR spacers pointed to a wealth of MGEs against which D. pigrum defends itself. These results reveal a role for horizontal gene transfer and mobile genetic elements in strain diversification while highlighting that in D. pigrum this occurs within the context of a highly stable chromosomal organization protected by a variety of defense mechanisms. IMPORTANCE Dolosigranulum pigrum is a candidate beneficial bacterium with potential for future therapeutic use. This is based on its positive associations with characteristics of health in multiple studies of human nasal microbiota across the span of human life. For example, high levels of D. pigrum nasal colonization in adults predicts the absence of Staphylococcus aureus nasal colonization. Also, D. pigrum nasal colonization in young children is associated with healthy control groups in studies of middle ear infections. Our analysis of 28 genomes revealed a remarkable stability of D. pigrum strains colonizing people in the United States across a 20-year span. We subsequently identified factors that can influence this stability, including genomic stability, phage predators, the role of MGEs in strain-level variation, and defenses against MGEs. Finally, these D. pigrum strains also lacked predicted virulence factors. Overall, these findings add additional support to the potential for D. pigrum as a therapeutic bacterium.

Interplay between the cell envelope and mobile genetic elements shapes gene flow in populations of the nosocomial pathogen Klebsiella pneumoniae

Mobile genetic elements (MGEs) drive genetic transfers between bacteria using mechanisms that require a physical interaction with the cellular envelope. In the high-priority multidrug-resistant nosocomial pathogens (ESKAPE), the first point of contact between the cell and virions or conjugative pili is the capsule. While the capsule can be a barrier to MGEs, it also evolves rapidly by horizontal gene transfer (HGT). Here, we aim at understanding this apparent contradiction by studying the covariation between the repertoire of capsule genes and MGEs in approximately 4,000 genomes of Klebsiella pneumoniae (Kpn). We show that capsules drive phage-mediated gene flow between closely related serotypes. Such serotype-specific phage predation also explains the frequent inactivation of capsule genes, observed in more than 3% of the genomes. Inactivation is strongly epistatic, recapitulating the capsule biosynthetic pathway. We show that conjugative plasmids are acquired at higher rates in natural isolates lacking a functional capsular locus and confirmed experimentally this result in capsule mutants. This suggests that capsule inactivation by phage pressure facilitates its subsequent reacquisition by conjugation. Accordingly, capsule reacquisition leaves long recombination tracts around the capsular locus. The loss and regain process rewires gene flow toward other lineages whenever it leads to serotype swaps. Such changes happen preferentially between chemically related serotypes, hinting that the fitness of serotype-swapped strains depends on the host genetic background. These results enlighten the bases of trade-offs between the evolution of virulence and multidrug resistance and caution that some alternatives to antibiotics by selecting for capsule inactivation may facilitate the acquisition of antibiotic resistance genes (ARGs).

A study of how the complex interaction between capsules and mobile genetic elements shapes gene flow in populations of Klebsiella pneumoniae reveals that capsule inactivation by phage pressure facilitates its subsequent re-acquisition by conjugation, and this loss and re-gain process influences the gene flow towards other lineages whenever it leads to serotype changes.