Robust analysis of prokaryotic pangenome gene gain and loss rates with Panstripe.

Horizontal gene transfer (HGT) plays a critical role in the evolution and diversification of many microbial species. The resulting dynamics of gene gain and loss can have important implications for the development of antibiotic resistance and the design of vaccine and drug interventions. Methods for the analysis of gene presence/absence patterns typically do not account for errors introduced in the automated annotation and clustering of gene sequences. In particular, methods adapted from ecological studies, including the pangenome gene accumulation curve, can be misleading as they may reflect the underlying diversity in the temporal sampling of genomes rather than a difference in the dynamics of HGT. Here, we introduce Panstripe, a method based on generalized linear regression that is robust to population structure, sampling bias, and errors in the predicted presence/absence of genes. We show using simulations that Panstripe can effectively identify differences in the rate and number of genes involved in HGT events, and illustrate its capability by analyzing several diverse bacterial genome data sets representing major human pathogens.

Dissemination of extensively drug-resistant NDM-producing Providencia stuartii in Europe linked to patients transferred from Ukraine, March 2022 to March 2023.

BackgroundThe war in Ukraine led to migration of Ukrainian people. Early 2022, several European national surveillance systems detected multidrug-resistant (MDR) bacteria related to Ukrainian patients.AimTo investigate the genomic epidemiology of New Delhi metallo-ß-lactamase (NDM)-producing Providencia stuartii from Ukrainian patients among European countries.MethodsWhole-genome sequencing of 66 isolates sampled in 2022-2023 in 10 European countries enabled whole-genome multilocus sequence typing (wgMLST), identification of resistance genes, replicons, and plasmid reconstructions. Five bla NDM-1-carrying-P. stuartii isolates underwent antimicrobial susceptibility testing (AST). Transferability to Escherichia coli of a bla NDM-1-carrying plasmid from a patient strain was assessed. Epidemiological characteristics of patients with NDM-producing P. stuartii were gathered by questionnaire.ResultswgMLST of the 66 isolates revealed two genetic clusters unrelated to Ukraine and three linked to Ukrainian patients. Of these three, two comprised bla NDM-1-carrying-P. stuartii and the third bla NDM-5-carrying-P. stuartii. The bla NDM-1 clusters (PstCluster-001, n = 22 isolates; PstCluster-002, n = 8 isolates) comprised strains from seven and four countries, respectively. The bla NDM-5 cluster (PstCluster-003) included 13 isolates from six countries. PstCluster-001 and PstCluster-002 isolates carried an MDR plasmid harbouring bla NDM-1, bla OXA-10, bla CMY-16, rmtC and armA, which was transferrable in vitro and, for some Ukrainian patients, shared by other Enterobacterales. AST revealed PstCluster-001 isolates to be extensively drug-resistant (XDR), but susceptible to cefiderocol and aztreonam-avibactam. Patients with data on age (n = 41) were 19-74 years old; of 49 with information on sex, 38 were male.ConclusionXDR P. stuartii were introduced into European countries, requiring increased awareness and precautions when treating patients from conflict-affected areas.

Rapid cross-border emergence of NDM-5-producing Escherichia coli in the European Union/European Economic Area, 2012 to June 2022.

Whole genome sequencing data of 874 Escherichia coli isolates carrying bla NDM-5 from 13 European Union/European Economic Area countries between 2012 and June 2022 showed the predominance of sequence types ST167, ST405, ST410, ST361 and ST648, and an increasing frequency of detection. Nearly a third (30.6%) of these isolates were associated with infections and more than half (58.2%) were predicted to be multidrug-resistant. Further spread of E. coli carrying bla NDM-5 would leave limited treatment options for serious E. coli infections.

Adaptive cell wall thickening in Enterococcus faecalis is associated with decreased vancomycin susceptibility.

