Rates Of Horizontal Gene Transfer Research Articles

Defence systems and horizontal gene transfer in bacteria.

Horizontal gene transfer (HGT) is a fundamental process in prokaryotic evolution, contributing significantly to diversification and adaptation. HGT is typically facilitated by mobile genetic elements (MGEs), such as conjugative plasmids and phages, which often impose fitness costs on their hosts. However, a considerable number of bacterial genes are involved in defence mechanisms that limit the propagation of MGEs, suggesting they may actively restrict HGT. In our study, we investigated whether defence systems limit HGT by examining the relationship between the HGT rate and the presence of 73 defence systems across 12 bacterial species. We discovered that only six defence systems, three of which were different CRISPR-Cas subtypes, were associated with a reduced gene gain rate at the species evolution scale. Hosts of these defence systems tend to have a smaller pangenome size and fewer phage-related genes compared to genomes without these systems. This suggests that these defence mechanisms inhibit HGT by limiting prophage integration. We hypothesize that the restriction of HGT by defence systems is species-specific and depends on various ecological and genetic factors, including the burden of MGEs and the fitness effect of HGT in bacterial populations.

Horizontal gene transfer facilitates the molecular reverse-evolution of antibiotic sensitivity in experimental populations of H. pylori.

Reversing the evolution of traits harmful to humans, such as antimicrobial resistance, is a key ambition of applied evolutionary biology. A major impediment to reverse evolution is the relatively low spontaneous mutation rates that revert evolved genotypes back to their ancestral state. However, the repeated re-introduction of ancestral alleles by horizontal gene transfer (HGT) could make reverse evolution likely. Here we evolve populations of an antibiotic-resistant strain of Helicobacter pylori in growth conditions without antibiotics while introducing an ancestral antibiotic-sensitive allele by HGT. We evaluate reverse evolution using DNA sequencing and find that HGT facilitates the molecular reverse evolution of the antibiotic resistance allele, and that selection for high rates of HGT drives the evolution of increased HGT rates in low-HGT treatment populations. Finally, we use a theoretical model and carry out simulations to infer how the fitness costs of antibiotic resistance, rates of HGT and effects of genetic drift interact to determine the probability and predictability of reverse evolution.

Metacommunity structure preserves genome diversity in the presence of gene-specific selective sweeps under moderate rates of horizontal gene transfer.

The horizontal transfer of genes is fundamental for the eco-evolutionary dynamics of microbial communities, such as oceanic plankton, soil, and the human microbiome. In the case of an acquired beneficial gene, classic population genetics would predict a genome-wide selective sweep, whereby the genome spreads clonally within the community and together with the beneficial gene, removing genome diversity. Instead, several sources of metagenomic data show the existence of "gene-specific sweeps", whereby a beneficial gene spreads across a bacterial community, maintaining genome diversity. Several hypotheses have been proposed to explain this process, including the decreasing gene flow between ecologically distant populations, frequency-dependent selection from linked deleterious allelles, and very high rates of horizontal gene transfer. Here, we propose an additional possible scenario grounded in eco-evolutionary principles. Specifically, we show by a mathematical model and simulations that a metacommunity where species can occupy multiple patches, acting together with a realistic (moderate) HGT rate, helps maintain genome diversity. Assuming a scenario of patches dominated by single species, our model predicts that diversity only decreases moderately upon the arrival of a new beneficial gene, and that losses in diversity can be quickly restored. We explore the generic behaviour of diversity as a function of three key parameters, frequency of insertion of new beneficial genes, migration rates and horizontal transfer rates.Our results provides a testable explanation for how diversity can be maintained by gene-specific sweeps even in the absence of high horizontal gene transfer rates.

Pangenome insights into the diversification and disease specificity of worldwide Xanthomonas outbreaks.

