Bacterial Genome Evolution Research Articles

Comparative genomics analysis to explore the biodiversity and mining novel target genes of Listeria monocytogenes strains from different regions.

As a common foodborne pathogen, infection with L. monocytogenes poses a significant threat to human life and health. The objective of this study was to employ comparative genomics to unveil the biodiversity and evolutionary characteristics of L. monocytogenes strains from different regions, screening for potential target genes and mining novel target genes, thus providing significant reference value for the specific molecular detection and therapeutic targets of L. monocytogenes strains. Pan-genomic analysis revealed that L. monocytogenes from different regions have open genomes, providing a solid genetic basis for adaptation to different environments. These strains contain numerous virulence genes that contribute to their high pathogenicity. They also exhibit relatively high resistance to phosphonic acid, glycopeptide, lincosamide, and peptide antibiotics. The results of mobile genetic elements indicate that, despite being located in different geographical locations, there is a certain degree of similarity in bacterial genome evolution and adaptation to specific environmental pressures. The potential target genes identified through pan-genomics are primarily associated with the fundamental life activities and infection invasion of L. monocytogenes, including known targets such as inlB, which can be utilized for molecular detection and therapeutic purposes. After screening a large number of potential target genes, we further screened them using hub gene selection methods to mining novel target genes. The present study employed eight different hub gene screening methods, ultimately identifying ten highly connected hub genes (bglF_1, davD, menE_1, tilS, dapX, iolC, gshAB, cysG, trpA, and hisC), which play crucial roles in the pathogenesis of L. monocytogenes. The results of pan-genomic analysis showed that L. monocytogenes from different regions exhibit high similarity in bacterial genome evolution. The PCR results demonstrated the excellent specificity of the bglF_1 and davD genes for L. monocytogenes. Therefore, the bglF_1 and davD genes hold promise as specific molecular detection and therapeutic targets for L. monocytogenes strains from different regions.

Reversions mask the contribution of adaptive evolution in microbiomes.

When examining bacterial genomes for evidence of past selection, the results obtained depend heavily on the mutational distance between chosen genomes. Even within a bacterial species, genomes separated by larger mutational distances exhibit stronger evidence of purifying selection as assessed by , the normalized ratio of nonsynonymous to synonymous mutations. Here, we show that the classical interpretation of this scale-dependence, weak purifying selection, leads to problematic mutation accumulation when applied to available gut microbiome data. We propose an alternative, adaptive reversion model with exactly opposite implications for dynamical intuition and applications of . Reversions that occur and sweep within-host populations are nearly guaranteed in microbiomes due to large population sizes, short generation times, and variable environments. Using analytical and simulation approaches, we show that adaptive reversion can explain the decay given only dozens of locally-fluctuating selective pressures, which is realistic in the context of Bacteroides genomes. The success of the adaptive reversion model argues for interpreting low values of obtained from long-time scales with caution, as they may emerge even when adaptive sweeps are frequent. Our work thus inverts the interpretation of an old observation in bacterial evolution, illustrates the potential of mutational reversions to shape genomic landscapes over time, and highlights the importance of studying bacterial genomic evolution on short time scales.

Homologous Recombination Shapes the Architecture and Evolution of Bacterial Genomes.

Homologous recombination is a key evolutionary force that varies considerably across bacterial species. However, how the landscape of homologous recombination varies across genes and within individual genomes has only been studied in a few species. Here, we used Approximate Bayesian Computation to estimate the recombination rate along the genomes of 145 bacterial species. Our results show that homologous recombination varies greatly along bacterial genomes and shapes many aspects of genome architecture and evolution. The genomic landscape of recombination presents several key signatures: rates are highest near the origin of replication in most species, patterns of recombination generally appear symmetrical in both replichores (i.e. replicational halves of circular chromosomes) and most species have genomic hotpots of recombination. Furthermore, many closely related species share conserved landscapes of recombination across orthologs indicating that recombination landscapes are conserved over significant evolutionary distances. We show evidence that recombination drives the evolution of GC-content through increasing the effectiveness of selection and not through biased gene conversion, thereby contributing to an ongoing debate. Finally, we demonstrate that the rate of recombination varies across gene function and that many hotspots of recombination are associated with adaptive and mobile regions often encoding genes involved in pathogenicity.

