The chromosomal accommodation and domestication of mobile genetic elements

Dinotefuran exposure alters biochemical, metabolomic, gut microbiome, and growth responses in decapoda pacific white shrimp Penaeus vannamei

Dinotefuran, a neonicotinoid insecticide, may impact nontarget organisms such as Decapoda P. vannamei shrimp with nervous systems similar to insects. Exposing shrimp to low dinotefuran concentrations (6, 60, and 600 μg/L) for 21 days affected growth, hepatosomatic index, and survival. Biomarkers erythromycin-N-demethylase, alanine aminotransferase, and catalase increased in all exposed groups, while glutathione S-transferase is the opposite; aminopyrin-N-demethylase, malondialdehyde, and aspartate aminotransferase increased at 60 and 600 μg/L. Concentration-dependent effects on gut microbiota altered the abundance of bacterial groups, increased potentially pathogenic and oxidative stress-resistant phenotypes, and decreased biofilm formation. Gram-positive/negative microbiota changed significantly. Metabolite differences between the exposed and control groups were identified using mass spectrometry and KEGG pathway enrichment. N-acetylcystathionine showed potential as a reliable dinotefuran metabolic marker. Weighted correlation network analysis (WGCNA) results indicated high connectivity of cruecdysone in the metabolite network and significant enrichment at 600 μg/L dinotefuran. The WGCNA results revealed a highly significant negative correlation between two key metabolites, caldine and indican, and the gut microbiota within co-expression modules. Overall, the risk of dinotefuran exposure to non-target organisms in aquatic environments still requires further attention.

Prophage maintenance is determined by environment-dependent selective sweeps rather than mutational availability

Prophages, viral sequences integrated into bacterial genomes, can be beneficial and costly. Despite the risk of prophage activation and subsequent bacterial death, active prophages are present in most bacterial genomes. However, our understanding of the selective forces that maintain prophages in bacterial populations is limited. Combining experimental evolution with stochastic modeling, we show that prophage maintenance and loss are primarily determined by environmental conditions that alter the net fitness effect of a prophage on its bacterial host. When prophages are too costly, they are rapidly lost through environment-specific sequences of selective sweeps. Conflicting selection pressures that select against the prophage but for a prophage-encoded accessory gene can maintain prophages. The dynamics of prophage maintenance additionally depend on the sociality of this accessory gene. Prophage-encoded genes that exclusively benefit the lysogen maintain prophages at higher frequencies compared with genes that benefit the entire population. That is because the latter can protect phage-free "cheaters," reducing the benefit of maintaining the prophage. Our simulations suggest that environmental variation plays a larger role than mutation rates in determining prophage maintenance. These findings highlight the complexity of selection pressures that act on mobile genetic elements and challenge our understanding of the role of environmental factors relative to random chance events in shaping the evolutionary trajectory of bacterial populations. By shedding light on the key factors that shape microbial populations in the face of environmental changes, our study significantly advances our understanding of the complex dynamics of microbial evolution and diversification.

Structure and dispersion of the conjugative mobilome in surface ocean bacterioplankton

Mobile genetic elements (MGEs), collectively referred to as the "mobilome", can have a significant impact on the fitness of microbial communities and therefore on ecological processes. Marine MGEs have mainly been associated with wide geographical and phylogenetic dispersal of adaptative traits. However, whether the structure of this mobilome exhibits deterministic patterns in the natural community is still an open question. The aim of this study was to characterize the structure of the conjugative mobilome in the ocean surface bacterioplankton by searching the publicly available marine metagenomes from the TARA Oceans survey, together with molecular markers, such as relaxases and type IV coupling proteins of the type IV secretion system (T4SS). The T4SS machinery was retrieved in more abundance than relaxases in the surface marine bacterioplankton. Moreover, among the identified MGEs, mobilizable elements were the most abundant, outnumbering self-conjugative sequences. Detection of a high number of incomplete T4SSs provides insight into possible strategies related to trans-acting activity between MGEs, and accessory functions of the T4SS (e.g. protein secretion), allowing the host to maintain a lower metabolic burden in the highly dynamic marine system. Additionally, the results demonstrate a wide geographical dispersion of MGEs throughout oceanic regions, while the Southern Ocean appears segregated from other regions. The marine mobilome also showed a high similarity of functions present in known plasmid databases. Moreover, cargo genes were mostly related to DNA processing, but scarcely associated with antibiotic resistance. Finally, within the MGEs, integrative and conjugative elements showed wider marine geographic dispersion than plasmids.

Evolutionary and ecological role of extracellular contractile injection systems: from threat to weapon

Contractile injection systems (CISs) are phage tail-related structures that are encoded in many bacterial genomes. These devices encompass the cell-based type VI secretion systems (T6SSs) as well as extracellular CISs (eCISs). The eCISs comprise the R-tailocins produced by various bacterial species as well as related phage tail-like structures such as the antifeeding prophages (Afps) of Serratia entomophila, the Photorhabdus virulence cassettes (PVCs), and the metamorphosis-associated contractile structures (MACs) of Pseudoalteromonas luteoviolacea. These contractile structures are released into the extracellular environment upon suicidal lysis of the producer cell and play important roles in bacterial ecology and evolution. In this review, we specifically portray the eCISs with a focus on the R-tailocins, sketch the history of their discovery and provide insights into their evolution within the bacterial host, their structures and how they are assembled and released. We then highlight ecological and evolutionary roles of eCISs and conceptualize how they can influence and shape bacterial communities. Finally, we point to their potential for biotechnological applications in medicine and agriculture.

One Earth: The Equilibrium between the Human and the Bacterial Worlds

Misuse and abuse of antibiotics on humans, cattle, and crops have led to the selection of multi-resistant pathogenic bacteria, the most feared 'superbugs'. Infections caused by superbugs are progressively difficult to treat, with a subsequent increase in lethality: the toll on human lives is predicted to reach 10 million by 2050. Here we review three concepts linked to the growing resistance to antibiotics, namely (i) the Resistome, which refers to the collection of bacterial genes that confer resistance to antibiotics, (ii) the Mobilome, which includes all the mobile genetic elements that participate in the spreading of antibiotic resistance among bacteria by horizontal gene transfer processes, and (iii) the Nichome, which refers to the set of genes that are expressed when bacteria try to colonize new niches. We also discuss the strategies that can be used to tackle bacterial infections and propose an entente cordiale with the bacterial world so that instead of war and destruction of the 'fierce enemy' we can achieve a peaceful coexistence (the One Earth concept) between the human and the bacterial worlds. This, in turn, will contribute to microbial biodiversity, which is crucial in a globally changing climate due to anthropogenic activities.

Modularity and diversity of target selectors in Tn7 transposons

To spread, transposons must integrate into target sites without disruption of essential genes while avoiding host defense systems. Tn7-like transposons employ multiple mechanisms for target-site selection, including protein-guided targeting and, in CRISPR-associated transposons (CASTs), RNA-guided targeting. Combining phylogenomic and structural analyses, we conducted a broad survey of target selectors, revealing diverse mechanisms used by Tn7 to recognize target sites, including previously uncharacterized target-selector proteins found in newly discovered transposable elements (TEs). We experimentally characterized a CAST I-D system and a Tn6022-like transposon that uses TnsF, which contains an inactivated tyrosine recombinase domain, to target the comM gene. Additionally, we identified a non-Tn7 transposon, Tsy, encoding a homolog of TnsF with an active tyrosine recombinase domain, which we show also inserts into comM. Our findings show that Tn7 transposons employ modular architecture and co-opt target selectors from various sources to optimize target selection and drive transposon spread.

