PANINI: Pangenome Neighbour Identification for Bacterial Populations

Genomic epidemiology of extended-spectrum beta-lactamase-producing Escherichia coli from humans and a river in Aotearoa New Zealand

In Aotearoa New Zealand, urinary tract infections in humans are commonly caused by extended-spectrum beta-lactamase (ESBL)-producing Escherichia coli. This group of antimicrobial-resistant bacteria are often multidrug resistant. However, there is limited information on ESBL-producing E. coli found in the environment and their link with human clinical isolates. In this study, we examined the genetic relationship between environmental and human clinical ESBL-producing E. coli and isolates collected in parallel within the same area over 14 months. Environmental samples were collected from treated effluent, stormwater and multiple locations along an Aotearoa New Zealand river. Treated effluent, stormwater and river water sourced downstream of the treated effluent outlet were the main samples that were positive for ESBL-producing E. coli (7/14 samples, 50.0%; 3/6 samples, 50%; and 15/28 samples, 54%, respectively). Whole-genome sequence comparison was carried out on 307 human clinical and 45 environmental ESBL-producing E. coli isolates. Sequence type 131 was dominant for both clinical (147/307, 47.9%) and environmental isolates (11/45, 24.4%). Only one ESBL gene was detected in each isolate. Among the clinical isolates, the most prevalent ESBL genes were bla CTX-M-27 (134/307, 43.6%) and bla CTX-M-15 (134/307, 43.6%). Among the environmental isolates, bla CTX-M-15 (28/45, 62.2%) was the most prevalent gene. A core SNP analysis of these isolates suggested that some strains were shared between humans and the local river. These results highlight the importance of understanding different transmission pathways for the spread of ESBL-producing E. coli.

Pan-genomics: Insight into the Functional Genome, Applications, Advancements, and Challenges

A pan-genome is a compilation of the common and unique genomes found in a given species. It incorporates the genetic information from all of the genomes sampled, producing a big and diverse set of genetic material. Pan-genomic analysis has various advantages over typical genomics research. It creates a vast and varied spectrum of genetic material by combining the genetic data from all the sampled genomes. Comparing pan-genomics analysis to conventional genomic research, there are a number of benefits. Although the most recent era of pan-genomic studies has used cutting-edge sequencing technology to shed fresh light on biological variety and improvement, the potential uses of pan-genomics in improvement have not yet been fully realized. Pan-genome research in various organisms has demonstrated that missing genetic components and the detection of significant Structural Variants (SVs) can be investigated using pan-genomic methods. Many individual-specific sequences have been linked to biological adaptability, phenotypic, and key economic attributes. This study aims to focus on how pangenome analysis uncovers genetic differences in various organisms, including human, and their effects on phenotypes, as well as how this might help us comprehend the diversity of species. The review also concentrated on potential problems and the prospects for future pangenome research.

Comparative pangenome analysis of Enterococcus faecium and Enterococcus lactis provides new insights into the adaptive evolution by horizontal gene acquisitions

BACKGROUND

Enterococcus faecium and E. lactis are phylogenetically closely related lactic acid bacteria that are ubiquitous in nature and are known to be beneficial or pathogenic. Despite their considerable industrial and clinical importance, comprehensive studies on their evolutionary relationships and genomic, metabolic, and pathogenic traits are still lacking. Therefore, we conducted comparative pangenome analyses using all available dereplicated genomes of these species.

RESULTS

E. faecium was divided into two subclades: subclade I, comprising strains derived from humans, animals, and food, and the more recent phylogenetic subclade II, consisting exclusively of human-derived strains. In contrast, E. lactis strains, isolated from diverse sources including foods, humans, animals, and the environment, did not display distinct clustering based on their isolation sources. Despite having similar metabolic features, noticeable genomic differences were observed between E. faecium subclades I and II, as well as E. lactis. Notably, E. faecium subclade II strains exhibited significantly larger genome sizes and higher gene counts compared to both E. faecium subclade I and E. lactis strains. Furthermore, they carried a higher abundance of antibiotic resistance, virulence, bacteriocin, and mobile element genes. Phylogenetic analysis of antibiotic resistance and virulence genes suggests that E. faecium subclade II strains likely acquired these genes through horizontal gene transfer, facilitating their effective adaptation in response to antibiotic use in humans.