Moving toward the Inclusion of Epigenomics in Bacterial Genome Evolution: Perspectives and Challenges.

The universality of DNA methylation as an epigenetic regulatory mechanism belongs to all biological kingdoms. However, while eukaryotic systems have been the primary focus of DNA methylation studies, the molecular mechanisms in prokaryotes are less known. Nevertheless, DNA methylation in prokaryotes plays a pivotal role in many cellular processes such as defense systems against exogenous DNA, cell cycle dynamics, and gene expression, including virulence. Thanks to single-molecule DNA sequencing technologies, genome-wide identification of methylated DNA is becoming feasible on a large scale, providing the possibility to investigate more deeply the presence, variability, and roles of DNA methylation. Here, we present an overview of the multifaceted roles of DNA methylation in prokaryotes and suggest research directions and tools which can enable us to better understand the contribution of DNA methylation to prokaryotic genome evolution and adaptation. In particular, we emphasize the need to understand the presence and role of transgenerational inheritance, as well as the impact of epigenomic signatures on adaptation and genome evolution. Research directions and the importance of novel computational tools are underlined.

Environmental Compatibility and Genome Flexibility of Klebsiella oxytoca Isolated from Eight Species of Aquatic Animals

Klebsiella oxytoca is an emerging pathogen that can cause life-threatening infectious diseases in humans. Recently, we firstly reported for the first time the presence of K. oxytoca in edible aquatic animals. In this study, we further investigated its bacterial environmental fitness and genome evolution signatures. The results revealed that K. oxytoca isolates (n = 8), originating from eight species of aquatic animals, were capable of growing under a broad spectrum of environmental conditions (pH 4.5–8.5, 0.5–6.5% NaCl), with different biofilm formation and swimming mobility profiles. The genome sequences of the K. oxytoca isolates were determined (5.84–6.02 Mb, 55.07–56.06% GC content). Strikingly, numerous putative mobile genetic elements (MGEs), particularly genomic islands (GIs, n = 105) and prophages (n = 24), were found in the K. oxytoca genomes, which provided the bacterium with specific adaptation traits, such as resistance, virulence, and material metabolism. Interestingly, the identified prophage-related clusters were derived from Burkholderia spp., Enterobacter spp., Klebsiella spp., Pseudomonas spp., and Haemophilus spp., suggesting phage transmission across Klebsiella and the other four genera. Many strain-specific (n = 10–447) genes were present in the K. oxytoca genomes, whereas the CRISPR-Cas protein-encoding gene was absent, indicating likely active horizontal gene transfer (HGT) and considerable genome variation in K. oxytoca evolution. Overall, the results of this study are the first to demonstrate the environmental compatibility and genome flexibility of K. oxytoca of aquatic animal origins.

Detecting horizontal gene transfer among microbiota: an innovative pipeline for identifying co-shared genes within the mobilome through advanced comparative analysis.

Horizontal gene transfer (HGT) is a key driver in the evolution of bacterial genomes. The acquisition of genes mediated by HGT may enable bacteria to adapt to ever-changing environmental conditions. Long-term application of antibiotics in intensive agriculture is associated with the dissemination of antibiotic resistance genes among bacteria with the consequences causing public health concern. Commensal farm-animal-associated gut microbiota are considered the reservoir of the resistance genes. Therefore, in this study, we identified known and not-yet characterized mobilized genes originating from chicken and porcine fecal samples using our innovative pipeline followed by network analysis to provide appropriate visualization to support proper interpretation.

Complex interactions between diverse mobile genetic elements drive the evolution of metal-resistant bacterial genomes.

