Diversification of bacterial genome content through distinct mechanisms over different timescales

Comprehensive analysis of genomic variation, pan-genome and biosynthetic potential of Corynebacterium glutamicum strains

Corynebacterium glutamicum is a non-pathogenic species of the Corynebacteriaceae family. It has been broadly used in industrial biotechnology for the production of valuable products. Though it is widely accepted at the industrial level, knowledge about the genomic diversity of the strains is limited. Here, we investigated the comparative genomic features of the strains and pan-genomic characteristics. We also observed phylogenetic relationships among the strains based on average nucleotide identity (ANI). We found diversity between strains at the genomic and pan-genomic levels. Less than one-third of the C. glutamicum pan-genome consists of core genes and soft-core genes. Whereas, a large number of strain-specific genes covered about half of the total pan-genome. Besides, C. glutamicum pan-genome is open and expanding, which indicates the possible addition of new gene families to the pan-genome. We also investigated the distribution of biosynthetic gene clusters (BGCs) among the strains. We discovered slight variations of BGCs at the strain level. Several BGCs with the potential to express novel bioactive secondary metabolites have been identified. Therefore, by utilizing the characteristic advantages of C. glutamicum, different strains can be potential applicants for natural drug discovery.

Analysis of the limited M. tuberculosis accessory genome reveals potential pitfalls of pan-genome analysis approaches

Pan-genome analysis is a fundamental tool for studying bacterial genome evolution; however, the variety of methods used to define and measure the pan-genome poses challenges to the interpretation and reliability of results. To quantify sources of bias and error related to common pan-genome analysis approaches, we evaluated different approaches applied to curated collection of 151 Mycobacterium tuberculosis ( Mtb ) isolates. Mtb is characterized by its clonal evolution, absence of horizontal gene transfer, and limited accessory genome, making it an ideal test case for this study. Using a state-of-the-art graph-genome approach, we found that a majority of the structural variation observed in Mtb originates from rearrangement, deletion, and duplication of redundant nucleotide sequences. In contrast, we found that pan-genome analyses that focus on comparison of coding sequences (at the amino acid level) can yield surprisingly variable results, driven by differences in assembly quality and the softwares used. Upon closer inspection, we found that coding sequence annotation discrepancies were a major contributor to inflated Mtb accessory genome estimates. To address this, we developed panqc, a software that detects annotation discrepancies and collapses nucleotide redundancy in pan-genome estimates. When applied to Mtb and E. coli pan-genomes, panqc exposed distinct biases influenced by the genomic diversity of the population studied. Our findings underscore the need for careful methodological selection and quality control to accurately map the evolutionary dynamics of a bacterial species.

Evolutionary dynamics of the accessory genomes of Staphylococcus aureus

Staphylococcus aureus is a ubiquitous commensal and opportunistic bacterial pathogen that can cause a wide gamut of infections, which are exacerbated by the presence of multidrug-resistant and methicillin-resistant S. aureus. S. aureus is genetically heterogeneous and consists of numerous distinct lineages. Using 558 complete genomes of S. aureus, we aim to determine how the accessory genome content among phylogenetic lineages of S. aureus is structured and has evolved. Bayesian hierarchical clustering identified 10 sequence clusters, of which seven contained major sequence types (ST 1, 5, 8, 30, 59, 239, and 398). The seven sequence clusters differed in their accessory gene content, including genes associated with antimicrobial resistance and virulence. Focusing on the two largest clusters, BAPS8 and BAPS10, and each consisting mostly of ST5 and ST8, respectively, we found that the structure and connected components in the co-occurrence networks of accessory genomes varied between them. These differences are explained, in part, by the variation in the rates at which the two sequence clusters gained and lost accessory genes, with the highest rate of gene accumulation occurring recently in their evolutionary histories. We also identified a divergent group within BAPS10 that has experienced high gene gain and loss early in its history. Together, our results show highly variable and dynamic accessory genomes in S. aureus that are structured by the history of the specific lineages that carry them.IMPORTANCEStaphylococcus aureus is an opportunistic, multi-host pathogen that can cause a variety of benign and life-threatening infections. Our results revealed considerable differences in the structure and evolution of the accessory genomes of major lineages within S. aureus. Such genomic variation within a species can have important implications on disease epidemiology, pathogenesis of infection, and interactions with the vertebrate host. Our findings provide important insights into the underlying genetic basis for the success of S. aureus as a highly adaptable and resistant pathogen, which will inform current efforts to control and treat staphylococcal diseases.

The metabolic, virulence and antimicrobial resistance profiles of colonising Streptococcus pneumoniae shift after PCV13 introduction in urban Malawi

Streptococcus pneumoniae causes substantial mortality among children under 5-years-old worldwide. Polysaccharide conjugate vaccines (PCVs) are highly effective at reducing vaccine serotype disease, but emergence of non-vaccine serotypes and persistent nasopharyngeal carriage threaten this success. We investigated the hypothesis that following vaccine, adapted pneumococcal genotypes emerge with the potential for vaccine escape. We genome sequenced 2804 penumococcal isolates, collected 4-8 years after introduction of PCV13 in Blantyre, Malawi. We developed a pipeline to cluster the pneumococcal population based on metabolic core genes into "Metabolic genotypes" (MTs). We show that S. pneumoniae population genetics are characterised by emergence of MTs with distinct virulence and antimicrobial resistance (AMR) profiles. Preliminary in vitro and murine experiments revealed that representative isolates from emerging MTs differed in growth, haemolytic, epithelial infection, and murine colonisation characteristics. Our results suggest that in the context of PCV13 introduction, pneumococcal population dynamics had shifted, a phenomenon that could further undermine vaccine control and promote spread of AMR.

Diverse regulatory pathways modulate bet hedging of competence induction in epigenetically-differentiated phase variants of Streptococcus pneumoniae

Despite enabling Streptococcus pneumoniae to acquire antibiotic resistance and evade vaccine-induced immunity, transformation occurs at variable rates across pneumococci. Phase variants of isolate RMV7, distinguished by altered methylation patterns driven by the translocating variable restriction-modification (tvr) locus, differed significantly in their transformation efficiencies and biofilm thicknesses. These differences were replicated when the corresponding tvr alleles were introduced into an RMV7 derivative lacking the locus. RNA-seq identified differential expression of the type 1 pilus, causing the variation in biofilm formation, and inhibition of competence induction in the less transformable variant, RMV7domi. This was partly attributable to RMV7domi's lower expression of ManLMN, which promoted competence induction through importing N-acetylglucosamine. This effect was potentiated by analogues of some proteobacterial competence regulatory machinery. Additionally, one of RMV7domi's phage-related chromosomal island was relatively active, which inhibited transformation by increasing expression of the stress response proteins ClpP and HrcA. However, HrcA increased competence induction in the other variant, with its effects depending on Ca2+ supplementation and heat shock. Hence the heterogeneity in transformation efficiency likely reflects the diverse signalling pathways by which it is affected. This regulatory complexity will modulate population-wide responses to synchronising quorum sensing signals to produce co-ordinated yet stochastic bet hedging behaviour.

Bactabolize is a tool for high-throughput generation of bacterial strain-specific metabolic models

Metabolic capacity can vary substantially within a bacterial species, leading to ecological niche separation, as well as differences in virulence and antimicrobial susceptibility. Genome-scale metabolic models are useful tools for studying the metabolic potential of individuals, and with the rapid expansion of genomic sequencing there is a wealth of data that can be leveraged for comparative analysis. However, there exist few tools to construct strain-specific metabolic models at scale. Here, we describe Bactabolize, a reference-based tool which rapidly produces strain-specific metabolic models and growth phenotype predictions. We describe a pan reference model for the priority antimicrobial-resistant pathogen, Klebsiella pneumoniae, and a quality control framework for using draft genome assemblies as input for Bactabolize. The Bactabolize-derived model for K. pneumoniae reference strain KPPR1 performed comparatively or better than currently available automated approaches CarveMe and gapseq across 507 substrate and 2317 knockout mutant growth predictions. Novel draft genomes passing our systematically defined quality control criteria resulted in models with a high degree of completeness (≥99% genes and reactions captured compared to models derived from matched complete genomes) and high accuracy (mean 0.97, n=10). We anticipate the tools and framework described herein will facilitate large-scale metabolic modelling analyses that broaden our understanding of diversity within bacterial species and inform novel control strategies for priority pathogens.