Systematic analysis of prophage elements in actinobacterial genomes reveals a remarkable phylogenetic diversity

Actinobacteria represent one of the largest bacterial phyla harboring many species of high medical, biotechnological and ecological relevance. Prophage elements are major contributors to bacterial genome diversity and were shown to significantly shape bacterial fitness and host-microbe interactions. In this study, we performed a systematic analysis of prophage elements in 2406 complete actinobacterial genomes. Overall, 2106 prophage elements were predicted to be present in about 50% (1172/2406) of the analyzed datasets. Interestingly, these identified sequences compose a high prevalence of cryptic prophage elements, indicating genetic decay and domestication. Analysis of the sequence relationship of predicted prophages with known actinobacteriophage genomes revealed an exceptional high phylogenetic diversity of prophage elements. As a trend, we observed a higher prevalence of prophage elements in vicinity to the terminus. Analysis of the prophage-encoded gene functions revealed that prophage sequences significantly contribute to the bacterial antiviral immune system, but no biosynthetic gene clusters involved in the synthesis of known antiphage molecules were identified in prophage genomes. Overall, the current study highlights the remarkable diversity of prophages in actinobacterial genomes, with highly divergent prophages in actinobacterial genomes and thus provides an important basis for further investigation of phage-host interactions in this important bacterial phylum.

Sending out an SOS - the bacterial DNA damage response

The term “SOS response” was first coined by Radman in 1974, in an intellectual effort to put together the data suggestive of a concerted gene expression program in cells undergoing DNA damage. A large amount of information about this cellular response has been collected over the following decades. In this review, we will focus on a few of the relevant aspects about the SOS response: its mechanism of control and the stressors which activate it, the diversity of regulated genes in different species, its role in mutagenesis and evolution including the development of antimicrobial resistance, and its relationship with mobile genetic elements.

Analysis of intact prophages in genomes of Paenibacillus larvae: An important pathogen for bees

Paenibacillus larvae is the etiological agent of American Foulbrood (AFB), a highly contagious and worldwide spread bacterial disease that affects honeybee brood. In this study, all complete P. larvae genomes available on the NCBI database were analyzed in order to detect presence of prophages using the PHASTER software. A total of 55 intact prophages were identified in 11 P. larvae genomes (5.0 ± 2.3 per genome) and were further investigated for the presence of genes encoding relevant traits related to P. larvae. A closer look at the prophage genomes revealed the presence of several putative genes such as metabolic and antimicrobial resistance genes, toxins or bacteriocins, potentially influencing host performance. Some of the coding DNA sequences (CDS) were present in all ERIC-genotypes, while others were only found in a specific genotype. While CDS encoding toxins and antitoxins such as HicB and MazE were found in prophages of all bacterial genotypes, others, from the same category, were provided by prophages particularly to ERIC I (enhancin-like toxin), ERIC II (antitoxin SocA) and ERIC V strains (subunit of Panton-Valentine leukocidin system (PVL) LukF-PV). This is the first in-depth analysis of P. larvae prophages. It provides better knowledge on their impact in the evolution of virulence and fitness of P. larvae, by discovering new features assigned by the viruses.

Phage-host coevolution in natural populations

Coevolution between bacteriophages (phages) and their bacterial hosts occurs through changes in resistance and counter-resistance mechanisms. To assess phage-host evolution in wild populations, we isolated 195 Vibrio crassostreae strains and 243 vibriophages during a 5-month time series from an oyster farm and combined these isolates with existing V. crassostreae and phage isolates. Cross-infection studies of 81,926 host-phage pairs delineated a modular network where phages are best at infecting co-occurring hosts, indicating local adaptation. Successful propagation of phage is restricted by the ability to adsorb to closely related bacteria and further constrained by strain-specific defence systems. These defences are highly diverse and predominantly located on mobile genetic elements, and multiple defences are active within a single genome. We further show that epigenetic and genomic modifications enable phage to adapt to bacterial defences and alter host range. Our findings reveal that the evolution of bacterial defences and phage counter-defences is underpinned by frequent genetic exchanges with, and between, mobile genetic elements.

Phages and their satellites encode hotspots of antiviral systems

Bacteria carry diverse genetic systems to defend against viral infection, some of which are found within prophages where they inhibit competing viruses. Phage satellites pose additional pressures on phages by hijacking key viral elements to their own benefit. Here, we show that E. coli P2-like phages and their parasitic P4-like satellites carry hotspots of genetic variation containing reservoirs of anti-phage systems. We validate the activity of diverse systems and describe PARIS, an abortive infection system triggered by a phage-encoded anti-restriction protein. Antiviral hotspots participate in inter-viral competition and shape dynamics between the bacterial host, P2-like phages, and P4-like satellites. Notably, the anti-phage activity of satellites can benefit the helper phage during competition with virulent phages, turning a parasitic relationship into a mutualistic one. Anti-phage hotspots are present across distant species and constitute a substantial source of systems that participate in the competition between mobile genetic elements.

Microbial defenses against mobile genetic elements and viruses: Who defends whom from what?

Prokaryotes have numerous mobile genetic elements (MGEs) that mediate horizontal gene transfer (HGT) between cells. These elements can be costly, even deadly, and cells use numerous defense systems to filter, control, or inactivate them. Recent studies have shown that prophages, conjugative elements, their parasites (phage satellites and mobilizable elements), and other poorly described MGEs encode defense systems homologous to those of bacteria. These constitute a significant fraction of the repertoire of cellular defense genes. As components of MGEs, these defense systems have presumably evolved to provide them, not the cell, adaptive functions. While the interests of the host and MGEs are aligned when they face a common threat such as an infection by a virulent phage, defensive functions carried by MGEs might also play more selfish roles to fend off other antagonistic MGEs or to ensure their maintenance in the cell. MGEs are eventually lost from the surviving host genomes by mutational processes and their defense systems can be co-opted when they provide an advantage to the cell. The abundance of defense systems in MGEs thus sheds new light on the role, effect, and fate of the so-called "cellular defense systems," whereby they are not only merely microbial defensive weapons in a 2-partner arms race, but also tools of intragenomic conflict between multiple genetic elements with divergent interests that shape cell fate and gene flow at the population level.

The role of interspecies recombination in the evolution of antibiotic-resistant pneumococci

Multidrug-resistant Streptococcus pneumoniae emerge through the modification of core genome loci by interspecies homologous recombinations, and acquisition of gene cassettes. Both occurred in the otherwise contrasting histories of the antibiotic-resistant S. pneumoniae lineages PMEN3 and PMEN9. A single PMEN3 clade spread globally, evading vaccine-induced immunity through frequent serotype switching, whereas locally circulating PMEN9 clades independently gained resistance. Both lineages repeatedly integrated Tn916-type and Tn1207.1-type elements, conferring tetracycline and macrolide resistance, respectively, through homologous recombination importing sequences originating in other species. A species-wide dataset found over 100 instances of such interspecific acquisitions of resistance cassettes and flanking homologous arms. Phylodynamic analysis of the most commonly sampled Tn1207.1-type insertion in PMEN9, originating from a commensal and disrupting a competence gene, suggested its expansion across Germany was driven by a high ratio of macrolide-to-β-lactam consumption. Hence, selection from antibiotic consumption was sufficient for these atypically large recombinations to overcome species boundaries across the pneumococcal chromosome.

Live cell dynamics of production, explosive release and killing activity of phage tail-like weapons for Pseudomonas kin exclusion