In this study, we compared the genomes of three metal-resistant bacteria isolated from mercury-contaminated soil. We identified diverse and novel MGEs with evidence of multiple LGT events shaping their genomic structure and heavy metal resistance. Among the three metal-resistant strains, Sphingobium sp SA2 and Sphingopyxis sp SE2 were resistant to multiple metals including mercury, cadmium, copper, zinc and lead. Pseudoxanthomonas sp SE1 showed resistance to mercury only. Whole genome sequencing by Illumina and Oxford Nanopore technologies was undertaken to obtain comprehensive genomic data. The Sphingobium and Sphingopyxis strains contained multiple chromosomes and plasmids, whereas the Pseudoxanthomonas strain contained one circular chromosome. Consistent with their metal resistance profiles, the strains of Sphingobium and Sphingopyxis contained a higher quantity of diverse metal resistance genes across their chromosomes and plasmids compared to the single-metal resistant Pseudoxanthomonas SE1. In all three strains, metal resistance genes were principally associated with various novel MGEs including genomic islands (GIs), integrative conjugative elements (ICEs), transposons, insertion sequences (IS), recombinase in trio (RIT) elements and group II introns, indicating their importance in facilitating metal resistance adaptation in a contaminated environment. In the Pseudoxanthomonas strain, metal resistance regions were largely situated on a GI. The chromosomes of the strains of Sphingobium and Sphingopyxis contained multiple metal resistance regions, which were likely acquired by several GIs, ICEs, numerous IS elements, several Tn3 family transposons and RIT elements. Two of the plasmids of Sphingobium were impacted by Tn3 family transposons and ISs likely integrating metal resistance genes. The two plasmids of Sphingopyxis harboured transposons, IS elements, an RIT element and a group II intron. This study provides a comprehensive annotation of complex genomic regions of metal resistance associated with novel MGEs. It highlights the critical importance of LGT in the evolution of metal resistance of bacteria in contaminated environments.

Differential adenine methylation analysis reveals increased variability in 6mA in the absence of methyl-directed mismatch repair.

Methylation greatly influences the bacterial genome by guiding DNA repair and regulating pathogenic and stress-response phenotypes. But, the rate of epigenetic changes and their consequences on molecular phenotypes are underexplored. Through a detailed characterization of genome-wide adenine methylation in a commonly used laboratory strain of Escherichia coli, we reveal that mismatch repair deficient populations experience an increase in epimutations resulting in a genome-wide reduction of 6mA methylation in a manner consistent with genetic drift. Our findings highlight how methylation patterns evolve and the constraints on epigenetic evolution due to post-replicative DNA repair, contributing to a deeper understanding of bacterial genome evolution and how epimutations may introduce semi-permanent variation that can influence adaptation.

Genome content predicts the carbon catabolic preferences of heterotrophic bacteria.

Heterotrophic bacteria-bacteria that utilize organic carbon sources-are taxonomically and functionally diverse across environments. It is challenging to map metabolic interactions and niches within microbial communities due to the large number of metabolites that could serve as potential carbon and energy sources for heterotrophs. Whether their metabolic niches can be understood using general principles, such as a small number of simplified metabolic categories, is unclear. Here we perform high-throughput metabolic profiling of 186 marine heterotrophic bacterial strains cultured in media containing one of 135 carbon substrates to determine growth rates, lag times and yields. We show that, despite high variability at all levels of taxonomy, the catabolic niches of heterotrophic bacteria can be understood in terms of their preference for either glycolytic (sugars) or gluconeogenic (amino and organic acids) carbon sources. This preference is encoded by the total number of genes found in pathways that feed into the two modes of carbon utilization and can be predicted using a simple linear model based on gene counts. This allows for coarse-grained descriptions of microbial communities in terms of prevalent modes of carbon catabolism. The sugar-acid preference is also associated with genomic GC content and thus with the carbon-nitrogen requirements of their encoded proteome. Our work reveals how the evolution of bacterial genomes is structured by fundamental constraints rooted in metabolism.

Empirical evidence that complexity limits horizontal gene transfer.

Horizontal gene transfer (HGT) is a major contributor to bacterial genome evolution, generating phenotypic diversity, driving the expansion of protein families, and facilitating the evolution of new phenotypes, new metabolic pathways, and new species. Comparative studies of gene gain in bacteria suggest that the frequency with which individual genes successfully undergo HGT varies considerably and may be associated with the number of protein-protein interactions in which the gene participates, i.e., its connectivity. Two non-exclusive hypotheses have emerged to explain why transferability should decrease with connectivity: the Complexity Hypothesis (Jain, Rivera, & Lake, 1999) and the Balance Hypothesis (Papp, Pál, & Hurst, 2003). These hypotheses predict that the functional costs of HGT arise from a failure of divergent homologues to make normal protein-protein interactions or from gene mis-expression, respectively. Here we describe genome-wide assessments of these hypotheses in which we used 74 existing prokaryotic whole genome shotgun libraries to estimate rates of horizontal transfer of genes from taxonomically diverse prokaryotic donors into E. coli. We show that 1) transferability declines as connectivity increases, 2) transferability declines as the divergence between donor and recipient orthologs increases, and that 3) the magnitude of this negative effect of divergence on transferability increases with connectivity. These effects are particularly robust among the translational proteins, which span the widest range of connectivities. Whereas the Complexity Hypotheses explains all three of these observations, the Balance Hypothesis explains only the first one.

