The Ecology and Evolution of Pangenomes

Pandora: nucleotide-resolution bacterial pan-genomics with reference graphs

We present pandora, a novel pan-genome graph structure and algorithms for identifying variants across the full bacterial pan-genome. As much bacterial adaptability hinges on the accessory genome, methods which analyze SNPs in just the core genome have unsatisfactory limitations. Pandora approximates a sequenced genome as a recombinant of references, detects novel variation and pan-genotypes multiple samples. Using a reference graph of 578 Escherichia coli genomes, we compare 20 diverse isolates. Pandora recovers more rare SNPs than single-reference-based tools, is significantly better than picking the closest RefSeq reference, and provides a stable framework for analyzing diverse samples without reference bias.

Gene-gene relationships in an Escherichia coli accessory genome are linked to function and mobility

The pangenome contains all genes encoded by a species, with the core genome present in all strains and the accessory genome in only a subset. Coincident gene relationships are expected within the accessory genome, where the presence or absence of one gene is influenced by the presence or absence of another. Here, we analysed the accessory genome of an Escherichia coli pangenome consisting of 400 genomes from 20 sequence types to identify genes that display significant co-occurrence or avoidance patterns with one another. We present a complex network of genes that are either found together or that avoid one another more often than would be expected by chance, and show that these relationships vary by lineage. We demonstrate that genes co-occur by function, and that several highly connected gene relationships are linked to mobile genetic elements. We find that genes are more likely to co-occur with, rather than avoid, another gene in the accessory genome. This work furthers our understanding of the dynamic nature of prokaryote pangenomes and implicates both function and mobility as drivers of gene relationships.

The dilution effect limits plasmid horizontal transmission in multispecies bacterial communities

By transferring ecologically important traits between species, plasmids drive genomic divergence and evolutionary innovation in their bacterial hosts. Bacterial communities are often diverse and contain multiple coexisting plasmids, but the dynamics of plasmids in multi-species communities are poorly understood. Here, we show, using experimental multi-species communities containing two plasmids, that bacterial diversity limits the horizontal transmission of plasmids due to the ‘dilution effect’; this is an epidemiological phenomenon whereby living alongside less proficient host species reduces the expected infection risk for a focal host species. In addition, plasmid horizontal transmission was also affected by plasmid diversity, such that the rate of plasmid conjugation was reduced from co-infected host cells carrying both plasmids. In diverse microbial communities, plasmid spread may be limited by the dilution effect and plasmid–plasmid interactions, reducing the rate of horizontal transmission.

The global population structure and beta-lactamase repertoire of the opportunistic pathogen Serratia marcescens

Serratia marcescens is a global spread nosocomial pathogen. This rod-shaped bacterium displays a broad host range and worldwide geographical distribution. Here we analyze an international collection of this multidrug-resistant, opportunistic pathogen from 35 countries to infer its population structure. We show that S. marcescens comprises 12 lineages; Sm1, Sm4, and Sm10 harbor 78.3% of the known environmental strains. Sm5, Sm6, and Sm7 comprise only human-associated strains which harbor smallest pangenomes, genomic fluidity and lowest levels of core recombination, indicating niche specialization. Sm7 and Sm9 lineages exhibit the most concerning resistome; bla_KPC-2 plasmid is widespread in Sm7, whereas Sm9, also an anthropogenic-exclusive lineage, presents highest plasmid/lineage size ratio and plasmid-diversity encoding metallo-beta-lactamases comprising bla_NDM-1. The heterogeneity of resistance patterns of S. marcescens lineages elucidated herein highlights the relevance of surveillance programs, using whole-genome sequencing, to provide insights into the molecular epidemiology of carbapenemase producing strains of this species.

Genetic innovations in animal-microbe symbioses

Animal hosts have initiated myriad symbiotic associations with microorganisms and often have maintained these symbioses for millions of years, spanning drastic changes in ecological conditions and lifestyles. The establishment and persistence of these relationships require genetic innovations on the parts of both symbionts and hosts. The nature of symbiont innovations depends on their genetic population structure, categorized here as open, closed or mixed. These categories reflect modes of inter-host transmission that result in distinct genomic features, or genomic syndromes, in symbionts. Although less studied, hosts also innovate in order to preserve and control symbiotic partnerships. New capabilities to sequence host-associated microbial communities and to experimentally manipulate both hosts and symbionts are providing unprecedented insights into how genetic innovations arise under different symbiont population structures and how these innovations function to support symbiotic relationships.

Integrated mass spectrometry-based multi-omics for elucidating mechanisms of bacterial virulence

Despite being considered the simplest form of life, bacteria remain enigmatic, particularly in light of pathogenesis and evolving antimicrobial resistance. After three decades of genomics, we remain some way from understanding these organisms, and a substantial proportion of genes remain functionally unknown. Methodological advances, principally mass spectrometry (MS), are paving the way for parallel analysis of the proteome, metabolome and lipidome. Each provides a global, complementary assay, in addition to genomics, and the ability to better comprehend how pathogens respond to changes in their internal (e.g. mutation) and external environments consistent with infection-like conditions. Such responses include accessing necessary nutrients for survival in a hostile environment where co-colonizing bacteria and normal flora are acclimated to the prevailing conditions. Multi-omics can be harnessed across temporal and spatial (sub-cellular) dimensions to understand adaptation at the molecular level. Gene deletion libraries, in conjunction with large-scale approaches and evolving bioinformatics integration, will greatly facilitate next-generation vaccines and antimicrobial interventions by highlighting novel targets and pathogen-specific pathways. MS is also central in phenotypic characterization of surface biomolecules such as lipid A, as well as aiding in the determination of protein interactions and complexes. There is increasing evidence that bacteria are capable of widespread post-translational modification, including phosphorylation, glycosylation and acetylation; with each contributing to virulence. This review focuses on the bacterial genotype to phenotype transition and surveys the recent literature showing how the genome can be validated at the proteome, metabolome and lipidome levels to provide an integrated view of organism response to host conditions.

Genome-Inferred Correspondence between Phylogeny and Metabolic Traits in the Wild Drosophila Gut Microbiome

Annotated genome sequences provide valuable insight into the functional capabilities of members of microbial communities. Nevertheless, most studies on the microbiome in animal guts use metagenomic data, hampering the assignment of genes to specific microbial taxa. Here, we make use of the readily culturable bacterial communities in the gut of the fruit fly Drosophila melanogaster to obtain draft genome sequences for 96 isolates from wild flies. These include 81 new de novo assembled genomes, assigned to three orders (Enterobacterales, Lactobacillales, and Rhodospirillales) with 80% of strains identified to species level using average nucleotide identity and phylogenomic reconstruction. Based on annotations by the RAST pipeline, among-isolate variation in metabolic function partitioned strongly by bacterial order, particularly by amino acid metabolism (Rhodospirillales), fermentation, and nucleotide metabolism (Lactobacillales) and arginine, urea, and polyamine metabolism (Enterobacterales). Seven bacterial species, comprising 2-3 species in each order, were well-represented among the isolates and included ≥5 strains, permitting analysis of metabolic functions in the accessory genome (i.e., genes not present in every strain). Overall, the metabolic function in the accessory genome partitioned by bacterial order. Two species, Gluconobacter cerinus (Rhodospirillales) and Lactiplantibacillus plantarum (Lactobacillales) had large accessory genomes, and metabolic functions were dominated by amino acid metabolism (G. cerinus) and carbohydrate metabolism (La. plantarum). The patterns of variation in metabolic capabilities at multiple phylogenetic scales provide the basis for future studies of the ecological and evolutionary processes shaping the diversity of microorganisms associated with natural populations of Drosophila.

Techniques in bacterial strain typing: past, present, and future

PURPOSE OF REVIEW