CONCLUSIONS

Our study offers valuable insights into the adaptive evolution of E. faecium strains, enabling their survival as pathogens in the human environment through horizontal gene acquisitions.

Staphylococcus aureus populations from the gut and the blood are not distinguished by virulence traits-a critical role of host barrier integrity

BACKGROUND

The opportunistic pathogen Staphylococcus aureus is an asymptomatically carried member of the microbiome of about one third of the human population at any given point in time. Body sites known to harbor S. aureus are the skin, nasopharynx, and gut. In particular, the mechanisms allowing S. aureus to pass the gut epithelial barrier and to invade the bloodstream were so far poorly understood. Therefore, the objective of our present study was to investigate the extent to which genetic differences between enteric S. aureus isolates and isolates that caused serious bloodstream infections contribute to the likelihood of invasive disease.

RESULTS

Here, we present genome-wide association studies (GWAS) that compare the genome sequences of 69 S. aureus isolates from enteric carriage by healthy volunteers and 95 isolates from bloodstream infections. We complement our GWAS results with a detailed characterization of the cellular and extracellular proteomes of the representative gut and bloodstream isolates, and by assaying the virulence of these isolates with infection models based on human gut epithelial cells, human blood cells, and a small animal infection model. Intriguingly, our results show that enteric and bloodstream isolates with the same sequence type (ST1 or ST5) are very similar to each other at the genomic and proteomic levels. Nonetheless, bloodstream isolates are not necessarily associated with an invasive profile. Furthermore, we show that the main decisive factor preventing infection of gut epithelial cells in vitro is the presence of a tight barrier.

CONCLUSIONS

Our data show that virulence is a highly variable trait, even within a single clone. Importantly, however, there is no evidence that blood stream isolates possess a higher virulence potential than those from the enteric carriage. In fact, some gut isolates from healthy carriers were more virulent than bloodstream isolates. Based on our present observations, we propose that the integrity of the gut epithelial layer, rather than the pathogenic potential of the investigated enteric S. aureus isolates, determines whether staphylococci from the gut microbiome will become invasive pathogens. Video Abstract.

Dynamics of extended-spectrum cephalosporin resistance genes in Escherichia coli from Europe and North America

Extended-spectrum cephalosporins (ESCs) are critically important antimicrobial agents for human and veterinary medicine. ESC resistance (ESC-R) genes have spread worldwide through plasmids and clonal expansion, yet the distribution and dynamics of ESC-R genes in different ecological compartments are poorly understood. Here we use whole genome sequence data of Enterobacterales isolates of human and animal origin from Europe and North America and identify contrasting temporal dynamics. AmpC β-lactamases were initially more dominant in North America in humans and farm animals, only later emerging in Europe. In contrast, specific extended-spectrum β-lactamases (ESBLs) were initially common in animals from Europe and later emerged in North America. This study identifies differences in the relative importance of plasmids and clonal expansion across different compartments for the spread of different ESC-R genes. Understanding the mechanisms of transmission will be critical in the design of interventions to reduce the spread of antimicrobial resistance.

Comparative genomics and pangenomics of vancomycin-resistant and susceptible Enterococcus faecium from Irish hospitals

Introduction. Enterococcus faecium has emerged as an important nosocomial pathogen, which is increasingly difficult to treat due to the genetic acquisition of vancomycin resistance. Ireland has a recalcitrant vancomycin-resistant bloodstream infection rate compared to other developed countries.

Hypothesis/Gap statement. Vancomycin resistance rates persist amongst E. faecium isolates from Irish hospitals. The evolutionary genomics governing these trends have not been fully elucidated.