Survival and Genome Evolution Signatures of Klebsiella pneumoniae Isolates Originated in Seven Species of Aquatic Animals

Klebsiella pneumoniae can cause life-threatening pneumonia in humans. The bacterium is also the causative agent of nosocomial infection diseases. In our recent research, we reported, for the first time, the presence of K. pneumoniae in fourteen species of aquatic animals sampled in Shanghai, China. Here, we further investigated the bacterial survival and genome evolution traits. The results revealed that K. pneumoniae isolates (n = 7), recovered from 7 species of commonly consumed aquatic animals, had multiple antibiotic and heavy metal resistance profiles. The isolates were capable of growing vigorously at pH 4.5−7.5 and 0.5−1.0% NaCl in TSB medium at 37 °C. Draft genome sequences of the K. pneumoniae isolates were determined (5,256,522−5,857,823 bp, 56.35–57.81% GC contents), which carried many mobile genetic elements, including genomic islands (n = 87), prophages (n = 14), integrons (n = 4), and insertion sequences (n = 22), indicating possible active horizontal gene transfer during the genome evolution. Meanwhile, numerous strain-specific (n = 199−605) genes, antibiotic resistance (n = 20−35, e.g., β-lactamase) genes, and virulence (n = 43−59, e.g., enterobactin)-related genes, were also identified, demonstrating considerable genome variation in the K. pneumoniae isolates. Overall, the results of this study fill prior gaps in understanding the K. pneumoniae genomes derived from aquatic animals.

Reduction, evolutionary pattern and positive selection of genes encoding formate dehydrogenase in Wood-Ljungdahl pathway of gastrointestinal acetogens suggests their adaptation to formate-rich habitats.

Acetogens are anaerobes using Wood-Ljungdahl pathway (WLP) as the terminal electron acceptor for both assimilation and dissimilation of CO2 and widely distributed in diverse habitats. However, their habitat adaptation is often unclear. Given that bacterial genome evolution is often the result of environmental selective pressure, hereby we analysed gene copy number, phylogeny and selective pressure of genes involved in WLP within known genomes of 43 species to study the habitat adaption of gastrointestinal acetogens. The gene copy number of formate dehydrogenase (FDH) in gastrointestinal acetogens was much lower than that of non-gastrointestinal acetogens, and in five cases, no FDH genes were found in the genomes of five gastrointestinal acetogens, but that of the other WLP genes showed no difference. The evolutionary pattern of FDH genes was significantly different from that of the other enzymes. Additionally, seven positively selected sites were only identified in the fdhF genes, which means fdhF mutations favoured their adaptation. Collectively, reduction or loss of FDH genes and their evolutionary pattern as well as positive selection in gastrointestinal acetogens indicated their adaptation to formate-rich habitats, implying that FDH genes catalysing CO2 reduction to formate as the first step of methyl branch of WLP may have evolved independently.

Evolutionary Responses to Acquiring a Multidrug Resistance Plasmid Are Dominated by Metabolic Functions across Diverse Escherichia coli Lineages.