The bacterial genus Xanthomonas is responsible for disease outbreaks in several hundred plant species, many of them economically important crops. In the era of next-generation sequencing, thousands of strains from this genus have now been sequenced as part of isolated studies that focus on outbreak characterization, host range, diversity, and virulence factor identification. However, these data have not been synthesized and we lack a comprehensive phylogeny for the genus, with some species designations in public databases still relying on phenotypic similarities and representative sequence typing. The extent of genetic cohesiveness among Xanthomonas strains, the distribution of virulence factors across strains, and the impact of evolutionary history on host range across the genus are also poorly understood. In this study, we present a pangenome analysis of 1,910 diverse Xanthomonas genomes, highlighting their evolutionary relationships, the distribution of virulence-associated genes across strains, and rates of horizontal gene transfer. We find a number of broadly conserved classes of virulence factors and considerable diversity in the Type 3 Secretion Systems (T3SSs) and Type 3 Secreted Effector (T3SE) repertoires of different Xanthomonas species. We also use these data to re-assign incorrectly classified strains to phylogenetically informed species designations and find evidence of both monophyletic host specificity and convergent evolution of phylogenetically distant strains to the same host. Finally, we explore the role of recombination in maintaining genetic cohesion within the Xanthomonas genus as a result of both ancestral and recent recombination events. Understanding the evolutionary history of Xanthomonas species and the relationship of key virulence factors with host-specificity provides valuable insight into the mechanisms through which Xanthomonas species shift between hosts and will enable us to develop more robust resistance strategies against these highly virulent pathogens.

Biofilm Lifestyle in Recurrent Urinary Tract Infections.

Urinary tract infections (UTIs) represent one of the most common infections that are frequently encountered in health care facilities. One of the main mechanisms used by bacteria that allows them to survive hostile environments is biofilm formation. Biofilms are closed bacterial communities that offer protection and safe hiding, allowing bacteria to evade host defenses and hide from the reach of antibiotics. Inside biofilm communities, bacteria show an increased rate of horizontal gene transfer and exchange of resistance and virulence genes. Additionally, bacterial communication within the biofilm allows them to orchestrate the expression of virulence genes, which further cements the infestation and increases the invasiveness of the infection. These facts stress the necessity of continuously updating our information and understanding of the etiology, pathogenesis, and eradication methods of this growing public health concern. This review seeks to understand the role of biofilm formation in recurrent urinary tact infections by outlining the mechanisms underlying biofilm formation in different uropathogens, in addition to shedding light on some biofilm eradication strategies.

Transposable elements drive intron gain in diverse eukaryotes

There is massive variation in intron numbers across eukaryotic genomes, yet the major drivers of intron content during evolution remain elusive. Rapid intron loss and gain in some lineages contrast with long-term evolutionary stasis in others. Episodic intron gain could be explained by recently discovered specialized transposons called Introners, but so far Introners are only known from a handful of species. Here, we performed a systematic search across 3,325 eukaryotic genomes and identified 27,563 Introner-derived introns in 175 genomes (5.2%). Species with Introners span remarkable phylogenetic diversity, from animals to basal protists, representing lineages whose last common ancestor dates to over 1.7 billion years ago. Aquatic organisms were 6.5 times more likely to contain Introners than terrestrial organisms. Introners exhibit mechanistic diversity but most are consistent with DNA transposition, indicating that Introners have evolved convergently hundreds of times from nonautonomous transposable elements. Transposable elements and aquatic taxa are associated with high rates of horizontal gene transfer, suggesting that this combination of factors may explain the punctuated and biased diversity of species containing Introners. More generally, our data suggest that Introners may explain the episodic nature of intron gain across the eukaryotic tree of life. These results illuminate the major source of ongoing intron creation in eukaryotic genomes.

Horizontal gene transfer enables programmable gene stability in synthetic microbiota.