The advancement of molecular techniques such as whole-genome sequencing (WGS) has revolutionized the field of bacterial strain typing, with important implications for epidemiological surveillance and outbreak investigations. This review summarizes state-of-the-art techniques in strain typing and examines barriers faced by clinical and public health laboratories in implementing these new methodologies.

RECENT FINDINGS

WGS-based methodologies are on track to become the new 'gold standards' in bacterial strain typing, replacing traditional methods like pulsed-field gel electrophoresis and multilocus sequence typing. These new techniques have an improved ability to identify genetic relationships among organisms of interest. Further, advances in long-read sequencing approaches will likely provide a highly discriminatory tool to perform pangenome analyses and characterize relevant accessory genome elements, including mobile genetic elements carrying antibiotic resistance determinants in real time. Barriers to widespread integration of these approaches include a lack of standardized workflows and technical training.

SUMMARY

Genomic bacterial strain typing has facilitated a paradigm shift in clinical and molecular epidemiology. The increased resolution that these new techniques provide, along with epidemiological data, will facilitate the rapid identification of transmission routes with high confidence, leading to timely and effective deployment of infection control and public health interventions in outbreak settings.

CRISPR-Cas System: The Powerful Modulator of Accessory Genomes in Prokaryotes

Being frequently exposed to foreign nucleic acids, bacteria and archaea have developed an ingenious adaptive defense system, called CRISPR-Cas. The system is composed of the Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) array, together with CRISPR (cas)-associated genes. This system consists of a complex machinery that integrates fragments of foreign nucleic acids from viruses and mobile genetic elements (MGEs), into CRISPR arrays. The inserted segments (spacers) are transcribed and then used by cas proteins as guide RNAs for recognition and inactivation of the targets. Different types and families of CRISPR-Cas systems consist of distinct adaptation and effector modules with evolutionary trajectories, partially independent. The origin of the effector modules and the mechanism of spacer integration/deletion is far less clear. A review of the most recent data regarding the structure, ecology, and evolution of CRISPR-Cas systems and their role in the modulation of accessory genomes in prokaryotes is proposed in this article. The CRISPR-Cas system's impact on the physiology and ecology of prokaryotes, modulation of horizontal gene transfer events, is also discussed here. This system gained popularity after it was proposed as a tool for plant and animal embryo editing, in cancer therapy, as antimicrobial against pathogenic bacteria, and even for combating the novel coronavirus - SARS-CoV-2; thus, the newest and promising applications are reviewed as well.

Pangenomics reveals alternative environmental lifestyles among chlamydiae

Chlamydiae are highly successful strictly intracellular bacteria associated with diverse eukaryotic hosts. Here we analyzed metagenome-assembled genomes of the "Genomes from Earth's Microbiomes" initiative from diverse environmental samples, which almost double the known phylogenetic diversity of the phylum and facilitate a highly resolved view at the chlamydial pangenome. Chlamydiae are defined by a relatively large core genome indicative of an intracellular lifestyle, and a highly dynamic accessory genome of environmental lineages. We observe chlamydial lineages that encode enzymes of the reductive tricarboxylic acid cycle and for light-driven ATP synthesis. We show a widespread potential for anaerobic energy generation through pyruvate fermentation or the arginine deiminase pathway, and we add lineages capable of molecular hydrogen production. Genome-informed analysis of environmental distribution revealed lineage-specific niches and a high abundance of chlamydiae in some habitats. Together, our data provide an extended perspective of the variability of chlamydial biology and the ecology of this phylum of intracellular microbes.

Interplay of Various Evolutionary Modes in Genome Diversification and Adaptive Evolution of the Family Sulfolobaceae

Sulfolobaceae family, comprising diverse thermoacidophilic and aerobic sulfur-metabolizing Archaea from various geographical locations, offers an ideal opportunity to infer the evolutionary dynamics across the members of this family. Comparative pan-genomics coupled with evolutionary analyses has revealed asymmetric genome evolution within the Sulfolobaceae family. The trend of genome streamlining followed by periods of differential gene gains resulted in an overall genome expansion in some species of this family, whereas there was reduction in others. Among the core genes, both Sulfolobus islandicus and Saccharolobus solfataricus showed a considerable fraction of positively selected genes and also higher frequencies of gene acquisition. In contrast, Sulfolobus acidocaldarius genomes experienced substantial amount of gene loss and strong purifying selection as manifested by relatively lower genome size and higher genome conservation. Central carbohydrate metabolism and sulfur metabolism coevolved with the genome diversification pattern of this archaeal family. The autotrophic CO2 fixation with three significant positively selected enzymes from S. islandicus and S. solfataricus was found to be more imperative than heterotrophic CO2 fixation for Sulfolobaceae. Overall, our analysis provides an insight into the interplay of various genomic adaptation strategies including gene gain-loss, mutation, and selection influencing genome diversification of Sulfolobaceae at various taxonomic levels and geographical locations.

Multi-Omics Approach Reveals the Potential Core Vaccine Targets for the Emerging Foodborne Pathogen Campylobacter jejuni

Campylobacter jejuni is a leading cause of bacterial gastroenteritis in humans around the world. The emergence of bacterial resistance is becoming more serious; therefore, development of new vaccines is considered to be an alternative strategy against drug-resistant pathogen. In this study, we investigated the pangenome of 173 C. jejuni strains and analyzed the phylogenesis and the virulence factor genes. In order to acquire a high-quality pangenome, genomic relatedness was firstly performed with average nucleotide identity (ANI) analyses, and an open pangenome of 8,041 gene families was obtained with the correct taxonomy genomes. Subsequently, the virulence property of the core genome was analyzed and 145 core virulence factor (VF) genes were obtained. Upon functional genomics and immunological analyses, five core VF proteins with high antigenicity were selected as potential core vaccine targets for humans. Furthermore, functional annotations indicated that these proteins are involved in important molecular functions and biological processes, such as adhesion, regulation, and secretion. In addition, transcriptome analysis in human cells and pig intestinal loop proved that these vaccine target genes are important in the virulence of C. jejuni in different hosts. Comprehensive pangenome and relevant animal experiments will facilitate discovering the potential core vaccine targets with improved efficiency in reverse vaccinology. Likewise, this study provided some insights into the genetic polymorphism and phylogeny of C. jejuni and discovered potential vaccine candidates for humans. Prospective development of new vaccines using the targets will be an alternative to the use of antibiotics and prevent the development of multidrug-resistant C. jejuni in humans and even other animals.

Local adaptation, geographical distance and phylogenetic relatedness: Assessing the drivers of siderophore-mediated social interactions in natural bacterial communities

In heterogenous, spatially structured habitats, individuals within populations can become adapted to the prevailing conditions in their local environment. Such local adaptation has been reported for animals and plants, and for pathogens adapting to hosts. There is increasing interest in applying the concept of local adaptation to microbial populations, especially in the context of microbe-microbe interactions. Here, we tested whether cooperation and cheating on cooperation can spur patterns of local adaptation in soil and pond communities of Pseudomonas bacteria, collected across a geographical scale of 0.5 to 50 m. We focussed on the production of pyoverdines, a group of secreted iron-scavenging siderophores that often differ among pseudomonads in their chemical structure and the receptor required for their uptake. A combination of supernatant-feeding and competition assays between isolates from four distance categories revealed tremendous variation in the extent to which pyoverdine non- and low-producers can benefit from pyoverdines secreted by producers. However, this variation was not explained by geographical distance, but primarily depended on the phylogenetic relatedness between interacting isolates. A notable exception occurred in local pond communities, where the effect of phylogenetic relatedness was eroded in supernatant assays, probably due to the horizontal transfer of receptor genes. While the latter result could be a signature of local adaptation, our results overall indicate that common ancestry and not geographical distance is the main predictor of siderophore-mediated social interactions among pseudomonads.