OBJECTIVES: Enterococcus faecalis can adopt both a commensal and a nosocomial lifestyle, resisting numerous antibiotics. In this study, we aim to investigate the relationship between the cell wall (CW) thickness and decreased susceptibility to vancomycin (VD) in van-gene negative clinical isolates of E. faecalis (nMIC 8 = 2, nMIC 4 = 3, ST30, ST40, and ST59). METHODS: The CW thickness was assessed in VD strains and compared with vancomycin susceptible isolates of the same sequence type (ST) (Vancomycin susceptible [VS]; nMIC 2 = 5). The VD and VS strains were subjected to serial passage (evolved [ev]) with and without vancomycin selection. Subsequent measurements of CW thickness and vancomycin MICs were performed. RESULTS: The VD strains exhibited increased CW thickness when compared with ST-related VS strains (ΔCW thickness VD vs. VS ST30 25 nm, ST59 15 nm, and ST40 1 nm). Serial passages without vancomycin selection led to a decrease in CW thickness and vancomycin MIC in VD strains (ΔCW thickness VD vs. evVD ST30 22 nm, ST59 3 nm, and ST40 2 nm). Serial passages with vancomycin selection caused an increase in CW thickness and vancomycin MIC in ST-related VS strains (ΔCW thickness VS vs. evVS ST30 22 nm, ST59 16 nm, and ST40 1 nm). DISCUSSION: Adaptive changes in CW thickness were observed in response to vancomycin exposure. Increased CW thickness correlated with decreased vancomycin susceptibility, whereas decreased CW thickness correlated with increased vancomycin susceptibility. Core single nucleotide polymorphisms in the evolved mutants were mostly found in genes encoding proteins associated with the cytoplasm or the cytoplasmic membrane. The potential relevance of these adaptive changes is underlined by the observed phenotypes in clinical isolates. Our findings emphasize the importance of monitoring adaptive changes, as vancomycin-resistant enterococci infections are a growing concern.

The first tigecycline resistant Enterococcus faecium in Norway was related to tigecycline exposure.

OBJECTIVES: We describe the first tigecycline resistant enterococcal isolate in Norway and the mechanisms involved. MATERIAL AND METHODS: The Norwegian National Advisory Unit on Detection of Antimicrobial Resistance (K-res). received in 2022 an Enterococcus faecium blood culture isolate with decreased susceptibility to tigecycline from a hospitalized patient in the South-Eastern Norway Health region for confirmatory testing. K-res verified a tigecycline-resistant E. faecium (TigR) with broth microdilution MIC of 0.5 mg/L. The patient had received treatment with tigecycline because of an infection with a linezolid- and vancomycin-resistant but tigecycline susceptible E. faecium (TigS) 47 days prior to the detection of the corresponding tigecycline-resistant isolate. Whole-genome comparisons, cgMLST and SNP analyses revealed that the two ST117 strains were closely related. RESULTS: The TigR isolate showed a novel deletion of 2 amino acids (K57Y58) in a polymorphic region of ribosomal protein S10 previously associated with tigecycline resistance and a deletion of the tet(M) leader peptide previously related to increased expression of tet(M) and tigecycline resistance in enterococci. CONCLUSIONS: Genomic and epidemiological analyses confirm that the two E. faecium (TigR and TigS) are closely related isolates of the same strain and that the two deletions (in rpsJ and of tet(M) leader peptide) account for the tigecycline resistance in TigR.

Detecting co-selection through excess linkage disequilibrium in bacterial genomes.

Population genomics has revolutionized our ability to study bacterial evolution by enabling data-driven discovery of the genetic architecture of trait variation. Genome-wide association studies (GWAS) have more recently become accompanied by genome-wide epistasis and co-selection (GWES) analysis, which offers a phenotype-free approach to generating hypotheses about selective processes that simultaneously impact multiple loci across the genome. However, existing GWES methods only consider associations between distant pairs of loci within the genome due to the strong impact of linkage-disequilibrium (LD) over short distances. Based on the general functional organisation of genomes it is nevertheless expected that majority of co-selection and epistasis will act within relatively short genomic proximity, on co-variation occurring within genes and their promoter regions, and within operons. Here, we introduce LDWeaver, which enables an exhaustive GWES across both short- and long-range LD, to disentangle likely neutral co-variation from selection. We demonstrate the ability of LDWeaver to efficiently generate hypotheses about co-selection using large genomic surveys of multiple major human bacterial pathogen species and validate several findings using functional annotation and phenotypic measurements. Our approach will facilitate the study of bacterial evolution in the light of rapidly expanding population genomic data.

Deep sequencing of Escherichia coli exposes colonisation diversity and impact of antibiotics in Punjab, Pakistan.

Multi-drug resistant (MDR) E. coli constitute a major public health burden globally, reaching the highest prevalence in the global south yet frequently flowing with travellers to other regions. However, our comprehension of the entire genetic diversity of E. coli colonising local populations remains limited. We quantified this diversity, its associated antimicrobial resistance (AMR), and assessed the impact of antibiotic use by recruiting 494 outpatients and 423 community dwellers in the Punjab province, Pakistan. Rectal swab and stool samples were cultured on CLED agar and DNA extracted from plate sweeps was sequenced en masse to capture both the genetic and AMR diversity of E. coli. We assembled 5,247 E. coli genomes from 1,411 samples, displaying marked genetic diversity in gut colonisation. Compared with high income countries, the Punjabi population generally showed a markedly different distribution of genetic lineages and AMR determinants, while use of antibiotics elevated the prevalence of well-known globally circulating MDR clinical strains. These findings implicate that longitudinal multi-regional genomics-based surveillance of both colonisation and infections is a prerequisite for developing mechanistic understanding of the interplay between ecology and evolution in the maintenance and dissemination of (MDR) E. coli.