Multidrug resistance (MDR) plasmids drive the spread of antibiotic resistance between bacterial lineages. The immediate impact of MDR plasmid acquisition on fitness and cellular processes varies among bacterial lineages, but how the evolutionary processes enabling the genomic integration of MDR plasmids vary is less well understood, particularly in clinical pathogens. Using diverse Escherichia coli lineages experimentally evolved for ~700 generations, we show that the evolutionary response to gaining the MDR plasmid pLL35 was dominated by chromosomal mutations affecting metabolic and regulatory functions, with both strain-specific and shared mutational targets. The expression of several of these functions, such as anaerobic metabolism, is known to be altered upon acquisition of pLL35. Interactions with resident mobile genetic elements, notably several IS-elements, potentiated parallel mutations, including insertions upstream of hns that were associated with its upregulation and the downregulation of the plasmid-encoded extended-spectrum beta-lactamase gene. Plasmid parallel mutations targeted conjugation-related genes, whose expression was also commonly downregulated in evolved clones. Beyond their role in horizontal gene transfer, plasmids can be an important selective force shaping the evolution of bacterial chromosomes and core cellular functions. IMPORTANCE Plasmids drive the spread of antimicrobial resistance genes between bacterial genomes. However, the evolutionary processes allowing plasmids to be assimilated by diverse bacterial genomes are poorly understood, especially in clinical pathogens. Using experimental evolution with diverse E. coli lineages and a clinical multidrug resistance plasmid, we show that although plasmids drove unique evolutionary paths per lineage, there was a surprising degree of convergence in the functions targeted by mutations across lineages, dominated by metabolic functions. Remarkably, these same metabolic functions show higher evolutionary rates in MDR-lineages in nature and in some cases, like anaerobic metabolism, their expression is directly manipulated by the plasmid. Interactions with other mobile elements resident in the genomes accelerated adaptation by disrupting genes and regulatory sequences that they inserted into. Beyond their role in horizontal gene transfer, plasmids are an important selective force driving the evolution of bacterial genomes and core cellular functions.

Streptococcus pyogenes Φ1207.3 Is a Temperate Bacteriophage Carrying the Macrolide Resistance Gene Pair mef(A)-msr(D) and Capable of Lysogenizing Different Streptococci.

Streptococcus pyogenes prophage Φ1207.3 (formerly Tn1207.3) carries the mef(A)-msr(D) resistance genes, responsible for type M macrolide resistance. To investigate if Φ1207.3 is a functional bacteriophage, we transferred the element from the original S. pyogenes host in a prophage-free and competence-deficient S. pneumoniae strain. Pneumococcal cultures of the Φ1207.3-carrying lysogen were treated with mitomycin C to assess if Φ1207.3 enters the lytic cycle. Mitomycin C induced a limited phage burst and a growth impairment, resulting in early entrance into the stationary phase. To determine if Φ1207.3 is able to produce mature phage particles, we prepared concentrated supernatants recovered from a mitomycin C-induced pneumococcal culture by sequential centrifugation and ultracentrifugation steps. Negative-staining transmission electron microscopy (TEM) of supernatants revealed the presence of phage particles with an icosahedral, electron-dense capsid and a long, noncontractile tail, typical of a siphovirus. Quantification of Φ1207.3 was performed by quantitative PCR (qPCR) and semiquantitatively by TEM. PCR quantified 3.34 × 104 and 6.06 × 104 excised forms of phage genome per milliliter of supernatant obtained from the untreated and mitomycin C-treated cultures, respectively. By TEM, we estimated 3.02 × 103 and 7.68 × 103 phage particles per milliliter of supernatant. The phage preparations of Φ1207.3 infected and lysogenized pneumococcal recipient strains at a frequency of 7.5 × 10-6 lysogens/recipient but did not show sufficient lytic activity to form plaques. Phage lysogenization efficiently occurred after 30 min of contact of the phages with the recipient cells and required a minimum of 103 phage particles. IMPORTANCE Bacteriophages play an important role in bacterial physiology and genome evolution. The widespread use of genome sequencing revealed that bacterial genomes can contain several different integrated temperate bacteriophages, which can constitute up to 20% of the genome. Most of these bacteriophages are only predicted in silico and are never shown to be functional. In fact, it is often difficult to induce the lytic cycle of temperate bacteriophages. In this work, we show that Φ1207.3, a peculiar bacteriophage originally from Streptococcus pyogenes, which can lysogenize different streptococci and carries the macrolide resistance mef(A)-msr(D) gene pair, is capable of producing mature virions, but only at a low level, while not being able to produce plaques. This temperate phage is probably a partially functional phage, which seems to have lost lytic characteristics to specialize in lysogenization. While we are not used to conceiving phages separately from lysis, this behavior could actually be more frequent than expected.

Gene flow and introgression are pervasive forces shaping the evolution of bacterial species