Plant NLR diversity: the known unknowns of pan-NLRomes

Plants and pathogens constantly adapt to each other. As a consequence, many members of the plant immune system, and especially the intracellular nucleotide-binding site leucine-rich repeat receptors, also known as NOD-like receptors (NLRs), are highly diversified, both among family members in the same genome, and between individuals in the same species. While this diversity has long been appreciated, its true extent has remained unknown. With pan-genome and pan-NLRome studies becoming more and more comprehensive, our knowledge of NLR sequence diversity is growing rapidly, and pan-NLRomes provide powerful platforms for assigning function to NLRs. These efforts are an important step toward the goal of comprehensively predicting from sequence alone whether an NLR provides disease resistance, and if so, to which pathogens.

A tri-tuple coordinate system derived for fast and accurate analysis of the colored de Bruijn graph-based pangenomes

Background

With the rapid development of accurate sequencing and assembly technologies, an increasing number of high-quality chromosome-level and haplotype-resolved assemblies of genomic sequences have been derived, from which there will be great opportunities for computational pangenomics. Although genome graphs are among the most useful models for pangenome representation, their structural complexity makes it difficult to present genome information intuitively, such as the linear reference genome. Thus, efficiently and accurately analyzing the genome graph spatial structure and coordinating the information remains a substantial challenge.

Results

We developed a new method, a colored superbubble (cSupB), that can overcome the complexity of graphs and organize a set of species- or population-specific haplotype sequences of interest. Based on this model, we propose a tri-tuple coordinate system that combines an offset value, topological structure and sample information. Additionally, cSupB provides a novel method that utilizes complete topological information and efficiently detects small indels (< 50 bp) for highly similar samples, which can be validated by simulated datasets. Moreover, we demonstrated that cSupB can adapt to the complex cycle structure.

Conclusions

Although the solution is made suitable for increasingly complex genome graphs by relaxing the constraint, the directed acyclic graph, the motif cSupB and the cSupB method can be extended to any colored directed acyclic graph. We anticipate that our method will facilitate the analysis of individual haplotype variants and population genomic diversity. We have developed a C +  + program for implementing our method that is available at https://github.com/eggleader/cSupB.

Supplementary information

The online version contains supplementary material available at 10.1186/s12859-021-04149-w.

Microevolution in the pansecondary metabolome of Aspergillus flavus and its potential macroevolutionary implications for filamentous fungi

Fungi produce a wealth of pharmacologically bioactive secondary metabolites (SMs) from biosynthetic gene clusters (BGCs). It is common practice for drug discovery efforts to treat species' secondary metabolomes as being well represented by a single or a small number of representative genomes. However, this approach misses the possibility that intraspecific population dynamics, such as adaptation to environmental conditions or local microbiomes, may harbor novel BGCs that contribute to the overall niche breadth of species. Using 94 isolates of Aspergillus flavus, a cosmopolitan model fungus, sampled from seven states in the United States, we dereplicate 7,821 BGCs into 92 unique BGCs. We find that more than 25% of pangenomic BGCs show population-specific patterns of presence/absence or protein divergence. Population-specific BGCs make up most of the accessory-genome BGCs, suggesting that different ecological forces that maintain accessory genomes may be partially mediated by population-specific differences in secondary metabolism. We use ultra-high-performance high-resolution mass spectrometry to confirm that these genetic differences in BGCs also result in chemotypic differences in SM production in different populations, which could mediate ecological interactions and be acted on by selection. Thus, our results suggest a paradigm shift that previously unrealized population-level reservoirs of SM diversity may be of significant evolutionary, ecological, and pharmacological importance. Last, we find that several population-specific BGCs from A. flavus are present in Aspergillus parasiticus and Aspergillus minisclerotigenes and discuss how the microevolutionary patterns we uncover inform macroevolutionary inferences and help to align fungal secondary metabolism with existing evolutionary theory.

Positive Selection Inhibits Plasmid Coexistence in Bacterial Genomes

Plasmids play an important role in bacterial evolution by transferring niche-adaptive functional genes between lineages, thus driving genomic diversification. Bacterial genomes commonly contain multiple, coexisting plasmid replicons, which could fuel adaptation by increasing the range of gene functions available to selection and allowing their recombination. However, plasmid coexistence is difficult to explain because the acquisition of plasmids typically incurs high fitness costs for the host cell. Here, we show that plasmid coexistence was stably maintained without positive selection for plasmid-borne gene functions and was associated with compensatory evolution to reduce fitness costs. In contrast, with positive selection, plasmid coexistence was unstable despite compensatory evolution. Positive selection discriminated between differential fitness benefits of functionally redundant plasmid replicons, retaining only the more beneficial plasmid. These data suggest that while the efficiency of negative selection against plasmid fitness costs declines over time due to compensatory evolution, positive selection to maximize plasmid-derived fitness benefits remains efficient. Our findings help to explain the forces structuring bacterial genomes: coexistence of multiple plasmids in a genome is likely to require either rare positive selection in nature or nonredundancy of accessory gene functions among the coexisting plasmids.

Comparative phylo-pangenomics reveals generalist lifestyles in representative Acinetobacter species and proposes candidate gene markers for species identification

Acinetobacter species have the potential to invade and colonize immunocompromised patients, therefore being well-known as opportunistic pathogens. Among these bacteria, the species of the Acinetobacter calcoaceticus-Acinetobacter baumannii "complex" (Acb members) emerge as the main often isolated bacteria in clinical specimens. The unequivocal taxonomy is crucial to correctly identify these species and associated with comparative genomic analyses aids to understand their life-styles as well. In this study, all publicly available Acinetobacter species at the date of this study preparation were analyzed. The results revealed that the Acb members are in fact a complex when phenotypic methods are confronted, while for comparative and phylogenomics analyses this term is misleading, since they composed a monophyletic group instead. Nine best gene markers (response regulator, recJ, recG, phosphomannomutase, pepSY, monovalent cation/H + antiporter subunit D, mnmE, glnE, and bamA) were selected for identification of Acinetobacter species. Moreover, representative strains of each species were split according their isolation sources in the categories: environmental, human, insect and non-human vertebrate. Neither niche-specific genome signature nor niche-associated functional and pathogenic potential were associated with their isolation source, meaning it is not the main force acting on Acinetobacter adaptation in a given niche and corroborating that their ubiquitous distribution is a reflex of their generalist life-styles.

CRISPR-Cas systems restrict horizontal gene transfer in Pseudomonas aeruginosa

CRISPR-Cas systems provide bacteria and archaea with an adaptive immune system that targets foreign DNA. However, the xenogenic nature of immunity provided by CRISPR-Cas raises the possibility that these systems may constrain horizontal gene transfer. Here we test this hypothesis in the opportunistic pathogen Pseudomonas aeruginosa, which has emerged as an important model system for understanding CRISPR-Cas function. Across the diversity of P. aeruginosa, active CRISPR-Cas systems are associated with smaller genomes and higher GC content, suggesting that CRISPR-Cas inhibits the acquisition of foreign DNA. Although phage is the major target of CRISPR-Cas spacers, more than 80% of isolates with an active CRISPR-Cas system have spacers that target integrative conjugative elements (ICE) or the conserved conjugative transfer machinery used by plasmids and ICE. Consistent with these results, genomes containing active CRISPR-Cas systems harbour a lower abundance of both prophage and ICE. Crucially, spacers in genomes with active CRISPR-Cas systems map to ICE and phage that are integrated into the chromosomes of closely related genomes lacking CRISPR-Cas immunity. We propose that CRISPR-Cas acts as an important constraint to horizontal gene transfer, and the evolutionary mechanisms that ensure its maintenance or drive its loss are key to the ability of this pathogen to adapt to new niches and stressors.

Limited and Strain-Specific Transcriptional and Growth Responses to Acquisition of a Multidrug Resistance Plasmid in Genetically Diverse Escherichia coli Lineages

