Search

Neisseria gonorrhoeae, the bacterium responsible for the sexually transmitted disease gonorrhoea, has shown an extraordinary ability to develop antimicrobial resistance (AMR) to multiple classes of antimicrobials. With no available vaccine, managingN. gonorrhoeaeinfections demands effective preventive measures, antibiotic treatment and epidemiological surveillance. The latter two are progressively being supported by the generation of whole-genome sequencing (WGS) data on behalf of national and international surveillance programmes. In this context, this study aims to performN. gonorrhoeaeclustering into genogroups based on WGS data, for enhanced prospective laboratory surveillance. Particularly, it aims to identify the major circulating WGS-genogroups in Europe and to establish a relationship between these and AMR. Ultimately, it enriches public databases by contributing with WGS data from Portuguese isolates spanning 15 years of surveillance. A total of 3791 carefully inspectedN. gonorrhoeaegenomes from isolates collected across Europe were analysed using a gene-by-gene approach (i.e. using cgMLST). Analysis of cluster composition and stability allowed the classification of isolates into a two-step hierarchical genogroup level determined by two allelic distance thresholds revealing cluster stability. Genogroup clustering in general agreed with availableN. gonorrhoeaetyping methods [i.e. MLST (multilocus sequence typing), NG-MAST (N. gonorrhoeaemulti-antigen sequence typing) and PubMLST core-genome groups], highlighting the predominant genogroups circulating in Europe, and revealed that the vast majority of the genogroups present a dominant AMR profile. Additionally, a non-static gene-by-gene approach combined with a more discriminatory threshold for potential epidemiological linkage enabled us to match data with previous reports on outbreaks or transmission chains. In conclusion, this genogroup assignment allows a comprehensive analysis ofN. gonorrhoeaegenetic diversity and the identification of the WGS-based genogroups circulating in Europe, while facilitating the assessment (and continuous monitoring) of their frequency, geographical dispersion and potential association with specific AMR signatures. This strategy may benefit public-health actions through the prioritization of genogroups to be controlled, the identification of emerging resistance carriage, and the potential facilitation of data sharing and communication.

Actinobacteriaare a large and diverse phylum of bacteria that contains medically and ecologically relevant organisms. Many members are valuable sources of bioactive natural products and chemical precursors that are exploited in the clinic. These are made using the enzyme pathways encoded in their complex genomes. Whilst the number of sequenced genomes has increased rapidly in the last twenty years, the large size and complexity of many Actinobacterial genomes means that the sequences remain incomplete and consist of large numbers of contigs with poor annotation, which hinders large scale comparative genomics and evolutionary studies. To enable greater understanding and exploitation of Actinobacterial genomes, specialist genomic databases must be linked to high-quality genome sequences. Here we provide a curated database of 612 high-quality actinobacterial genomes from 80 genera, chosen to represent a broad phylogenetic group with equivalent genome reannotation. Utilising this database will provide researchers with a framework for evolutionary and metabolic studies, to enable a foundation for genome and metabolic engineering, to facilitate discovery of novel bioactive therapeutics and studies on gene family evolution.Significance as a bioresource to the communityTheActinobacteriaare a large diverse phylum of bacteria, often with large, complex genomes with a high G+C content. Sequence databases have great variation in the quality of sequences, equivalence of annotation and phylogenetic representation, which makes it challenging to undertake evolutionary and phylogenetic studies. To address this, we have assembled a curated, taxa-specific, non-redundant database to aid detailed comparative analysis of Actinobacteria. ActDES constitutes a novel resource for the community of Actinobacterial researchers that will be useful primarily for two types of analyses: (i) comparative genomic studies – facilitated by reliable identification of orthologs across a set of defined, phylogenetically-representative genomes, and (ii) phylogenomic studies which will be improved by identification of gene subsets at specified taxonomic level. These analyses can then act as a springboard for the studies of the evolution of virulence genes, the evolution of metabolism and identification of targets for metabolic engineering.Data summaryAll genome sequences used in this study can be found in the NCBI taxonomy browserhttps://www.ncbi.nlm.nih.gov/Taxonomy/Browser/www.tax.cgiand are summarised along with Accession numbers in Table S1All other data are available on Figsharehttps://doi.org/10.6084/m9.figshare.12167529andhttps://doi.org/10.5281/zenodo.3830391Perl script files available on GitHubhttps://github.com/nselem/ActDESincluding details of how to batch annotate genomes in RAST from the terminalhttps://github.com/nselem/myrastSupp. Table S1List of genomes from NCBI (Actinobacteria database.xlsx)https://doi.org/10.6084/m9.figshare.12167529CVS genome annotation files including the FASTA files of nucleotide and amino acids sequences (individual .cvs files)https://doi.org/10.6084/m9.figshare.12167880BLAST nucleotide database (.fasta file)https://doi.org/10.6084/m9.figshare.12167724BLAST protein database (.fasta file)https://doi.org/10.6084/m9.figshare.12167724Supp. Table S2 Expansion table genus level (Expansion table.xlsx Tab Genus level)https://doi.org/10.6084/m9.figshare.12167529Supp. Table S2 Expansion table species level (Expansion table.xlsx Tab species level)https://doi.org/10.6084/m9.figshare.12167529All GlcP and Glk data – blast hits from ActDES database, MUSCLE Alignment files and .nwk tree files can be found athttps://doi.org/10.6084/m9.figshare.12167529Interactive trees in Microreact for Glk treehttps://microreact.org/project/w_KDfn1xA/90e6759eand associated files can be found athttps://doi.org/10.6084/m9.figshare.12326441.v1Interactive trees in Microreact for GlcP treehttps://microreact.org/project/VBUdiQ5_k/0fc4622band associated files can be found athttps://doi.org/10.6084/m9.figshare.12326441.v1

BackgroundCarbapenem-resistant Klebsiella pneumoniae (CRKP) remains a major clinical pathogen and public health threat with few therapeutic options. The mobilome, resistome, methylome, virulome and phylogeography of CRKP were characterised.MethodsCRKP collected in 2018 were subjected to antimicrobial susceptibility testing, screening by multiplex-PCR, genotyping by Repetitive Element Palindromic-Polymerase Chain Reaction (REP-PCR), plasmid size, number, incompatibility, and mobility analyses, and PacBio’s SMRT sequencing (n=6).Results & conclusionThere were 56 multidrug-resistant CRKP, having blaOXA-48-like and blaNDM-1/7 carbapenemases on self-transmissible IncF, A/C, IncL/M and IncX3 plasmids endowed with prophages, traT, resistance islands and type I and II restriction modification systems (RMS). These plasmids were of close evolutionary relationship to several plasmids globally whilst the strains also clustered with several global clades, evincing transboundary horizontal and vertical dissemination. Reduced susceptibility to colistin occurred in 23 strains. Common clones included ST307, ST607, ST17, ST39, and ST3559. IncFIIk virulent plasmid replicon was present in 56 strains. The six strains contained at least 41 virulence genes and four different K- and O-loci types: KL2, KL25, KL27, KL102, O1, O2, O4 and O5. Types I, II, and III RMS, conferring m6A (GATC, GATGNNNNNNTTG, CAANNNNNNCATC motifs) and m4C (CCWGG) modifications on chromosomes and plasmids, were found.There is plasmid-mediated, clonal, and multiclonal dissemination of blaOXA-48-like and blaNDM-1 in South Africa, mirroring international epidemiology of similar clones and plasmids. Plasmid-mediated transmission of RMS, virulome and prophages influence bacterial evolution, epidemiology, pathogenicity, and resistance, threatening infection treatment. RMS influence on antimicrobial and bacteriophage therapy needs urgent investigation.Highlights/ImportanceK. pneumoniae is a major pathogen implicated in numerous nosocomial infections. Worryingly, we show that K. pneumoniae isolates from South Africa, Africa and globally are endowed with rich resistomes and mobilomes that make them almost pandrug resistant. The isolates in this study contained rich virulomes and prophages on both chromosomes and plasmids, with close evolutionary kith or kin to other plasmids identified worldwide. There was a rich diversity of restriction modification systems that regulate virulence, transcription, and plasmid mobility in bacteria, facilitating the epidemiology, resistance, pathogenicity and genomic evolution of the strains, and threatening antimicrobial and bacteriophage therapy.

Mucormycoses are invasive infections by Rhizopus species and other Mucorales. Over 10 months, four solid organ transplant (SOT) recipients at our centre developed mucormycosis due to Rhizopus microsporus (n=2), R. arrhizus (n=1) or Lichtheimia corymbifera (n=1), at a median 31.5 days (range: 13–34) post-admission. We performed whole genome sequencing (WGS) on 72 Mucorales isolates (45 R. arrhizus, 19 R. delemar, six R. microsporus, two Lichtheimia species) from these patients, from five patients with community-acquired mucormycosis, and from hospital and regional environments. Isolates were compared by core protein phylogeny and global genomic features, including genome size, guanine–cytosine percentages, shared protein families and paralogue expansions. Patient isolates fell into six core phylogenetic lineages (clades). Phylogenetic and genomic similarities of R. microsporus isolates recovered 7 months apart from two SOT recipients in adjoining hospitals suggested a potential common source exposure. However, isolates from other patients and environmental sites had unique genomes. Many isolates that were indistinguishable by core phylogeny were distinct by one or more global genomic comparisons. Certain clades were recovered throughout the study period, whereas others were found at particular time points. In conclusion, mucormycosis cases could not be genetically linked to a definitive environmental source. Comprehensive genomic analyses eliminated false associations between Mucorales isolates that would have been assigned using core phylogenetic or less extensive genomic comparisons. The genomic diversity of Mucorales mandates that multiple isolates from individual patients and environmental sites undergo WGS during epidemiological investigations. However, exhaustive surveillance of fungal populations in a hospital and surrounding community is probably infeasible.