BackgroundAlthough originally thought to evolve clonally, studies have revealed that most bacteria exchange DNA. However, it remains unclear to what extent gene flow shapes the evolution of bacterial genomes and maintains the cohesion of species.ResultsHere, we analyze the patterns of gene flow within and between >2600 bacterial species. Our results show that fewer than 10% of bacterial species are truly clonal, indicating that purely asexual species are rare in nature. We further demonstrate that the taxonomic criterion of ~95% genome sequence identity routinely used to define bacterial species does not accurately represent a level of divergence that imposes an effective barrier to gene flow across bacterial species. Interruption of gene flow can occur at various sequence identities across lineages, generally from 90 to 98% genome identity. This likely explains why a ~95% genome sequence identity threshold has empirically been judged as a good approximation to define bacterial species. Our results support a universal mechanism where the availability of identical genomic DNA segments required to initiate homologous recombination is the primary determinant of gene flow and species boundaries in bacteria. We show that these barriers of gene flow remain porous since many distinct species maintain some level of gene flow, similar to introgression in sexual organisms.ConclusionsOverall, bacterial evolution and speciation are likely shaped by similar forces driving the evolution of sexual organisms. Our findings support a model where the interruption of gene flow—although not necessarily the initial cause of speciation—leads to the establishment of permanent and irreversible species borders.

The effect of bacterial growth strategies on plasmid transfer and naphthalene degradation for bioremediation

Mobilizable plasmids are extra-chromosomal, circular DNA that have contributed to the rapid evolution of bacterial genomes and have been used in environmental, biotechnological, and medicinal applications. Degradative plasmids with genetic capabilities to degrade organic contaminants, such as polycyclic aromatic hydrocarbons (PAHs), have the potential to be useful for more environmentally friendly and cost-effective remediation technologies compared to existing physical remediation methods. Genetic bioaugmentation, the addition of catabolic genes into well-adapted communities via plasmid transfer (conjugation), is being explored as a remediation approach that is sustainable and long-lasting. Here, we explored the effect of the ecological growth strategies of plasmid donors and recipients on conjugation and naphthalene degradation of two PAH-degrading plasmids, pNL1 and NAH7. Overall, both pNL1 and NAH7 showed conjugation preferences towards a slow-growing ecological growth strategy, except when NAH7 was in a mixed synthetic community. These conjugation preferences were partially described by a combination of growth strategy, GC content, and phylogenetic relatedness. Further, removal of naphthalene via plasmid-mediated degradation was consistently higher in a community consisting of recipients with a slow-growing ecological growth strategy compared to a mixed community or a community consisting of fast-growing ecological growth strategy. Understanding plasmid conjugation and degradative preferences has the capacity to influence future remediation technology design and has broad implications in biomedical, environmental, and health fields.

Fitness costs of CRISPR-Cas systems in bacteria.

CRISPR-Cas systems provide bacteria with both specificity and adaptability in defence against invading genetic elements. From a theoretical perspective, CRISPR-Cas systems confer many benefits. However, they are observed at an unexpectedly low prevalence across the bacterial domain. While these defence systems can be gained horizontally, fitness costs may lead to selection against their carriage. Understanding the source of CRISPR-related fitness costs will help us to understand the evolutionary dynamics of CRISPR-Cas systems and their role in shaping bacterial genome evolution. Here, we review our current understanding of the potential fitness costs associated with CRISPR-Cas systems. In addition to potentially restricting the acquisition of genetic material that could confer fitness benefits, we explore five alternative biological factors that from a theoretical perspective may influence the fitness costs associated with CRISPR-Cas system carriage: (1) the repertoire of defence mechanisms a bacterium has available to it, (2) the potential for a metabolic burden, (3) larger-scale population and environmental factors, (4) the phenomenon of self-targeting spacers, and (5) alternative non-defence roles for CRISPR-Cas.

Dissemination of virulence and resistance genes among Klebsiella pneumoniae via outer membrane vesicle: An important plasmid transfer mechanism to promote the emergence of carbapenem-resistant hypervirulent Klebsiella pneumoniae.