Multidrug-resistant (MDR) Escherichia coli strains are a major global threat to human health, wherein multidrug resistance is primarily spread by MDR plasmid acquisition. MDR plasmids are not widely distributed across the entire E. coli species, but instead are concentrated in a small number of clones. Here, we test if diverse E. coli strains vary in their ability to acquire and maintain MDR plasmids and if this relates to their transcriptional response following plasmid acquisition. We used strains from across the diversity of E. coli strains, including the common MDR lineage sequence type 131 (ST131) and the IncF plasmid pLL35, carrying multiple antibiotic resistance genes. Strains varied in their ability to acquire pLL35 by conjugation, but all were able to stably maintain the plasmid. The effects of pLL35 acquisition on cefotaxime resistance and growth also varied among strains, with growth responses ranging from a small decrease to a small increase in growth of the plasmid carrier relative to the parental strain. Transcriptional responses to pLL35 acquisition were limited in scale and highly strain specific. We observed transcriptional responses at the operon or regulon level-possibly due to stress responses or interactions with resident mobile genetic elements (MGEs). Subtle transcriptional responses consistent across all strains were observed affecting functions, such as anaerobic metabolism, previously shown to be under negative frequency-dependent selection in MDR E. coli Overall, there was no correlation between the magnitudes of the transcriptional and growth responses across strains. Together, these data suggest that fitness costs arising from transcriptional disruption are unlikely to act as a barrier to dissemination of this MDR plasmid in E. coli IMPORTANCE Plasmids play a key role in bacterial evolution by transferring adaptive functions between lineages that often enable invasion of new niches, including driving the spread of antibiotic resistance genes. Fitness costs of plasmid acquisition arising from the disruption of cellular processes could limit the spread of multidrug resistance plasmids. However, the impacts of plasmid acquisition are typically measured in lab-adapted strains rather than natural isolates, which act as reservoirs for the maintenance and transmission of plasmids to clinically relevant strains. Using a clinical multidrug resistance plasmid and a diverse collection of E. coli strains isolated from clinical infections and natural environments, we show that plasmid acquisition had only limited and highly strain-specific effects on bacterial growth and transcription under laboratory conditions. These findings suggest that fitness costs arising from transcriptional disruption are unlikely to act as a barrier to transmission of this plasmid in natural populations of E. coli.

The proficiency of the original host species determines community-level plasmid dynamics

Plasmids are common in natural bacterial communities, facilitating bacterial evolution via horizontal gene transfer. Bacterial species vary in their proficiency to host plasmids: whereas plasmids are stably maintained in some species regardless of selection for plasmid-encoded genes, in other species, even beneficial plasmids are rapidly lost. It is, however, unclear how this variation in host proficiency affects plasmid persistence in communities. Here, we test this using multispecies bacterial soil communities comprising species varying in their proficiency to host a large conjugative mercury resistance plasmid, pQBR103. The plasmid reached higher community-level abundance where beneficial and when introduced to the community in a more proficient host species. Proficient plasmid host species were also better able to disseminate the plasmid to a wider diversity of host species. These findings suggest that the dynamics of plasmids in natural bacterial communities depend not only upon the plasmid's attributes and the selective environment but also upon the proficiency of their host species.

Pangenome analyses of LuxS-coding genes and enzymatic repertoires in cocoa-related lactic acid bacteria

Lactobacillaceae presents potential for interspecific Quorum Sensing (QS) in spontaneous cocoa fermentation, correlated with high abundance of luxS. Three Brazilian isolates from cocoa fermentation were characterized by Whole Genome Sequencing and luxS gene was surveyed in their genomes, in comparison with public databases. They were classified as Lactiplantibacillus plantarum, Limosilactobacillus fermentum and Pediococcus acidilactici. LuxS genes were conserved in core genomes of the novel isolates, but in some non-cocoa related Lactic Acid Bacteria (LAB) it was accessory and plasmid-borne. The conservation and horizontal acquisition of luxS reinforces that QS is determinant for bacterial adaptation in several environments, especially taking into account the luxS has been correlated with modulation of bacteriocin production, stress tolerance and biofilm formation. Therefore, in this paper, new clade and species-specific primers were designed for future application for screening of luxS gene in LAB to evaluate the adaptive potential to diverse food fermentations.

Niche and local geography shape the pangenome of wastewater- and livestock-associated Enterobacteriaceae

Escherichia coli and other Enterobacteriaceae are diverse species with "open" pangenomes, where genes move intra- and interspecies via horizontal gene transfer. However, most analyses focus on clinical isolates. The pangenome dynamics of natural populations remain understudied, despite their suggested role as reservoirs for antimicrobial resistance (AMR) genes. Here, we analyze near-complete genomes for 827 Enterobacteriaceae (553 Escherichia and 274 non-Escherichia spp.) with 2292 circularized plasmids in total, collected from 19 locations (livestock farms and wastewater treatment works in the United Kingdom) within a 30-km radius at three time points over a year. We find different dynamics for chromosomal and plasmid-borne genes. Plasmids have a higher burden of AMR genes and insertion sequences, and AMR-gene-carrying plasmids show evidence of being under stronger selective pressure. Environmental niche and local geography both play a role in shaping plasmid dynamics. Our results highlight the importance of local strategies for controlling the spread of AMR.

Rhizobiales-Specific RirA Represses a Naturally "Synthetic" Foreign Siderophore Gene Cluster To Maintain Sinorhizobium-Legume Mutualism

Iron homeostasis is strictly regulated in cellular organisms. The Rhizobiales order enriched with symbiotic and pathogenic bacteria has evolved a lineage-specific regulator, RirA, responding to iron fluctuations. However, the regulatory role of RirA in bacterium-host interactions remains largely unknown. Here, we report that RirA is essential for mutualistic interactions of Sinorhizobium fredii with its legume hosts by repressing a gene cluster directing biosynthesis and transport of petrobactin siderophore. Genes encoding an inner membrane ABC transporter (fat) and the biosynthetic machinery (asb) of petrobactin siderophore are sporadically distributed in Gram-positive and Gram-negative bacteria. An outer membrane siderophore receptor gene (fprA) was naturally assembled with asb and fat, forming a long polycistron in S. fredii. An indigenous regulation cascade harboring an inner membrane protease (RseP), a sigma factor (FecI), and its anti-sigma protein (FecR) were involved in direct activation of the fprA-asb-fat polycistron. Operons harboring fecI and fprA-asb-fat, and those encoding the indigenous TonB-ExbB-ExbD complex delivering energy to the outer membrane transport activity, were directly repressed by RirA under iron-replete conditions. The rirA deletion led to upregulation of these operons and iron overload in nodules, impaired intracellular persistence, and symbiotic nitrogen fixation of rhizobia. Mutualistic defects of the rirA mutant can be rescued by blocking activities of this naturally "synthetic" circuit for siderophore biosynthesis and transport. These findings not only are significant for understanding iron homeostasis of mutualistic interactions but also provide insights into assembly and integration of foreign machineries for biosynthesis and transport of siderophores, horizontal transfer of which is selected in microbiota. IMPORTANCE Iron is a public good explored by both eukaryotes and prokaryotes. The abundant ferric form is insoluble under neutral and basic pH conditions, and many bacteria secrete siderophores forming soluble ferric siderophore complexes, which can be then taken up by specific receptors and transporters. Siderophore biosynthesis and uptake machineries can be horizontally transferred among bacteria in nature. Despite increasing attention on the importance of siderophores in host-microbiota interactions, the regulatory integration process of transferred siderophore biosynthesis and transport genes is poorly understood in an evolutionary context. By focusing on the mutualistic rhizobium-legume symbiosis, here, we report how a naturally synthetic foreign siderophore gene cluster was integrated with the rhizobial indigenous regulation cascade, which is essential for maintaining mutualistic interactions.

