Reference-free structural variant detection in microbiomes via long-read co-assembly graphs.

MOTIVATION: The study of bacterial genome dynamics is vital for understanding the mechanisms underlying microbial adaptation, growth, and their impact on host phenotype. Structural variants (SVs), genomic alterations of 50 base pairs or more, play a pivotal role in driving evolutionary processes and maintaining genomic heterogeneity within bacterial populations. While SV detection in isolate genomes is relatively straightforward, metagenomes present broader challenges due to the absence of clear reference genomes and the presence of mixed strains. In response, our proposed method rhea, forgoes reference genomes and metagenome-assembled genomes (MAGs) by encompassing all metagenomic samples in a series (time or other metric) into a single co-assembly graph. The log fold change in graph coverage between successive samples is then calculated to call SVs that are thriving or declining. RESULTS: We show rhea to outperform existing methods for SV and horizontal gene transfer (HGT) detection in two simulated mock metagenomes, particularly as the simulated reads diverge from reference genomes and an increase in strain diversity is incorporated. We additionally demonstrate use cases for rhea on series metagenomic data of environmental and fermented food microbiomes to detect specific sequence alterations between successive time and temperature samples, suggesting host advantage. Our approach leverages previous work in assembly graph structural and coverage patterns to provide versatility in studying SVs across diverse and poorly characterized microbial communities for more comprehensive insights into microbial gene flux. AVAILABILITY AND IMPLEMENTATION: rhea is open source and available at: https://github.com/treangenlab/rhea.

Type IV-A3 CRISPR-Cas systems drive inter-plasmid conflicts by acquiring spacers in trans.

Plasmid-encoded type IV-A CRISPR-Cas systems lack an acquisition module, feature a DinG helicase instead of a nuclease, and form ribonucleoprotein complexes of unknown biological functions. Type IV-A3 systems are carried by conjugative plasmids that often harbor antibiotic-resistance genes and their CRISPR array contents suggest a role in mediating inter-plasmid conflicts, but this function remains unexplored. Here, we demonstrate that a plasmid-encoded type IV-A3 system co-opts the type I-E adaptation machinery from its host, Klebsiella pneumoniae (K. pneumoniae), to update its CRISPR array. Furthermore, we reveal that robust interference of conjugative plasmids and phages is elicited through CRISPR RNA-dependent transcriptional repression. By silencing plasmid core functions, type IV-A3 impacts the horizontal transfer and stability of targeted plasmids, supporting its role in plasmid competition. Our findings shed light on the mechanisms and ecological function of type IV-A3 systems and demonstrate their practical efficacy for countering antibiotic resistance in clinically relevant strains.

Spread of carbapenemase-producing Morganella spp from 2013 to 2021: a comparative genomic study.

BACKGROUND: Morganella spp are opportunistic pathogens involved in various infections. Intrinsic resistance to multiple antibiotics (including colistin) combined with the emergence of carbapenemase producers reduces the number of active antimicrobials. The aim of this study was to characterise genetic features related to the spread of carbapenem-resistant Morganella spp. METHODS: This comparative genomic study included extensively drug-resistant Morganella spp isolates collected between Jan 1, 2013, and March 1, 2021, by the French National Reference Center (NRC; n=68) and European antimicrobial resistance reference centres in seven European countries (n=104), as well as one isolate from Canada, two reference strains from the Pasteur Institute collection (Paris, France), and two colistin-susceptible isolates from Bicêtre Hospital (Kremlin-Bicêtre, France). The isolates were characterised by whole-genome sequencing, antimicrobial susceptibility testing, and biochemical tests. Complete genomes from GenBank (n=103) were also included for genomic analysis, including phylogeny and determination of core genomes and resistomes. Genetic distance between different species or subspecies was performed using average nucleotide identity (ANI). Intrinsic resistance mechanisms to polymyxins were investigated by combining genetic analysis with mass spectrometry on lipid A. FINDINGS: Distance analysis by ANI of 275 isolates identified three groups: Morganella psychrotolerans, Morganella morganii subspecies sibonii, and M morganii subspecies morganii, and a core genome maximum likelihood phylogenetic tree showed that the M morganii isolates can be separated into four subpopulations. On the basis of these findings and of phenotypic divergences between isolates, we propose a modified taxonomy for the Morganella genus including four species, Morganella psychrotolerans, Morganella sibonii, Morganella morganii, and a new species represented by a unique environmental isolate. We propose that M morganii include two subspecies: M morganii subspecies morganii (the most prevalent) and M morganii subspecies intermedius. This modified taxonomy was supported by a difference in intrinsic resistance to tetracycline and conservation of metabolic pathways such as trehalose assimilation, both only present in M sibonii. Carbapenemase producers were mostly identified among five high-risk clones of M morganii subspecies morganii. The most prevalent carbapenemase corresponded to NDM-1, followed by KPC-2, and OXA-48. A cefepime-zidebactam combination was the most potent antimicrobial against the 172 extensively drug-resistant Morganella spp isolates in our collection from different European countries, which includes metallo-ß-lactamase producers. Lipid A analysis showed that the intrinsic resistance to colistin was associated with the presence of L-ARA4N on lipid A. INTERPRETATION: This global characterisation of, to our knowledge, the widest collection of extensively drug-resistant Morganella spp highlights the need to clarify the taxonomy and decipher intrinsic resistance mechanisms, and paves the way for further genomic comparisons. FUNDING: None.