Modulation of multidrug-resistant clone success in Escherichia coli populations: a longitudinal, multi-country, genomic and antibiotic usage cohort study.

BACKGROUND: The effect of antibiotic usage on the success of multidrug-resistant (MDR) clones in a population remains unclear. With this genomics-based molecular epidemiology study, we aimed to investigate the contribution of antibiotic use to Escherichia coli clone success, relative to intra-strain competition for colonisation and infection. METHODS: We sequenced all the available E coli bloodstream infection isolates provided by the British Society for Antimicrobial Chemotherapy (BSAC) from 2012 to 2017 (n=718) and combined these with published data from the UK (2001-11; n=1090) and Norway (2002-17; n=3254). Defined daily dose (DDD) data from the European Centre for Disease Prevention and Control (retrieved on Sept 21, 2021) for major antibiotic classes (ß-lactam, tetracycline, macrolide, sulfonamide, quinolone, and non-penicillin ß-lactam) were used together with sequence typing, resistance profiling, regression analysis, and non-neutral Wright-Fisher simulation-based modelling to enable systematic comparison of resistance levels, clone success, and antibiotic usage between the UK and Norway. FINDINGS: Sequence type (ST)73, ST131, ST95, and ST69 accounted for 892 (49·3%) of 1808 isolates in the BSAC collection. In the UK, the proportion of ST69 increased between 2001-10 and 2011-17 (p=0·0004), whereas the proportions of ST73 and ST95 did not vary between periods. ST131 expanded quickly after its emergence in 2003 and its prevalence remained consistent throughout the study period (apart from a brief decrease in 2009-10). The extended-spectrum ß-lactamase (ESBL)-carrying, globally disseminated MDR clone ST131-C2 showed overall greater success in the UK (154 [56·8%] of 271 isolates in 2003-17) compared with Norway (51 [18·3%] of 278 isolates in 2002-17; p<0·0001). DDD data indicated higher total use of antimicrobials in the UK, driven mainly by the class of non-penicillin ß-lactams, which were used between 2·7-times and 5·1-times more in the UK per annum (ratio mean 3·7 [SD 0·8]). This difference was associated with the higher success of the MDR clone ST131-C2 (pseudo-R2 69·1%). A non-neutral Wright-Fisher model replicated the observed expansion of non-MDR and MDR sequence types under higher DDD regimes. INTERPRETATION: Our study indicates that resistance profiles of contemporaneously successful clones can vary substantially, warranting caution in the interpretation of correlations between aggregate measures of resistance and antibiotic usage. Our study further suggests that in countries with low-to-moderate use of antibiotics, such as the UK and Norway, the extent of non-penicillin ß-lactam use modulates rather than determines the success of widely disseminated MDR ESBL-carrying E coli clones. Detailed understanding of underlying causal drivers of success is important for improved control of resistant pathogens. FUNDING: Trond Mohn Foundation, Marie Sklodowska-Curie Actions, European Research Council, Royal Society, and Wellcome Trust.

Mge-cluster: a reference-free approach for typing bacterial plasmids.

Extrachromosomal elements of bacterial cells such as plasmids are notorious for their importance in evolution and adaptation to changing ecology. However, high-resolution population-wide analysis of plasmids has only become accessible recently with the advent of scalable long-read sequencing technology. Current typing methods for the classification of plasmids remain limited in their scope which motivated us to develop a computationally efficient approach to simultaneously recognize novel types and classify plasmids into previously identified groups. Here, we introduce mge-cluster that can easily handle thousands of input sequences which are compressed using a unitig representation in a de Bruijn graph. Our approach offers a faster runtime than existing algorithms, with moderate memory usage, and enables an intuitive visualization, classification and clustering scheme that users can explore interactively within a single framework. Mge-cluster platform for plasmid analysis can be easily distributed and replicated, enabling a consistent labelling of plasmids across past, present, and future sequence collections. We underscore the advantages of our approach by analysing a population-wide plasmid data set obtained from the opportunistic pathogen Escherichia coli, studying the prevalence of the colistin resistance gene mcr-1.1 within the plasmid population, and describing an instance of resistance plasmid transmission within a hospital environment.