Interference competition among bacteria requires a highly specialized, narrow-spectrum weaponry when targeting closely-related competitors while sparing individuals from the same clonal population. Here we investigated mechanisms by which environmentally important Pseudomonas bacteria with plant-beneficial activity perform kin interference competition. We show that killing between phylogenetically closely-related strains involves contractile phage tail-like devices called R-tailocins that puncture target cell membranes. Using live-cell imaging, we evidence that R-tailocins are produced at the cell center, transported to the cell poles and ejected by explosive cell lysis. This enables their dispersal over several tens of micrometers to reach targeted cells. We visualize R-tailocin-mediated competition dynamics between closely-related Pseudomonas strains at the single-cell level, both in non-induced condition and upon artificial induction. We document the fatal impact of cellular self-sacrifice coupled to deployment of phage tail-like weaponry in the microenvironment of kin bacterial competitors, emphasizing the necessity for microscale assessment of microbial competitions.

Microbiome Engineering: Synthetic Biology of Plant-Associated Microbiomes in Sustainable Agriculture

To support an ever-increasing population, modern agriculture faces numerous challenges that pose major threats to global food and energy security. Plant-associated microbes, with their many plant growth-promoting (PGP) traits, have enormous potential in helping to solve these challenges. However, the results of their use in agriculture have been variable, probably because of poor colonization. Phytomicrobiome engineering is an emerging field of synthetic biology that may offer ways to alleviate this limitation. This review highlights recent advances in both bottom-up and top-down approaches to engineering non-model bacteria and microbiomes to promote beneficial plant-microbe interactions, as well as advances in strategies to evaluate these interactions. Biosafety, biosecurity, and biocontainment strategies to address the environmental concerns associated with field use of synthetic microbes are also discussed.

Modular prophage interactions driven by capsule serotype select for capsule loss under phage predation

Klebsiella species are able to colonize a wide range of environments and include worrisome nosocomial pathogens. Here, we sought to determine the abundance and infectivity of prophages of Klebsiella to understand how the interactions between induced prophages and bacteria affect population dynamics and evolution. We identified many prophages in the species, placing these taxa among the top 5% of the most polylysogenic bacteria. We selected 35 representative strains of the Klebsiella pneumoniae species complex to establish a network of induced phage-bacteria interactions. This revealed that many prophages are able to enter the lytic cycle, and subsequently kill or lysogenize closely related Klebsiella strains. Although 60% of the tested strains could produce phages that infect at least one other strain, the interaction network of all pairwise cross-infections is very sparse and mostly organized in modules corresponding to the strains' capsule serotypes. Accordingly, capsule mutants remain uninfected showing that the capsule is a key factor for successful infections. Surprisingly, experiments in which bacteria are predated by their own prophages result in accelerated loss of the capsule. Our results show that phage infectiousness defines interaction modules between small subsets of phages and bacteria in function of capsule serotype. This limits the role of prophages as competitive weapons because they can infect very few strains of the species complex. This should also restrict phage-driven gene flow across the species. Finally, the accelerated loss of the capsule in bacteria being predated by their own phages, suggests that phages drive serotype switch in nature.

The Ingenuity of Bacterial Genomes

The genomes of bacteria contain fewer genes and substantially less noncoding DNA than those of eukaryotes, and as a result, they have much less raw material to invent new traits. Yet, bacteria are vastly more taxonomically diverse, numerically abundant, and globally successful in colonizing new habitats compared to eukaryotes. Although bacterial genomes are generally considered to be optimized for efficient growth and rapid adaptation, nonadaptive processes have played a major role in shaping the size, contents, and compact organization of bacterial genomes and have allowed the establishment of deleterious traits that serve as the raw materials for genetic innovation.

Phylogenetic background and habitat drive the genetic diversification of Escherichia coli

Escherichia coli is mostly a commensal of birds and mammals, including humans, where it can act as an opportunistic pathogen. It is also found in water and sediments. We investigated the phylogeny, genetic diversification, and habitat-association of 1,294 isolates representative of the phylogenetic diversity of more than 5,000 isolates from the Australian continent. Since many previous studies focused on clinical isolates, we investigated mostly other isolates originating from humans, poultry, wild animals and water. These strains represent the species genetic diversity and reveal widespread associations between phylogroups and isolation sources. The analysis of strains from the same sequence types revealed very rapid change of gene repertoires in the very early stages of divergence, driven by the acquisition of many different types of mobile genetic elements. These elements also lead to rapid variations in genome size, even if few of their genes rise to high frequency in the species. Variations in genome size are associated with phylogroup and isolation sources, but the latter determine the number of MGEs, a marker of recent transfer, suggesting that gene flow reinforces the association of certain genetic backgrounds with specific habitats. After a while, the divergence of gene repertoires becomes linear with phylogenetic distance, presumably reflecting the continuous turnover of mobile element and the occasional acquisition of adaptive genes. Surprisingly, the phylogroups with smallest genomes have the highest rates of gene repertoire diversification and fewer but more diverse mobile genetic elements. This suggests that smaller genomes are associated with higher, not lower, turnover of genetic information. Many of these genomes are from freshwater isolates and have peculiar traits, including a specific capsule, suggesting adaptation to this environment. Altogether, these data contribute to explain why epidemiological clones tend to emerge from specific phylogenetic groups in the presence of pervasive horizontal gene transfer across the species.

Previous large scale studies on the evolution of E. coli focused on clinical isolates emphasizing virulence and antibiotic resistance in medically important lineages. Yet, most E. coli strains are either human commensals or not associated with humans at all. Here, we analyzed a large collection of non-clinical isolates of the species to assess the mechanisms of gene repertoire diversification in the light of isolation sources and phylogeny. We show that gene repertoires evolve so rapidly by the high turnover of mobile genetic elements that epidemiologically indistinguishable strains can be phenotypically extremely heterogeneous, illustrating the velocity of bacterial adaptation and the importance of accounting for the information on the whole genome at the epidemiological scale. Phylogeny and habitat shape the genetic diversification of E. coli to similar extents. Surprisingly, freshwater strains seem specifically adapted to this environment, breaking the paradigm that E. coli environmental isolates are systematically fecal contaminations. As a consequence, the evolution of this species is also shaped by environmental habitats, and it may diversify by acquiring genes and mobile elements from environmental bacteria (and not just from gut bacteria). This may facilitate the acquisition of virulence factors and antibiotic resistance in the strains that become pathogenic.

The Role of Integrative and Conjugative Elements in Antibiotic Resistance Evolution

Mobile genetic elements (MGEs), such as plasmids and integrative and conjugative elements (ICEs), are main drivers for the spread of antibiotic resistance (AR). Coevolution between bacteria and plasmids shapes the transfer and stability of plasmids across bacteria. Although ICEs outnumber conjugative plasmids, the dynamics of ICE-bacterium coevolution, ICE transfer rates, and fitness costs are as yet largely unexplored. Conjugative plasmids and ICEs are both transferred by type IV secretion systems, but ICEs are typically immune to segregational loss, suggesting that the evolution of ICE-bacterium associations varies from that of plasmid-bacterium associations. Considering the high abundance of ICEs among bacteria, ICE-bacterium dynamics represent a promising challenge for future research that will enhance our understanding of AR spread in human pathogens.

Early fate of exogenous promoters in E. coli