Methodology. A set of 28 vancomycin-resistant isolates was sequenced to construct a dataset alongside 61 other publicly available Irish genomes. This dataset was extensively analysed using in silico methodologies (comparative genomics, pangenomics, phylogenetics, genotypics and comparative functional analyses) to uncover distinct evolutionary, coevolutionary and clinically relevant population trends.

Results. These results suggest that a stable (in terms of genome size, GC% and number of genes), yet genetically diverse population (in terms of gene content) of E. faecium persists in Ireland with acquired resistance arising via plasmid acquisition (vanA) or, to a lesser extent, chromosomal recombination (vanB). Population analysis revealed five clusters with one cluster partitioned into four clades which transcend isolation dates. Pangenomic and recombination analyses revealed an open (whole genome and chromosomal specific) pangenome illustrating a rampant evolutionary pattern. Comparative resistomics and virulomics uncovered distinct chromosomal and mobilomal propensity for multidrug resistance, widespread chromosomal point-mutation-mediated resistance and chromosomally harboured arsenals of virulence factors. Interestingly, a potential difference in biofilm formation strategies was highlighted by coevolutionary analysis, suggesting differential biofilm genotypes between vanA and vanB isolates.

Conclusions. These results highlight the evolutionary history of Irish E. faecium isolates and may provide insight into underlying infection dynamics in a clinical setting. Due to the apparent ease of vancomycin resistance acquisition over time, susceptible E. faecium should be concurrently reduced in Irish hospitals to mitigate potential resistant infections.

Whole-Genome Sequence Comparisons of Listeria monocytogenes Isolated from Meat and Fish Reveal High Inter- and Intra-Sample Diversity

Interpretation of whole-genome sequencing (WGS) data for foodborne outbreak investigations is complex, as the genetic diversity within processing plants and transmission events need to be considered. In this study, we analyzed 92 food-associated Listeria monocytogenes isolates by WGS-based methods. We aimed to examine the genetic diversity within meat and fish production chains and to assess the applicability of suggested thresholds for clustering of potentially related isolates. Therefore, meat-associated isolates originating from the same samples or processing plants as well as fish-associated isolates were analyzed as distinct sets. In silico serogrouping, multilocus sequence typing (MLST), core genome MLST (cgMLST), and pangenome analysis were combined with screenings for prophages and genetic traits. Isolates of the same subtypes (cgMLST types (CTs) or MLST sequence types (STs)) were additionally compared by SNP calling. This revealed the occurrence of more than one CT within all three investigated plants and within two samples. Analysis of the fish set resulted in predominant assignment of isolates from pangasius catfish and salmon to ST2 and ST121, respectively, potentially indicating persistence within the respective production chains. The approach not only allowed the detection of distinct subtypes but also the determination of differences between closely related isolates, which need to be considered when interpreting WGS data for surveillance.

Pangenomics in Microbial and Crop Research: Progress, Applications, and Perspectives

Advances in sequencing technologies and bioinformatics tools have fueled a renewed interest in whole genome sequencing efforts in many organisms. The growing availability of multiple genome sequences has advanced our understanding of the within-species diversity, in the form of a pangenome. Pangenomics has opened new avenues for future research such as allowing dissection of complex molecular mechanisms and increased confidence in genome mapping. To comprehensively capture the genetic diversity for improving plant performance, the pangenome concept is further extended from species to genus level by the inclusion of wild species, constituting a super-pangenome. Characterization of pangenome has implications for both basic and applied research. The concept of pangenome has transformed the way biological questions are addressed. From understanding evolution and adaptation to elucidating host–pathogen interactions, finding novel genes or breeding targets to aid crop improvement to design effective vaccines for human prophylaxis, the increasing availability of the pangenome has revolutionized several aspects of biological research. The future availability of high-resolution pangenomes based on reference-level near-complete genome assemblies would greatly improve our ability to address complex biological problems.

International links between Streptococcus pneumoniae vaccine serotype 4 sequence type (ST) 801 in Northern European shipyard outbreaks of invasive pneumococcal disease