Belt and braces: Two escape ways to maintain the cassette reservoir of large chromosomal integrons.

Integrons are adaptive devices that capture, stockpile, shuffle and express gene cassettes thereby sampling combinatorial phenotypic diversity. Some integrons called sedentary chromosomal integrons (SCIs) can be massive structures containing hundreds of cassettes. Since most of these cassettes are non-expressed, it is not clear how they remain stable over long evolutionary timescales. Recently, it was found that the experimental inversion of the SCI of Vibrio cholerae led to a dramatic increase of the cassette excision rate associated with a fitness defect. Here, we question the evolutionary sustainability of this apparently counter selected genetic context. Through experimental evolution, we find that the integrase is rapidly inactivated and that the inverted SCI can recover its original orientation by homologous recombination between two insertion sequences (ISs) present in the array. These two outcomes of SCI inversion restore the normal growth and prevent the loss of cassettes, enabling SCIs to retain their roles as reservoirs of functions. These results illustrate a nice interplay between gene orientation, genome rearrangement, bacterial fitness and demonstrate how integrons can benefit from their embedded ISs.

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.

Perturbation and resilience of the gut microbiome up to 3 months after ß-lactams exposure in healthy volunteers suggest an important role of microbial ß-lactamases.

BACKGROUND: Antibiotics notoriously perturb the gut microbiota. We treated healthy volunteers either with cefotaxime or ceftriaxone for 3 days, and collected in each subject 12 faecal samples up to day 90. Using untargeted and targeted phenotypic and genotypic approaches, we studied the changes in the bacterial, phage and fungal components of the microbiota as well as the metabolome and the ß-lactamase activity of the stools. This allowed assessing their degrees of perturbation and resilience. RESULTS: While only two subjects had detectable concentrations of antibiotics in their faeces, suggesting important antibiotic degradation in the gut, the intravenous treatment perturbed very significantly the bacterial and phage microbiota, as well as the composition of the metabolome. In contrast, treatment impact was relatively low on the fungal microbiota. At the end of the surveillance period, we found evidence of resilience across the gut system since most components returned to a state like the initial one, even if the structure of the bacterial microbiota changed and the dynamics of the different components over time were rarely correlated. The observed richness of the antibiotic resistance genes repertoire was significantly reduced up to day 30, while a significant increase in the relative abundance of ß-lactamase encoding genes was observed up to day 10, consistent with a concomitant increase in the ß-lactamase activity of the microbiota. The level of ß-lactamase activity at baseline was positively associated with the resilience of the metabolome content of the stools. CONCLUSIONS: In healthy adults, antibiotics perturb many components of the microbiota, which return close to the baseline state within 30 days. These data suggest an important role of endogenous ß-lactamase-producing anaerobes in protecting the functions of the microbiota by de-activating the antibiotics reaching the colon. Video Abstract.

Reference-free Structural Variant Detection in Microbiomes via Long-read Coassembly Graphs.