The functions of many microbial communities exhibit remarkable stability despite fluctuations in the compositions of these communities. To date, a mechanistic understanding of this function-composition decoupling is lacking. Statistical mechanisms have been commonly hypothesized to explain such decoupling. Here, we proposed that dynamic mechanisms, mediated by horizontal gene transfer (HGT), also enable the independence of functions from the compositions of microbial communities. We combined theoretical analysis with numerical simulations to illustrate that HGT rates can determine the stability of gene abundance in microbial communities. We further validated these predictions using engineered microbial consortia of different complexities transferring one or more than a dozen clinically isolated plasmids, as well as through the reanalysis of data from the literature. Our results demonstrate a generalizable strategy to program the gene stability of microbial communities.

Organic fertilizer potentiates the transfer of typical antibiotic resistance gene among special bacterial species

The propagation of antibiotic resistance genes (ARGs) in environments has evoked many attentions, however, how to identify their host pathogenic bacteria in situ remains a great challenge. Here we explored the bacterial host distribution and dissemination of a typical ARG, sul1 gene, in agricultural soils through the simultaneous detection of sul1 and its host 16S rRNA gene by emulsion paired isolation and concatenation PCR (epicPCR). Compared to chemical fertilizer, organic fertilizer (chicken manure) led to a higher prevalence of sul1 gene in the soil, and dominant bacterial hosts of sul1 gene were classified into Proteobacteria and Bacteroidetes phyla. Additionally, significant higher diversity of antibiotic resistance bacteria (ARB), higher rate of horizontal gene transfer (HGT), higher rate of mobile genetic elements (MGE) and higher proportion of pathogens were all observed in the treatment of organic fertilizer. This study alerts potential health risks of manure applications in agricultural soils.

Bacterial chromosomal mobility via lateral transduction exceeds that of classical mobile genetic elements

It is commonly assumed that the horizontal transfer of most bacterial chromosomal genes is limited, in contrast to the frequent transfer observed for typical mobile genetic elements. However, this view has been recently challenged by the discovery of lateral transduction in Staphylococcus aureus, where temperate phages can drive the transfer of large chromosomal regions at extremely high frequencies. Here, we analyse previously published as well as new datasets to compare horizontal gene transfer rates mediated by different mechanisms in S. aureus and Salmonella enterica. We find that the horizontal transfer of core chromosomal genes via lateral transduction can be more efficient than the transfer of classical mobile genetic elements via conjugation or generalized transduction. These results raise questions about our definition of mobile genetic elements, and the potential roles played by lateral transduction in bacterial evolution.

Genome and Methylome analysis of a phylogenetic novel Campylobacter coli cluster with C. jejuni introgression.

The intriguing recent discovery of Campylobacter coli strains, especially of clade 1, that (i) possess mosaic C. coli / C. jejuni alleles, (ii) demonstrate mixed multilocus sequence types (MLSTs) and (iii) have undergone genome-wide introgression has led to the speculation that these two species may be involved in an accelerated rate of horizontal gene transfer that is progressively leading to the merging of both species in a process coined ‘despeciation’. In an MLST-based neighbour-joining tree of a number of C. coli and C. jejuni isolates of different clades, three prominent Campylobacter isolates formed a seemingly separate cluster besides the previously described C. coli and C. jejuni clades. In the light of the suspected, ongoing genetic introgression between the C. coli and C. jejuni species, this cluster of Campylobacter isolates is proposed to present one of the hybrid clonal complexes in the despeciation process of the genus. Specific DNA methylation as well as restriction modification systems are known to be involved in selective uptake of external DNA and their role in such genetic introgression remains to be further investigated. In this study, the phylogeny and DNA methylation of these putative C. coli / C. jejuni hybrid strains were explored, their genomic mosaic structure caused by C. jejuni introgression was demonstrated and basic phenotypic assays were used to characterize these isolates. The genomes of the three hybrid Campylobacter strains were sequenced using PacBio SMRT sequencing, followed by methylome analysis by Restriction-Modification Finder and genome analysis by Parsnp, Smash++ and blast. Additionally, the strains were phenotypically characterized with respect to growth behaviour, motility, eukaryotic cell invasion and adhesion, autoagglutination, biofilm formation, and water survival ability. Our analyses show that the three hybrid Campylobacter strains are clade 1 C . coli strains, which have acquired between 8.1 and 9.1 % of their genome from C. jejuni . The C. jejuni genomic segments acquired are distributed over the entire genome and do not form a coherent cluster. Most of the genes originating from C. jejuni are involved in chemotaxis and motility, membrane transport, cell signalling, or the resistance to toxic compounds such as bile acids. Interspecies gene transfer from C. jejuni has contributed 8.1–9.1% to the genome of three C. coli isolates and initiated the despeciation between C. jejuni and C. coli . Based on their functional annotation, the genes originating from C. jejuni enable the adaptation of the three strains to an intra-intestinal habitat. The transfer of a fused type II restriction-modification system that recognizes the CAYNNNNNCTC/GAGNNNNNRTG motif seems to be the key for the recombination of the C. jejuni genetic material with C. coli genomes.