The ESKAPE mobilome contributes to the spread of antimicrobial resistance and CRISPR-mediated conflict between mobile genetic elements

Mobile genetic elements (MGEs) mediate the shuffling of genes among organisms. They contribute to the spread of virulence and antibiotic resistance (AMR) genes in human pathogens, such as the particularly problematic group of ESKAPE pathogens. Here, we performed the first systematic analysis of MGEs, including plasmids, prophages, and integrative and conjugative/mobilizable elements (ICEs/IMEs), across all ESKAPE pathogens. We found that different MGE types are asymmetrically distributed across these pathogens, and that most horizontal gene transfer (HGT) events are restricted by phylum or genus. We show that the MGEs proteome is involved in diverse functional processes and distinguish widespread proteins within the ESKAPE context. Moreover, anti-CRISPRs and AMR genes are overrepresented in the ESKAPE mobilome. Our results also underscore species-specific trends shaping the number of MGEs, AMR, and virulence genes across pairs of conspecific ESKAPE genomes with and without CRISPR-Cas systems. Finally, we observed that CRISPR spacers found on prophages, ICEs/IMEs, and plasmids have different targeting biases: while plasmid and prophage CRISPRs almost exclusively target other plasmids and prophages, respectively, ICEs/IMEs CRISPRs preferentially target prophages. Overall, our study highlights the general importance of the ESKAPE mobilome in contributing to the spread of AMR and mediating conflict among MGEs.

Post-vaccine epidemiology of serotype 3 pneumococci identifies transformation inhibition through prophage-driven alteration of a non-coding RNA

BACKGROUND

The respiratory pathogen Streptococcus pneumoniae (the pneumococcus) is a genetically diverse bacterium associated with over 101 immunologically distinct polysaccharide capsules (serotypes). Polysaccharide conjugate vaccines (PCVs) have successfully eliminated multiple targeted serotypes, yet the mucoid serotype 3 has persisted despite its inclusion in PCV13. This capsule type is predominantly associated with a single globally disseminated strain, GPSC12 (clonal complex 180).

METHODS

A genomic epidemiology study combined previous surveillance datasets of serotype 3 pneumococci to analyse the population structure, dynamics, and differences in rates of diversification within GPSC12 during the period of PCV introductions. Transcriptomic analyses, whole genome sequencing, mutagenesis, and electron microscopy were used to characterise the phenotypic impact of loci hypothesised to affect this strain's evolution.

RESULTS

GPSC12 was split into clades by a genomic analysis. Clade I, the most common, rarely underwent transformation, but was typically infected with the prophage ϕOXC141. Prior to the introduction of PCV13, this clade's composition shifted towards a ϕOXC141-negative subpopulation in a systematically sampled UK collection. In the post-PCV13 era, more rapidly recombining non-Clade I isolates, also ϕOXC141-negative, have risen in prevalence. The low in vitro transformation efficiency of a Clade I isolate could not be fully explained by the ~100-fold reduction attributable to the serotype 3 capsule. Accordingly, prophage ϕOXC141 was found to modify csRNA3, a non-coding RNA that inhibits the induction of transformation. This alteration was identified in ~30% of all pneumococci and was particularly common in the unusually clonal serotype 1 GPSC2 strain. RNA-seq and quantitative reverse transcriptase PCR experiments using a genetically tractable pneumococcus demonstrated the altered csRNA3 was more effective at inhibiting production of the competence-stimulating peptide pheromone. This resulted in a reduction in the induction of competence for transformation.

CONCLUSION

This interference with the quorum sensing needed to induce competence reduces the risk of the prophage being deleted by homologous recombination. Hence the selfish prophage-driven alteration of a regulatory RNA limits cell-cell communication and horizontal gene transfer, complicating the interpretation of post-vaccine population dynamics.

Epidemiological dynamics of bacteriocin competition and antibiotic resistance

Bacteriocins, toxic peptides involved in the competition between bacterial strains, are extremely diverse. Previous work on bacteriocin dynamics has highlighted the role of non-transitive 'rock-paper-scissors' competition in maintaining the coexistence of different bacteriocin profiles. The focus to date has primarily been on bacteriocin interactions at the within-host scale (i.e. within a single bacterial population). Yet in species such as Streptococcus pneumoniae, with relatively short periods of colonization and limited within-host diversity, ecological outcomes are also shaped by processes at the epidemiological (between-host) scale. Here, we first investigate bacteriocin dynamics and diversity in epidemiological models. We find that in these models, bacteriocin diversity is more readily maintained than in within-host models, and with more possible combinations of coexisting bacteriocin profiles. Indeed, maintenance of diversity in epidemiological models does not require rock-paper-scissors dynamics; it can also occur through a competition-colonization trade-off. Second, we investigate the link between bacteriocin diversity and diversity at antibiotic resistance loci. Previous work has proposed that bacterial duration of colonization modulates the fitness of antibiotic resistance. Due to their inhibitory effects, bacteriocins are a plausible candidate for playing a role in the duration of colonization episodes. We extend the epidemiological model of bacteriocin dynamics to incorporate an antibiotic resistance locus and demonstrate that bacteriocin diversity can indeed maintain the coexistence of antibiotic-sensitive and -resistant strains.

Detecting bacterial adaptation within individual microbiomes

The human microbiome harbours a large capacity for within-person adaptive mutations. Commensal bacterial strains can stably colonize a person for decades, and billions of mutations are generated daily within each person's microbiome. Adaptive mutations emerging during health might be driven by selective forces that vary across individuals, vary within an individual, or are completely novel to the human population. Mutations emerging within individual microbiomes might impact the immune system, the metabolism of nutrients or drugs, and the stability of the community to perturbations. Despite this potential, relatively little attention has been paid to the possibility of adaptive evolution within complex human-associated microbiomes. This review discusses the promise of studying within-microbiome adaptation, the conceptual and technical limitations that may have contributed to an underappreciation of adaptive de novo mutations occurring within microbiomes to date, and methods for detecting recent adaptive evolution. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

Metagenomic methylation patterns resolve bacterial genomes of unusual size and structural complexity

The plasticity of bacterial and archaeal genomes makes examining their ecological and evolutionary dynamics both exciting and challenging. The same mechanisms that enable rapid genomic change and adaptation confound current approaches for recovering complete genomes from metagenomes. Here, we use strain-specific patterns of DNA methylation to resolve complex bacterial genomes from long-read metagenomic data of a marine microbial consortium, the “pink berries” of the Sippewissett Marsh (USA). Unique combinations of restriction-modification (RM) systems encoded by the bacteria produced distinctive methylation profiles that were used to accurately bin and classify metagenomic sequences. Using this approach, we finished the largest and most complex circularized bacterial genome ever recovered from a metagenome (7.9 Mb with >600 transposons), the finished genome of Thiohalocapsa sp. PB-PSB1 the dominant bacteria in the consortia. From genomes binned by methylation patterns, we identified instances of horizontal gene transfer between sulfur-cycling symbionts (Thiohalocapsa sp. PB-PSB1 and Desulfofustis sp. PB-SRB1), phage infection, and strain-level structural variation. We also linked the methylation patterns of each metagenome-assembled genome with encoded DNA methyltransferases and discovered new RM defense systems, including novel associations of RM systems with RNase toxins.

Comparative genome analysis of four Leuconostoc strains with a focus on carbohydrate-active enzymes and oligosaccharide utilization pathways

•

Comparative genomic analysis of four Leuconostoc strains was performed.

•

Leuconostoc spp. shared genomic similarity, but their genetic content differed.

•

Leuconostoc spp. showed different genes encoding CAZymes.

•

Oligosaccharide’s utilization and folate biosynthesis pathways were investigated.