Bacterial genome dynamics are vital for understanding the mechanisms underlying microbial adaptation, growth, and their broader impact on host phenotype. Structural variants (SVs), genomic alterations of 10 base pairs or more, play a pivotal role in driving evolutionary processes and maintaining genomic heterogeneity within bacterial populations. While SV detection in isolate genomes is relatively straightforward, metagenomes present broader challenges due to absence of clear reference genomes and presence of mixed strains. In response, our proposed method rhea, forgoes reference genomes and metagenome-assembled genomes (MAGs) by encompassing a single metagenome coassembly graph constructed from all samples in a series. The log fold change in graph coverage between subsequent samples is then calculated to call SVs that are thriving or declining throughout the series. We show rhea to outperform existing methods for SV and horizontal gene transfer (HGT) detection in two simulated mock metagenomes, which is particularly noticeable as the simulated reads diverge from reference genomes and an increase in strain diversity is incorporated. We additionally demonstrate use cases for rhea on series metagenomic data of environmental and fermented food microbiomes to detect specific sequence alterations between subsequent time and temperature samples, suggesting host advantage. Our innovative approach leverages raw read patterns rather than references or MAGs to include all sequencing reads in analysis, and thus provide versatility in studying SVs across diverse and poorly characterized microbial communities for more comprehensive insights into microbial genome dynamics.

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.

Cassette recombination dynamics within chromosomal integrons are regulated by toxin-antitoxin systems.

Integrons are adaptive bacterial devices that rearrange promoter-less gene cassettes into variable ordered arrays under stress conditions, thereby sampling combinatorial phenotypic diversity. Chromosomal integrons often carry hundreds of silent gene cassettes, with integrase-mediated recombination leading to rampant DNA excision and integration, posing a potential threat to genome integrity. How this activity is regulated and controlled, particularly through selective pressures, to maintain such large cassette arrays is unknown. Here, we show a key role of promoter-containing toxin-antitoxin (TA) cassettes as systems that kill the cell when the overall cassette excision rate is too high. These results highlight the importance of TA cassettes regulating the cassette recombination dynamics and provide insight into the evolution and success of integrons in bacterial genomes.

Integron cassettes integrate into bacterial genomes via widespread non-classical attG sites.

Integrons are genetic elements involved in bacterial adaptation which capture, shuffle and express genes encoding adaptive functions embedded in cassettes. These events are governed by the integron integrase through site-specific recombination between attC and attI integron sites. Using computational and molecular genetic approaches, here we demonstrate that the integrase also catalyses cassette integration into bacterial genomes outside of its known att sites. Once integrated, these cassettes can be expressed if located near bacterial promoters and can be excised at the integration point or outside, inducing chromosomal modifications in the latter case. Analysis of more than 5 × 105 independent integration events revealed a very large genomic integration landscape. We identified consensus recombination sequences, named attG sites, which differ greatly in sequence and structure from classical att sites. These results unveil an alternative route for dissemination of adaptive functions in bacteria and expand the role of integrons in bacterial evolution.

IMG/PR: a database of plasmids from genomes and metagenomes with rich annotations and metadata.

Plasmids are mobile genetic elements found in many clades of Archaea and Bacteria. They drive horizontal gene transfer, impacting ecological and evolutionary processes within microbial communities, and hold substantial importance in human health and biotechnology. To support plasmid research and provide scientists with data of an unprecedented diversity of plasmid sequences, we introduce the IMG/PR database, a new resource encompassing 699 973 plasmid sequences derived from genomes, metagenomes and metatranscriptomes. IMG/PR is the first database to provide data of plasmid that were systematically identified from diverse microbiome samples. IMG/PR plasmids are associated with rich metadata that includes geographical and ecosystem information, host taxonomy, similarity to other plasmids, functional annotation, presence of genes involved in conjugation and antibiotic resistance. The database offers diverse methods for exploring its extensive plasmid collection, enabling users to navigate plasmids through metadata-centric queries, plasmid comparisons and BLAST searches. The web interface for IMG/PR is accessible at https://img.jgi.doe.gov/pr. Plasmid metadata and sequences can be downloaded from https://genome.jgi.doe.gov/portal/IMG_PR.

Antimicrobial resistance level and conjugation permissiveness shape plasmid distribution in clinical enterobacteria.