Gene gain by horizontal gene transfer is a major pathway of genome innovation in bacteria. The current view posits that acquired genes initially need to be silenced and that a bacterial chromatin protein, H-NS, plays a role in this silencing. However, we lack direct observation of the early fate of a horizontally transferred gene to prove this theory. We combine sequencing, flow cytometry and sorting, followed by microscopy to monitor gene expression and its variability after large-scale random insertions of a reporter gene in a population of Escherichia coli bacteria. We find that inserted promoters have a wide range of gene-expression variability related to their location. We find that high-expression clones carry insertions that are not correlated with H-NS binding. Conversely, binding of H-NS correlates with silencing. Finally, while most promoters show a common level of extrinsic noise, some insertions show higher noise levels. Analysis of these high-noise clones supports a scenario of switching due to transcriptional interference from divergent ribosomal promoters. Altogether, our findings point to evolutionary pathways where newly-acquired genes are not necessarily silenced, but may immediately explore a wide range of expression levels to probe the optimal ones.

Atypical organizations and epistatic interactions of CRISPRs and cas clusters in genomes and their mobile genetic elements

Prokaryotes use CRISPR-Cas systems for adaptive immunity, but the reasons for the frequent existence of multiple CRISPRs and cas clusters remain poorly understood. Here, we analysed the joint distribution of CRISPR and cas genes in a large set of fully sequenced bacterial genomes and their mobile genetic elements. Our analysis suggests few negative and many positive epistatic interactions between Cas subtypes. The latter often result in complex genetic organizations, where a locus has a single adaptation module and diverse interference mechanisms that might provide more effective immunity. We typed CRISPRs that could not be unambiguously associated with a cas cluster and found that such complex loci tend to have unique type I repeats in multiple CRISPRs. Many chromosomal CRISPRs lack a neighboring Cas system and they often have repeats compatible with the Cas systems encoded in trans. Phages and 25 000 prophages were almost devoid of CRISPR-Cas systems, whereas 3% of plasmids had CRISPR-Cas systems or isolated CRISPRs. The latter were often compatible with the chromosomal cas clusters, suggesting that plasmids can co-opt the latter. These results highlight the importance of interactions between CRISPRs and cas present in multiple copies and in distinct genomic locations in the function and evolution of bacterial immunity.

Synergistic Activity of Mobile Genetic Element Defences in Streptococcus pneumoniae

A diverse set of mobile genetic elements (MGEs) transmit between Streptococcus pneumoniae cells, but many isolates remain uninfected. The best-characterised defences against horizontal transmission of MGEs are restriction-modification systems (RMSs), of which there are two phase-variable examples in S. pneumoniae. Additionally, the transformation machinery has been proposed to limit vertical transmission of chromosomally integrated MGEs. This work describes how these mechanisms can act in concert. Experimental data demonstrate RMS phase variation occurs at a sub-maximal rate. Simulations suggest this may be optimal if MGEs are sometimes vertically inherited, as it reduces the probability that an infected cell will switch between RMS variants while the MGE is invading the population, and thereby undermine the restriction barrier. Such vertically inherited MGEs can be deleted by transformation. The lack of between-strain transformation hotspots at known prophage att sites suggests transformation cannot remove an MGE from a strain in which it is fixed. However, simulations confirmed that transformation was nevertheless effective at preventing the spread of MGEs into a previously uninfected cell population, if a recombination barrier existed between co-colonising strains. Further simulations combining these effects of phase variable RMSs and transformation found they synergistically inhibited MGEs spreading, through limiting both vertical and horizontal transmission.

May the Phage be With You? Prophage-Like Elements in the Genomes of Soft Rot Pectobacteriaceae: Pectobacterium spp. and Dickeya spp

Soft Rot Pectobacteriaceae (SRP; Pectobacterium spp. and Dickeya spp., formerly known as pectinolytic Erwinia spp.) are necrotrophic bacterial pathogens infecting a large number of plant species worldwide, including agriculturally-important crops. Despite the SRP importance in agriculture, little is known about the bacteriophages infecting them, and even less about the prophages present in their genomes. Prophages are recognized as factors underlying bacterial virulence, genomic diversification and ecological fitness that contribute to the novel phenotypic properties of bacterial hosts. Likewise, they are recognized as a driving force of bacterial evolution. In this study, 57 complete genomes of Pectobacterium spp. and Dickeya spp. deposited in NCBI GenBank, were analyzed for the presence of prophage-like elements. Viral sequences were discovered in 95% of bacterial genomes analyzed with the use of PHASTER, PhiSpy, and manual curation of the candidate sequences using NCBI BLAST. In total 37 seemingly intact and 48 putatively defective prophages were found. The 37 seemingly intact prophages (27 sequences in Dickeya spp. genomes and 10 sequences in Pectobacterium spp. genomes) were annotated using RAST. Analysis of the prophage genes encoding viral structural proteins allowed classification of these prophages into different families of the order Caudovirales (tailed bacteriophages) with the SRP prophages of the Myoviridae family (81% of found prophages) being the most abundant. The phylogenetic relationships between prophages were analyzed using amino acid sequences of terminase large subunit (gene terL), integrase (gene int), holin (gene hol), and lysin (gene lys). None of these markers however proved fully useful for clear phylogenetic separation of prophages of SRP into distinct clades. Comparative analyses of prophage proteomes revealed six clusters: five present in Dickeya spp. and one within Pectobacterium spp. When screened for the presence of bacterial genes in the genomes of intact prophages, only one prophage did not contain any ORFs of bacterial origin, the other prophages contained up to 23 genes acquired from bacterial hosts. The bacterial genes present in prophages could possibly affect fitness and virulence of their hosts. The implication of prophage presence in the genomes of Pectobacterium spp. and Dickeya spp. is discussed.

Into the Thermus Mobilome: Presence, Diversity and Recent Activities of Insertion Sequences Across Thermus spp

A high level of transposon-mediated genome rearrangement is a common trait among microorganisms isolated from thermal environments, probably contributing to the extraordinary genomic plasticity and horizontal gene transfer (HGT) observed in these habitats. In this work, active and inactive insertion sequences (ISs) spanning the sequenced members of the genus Thermus were characterized, with special emphasis on three T. thermophilus strains: HB27, HB8, and NAR1. A large number of full ISs and fragments derived from different IS families were found, concentrating within megaplasmids present in most isolates. Potentially active ISs were identified through analysis of transposase integrity, and domestication-related transposition events of ISTth7 were identified in laboratory-adapted HB27 derivatives. Many partial copies of ISs appeared throughout the genome, which may serve as specific targets for homologous recombination contributing to genome rearrangement. Moreover, recruitment of IS1000 32 bp segments as spacers for CRISPR sequence was identified, pointing to the adaptability of these elements in the biology of these thermophiles. Further knowledge about the activity and functional diversity of ISs in this genus may contribute to the generation of engineered transposons as new genetic tools, and enrich our understanding of the outstanding plasticity shown by these thermophiles.

The depths of virus exaptation

Viruses are ubiquitous parasites of cellular life forms and the most abundant biological entities on earth. The relationships between viruses and their hosts involve the continuous arms race but are by no account limited to it. Growing evidence shows that, in the course of evolution, viruses and their components are repeatedly recruited (exapted) for host functions. The functions of exapted viruses typically involve either defense from other viruses or cellular competitors or transfer of nucleic acids between cells, or storage functions. Virus exaptation can reach different depths, from recruitment of a fully functional virus to exploitation of defective, partially degraded viruses, to utilization of individual virus proteins.

Chromosome architecture constrains horizontal gene transfer in bacteria