Leuconostoc is mostly found in food, plants, and dairy products. Due to their innate genomic features, such as the presence of carbohydrate-active enzymes, bacteriocins, and plasmids, Leuconostoc spp. have great biotechnological potential. In this study, four strains were isolated and identified as Leuconostoc mesenteroides SG315 (LA), L. citreum SG255 (LB), L. lactis CCK940 (LC), and L. lactis SBC001 (LD). Comparative analysis was performed using their draft genome sequences. Differences among the four strains were analyzed using the average nucleotide identity, dot plot, and multiple alignments of conserved genomic sequences. Functional profiling revealed 2134, 1917, 1751, and 1816 open reading frames; 2023, 1823, 1655, and 1699 protein-coding genes; 60, 57, 83, and 82 RNA-coding genes; and GC content of 37.5 %, 38.8 %, 43.3 %, and 43.2 %, in LA, LB, LC, and LD, respectively. The total number of genes encoding carbohydrate-active enzymes was 76 (LA), 73 (LB), 57 (LC), and 67 (LD). These results indicate that the four strains shared a large number of genes, but their gene content is different. Furthermore, most genes with unknown functions were observed in the prophage regions of the genome. This study also elucidated the oligosaccharide utilization and folate biosynthesis pathways in Leuconostoc spp. Taken together, our findings provide useful information on the genomic diversity of CAZymes in the four Leuconostoc strains and suggest that these species could be used for potent exploitation.

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 to core genes, accessory genes are enriched in genome plasticity functions. We found no evidence for recent horizontal gene transfer 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.​.

Analysing pneumococcal invasiveness using Bayesian models of pathogen progression rates