Diverse Viruses in Deep-Sea Hydrothermal Vent Fluids Have Restricted Dispersal across Ocean Basins.

ABSTRACTIn the ocean, viruses impact microbial mortality, regulate biogeochemical cycling, and alter the metabolic potential of microbial lineages. At deep-sea hydrothermal vents, abundant viruses infect a wide range of hosts among the archaea and bacteria that inhabit these dynamic habitats. However, little is known about viral diversity, host range, and biogeography across different vent ecosystems, which has important implications for how viruses manipulate microbial function and evolution. Here, we examined viral diversity, viral and host distribution, and virus-host interactions in microbial metagenomes generated from venting fluids from several vent sites within three different geochemically and geographically distinct hydrothermal systems: Piccard and Von Damm vent fields at the Mid-Cayman Rise in the Caribbean Sea, and at several vent sites within Axial Seamount in the Pacific Ocean. Analysis of viral sequences and clustered regularly interspaced short palindromic repeat (CRISPR) spacers revealed highly diverse viral assemblages and evidence of active infection. Network analysis revealed that viral host range was relatively narrow, with very few viruses infecting multiple microbial lineages. Viruses were largely endemic to individual vent sites, indicating restricted dispersal, and in some cases, viral assemblages persisted over time. Thus, we show that hydrothermal vent fluids are home to novel, diverse viral assemblages that are highly localized to specific regions and taxa.IMPORTANCE Viruses play important roles in manipulating microbial communities and their evolution in the ocean, yet not much is known about viruses in deep-sea hydrothermal vents. However, viral ecology and evolution are of particular interest in hydrothermal vent habitats because of their unique nature: previous studies have indicated that most viruses in hydrothermal vents are temperate rather than lytic, and it has been established that rates of horizontal gene transfer (HGT) are particularly high among thermophilic vent microbes, and viruses are common vectors for HGT. If viruses have broad host range or are widespread across vent sites, they have increased potential to act as gene-sharing “highways” between vent sites. By examining viral diversity, distribution, and infection networks across disparate vent sites, this study provides the opportunity to better characterize and constrain the viral impact on hydrothermal vent microbial communities. We show that viruses in hydrothermal vents are diverse and apparently active, but most have restricted host range and are not widely distributed among vent sites. Thus, the impacts of viral infection are likely to be highly localized and constrained to specific taxa in these habitats.

Expanding Diversity of Firmicutes Single-Strand Annealing Proteins: A Putative Role of Bacteriophage-Host Arms Race.