Despite significant frequencies of lateral gene transfer between species, higher taxonomic groups of bacteria show ecological and phenotypic cohesion. This suggests that barriers prevent panmictic dissemination of genes via lateral gene transfer. We have proposed that most bacterial genomes have a functional architecture imposed by Architecture IMparting Sequences (AIMS). AIMS are defined as 8 base pair sequences preferentially abundant on leading strands, whose abundance and strand-bias are positively correlated with proximity to the replication terminus. We determined that inversions whose endpoints lie within a single chromosome arm, which would reverse the polarity of AIMS in the inverted region, are both shorter and less frequent near the replication terminus. This distribution is consistent with the increased selection on AIMS function in this region, thus constraining DNA rearrangement. To test the hypothesis that AIMS also constrain DNA transfer between genomes, AIMS were identified in genomes while ignoring atypical, potentially laterally-transferred genes. The strand-bias of AIMS within recently acquired genes was negatively correlated with the distance of those genes from their genome’s replication terminus. This suggests that selection for AIMS function prevents the acquisition of genes whose AIMS are not found predominantly in the permissive orientation. This constraint has led to the loss of at least 18% of genes acquired by transfer in the terminus-proximal region. We used completely sequenced genomes to produce a predictive road map of paths of expected horizontal gene transfer between species based on AIMS compatibility between donor and recipient genomes. These results support a model whereby organisms retain introgressed genes only if the benefits conferred by their encoded functions outweigh the detriments incurred by the presence of foreign DNA lacking genome-wide architectural information.

The potential success of horizontal gene transfer events is historically equated to the benefits conferred by encoded products. Here we show that gene transfer events are observed less frequently if the introduced genes disrupt important patterns of genomic information, suggesting that this disruption would confer an unacceptable cost. As a result, gene transfer events are less likely to be successful if the potential donor genomes have incompatible genome architecture. Because more distantly-related genes are less compatible, chromosome architecture serves as a mechanism to bias gene transfer events to those involving closer relatives, thereby providing a mechanism for the genotypic and phenotypic cohesion of higher taxonomic groups.

Assessing the benefits of horizontal gene transfer by laboratory evolution and genome sequencing

BACKGROUND

Recombination is widespread across the tree of life, because it helps purge deleterious mutations and creates novel adaptive traits. In prokaryotes, it often takes the form of horizontal gene transfer from a donor to a recipient bacterium. While such transfer is widespread in natural communities, its immediate fitness benefits are usually unknown. We asked whether any such benefits depend on the environment, and on the identity of donor and recipient strains. To this end, we adapted Escherichia coli to two novel carbon sources over several hundred generations of laboratory evolution, exposing evolving populations to various DNA donors.

RESULTS

At the end of these experiments, we measured fitness and sequenced the genomes of 65 clones from 34 replicate populations to study the genetic changes associated with adaptive evolution. Furthermore, we identified candidate de novo beneficial mutations. During adaptive evolution on the first carbon source, 4-Hydroxyphenylacetic acid (HPA), recombining populations adapted better, which was likely mediated by acquiring the hpa operon from the donor. In contrast, recombining populations did not adapt better to the second carbon source, butyric acid, even though they suffered fewer extinctions than non-recombining populations. The amount of DNA transferred, but not its benefit, strongly depended on the donor-recipient strain combination.

CONCLUSIONS

To our knowledge, our study is the first to investigate the genomic consequences of prokaryotic recombination and horizontal gene transfer during laboratory evolution. It shows that the benefits of recombination strongly depend on the environment and the foreign DNA donor.

Integrative and conjugative elements and their hosts: composition, distribution and organization

Conjugation of single-stranded DNA drives horizontal gene transfer between bacteria and was widely studied in conjugative plasmids. The organization and function of integrative and conjugative elements (ICE), even if they are more abundant, was only studied in a few model systems. Comparative genomics of ICE has been precluded by the difficulty in finding and delimiting these elements. Here, we present the results of a method that circumvents these problems by requiring only the identification of the conjugation genes and the species’ pan-genome. We delimited 200 ICEs and this allowed the first large-scale characterization of these elements. We quantified the presence in ICEs of a wide set of functions associated with the biology of mobile genetic elements, including some that are typically associated with plasmids, such as partition and replication. Protein sequence similarity networks and phylogenetic analyses revealed that ICEs are structured in functional modules. Integrases and conjugation systems have different evolutionary histories, even if the gene repertoires of ICEs can be grouped in function of conjugation types. Our characterization of the composition and organization of ICEs paves the way for future functional and evolutionary analyses of their cargo genes, composed of a majority of unknown function genes.

How a Genetically Stable Extremophile Evolves: Modes of Genome Diversification in the Archaeon Sulfolobus acidocaldarius

In order to analyze in molecular terms how Sulfolobus genomes diverge, damage-induced mutations and natural polymorphisms (PMs) were identified in laboratory constructs and wild-type isolates, respectively, of Sulfolobus acidocaldarius Among wild-type isolates drawn from one local population, pairwise nucleotide divergence averaged 4 × 10^-6, which is about 0.15% of the corresponding divergence reported for Sulfolobus islandicus The most variable features of wild-type S. acidocaldarius genomes were homopolymer (mononucleotide) tracts and longer tandem repeats, consistent with the spontaneous mutations that occur under laboratory conditions. Natural isolates, however, also revealed large insertions/deletions and inversions, which did not occur in any of the laboratory-manipulated strains. Several of the large insertions/deletions could be attributed to the integration or excision of mobile genetic elements (MGEs), and each MGE represented a distinct system of site-specific recombination. The mode of recombination associated with one MGE, a provirus related to Sulfolobus turreted icosahedral virus, was also seen in certain chromosomal inversions. Artificially induced mutations, non-MGE insertions/deletions, and small PMs exhibited different distributions over the genome, suggesting that large-scale patterning of Sulfolobus genomes begins early in the divergence process. Unlike induced mutations, natural base pair substitutions occurred in clusters, and one cluster exhibited properties expected of nonreciprocal recombination (gene conversion) between dispersed imperfect repeats. Taken together, the results identify simple replication errors, slipped-strand events promoted by tandem repeats, homologous recombination, and rearrangements promoted by MGEs as the primary sources of genetic variation for this extremely acidophilic archaeon in its geothermal environment.IMPORTANCE The optimal growth temperatures of hyperthermophilic archaea accelerate DNA decomposition, which is expected to make DNA repair especially important for their genetic stability, yet these archaea lack certain broadly conserved types of DNA repair proteins. In this study, the genome of the extreme thermoacidophile Sulfolobus acidocaldarius was found to be remarkably stable, accumulating few mutations in many (though not all) laboratory manipulations and in natural populations. Furthermore, all the genetic processes that were inferred to diversify these genomes also operate in mesophilic bacteria and eukaryotes. This suggests that a common set of mechanisms produces most of the genetic variation in all microorganisms, despite the fundamental differences in physiology, DNA repair systems, and genome structure represented in the three domains of life.

Crosstalk between vertical and horizontal gene transfer: plasmid replication control by a conjugative relaxase

Horizontal gene transfer is a key process in the evolution of bacteria and also represents a source of genetic variation in eukaryotes. Among elements participating in gene transfer, thousands of small (<10 kb) mobile bacterial plasmids that replicate by the rolling circle mechanism represent a driving force in the spread of antibiotic resistances. In general, these plasmids are built as genetic modules that encode a replicase, an antibiotic-resistance determinant, and a relaxase that participates in their conjugative mobilization. Further, they control their relatively high copy number (∼30 copies per genome equivalent) by antisense RNAs alone or combined with a repressor protein. We report here that the MobM conjugative relaxase encoded by the promiscuous plasmid pMV158 participates in regulation of the plasmid copy number by transcriptional repression of the antisense RNA, thus increasing the number of plasmid molecules ready to be horizontally transferred (mobilization) and/or vertically inherited (replication). This type of crosstalk between genetic modules involved in vertical and horizontal gene flow has not been reported before.