The disease burden attributable to opportunistic pathogens depends on their prevalence in asymptomatic colonisation and the rate at which they progress to cause symptomatic disease. Increases in infections caused by commensals can result from the emergence of “hyperinvasive” strains. Such pathogens can be identified through quantifying progression rates using matched samples of typed microbes from disease cases and healthy carriers. This study describes Bayesian models for analysing such datasets, implemented in an RStan package (https://github.com/nickjcroucher/progressionEstimation). The models converged on stable fits that accurately reproduced observations from meta-analyses of Streptococcus pneumoniae datasets. The estimates of invasiveness, the progression rate from carriage to invasive disease, in cases per carrier per year correlated strongly with the dimensionless values from meta-analysis of odds ratios when sample sizes were large. At smaller sample sizes, the Bayesian models produced more informative estimates. This identified historically rare but high-risk S. pneumoniae serotypes that could be problematic following vaccine-associated disruption of the bacterial population. The package allows for hypothesis testing through model comparisons with Bayes factors. Application to datasets in which strain and serotype information were available for S. pneumoniae found significant evidence for within-strain and within-serotype variation in invasiveness. The heterogeneous geographical distribution of these genotypes is therefore likely to contribute to differences in the impact of vaccination in between locations. Hence genomic surveillance of opportunistic pathogens is crucial for quantifying the effectiveness of public health interventions, and enabling ongoing meta-analyses that can identify new, highly invasive variants.

Characterization of the Type I Restriction Modification System Broadly Conserved among Group A Streptococci

Although prokaryotic DNA methylation investigations have long focused on immunity against exogenous DNA, it has been recently recognized that DNA methylation impacts gene expression and phase variation in Streptococcus pneumoniae and Streptococcus suis. A comprehensive analysis of DNA methylation is lacking for beta-hemolytic streptococci, and thus we sought to examine DNA methylation in the major human pathogen group A Streptococcus (GAS). Using a database of 224 GAS genomes encompassing 80 emm types, we found that nearly all GAS strains encode a type I restriction modification (RM) system that lacks the hsdS′ alleles responsible for impacting gene expression in S. pneumoniae and S. suis. The GAS type I system is located on the core chromosome, while sporadically present type II orphan methyltransferases were identified on prophages. By combining single-molecule real-time (SMRT) analyses of 10 distinct emm types along with phylogenomics of 224 strains, we were able to assign 13 methylation patterns to the GAS population. Inactivation of the type I RM system, occurring either naturally through phage insertion or through laboratory-induced gene deletion, abrogated DNA methylation detectable via either SMRT or MinION sequencing. Contrary to a previous report, inactivation of the type I system did not impact transcript levels of the gene (mga) encoding the key multigene activator protein (Mga) or Mga-regulated genes. Inactivation of the type I system significantly increased plasmid transformation rates. These data delineate the breadth of the core chromosomal type I RM system in the GAS population and clarify its role in immunity rather than impacting Mga regulon expression.

IMPORTANCE The advent of whole-genome approaches capable of detecting DNA methylation has markedly expanded appreciation of the diverse roles of epigenetic modification in prokaryotic physiology. For example, recent studies have suggested that DNA methylation impacts gene expression in some streptococci. The data described herein are from the first systematic analysis of DNA methylation in a beta-hemolytic streptococcus and one of the few analyses to comprehensively characterize DNA methylation across hundreds of strains of the same bacterial species. We clarify that DNA methylation in group A Streptococcus (GAS) is primarily due to a type I restriction modification (RM) system present in the core genome and does not impact mga-regulated virulence gene expression, but does impact immunity against exogenous DNA. The identification of the DNA motifs recognized by each type I RM system may assist with optimizing methods for GAS genetic manipulation and help us understand how bacterial pathogens acquire exogenous DNA elements.

Evolutionary Pathways and Trajectories in Antibiotic Resistance

Evolution is the hallmark of life. Descriptions of the evolution of microorganisms have provided a wealth of information, but knowledge regarding "what happened" has precluded a deeper understanding of "how" evolution has proceeded, as in the case of antimicrobial resistance. The difficulty in answering the "how" question lies in the multihierarchical dimensions of evolutionary processes, nested in complex networks, encompassing all units of selection, from genes to communities and ecosystems. At the simplest ontological level (as resistance genes), evolution proceeds by random (mutation and drift) and directional (natural selection) processes; however, sequential pathways of adaptive variation can occasionally be observed, and under fixed circumstances (particular fitness landscapes), evolution is predictable. At the highest level (such as that of plasmids, clones, species, microbiotas), the systems' degrees of freedom increase dramatically, related to the variable dispersal, fragmentation, relatedness, or coalescence of bacterial populations, depending on heterogeneous and changing niches and selective gradients in complex environments. Evolutionary trajectories of antibiotic resistance find their way in these changing landscapes subjected to random variations, becoming highly entropic and therefore unpredictable. However, experimental, phylogenetic, and ecogenetic analyses reveal preferential frequented paths (highways) where antibiotic resistance flows and propagates, allowing some understanding of evolutionary dynamics, modeling and designing interventions. Studies on antibiotic resistance have an applied aspect in improving individual health, One Health, and Global Health, as well as an academic value for understanding evolution. Most importantly, they have a heuristic significance as a model to reduce the negative influence of anthropogenic effects on the environment.

Heuristic and Hierarchical-Based Population Mining of Salmonella enterica Lineage I Pan-Genomes as a Platform to Enhance Food Safety

The recent incorporation of bacterial whole-genome sequencing (WGS) into Public Health laboratories has enhanced foodborne outbreak detection and source attribution. As a result, large volumes of publicly available datasets can be used to study the biology of foodborne pathogen populations at an unprecedented scale. To demonstrate the application of a heuristic and agnostic hierarchical population structure guided pan-genome enrichment analysis (PANGEA), we used populations of S. enterica lineage I to achieve two main objectives: (i) show how hierarchical population inquiry at different scales of resolution can enhance ecological and epidemiological inquiries; and (ii) identify population-specific inferable traits that could provide selective advantages in food production environments. Publicly available WGS data were obtained from NCBI database for three serovars of Salmonella enterica subsp. enterica lineage I (S. Typhimurium, S. Newport, and S. Infantis). Using the hierarchical genotypic classifications (Serovar, BAPS1, ST, cgMLST), datasets from each of the three serovars showed varying degrees of clonal structuring. When the accessory genome (PANGEA) was mapped onto these hierarchical structures, accessory loci could be linked with specific genotypes. A large heavy-metal resistance mobile element was found in the Monophasic ST34 lineage of S. Typhimurium, and laboratory testing showed that Monophasic isolates have on average a higher degree of copper resistance than the Biphasic ones. In S. Newport, an extra sugE gene copy was found among most isolates of the ST45 lineage, and laboratory testing of multiple isolates confirmed that isolates of S. Newport ST45 were on average less sensitive to the disinfectant cetylpyridimium chloride than non-ST45 isolates. Lastly, data-mining of the accessory genomic content of S. Infantis revealed two cryptic Ecotypes with distinct accessory genomic content and distinct ecological patterns. Poultry appears to be the major reservoir for Ecotype 1, and temporal analysis further suggested a recent ecological succession, with Ecotype 2 apparently being displaced by Ecotype 1. Altogether, the use of a heuristic hierarchical-based population structure analysis that includes bacterial pan-genomes (core and accessory genomes) can (1) improve genomic resolution for mapping populations and accessing epidemiological patterns; and (2) define lineage-specific informative loci that may be associated with survival in the food chain.

Evolutionary Ecology and Interplay of Prokaryotic Innate and Adaptive Immune Systems

Like many organisms, bacteria and archaea have both innate and adaptive immune systems to defend against infection by viruses and other parasites. Innate immunity most commonly relies on the endonuclease-mediated cleavage of any incoming DNA that lacks a specific epigenetic modification, through a system known as restriction-modification. CRISPR-Cas-mediated adaptive immunity relies on the insertion of short DNA sequences from parasite genomes into CRISPR arrays on the host genome to provide sequence-specific protection. The discovery of each of these systems has revolutionised our ability to carry out genetic manipulations, and, as a consequence, the enzymes involved have been characterised in exquisite detail. In comparison, much less is known about the importance of these two arms of the defence for the ecology and evolution of prokaryotes and their parasites. Here, we review our current ecological and evolutionary understanding of these systems in isolation, and discuss the need to study how innate and adaptive immune responses are integrated when they coexist in the same cell.

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

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

Distribution of Antimicrobial Resistance and Virulence Genes within the Prophage-Associated Regions in Nosocomial Pathogens

Prophages are often involved in host survival strategies and contribute toward increasing the genetic diversity of the host genome. Prophages also drive horizontal propagation of various genes as vehicles. However, there are few retrospective studies contributing to the propagation of antimicrobial resistance (AMR) and virulence factor (VF) genes by prophage. We extracted the complete genome sequences of seven pathogens, including ESKAPE bacteria and Escherichia coli from a public database, and examined the distribution of both the AMR and VF genes in prophage-like regions. We found that the ratios of AMR and VF genes greatly varied among the seven species. More than 70% of Enterobacter cloacae strains had VF genes, but only 1.2% of Klebsiella pneumoniae strains had VF genes from prophages. AMR and VF genes are unlikely to exist together in the same prophage region except in E. coli and Staphylococcus aureus, and the distribution patterns of prophage types containing AMR genes are distinct from those of VF gene-carrying prophage types. AMR genes in the prophage were located near transposase and/or integrase. The prophage containing class 1 integrase possessed a significantly greater number of AMR genes than did prophages with no class 1 integrase. The results of this study present a comprehensive picture of AMR and VF genes present within, or close to, prophage-like elements and different prophage patterns between AMR- or VF-encoding prophage-like elements.

IMPORTANCE Although we believe phages play an important role in horizontal gene transfer in exchanging genetic material, we do not know the distribution of the antimicrobial resistance (AMR) and/or virulence factor (VF) genes in prophages. We collected different prophage elements from the complete genome sequences of seven species—Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, Enterobacter cloacae, and Escherichia coli—and characterized the distribution of antimicrobial resistance and virulence genes located in the prophage region. While virulence genes in prophage were species specific, antimicrobial resistance genes in prophages were highly conserved in various species. An integron structure was detected within specific prophage regions such as P1-like prophage element. Maximum of 10 antimicrobial resistance genes were found in a single prophage region, suggesting that prophages act as a reservoir for antimicrobial resistance genes. The results of this study show the different characteristic structures between AMR- or VF-encoding prophages.

ProkEvo: an automated, reproducible, and scalable framework for high-throughput bacterial population genomics analyses

Whole Genome Sequence (WGS) data from bacterial species is used for a variety of applications ranging from basic microbiological research, diagnostics, and epidemiological surveillance. The availability of WGS data from hundreds of thousands of individual isolates of individual microbial species poses a tremendous opportunity for discovery and hypothesis-generating research into ecology and evolution of these microorganisms. Flexibility, scalability, and user-friendliness of existing pipelines for population-scale inquiry, however, limit applications of systematic, population-scale approaches. Here, we present ProkEvo, an automated, scalable, reproducible, and open-source framework for bacterial population genomics analyses using WGS data. ProkEvo was specifically developed to achieve the following goals: (1) Automation and scaling of complex combinations of computational analyses for many thousands of bacterial genomes from inputs of raw Illumina paired-end sequence reads; (2) Use of workflow management systems (WMS) such as Pegasus WMS to ensure reproducibility, scalability, modularity, fault-tolerance, and robust file management throughout the process; (3) Use of high-performance and high-throughput computational platforms; (4) Generation of hierarchical-based population structure analysis based on combinations of multi-locus and Bayesian statistical approaches for classification for ecological and epidemiological inquiries; (5) Association of antimicrobial resistance (AMR) genes, putative virulence factors, and plasmids from curated databases with the hierarchically-related genotypic classifications; and (6) Production of pan-genome annotations and data compilation that can be utilized for downstream analysis such as identification of population-specific genomic signatures. The scalability of ProkEvo was measured with two datasets comprising significantly different numbers of input genomes (one with ~2,400 genomes, and the second with ~23,000 genomes). Depending on the dataset and the computational platform used, the running time of ProkEvo varied from ~3-26 days. ProkEvo can be used with virtually any bacterial species, and the Pegasus WMS uniquely facilitates addition or removal of programs from the workflow or modification of options within them. To demonstrate versatility of the ProkEvo platform, we performed a hierarchical-based population structure analyses from available genomes of three distinct pathogenic bacterial species as individual case studies. The specific case studies illustrate how hierarchical analyses of population structures, genotype frequencies, and distribution of specific gene functions can be integrated into an analysis. Collectively, our study shows that ProkEvo presents a practical viable option for scalable, automated analyses of bacterial populations with direct applications for basic microbiology research, clinical microbiological diagnostics, and epidemiological surveillance.

Global genomic pathogen surveillance to inform vaccine strategies: a decade-long expedition in pneumococcal genomics

Vaccines are powerful agents in infectious disease prevention but often designed to protect against some strains that are most likely to spread and cause diseases. Most vaccines do not succeed in eradicating the pathogen and thus allow the potential emergence of vaccine evading strains. As with most evolutionary processes, being able to capture all variations across the entire genome gives us the best chance of monitoring and understanding the processes of vaccine evasion. Genomics is being widely adopted as the optimum approach for pathogen surveillance with the potential for early and precise identification of high-risk strains. Given sufficient longitudinal data, genomics also has the potential to forecast the emergence of such strains enabling immediate or pre-emptive intervention. In this review, we consider the strengths and challenges for pathogen genomic surveillance using the experience of the Global Pneumococcal Sequencing (GPS) project as an early example. We highlight the multifaceted nature of genome data and recent advances in genome-based tools to extract useful information relevant to inform vaccine strategies and treatment options. We conclude with future perspectives for genomic pathogen surveillance.

Global population structure and genotyping framework for genomic surveillance of the major dysentery pathogen, Shigella sonnei

Shigella sonnei is the most common agent of shigellosis in high-income countries, and causes a significant disease burden in low- and middle-income countries. Antimicrobial resistance is increasingly common in all settings. Whole genome sequencing (WGS) is increasingly utilised for S. sonnei outbreak investigation and surveillance, but comparison of data between studies and labs is challenging. Here, we present a genomic framework and genotyping scheme for S. sonnei to efficiently identify genotype and resistance determinants from WGS data. The scheme is implemented in the software package Mykrobe and tested on thousands of genomes. Applying this approach to analyse >4,000 S. sonnei isolates sequenced in public health labs in three countries identified several common genotypes associated with increased rates of ciprofloxacin resistance and azithromycin resistance, confirming intercontinental spread of highly-resistant S. sonnei clones and demonstrating the genomic framework can facilitate monitoring the spread of resistant clones, including those that have recently emerged, at local and global scales.

Negative frequency-dependent selection and asymmetrical transformation stabilise multi-strain bacterial population structures

Streptococcus pneumoniae can be divided into many strains, each a distinct set of isolates sharing similar core and accessory genomes, which co-circulate within the same hosts. Previous analyses suggested the short-term vaccine-associated dynamics of S. pneumoniae strains may be mediated through multi-locus negative frequency-dependent selection (NFDS), which maintains accessory loci at equilibrium frequencies. Long-term simulations demonstrated NFDS stabilised clonally-evolving multi-strain populations through preventing the loss of variation through drift, based on polymorphism frequencies, pairwise genetic distances and phylogenies. However, allowing symmetrical recombination between isolates evolving under multi-locus NFDS generated unstructured populations of diverse genotypes. Replication of the observed data improved when multi-locus NFDS was combined with recombination that was instead asymmetrical, favouring deletion of accessory loci over insertion. This combination separated populations into strains through outbreeding depression, resulting from recombinants with reduced accessory genomes having lower fitness than their parental genotypes. Although simplistic modelling of recombination likely limited these simulations' ability to maintain some properties of genomic data as accurately as those lacking recombination, the combination of asymmetrical recombination and multi-locus NFDS could restore multi-strain population structures from randomised initial populations. As many bacteria inhibit insertions into their chromosomes, this combination may commonly underlie the co-existence of strains within a niche.

Genome Informatics and Machine Learning-Based Identification of Antimicrobial Resistance-Encoding Features and Virulence Attributes in Escherichia coli Genomes Representing Globally Prevalent Lineages, Including High-Risk Clonal Complexes

Escherichia coli, a ubiquitous commensal/pathogenic member from the Enterobacteriaceae family, accounts for high infection burden, morbidity, and mortality throughout the world. With emerging multidrug resistance (MDR) on a massive scale, E. coli has been listed as one of the Global Antimicrobial Resistance and Use Surveillance System (GLASS) priority pathogens. Understanding the resistance mechanisms and underlying genomic features appears to be of utmost importance to tackle further spread of these multidrug-resistant superbugs. While a few of the globally prevalent sequence types (STs) of E. coli, such as ST131, ST69, ST405, and ST648, have been previously reported to be highly virulent and harboring MDR, there is no clarity if certain ST lineages have a greater propensity to acquire MDR. In this study, large-scale comparative genomics of a total of 5,653 E. coli genomes from 19 ST lineages revealed ST-wide prevalence patterns of genomic features, such as antimicrobial resistance (AMR)-encoding genes/mutations, virulence genes, integrons, and transposons. Interpretation of the importance of these features using a Random Forest Classifier trained with 11,988 genomic features from whole-genome sequence data identified ST-specific or phylogroup-specific signature proteins mostly belonging to different protein superfamilies, including the toxin-antitoxin systems. Our study provides a comprehensive understanding of a myriad of genomic features, ST-specific proteins, and resistance mechanisms entailing different lineages of E. coli at the level of genomes; this could be of significant downstream importance in understanding the mechanisms of AMR, in clinical discovery, in epidemiology, and in devising control strategies. IMPORTANCE With the leap in whole-genome data being generated, the application of relevant methods to mine biologically significant information from microbial genomes is of utmost importance to public health genomics. Machine-learning methods have been used not only to mine, curate, or classify the data but also to identify the relevant features that could be linked to a particular class/target. This is perhaps one of the pioneering studies that has attempted to classify a large repertoire of E. coli genome data sets (5,653 genomes) belonging to 19 different STs (including well-studied as well as understudied STs) using machine learning approaches. Important features identified by these approaches have revealed ST-specific signature proteins, which could be further studied to predict possible associations with the phenotypic profiles, thereby providing a better understanding of virulence and the resistance mechanisms among different clonal lineages of E. coli.

Time-calibrated genomic evolution of a monomorphic bacterium during its establishment as an endemic crop pathogen

Horizontal gene transfer is of major evolutionary importance as it allows for the redistribution of phenotypically important genes among lineages. Such genes with essential functions include those involved in resistance to antimicrobial compounds and virulence factors in pathogenic bacteria. Understanding gene turnover at microevolutionary scales is critical to assess the pace of this evolutionary process. Here, we characterized and quantified gene turnover for the epidemic lineage of a bacterial plant pathogen of major agricultural importance worldwide. Relying on a dense geographic sampling spanning 39 years of evolution, we estimated both the dynamics of single nucleotide polymorphism accumulation and gene content turnover. We identified extensive gene content variation among lineages even at the smallest phylogenetic and geographic scales. Gene turnover rate exceeded nucleotide substitution rate by three orders of magnitude. Accessory genes were found preferentially located on plasmids, but we identified a highly plastic chromosomal region hosting ecologically important genes such as transcription activator-like effectors. Whereas most changes in the gene content are probably transient, the rapid spread of a mobile element conferring resistance to copper compounds widely used for the management of plant bacterial pathogens illustrates how some accessory genes can become ubiquitous within a population over short timeframes.

Prevalence of phase variable epigenetic invertons among host-associated bacteria

Type I restriction-modification (R-M) systems consist of a DNA endonuclease (HsdR, HsdM and HsdS subunits) and methyltransferase (HsdM and HsdS subunits). The hsdS sequences flanked by inverted repeats (referred to as epigenetic invertons) in certain Type I R-M systems undergo invertase-catalyzed inversions. Previous studies in Streptococcus pneumoniae have shown that hsdS inversions within clonal populations produce subpopulations with profound differences in the methylome, cellular physiology and virulence. In this study, we bioinformatically identified six major clades of the tyrosine and serine family invertases homologs from 16 bacterial phyla, which potentially catalyze hsdS inversions in the epigenetic invertons. In particular, the epigenetic invertons are highly enriched in host-associated bacteria. We further verified hsdS inversions in the Type I R-M systems of four representative host-associated bacteria and found that each of the resultant hsdS allelic variants specifies methylation of a unique DNA sequence. In addition, transcriptome analysis revealed that hsdS allelic variations in Enterococcus faecalis exert significant impact on gene expression. These findings indicate that epigenetic switches driven by invertases in the epigenetic invertons broadly operate in the host-associated bacteria, which may broadly contribute to bacterial host adaptation and virulence beyond the role of the Type I R-M systems against phage infection.

Accurate reconstruction of bacterial pan- and core genomes with PEPPAN

Bacterial genomes can contain traces of a complex evolutionary history, including extensive homologous recombination, gene loss, gene duplications, and horizontal gene transfer. To reconstruct the phylogenetic and population history of a set of multiple bacteria, it is necessary to examine their pangenome, the composite of all the genes in the set. Here we introduce PEPPAN, a novel pipeline that can reliably construct pangenomes from thousands of genetically diverse bacterial genomes that represent the diversity of an entire genus. PEPPAN outperforms existing pangenome methods by providing consistent gene and pseudogene annotations extended by similarity-based gene predictions, and identifying and excluding paralogs by combining tree- and synteny-based approaches. The PEPPAN package additionally includes PEPPAN_parser, which implements additional downstream analyses, including the calculation of trees based on accessory gene content or allelic differences between core genes. To test the accuracy of PEPPAN, we implemented SimPan, a novel pipeline for simulating the evolution of bacterial pangenomes. We compared the accuracy and speed of PEPPAN with four state-of-the-art pangenome pipelines using both empirical and simulated data sets. PEPPAN was more accurate and more specific than any of the other pipelines and was almost as fast as any of them. As a case study, we used PEPPAN to construct a pangenome of approximately 40,000 genes from 3052 representative genomes spanning at least 80 species of Streptococcus The resulting gene and allelic trees provide an unprecedented overview of the genomic diversity of the entire Streptococcus genus.

Lysogeny in Streptococcus pneumoniae

Bacterial viruses, or bacteriophages, are major contributors to the evolution, pathogenesis and overall biology of their host bacteria. During their life cycle, temperate bacteriophages form stable associations with their host by integrating into the chromosome, a process called lysogeny. Isolates of the human pathogen Streptococcus pneumoniae are frequently lysogenic, and genomic studies have allowed the classification of these phages into distinct phylogenetic groups. Here, we review the recent advances in the characterization of temperate pneumococcal phages, with a focus on their genetic features and chromosomal integration loci. We also discuss the contribution of phages, and specific phage-encoded features, to colonization and virulence. Finally, we discuss interesting research perspectives in this field.

Recent advances in understanding the ecology of the lung microbiota and deciphering the gut-lung axis

A rapidly expanding new field of lung research has been produced by the emergence of culture-independent next-generation sequencing technologies. While pulmonary microbiome research lags behind the exploration of the microbiome in other organ systems, the field is maturing and has recently produced multiple exciting discoveries. In this mini-review, we will explore recent advances in our understanding of the lung microbiome and the gut-lung axis from an ecological perspective.

Frequency-dependent selection can forecast evolution in Streptococcus pneumoniae

Predicting how pathogen populations will change over time is challenging. Such has been the case with Streptococcus pneumoniae, an important human pathogen, and the pneumococcal conjugate vaccines (PCVs), which target only a fraction of the strains in the population. Here, we use the frequencies of accessory genes to predict changes in the pneumococcal population after vaccination, hypothesizing that these frequencies reflect negative frequency-dependent selection (NFDS) on the gene products. We find that the standardized predicted fitness of a strain, estimated by an NFDS-based model at the time the vaccine is introduced, enables us to predict whether the strain increases or decreases in prevalence following vaccination. Further, we are able to forecast the equilibrium post-vaccine population composition and assess the invasion capacity of emerging lineages. Overall, we provide a method for predicting the impact of an intervention on pneumococcal populations with potential application to other bacterial pathogens in which NFDS is a driving force.

Site-Specific Recombination - How Simple DNA Inversions Produce Complex Phenotypic Heterogeneity in Bacterial Populations

Many bacterial species generate phenotypically heterogeneous subpopulations as a strategy for ensuring the survival of the population as a whole - an environmental stress that eradicates one subpopulation may leave other phenotypic groups unharmed, allowing the lineage to continue. Phase variation, a process that functions as an ON/OFF switch for gene expression, is one way that bacteria achieve phenotypic heterogeneity. Phase variation occurs stochastically and reversibly, and in the presence of a selective pressure the advantageous phenotype(s) predominates in the population. Phase variation can occur through multiple genetic and epigenetic mechanisms. This review focuses on conservative site-specific recombination that generates reversible DNA inversions as a genetic mechanism mediating phase variation. Recent studies have sparked a renewed interest in phase variation mediated through DNA inversion, revealing a high level of complexity beyond simple ON/OFF switching, including unusual modes of gene regulation, and highlighting an underappreciation of the use of these mechanisms by bacteria.

Recent advances in understanding the ecology of the lung microbiota and deciphering the gut-lung axis

A rapidly expanding new field of lung research has been produced by the emergence of culture-independent next-generation sequencing technologies. While pulmonary microbiome research lags behind the exploration of the microbiome in other organ systems, the field is maturing and has recently produced multiple exciting discoveries. In this mini-review, we will explore recent advances in our understanding of the lung microbiome and the gut-lung axis from an ecological perspective.

Potential evolutionary impact of integrative and conjugative elements (ICEs) and genomic islands in the Ralstonia solanacearum species complex

Ralstonia solanacearum, a soil-borne plant pathogen, encompasses a large number of strains known as R. solanacearum species complex (RSSC). Although it has been suggested that mobile genetic elements (MGEs) may play an important role in the RSSC genome, the evolutionary impact of these elements remains unknown. Here, we identified and analysed Integrative and Conjugative Elements (ICEs) and Genomic Islands (GIs) in the 121 genomes published for Ralstonia spp., including RSSC strains and three other non-plant pathogenic Ralstonia spp. Our results provided a dataset of 12 ICEs and 31 GIs distributed throughout Ralstonia spp. Four novel ICEs in RSSC were found. Some of these elements cover 5% of the host genome and carry accessory genes with a potential impact on the fitness and pathogenicity of RSSC. In addition, phylogenetic analysis revealed that these MGEs clustered to the same species, but there is evidence of strains from different countries that host the same element. Our results provide novel insight into the RSSC adaptation, opening new paths to a better understanding of how these elements affect this soil-borne plant pathogen.

Producing polished prokaryotic pangenomes with the Panaroo pipeline

Population-level comparisons of prokaryotic genomes must take into account the substantial differences in gene content resulting from horizontal gene transfer, gene duplication and gene loss. However, the automated annotation of prokaryotic genomes is imperfect, and errors due to fragmented assemblies, contamination, diverse gene families and mis-assemblies accumulate over the population, leading to profound consequences when analysing the set of all genes found in a species. Here, we introduce Panaroo, a graph-based pangenome clustering tool that is able to account for many of the sources of error introduced during the annotation of prokaryotic genome assemblies. Panaroo is available at https://github.com/gtonkinhill/panaroo .

Dynamic changes in paediatric invasive pneumococcal disease after sequential switches of conjugate vaccine in Belgium: a national retrospective observational study

BACKGROUND

Ten-valent and 13-valent pneumococcal conjugate vaccines (PCVs) have shown important benefits by decreasing invasive pneumococcal disease caused by vaccine serotypes. Belgium had an uncommon situation with sequential use of PCV7, PCV13, and PCV10 in the childhood vaccination programmes between 2007 and 2018. We aimed to analyse the changes in incidence of invasive pneumococcal disease and serotype distribution in children throughout this period.

METHODS

Streptococcus pneumoniae isolates were obtained from patients with invasive pneumococcal disease in Belgium between 2007 and 2018 by the national laboratory-based surveillance. Paediatric invasive pneumococcal disease incidence, serotype distribution, and antimicrobial susceptibility were analysed in periods during which PCV7 (2009-10), PCV13 (2013-14), both PCV13 and PCV10 (2015-16), and PCV10 (2017-18) were used. Incidence rates and trends were compared. Vaccination status was collected. For a subset of serotype 19A isolates, multilocus sequence type was identified.

FINDINGS

After a decrease in PCV7 serotype invasive pneumococcal disease was observed during the PCV7 period, total paediatric invasive pneumococcal disease incidence significantly declined during the PCV13 period (-2·6% monthly, p<0·0001). During the PCV13-PCV10 period (2015-16), the lowest mean in paediatric invasive pneumococcal disease incidence was achieved, but the incidence increased again during the PCV10 period (2017-18), especially in children younger than 2 years (+1·7% monthly; p=0·028). This increase was mainly due to a significant rise in serotype 19A invasive pneumococcal disease incidence in the PCV10 period compared with the PCV13 period (p<0·0001), making serotype 19A the predominant serotype in paediatric invasive pneumococcal disease in the PCV10 period. Genetic diversity within the 2017-18 serotype 19A collection was seen, with two predominant clones, ST416 and ST994, that were infrequently observed before PCV10 introduction. In 2018, among children younger than 5 years with invasive pneumococcal disease who were correctly vaccinated, 37% (37 of 100) had PCV13 serotype invasive pneumococcal disease, all caused by serotype 19A and serotype 3.

INTERPRETATION

After a significant decrease during the PCV13 period, paediatric invasive pneumococcal disease incidence increased again during the PCV10 period. This observation mainly resulted from a significant increase of serotype 19A cases. During the PCV10 period, dominant serotype 19A clones differed from those detected during previous vaccine periods. Whether changes in epidemiology resulted from the vaccine switch or also from natural evolution remains to be further elucidated.

FUNDING

The Belgian National Reference is funded by the Belgian National Institute for Health and Disability Insurance and the whole genome sequencing by an investigator-initiated research grant from Pfizer.

Population genomics of Staphylococcus pseudintermedius in companion animals in the United States

Staphylococcus pseudintermedius is a commensal bacterium and a major opportunistic pathogen of dogs. The emergence of methicillin-resistant S. pseudintermedius (MRSP) is also becoming a serious concern. We carried out a population genomics study of 130 clinical S. pseudintermedius isolates from dogs and cats in the New England region of the United States. Results revealed the co-circulation of phylogenetically diverse lineages that have access to a large pool of accessory genes. Many MRSP and multidrug-resistant clones have emerged through multiple independent, horizontal acquisition of resistance determinants and frequent genetic exchange that disseminate DNA to the broader population. When compared to a Texas population, we found evidence of clonal expansion of MRSP lineages that have disseminated over large distances. These findings provide unprecedented insight into the diversification of a common cutaneous colonizer of man's oldest companion animal and the widespread circulation of multiple high-risk resistant clones.

Conserved DNA Methyltransferases: A Window into Fundamental Mechanisms of Epigenetic Regulation in Bacteria

An increasing number of studies have reported that bacterial DNA methylation has important functions beyond the roles in restriction-modification systems, including the ability of affecting clinically relevant phenotypes such as virulence, host colonization, sporulation, biofilm formation, among others. Although insightful, such studies have a largely ad hoc nature and would benefit from a systematic strategy enabling a joint functional characterization of bacterial methylomes by the microbiology community. In this opinion article, we propose that highly conserved DNA methyltransferases (MTases) represent a unique opportunity for bacterial epigenomic studies. These MTases are rather common in bacteria, span various taxonomic scales, and are present in multiple human pathogens. Apart from well-characterized core DNA MTases, like those from Vibrio cholerae, Salmonella enterica, Clostridioides difficile, or Streptococcus pyogenes, multiple highly conserved DNA MTases are also found in numerous human pathogens, including those belonging to the genera Burkholderia and Acinetobacter. We discuss why and how these MTases can be prioritized to enable a community-wide, integrative approach for functional epigenomic studies. Ultimately, we discuss how some highly conserved DNA MTases may emerge as promising targets for the development of novel epigenetic inhibitors for biomedical applications.

Designing ecologically optimized pneumococcal vaccines using population genomics

Streptococcus pneumoniae (the pneumococcus) is a common nasopharyngeal commensal that can cause invasive pneumococcal disease (IPD). Each component of current protein-polysaccharide conjugate vaccines (PCVs) generally induces immunity specific to one of the approximately 100 pneumococcal serotypes, and typically eliminates it from carriage and IPD through herd immunity. Overall carriage rates remain stable owing to replacement by non-PCV serotypes. Consequently, the net change in IPD incidence is determined by the relative invasiveness of the pre- and post-PCV-carried pneumococcal populations. In the present study, we identified PCVs expected to minimize the post-vaccine IPD burden by applying Bayesian optimization to an ecological model of serotype replacement that integrated epidemiological and genomic data. We compared optimal formulations for reducing infant-only or population-wide IPD, and identified potential benefits to including non-conserved pneumococcal carrier proteins. Vaccines were also devised to minimize IPD resistant to antibiotic treatment, despite the ecological model assuming that resistance levels in the carried population would be preserved. We found that expanding infant-administered PCV valency is likely to result in diminishing returns, and that complementary pairs of infant- and adult-administered vaccines could be a superior strategy. PCV performance was highly dependent on the circulating pneumococcal population, further highlighting the advantages of a diversity of anti-pneumococcal vaccination strategies.

Phase variation in pneumococcal populations during carriage in the human nasopharynx

Streptococcus pneumoniae is one of the world's leading bacterial pathogens, responsible for pneumonia, septicaemia and meningitis. Asymptomatic colonisation of the nasopharynx is considered to be a prerequisite for these severe infections, however little is understood about the biological changes that permit the pneumococcus to switch from asymptomatic coloniser to invasive pathogen. A phase variable type I restriction-modification (R-M) system (SpnIII) has been linked to a change in capsule expression and to the ability to successfully colonise the murine nasopharynx. Using our laboratory data, we have developed a Markov change model that allows prediction of the expected level of phase variation within a population, and as a result measures when populations deviate from those expected at random. Using this model, we have analysed samples from the Experimental Human Pneumococcal Carriage (EHPC) project. Here we show, through mathematical modelling, that the patterns of dominant SpnIII alleles expressed in the human nasopharynx are significantly different than those predicted by stochastic switching alone. Our inter-disciplinary work demonstrates that the expression of alternative methylation patterns should be an important consideration in studies of pneumococcal colonisation.

Phylogenomic Analyses of Members of the Widespread Marine Heterotrophic Genus Pseudovibrio Suggest Distinct Evolutionary Trajectories and a Novel Genus, Polycladidibacter gen. nov

Bacteria belonging to the Pseudovibrio genus are widespread, metabolically versatile, and able to thrive as both free-living and host-associated organisms. Although more than 50 genomes are available, a comprehensive comparative genomics study to resolve taxonomic inconsistencies is currently missing. We analyzed all available genomes and used 552 core genes to perform a robust phylogenomic reconstruction. This in-depth analysis revealed the divergence of two monophyletic basal lineages of strains isolated from polyclad flatworm hosts, namely, Pseudovibrio hongkongensis and Pseudovibrio stylochi These strains have reduced genomes and lack sulfur-related metabolisms and major biosynthetic gene clusters, and their environmental distribution appears to be tightly associated with invertebrate hosts. We showed experimentally that the divergent strains are unable to utilize various sulfur compounds that, in contrast, can be utilized by the type strain Pseudovibrio denitrificans Our analyses suggest that the lineage leading to these two strains has been subject to relaxed purifying selection resulting in great gene loss. Overall genome relatedness indices (OGRI) indicate substantial differences between the divergent strains and the rest of the genus. While 16S rRNA gene analyses do not support the establishment of a different genus for the divergent strains, their substantial genomic, phylogenomic, and physiological differences strongly suggest a divergent evolutionary trajectory and the need for their reclassification. Therefore, we propose the novel genus Polycladidibacter gen. nov.IMPORTANCE The genus Pseudovibrio is commonly associated with marine invertebrates, which are essential for ocean health and marine nutrient cycling. Traditionally, the phylogeny of the genus has been based on 16S rRNA gene analysis. The use of the 16S rRNA gene or any other single marker gene for robust phylogenetic placement has recently been questioned. We used a large set of marker genes from all available Pseudovibrio genomes for in-depth phylogenomic analyses. We identified divergent monophyletic basal lineages within the Pseudovibrio genus, including two strains isolated from polyclad flatworms. These strains showed reduced sulfur metabolism and biosynthesis capacities. The phylogenomic analyses revealed distinct evolutionary trajectories and ecological adaptations that differentiate the divergent strains from the other Pseudovibrio members and suggest that they fall into a novel genus. Our data show the importance of widening the use of phylogenomics for better understanding bacterial physiology, phylogeny, and evolution.

The Contribution of Genetic Variation of Streptococcus pneumoniae to the Clinical Manifestation of Invasive Pneumococcal Disease

Background

Different clinical manifestations of invasive pneumococcal disease (IPD) have thus far mainly been explained by patient characteristics. Here we studied the contribution of pneumococcal genetic variation to IPD phenotype.

Methods

The index cohort consisted of 349 patients admitted to 2 Dutch hospitals between 2000-2011 with pneumococcal bacteremia. We performed genome-wide association studies to identify pneumococcal lineages, genes, and allelic variants associated with 23 clinical IPD phenotypes. The identified associations were validated in a nationwide (n = 482) and a post-pneumococcal vaccination cohort (n = 121). The contribution of confirmed pneumococcal genotypes to the clinical IPD phenotype, relative to known clinical predictors, was tested by regression analysis.

Results

Among IPD patients, the presence of pneumococcal gene slaA was a nationwide confirmed independent predictor of meningitis (odds ratio [OR], 10.5; P = .001), as was sequence cluster 9 (serotype 7F: OR, 3.68; P = .057). A set of 4 pneumococcal genes co-located on a prophage was a confirmed independent predictor of 30-day mortality (OR, 3.4; P = .003). We could detect the pneumococcal variants of concern in these patients' blood samples.

Conclusions

In this study, knowledge of pneumococcal genotypic variants improved the clinical risk assessment for detrimental manifestations of IPD. This provides us with novel opportunities to target, anticipate, or avert the pathogenic effects related to particular pneumococcal variants, and indicates that information on pneumococcal genotype is important for the diagnostic and treatment strategy in IPD. Ongoing surveillance is warranted to monitor the clinical value of information on pneumococcal variants in dynamic microbial and susceptible host populations.

The Best of All Worlds: Streptococcus pneumoniae Conjunctivitis through the Lens of Community Ecology and Microbial Biogeography

The study of the forces which govern the geographical distributions of life is known as biogeography, a subject which has fascinated zoologists, botanists and ecologists for centuries. Advances in our understanding of community ecology and biogeography—supported by rapid improvements in next generation sequencing technology—have now made it possible to identify and explain where and why life exists as it does, including within the microbial world. In this review, we highlight how a unified model of microbial biogeography, one which incorporates the classic ecological principles of selection, diversification, dispersion and ecological drift, can be used to explain community dynamics in the settings of both health and disease. These concepts operate on a multiplicity of temporal and spatial scales, and together form a powerful lens through which to study microbial population structures even at the finest anatomical resolutions. When applied specifically to curious strains of conjunctivitis-causing, nonencapsulated Streptococcus pneumoniae, we show how this conceptual framework can be used to explain the possible evolutionary and disease-causing mechanisms which allowed these lineages to colonize and invade a separate biogeography. An intimate knowledge of this radical bifurcation in phylogeny, still the only known niche subspecialization for S. pneumoniae to date, is critical to understanding the pathogenesis of ocular surface infections, nature of host-pathogen interactions, and developing strategies to curb disease transmission.

Gene gain and loss push prokaryotes beyond the homologous recombination barrier and accelerate genome sequence divergence

Bacterial and archaeal evolution involve extensive gene gain and loss. Thus, phylogenetic trees of prokaryotes can be constructed both by traditional sequence-based methods (gene trees) and by comparison of gene compositions (genome trees). Comparing the branch lengths in gene and genome trees with identical topologies for 34 clusters of closely related bacterial and archaeal genomes, we show here that terminal branches of gene trees are systematically compressed compared to those of genome trees. Thus, sequence evolution is delayed compared to genome evolution by gene gain and loss. The extent of this delay differs widely among bacteria and archaea. Mathematical modeling shows that the divergence delay can result from sequence homogenization by homologous recombination. The model explains how homologous recombination maintains the cohesiveness of the core genome of a species while allowing extensive gene gain and loss within the accessory genome. Once evolving genomes become isolated by barriers impeding homologous recombination, gene and genome evolution processes settle into parallel trajectories, and genomes diverge, resulting in speciation.

A significant proportion of the molecular evolution of bacteria and archaea occurs through gene gain and loss. Here Iranzo et al. develop a mathematical model that explains observed differential patterns of sequence evolution vs. gene content evolution as a consequence of homologous recombination.

Immunization with a recombinant antigen composed of conserved blocks from TSA56 provides broad genotype protection against scrub typhus

Scrub typhus is an acute febrile disease caused by Orientia tsutsugamushi infection. Despite the wide range of approaches explored during the last seventy years, an effective prophylactic vaccine is not yet available. Here, we developed a novel recombinant antigen derived from conserved regions of 56 kDa type-specific antigen (TSA56), a major outer membrane protein responsible for genetic heterogeneity and antigenicity, and evaluated it as a protective vaccine antigen. Our findings demonstrate that immunization with conserved blocks of TSA56 (cTSA56) not only provides protective immunity against lethal challenges with the homologous genotype, but also confers significantly better protection against heterologous genotypes than TSA56. Adoptive transfer of CD4⁺ or CD8⁺ T cells from immunized mice provided significantly enhanced protection against lethal challenge, whereas immune B cells failed to do so, indicating that cellular immunity against the conserved epitopes plays a protective role. Moreover, immunization with a 10-mer peptide mixture, screened from CD8⁺ T cell epitopes within the conserved region of TSA56, provided enhanced protection against lethal challenge with O. tsutsugamushi. Therefore, this novel recombinant antigen is a promising candidate for scrub typhus vaccine against a wide range of O. tsutsugamushi genotypes.

Population Gene Introgression and High Genome Plasticity for the Zoonotic Pathogen Streptococcus agalactiae

The influence that bacterial adaptation (or niche partitioning) within species has on gene spillover and transmission among bacterial populations occupying different niches is not well understood. Streptococcus agalactiae is an important bacterial pathogen that has a taxonomically diverse host range making it an excellent model system to study these processes. Here, we analyze a global set of 901 genome sequences from nine diverse host species to advance our understanding of these processes. Bayesian clustering analysis delineated 12 major populations that closely aligned with niches. Comparative genomics revealed extensive gene gain/loss among populations and a large pan genome of 9,527 genes, which remained open and was strongly partitioned among niches. As a result, the biochemical characteristics of 11 populations were highly distinctive (significantly enriched). Positive selection was detected and biochemical characteristics of the dispensable genes under selection were enriched in ten populations. Despite the strong gene partitioning, phylogenomics detected gene spillover. In particular, tetracycline resistance (which likely evolved in the human-associated population) from humans to bovine, canines, seals, and fish, demonstrating how a gene selected in one host can ultimately be transmitted into another, and biased transmission from humans to bovines was confirmed with a Bayesian migration analysis. Our findings show high bacterial genome plasticity acting in balance with selection pressure from distinct functional requirements of niches that is associated with an extensive and highly partitioned dispensable genome, likely facilitating continued and expansive adaptation.

Genetic variation regulates the activation and specificity of Restriction-Modification systems in Neisseria gonorrhoeae

Restriction-Modification systems (RMS) are one of the main mechanisms of defence against foreign DNA invasion and can have an important role in the regulation of gene expression. The obligate human pathogen Neisseria gonorrhoeae carries one of the highest loads of RMS in its genome; between 13 to 15 of the three main types. Previous work has described their organization in the reference genome FA1090 and has inferred the associated methylated motifs. Here, we studied the structure of RMS and target methylated motifs in 25 gonococcal strains sequenced with Single Molecule Real-Time (SMRT) technology, which provides data on DNA modification. The results showed a variable picture of active RMS in different strains, with phase variation switching the activity of Type III RMS, and both the activity and specificity of a Type I RMS. Interestingly, the Dam methylase was found in place of the NgoAXI endonuclease in two of the strains, despite being previously thought to be absent in the gonococcus. We also identified the real methylation target of NgoAXII as 5′-GCAGA-3′, different from that previously described. Results from this work give further insights into the diversity and dynamics of RMS and methylation patterns in N. gonorrhoeae.

Synergistic Activity of Mobile Genetic Element Defences in Streptococcus pneumoniae

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

Methylation Warfare: Interaction of Pneumococcal Bacteriophages with Their Host

With antimicrobial drug resistance becoming an increasing burden on human health, much attention has been focused on the potential use of bacteriophages and their enzymes as therapeutics. However, the investigations into the physiology of the complex interactions of bacteriophages with their hosts have attracted far less attention, in comparison. This work describes the molecular characterization of the infectious cycle of a bacteriophage in the important human pathogen Streptococcus pneumoniae and explores the intricate relationship between phase-variable host defense mechanisms and the virus. This is the first report showing how a phase-variable type I restriction-modification system is involved in bacteriophage restriction while it also provides an additional level of infection control through abortive infection.

Virus-host interactions are regulated by complex coevolutionary dynamics. In Streptococcus pneumoniae, phase-variable type I restriction-modification (R-M) systems are part of the core genome. We hypothesized that the ability of the R-M systems to switch between six target DNA specificities also has a key role in preventing the spread of bacteriophages. Using the streptococcal temperate bacteriophage SpSL1, we show that the variants of both the SpnIII and SpnIV R-M systems are able to restrict invading bacteriophage with an efficiency approximately proportional to the number of target sites in the bacteriophage genome. In addition to restriction of lytic replication, SpnIII also led to abortive infection in the majority of host cells. During lytic infection, transcriptional analysis found evidence of phage-host interaction through the strong upregulation of the nrdR nucleotide biosynthesis regulon. During lysogeny, the phage had less of an effect on host gene regulation. This research demonstrates a novel combined bacteriophage restriction and abortive infection mechanism, highlighting the importance that the phase-variable type I R-M systems have in the multifunctional defense against bacteriophage infection in the respiratory pathogen S. pneumoniae.

IMPORTANCE With antimicrobial drug resistance becoming an increasing burden on human health, much attention has been focused on the potential use of bacteriophages and their enzymes as therapeutics. However, the investigations into the physiology of the complex interactions of bacteriophages with their hosts have attracted far less attention, in comparison. This work describes the molecular characterization of the infectious cycle of a bacteriophage in the important human pathogen Streptococcus pneumoniae and explores the intricate relationship between phase-variable host defense mechanisms and the virus. This is the first report showing how a phase-variable type I restriction-modification system is involved in bacteriophage restriction while it also provides an additional level of infection control through abortive infection.