Genomic Background Governs Opposing Responses to Nalidixic Acid upon Megaplasmid Acquisition in Pseudomonas

Horizontal gene transfer is a significant driver of evolutionary dynamics across microbial populations. Although the benefits of the acquisition of new genetic material are often quite clear, experiments across systems have demonstrated that gene transfer events can cause significant phenotypic changes and entail fitness costs in a way that is dependent on the genomic and environmental context. Here, we test for the generality of one previously identified cost, sensitization of cells to the antibiotic nalidixic acid after acquisition of an ∼1-Mb megaplasmid, across Pseudomonas strains and species. Overall, we find that the presence of this megaplasmid sensitizes many different Pseudomonas strains to nalidixic acid but that this same horizontal gene transfer event increases resistance of Pseudomonas putida KT2440 to nalidixic acid across assays as well as to ciprofloxacin under competitive conditions. These phenotypic results are not easily explained away as secondary consequences of overall fitness effects and appear to occur independently of another cost associated with this megaplasmid, sensitization to higher temperatures. Lastly, we draw parallels between these reported results and the phenomenon of sign epistasis for de novo mutations and explore how context dependence of effects of plasmid acquisition could impact overall evolutionary dynamics and the evolution of antimicrobial resistance.IMPORTANCE Numerous studies have demonstrated that gene transfer events (e.g., plasmid acquisition) can entail a variety of costs that arise as by-products of the incorporation of foreign DNA into established physiological and genetic systems. These costs can be ameliorated through evolutionary time by the occurrence of compensatory mutations, which stabilize the presence of a horizontally transferred region within the genome but which also may skew future adaptive possibilities for these lineages. Here, we demonstrate another possible outcome, that phenotypic changes arising as a consequence of the same horizontal gene transfer (HGT) event are costly to some strains but may actually be beneficial in other genomic backgrounds under the right conditions. These results provide a new viewpoint for considering conditions that promote plasmid maintenance and highlight the influence of genomic and environmental contexts when considering amelioration of fitness costs after HGT events.

Pangenomics of the Symbiotic Rhizobiales. Core and Accessory Functions Across a Group Endowed with High Levels of Genomic Plasticity

Pangenome analyses reveal major clues on evolutionary instances and critical genome core conservation. The order Rhizobiales encompasses several families with rather disparate ecological attitudes. Among them, Rhizobiaceae, Bradyrhizobiaceae, Phyllobacteriacreae and Xanthobacteriaceae, include members proficient in mutualistic symbioses with plants based on the bacterial conversion of N₂ into ammonia (nitrogen-fixation). The pangenome of 12 nitrogen-fixing plant symbionts of the Rhizobiales was analyzed yielding total 37,364 loci, with a core genome constituting 700 genes. The percentage of core genes averaged 10.2% over single genomes, and between 5% to 7% were found to be plasmid-associated. The comparison between a representative reference genome and the core genome subset, showed the core genome highly enriched in genes for macromolecule metabolism, ribosomal constituents and overall translation machinery, while membrane/periplasm-associated genes, and transport domains resulted under-represented. The analysis of protein functions revealed that between 1.7% and 4.9% of core proteins could putatively have different functions.

Assessment of assumptions underlying models of prokaryotic pangenome evolution

BACKGROUND

The genomes of bacteria and archaea evolve by extensive loss and gain of genes which, for any group of related prokaryotic genomes, result in the formation of a pangenome with the universal, asymmetrical U-shaped distribution of gene commonality. However, the evolutionary factors that define the specific shape of this distribution are not thoroughly understood.

RESULTS

We investigate the fit of simple models of genome evolution to the empirically observed gene commonality distributions and genome intersections for 33 groups of closely related bacterial genomes. A model with an infinite external gene pool available for gene acquisition and constant genome size (IGP-CGS model), and two gene turnover rates, one for slow- and the other one for fast-evolving genes, allows two approaches to estimate the parameters for gene content dynamics. One is by fitting the model prediction to the distribution of the number of genes shared by precisely k genomes (gene commonality distribution) and another by analyzing the distribution of the number of genes common for k genome sets (k-cores). Both approaches produce a comparable overall quality of fit, although the former significantly overestimates the number of the universally conserved genes, while the latter overestimates the number of singletons. We further explore the effect of dropping each of the assumptions of the IGP-CGS model on the fit to the gene commonality distributions and show that models with either a finite gene pool or unequal rates of gene loss and gain (greater gene loss rate) eliminate the overestimate of the number of singletons or the core genome size.

CONCLUSIONS

We examine the assumptions that are usually adopted for modeling the evolution of the U-shaped gene commonality distributions in prokaryote genomes, namely, those of infinitely many genes and constant genome size. The combined analysis of genome intersections and gene commonality suggests that at least one of these assumptions is invalid. The violation of both these assumptions reflects the limited ability of prokaryotes to gain new genes. This limitation seems to stem, at least partly, from the horizontal gene transfer barrier, i.e., the cost of accommodation of foreign genes by prokaryotes. Further development of models taking into account the complexity of microbial evolution is necessary for an improved understanding of the evolution of prokaryotes.

Mechanisms That Shape Microbial Pangenomes

Analyses of multiple whole-genome sequences from the same species have revealed that differences in gene content can be substantial, particularly in prokaryotes. Such variation has led to the recognition of pangenomes, the complete set of genes present in a species - consisting of core genes, present in all individuals, and accessory genes whose presence is variable. Questions now arise about how pangenomes originate and evolve. We describe how gene content variation can arise as a result of the combination of several processes, including random drift, selection, gain/loss balance, and the influence of ecological and epistatic interactions. We believe that identifying the contributions of these processes to pangenomes will need novel theoretical approaches and empirical data.

Integration of diverse DNA substrates by a casposase can be targeted to R-loops in vitro by its fusion to Cas9

CRISPR systems build adaptive immunity against mobile genetic elements by DNA capture and integration catalysed by Cas1–Cas2 protein complexes. Recent studies suggested that CRISPR repeats and adaptation module originated from a novel type of DNA transposons called casposons. Casposons encode a Cas1 homologue called casposase that alone integrates into target molecules single and double-stranded DNA containing terminal inverted repeats (TIRs) from casposons. A recent study showed Methanosarcina mazei casposase is able to integrate random DNA oligonucleotides, followed up in this work using Acidoprofundum boonei casposase, from which we also observe promiscuous substrate integration. Here we first show that the substrate flexibility of Acidoprofundum boonei casposase extends to random integration of DNA without TIRs, including integration of a functional gene. We then used this to investigate targeting of the casposase-catalysed DNA integration reactions to specific DNA sites that would allow insertion of defined DNA payloads. Casposase–Cas9 fusions were engineered that were catalytically proficient in vitro and generated RNA-guided DNA integration products from short synthetic DNA or a gene, with or without TIRs. However, DNA integration could still occur unguided due to the competing background activity of the casposase moiety. Expression of Casposase-dCas9 in Escherichia coli cells effectively targeted chromosomal and plasmid lacZ revealed by reduced β-galactosidase activity but DNA integration was not detected. The promiscuous substrate integration properties of casposases make them potential DNA insertion tools. The Casposase–dCas9 fusion protein may serves as a prototype for development in genetic editing for DNA insertion that is independent of homology-directed DNA repair.

Coexistence of cry9 with the vip3A Gene in an Identical Plasmid of Bacillus thuringiensis Indicates Their Synergistic Insecticidal Toxicity