Bacteriophage-encoded single strand annealing proteins (SSAPs) are recombinases which can substitute the classical, bacterial RecA and manage the DNA metabolism at different steps of phage propagation. SSAPs have been shown to efficiently promote recombination between short and rather divergent DNA sequences and were exploited for in vivo genetic engineering mainly in Gram-negative bacteria. In opposition to the conserved and almost universal bacterial RecA protein, SSAPs display great sequence diversity. The importance for SSAPs in phage biology and phage-bacteria evolution is underlined by their role as key players in events of horizontal gene transfer (HGT). All of the above provoke a constant interest for the identification and study of new phage recombinase proteins in vivo, in vitro as well as in silico. Despite this, a huge body of putative ssap genes escapes conventional classification, as they are not properly annotated. In this work, we performed a wide-scale identification, classification and analysis of SSAPs encoded by the Firmicutes bacteria and their phages. By using sequence similarity network and gene context analyses, we created a new high quality dataset of phage-related SSAPs, substantially increasing the number of annotated SSAPs. We classified the identified SSAPs into seven distinct families, namely RecA, Gp2.5, RecT/Redβ, Erf, Rad52/22, Sak3, and Sak4, organized into three superfamilies. Analysis of the relationships between the revealed protein clusters led us to recognize Sak3-like proteins as a new distinct SSAP family. Our analysis showed an irregular phylogenetic distribution of ssap genes among different bacterial phyla and specific phages, which can be explained by the high rates of ssap HGT. We propose that the evolution of phage recombinases could be tightly linked to the dissemination of bacterial phage-resistance mechanisms (e.g., abortive infection and CRISPR/Cas systems) targeting ssap genes and be a part of the constant phage-bacteria arms race.

Assumption and Criticism on RNA World Hypothesis from Ribozymes to Functional Cells

The origin of life has long baffled scientists, trying to shed some light on the universe of the unknown we live in. The RNA world hypothesis is merely an idea that tries to grasp the root of life by going beyond phylogenetic trees, and the ‘Drawinian threshold'. The birth of this hypothesis came with the discovery of ribozymes and certain relic molecules left behind by an RNA ruled world about -3.8 billion years ago. These primitive cells were the most minimal forms of life and experiments are trying to show this by mimicking primitive earth environments. Nucleotides have been shown to polymerize into random chains easily on clay or in ice and the chemical synthesis of RNA is possible although extremely challenging without the modern cellular machinery. The RNA world is untraceable because of horizontal gene transfer, which is highly criticized by researches, and still remains a widely researched and studied topic as it could be the heart of life on other earth like planets in the universe. The state of the model needs improvements in analytical techniques and experiments for the rate of Horizontal Gene Transfer on loosely bound components, and continued comparison of the different types of phylogenetic trees to widen the evolution of the genetic core. Furthermore the development of analytical techniques used for the monitoring of the Donan effect may show more insight into the division and growth of primitive cellular membranes and the continued research and exploration in the fields of Ribozymes, with a goal to reach a higher catalytic range with high fidelity and rate.

Elevated rates of horizontal gene transfer in the industrialized human microbiome

Industrialization has impacted the human gut ecosystem, resulting in altered microbiome composition and diversity. Whether bacterial genomes may also adapt to the industrialization of their host populations remains largely unexplored. Here, we investigate the extent to which the rates and targets of horizontal gene transfer (HGT) vary across thousands of bacterial strains from 15 human populations spanning a range of industrialization. We show that HGTs have accumulated in the microbiome over recent host generations and that HGT occurs at high frequency within individuals. Comparison across human populations reveals that industrialized lifestyles are associated with higher HGT rates and that the functions of HGTs are related to the level of host industrialization. Our results suggest that gut bacteria continuously acquire new functionality based on host lifestyle and that high rates of HGT may be a recent development in human history linked to industrialization.

Investigating the effects of municipal and hospital wastewaters on horizontal gene transfer.