The hidden life of integrative and conjugative elements

Integrative and conjugative elements (ICEs) are widespread mobile DNA that transmit both vertically, in a host-integrated state, and horizontally, through excision and transfer to new recipients. Different families of ICEs have been discovered with more or less restricted host ranges, which operate by similar mechanisms but differ in regulatory networks, evolutionary origin and the types of variable genes they contribute to the host. Based on reviewing recent experimental data, we propose a general model of ICE life style that explains the transition between vertical and horizontal transmission as a result of a bistable decision in the ICE-host partnership. In the large majority of cells, the ICE remains silent and integrated, but hidden at low to very low frequencies in the population specialized host cells appear in which the ICE starts its process of horizontal transmission. This bistable process leads to host cell differentiation, ICE excision and transfer, when suitable recipients are present. The ratio of ICE bistability (i.e. ratio of horizontal to vertical transmission) is the outcome of a balance between fitness costs imposed by the ICE horizontal transmission process on the host cell, and selection for ICE distribution (i.e. ICE 'fitness'). From this emerges a picture of ICEs as elements that have adapted to a mostly confined life style within their host, but with a very effective and dynamic transfer from a subpopulation of dedicated cells.

Long-term evolution of viruses: A Janus-faced balance

The popular textbook image of viruses as noxious and selfish genetic parasites greatly underestimates the beneficial contributions of viruses to the biosphere. Given the crucial dependency of viruses to reproduce in an intracellular environment, viruses that engage in excessive killing (lysis) can drive their cellular hosts to extinction and will not survive. The lytic mode of virus propagation must, therefore, be tempered and balanced by non-lytic modes of virus latency and symbiosis. Here, we review recent bioinformatics and metagenomic studies to argue that viral endogenization and domestication may be more frequent mechanisms of virus persistence than lysis. We use a triangle diagram to explain the three major virus persistence strategies that explain the global scope of virus-cell interactions including lysis, latency and virus-cell symbiosis. This paradigm can help identify novel directions in virology research where scientists could artificially gain control over switching lytic and beneficial viral lifestyles. Also see the Video Abstract: http://youtu.be/GwXWz4N8o8.

Temperate bacteriophages as regulators of host behavior

Bacteriophages are ubiquitous and affect most facets of life, from evolution of bacteria, through ecology and global biochemical cycling to human health. The interactions between phages and bacteria often lead to biological novelty and an important milestone in this process is the ability of phages to regulate their host's behavior. In this review article, we will focus on newly reported cases that demonstrate how temperate phages regulate bacterial gene expression and behavior in a variety of bacterial species, pathogenic and environmental. This regulation is mediated by diverse mechanisms such as transcription factors, sRNAs, DNA rearrangements, and even controlled bacterial lysis. The outcome is mutualistic relationships that enable adaptively enhanced communal phage-host fitness under specific conditions.

Embracing the enemy: the diversification of microbial gene repertoires by phage-mediated horizontal gene transfer

Bacteriophages and archaeal viruses contribute, through lysogenic conversion or transduction, to the horizontal transfer of genetic material between microbial genomes. Recent genomics, metagenomics, and single cell studies have shown that lysogenic conversion is widespread and provides hosts with adaptive traits often associated with biotic interactions. The quantification of the evolutionary impact of transduction has lagged behind and requires further theoretical and experimental work. Nevertheless, recent studies suggested that generalized transduction plays a role in the transfer of antibiotic resistance genes and in the acquisition of novel genes during intra-specific bacterial competition. The characteristics of transduction and lysogenic conversion complement those of other mechanisms of transfer, and could play a key role in the spread of adaptive genes between communities.

In and out-contribution of natural transformation to the shuffling of large genomic regions

Naturally competent bacteria can pull free DNA from their surroundings. This incoming DNA can serve various purposes, ranging from acting as a source of nutrients or DNA stretches for repair to the acquisition of novel genetic information. The latter process defines the natural competence for transformation as a mode of horizontal gene transfer (HGT) and led to its discovery almost a century ago. However, although it is widely accepted that natural competence can contribute to the spread of genetic material among prokaryotes, the question remains whether this mode of HGT can foster the transfer of larger DNA regions or only transfers shorter fragments, given that extracellular DNA is often heavily fragmented. Here, I outline examples of competence-mediated movement of large genomic segments. Moreover, I discuss a recent proposition that transformation is used to cure bacteria of selfish mobile genetic elements. Such a transformation-mediated genome maintenance mechanism could indeed be an important and underappreciated function of natural competence.

The Streptococcus pneumoniae pezAT Toxin-Antitoxin System Reduces β-Lactam Resistance and Genetic Competence

Chromosomally encoded Type II Toxin–Antitoxin operons are ubiquitous in bacteria and archaea. Antitoxins neutralize the toxic effect of cognate Toxins by protein–protein interactions and sequestering the active residues of the Toxin. Toxins target essential bacterial processes, mostly translation and replication. However, one class apart is constituted by the PezAT pair because the PezT toxin target cell wall biosynthesis. Here, we have examined the role of the pezAT toxin–antitoxin genes in its natural host, the pathogenic bacterium Streptococcus pneumoniae. The pezAT operon on Pneumococcal Pathogenicity Island 1 was deleted from strain R6 and its phenotypic traits were compared with those of the wild type. The mutant cells formed shorter chains during exponential phase, leading to increased colony-forming units. At stationary phase, the mutant was more resilient to lysis. Importantly, the mutant exhibited higher resistance to antibiotics targeting cell walls (β-lactams), but not to antibiotics acting at other levels. In addition, the mutants also showed enhanced genetic competence. We suggest that PezAT participates in a subtle equilibrium between loss of functions (resistance to β-lactams and genetic competence) and gain of other traits (virulence).

IncA/C Conjugative Plasmids Mobilize a New Family of Multidrug Resistance Islands in Clinical Vibrio cholerae Non-O1/Non-O139 Isolates from Haiti

Mobile genetic elements play a pivotal role in the adaptation of bacterial populations, allowing them to rapidly cope with hostile conditions, including the presence of antimicrobial compounds. IncA/C conjugative plasmids (ACPs) are efficient vehicles for dissemination of multidrug resistance genes in a broad range of pathogenic species of Enterobacteriaceae. ACPs have sporadically been reported in Vibrio cholerae, the infectious agent of the diarrheal disease cholera. The regulatory network that controls ACP mobility ultimately depends on the transcriptional activation of multiple ACP-borne operons by the master activator AcaCD. Beyond ACP conjugation, AcaCD has also recently been shown to activate the expression of genes located in the Salmonella genomic island 1 (SGI1). Here, we describe MGIVchHai6, a novel and unrelated mobilizable genomic island (MGI) integrated into the 3′ end of trmE in chromosome I of V. cholerae HC-36A1, a non-O1/non-O139 multidrug-resistant clinical isolate recovered from Haiti in 2010. MGIVchHai6 contains a mercury resistance transposon and an integron In104-like multidrug resistance element similar to the one of SGI1. We show that MGIVchHai6 excises from the chromosome in an AcaCD-dependent manner and is mobilized by ACPs. Acquisition of MGIVchHai6 confers resistance to β-lactams, sulfamethoxazole, tetracycline, chloramphenicol, trimethoprim, and streptomycin/spectinomycin. In silico analyses revealed that MGIVchHai6-like elements are carried by several environmental and clinical V. cholerae strains recovered from the Indian subcontinent, as well as from North and South America, including all non-O1/non-O139 clinical isolates from Haiti.