Bacillus thuringiensis (Bt) strains may express several insecticidal proteins with synergistic features, achieving high insecticidal toxicity and delaying development of resistance in insect pests. Previous work showed that Cry9Aa and Vip3Aa proteins present synergistic activity against Chilo suppressalis. In this study, genome-wide analysis of 489 Bt genomes revealed that cry9A was associated with the vip3A gene in seven Bt strains. Among all Bt genomes analyzed, not a single strain was found to have the cry9A gene alone without the presence of the vip3A gene. The complete genome sequencing of two Bt strains, 4AP1 and 4AO1, revealed that cry9A and vip3A genes were located in the same plasmid in both strains. The genome context analysis suggested a recombination mechanism responsible for the insertion of the cry9A gene into the plasmid containing vip3A. The coexistence of Cry9A with Vip3A proteins in strain 4AP1 was confirmed by liquid chromatography-tandem mass spectrometry and western blot analyses. Furthermore, another Cry9 protein codified by the gene in the identical plasmid also showed synergistic activity with the Vip3A protein. Overall, our results support that cry9 genes coexisted with vip3A and that complete genome sequencing combined with protein expression analysis may be used to identify associations of insecticidal proteins with potential synergistic toxicity.

Ecogenetics of antibiotic resistance in Listeria monocytogenes

The acquisition process of antibiotic resistance in an otherwise susceptible organism is shaped by the ecology of the species. Unlike other relevant human pathogens, Listeria monocytogenes has maintained a high rate of susceptibility to the antibiotics used for decades to treat human and animal infections. However, L. monocytogenes can acquire antibiotic resistance genes from other organisms' plasmids and conjugative transposons. Ecological factors could account for its susceptibility. L. monocytogenes is ubiquitous in nature, most frequently including reservoirs unexposed to antibiotics, including intracellular sanctuaries. L. monocytogenes has a remarkably closed genome, reflecting limited community interactions, small population sizes and high niche specialization. The L. monocytogenes species is divided into variants that are specialized in small specific niches, which reduces the possibility of coexistence with potential donors of antibiotic resistance. Interactions with potential donors are also hampered by interspecies antagonism. However, occasional increases in population sizes (and thus the possibility of acquiring antibiotic resistance) can derive from selection of the species based on intrinsic or acquired resistance to antibiotics, biocides, heavy metals or by a natural tolerance to extreme conditions. High-quality surveillance of the emergence of resistance to the key drugs used in primary therapy is mandatory.

A spectrum of verticality across genes

Lateral gene transfer (LGT) has impacted prokaryotic genome evolution, yet the extent to which LGT compromises vertical evolution across individual genes and individual phyla is unknown, as are the factors that govern LGT frequency across genes. Estimating LGT frequency from tree comparisons is problematic when thousands of genomes are compared, because LGT becomes difficult to distinguish from phylogenetic artefacts. Here we report quantitative estimates for verticality across all genes and genomes, leveraging a well-known property of phylogenetic inference: phylogeny works best at the tips of trees. From terminal (tip) phylum level relationships, we calculate the verticality for 19,050,992 genes from 101,422 clusters in 5,655 prokaryotic genomes and rank them by their verticality. Among functional classes, translation, followed by nucleotide and cofactor biosynthesis, and DNA replication and repair are the most vertical. The most vertically evolving lineages are those rich in ecological specialists such as Acidithiobacilli, Chlamydiae, Chlorobi and Methanococcales. Lineages most affected by LGT are the α-, β-, γ-, and δ- classes of Proteobacteria and the Firmicutes. The 2,587 eukaryotic clusters in our sample having prokaryotic homologues fail to reject eukaryotic monophyly using the likelihood ratio test. The low verticality of α-proteobacterial and cyanobacterial genomes requires only three partners-an archaeal host, a mitochondrial symbiont, and a plastid ancestor-each with mosaic chromosomes, to directly account for the prokaryotic origin of eukaryotic genes. In terms of phylogeny, the 100 most vertically evolving prokaryotic genes are neither representative nor predictive for the remaining 97% of an average genome. In search of factors that govern LGT frequency, we find a simple but natural principle: Verticality correlates strongly with gene distribution density, LGT being least likely for intruding genes that must replace a preexisting homologue in recipient chromosomes. LGT is most likely for novel genetic material, intruding genes that encounter no competing copy.

Gene Loss and Acquisition in Lineages of Pseudomonas aeruginosa Evolving in Cystic Fibrosis Patient Airways

Bacterial airway infections, predominantly caused by P. aeruginosa, are a major cause of mortality and morbidity of CF patients. While short insertions and deletions as well as point mutations occurring during infection are well studied, there is a lack of understanding of how gene loss and acquisition play roles in bacterial adaptation to the human airways. Here, we investigated P. aeruginosa within-host evolution with regard to gene loss and acquisition. We show that during long-term infection P. aeruginosa genomes tend to lose genes, in particular, genes related to virulence. This adaptive strategy allows reduction of the genome size and evasion of the host’s immune response. This knowledge is crucial to understand the basic mutational steps that, on the timescale of years, diversify lineages and adds to the identification of bacterial genetic determinants that have implications for CF disease.

Genome analyses have documented that there are differences in gene repertoire between evolutionary distant lineages of the same bacterial species; however, less is known about microevolutionary dynamics of gene loss and acquisition within bacterial lineages as they evolve over years. Here, we analyzed the genomes of 45 Pseudomonas aeruginosa lineages evolving in the lungs of cystic fibrosis (CF) patients to identify genes that are lost or acquired during the first years of infection. On average, lineage genome content changed by 88 genes (range, 0 to 473). Genes were more often lost than acquired, and prophage genes were more variable than bacterial genes. We identified convergent loss or acquisition of the same genes across lineages, suggesting selection for loss and acquisition of certain genes in the host environment. We found that a notable proportion of such genes are associated with virulence; a trait previously shown to be important for adaptation. Furthermore, we also compared the genomes across lineages to show that the within-lineage variable genes (i.e., genes that had been lost or acquired during the infection) often belonged to genomic content not shared across all lineages. In sum, our analysis adds to the knowledge on the pace and drivers of gene loss and acquisition in bacteria evolving over years in a human host environment and provides a basis to further understand how gene loss and acquisition play roles in lineage differentiation and host adaptation.

Unbiased homeologous recombination during pneumococcal transformation allows for multiple chromosomal integration events

The spread of antimicrobial resistance and vaccine escape in the human pathogen Streptococcus pneumoniae can be largely attributed to competence-induced transformation. Here, we studied this process at the single-cell level. We show that within isogenic populations, all cells become naturally competent and bind exogenous DNA. We find that transformation is highly efficient and that the chromosomal location of the integration site or whether the transformed gene is encoded on the leading or lagging strand has limited influence on recombination efficiency. Indeed, we have observed multiple recombination events in single recipients in real-time. However, because of saturation and because a single-stranded donor DNA replaces the original allele, transformation efficiency has an upper threshold of approximately 50% of the population. The fixed mechanism of transformation results in a fail-safe strategy for the population as half of the population generally keeps an intact copy of the original genome.

A Genome-Based Model to Predict the Virulence of Pseudomonas aeruginosa Isolates