Horizontal gene transfer (HGT) plays an important role in the dissemination of antibiotic resistance genes. In sewer systems, human-associated and environmental bacteria are mixed together and exposed to many substances known to increase HGT, including various antibacterial compounds. In wastewaters, those substances are most often detected below concentrations known to induce HGT individually. Still, it is possible that such wastewaters induce HGT, for example via mixture effects. Here, a panel of antibiotics, biocides and other pharmaceuticals was measured in filter-sterilized municipal and hospital wastewater samples from Gothenburg, Sweden. The effects on HGT of the chemical mixtures in these samples were investigated by exposing a complex bacterial donor community together with a GFP-tagged E.coli recipient strain. Recipients that captured sulfonamide resistance-conferring mobile genetic elements (MGEs) from the bacterial community were enumerated and characterized by replicon typing, antibiotic susceptibility testing and long read sequencing. While exposure to municipal wastewater did not result in any detectable change in HGT rates, exposure to hospital wastewater was associated with an increase in the proportion of recipients that acquired sulfonamide resistance but also a drastic decrease in the total number of recipients. Although, concentrations were generally higher in hospital than municipal wastewater, none of the measured substances could individually explain the observed effects of hospital wastewater. The great majority of the MGEs captured were IncN plasmids, and resistance to several antibiotics was co-transferred in most cases. Taken together, the data show no evidence that chemicals present in the studied municipal wastewater induce HGT. Still, the increased relative abundance of transconjugants after exposure to hospital wastewater could have implications for the risks of both emergence and transmission of resistant bacteria.

Deciphering functional redundancy in the human microbiome

Although the taxonomic composition of the human microbiome varies tremendously across individuals, its gene composition or functional capacity is highly conserved — implying an ecological property known as functional redundancy. Such functional redundancy has been hypothesized to underlie the stability and resilience of the human microbiome, but this hypothesis has never been quantitatively tested. The origin of functional redundancy is still elusive. Here, we investigate the basis for functional redundancy in the human microbiome by analyzing its genomic content network — a bipartite graph that links microbes to the genes in their genomes. We find that this network exhibits several topological features that favor high functional redundancy. Furthermore, we develop a simple genome evolution model to generate genomic content network, finding that moderate selection pressure and high horizontal gene transfer rate are necessary to generate genomic content networks with key topological features that favor high functional redundancy. Finally, we analyze data from two published studies of fecal microbiota transplantation (FMT), finding that high functional redundancy of the recipient’s pre-FMT microbiota raises barriers to donor microbiota engraftment. This work elucidates the potential ecological and evolutionary processes that create and maintain functional redundancy in the human microbiome and contribute to its resilience.

A New Phylogenomic Approach For Quantifying Horizontal Gene Transfer Trends in Prokaryotes

It is well established nowadays that among prokaryotes, various families of orthologous genes exhibit conflicting evolutionary history. A prime factor for this conflict is horizontal gene transfer (HGT) - the transfer of genetic material not via vertical descent. Thus, the prevalence of HGT is challenging the meaningfulness of the classical Tree of Life concept. Here we present a comprehensive study of HGT representing the entire prokaryotic world. We mainly rely on a novel analytic approach for analyzing an aggregate of gene histories, by means of the quartet plurality distribution (QPD) that we develop. Through the analysis of real and simulated data, QPD is used to reveal evidence of a barrier against HGT, separating the archaea from the bacteria and making HGT between the two domains, in general, quite rare. In contrast, bacteria’s confined HGT is substantially more frequent than archaea’s. Our approach also reveals that despite intensive HGT, a strong tree-like signal can be extracted, corroborating several previous works. Thus, QPD, which enables one to analytically combine information from an aggregate of gene trees, can be used for understanding patterns and rates of HGT in prokaryotes, as well as for validating or refuting models of horizontal genetic transfers and evolution in general.

Population genomics of Vibrionaceae isolated from an endangered oasis reveals local adaptation after an environmental perturbation