Vibrio cholerae, the causative agent of cholera, remains a global public health threat. Seventh-pandemic V. cholerae acquired multidrug resistance genes primarily through circulation of SXT/R391 integrative and conjugative elements. IncA/C conjugative plasmids have sporadically been reported to mediate antimicrobial resistance in environmental and clinical V. cholerae isolates. Our results showed that while IncA/C plasmids are rare in V. cholerae populations, they play an important yet insidious role by specifically propagating a new family of genomic islands conferring resistance to multiple antibiotics. These results suggest that nonepidemic V. cholerae non-O1/non-O139 strains bearing these genomic islands constitute a reservoir of transmissible resistance genes that can be propagated by IncA/C plasmids to V. cholerae populations in epidemic geographical areas as well to pathogenic species of Enterobacteriaceae. We recommend future epidemiological surveys take into account the circulation of these genomic islands.

Functional and Evolutionary Characterization of a Gene Transfer Agent's Multilocus "Genome"

Gene transfer agents (GTAs) are phage-like particles that can package and transfer a random piece of the producing cell’s genome, but are unable to transfer all the genes required for their own production. As such, GTAs represent an evolutionary conundrum: are they selfish genetic elements propagating through an unknown mechanism, defective viruses, or viral structures “repurposed” by cells for gene exchange, as their name implies? In Rhodobacter capsulatus, production of the R. capsulatus GTA (RcGTA) particles is associated with a cluster of genes resembling a small prophage. Utilizing transcriptomic, genetic and biochemical approaches, we report that the RcGTA “genome” consists of at least 24 genes distributed across five distinct loci. We demonstrate that, of these additional loci, two are involved in cell recognition and binding and one in the production and maturation of RcGTA particles. The five RcGTA “genome” loci are widespread within Rhodobacterales, but not all loci have the same evolutionary histories. Specifically, two of the loci have been subject to frequent, probably virus-mediated, gene transfer events. We argue that it is unlikely that RcGTA is a selfish genetic element. Instead, our findings are compatible with the scenario that RcGTA is a virus-derived element maintained by the producing organism due to a selective advantage of within-population gene exchange. The modularity of the RcGTA “genome” is presumably a result of selection on the host organism to retain GTA functionality.

Using Sex to Cure the Genome

The diversification of prokaryotes is accelerated by their ability to acquire DNA from other genomes. However, the underlying processes also facilitate genome infection by costly mobile genetic elements. The discovery that cells can uptake DNA by natural transformation was instrumental to the birth of molecular biology nearly a century ago. Surprisingly, a new study shows that this mechanism could efficiently cure the genome of mobile elements acquired through previous sexual exchanges.

Natural transformation was thought to provide new genetic information to bacteria. Instead, a new study suggests it cures the genome of deleterious mobile elements.

The evolutionary dynamics of integrons in changing environments

Integrons are genetic elements that are common in bacteria and are hotspots for genome evolution. They facilitate the acquisition and reassembly of gene cassettes encoding a variety of functions, including drug resistance. Despite their importance in clinical settings, the selective forces responsible for the evolution and maintenance of integrons are poorly understood. We present a mathematical model of integron evolution within bacterial populations subject to fluctuating antibiotic exposures. Bacteria carrying a functional integrase that mediates reshuffling of cassette genes and thereby modulates gene expression patterns compete with bacteria without a functional integrase. Our results indicate that for a wide range of parameters, the functional integrase can be stably maintained in the population despite substantial fitness costs. This selective advantage arises because gene-cassette shuffling generates genetic diversity, thus enabling the population to respond rapidly to changing selective pressures. We also show that horizontal gene transfer promotes stable maintenance of the integrase and can also lead to de novo assembly of integrons. Our model generates testable predictions for integron evolution, including loss of functional integrases in stable environments and selection for intermediate gene-shuffling rates in changing environments. Our results highlight the need for experimental studies of integron population biology.

Coevolution of the Organization and Structure of Prokaryotic Genomes

The cytoplasm of prokaryotes contains many molecular machines interacting directly with the chromosome. These vital interactions depend on the chromosome structure, as a molecule, and on the genome organization, as a unit of genetic information. Strong selection for the organization of the genetic elements implicated in these interactions drives replicon ploidy, gene distribution, operon conservation, and the formation of replication-associated traits. The genomes of prokaryotes are also very plastic with high rates of horizontal gene transfer and gene loss. The evolutionary conflicts between plasticity and organization lead to the formation of regions with high genetic diversity whose impact on chromosome structure is poorly understood. Prokaryotic genomes are remarkable documents of natural history because they carry the imprint of all of these selective and mutational forces. Their study allows a better understanding of molecular mechanisms, their impact on microbial evolution, and how they can be tinkered in synthetic biology.

Co-occurrence of resistance genes to antibiotics, biocides and metals reveals novel insights into their co-selection potential

BACKGROUND

Antibacterial biocides and metals can co-select for antibiotic resistance when bacteria harbour resistance or tolerance genes towards both types of compounds. Despite numerous case studies, systematic and quantitative data on co-occurrence of such genes on plasmids and chromosomes is lacking, as is knowledge on environments and bacterial taxa that tend to carry resistance genes to such compounds. This effectively prevents identification of risk scenarios. Therefore, we aimed to identify general patterns for which biocide/metal resistance genes (BMRGs) and antibiotic resistance genes (ARGs) that tend to occur together. We also aimed to quantify co-occurrence of resistance genes in different environments and taxa, and investigate to what extent plasmids carrying both types of genes are conjugative and/or are carrying toxin-antitoxin systems.

RESULTS

Co-occurrence patterns of resistance genes were derived from publicly available, fully sequenced bacterial genomes (n = 2522) and plasmids (n = 4582). The only BMRGs commonly co-occurring with ARGs on plasmids were mercury resistance genes and the qacE∆1 gene that provides low-level resistance to quaternary ammonium compounds. Novel connections between cadmium/zinc and macrolide/aminoglycoside resistance genes were also uncovered. Several clinically important bacterial taxa were particularly prone to carry both BMRGs and ARGs. Bacteria carrying BMRGs more often carried ARGs compared to bacteria without (p < 0.0001). BMRGs were found in 86 % of bacterial genomes, and co-occurred with ARGs in 17 % of the cases. In contrast, co-occurrences of BMRGs and ARGs were rare on plasmids from all external environments (<0.7 %) but more common on those of human and domestic animal origin (5 % and 7 %, respectively). Finally, plasmids with both BMRGs and ARGs were more likely to be conjugative (p < 0.0001) and carry toxin-antitoxin systems (p < 0.0001) than plasmids without resistance genes.

CONCLUSIONS

This is the first large-scale identification of compounds, taxa and environments of particular concern for co-selection of resistance against antibiotics, biocides and metals. Genetic co-occurrences suggest that plasmids provide limited opportunities for biocides and metals to promote horizontal transfer of antibiotic resistance through co-selection, whereas ample possibilities exist for indirect selection via chromosomal BMRGs. Taken together, the derived patterns improve our understanding of co-selection potential between biocides, metals and antibiotics, and thereby provide guidance for risk-reducing actions.

Co-occurrence of resistance genes to antibiotics, biocides and metals reveals novel insights into their co-selection potential

BACKGROUND

Antibacterial biocides and metals can co-select for antibiotic resistance when bacteria harbour resistance or tolerance genes towards both types of compounds. Despite numerous case studies, systematic and quantitative data on co-occurrence of such genes on plasmids and chromosomes is lacking, as is knowledge on environments and bacterial taxa that tend to carry resistance genes to such compounds. This effectively prevents identification of risk scenarios. Therefore, we aimed to identify general patterns for which biocide/metal resistance genes (BMRGs) and antibiotic resistance genes (ARGs) that tend to occur together. We also aimed to quantify co-occurrence of resistance genes in different environments and taxa, and investigate to what extent plasmids carrying both types of genes are conjugative and/or are carrying toxin-antitoxin systems.