Variation in the genome of Pseudomonas aeruginosa, an important pathogen, can have dramatic impacts on the bacterium's ability to cause disease. We therefore asked whether it was possible to predict the virulence of P. aeruginosa isolates based on their genomic content. We applied a machine learning approach to a genetically and phenotypically diverse collection of 115 clinical P. aeruginosa isolates using genomic information and corresponding virulence phenotypes in a mouse model of bacteremia. We defined the accessory genome of these isolates through the presence or absence of accessory genomic elements (AGEs), sequences present in some strains but not others. Machine learning models trained using AGEs were predictive of virulence, with a mean nested cross-validation accuracy of 75% using the random forest algorithm. However, individual AGEs did not have a large influence on the algorithm's performance, suggesting instead that virulence predictions are derived from a diffuse genomic signature. These results were validated with an independent test set of 25 P. aeruginosa isolates whose virulence was predicted with 72% accuracy. Machine learning models trained using core genome single-nucleotide variants and whole-genome k-mers also predicted virulence. Our findings are a proof of concept for the use of bacterial genomes to predict pathogenicity in P. aeruginosa and highlight the potential of this approach for predicting patient outcomes.IMPORTANCE Pseudomonas aeruginosa is a clinically important Gram-negative opportunistic pathogen. P. aeruginosa shows a large degree of genomic heterogeneity both through variation in sequences found throughout the species (core genome) and through the presence or absence of sequences in different isolates (accessory genome). P. aeruginosa isolates also differ markedly in their ability to cause disease. In this study, we used machine learning to predict the virulence level of P. aeruginosa isolates in a mouse bacteremia model based on genomic content. We show that both the accessory and core genomes are predictive of virulence. This study provides a machine learning framework to investigate relationships between bacterial genomes and complex phenotypes such as virulence.

The Impact of Mercury Selection and Conjugative Genetic Elements on Community Structure and Resistance Gene Transfer

Carriage of resistance genes can underpin bacterial survival, and by spreading these genes between species, mobile genetic elements (MGEs) can potentially protect diversity within microbial communities. The spread of MGEs could be affected by environmental factors such as selection for resistance, and biological factors such as plasmid host range, with consequences for individual species and for community structure. Here we cultured a focal bacterial strain, Pseudomonas fluorescens SBW25, embedded within a soil microbial community, with and without mercury selection, and with and without mercury resistance plasmids (pQBR57 or pQBR103), to investigate the effects of selection and resistance gene introduction on (1) the focal species; (2) the community as a whole; (3) the spread of the introduced mer resistance operon. We found that P. fluorescens SBW25 only escaped competitive exclusion by other members of community under mercury selection, even when it did not begin with a mercury resistance plasmid, due to its propensity to acquire resistance from the community by horizontal gene transfer. Mercury pollution had a significant effect on community structure, decreasing alpha diversity within communities while increasing beta diversity between communities, a pattern that was not affected by the introduction of mercury resistance plasmids by P. fluorescens SBW25. Nevertheless, the introduced merA gene spread to a phylogenetically diverse set of recipients over the 5 weeks of the experiment, as assessed by epicPCR. Our data demonstrates how the effects of MGEs can be experimentally assessed for individual lineages, the wider community, and for the spread of adaptive traits.

Bacterial dispersal and biogeography as underappreciated influences on phytobiomes

Bacterial strains are not distributed evenly throughout the environment. Here I explore how differential distribution and dispersal patterns of bacteria could affect interactions and coevolutionary dynamics with plants, and highlight ways that variation could be taken advantage of to develop robust and effective microbial consortia to inoculate crops. Questions about biogeographical patterns in viruses, fungi, and other eukaryotes are equally as prevalent and important for agriculture, and are in some cases more thoroughly explored. For simplicity as well as to bring attention to bacterial biogeography and dispersal in the context of plant interactions, I focus solely on bacterial patterns and questions for this article. The next few years will no doubt bring great advances in our understanding of dispersal capabilities and population dynamics for many plant-associated bacteria, and one of the next looming challenges will be learning to harvest this diversity in ways that can benefit agriculture.

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.

Genome Analysis of Endobacterium cerealis, a Novel Genus and Species Isolated from Zea mays Roots in North Spain

In the present work, we analyse the genomic and phenotypic characteristics of a strain named RZME27^T isolated from roots of a Zea mays plant grown in Spain. The phylogenetic analyses of 16S rRNA gene and whole genome sequences showed that the strain RZME27^T clustered with the type strains of Neorhizobium galegae and Pseudorhizobium pelagicum from the family Rhizobiaceae. This family encompasses several genera establishing symbiosis with legumes, but the genes involved in nodulation and nitrogen fixation are absent in its genome. Nevertheless, genes related to plant colonization, such as those involved in motility, chemotaxis, quorum sensing, exopolysaccharide biosynthesis and hydrolytic enzymes production were found. The comparative pangenomic analyses showed that 78 protein clusters present in the strain RZME27^T were not found in the type strains of its closest relatives N. galegae and P. pelagicum. The calculated average nucleotide identity (ANI) values between the strain RZME27^T and the type strains of N. galegae and P. pelagicum were 75.61% and 75.1%, respectively, similar or lower than those found for other genera from family Rhizobiaceae. Several phenotypic differences were also found, highlighting the absence of the fatty acid C_19:0 cyclo ω8c and propionate assimilation. These results support the definition of a novel genus and species named Endobacterium cerealis gen. nov. sp. nov. whose type strain is RZME27^T.

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.

Vast Differences in Strain-Level Diversity in the Gut Microbiota of Two Closely Related Honey Bee Species

Most bacterial species encompass strains with vastly different gene content. Strain diversity in microbial communities is therefore considered to be of functional importance. Yet little is known about the extent to which related microbial communities differ in diversity at this level and which underlying mechanisms may constrain and maintain strain-level diversity. Here, we used shotgun metagenomics to characterize and compare the gut microbiota of two honey bee species, Apis mellifera and Apis cerana, which diverged about 6 mya. Although the host species are colonized largely by the same bacterial 16S rRNA phylotypes, we find that their communities are host specific when analyzed with genomic resolution. Moreover, despite their similar ecology, A. mellifera displayed a much higher diversity of strains and functional gene content in the microbiota compared to A. cerana, both per colony and per individual bee. In particular, the gene repertoire for polysaccharide degradation was massively expanded in the microbiota of A. mellifera relative to A. cerana. Bee management practices, divergent ecological adaptation, or habitat size may have contributed to the observed differences in microbiota genomic diversity of these key pollinator species. Our results illustrate that the gut microbiota of closely related animal hosts can differ vastly in genomic diversity while displaying similar levels of diversity based on the 16S rRNA gene. Such differences are likely to have consequences for gut microbiota functioning and host-symbiont interactions, highlighting the need for metagenomic studies to understand the ecology and evolution of microbial communities.

•

Metagenomics reveals differences in gut microbiota diversity beyond the 16S rRNA gene

•

Apis cerana and Apis mellifera harbor distinct species and strains in their gut

•

Diversity is much higher in A. mellifera per individual bee and within colonies

•

Major differences in functions are related to polysaccharide degradation

Genomic variation among closely related Vibrio alginolyticus strains is located on mobile genetic elements

BACKGROUND

Species of the genus Vibrio, one of the most diverse bacteria genera, have undergone niche adaptation followed by clonal expansion. Niche adaptation and ultimately the formation of ecotypes and speciation in this genus has been suggested to be mainly driven by horizontal gene transfer (HGT) through mobile genetic elements (MGEs). Our knowledge about the diversity and distribution of Vibrio MGEs is heavily biased towards human pathogens and our understanding of the distribution of core genomic signatures and accessory genes encoded on MGEs within specific Vibrio clades is still incomplete. We used nine different strains of the marine bacterium Vibrio alginolyticus isolated from pipefish in the Kiel-Fjord to perform a multiscale-comparative genomic approach that allowed us to investigate [1] those genomic signatures that characterize a habitat-specific ecotype and [2] the source of genomic variation within this ecotype.

RESULTS

We found that the nine isolates from the Kiel-Fjord have a closed-pangenome and did not differ based on core-genomic signatures. Unique genomic regions and a unique repertoire of MGEs within the Kiel-Fjord isolates suggest that the acquisition of gene-blocks by HGT played an important role in the evolution of this ecotype. Additionally, we found that ~ 90% of the genomic variation among the nine isolates is encoded on MGEs, which supports ongoing theory that accessory genes are predominately located on MGEs and shared by HGT. Lastly, we could show that these nine isolates share a unique virulence and resistance profile which clearly separates them from all other investigated V. alginolyticus strains and suggests that these are habitat-specific genes, required for a successful colonization of the pipefish, the niche of this ecotype.