Conjugative plasmids play a key role in the dissemination of antimicrobial resistance (AMR) genes across bacterial pathogens. AMR plasmids are widespread in clinical settings, but their distribution is not random, and certain associations between plasmids and bacterial clones are particularly successful. For example, the globally spread carbapenem resistance plasmid pOXA-48 can use a wide range of enterobacterial species as hosts, but it is usually associated with a small number of specific Klebsiella pneumoniae clones. These successful associations represent an important threat for hospitalized patients. However, knowledge remains limited about the factors determining AMR plasmid distribution in clinically relevant bacteria. Here, we combined in vitro and in vivo experimental approaches to analyze pOXA-48-associated AMR levels and conjugation dynamics in a collection of wild-type enterobacterial strains isolated from hospitalized patients. Our results revealed significant variability in these traits across different bacterial hosts, with Klebsiella spp. strains showing higher pOXA-48-mediated AMR and conjugation frequencies than Escherichia coli strains. Using experimentally determined parameters, we developed a simple mathematical model to interrogate the contribution of AMR levels and conjugation permissiveness to plasmid distribution in bacterial communities. The simulations revealed that a small subset of clones, combining high AMR levels and conjugation permissiveness, play a critical role in stabilizing the plasmid in different polyclonal microbial communities. These results help to explain the preferential association of plasmid pOXA-48 with K. pneumoniae clones in clinical settings. More generally, our study reveals that species- and strain-specific variability in plasmid-associated phenotypes shape AMR evolution in clinically relevant bacterial communities.

Ongoing shuffling of protein fragments diversifies core viral functions linked to interactions with bacterial hosts.

Biological modularity enhances evolutionary adaptability. This principle is vividly exemplified by bacterial viruses (phages), which display extensive genomic modularity. Phage genomes are composed of independent functional modules that evolve separately and recombine in various configurations. While genomic modularity in phages has been extensively studied, less attention has been paid to protein modularity-proteins consisting of distinct building blocks that can evolve and recombine, enhancing functional and genetic diversity. Here, we use a set of 133,574 representative phage proteins and highly sensitive homology detection to capture instances of domain mosaicism, defined as fragment sharing between two otherwise unrelated proteins, and to understand its relationship with functional diversity in phage genomes. We discover that unrelated proteins from diverse functional classes frequently share homologous domains. This phenomenon is particularly pronounced within receptor-binding proteins, endolysins, and DNA polymerases. We also identify multiple instances of recent diversification via domain shuffling in receptor-binding proteins, neck passage structures, endolysins and some members of the core replication machinery, often transcending distant taxonomic and ecological boundaries. Our findings suggest that ongoing diversification via domain shuffling is reflective of a co-evolutionary arms race, driven by the need to overcome various bacterial resistance mechanisms against phages.

Chance Favors the Prepared Genomes: Horizontal Transfer Shapes the Emergence of Antibiotic Resistance Mutations in Core Genes.

Bacterial lineages acquire novel traits at diverse rates in part because the genetic background impacts the successful acquisition of novel genes by horizontal transfer. Yet, how horizontal transfer affects the subsequent evolution of core genes remains poorly understood. Here, we studied the evolution of resistance to quinolones in Escherichia coli accounting for population structure. We found 60 groups of genes whose gain or loss induced an increase in the probability of subsequently becoming resistant to quinolones by point mutations in the gyrase and topoisomerase genes. These groups include functions known to be associated with direct mitigation of the effect of quinolones, with metal uptake, cell growth inhibition, biofilm formation, and sugar metabolism. Many of them are encoded in phages or plasmids. Although some of the chronologies may reflect epidemiological trends, many of these groups encoded functions providing latent phenotypes of antibiotic low-level resistance, tolerance, or persistence under quinolone treatment. The mutations providing resistance were frequent and accumulated very quickly. Their emergence was found to increase the rate of acquisition of other antibiotic resistances setting the path for multidrug resistance. Hence, our findings show that horizontal gene transfer shapes the subsequent emergence of adaptive mutations in core genes. In turn, these mutations further affect the subsequent evolution of resistance by horizontal gene transfer. Given the substantial gene flow within bacterial genomes, interactions between horizontal transfer and point mutations in core genes may be a key to the success of adaptation processes.

Emergence of novel non-aggregative variants under negative frequency-dependent selection in Klebsiella variicola.