RESULTS

Co-occurrence patterns of resistance genes were derived from publicly available, fully sequenced bacterial genomes (n = 2522) and plasmids (n = 4582). The only BMRGs commonly co-occurring with ARGs on plasmids were mercury resistance genes and the qacE∆1 gene that provides low-level resistance to quaternary ammonium compounds. Novel connections between cadmium/zinc and macrolide/aminoglycoside resistance genes were also uncovered. Several clinically important bacterial taxa were particularly prone to carry both BMRGs and ARGs. Bacteria carrying BMRGs more often carried ARGs compared to bacteria without (p < 0.0001). BMRGs were found in 86 % of bacterial genomes, and co-occurred with ARGs in 17 % of the cases. In contrast, co-occurrences of BMRGs and ARGs were rare on plasmids from all external environments (<0.7 %) but more common on those of human and domestic animal origin (5 % and 7 %, respectively). Finally, plasmids with both BMRGs and ARGs were more likely to be conjugative (p < 0.0001) and carry toxin-antitoxin systems (p < 0.0001) than plasmids without resistance genes.

CONCLUSIONS

This is the first large-scale identification of compounds, taxa and environments of particular concern for co-selection of resistance against antibiotics, biocides and metals. Genetic co-occurrences suggest that plasmids provide limited opportunities for biocides and metals to promote horizontal transfer of antibiotic resistance through co-selection, whereas ample possibilities exist for indirect selection via chromosomal BMRGs. Taken together, the derived patterns improve our understanding of co-selection potential between biocides, metals and antibiotics, and thereby provide guidance for risk-reducing actions.

The dualistic nature of integrative and conjugative elements

Integrative and conjugative elements (ICEs) are mobile genetic elements that play a key role in bacterial adaptation. Such elements are found in almost every bacterial genera and species, and often code for adaptive traits conferring selective advantages to their host. ICEs maintain by integrating into and replicating along with a replicon of the host genome. ICEs can propagate by conjugative transfer toward a recipient cell following excision from the replicon as a circular covalently-closed molecule. For a long time, the excised form of ICEs was assumed to be non-replicative. This assumption predicts that excised ICEs are sensitive to loss during cell division, unless they carry stabilization systems such as addiction modules or antibiotic resistance genes. Over the past few years, growing evidence have been presented that support conditional replication of the circular intermediate as an intrinsic feature of ICEs. We recently confirmed this feature in the large family of SXT/R391 ICEs, which thrive in several species of Enterobacteriaceae and Vibrionaceae. Furthermore, we demonstrated that SXT/R391 ICEs encode a functional plasmid-like type II partition system that enhances their stability, such systems being probably encoded by other ICEs. The lifecycle of ICEs is therefore much more complex than initially thought as many ICEs may use plasmid-like features to improve their stability and dissemination.

Different Ancestries of R Tailocins in Rhizospheric Pseudomonas Isolates

Bacterial genomes accommodate a variety of mobile genetic elements, including bacteriophage-related clusters that encode phage tail-like protein complexes playing a role in interactions with eukaryotic or prokaryotic cells. Such tailocins are unable to replicate inside target cells due to the lack of a phage head with associated DNA. A subset of tailocins mediate antagonistic activities with bacteriocin-like specificity. Functional characterization of bactericidal tailocins of two Pseudomonas putida rhizosphere isolates revealed not only extensive similarity with the tail assembly module of the Pseudomonas aeruginosa R-type pyocins but also differences in genomic integration site, regulatory genes, and lytic release modules. Conversely, these three features are quite similar between strains of the P. putida and Pseudomonas fluorescens clades, although phylogenetic analysis of tail genes suggests them to have evolved separately. Unlike P. aeruginosa R pyocin elements, the tailocin gene clusters of other pseudomonads frequently carry cargo genes, including bacteriocins. Compared with P. aeruginosa, the tailocin tail fiber sequences that act as specificity determinants have diverged much more extensively among the other pseudomonad species, mostly isolates from soil and plant environments. Activity of the P. putida antibacterial particles requires a functional lipopolysaccharide layer on target cells, but contrary to R pyocins from P. aeruginosa, strain susceptibilities surpass species boundaries.

Independent Co-Option of a Tailed Bacteriophage into a Killing Complex in Pseudomonas

Competition between microbes is widespread in nature, especially among those that are closely related. To combat competitors, bacteria have evolved numerous protein-based systems (bacteriocins) that kill strains closely related to the producer. In characterizing the bacteriocin complement and killing spectra for the model strain Pseudomonas syringae B728a, we discovered that its activity was not linked to any predicted bacteriocin but is derived from a prophage. Instead of encoding an active prophage, this region encodes a bacteriophage-derived bacteriocin, termed an R-type syringacin. This R-type syringacin is striking in its convergence with the well-studied R-type pyocin of P. aeruginosa in both genomic location and molecular function. Genomic alignment, amino acid percent sequence identity, and phylogenetic inference all support a scenario where the R-type syringacin has been co-opted independently of the R-type pyocin. Moreover, the presence of this region is conserved among several other Pseudomonas species and thus is likely important for intermicrobial interactions throughout this important genus.

Evolutionary innovation is often achieved through modification of complexes or processes for alternate purposes, termed co-option. Notable examples include the co-option of a structure functioning in locomotion (bacterial flagellum) to one functioning in protein secretion (type three secretion system). Similar co-options can occur independently in distinct lineages. We discovered a genomic region in the plant pathogen Pseudomonas syringae that consists of a fragment of a bacteriophage genome. The fragment encodes only the tail of the bacteriophage, which is lethal toward strains of this species. This structure is similar to a previously described structure produced by the related species Pseudomonas aeruginosa. The two structures, however, are not derived from the same evolutionary event. Thus, they represent independent bacteriophage co-options. The co-opted bacteriophage from P. syringae is found in the genomes of many other Pseudomonas species, suggesting ecological importance across this genus.

Functionally Structured Genomes in Lactobacillus kunkeei Colonizing the Honey Crop and Food Products of Honeybees and Stingless Bees

Lactobacillus kunkeei is the most abundant bacterial species in the honey crop and food products of honeybees. The 16 S rRNA genes of strains isolated from different bee species are nearly identical in sequence and therefore inadequate as markers for studies of coevolutionary patterns. Here, we have compared the 1.5 Mb genomes of ten L. kunkeei strains isolated from all recognized Apis species and another two strains from Meliponini species. A gene flux analysis, including previously sequenced Lactobacillus species as outgroups, indicated the influence of reductive evolution. The genome architecture is unique in that vertically inherited core genes are located near the terminus of replication, whereas genes for secreted proteins and putative host-adaptive traits are located near the origin of replication. We suggest that these features have resulted from a genome-wide loss of genes, with integrations of novel genes mostly occurring in regions flanking the origin of replication. The phylogenetic analyses showed that the bacterial topology was incongruent with the host topology, and that strains of the same microcluster have recombined frequently across the host species barriers, arguing against codiversification. Multiple genotypes were recovered in the individual hosts and transfers of mobile elements could be demonstrated for strains isolated from the same host species. Unlike other bacteria with small genomes, short generation times and multiple rRNA operons suggest that L. kunkeei evolves under selection for rapid growth in its natural growth habitat. The results provide an extended framework for reductive genome evolution and functional genome organization in bacteria.