Sellimonas intestinalis is a Gram-positive and anaerobic bacterial species previously considered as uncultivable. Although little is known about this Lachnospiraceae family member, its increased abundance has been reported in patients who have recovered from intestinal homeostasis after dysbiosis events. In this context, the aim of the present study was to take advantage of a massive in vitro culture protocol that allowed the recovery of extremely oxygen-sensitive species from faecal samples, which led to isolation of S. intestinalis . Whole genome analyses of 11 S . intestinalis genomes revealed that this species has a highly conserved genome with 99.7 % 16S rRNA gene sequence similarity, average nucleotide polymorphism results >95, and 50.1 % of its coding potential being part of the core genome. Despite this, the variable portion of its genome was informative enough to reveal the existence of three lineages (lineage-I including isolates from Chile and France, lineage-II from South Korea and Finland, and lineage-III from China and one isolate from the USA) and evidence of some recombination signals. The identification of a cluster of orthologous groups revealed a high number of genes involved in metabolism, including amino acid and carbohydrate transport as well as energy production and conversion, which matches with the metabolic profile previously reported for microbiota from healthy individuals. Additionally, virulence factors and antimicrobial resistance genes were found (mainly in lineage-III), which could favour their survival during antibiotic-induced dysbiosis. These findings provide the basis of knowledge about the potential of S. intestinalis as a bioindicator of intestinal homeostasis recovery and contribute to advancing the characterization of gut microbiota members with beneficial potential.

Urinary tract infections (UTIs) are the most common bacterial infections requiring medical attention and a leading justification for antibiotic prescription. Trimethoprim is prescribed empirically for uncomplicated cases. UTIs are primarily caused by extraintestinal pathogenic Escherichia coli (ExPEC) and ExPEC strains play a central role in disseminating antimicrobial-resistance genes worldwide. Here, we describe the whole-genome sequences of trimethoprim-resistant ExPEC and/or ExPEC from recurrent UTIs (67 in total) from patients attending a regional Australian hospital from 2006 to 2008. Twenty-three sequence types (STs) were observed, with ST131 predominating (28 %), then ST69 and ST73 (both 7 %). Co-occurrence of trimethoprim-resistance genes with genes conferring resistance to extended-spectrum β-lactams, heavy metals and quaternary ammonium ions was a feature of the ExPEC described here. Seven trimethoprim-resistance genes were identified, most commonly dfrA17 (38 %) and dfrA12 (18 %). An uncommon dfrB4 variant was also observed. Two blaCTX-M variants were identified – blaCTX-M-15 (16 %) and blaCTX-M-14 (10 %). The former was always associated with dfrA12, the latter with dfrA17, and all blaCTX-M genes co-occurred with chromate-resistance gene chrA. Eighteen class 1 integron structures were characterized, and chrA featured in eight structures; dfrA genes featured in seventeen. ST131 H30Rx isolates possessed distinct antimicrobial gene profiles comprising aac(3)-IIa, aac(6)-Ib-cr, aph(3′)-Ia, aadA2, blaCTX-M-15, blaOXA-1 and dfrA12. The most common virulence-associated genes (VAGs) were fimH, fyuA, irp2 and sitA (all 91 %). Virulence profile clustering showed ST131 H30 isolates carried similar VAGs to ST73, ST405, ST550 and ST1193 isolates. The sole ST131 H27 isolate carried molecular predictors of enteroaggregative E. coli /ExPEC hybrid strains (aatA, aggR, fyuA). Seven isolates (10 %) carried VAGs suggesting ColV plasmid carriage. Finally, SNP analysis of serial UTI patients experiencing worsening sequelae demonstrated a high proportion of point mutations in virulence factors.

Epidemic and pandemic clones of bacterial pathogens with distinct characteristics continually emerge, replacing those previously dominant through mechanisms that remain poorly characterized. Here, whole-genome-sequencing-powered epidemiology linked horizontal transfer of a virulence gene, sopE, to the emergence and clonal expansion of a new epidemic Salmonella enterica serovar Typhimurium (S. Typhimurium) clone. The sopE gene is sporadically distributed within the genus Salmonella and rare in S . enterica Typhimurium lineages, but was acquired multiple times during clonal expansion of the currently dominant pandemic monophasic S. Typhimurium sequence type (ST) 34 clone. Ancestral state reconstruction and time-scaled phylogenetic analysis indicated that sopE was not present in the common ancestor of the epidemic clade, but later acquisition resulted in increased clonal expansion of sopE-containing clones that was temporally associated with emergence of the epidemic, consistent with increased fitness. The sopE gene was mainly associated with a temperate bacteriophage mTmV, but recombination with other bacteriophage and apparent horizontal gene transfer of the sopE gene cassette resulted in distribution among at least four mobile genetic elements within the monophasic S . enterica Typhimurium ST34 epidemic clade. The mTmV prophage lysogenic transfer to other S. enterica serovars in vitro was limited, but included the common pig-associated S . enterica Derby (S. Derby). This may explain mTmV in S. Derby co-circulating on farms with monophasic S. Typhimurium ST34, highlighting the potential for further transfer of the sopE virulence gene in nature. We conclude that whole-genome epidemiology pinpoints potential drivers of evolutionary and epidemiological dynamics during pathogen emergence, and identifies targets for subsequent research in epidemiology and bacterial pathogenesis.

Ovine and caprine brucellosis, caused by Brucella melitensis , is one of the world’s most widespread zoonoses and is a major cause of economic losses in domestic ruminant production. In Italy, the disease remains endemic in several southern provinces, despite an ongoing brucellosis eradication programme. In this study, we used whole-genome sequencing to detail the genetic diversity of circulating strains, and to examine the origins of the predominant sub-lineages of B. melitensis in Italy. We reconstructed a global phylogeny of B. melitensis , strengthened by 339 new whole-genome sequences, from Italian isolates collected from 2011 to 2018 as part of a national livestock surveillance programme. All Italian strains belonged to the West Mediterranean lineage, which further divided into two major clades that diverged roughly between the 5th and 7th centuries. We observed that Sicily serves as a brucellosis burden hotspot, giving rise to several distinct sub-lineages. More than 20 putative outbreak clusters of ovine and caprine brucellosis were identified, several of which persisted over the 8 year survey period despite an aggressive brucellosis eradication campaign. While the outbreaks in Central and Northern Italy were generally associated with introductions of single clones of B. melitensis and their subsequent dissemination within neighbouring territories, we observed weak geographical segregation of genotypes in the southern regions. Biovar determination, recommended in routine analysis of all Brucella strains by the World Organisation for Animal Health (OIE), could not discriminate among the four main global clades. This demonstrates a need for updating the guidelines used for monitoring B. melitensis transmission and spread, both at the national and international level, and to include whole-genome-based typing as the principal method for identification and tracing of brucellosis outbreaks.

Although several studies have investigated genetic diversity of Leishmania infantum in North Africa, genome-wide analyses are lacking. Here, we conducted comparative analyses of nuclear and mitochondrial genomes of seven L. infantum isolates from Tunisia with the aim to gain insight into factors that drive genomic and phenotypic adaptation. Isolates were from cured (n=4) and recurrent (n=3) visceral leishmaniasis (VL) cases, originating from northern (n=2) and central (n=5) Tunisia, where respectively stable and emerging VL foci are observed. All isolates from relapsed patients were from Kairouan governorate (Centre); one showing resistance to the anti-leishmanial drug Meglumine antimoniate. Nuclear genome diversity of the isolates was analysed by comparison to the L. infantum JPCM5 reference genome. Kinetoplast maxi and minicircle sequences (1 and 59, respectively) were extracted from unmapped reads and identified by blast analysis against public data sets. The genome variation analysis grouped together isolates from the same geographical origins. Strains from the North were very different from the reference showing more than 34 587 specific single nucleotide variants, with one isolate representing a full genetic hybrid as judged by variant frequency. Composition of minicircle classes within isolates corroborated this geographical population structure. Read depth analysis revealed several significant gene copy number variations correlating with either geographical origin (amastin and Hsp33 genes) or relapse (CLN3 gene). However, no specific gene copy number variation was found in the drug-resistant isolate. In contrast, resistance was associated with a specific minicircle pattern suggesting Leishmania mitochondrial DNA as a potential novel source for biomarker discovery.

Klebsiella pneumoniae (Kp) reference strain Kp52.145 is widely used in experimental Klebsiella pathophysiology. Since 1935, only one other strain of the same sublineage (sequence type ST66, capsular serotype K2) was isolated (AJ210, Australia). Here, we describe a community-acquired invasive infection caused by a ST66-K2 Kp strain in France. Four hypermucoviscous Kp isolates responsible for acute otitis media, meningitis, bacteraemia and bacteriuria, respectively, were obtained from a patient with a history of chronic alcoholism and diabetes mellitus, and infected with HIV. The isolates were characterized by phenotypic and genomic methods. The four genetically identical ST66-K2 isolates presented a full antimicrobial susceptibility profile, including to ampicillin, corresponding to a single strain (SB5881), which was more closely related to AJ210 (135 SNPs) than to Kp52.145 (388 SNPs). Colibactin and yersiniabactin gene clusters were present on the integrative and conjugative element ICEKp10 in the chromosome. The two plasmids from Kp52.145 were detected in SB5881. In addition to carrying genes for virulence factors RmpA, aerobactin and salmochelin, plasmid II has acquired in SB5881, the conjugation machinery gene cluster from plasmid I. We report the first case of community-acquired infection caused by a hypervirulent ST66-K2 Kp strain in Europe. This demonstrates the long-term persistence of the high-virulence and laboratory model ST66-K2 sublineage. The combination of a conjugative apparatus and major virulence genes on a single plasmid may contribute to the co-occurrence of hypervirulence and multidrug resistance in single Kp strains.