Klebsiella pneumoniae is well-known opportunistic enterobacteria involved in complex clinical infections in humans and animals. The domestic animals might be a source of the multidrug-resistant virulent K. pneumoniae to humans. K. pneumoniae infections in domestic animals are considered as an emergent global concern. The horizontal gene transfer plays essential roles in bacterial genome evolution by spread of virulence and resistance determinants. However, the virulence genes can be transferred horizontally via K. pneumoniae-derived outer membrane vesicles (OMVs) remains to be unreported. In this study, we performed complete genome sequencing of two K. pneumoniae HvK2115 and CRK3022 with hypervirulent or carbapenem-resistant traits. OMVs from K. pneumoniae HvK2115 and CRK3022 were purified and observed. The carriage of virulence or resistance genes in K. pneumoniae OMVs was identified. The influence of OMVs on the horizontal transfer of virulence-related or drug-resistant plasmids among K. pneumoniae strains was evaluated thoroughly. The plasmid transfer to recipient bacteria through OMVs was identified by polymerase chain reaction, pulsed field gel electrophoresis and Southern blot. This study revealed that OMVs could mediate the intraspecific and interspecific horizontal transfer of the virulence plasmid phvK2115. OMVs could simultaneously transfer two resistance plasmids into K. pneumoniae and Escherichia coli recipient strains. OMVs-mediated horizontal transfer of virulence plasmid phvK2115 could significantly enhance the pathogenicity of human carbapenem-resistant K. pneumoniae CRK3022. The CRK3022 acquired the virulence plasmid phvK2115 could become a CR-hvKp strain. It was critically important that OMVs-mediated horizontal transfer of phvK2115 lead to the coexistence of virulence and carbapenem-resistance genes in K. pneumoniae, resulting in the emerging of carbapenem-resistant hypervirulent K. pneumoniae.

Pangenome Evolution in Environmentally Transmitted Symbionts of Deep-Sea Mussels Is Governed by Vertical Inheritance.

Microbial pangenomes vary across species; their size and structure are determined by genetic diversity within the population and by gene loss and horizontal gene transfer (HGT). Many bacteria are associated with eukaryotic hosts where the host colonization dynamics may impact bacterial genome evolution. Host-associated lifestyle has been recognized as a barrier to HGT in parentally transmitted bacteria. However, pangenome evolution of environmentally acquired symbionts remains understudied, often due to limitations in symbiont cultivation. Using high-resolution metagenomics, here we study pangenome evolution of two co-occurring endosymbionts inhabiting Bathymodiolus brooksi mussels from a single cold seep. The symbionts, sulfur-oxidizing (SOX) and methane-oxidizing (MOX) gamma-proteobacteria, are environmentally acquired at an early developmental stage and individual mussels may harbor multiple strains of each symbiont species. We found differences in the accessory gene content of both symbionts across individual mussels, which are reflected by differences in symbiont strain composition. Compared with core genes, accessory genes are enriched in genome plasticity functions. We found no evidence for recent HGT between both symbionts. A comparison between the symbiont pangenomes revealed that the MOX population is less diverged and contains fewer accessory genes, supporting that the MOX association with B. brooksi is more recent in comparison to that of SOX. Our results show that the pangenomes of both symbionts evolved mainly by vertical inheritance. We conclude that genome evolution of environmentally transmitted symbionts that associate with individual hosts over their lifetime is affected by a narrow symbiosis where the frequency of HGT is constrained.

Plasmid fitness costs are caused by specific genetic conflicts enabling resolution by compensatory mutation.

Plasmids play an important role in bacterial genome evolution by transferring genes between lineages. Fitness costs associated with plasmid carriage are expected to be a barrier to gene exchange, but the causes of plasmid fitness costs are poorly understood. Single compensatory mutations are often sufficient to completely ameliorate plasmid fitness costs, suggesting that such costs are caused by specific genetic conflicts rather than generic properties of plasmids, such as their size, metabolic burden, or gene expression level. By combining the results of experimental evolution with genetics and transcriptomics, we show here that fitness costs of 2 divergent large plasmids in Pseudomonas fluorescens are caused by inducing maladaptive expression of a chromosomal tailocin toxin operon. Mutations in single genes unrelated to the toxin operon, and located on either the chromosome or the plasmid, ameliorated the disruption associated with plasmid carriage. We identify one of these compensatory loci, the chromosomal gene PFLU4242, as the key mediator of the fitness costs of both plasmids, with the other compensatory loci either reducing expression of this gene or mitigating its deleterious effects by up-regulating a putative plasmid-borne ParAB operon. The chromosomal mobile genetic element Tn6291, which uses plasmids for transmission, remained up-regulated even in compensated strains, suggesting that mobile genetic elements communicate through pathways independent of general physiological disruption. Plasmid fitness costs caused by specific genetic conflicts are unlikely to act as a long-term barrier to horizontal gene transfer (HGT) due to their propensity for amelioration by single compensatory mutations, helping to explain why plasmids are so common in bacterial genomes.