CONCLUSION

We conclude that all nine V. alginolyticus strains from the Kiel-Fjord belong to a unique ecotype, which we named the Kiel-alginolyticus ecotype. The low sequence variation of the core-genome in combination with the presence of MGE encoded relevant traits, as well as the presence of a suitable niche (here the pipefish), suggest, that this ecotype might have evolved from a clonal expansion following HGT driven niche-adaptation.

Coinfinder: detecting significant associations and dissociations in pangenomes

The accessory genes of prokaryote and eukaryote pangenomes accumulate by horizontal gene transfer, differential gene loss, and the effects of selection and drift. We have developed Coinfinder, a software program that assesses whether sets of homologous genes (gene families) in pangenomes associate or dissociate with each other (i.e. are 'coincident') more often than would be expected by chance. Coinfinder employs a user-supplied phylogenetic tree in order to assess the lineage-dependence (i.e. the phylogenetic distribution) of each accessory gene, allowing Coinfinder to focus on coincident gene pairs whose joint presence is not simply because they happened to appear in the same clade, but rather that they tend to appear together more often than expected across the phylogeny. Coinfinder is implemented in C++, Python3 and R and is freely available under the GNU license from https://github.com/fwhelan/coinfinder.

PPanGGOLiN: Depicting microbial diversity via a partitioned pangenome graph

The use of comparative genomics for functional, evolutionary, and epidemiological studies requires methods to classify gene families in terms of occurrence in a given species. These methods usually lack multivariate statistical models to infer the partitions and the optimal number of classes and don't account for genome organization. We introduce a graph structure to model pangenomes in which nodes represent gene families and edges represent genomic neighborhood. Our method, named PPanGGOLiN, partitions nodes using an Expectation-Maximization algorithm based on multivariate Bernoulli Mixture Model coupled with a Markov Random Field. This approach takes into account the topology of the graph and the presence/absence of genes in pangenomes to classify gene families into persistent, cloud, and one or several shell partitions. By analyzing the partitioned pangenome graphs of isolate genomes from 439 species and metagenome-assembled genomes from 78 species, we demonstrate that our method is effective in estimating the persistent genome. Interestingly, it shows that the shell genome is a key element to understand genome dynamics, presumably because it reflects how genes present at intermediate frequencies drive adaptation of species, and its proportion in genomes is independent of genome size. The graph-based approach proposed by PPanGGOLiN is useful to depict the overall genomic diversity of thousands of strains in a compact structure and provides an effective basis for very large scale comparative genomics. The software is freely available at https://github.com/labgem/PPanGGOLiN.

Toward a high-quality pan-genome landscape of Bacillus subtilis by removal of confounding strains

Pan-genome analysis is widely used to study the evolution and genetic diversity of species, particularly in bacteria. However, the impact of strain selection on the outcome of pan-genome analysis is poorly understood. Furthermore, a standard protocol to ensure high-quality pan-genome results is lacking. In this study, we carried out a series of pan-genome analyses of different strain sets of Bacillus subtilis to understand the impact of various strains on the performance and output quality of pan-genome analyses. Consequently, we found that the results obtained by pan-genome analyses of B. subtilis can be influenced by the inclusion of incorrectly classified Bacillus subspecies strains, phylogenetically distinct strains, engineered genome-reduced strains, chimeric strains, strains with a large number of unique genes or a large proportion of pseudogenes, and multiple clonal strains. Since the presence of these confounding strains can seriously affect the quality and true landscape of the pan-genome, we should remove these deviations in the process of pan-genome analyses. Our study provides new insights into the removal of biases from confounding strains in pan-genome analyses at the beginning of data processing, which enables the achievement of a closer representation of a high-quality pan-genome landscape of B. subtilis that better reflects the performance and credibility of the B. subtilis pan-genome. This procedure could be added as an important quality control step in pan-genome analyses for improving the efficiency of analyses, and ultimately contributing to a better understanding of genome function, evolution and genome-reduction strategies for B. subtilis in the future.

Population genomics of Klebsiella pneumoniae

Klebsiella pneumoniae is a common cause of antimicrobial-resistant opportunistic infections in hospitalized patients. The species is naturally resistant to penicillins, and members of the population often carry acquired resistance to multiple antimicrobials. However, knowledge of K. pneumoniae ecology, population structure or pathogenicity is relatively limited. Over the past decade, K. pneumoniae has emerged as a major clinical and public health threat owing to increasing prevalence of healthcare-associated infections caused by multidrug-resistant strains producing extended-spectrum β-lactamases and/or carbapenemases. A parallel phenomenon of severe community-acquired infections caused by 'hypervirulent' K. pneumoniae has also emerged, associated with strains expressing acquired virulence factors. These distinct clinical concerns have stimulated renewed interest in K. pneumoniae research and particularly the application of genomics. In this Review, we discuss how genomics approaches have advanced our understanding of K. pneumoniae taxonomy, ecology and evolution as well as the diversity and distribution of clinically relevant determinants of pathogenicity and antimicrobial resistance. A deeper understanding of K. pneumoniae population structure and diversity will be important for the proper design and interpretation of experimental studies, for interpreting clinical and public health surveillance data and for the design and implementation of novel control strategies against this important pathogen.

Disentangling the impact of environmental and phylogenetic constraints on prokaryotic within-species diversity

Microbial organisms inhabit virtually all environments and encompass a vast biological diversity. The pangenome concept aims to facilitate an understanding of diversity within defined phylogenetic groups. Hence, pangenomes are increasingly used to characterize the strain diversity of prokaryotic species. To understand the interdependence of pangenome features (such as the number of core and accessory genes) and to study the impact of environmental and phylogenetic constraints on the evolution of conspecific strains, we computed pangenomes for 155 phylogenetically diverse species (from ten phyla) using 7,000 high-quality genomes to each of which the respective habitats were assigned. Species habitat ubiquity was associated with several pangenome features. In particular, core-genome size was more important for ubiquity than accessory genome size. In general, environmental preferences had a stronger impact on pangenome evolution than phylogenetic inertia. Environmental preferences explained up to 49% of the variance for pangenome features, compared with 18% by phylogenetic inertia. This observation was robust when the dataset was extended to 10,100 species (59 phyla). The importance of environmental preferences was further accentuated by convergent evolution of pangenome features in a given habitat type across different phylogenetic clades. For example, the soil environment promotes expansion of pangenome size, while host-associated habitats lead to its reduction. Taken together, we explored the global principles of pangenome evolution, quantified the influence of habitat, and phylogenetic inertia on the evolution of pangenomes and identified criteria governing species ubiquity and habitat specificity.

The Role of Gram-Negative Bacteria in Urinary Tract Infections: Current Concepts and Therapeutic Options

Urinary tract infections (UTIs) are some of the most common infections in human medicine worldwide, recognized as an important public health concern to healthcare systems around the globe. In addition, urine specimens are one of the most frequently submitted samples for culture to the clinical microbiology laboratory, exceeding the number of most of the other sample types. The epidemiology, species-distribution and susceptibility-patterns of uropathogens vary greatly in a geographical and time-dependent manner and it also strongly correlated with the reported patient population studied. Nevertheless, many studies highlight the fact that the etiological agents in UTIs have changed considerably, both in nosocomial and community settings, with a shift towards "less common" microorganisms having more pronounced roles. There is increasing demand for further research to advance diagnostics and treatment options, and to improve care of the patients. The aim of this review paper was to summarize current developments in the global burden of UTI, the diagnostic aspects of these infectious pathologies, the possible etiological agents and their virulence determinants (with a special focus on the members of the Enterobacterales order), current guidelines and quality indicators in the therapy of UTIs and the emergence of multidrug resistance in urinary pathogens.