Clostridioides difficile is the primary infectious cause of antibiotic-associated diarrhea. Local transmissions and international outbreaks of this pathogen have been previously elucidated by bacterial whole-genome sequencing, but comparative genomic analyses at the global scale were hampered by the lack of specific bioinformatic tools. Here we introduce a publicly accessible database within EnteroBase (http://enterobase.warwick.ac.uk) that automatically retrieves and assembles C. difficile short-reads from the public domain, and calls alleles for core-genome multilocus sequence typing (cgMLST). We demonstrate that comparable levels of resolution and precision are attained by EnteroBase cgMLST and single-nucleotide polymorphism analysis. EnteroBase currently contains 18 254 quality-controlled C. difficile genomes, which have been assigned to hierarchical sets of single-linkage clusters by cgMLST distances. This hierarchical clustering is used to identify and name populations of C. difficile at all epidemiological levels, from recent transmission chains through to epidemic and endemic strains. Moreover, it puts newly collected isolates into phylogenetic and epidemiological context by identifying related strains among all previously published genome data. For example, HC2 clusters (i.e. chains of genomes with pairwise distances of up to two cgMLST alleles) were statistically associated with specific hospitals (P

Groundbreaking studies conducted in the mid-1980s demonstrated the possibility of sequencing ancient DNA (aDNA), which has allowed us to answer fundamental questions about the human past. Microbiologists were thus given a powerful tool to glimpse directly into inscrutable bacterial history, hitherto inaccessible due to a poor fossil record. Initially plagued by concerns regarding contamination, the field has grown alongside technical progress, with the advent of high-throughput sequencing being a breakthrough in sequence output and authentication. Albeit burdened with challenges unique to the analysis of bacteria, a growing number of viable sources for aDNA has opened multiple avenues of microbial research. Ancient pathogens have been extracted from bones, dental pulp, mummies and historical medical specimens and have answered focal historical questions such as identifying the aetiological agent of the black death as Yersinia pestis . Furthermore, ancient human microbiomes from fossilized faeces, mummies and dental plaque have shown shifts in human commensals through the Neolithic demographic transition and industrial revolution, whereas environmental isolates stemming from permafrost samples have revealed signs of ancient antimicrobial resistance. Culminating in an ever-growing repertoire of ancient genomes, the quickly expanding body of bacterial aDNA studies has also enabled comparisons of ancient genomes to their extant counterparts, illuminating the evolutionary history of bacteria. In this review we summarize the present avenues of research and contextualize them in the past of the field whilst also pointing towards questions still to be answered.

Verotoxigenic Escherichia coli (VTEC) are food- and water-borne pathogens associated with both sporadic illness and outbreaks of enteric disease. While it is known that cattle are reservoirs of VTEC, little is known about the genomic variation of VTEC in cattle, and whether the variation in genomes reported for human outbreak strains is consistent with individual animal or group/herd sources of infection. A previous study of VTEC prevalence identified serotypes carried persistently by three consecutive cohorts of heifers within a closed herd of cattle. This present study aimed to: (i) determine whether the genomic relatedness of bovine isolates is similar to that reported for human strains associated with single source outbreaks, (ii) estimate the rates of genome change among dominant serotypes over time within a cattle herd, and (iii) identify genomic features of serotypes associated with persistence in cattle. Illumina MiSeq genome sequencing and genotyping based on allelic and single nucleotide variations were completed, while genome change over time was measured using Bayesian evolutionary analysis sampling trees. The accessory genome, including the non-protein-encoding intergenic regions (IGRs), virulence factors, antimicrobial-resistance genes and plasmid gene content of representative persistent and sporadic cattle strains were compared using Fisher’s exact test corrected for multiple comparisons. Herd strains from serotypes O6:H34 (n=22), O22:H8 (n=30), O108:H8 (n=39), O139:H19 (n=44) and O157:H7 (n=106) were readily distinguishable from epidemiologically unrelated strains of the same serotype using a similarity threshold of 10 or fewer allele differences between adjacent nodes. Temporal-cohort clustering within each serotype was supported by date randomization analysis. Substitutions per site per year were consistent with previously reported values for E. coli ; however, there was low branch support for these values. Acquisition of the phage-encoded Shiga toxin 2 gene in serotype O22:H8 was observed. Pan-genome analyses identified accessory regions that were more prevalent in persistent serotypes (P≤0.05) than in sporadic serotypes. These results suggest that VTEC serotypes from a specific cattle population are highly clonal with a similar level of relatedness as human single-source outbreak-associated strains, but changes in the genome occur gradually over time. Additionally, elements in the accessory genomes may provide a selective advantage for persistence of VTEC within cattle herds.

Mobile genetic elements (MGEs) are agents of bacterial evolution and adaptation. Genome sequencing provides an unbiased approach that has revealed an abundance of MGEs in prokaryotes, mainly plasmids and integrative conjugative elements. Nevertheless, many mobilomes, particularly those from environmental bacteria, remain underexplored despite their representing a reservoir of genes that can later emerge in the clinic. Here, we explored the mobilome of the Mycobacteriaceae family, focusing on strains from Brazilian Atlantic Forest soil. Novel Mycolicibacterium and Mycobacteroides strains were identified, with the former ones harbouring linear and circular plasmids encoding the specialized type-VII secretion system (T7SS) and mobility-associated genes. In addition, we also identified a T4SS-mediated integrative conjugative element (ICEMyc226) encoding two T7SSs and a number of xenobiotic degrading genes. Our study uncovers the diversity of the Mycobacteriaceae mobilome, providing the evidence of an ICE in this bacterial family. Moreover, the presence of T7SS genes in an ICE, as well as plasmids, highlights the role of these mobile genetic elements in the dispersion of T7SS.

Klebsiella pneumoniae is the clinically most important species within the genus Klebsiella and, as a result of the continuous emergence of multi-drug resistant (MDR) strains, the cause of severe nosocomial infections. The decline in the effectiveness of antibiotic treatments for infections caused by MDR bacteria has generated particular interest in the study of bacteriophages. In this study, we characterized a total of 40 temperate bacteriophages (prophages) with a genome range of 11.454–84.199 kb, predicted from 16 carbapenemase-producing clinical strains of K. pneumoniae belonging to different sequence types, previously identified by multilocus sequence typing. These prophages were grouped into the three families in the order Caudovirales (27 prophages belonging to the family Myoviridae, 10 prophages belonging to the family Siphoviridae and 3 prophages belonging to the family Podoviridae). Genomic comparison of the 40 prophage genomes led to the identification of four prophages isolated from different strains and of genome sizes of around 33.3, 36.1, 39.6 and 42.6 kb. These prophages showed sequence similarities (query cover >90 %, identity >99.9 %) with international Microbe Versus Phage (MVP) (http://mvp.medgenius.info/home) clusters 4762, 4901, 3499 and 4280, respectively. Phylogenetic analysis revealed the evolutionary proximity among the members of the four groups of the most frequently identified prophages in the bacterial genomes studied (33.3, 36.1, 39.6 and 42.6 kb), with bootstrap values of 100 %. This allowed the prophages to be classified into three clusters: A, B and C. Interestingly, these temperate bacteriophages did not infect the highest number of strains as indicated by a host-range assay, these results could be explained by the development of superinfection exclusion mechanisms. In addition, bioinformatic analysis of the 40 identified prophages revealed the presence of 2363 proteins. In total, 59.7 % of the proteins identified had a predicted function, mainly involving viral structure, transcription, replication and regulation (lysogenic/lysis). Interestingly, some proteins had putative functions associated with bacterial virulence (toxin expression and efflux pump regulators), phage defence profiles such as toxin–antitoxin modules, an anti-CRISPR/Cas9 protein, TerB protein (from terZABCDE operon) and methyltransferase proteins.

This study sought to assess the genetic variability of Escherichia coli isolated from bloodstream infections (BSIs) presenting at Concord Hospital, Sydney during 2013-2016. Whole-genome sequencing was used to characterize 81 E. coli isolates sourced from community-onset (CO) and hospital-onset (HO) BSIs. The cohort comprised 64 CO and 17 HO isolates, including 35 multidrug-resistant (MDR) isolates exhibiting phenotypic resistance to three or more antibiotic classes. Phylogenetic analysis identified two major ancestral clades. One was genetically diverse with 25 isolates distributed in 16 different sequence types (STs) representing phylogroups A, B1, B2, C and F, while the other comprised phylogroup B2 isolates in subclades representing the ST131, ST73 and ST95 lineages. Forty-seven isolates contained a class 1 integron, of which 14 carried bla CTX -M-gene. Isolates with a class 1 integron carried more antibiotic resistance genes than isolates without an integron and, in most instances, resistance genes were localized within complex resistance loci (CRL). Resistance to fluoroquinolones could be attributed to point mutations in chromosomal parC and gyrB genes and, in addition, two isolates carried a plasmid-associated qnrB4 gene. Co-resistance to fluoroquinolone and broad-spectrum beta-lactam antibiotics was associated with ST131 (HO and CO), ST38 (HO), ST393 (CO), ST2003 (CO) and ST8196 (CO and HO), a novel ST identified in this study. Notably, 10/81 (12.3 %) isolates with ST95 (5 isolates), ST131 (2 isolates), ST88 (2 isolates) and a ST540 likely carry IncFII-IncFIB plasmid replicons with a full spectrum of virulence genes consistent with the carriage of ColV-like plasmids. Our data indicate that IncF plasmids play an important role in shaping virulence and resistance gene carriage in BSI E. coli in Australia.

Knowledge of pneumococcal lineages, their geographic distribution and antibiotic resistance patterns, can give insights into global pneumococcal disease. We provide interactive bioinformatic outputs to explore such topics, aiming to increase dissemi-nation of genomic insights to the wider community, without the need for specialist training. We prepared 12 country-specific phylogenetic snapshots, and international phylogenetic snapshots of 73 common Global Pneumococcal Sequence Clusters (GPSCs) previously defined using PopPUNK, and present them in Microreact. Gene presence and absence defined using Roary, and recombination profiles derived from Gubbins are presented in Phandango for each GPSC. Temporal phylogenetic signal was assessed for each GPSC using BactDating. We provide examples of how such resources can be used. In our example use of a country-specific phylogenetic snapshot we determined that serotype 14 was observed in nine unrelated genetic backgrounds in South Africa. The international phylogenetic snapshot of GPSC9, in which most serotype 14 isolates from South Africa were observed, highlights that there were three independent sub-clusters represented by South African serotype 14 isolates. We estimated from the GPSC9-dated tree that the sub-clusters were each established in South Africa during the 1980s. We show how recombination plots allowed the identification of a 20 kb recombination spanning the capsular polysaccharide locus within GPSC97. This was consistent with a switch from serotype 6A to 19A estimated to have occured in the 1990s from the GPSC97-dated tree. Plots of gene presence/absence of resistance genes (tet, erm, cat) across the GPSC23 phylogeny were consistent with acquisition of a composite transposon. We estimated from the GPSC23-dated tree that the acquisition occurred between 1953 and 1975. Finally, we demonstrate the assignment of GPSC31 to 17 externally generated pneumococcal serotype 1 assemblies from Utah via Pathogenwatch. Most of the Utah isolates clustered within GPSC31 in a USA-specific clade with the most recent common ancestor estimated between 1958 and 1981. The resources we have provided can be used to explore to data, test hypothesis and generate new hypotheses. The accessible assignment of GPSCs allows others to contextualize their own collections beyond the data presented here.

Lawsonia intracellularis is a Gram-negative obligate intracellular bacterium that is the aetiological agent of proliferative enteropathy (PE), a common intestinal disease of major economic importance in pigs and other animal species. To date, progress in understanding the biology of L. intracellularis for improved disease control has been hampered by the inability to culture the organism in vitro. In particular, our understanding of the genomic diversity and population structure of clinical L. intercellularis is very limited. Here, we utilized a metagenomic shotgun approach to directly sequence and assemble 21 L . intracellularis genomes from faecal and ileum samples of infected pigs and horses across three continents. Phylogenetic analysis revealed a genetically monomorphic clonal lineage responsible for infections in pigs, with distinct subtypes associated with infections in horses. The genome was highly conserved, with 94 % of genes shared by all isolates and a very small accessory genome made up of only 84 genes across all sequenced strains. In part, the accessory genome was represented by regions with a high density of SNPs, indicative of recombination events importing novel gene alleles. In summary, our analysis provides the first view of the population structure for L. intracellularis , revealing a single major lineage associated with disease of pigs. The limited diversity and broad geographical distribution suggest the recent emergence and clonal expansion of an important livestock pathogen.

Bacterial sequencing will become increasingly adopted in routine microbiology laboratories. Here, we report the findings of a technical evaluation of almost 800 clinical methicillin-resistant Staphylococcus aureus (MRSA) isolates, in which we sought to define key quality metrics to support MRSA sequencing in clinical practice. We evaluated the accuracy of mapping to a generic reference versus clonal complex (CC)-specific mapping, which is more computationally challenging. Focusing on isolates that were genetically related (50 bp apart to identify same-species contamination for MRSA. These metrics were combined into a quality-control (QC) flowchart to determine whether sequence runs and individual clinical isolates passed QC, which could be adapted by future automated analysis systems to enable rapid hands-off sequence analysis by clinical laboratories.

Haemophilus influenzae causes common and sometimes severe adult and pediatric disease including chronic obstructive respiratory disease, otitis media and infections of the central nervous system. Serotype b strains, with a b-type capsule, have been the historical cause of invasive disease, and the introduction of a serotype b-specific vaccine has led to their decline. However, unencapsulated or non-b-type H. influenzae infections are not prevented by the vaccine and appear to be increasing in frequency. Here we report two pediatric cases of severe central nervous system H. influenzae infection presenting to the same hospital in San Diego, California during the same week in January 2016. Due to good vaccine coverage in this part of the world, H. influenzae cases are normally rare and seeing two cases in the same week was unexpected. We thus suspected a recent transmission chain, and possible local outbreak. To test this hypothesis, we isolated and sequenced whole genomes from each patient and placed them in a phylogenetic tree spanning the known diversity of H. influenzae . Surprisingly, we found that the two isolates (SD2016_1 and SD2016_2) belonged to distantly related lineages, suggesting two independent transmission events and ruling out a local outbreak. Despite being distantly related, the two isolates belong to two different lineages that have exchanged capsule loci in the recent past. Therefore, as in other bacterial pathogens, capsule switching by horizontal gene transfer may be an important evolutionary mechanism of vaccine evasion in H. influenzae .

Linezolid-resistant Enterococcus faecalis (LREfs) carrying optrA are increasingly reported globally from multiple sources, but we lack a comprehensive analysis of human and animal optrA-LREfs strains. To assess if optrA is dispersed in isolates with varied genetic backgrounds or with common genetic features, we investigated the phylogenetic structure, genetic content [antimicrobial resistance (AMR), virulence, prophages, plasmidome] and optrA-containing platforms of 27 publicly available optrA-positive E. faecalis genomes from different hosts in seven countries. At the genome-level analysis, an in-house database with 64 virulence genes was tested for the first time. Our analysis showed a diversity of clones and adaptive gene sequences related to a wide range of genera from Firmicutes . Phylogenies of core and accessory genomes were not congruent, and at least PAI-associated and prophage genes contribute to such differences. Epidemiologically unrelated clones (ST21, ST476-like and ST489) obtained from human clinical and animal hosts in different continents over eight years (2010–2017) could be phylogenetically related (3–126 SNPs difference). optrA was located on the chromosome within a Tn6674-like element (n=10) or on medium-size plasmids (30–60 kb; n=14) belonging to main plasmid families (RepA_N/Inc18/Rep_3). In most cases, the immediate gene vicinity of optrA was generally identical in chromosomal (Tn6674) or plasmid (impB-fexA-optrA) backbones. Tn6674 was always inserted into the same ∆radC integration site and embedded in a 32 kb chromosomal platform common to strains from different origins (patients, healthy humans, and animals) in Europe, Africa, and Asia during 2012–2017. This platform is conserved among hundreds of E. faecalis genomes and proposed as a chromosomal hotspot for optrA integration. The finding of optrA in strains sharing common adaptive features and genetic backgrounds across different hosts and countries suggests the occurrence of common and independent genetic events occurring in distant regions and might explain the easy de novo generation of optrA-positive strains. It also anticipates a dramatic increase of optrA carriage and spread with a serious impact on the efficacy of linezolid for the treatment of Gram-positive infections.

In this work we compared the most frequently used Klebsiella pneumoniae typing methods: PFGE, cgMLST and coreSNP. We evaluated the discriminatory power of the three methods to confirm or exclude nosocomial transmission on K. pneumoniae strains isolated from January to December 2017, in the framework of the routine surveillance for multidrug-resistant organisms at the San Raffaele Hospital, in Milan. We compared the results of the different methods to the results of epidemiological investigation. Our results showed that cgMLST and coreSNP are more discriminant than PFGE, and that both approaches are suitable for transmission analyses. cgMLST appeared to be inferior to coreSNP in the K. pneumoniae CG258 phylogenetic reconstruction. Indeed, we found that the phylogenetic reconstruction based on cgMLST genes wrongly clustered ST258 clade1 and clade2 strains, conversely properly assigned by coreSNP approach. In conclusion, this study provides evidences supporting the reliability of both cgMLST and coreSNP for hospital surveillance programs and highlights the limits of cgMLST scheme genes for phylogenetic reconstructions.

Several members of the Gram-negative environmental bacterial genus Achromobacter are associated with serious infections, with Achromobacter xylosoxidans being the most common. Despite their pathogenic potential, little is understood about these intrinsically drug-resistant bacteria and their role in disease, leading to suboptimal diagnosis and management. Here, we performed comparative genomics for 158 Achromobacter spp. genomes to robustly identify species boundaries, reassign several incorrectly speciated taxa and identify genetic sequences specific for the genus Achromobacter and for A. xylosoxidans . Next, we developed a Black Hole Quencher probe-based duplex real-time PCR assay, Ac-Ax, for the rapid and simultaneous detection of Achromobacter spp. and A. xylosoxidans from both purified colonies and polymicrobial clinical specimens. Ac-Ax was tested on 119 isolates identified as Achromobacter spp. using phenotypic or genotypic methods. In comparison to these routine diagnostic methods, the duplex assay showed superior identification of Achromobacter spp. and A. xylosoxidans , with five Achromobacter isolates failing to amplify with Ac-Ax confirmed to be different genera according to 16S rRNA gene sequencing. Ac-Ax quantified both Achromobacter spp. and A. xylosoxidans down to ~110 genome equivalents and detected down to ~12 and ~1 genome equivalent(s), respectively. Extensive in silico analysis, and laboratory testing of 34 non- Achromobacter isolates and 38 adult cystic fibrosis sputa, confirmed duplex assay specificity and sensitivity. We demonstrate that the Ac-Ax duplex assay provides a robust, sensitive and cost-effective method for the simultaneous detection of all Achromobacter spp. and A. xylosoxidans and will facilitate the rapid and accurate diagnosis of this important group of pathogens.

In this work, we used a whole-genome sequencing (WGS) approach to study the features of KPC-producing Klebsiella pneumoniae (KPC-Kp) spreading in a large Italian long-term acute-care rehabilitation facility (LTACRF), and to track the dynamics of dissemination within this setting. Thirty-eight, non-replicated, KPC-Kp isolates from colonized patients (either already colonized at admission or colonized during admission), collected during 2016, were subjected to antimicrobial-susceptibility testing and WGS. All isolates were resistant to β-lactams, with the exception of ceftazidime/avibactam (97.4 % susceptible). The second most effective agent was fosfomycin, followed by colistin, trimethoprim/sulfamethoxazole, gentamicin and amikacin (92.1, 86.8, 60.5, 44.7 and 50 % of susceptibility, respectively). A large proportion of isolates (n=18/38, 47.4%) belonged to clonal group (CG) 101, and most of them (n=15) to a new sequence type (ST) designated as ST2502. All the CG101 isolates had a capsule locus type KL17. The ST2502 harboured the genes encoding for the yersiniabactin siderophore and the ArmA methylase, conferring high-level resistance to aminoglycosides. The second most represented lineage of isolates (16/38, 42.1%) belonged to ST512 of CG258. Analysing WGS data, we were able to ascertain the common origin of some isolates imported from other hospitals, and to track several clusters of in-LTACRF cross-transmissions. The results revealed that, in peculiar epidemiological settings such as LTACRF, new KPC-Kp clones different from those prevailing in acute-care hospitals and associated with uncommon resistance and virulence determinants can successfully emerge and disseminate.

This study presents the nationwide epidemiology of Neisseria gonorrhoeae , using whole-genome sequencing of all culture-positive cases, which comprise roughly 40 % of all cases of gonorrhea reported in Norway from 2016 to 2017. Isolates were assigned to sequence types and Bayesian analysis clusters and variation in genes coding for antibiotic resistance was linked to phenotypic resistance data. The study also included isolates taken from the same patients from different anatomical sites at one or more time points. Comparing these isolates allows for observation of patterns of infections, i.e. multiple reinfections of genetically related clones vs. reinfections of genetically distant clones, and quantification of the genomic variation of closely related isolates from samples taken from a patient within the same day. Demographically, the patients in the study could be split into two groups; one group of patients from the capital with a high proportion of men who have sex with men (MSM), and another consisting of young adults with transmission primarily between males and females from outside the capital. Some clusters of N. gonorrhoeae were restricted to one of these two demographic groups. Pairwise comparison of multiple isolates from the same patients revealed that most were reinfected with different clones. Observations of frequent reinfections in patients is a concern and should be taken into account in the development of improved information and treatment guidelines.

Resistance to meticillin and vancomycin in Staphylococcus aureus significantly complicates the management of severe infections like bacteraemia, endocarditis or osteomyelitis. Here, we review the molecular mechanisms and genomic epidemiology of resistance to these agents, with a focus on how genomics has provided insights into the emergence and evolution of major meticillin-resistant S. aureus clones. We also provide insights on the use of bacterial whole-genome sequencing to inform management of S. aureus infections and for control of transmission at the hospital and in the community.

The heterogeneous and highly recombinogenic genus Haemophilus comprises several species, some of which are pathogenic to humans. All share an absolute requirement for blood-derived factors during growth. Certain species, such as the pathogen Haemophilus influenzae and the commensal Haemophilus haemolyticus, are thought to require both haemin (X-factor) and nicotinamide adenine dinucleotide (NAD, V-factor), whereas others, such as the informally classified 'Haemophilus intermedius subsp. intermedius', and Haemophilus parainfluenzae, only require V-factor. These differing growth requirements are commonly used for species differentiation, although a number of studies are now revealing issues with this approach. Here, we perform large-scale phylogenomics of 240 Haemophilus spp. genomes, including five 'H. intermedius' genomes generated in the current study, to reveal that strains of the 'H. intermedius' group are in fact haemin-independent H. haemolyticus (hiHh). Closer examination of these hiHh strains revealed that they encode an intact haemin biosynthesis pathway, unlike haemin-dependent H. haemolyticus and H. influenzae, which lack most haemin biosynthesis genes. Our results suggest that the common ancestor of modern-day H. haemolyticus and H. influenzae lost key haemin biosynthesis loci, likely as a consequence of specialized adaptation to otorhinolaryngeal and respiratory niches during their divergence from H. parainfluenzae. Genetic similarity analysis demonstrated that the haemin biosynthesis loci acquired in the hiHh lineage were likely laterally transferred from a H. parainfluenzae ancestor, and that this event probably occurred only once in hiHh. This study further challenges the validity of phenotypic methods for differentiating among Haemophilus species, and highlights the need for whole-genome sequencing for accurate characterization of species within this taxonomically challenging genus.

Vibrio parahaemolyticus is an important cause of foodborne gastroenteritis globally. Thermostable direct haemolysin (TDH) and the TDH-related haemolysin are the two key virulence factors in V. parahaemolyticus. Vibrio pathogenicity islands harbour the genes encoding these two haemolysins. The serotyping of V. parahaemolyticus is based on the combination of O and K antigens. Frequent recombination has been observed in V. parahaemolyticus , including in the genomic regions encoding the O and K antigens. V. parahaemolyticus serotype O4:K12 has caused gastroenteritis outbreaks in the USA and Spain. Recently, outbreaks caused by this serotype of V. parahaemolyticus have been reported in China. However, the relationships among this serotype of V. parahaemolyticus strains isolated in different regions have not been addressed. Here, we investigated the genome variation of the V. parahaemolyticus serotype O4:K12 using the whole-genome sequences of 29 isolates. We determined five distinct lineages in this strain collection. We observed frequent recombination among different lineages. In contrast, little recombination was observed within each individual lineage. We showed that the lineage of this serotype of V. parahaemolyticus isolated in America was different from those isolated in Asia and identified genes that exclusively existed in the strains isolated in America. Pan-genome analysis showed that strain-specific and cluster-specific genes were mostly located in the genomic islands. Pan-genome analysis also showed that the vast majority of the accessory genes in the O4:K12 serotype of V. parahaemolyticus were acquired from within the genus Vibrio . Hence, we have shown that multiple distinct lineages exist in V. parahaemolyticus serotype O4:K12 and have provided more evidence about the gene segregation found in V. parahaemolyticus isolated in different continents.

Stenotrophomonas maltophilia is emerging as an important cause of disease in nosocomial and community-acquired settings, including bloodstream, wound and catheter-associated infections. Cystic fibrosis (CF) airways also provide optimal growth conditions for various opportunistic pathogens with high antibiotic tolerance, including S. maltophilia . Currently, there is no rapid, cost-effective and accurate molecular method for detecting this potentially life-threatening pathogen, particularly in polymicrobial specimens, suggesting that its true prevalence is underestimated. Here, we used large-scale comparative genomics to identify a specific genetic target for S. maltophilia , with subsequent development and validation of a real-time PCR assay for its detection. Analysis of 167 Stenotrophomonas spp. genomes identified a conserved 4 kb region in S. maltophilia , which was targeted for Black Hole Quencher assay design. Our assay yielded the positive detection of 89 of 89 (100%) clinical S. maltophilia strains, and no amplification of 23 non- S. maltophilia clinical isolates. S. maltophilia was detected in 10 of 16 CF sputa, demonstrating the assay's utility for direct detection in respiratory specimens. The assay demonstrated good sensitivity, with limits of detection and quantitation on pure culture of ~10 and ~100 genome equivalents, respectively. Our assay provides a highly specific, sensitive and cost-effective method for the accurate identification of S. maltophilia , and will improve the diagnosis and treatment of this under-recognized pathogen by enabling its accurate and rapid detection from polymicrobial clinical and environmental samples.

Most pneumococci express a polysaccharide capsule, a key virulence factor and target for pneumococcal vaccines. However, pneumococci showing no serological evidence of capsule expression [nontypeable pneumococci (NTPn)] are more frequently isolated from carriage studies than in invasive disease. Limited data exist about the population structure of carriage NTPn from the African continent. We aimed to characterize carriage NTPn and compare them to previously described invasive NTPn. Carriage and invasive NTPn isolates were obtained from South African cross-sectional studies (2009 and 2012) and laboratory-based surveillance for invasive pneumococcal disease (2003–2013), respectively. Isolates were characterized by capsular locus sequence analysis, multilocus sequence typing, antimicrobial non-susceptibility patterns and phylogenetic analysis. NTPn represented 3.7 % (137/3721) of carriage isolates compared to 0.1 % (39/32 824) of invasive isolates (P

The interplay between food production animals, humans and the environment with respect to the transmission of drug-resistant pathogens is widely debated and poorly understood. Pandemic uropathogenic Escherichia coli ST131-H30Rx, with conserved fluoroquinolone and cephalosporin resistance, are not frequently identified in animals. However, the phylogenetic precursor lineage ST131-H22 in animals and associated meat products is being reported with increasing frequency. Here we characterized two highly related ST131-H22 strains, one from a healthy pig and the other from a human infection (in 2007 and 2009, respectively). We used both long and short genome sequencing and compared them to ST131-H22 genome sequences available in public repositories. Even within the context of H22 strains, the two strains in question were highly related, separated by only 20 core SNPs. Furthermore, they were closely related to a faecal strain isolated in 2010 from a geographically distinct, healthy human in New South Wales, Australia. The porcine and hospital strains carried highly similar HI2-ST3 multidrug resistant plasmids with differences in the hospital strain arising due to IS-mediated insertions and rearrangements. Near identical ColV plasmids were also present in both strains, further supporting their shared evolutionary history. This work highlights the importance of adopting a One Health approach to genomic surveillance to gain insights into pathogen evolution and spread.

Malaria was present in most of Europe until the second half of the 20th century, when it was eradicated through a combination of increased surveillance and mosquito control strategies, together with cross-border and political collaboration. Despite the severe burden of malaria on human populations, it remains contentious how the disease arrived and spread in Europe. Here, we report a partial Plasmodium falciparum nuclear genome derived from a set of antique medical slides stained with the blood of malaria-infected patients from Spain’s Ebro Delta, dating to the 1940s. Our analyses of the genome of this now eradicated European P. falciparum strain confirms stronger phylogeographical affinity to present-day strains in circulation in central south Asia, rather than to those in Africa. This points to a longitudinal, rather than a latitudinal, spread of malaria into Europe. In addition, this genome displays two derived alleles in the pfmrp1 gene that have been associated with drug resistance. Whilst this could represent standing variation in the ancestral P. falciparum population, these mutations may also have arisen due to the selective pressure of quinine treatment, which was an anti-malarial drug already in use by the time the sample we sequenced was mounted on a slide., C. L. -F. is supported by Obra Social ‘La Caixa’, Secretaria d’Universitats i Recerca Programme del Departament d’Economia i Coneixement de la Generalitat de Catalunya (GRC 2017 SGR 880) and by FEDER-MINECO (PGC2018-095931-B-100). L. v. D. and F. B. acknowledge financial support from the MRC Newton Fund UK–China NSFC initiative (grants MR/P007597/1 and 81661138006).

Burkholderia pseudomallei is a Gram-negative saprophytic bacillus and the aetiological agent of melioidosis, a disease of public-health importance throughout Southeast Asia and northern Australia. Infection can occur in humans and a wide array of animal species, though zoonotic transmission and case clusters are rare. Despite its highly plastic genome and extensive strain diversity, fine-scale investigations into the population structure of B. pseudomallei indicate there is limited geographical dispersal amongst sequence types (STs). In the ‘Top End’ of northern Australia, five STs comprise 90 % of the overall abundance, the most prevalent and widespread of which is ST-109. In May 2016, ST-109 was implicated in two fatal cases of melioidosis in juvenile saltwater crocodiles at a wildlife park near Darwin, Australia. To determine the probable source of infection, we sampled the crocodile enclosures and analysed the phylogenetic relatedness of crocodile and culture-positive ST-109 environmental park isolates against an additional 135 ST-109 B. pseudomallei isolates from the Top End. Collectively, our whole-genome sequencing (WGS) and pathology findings confirmed B. pseudomallei detected in the hatchling incubator as the likely source of infection, with zero SNPs identified between clinical and environmental isolates. Our results also demonstrate little variation across the ST-109 genome, with SNPs in recombinogenic regions and one suspected case of ST homoplasy accounting for nearly all observed diversity. Collectively, this study supports the use of WGS for outbreak source attribution in highly recombinogenic pathogens, and confirms the epidemiological and phylogenetic insights that can be gained from high-resolution sequencing platforms.

We undertook a comprehensive comparative analysis of a collection of 30 small (

Bacteriophages are the most prevalent biological entities impacting on the ecosystem and are characterized by their extensive diversity. However, there are two aspects of phages that have remained largely unexplored: genetic flux by recombination between phage populations and characterization of specific phages in terms of the pan-genome. Here, we examined the recombination and pan-genome in Helicobacter pylori prophages at both the genome and gene level. In the genome-level analysis, we applied, for the first time, chromosome painting and fineSTRUCTURE algorithms to a phage species, and showed novel trends in inter-population genetic flux. Notably, hpEastAsia is a phage population that imported a higher proportion of DNA fragments from other phages, whereas the hpSWEurope phages showed weaker signatures of inter-population recombination, suggesting genetic isolation. The gene-level analysis showed that, after parameter tuning of the prokaryote pan-genome analysis program, H. pylori phages have a pan-genome consisting of 75 genes and a soft-core genome of 10 genes, which includes genes involved in the lytic and lysogenic life cycles. Quantitative analysis of recombination events of the soft-core genes showed no substantial variation in the intensity of recombination across the genes, but rather equally frequent recombination among housekeeping genes that were previously reported to be less prone to recombination. The signature of frequent recombination appears to reflect the host–phage evolutionary arms race, either by contributing to escape from bacterial immunity or by protecting the host by producing defective phages.

The pneumococcus produces a polysaccharide capsule, encoded by the cps locus, that provides protection against phagocytosis and determines serotype. Nearly 100 serotypes have been identified with new serotypes still being discovered, especially in previously understudied regions. Here we present an analysis of the cps loci of more than 18 000 genomes from the Global Pneumococcal Sequencing (GPS) project with the aim of identifying novel cps loci with the potential to produce previously unrecognized capsule structures. Serotypes were assigned using whole genome sequence data and 66 of the approximately 100 known serotypes were included in the final dataset. Closer examination of each serotype’s sequences identified nine putative novel cps loci (9X, 11X, 16X, 18X1, 18X2, 18X3, 29X, 33X and 36X) found in ~2.6 % of the genomes. The large number and global distribution of GPS genomes provided an unprecedented opportunity to identify novel cps loci and consider their phylogenetic and geographical distribution. Nine putative novel cps loci were identified and examples of each will undergo subsequent structural and immunological analysis. Fil: van Tonder, Andries. Wellcome Genome Campus; Fil: Gladstone, Rebecca A.. Wellcome Genome Campus; Fil: Lo, Stephanie W.. Wellcome Genome Campus; Fil: Nahm, Moon H.. University Of Alabama At Birmingahm. School Of Medicine; Estados Unidos Fil: du Plessis, Mignon. National Institute For Communicable Diseases; Sudáfrica. University of the Witwatersrand; Sudáfrica Fil: Cornick, Jennifer. Malawi Liverpool Wellcome Trust; Malaui Fil: Kwambana Adams, Brenda. Vaccines And Immunity Theme; Gambia Fil: Madhi, Shabir A.. University of the Witwatersrand; Sudáfrica Fil: Hawkins, Paulina A.. University Of Emory. Rollins School Of Public Health; Estados Unidos Fil: Benisty, Rachel. Ben Gurion University of the Negev; Israel Fil: Dagan, Ron. Ben Gurion University of the Negev; Israel Fil: Everett, Dean. University of Edinburgh; Reino Unido Fil: Antonio, Martin. Vaccines And Immunity Theme; Gambia Fil: Klugman, Keith P.. University Of Emory. Rollins School Of Public Health; Estados Unidos Fil: von Gottberg, Anne. National Institute For Communicable Diseases; Sudáfrica. University of the Witwatersrand; Sudáfrica Fil: Breiman, Robert F.. University Of Emory. Rollins School Of Public Health; Estados Unidos. University of Edinburgh; Reino Unido Fil: McGee, Lesley. Centers for Disease Control and Prevention; Estados Unidos Fil: Bentley, Stephen D.. Wellcome Genome Campus; . University Of Emory. Emory College Of Arts And Sciences.; Estados Unidos Fil: Faccone, Diego Francisco. Administracion Nacional de Laboratorios E Institutos de Salud "dr. Carlos G. Malbran". Instituto Nacional de Epidemiologia. Departamento de Investigacion.; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina

Salmonella enterica serotype Kentucky can be a common causative agent of salmonellosis, usually associated with consumption of contaminated poultry. Antimicrobial resistance (AMR) to multiple drugs, including ciprofloxacin, is an emerging problem within this serotype. We used whole-genome sequencing (WGS) to investigate the phylogenetic structure and AMR content of 121 S. e nterica serotype Kentucky sequence type 198 isolates from five continents. Population structure was inferred using phylogenomic analysis and whole genomes were compared to investigate changes in gene content, with a focus on acquired AMR genes. Our analysis showed that multidrug-resistant (MDR) S. enterica serotype Kentucky isolates belonged to a single lineage, which we estimate emerged circa 1989 following the acquisition of the AMR-associated Salmonella genomic island (SGI) 1 (variant SGI1-K) conferring resistance to ampicillin, streptomycin, gentamicin, sulfamethoxazole and tetracycline. Phylogeographical analysis indicates this clone emerged in Egypt before disseminating into Northern, Southern and Western Africa, then to the Middle East, Asia and the European Union. The MDR clone has since accumulated various substitution mutations in the quinolone-resistance-determining regions (QRDRs) of DNA gyrase (gyrA) and DNA topoisomerase IV (parC), such that most strains carry three QRDR mutations which together confer resistance to ciprofloxacin. The majority of AMR genes in the S. e nterica serotype Kentucky genomes were carried either on plasmids or SGI structures. Remarkably, each genome of the MDR clone carried a different SGI1-K derivative structure; this variation could be attributed to IS26-mediated insertions and deletions, which appear to have hampered previous attempts to trace the clone’s evolution using sub-WGS resolution approaches. Several different AMR plasmids were also identified, encoding resistance to chloramphenicol, third-generation cephalosporins, carbapenems and/or azithromycin. These results indicate that most MDR S. e nterica serotype Kentucky circulating globally result from the clonal expansion of a single lineage that acquired chromosomal AMR genes 30 years ago, and has continued to diversify and accumulate additional resistances to last-line oral antimicrobials. This article contains data hosted by Microreact.

The majority of Acinetobacter baumannii isolates that are multiply, extensively and pan-antibiotic resistant belong to two globally disseminated clones, GC1 and GC2, that were first noticed in the 1970s. Here, we investigated microevolution and phylodynamics within GC1 via analysis of 45 whole-genome sequences, including 23 sequenced for this study. The most recent common ancestor of GC1 arose around 1960 and later diverged into two phylogenetically distinct lineages. In the 1970s, the main lineage acquired the AbaR resistance island, conferring resistance to older antibiotics, via a horizontal gene transfer event. We estimate a mutation rate of ∼5 SNPs genome− 1 year− 1 and detected extensive recombination within GC1 genomes, introducing nucleotide diversity into the population at >20 times the substitution rate (the ratio of SNPs introduced by recombination compared with mutation was 22). The recombination events were non-randomly distributed in the genome and created significant diversity within loci encoding outer surface molecules (including the capsular polysaccharide, the outer core lipooligosaccharide and the outer membrane protein CarO), and spread antimicrobial resistance-conferring mutations affecting the gyrA and parC genes and insertion sequence insertions activating the ampC gene. Both GC1 lineages accumulated resistance to newer antibiotics through various genetic mechanisms, including the acquisition of plasmids and transposons or mutations in chromosomal genes. Our data show that GC1 has diversified into multiple successful extensively antibiotic-resistant subclones that differ in their surface structures. This has important implications for all avenues of control, including epidemiological tracking, antimicrobial therapy and vaccination.

Enterococcus faeciumis a gut commensal of humans and animals. In addition, it has recently emerged as an important nosocomial pathogen through the acquisition of genetic elements that confer resistance to antibiotics and virulence. We performed a whole-genome sequencing-based study on 96 multidrug-resistantE. faeciumstrains that asymptomatically colonized five patients with the aim of describing the genome dynamics of this species. The patients were hospitalized on multiple occasions and isolates were collected over periods ranging from 15 months to 6.5 years. Ninety-five of the sequenced isolates belonged toE. faeciumclade A1, which was previously determined to be responsible for the vast majority of clinical infections. The clade A1 strains clustered into six clonal groups of highly similar isolates, three of which consisted entirely of isolates from a single patient. We also found evidence of concurrent colonization of patients by multiple distinct lineages and transfer of strains between patients during hospitalization. We estimated the evolutionary rate of two clonal groups that each colonized single patients at 12.6 and 25.2 single-nucleotide polymorphisms (SNPs)/genome/year. A detailed analysis of the accessory genome of one of the clonal groups revealed considerable variation due to gene gain and loss events, including the chromosomal acquisition of a 37 kbp prophage and the loss of an element containing carbohydrate metabolism-related genes. We determined the presence and location of 12 different insertion sequence (IS) elements, with ISEfa5showing a unique pattern of location in 24 of the 25 isolates, suggesting widespread ISEfa5excision and insertion into the genome during gut colonization. Our findings show that theE. faeciumgenome is highly dynamic during asymptomatic colonization of the human gut. We observed considerable genomic flexibility due to frequent horizontal gene transfer and recombination, which can contribute to the generation of genetic diversity within the species and, ultimately, can contribute to its success as a nosocomial pathogen.

Pseudomonas aeruginosa is a highly versatile, antibiotic-resistant Gram-negative bacterium known for causing opportunistic infections and contamination of industrial products. Despite extensive genomic analysis of clinical P. aeruginosa strains, no genomes exist for preservative-tolerant industrial strains. A unique collection of 69 industrial isolates was assembled and compared to clinical and environmental strains; 16 genetically distinct industrial strains were subjected to array tube genotyping, multilocus sequence typing and whole-genome sequencing. The industrial strains possessed high preservative tolerance and were dispersed widely across P. aeruginosa as a species, but recurrence of strains from the same lineage within specific industrial products and locations was identified. The industrial P. aeruginosa genomes (mean=7.0 Mb) were significantly larger than those of previously sequenced environmental (mean=6.5 Mb; n=19) and clinical (mean=6.6 Mb; n=66) strains. Complete sequencing of the P. aeruginosa industrial strain RW109, which encoded the largest genome (7.75 Mb), revealed a multireplicon structure including a megaplasmid (555 265 bp) and large plasmid (151 612 bp). The RW109 megaplasmid represented an emerging plasmid family conserved in seven industrial and two clinical P. aeruginosa strains, and associated with extremely stress-resilient phenotypes, including antimicrobial resistance and solvent tolerance. Here, by defining the detailed phylogenomics of P. aeruginosa industrial strains, we show that they uniquely possess multireplicon, megaplasmid-bearing genomes, and significantly greater genomic content worthy of further study.

Clostridioides difficileis a ubiquitous, diarrhoeagenic pathogen often associated with healthcare-acquired infections that can cause a range of symptoms from mild, self-limiting disease to toxic megacolon and death. Since the early 2000s, a large proportion ofC. difficilecases have been attributed to the ribotype 027 (RT027) lineage, which is associated with sequence type 1 (ST1) in theC. difficilemultilocus sequence typing scheme. The spread of ST1 has been attributed, in part, to resistance to fluoroquinolones used to treat unrelated infections, which creates conditions ideal forC. difficilecolonization and proliferation. In this study, we analysed 27 isolates from a healthcare network in northern Arizona, USA, and 1352 publicly available ST1 genomes to place locally sampled isolates into a global context. Whole genome, single nucleotide polymorphism analysis demonstrated that at least six separate introductions of ST1 were observed in healthcare facilities in northern Arizona over an 18-month sampling period. A reconstruction of transmission networks identified potential nosocomial transmission of isolates, which were only identified via whole genome sequence analysis. Antibiotic resistance heterogeneity was observed among ST1 genomes, including variability in resistance profiles among locally sampled ST1 isolates. To investigate why ST1 genomes are so common globally and in northern Arizona, we compared all high-qualityC. difficilegenomes and identified that ST1 genomes have gained and lost a number of genomic regions compared to all otherC. difficilegenomes; analyses of other toxigenicC. difficilesequence types demonstrate that this loss may be anomalous and could be related to niche specialization. These results suggest that a combination of antimicrobial resistance and gain and loss of specific genes may explain the prominent association of this sequence type withC. difficileinfection cases worldwide. The degree of genetic variability in ST1 suggests that classifying all ST1 genomes into a quinolone-resistant hypervirulent clone category may not be appropriate. Whole genome sequencing of clinicalC. difficileisolates provides a high-resolution surveillance strategy for monitoring persistence and transmission ofC. difficileand for assessing the performance of infection prevention and control strategies.

Shigella species are a major cause of gastroenteritis worldwide, and Shigella sonnei is the most common species isolated within the United States. Previous surveillance work in Pennsylvania documented increased antimicrobial resistance (AMR) in S. sonnei associated with reported illnesses. The present study examined a subset of these isolates by whole genome sequencing (WGS) to determine the relationship between domestic and international isolates, to identify genes that may be useful for identifying specific Global Lineages of S. sonnei and to test the accuracy of WGS for predicting AMR phenotype. A collection of 22 antimicrobial-resistant isolates from patients infected within the United States or while travelling internationally between 2009 and 2014 was chosen for WGS. Phylogenetic analysis revealed both international and domestic isolates were one of two previously defined Global Lineages of S. sonnei, designated Lineage II and Lineage III. Twelve of 17 alleles tested distinguish these two lineages. Lastly, genome analysis was used to identify AMR determinants. Genotypic analysis was concordant with phenotypic resistance for six of eight antibiotic classes. For aminoglycosides and trimethoprim, resistance genes were identified in two and three phenotypically sensitive isolates, respectively. This article contains data hosted by Microreact.

The standard workhorse for genomic analysis of the evolution of bacterial populations is phylogenetic modelling of mutations in the core genome. However, a notable amount of information about evolutionary and transmission processes in diverse populations can be lost unless the accessory genome is also taken into consideration. Here, we introduce PANINI (Pangenome Neighbour Identification for Bacterial Populations), a computationally scalable method for identifying the neighbours for each isolate in a data set using unsupervised machine learning with stochastic neighbour embedding based on the t-SNE (t-distributed stochastic neighbour embedding) algorithm. PANINI is browser-based and integrates with the Microreact platform for rapid online visualization and exploration of both core and accessory genome evolutionary signals, together with relevant epidemiological, geographical, temporal and other metadata. Several case studies with single- and multi-clone pneumococcal populations are presented to demonstrate the ability to identify biologically important signals from gene content data. PANINI is available at http://panini.pathogen.watch and code at http://gitlab.com/cgps/panini.

The 13-valent pneumococcal conjugate vaccine (PCV-13) was introduced in the United States in 2010. Using a large paediatric carriage sample collected from shortly after the introduction of PCV-7 to several years after the introduction of PCV-13, we investigate alterations in the composition of the pneumococcal population following the introduction of PCV-13, evaluating the extent to which the post-vaccination non-vaccine type (NVT) population mirrors that from prior to vaccine introduction and the effect of PCV-13 on vaccine type lineages. Draft genome assemblies from 736 newly sequenced and 616 previously published pneumococcal carriage isolates from children in Massachusetts between 2001 and 2014 were analysed. Isolates were classified into one of 22 sequence clusters (SCs) on the basis of their core genome sequence. We calculated the SC diversity for each sampling period as the probability that any two randomly drawn isolates from that period belong to different SCs. The sampling period immediately after the introduction of PCV-13 (2011) was found to have higher diversity than preceding (2007) or subsequent (2014) sampling periods {Simpson’s D 2007: 0.915 [95 % confidence interval (CI) 0.901, 0.929]; 2011: 0.935 [0.927, 0.942]; 2014 : 0.912 [0.901, 0.923]}. Amongst NVT isolates, we found the distribution of SCs in 2011 to be significantly different from that in 2007 or 2014 (Fisher’s exact test P=0.018, 0.0078), but did not find a difference comparing 2007 to 2014 (Fisher’s exact test P=0.24), indicating greater similarity between samples separated by a longer time period than between samples from closer time periods. We also found changes in the accessory gene content of the NVT population between 2007 and 2011 to have been reduced by 2014. Amongst the new serotypes targeted by PCV-13, four were present in our sample. The proportion of our sample composed of PCV-13-only vaccine serotypes 19A, 6C and 7F decreased between 2007 and 2014, but no such reduction was seen for serotype 3. We did, however, observe differences in the genetic composition of the pre- and post-PCV-13 serotype 3 population. Our isolates were collected during discrete sampling periods from a small geographical area, which may limit the generalizability of our findings. Pneumococcal diversity increased immediately following the introduction of PCV-13, but subsequently returned to pre-vaccination levels. This is reflected in the distribution of NVT lineages, and, to a lesser extent, their accessory gene frequencies. As such, there may be a period during which the population is particularly disrupted by vaccination before returning to a more stable distribution. The persistence and shifting genetic composition of serotype 3 is a concern and warrants further investigation.

Acinetobacter baumannii causes severe, fulminant, community-acquired pneumonia (CAP) in tropical and subtropical regions. We compared the population structure, virulence and antimicrobial resistance determinants of northern Australian community-onset A. baumannii strains with local and global strains. We performed whole-genome sequencing on 55 clinical and five throat colonization A. baumannii isolates collected in northern Australia between 1994 and 2016. Clinical isolates included CAP (n=41), healthcare-associated pneumonia (n=7) and nosocomial bloodstream (n=7) isolates. We also included 93 publicly available international A. baumannii genome sequences in the analyses. Patients with A. baumannii CAP were almost all critically unwell; 82 % required intensive care unit admission and 18 % died during their inpatient stay. Whole-genome phylogenetic analysis demonstrated that community-onset strains were not phylogenetically distinct from nosocomial strains. Some non-multidrug-resistant local strains were closely related to multidrug-resistant strains from geographically distant locations. Pasteur sequence type (ST)10 was the dominant ST and accounted for 31/60 (52 %) northern Australian strains; the remainder belonged to a diverse range of STs. The most recent common ancestor for ST10 was estimated to have occurred in 1738 (95 % highest posterior density, 1626–1826), with evidence of multiple introduction events between Australia and Southeast Asia between then and the present day. Virulence genes associated with biofilm formation and the type 6 secretion system (T6SS) were absent in many strains, and were not associated with in-hospital mortality. All strains were susceptible to gentamicin and meropenem; none carried an AbaR resistance island. Our results suggest that international dissemination of A. baumannii is occurring in the community on a contemporary timescale. Genes associated with biofilm formation and the T6SS may not be required for survival in community niches. The relative contributions of host and bacterial factors to the clinical severity of community-onset A. baumannii infection require further investigation.

Gonorrhoea infections are on the increase and strains that are resistant to all antimicrobials used to treat the disease have been found worldwide. These observations encouraged the World Health Organization to include Neisseria gonorrhoeae on their list of high-priority organisms in need of new treatments. Fortunately, concurrent resistance to both antimicrobials used in dual therapy is still rare. The fight against antimicrobial resistance (AMR) must begin from an understanding of how it evolves and spreads in sexual networks. Genome-based analyses have allowed the study of the gonococcal population dynamics and transmission, giving a novel perspective on AMR gonorrhoea. Here, we will review past, present and future treatment options for gonorrhoea and explain how genomics is helping to increase our understanding of the changing AMR and transmission landscape. This article contains data hosted by Microreact.

© 2019 The Authors. Avian pathogenic Escherichia coli (APEC) cause widespread economic losses in poultry production and are potential zoonotic pathogens. Genome sequences of 95 APEC from commercial poultry operations in four Australian states that carried the class 1 integrase gene intI1, a proxy for multiple drug resistance (MDR), were characterized. Sequence types ST117 (22/95), ST350 (10/95), ST429 and ST57 (each 9/95), ST95 (8/95) and ST973 (7/95) dominated, while 24 STs were represented by one or two strains. FII and FIB repA genes were the predominant (each 93/95, 98 %) plasmid incompatibility groups identified, but those of B/O/K/Z (25/95, 26 %) and I1 (24/95, 25 %) were also identified frequently. Virulence-associated genes (VAGs) carried by ColV and ColBM virulence plasmids, including those encoding protectins [iss (91/95, 96 %), ompT (91/95, 96 %) and traT (90/95, 95 %)], iron-acquisition systems [sitA (88/95, 93 %), etsA (87/95, 92 %), iroN (84/95, 89 %) and iucD/iutA (84/95, 89 %)] and the putative avian haemolysin hylF (91/95, 96 %), featured prominently. Notably, mobile resistance genes conferring resistance to fluoroquinolones, colistin, extended-spectrum b-lactams and carbapenems were not detected in the genomes of these 95 APEC but carriage of the sulphonamide resistance gene, sul1 (59/95, 63 %), the trimethoprim resistance gene cassettes dfrA5 (48/95, 50 %) and dfrA1 (25/95, 27 %), the tetracycline resistance determinant tet(A) (51/95, 55 %) and the ampicillin resistance genes bla TEM-1A/B/C (48/95, 52 %) was common. IS26 (77/95, 81 %), an insertion element known to capture and mobilize a wide spectrum of antimicrobial resistance genes, was also frequently identified. These studies provide a baseline snapshot of drug-resistant APEC in Australia and their role in the carriage of ColV-like virulence plasmids.

Campylobacter is the leading cause of bacterial enteritis in the developed world, and infections with the organism are largely sporadic in nature. Links between sporadic cases have not been established, with the majority of infections thought to be caused by genetically distinct isolates. Using a read-mapping approach, 158 clinical isolates collected during 2014 from the greater Nottinghamshire area were analysed to assess the local population structure and investigate potential case linkages between sporadic cases of campylobacteriosis. Four instances (2.5 %) of case linkage were observed across the dataset. This study demonstrates that case linkage does occur between sporadic Campylobacter infections, and provides evidence that a dual multi-locus sequence typing/within-lineage single nucleotide polymorphism typing approach to Campylobacter genomic epidemiology provides a benefit to public-health investigations.

The first extensively drug resistant (XDR) Neisseria gonorrhoeae strain with high resistance to the extended-spectrum cephalosporin ceftriaxone was identified in 2009 in Japan, but no other strain with this antimicrobial-resistance profile has been reported since. However, surveillance to date has been based on phenotypic methods and sequence typing, not genome sequencing. Therefore, little is known about the local population structure at the genomic level, and how resistance determinants and lineages are distributed and evolve. We analysed the whole-genome sequence data and the antimicrobial-susceptibility testing results of 204 strains sampled in a region where the first XDR ceftriaxone-resistant N. gonorrhoeae was isolated, complemented with 67 additional genomes from other time frames and locations within Japan. Strains resistant to ceftriaxone were not found, but we discovered a sequence type (ST)7363 sub-lineage susceptible to ceftriaxone and cefixime in which the mosaic penA allele responsible for reduced susceptibility had reverted to a susceptible allele by recombination. Approximately 85 % of isolates showed resistance to fluoroquinolones (ciprofloxacin) explained by linked amino acid substitutions at positions 91 and 95 of GyrA with 99 % sensitivity and 100 % specificity. Approximately 10 % showed resistance to macrolides (azithromycin), for which genetic determinants are less clear. Furthermore, we revealed different evolutionary paths of the two major lineages: single acquisition of penA X in the ST7363-associated lineage, followed by multiple independent acquisitions of the penA X and XXXIV in the ST1901-associated lineage. Our study provides a detailed picture of the distribution of resistance determinants and disentangles the evolution of the two major lineages spreading worldwide.

Orientia tsutsugamushi, formerly Rickettsia tsutsugamushi, is an obligate intracellular pathogen that causes scrub typhus, an underdiagnosed acute febrile disease with high morbidity. Scrub typhus is transmitted by the larval stage (chigger) of Leptotrombidium mites and is irregularly distributed across endemic regions of Asia, Australia and islands of the western Pacific Ocean. Previous work to understand population genetics in O. tsutsugamushi has been based on sub-genomic sampling methods and whole-genome characterization of two genomes. In this study, we compared 40 genomes from geographically dispersed areas and confirmed patterns of extensive homologous recombination likely driven by transposons, conjugative elements and repetitive sequences. High rates of lateral gene transfer (LGT) among O. tsutsugamushi genomes appear to have effectively eliminated a detectable clonal frame, but not our ability to infer evolutionary relationships and phylogeographical clustering. Pan-genomic comparisons using 31 082 high-quality bacterial genomes from 253 species suggests that genomic duplication in O. tsutsugamushi is almost unparalleled. Unlike other highly recombinant species where the uptake of exogenous DNA largely drives genomic diversity, the pan-genome of O. tsutsugamushi is driven by duplication and divergence. Extensive gene innovation by duplication is most commonly attributed to plants and animals and, in contrast with LGT, is thought to be only a minor evolutionary mechanism for bacteria. The near unprecedented evolutionary characteristics of O. tsutsugamushi, coupled with extensive intra-specific LGT, expand our present understanding of rapid bacterial evolutionary adaptive mechanisms.