BACKGROUND

Pneumococcal disease outbreaks of vaccine preventable serotype 4 sequence type (ST)801 in shipyards have been reported in several countries. We aimed to use genomics to establish any international links between them.

METHODS

Sequence data from ST801-related outbreak isolates from Norway (n = 17), Finland (n = 11) and Northern Ireland (n = 2) were combined with invasive pneumococcal disease surveillance from the respective countries, and ST801-related genomes from an international collection (n = 41 of > 40,000), totalling 106 genomes. Raw data were mapped and recombination excluded before phylogenetic dating.

RESULTS

Outbreak isolates were relatively diverse, with up to 100 SNPs (single nucleotide polymorphisms) and a common ancestor estimated around the year 2000. However, 19 Norwegian and Finnish isolates were nearly indistinguishable (0-2 SNPs) with the common ancestor dated around 2017.

CONCLUSION

The total diversity of ST801 within the outbreaks could not be explained by recent transmission alone, suggesting that harsh environmental and associated living conditions reported in the shipyards may facilitate invasion of colonising pneumococci. However, near identical strains in the Norwegian and Finnish outbreaks does suggest that transmission between international shipyards also contributed to those outbreaks. This indicates the need for improved preventative measures in this working population including pneumococcal vaccination.

Diversity of carbapenem-resistant Acinetobacter baumannii and bacteriophage-mediated spread of the Oxa23 carbapenemase

Carbapenem-resistant Acinetobacter baumannii are prevalent in low- and middle-income countries such as Egypt, but little is known about the molecular epidemiology and mechanisms of resistance in these settings. Here, we characterize carbapenem-resistant A. baumannii from Alexandria, Egypt, and place it in a regional context. Fifty-four carbapenem-resistant isolates from Alexandria Main University Hospital (AMUH), Alexandria, Egypt, collected between 2010 and 2015 were genome sequenced using Illumina technology. Genomes were de novo assembled and annotated. Genomes for 36 isolates from the Middle East region were downloaded from GenBank. The core-gene compliment was determined using Roary, and analyses of recombination were performed in Gubbins. Multilocus sequence typing (MLST) sequence type (ST) and antibiotic-resistance genes were identified. The majority of Egyptian isolates belonged to one of three major clades, corresponding to Pasteur MLST clonal complex (CC^PAS) 1, CC^PAS2 and ST^PAS158. Strains belonging to ST^PAS158 have been reported almost exclusively from North Africa, the Middle East and Pakistan, and may represent a region-specific lineage. All isolates carried an oxa23 gene, six carried bla _NDM-1 and one carried bla _NDM-2. The oxa23 gene was located on a variety of different mobile elements, with Tn2006 predominant in CC^PAS2 strains, and Tn2008 predominant in other lineages. Of particular concern, in 8 of the 13 CC^PAS1 strains, the oxa23 gene was located in a temperate bacteriophage phiOXA, previously identified only once before in a CC^PAS1 clone from the USA military. The carbapenem-resistant A. baumannii population in AMUH is very diverse, and indicates an endemic circulating population, including a region-specific lineage. A major mechanism for oxa23 dissemination in CC^PAS1 isolates appears to be a bacteriophage, presenting new concerns about the ability of these carbapenemases to spread throughout the bacterial population.

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.

Population Genomics of emm4 Group A Streptococcus Reveals Progressive Replacement with a Hypervirulent Clone in North America