BackgroundIn bacteria, pan-genomes are the result of an evolutionary “tug of war” between selection and horizontal gene transfer (HGT). High rates of HGT increase the genetic pool and the effective population size (Ne), resulting in open pan-genomes. In contrast, selective pressures can lead to local adaptation by purging the variation introduced by HGT and mutation, resulting in closed pan-genomes and clonal lineages. In this study, we explored both hypotheses, elucidating the pan-genome of Vibrionaceae isolates after a perturbation event in the endangered oasis of Cuatro Ciénegas Basin (CCB), Mexico, and looking for signals of adaptation to the environments in their genomes.ResultsWe obtained 42 genomes of Vibrionaceae distributed in six lineages, two of them did not showed any close reference strain in databases. Five of the lineages showed closed pan-genomes and were associated to either water or sediment environment; their high Ne estimates suggest that these lineages are not from a recent origin. The only clade with an open pan-genome was found in both environments and was formed by ten genetic groups with low Ne, suggesting a recent origin. The recombination and mutation estimators (r/m) ranged from 0.005 to 2.725, which are similar to oceanic Vibrionaceae estimations. However, we identified 367 gene families with signals of positive selection, most of them found in the core genome; suggesting that despite recombination, natural selection moves the Vibrionaceae CCB lineages to local adaptation, purging the genomes and keeping closed pan-genome patterns. Moreover, we identify 598 SNPs associated with an unstructured environment; some of the genes associated with these SNPs were related to sodium transport.ConclusionsDifferent lines of evidence suggest that the sampled Vibrionaceae, are part of the rare biosphere usually living under famine conditions. Two of these lineages were reported for the first time. Most Vibrionaceae lineages of CCB are adapted to their micro-habitats rather than to the sampled environments. This pattern of adaptation is concordant with the association of closed pan-genomes and local adaptation.

Horizontal gene transfer rate is not the primary determinant of observed antibiotic resistance frequencies in Streptococcus pneumoniae.

The extent to which evolution is constrained by the rate at which horizontal gene transfer (HGT) allows DNA to move between genetic lineages is an open question, which we address in the context of antibiotic resistance in Streptococcus pneumoniae. We analyze microbiological, genomic, and epidemiological data from the largest-to-date sequenced pneumococcal carriage study in 955 infants from a refugee camp on the Thailand-Myanmar border. Using a unified framework, we simultaneously test prior hypotheses on rates of HGT and a key evolutionary covariate (duration of carriage) as determinants of resistance frequencies. We conclude that in this setting, there is little evidence of HGT playing a major role in determining resistance frequencies. Instead, observed resistance frequencies are best explained as the outcome of selection acting on a pool of variants, irrespective of the rate at which resistance determinants move between genetic lineages.

Diversity and evolution of surface polysaccharide synthesis loci in Enterobacteriales

Bacterial capsules and lipopolysaccharides are diverse surface polysaccharides (SPs) that serve as the frontline for interactions with the outside world. While SPs can evolve rapidly, their diversity and evolutionary dynamics across different taxonomic scales has not been investigated in detail. Here, we focused on the bacterial order Enterobacteriales (including the medically relevant Enterobacteriaceae), to carry out comparative genomics of two SP locus synthesis regions, cps and kps, using 27,334 genomes from 45 genera. We identified high-quality cps loci in 22 genera and kps in 11 genera, around 4% of which were detected in multiple species. We found SP loci to be highly dynamic genetic entities: their evolution was driven by high rates of horizontal gene transfer (HGT), both of whole loci and component genes, and relaxed purifying selection, yielding large repertoires of SP diversity. In spite of that, we found the presence of (near-)identical locus structures in distant taxonomic backgrounds that could not be explained by recent exchange, pointing to long-term selective preservation of locus structures in some populations. Our results reveal differences in evolutionary dynamics driving SP diversity within different bacterial species, with lineages of Escherichia coli, Enterobacter hormaechei and Klebsiella aerogenes most likely to share SP loci via recent exchange; and lineages of Salmonella enterica, Cronobacter sakazakii and Serratia marcescens most likely to share SP loci via other mechanisms such as long-term preservation. Overall, the evolution of SP loci in Enterobacteriales is driven by a range of evolutionary forces and their dynamics and relative importance varies between different species.