The population-level impact of Enterococcus faecalis genetics on intestinal colonization and extraintestinal infection.

IMPORTANCE: Enterococcus faecalis causes life-threatening invasive hospital- and community-associated infections that are usually associated with multidrug resistance globally. Although E. faecalis infections cause opportunistic infections typically associated with antibiotic use, immunocompromised immune status, and other factors, they also possess an arsenal of virulence factors crucial for their pathogenicity. Despite this, the relative contribution of these virulence factors and other genetic changes to the pathogenicity of E. faecalis strains remain poorly understood. Here, we investigated whether specific genomic changes in the genome of E. faecalis isolates influence its pathogenicity-infection of hospitalized and nonhospitalized individuals and the propensity to cause extraintestinal infection and intestinal colonization. Our findings indicate that E. faecalis genetics partially influence the infection of hospitalized and nonhospitalized individuals and the propensity to cause extraintestinal infection, possibly due to gut-to-bloodstream translocation, highlighting the potential substantial role of host and environmental factors, including gut microbiota, on the opportunistic pathogenic lifestyle of this bacterium.

Evolutionary and functional history of the Escherichia coli K1 capsule.

Escherichia coli is a leading cause of invasive bacterial infections in humans. Capsule polysaccharide has an important role in bacterial pathogenesis, and the K1 capsule has been firmly established as one of the most potent capsule types in E. coli through its association with severe infections. However, little is known about its distribution, evolution and functions across the E. coli phylogeny, which is fundamental to elucidating its role in the expansion of successful lineages. Using systematic surveys of invasive E. coli isolates, we show that the K1-cps locus is present in a quarter of bloodstream infection isolates and has emerged in at least four different extraintestinal pathogenic E. coli (ExPEC) phylogroups independently in the last 500 years. Phenotypic assessment demonstrates that K1 capsule synthesis enhances E. coli survival in human serum independent of genetic background, and that therapeutic targeting of the K1 capsule re-sensitizes E. coli from distinct genetic backgrounds to human serum. Our study highlights that assessing the evolutionary and functional properties of bacterial virulence factors at population levels is important to better monitor and predict the emergence of virulent clones, and to also inform therapies and preventive medicine to effectively control bacterial infections whilst significantly lowering antibiotic usage.

Strong pathogen competition in neonatal gut colonisation.

Opportunistic bacterial pathogen species and their strains that colonise the human gut are generally understood to compete against both each other and the commensal species colonising this ecosystem. Currently we are lacking a population-wide quantification of strain-level colonisation dynamics and the relationship of colonisation potential to prevalence in disease, and how ecological factors might be modulating these. Here, using a combination of latest high-resolution metagenomics and strain-level genomic epidemiology methods we performed a characterisation of the competition and colonisation dynamics for a longitudinal cohort of neonatal gut microbiomes. We found strong inter- and intra-species competition dynamics in the gut colonisation process, but also a number of synergistic relationships among several species belonging to genus Klebsiella, which includes the prominent human pathogen Klebsiella pneumoniae. No evidence of preferential colonisation by hospital-adapted pathogen lineages in either vaginal or caesarean section birth groups was detected. Our analysis further enabled unbiased assessment of strain-level colonisation potential of extra-intestinal pathogenic Escherichia coli (ExPEC) in comparison with their propensity to cause bloodstream infections. Our study highlights the importance of systematic surveillance of bacterial gut pathogens, not only from disease but also from carriage state, to better inform therapies and preventive medicine in the future.

A high-throughput multiplexing and selection strategy to complete bacterial genomes.

BACKGROUND: Bacterial whole-genome sequencing based on short-read technologies often results in a draft assembly formed by contiguous sequences. The introduction of long-read sequencing technologies permits those contiguous sequences to be unambiguously bridged into complete genomes. However, the elevated costs associated with long-read sequencing frequently limit the number of bacterial isolates that can be long-read sequenced. Here we evaluated the recently released 96 barcoding kit from Oxford Nanopore Technologies (ONT) to generate complete genomes on a high-throughput basis. In addition, we propose an isolate selection strategy that optimizes a representative selection of isolates for long-read sequencing considering as input large-scale bacterial collections. RESULTS: Despite an uneven distribution of long reads per barcode, near-complete chromosomal sequences (assembly contiguity = 0.89) were generated for 96 Escherichia coli isolates with associated short-read sequencing data. The assembly contiguity of the plasmid replicons was even higher (0.98), which indicated the suitability of the multiplexing strategy for studies focused on resolving plasmid sequences. We benchmarked hybrid and ONT-only assemblies and showed that the combination of ONT sequencing data with short-read sequencing data is still highly desirable (i) to perform an unbiased selection of isolates for long-read sequencing, (ii) to achieve an optimal genome accuracy and completeness, and (iii) to include small plasmids underrepresented in the ONT library. CONCLUSIONS: The proposed long-read isolate selection ensures the completion of bacterial genomes that span the genome diversity inherent in large collections of bacterial isolates. We show the potential of using this multiplexing approach to close bacterial genomes on a high-throughput basis.