Clonal replacement is a major driver for changes in bacterial disease epidemiology. Recently, it has been proposed that episodic emergence of novel, hypervirulent clones of group A Streptococcus (GAS) results from acquisition of a 36-kb DNA region leading to increased expression of the cytotoxins Nga (NADase) and SLO (streptolysin O). We previously described a gene fusion event involving the gene encoding the GAS M protein (emm) and an adjacent M-like protein (enn) in the emm4 GAS population, a GAS emm type that lacks the hyaluronic acid capsule. Using whole-genome sequencing of a temporally and geographically diverse set of 1,126 isolates, we discovered that the North American emm4 GAS population has undergone clonal replacement with emergent GAS strains completely replacing historical isolates by 2017. Emergent emm4 GAS strains contained a handful of small genetic variations, including the emm-enn gene fusion, and showed a marked in vitro growth defect compared to historical strains. In contrast to other previously described GAS clonal replacement events, emergent emm4 GAS strains were not defined by acquisition of exogenous DNA and had no significant increase in transcript levels of nga and slo toxin genes via RNA sequencing and quantitative real-time PCR analysis relative to historic strains. Despite the in vitro growth differences, emergent emm4 GAS strains were hypervirulent in mice and ex vivo growth in human blood compared to historical strains. Thus, these data detail the emergence and dissemination of a hypervirulent acapsular GAS clone defined by small, endogenous genetic variation, thereby defining a novel model for GAS strain replacement.

IMPORTANCE Severe invasive infections caused by group A Streptococcus (GAS) result in substantial morbidity and mortality in children and adults worldwide. Previously, GAS clonal strain replacement has been attributed to acquisition of exogenous DNA leading to novel virulence gene acquisition or increased virulence gene expression. Our study of type emm4 GAS identified emergence of a hypervirulent GAS clade defined by variation in endogenous DNA content and lacking augmented toxin gene expression relative to replaced strains. These findings expand our understanding of the molecular mechanisms underlying bacterial clonal emergence.

Phylogeography and resistome of pneumococcal meningitis in West Africa before and after vaccine introduction