Klebsiella variicola is an emergent human pathogen causing diverse infections, some of which in the urinary tract. However, little is known about the evolution and maintenance of genetic diversity in this species, the molecular mechanisms and their population dynamics. Here, we characterized the emergence of a novel rdar-like (rough and dry) morphotype which is contingent both on the genetic background and the environment. We show that mutations in either the nitrogen assimilation control gene (nac) or the type III fimbriae regulator, mrkH, suffice to generate rdar-like colonies. These morphotypes are primarily selected for the reduced inter-cellular aggregation as a result of MrkH loss-of-function which reduces type 3 fimbriae expression. Additionally, these clones also display increased growth rate and reduced biofilm formation. Direct competitions between rdar and wild type clones show that mutations in mrkH provide large fitness advantages. In artificial urine, the morphotype is under strong negative frequency-dependent selection and can socially exploit wild type strains. An exhaustive search for mrkH mutants in public databases revealed that ca 8% of natural isolates analysed had a truncated mrkH gene many of which were due to insertions of IS elements, including a reported clinical isolate with rdar morphology. These strains were rarely hypermucoid and often isolated from human, mostly from urine and blood. The decreased aggregation of these mutants could have important clinical implications as we hypothesize that such clones could better disperse within the host allowing colonisation of other body sites and potentially leading to systemic infections.

Latent evolution of biofilm formation depends on life-history and genetic background.

Adaptation to one environment can often generate phenotypic and genotypic changes which impact the future ability of an organism to thrive in other environmental conditions. In the context of host-microbe interactions, biofilm formation can increase survival rates in vivo upon exposure to stresses, like the host's immune system or antibiotic therapy. However, how the generic process of adaptation impacts the ability to form biofilm and how it may change through time has seldomly been studied. To do so, we used a previous evolution experiment with three strains of the Klebsiella pneumoniae species complex, in which we specifically did not select for biofilm formation. We observed that changes in the ability to form biofilm happened very fast at first and afterwards reverted to ancestral levels in many populations. Biofilm changes were associated to changes in population yield and surface polysaccharide production. Genotypically, mutations in the tip adhesin of type III fimbriae (mrkD) or the fim switch of type I fimbriae were shaped by nutrient availability during evolution, and their impact on biofilm formation was dependent on capsule production. Analyses of natural isolates revealed similar mutations in mrkD, suggesting that such mutations also play an important role in adaptation outside the laboratory. Our work reveals that the latent evolution of biofilm formation, and its temporal dynamics, depend on nutrient availability, the genetic background and other intertwined phenotypic and genotypic changes. Ultimately, it suggests that small differences in the environment can alter an organism's fate in more complex niches like the host.

Restriction-modification systems have shaped the evolution and distribution of plasmids across bacteria.

Many novel traits such as antibiotic resistance are spread by plasmids between species. Yet plasmids have different host ranges. Restriction-modification systems (R-M systems) are by far the most abundant bacterial defense system and therefore represent one of the key barriers to plasmid spread. However, their effect on plasmid evolution and host range has been neglected. Here we analyse the avoidance of targets of the most abundant R-M systems (Type II) for complete genomes and plasmids across bacterial diversity. For the most common target length (6 bp) we show that target avoidance is strongly correlated with the taxonomic distribution of R-M systems and is greater in plasmid genes than core genes. We find stronger avoidance of R-M targets in plasmids which are smaller and have a broader host range. Our results suggest two different evolutionary strategies for plasmids: small plasmids primarily adapt to R-M systems by tuning their sequence composition, and large plasmids primarily adapt through the carriage of additional genes protecting from restriction. Our work provides systematic evidence that R-M systems are important barriers to plasmid transfer and have left their mark on plasmids over long evolutionary time.

The Capsule Increases Susceptibility to Last-Resort Polymyxins, but Not to Other Antibiotics, in Klebsiella pneumoniae.

The extracellular capsule is a virulence factor present in many facultative pathogens, but its role in antimicrobial resistance remains controversial. To shed light on this debate, we tested six antibiotics on four Klebsiella pneumoniae species complex strains. Noncapsulated strains exhibited increased tolerance to polymyxins, but not to other antibiotics, as measured using the MIC. Our results urge caution on the use of therapeutic agents that target the capsule and may result in selection for its inactivation.