Apparent nosocomial adaptation of Enterococcus faecalis predates the modern hospital era.

Enterococcus faecalis is a commensal and nosocomial pathogen, which is also ubiquitous in animals and insects, representing a classical generalist microorganism. Here, we study E. faecalis isolates ranging from the pre-antibiotic era in 1936 up to 2018, covering a large set of host species including wild birds, mammals, healthy humans, and hospitalised patients. We sequence the bacterial genomes using short- and long-read techniques, and identify multiple extant hospital-associated lineages, with last common ancestors dating back as far as the 19th century. We find a population cohesively connected through homologous recombination, a metabolic flexibility despite a small genome size, and a stable large core genome. Our findings indicate that the apparent hospital adaptations found in hospital-associated E. faecalis lineages likely predate the "modern hospital" era, suggesting selection in another niche, and underlining the generalist nature of this nosocomial pathogen.

Emergence and dissemination of antimicrobial resistance in Escherichia coli causing bloodstream infections in Norway in 2002-17: a nationwide, longitudinal, microbial population genomic study.

BACKGROUND: The clonal diversity underpinning trends in multidrug resistant Escherichia coli causing bloodstream infections remains uncertain. We aimed to determine the contribution of individual clones to resistance over time, using large-scale genomics-based molecular epidemiology. METHODS: This was a longitudinal, E coli population, genomic, cohort study that sampled isolates from 22 512 E coli bloodstream infections included in the Norwegian surveillance programme on resistant microbes (NORM) from 2002 to 2017. 15 of 22 laboratories were able to share their isolates, and the first 22·5% of isolates from each year were requested. We used whole genome sequencing to infer the population structure (PopPUNK), and we investigated the clade composition of the dominant multidrug resistant clonal complex (CC)131 using genetic markers previously reported for sequence type (ST)131, effective population size (BEAST), and presence of determinants of antimicrobial resistance (ARIBA, PointFinder, and ResFinder databases) over time. We compared these features between the 2002-10 and 2011-17 time periods. We also compared our results with those of a longitudinal study from the UK done between 2001 and 2011. FINDINGS: Of the 3500 isolates requested from the participating laboratories, 3397 (97·1%) were received, of which 3254 (95·8%) were successfully sequenced and included in the analysis. A significant increase in the number of multidrug resistant CC131 isolates from 71 (5·6%) of 1277 in 2002-10 to 207 (10·5%) of 1977 in 2011-17 (p<0·0001), was the largest clonal expansion. CC131 was the most common clone in extended-spectrum ß-lactamase (ESBL)-positive isolates (75 [58·6%] of 128) and fluoroquinolone non-susceptible isolates (148 [39·2%] of 378). Within CC131, clade A increased in prevalence from 2002, whereas the global multidrug resistant clade C2 was not observed until 2007. Multiple de-novo acquisitions of both blaCTX-M ESBL-encoding genes in clades A and C1 and gain of phenotypic fluoroquinolone non-susceptibility across the clade A phylogeny were observed. We estimated that exponential increases in the effective population sizes of clades A, C1, and C2 occurred in the mid-2000s, and in clade B a decade earlier. The rate of increase in the estimated effective population size of clade A (Ne=3147) was nearly ten-times that of C2 (Ne=345), with clade A over-represented in Norwegian CC131 isolates (75 [27·0%] of 278) compared with the UK study (8 [5·4%] of 147 isolates). INTERPRETATION: The early and sustained establishment of predominantly antimicrobial susceptible CC131 clade A isolates, relative to multidrug resistant clade C2 isolates, suggests that resistance is not necessary for clonal success. However, even in the low antibiotic use setting of Norway, resistance to important antimicrobial classes has rapidly been selected for in CC131 clade A isolates. This study shows the importance of genomic surveillance in uncovering the complex ecology underlying multidrug resistance dissemination and competition, which have implications for the design of strategies and interventions to control the spread of high-risk multidrug resistant clones. FUNDING: Trond Mohn Foundation, European Research Council, Marie Sklodowska-Curie Actions, and the Wellcome Trust.