Despite contributing to the large disease burden in West Africa, little is known about the genomic epidemiology of Streptococcus pneumoniae which cause meningitis among children under 5 years old in the region. We analysed whole-genome sequencing data from 185 S. pneumoniae isolates recovered from suspected paediatric meningitis cases as part of the World Health Organization (WHO) invasive bacterial diseases surveillance from 2010 to 2016. The phylogeny was reconstructed, accessory genome similarity was computed and antimicrobial-resistance patterns were inferred from the genome data and compared to phenotypic resistance from disc diffusion. We studied the changes in the distribution of serotypes pre- and post-pneumococcal conjugate vaccine (PCV) introduction in the Central and Western sub-regions separately. The overall distribution of non-vaccine, PCV7 (4, 6B, 9V, 14, 18C, 19F and 23F) and additional PCV13 serotypes (1, 3, 5, 6A, 19A and 7F) did not change significantly before and after PCV introduction in the Central region (Fisher's test P value 0.27) despite an increase in the proportion of non-vaccine serotypes to 40 % (n=6) in the post-PCV introduction period compared to 21.9 % (n=14). In the Western sub-region, PCV13 serotypes were more dominant among isolates from The Gambia following the introduction of PCV7, 81 % (n=17), compared to the pre-PCV period in neighbouring Senegal, 51 % (n=27). The phylogeny illustrated the diversity of strains associated with paediatric meningitis in West Africa and highlighted the existence of phylogeographical clustering, with isolates from the same sub-region clustering and sharing similar accessory genome content. Antibiotic-resistance genotypes known to confer resistance to penicillin, chloramphenicol, co-trimoxazole and tetracycline were detected across all sub-regions. However, there was no discernible trend linking the presence of resistance genotypes with the vaccine introduction period or whether the strain was a vaccine or non-vaccine serotype. Resistance genotypes appeared to be conserved within selected sub-clades of the phylogenetic tree, suggesting clonal inheritance. Our data underscore the need for continued surveillance on the emergence of non-vaccine serotypes as well as chloramphenicol and penicillin resistance, as these antibiotics are likely still being used for empirical treatment in low-resource settings. This article contains data hosted by Microreact.

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.

PanACoTA: a modular tool for massive microbial comparative genomics

The study of the gene repertoires of microbial species, their pangenomes, has become a key part of microbial evolution and functional genomics. Yet, the increasing number of genomes available complicates the establishment of the basic building blocks of comparative genomics. Here, we present PanACoTA (https://github.com/gem-pasteur/PanACoTA), a tool that allows to download all genomes of a species, build a database with those passing quality and redundancy controls, uniformly annotate and then build their pangenome, several variants of core genomes, their alignments and a rapid but accurate phylogenetic tree. While many programs building pangenomes have become available in the last few years, we have focused on a modular method, that tackles all the key steps of the process, from download to phylogenetic inference. While all steps are integrated, they can also be run separately and multiple times to allow rapid and extensive exploration of the parameters of interest. PanACoTA is built in Python3, includes a singularity container and features to facilitate its future development. We believe PanACoTa is an interesting addition to the current set of comparative genomics tools, since it will accelerate and standardize the more routine parts of the work, allowing microbial genomicists to more quickly tackle their specific questions.

Diversification of Colonization Factors in a Multidrug-Resistant Escherichia coli Lineage Evolving under Negative Frequency-Dependent Selection

Escherichia coli is a major cause of bloodstream and urinary tract infections globally. The wide dissemination of multidrug-resistant (MDR) strains of extraintestinal pathogenic E. coli (ExPEC) poses a rapidly increasing public health burden due to narrowed treatment options and increased risk of failure to clear an infection. Here, we present a detailed population genomic analysis of the ExPEC ST131 clone, in which we seek explanations for its success as an emerging pathogenic strain beyond the acquisition of antimicrobial resistance (AMR) genes. We show evidence for evolution toward separate ecological niches for the main clades of ST131 and differential evolution of anaerobic metabolism, key colonization, and virulence factors. We further demonstrate that negative frequency-dependent selection acting across accessory loci is a major mechanism that has shaped the population evolution of this pathogen.IMPORTANCE Infections with multidrug-resistant (MDR) strains of Escherichia coli are a significant global public health concern. To combat these pathogens, we need a deeper understanding of how they evolved from their background populations. By understanding the processes that underpin their emergence, we can design new strategies to limit evolution of new clones and combat existing clones. By combining population genomics with modelling approaches, we show that dominant MDR clones of E. coli are under the influence of negative frequency-dependent selection, preventing them from rising to fixation in a population. Furthermore, we show that this selection acts on genes involved in anaerobic metabolism, suggesting that this key trait, and the ability to colonize human intestinal tracts, is a key step in the evolution of MDR clones of E. coli.

Fast and flexible bacterial genomic epidemiology with PopPUNK

The routine use of genomics for disease surveillance provides the opportunity for high-resolution bacterial epidemiology. Current whole-genome clustering and multilocus typing approaches do not fully exploit core and accessory genomic variation, and they cannot both automatically identify, and subsequently expand, clusters of significantly similar isolates in large data sets spanning entire species. Here, we describe PopPUNK (Population Partitioning Using Nucleotide K -mers), a software implementing scalable and expandable annotation- and alignment-free methods for population analysis and clustering. Variable-length k-mer comparisons are used to distinguish isolates' divergence in shared sequence and gene content, which we demonstrate to be accurate over multiple orders of magnitude using data from both simulations and genomic collections representing 10 taxonomically widespread species. Connections between closely related isolates of the same strain are robustly identified, despite interspecies variation in the pairwise distance distributions that reflects species' diverse evolutionary patterns. PopPUNK can process 103-104 genomes in a single batch, with minimal memory use and runtimes up to 200-fold faster than existing model-based methods. Clusters of strains remain consistent as new batches of genomes are added, which is achieved without needing to reanalyze all genomes de novo. This facilitates real-time surveillance with consistent cluster naming between studies and allows for outbreak detection using hundreds of genomes in minutes. Interactive visualization and online publication is streamlined through the automatic output of results to multiple platforms. PopPUNK has been designed as a flexible platform that addresses important issues with currently used whole-genome clustering and typing methods, and has potential uses across bacterial genetics and public health research.