Search

The polysaccharide (PS) capsule is essential for immune evasion and virulence of Streptococcus pneumoniae. Existing pneumococcal vaccines are designed to elicit anticapsule antibodies; however, the effectiveness of these vaccines is being challenged by the emergence of new capsule types or variants. Herein, we characterize a newly discovered capsule type, 33E, that appears to have repeatedly emerged from vaccine type 33F via an inactivation mutation in the capsule glycosyltransferase gene, wciE. Structural analysis demonstrated that 33E and 33F share an identical repeat unit backbone [→5)-β-D-Galf2Ac-(1→3)-β-D-Galp-(1→3)-α-D-Galp-(1→3)-β-D-Galf-(1→3)-β-D-Glcp-(1→], except that a galactose (α-D-Galp) branch is present in 33F but not in 33E. Though the two capsule types were indistinguishable using conventional typing methods, the monoclonal antibody Hyp33FM1 selectively bound 33F but not 33E pneumococci. Further, we confirmed that wciE encodes a glycosyltransferase that catalyzes the addition of the branching α-D-Galp and that its inactivation in 33F strains results in the expression of the 33E capsule type. Though 33F and 33E share a structural and antigenic similarity, our pilot study suggested that immunization with a 23-valent pneumococcal PS vaccine containing 33F PS did not significantly elicit cross-opsonic antibodies to 33E. New conjugate vaccines that target capsule type 33F may not necessarily protect against 33E. Therefore, studies of new conjugate vaccines require knowledge of the newly identified capsule type 33E and reliable pneumococcal typing methods capable of distinguishing it from 33F., Competing Interests: Conflict of interest U. A. B. has Intellectual Property rights on some reagents used in the study. F. A. G., J. S. S., J. J. C., and M. H. N. are U. A. B. employees. The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Brucellosis remains one of the most significant zoonotic diseases globally, responsible for both considerable human morbidity and economic losses due to its impacts on livestock productivity. Despite this, there remain significant evidence gaps in many low- and middle-income countries, including those of sub-Saharan Africa. Here we report the first molecular characterisation of Brucella sp. from Ethiopia. Fifteen Brucella sp. isolates from an outbreak in cattle from a herd in central Ethiopia were identified as Brucella abortus , using bacterial culture and molecular methods. Sequencing of the Ethiopian B. abortus isolates allowed their phylogenetic comparison with 411 B. abortus strains of diverse geographical origins, using whole genome single nucleotide polymorphisms (wgSNP). The Ethiopian isolates belonged to an early-branching lineage (Lineage A) previously only represented by data from two strains, both of sub-Saharan African origin (Kenya and Mozambique). A second B. abortus lineage (Lineage B), also comprised solely of strains originating from sub-Saharan Africa, was identified. The majority of strains belonged to one of two lineages of strains originating from a much broader geographical range. Further analyses based on multi-locus sequence typing (MLST) and multi-locus variable-number tandem repeat analysis (MLVA) expanded the number of B. abortus strains available for comparison with the Ethiopian isolates and were consistent with the findings from wgSNP analysis. MLST profiles of the Ethiopian isolates expanded the sequence type (ST) diversity of the early branching lineage of B. abortus , equivalent to wgSNP Lineage A. A more diverse cluster of STs, equivalent to wgSNP Lineage B, was comprised solely of strains originating from sub-Saharan Africa. Similarly, analysis of B. abortus MLVA profiles ( n  = 1891) confirmed that the Ethiopian isolates formed a unique cluster, similar to only two existing strains, and distinct from the majority of other strains of sub-Saharan African origin. These findings expand the known diversity of an under-represented lineage of B. abortus and suggest a potential evolutionary origin for the species in East Africa. In addition to providing information concerning Brucella species extant within Ethiopia this work serves as the basis for further studies on the global population structure and evolutionary history of a major zoonotic pathogen., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Edao, Ameni, Berg, Tekle, Whatmore, Wood, van Tonder and Ashford.)

Background: Non-tuberculous mycobacteria (NTM) are environmental microorganisms and opportunistic pathogens in individuals with pre-existing lung conditions such as cystic fibrosis (CF) and non-CF bronchiectasis. While recent studies of Mycobacterium abscessus have identified transmission within single CF centres as well as nationally and globally, transmission of other NTM species is less well studied., Methods: To investigate the potential for transmission of the Mycobacterium avium complex (MAC) we sequenced 996 isolates from 354 CF and non-CF patients at the Royal Brompton Hospital (London, UK; collected 2013-2016) and analysed them in a global context. Epidemiological links were identified from patient records. Previously published genomes were used to characterise global population structures., Results: We identified putative transmission clusters in three MAC species, although few epidemiological links could be identified. For M. avium , lineages were largely limited to single countries, while for Mycobacterium chimaera , global transmission clusters previously associated with heater-cooler units (HCUs) were found. However, the immediate ancestor of the lineage causing the major HCU-associated outbreak was a lineage already circulating in patients., Conclusions: CF and non-CF patients shared transmission chains, although the lack of epidemiological links suggested that most transmission is indirect and may involve environmental intermediates or asymptomatic carriage in the wider population., Competing Interests: Conflict of interest: J. Parkhill received consulting fees from and holds stock in Next Gen Diagnostics LLC. M.R. Loebinger received consulting fees from Savara, AstraZeneca, Insmed, Grifols, Zambon, Armarta and 30T, and received honoraria from Insmed and Grifols. No other authors have potential competing interests., (Copyright ©The authors 2023.)

Streptococcus pneumoniae can produce a wide breadth of antigenically diverse capsule types, a fact that poses a looming threat to the success of vaccines that target pneumococcal polysaccharide (PS) capsule. Yet, many pneumococcal capsule types remain undiscovered and/or uncharacterized. Prior sequence analysis of pneumococcal capsule synthesis ( cps ) loci suggested the existence of capsule subtypes among isolates identified as "serotype 36" according to conventional capsule typing methods. We discovered these subtypes represent two antigenically similar but distinguishable pneumococcal capsule serotypes, 36A and 36B. Biochemical analysis of their capsule PS structure reveals that both have the shared repeat unit backbone [→5)-α-d-Gal f -(1→1)-d-Rib-ol-(5→P→6)-β-d-Man p NAc-(1→4)-β-d-Glc p -(1→] with two branching structures. Both serotypes have a β-d-Gal p branch to Ribitol. Serotypes 36A and 36B differ by the presence of a α-d-Glc p -(1→3)-β-d-Man p NAc or α-d-Gal p -(1→3)-β-d-Man p NAc branch, respectively. Comparison of the phylogenetically distant serogroup 9 and 36  cps loci, which all encode this distinguishing glycosidic bond, revealed that the incorporation of Glc p (in types 9N and 36A) versus Gal p (in types 9A, 9V, 9L, and 36B) is associated with the identity of four amino acids in the cps- encoded glycosyltransferase WcjA. Identifying functional determinants of cps -encoded enzymes and their impact on capsule PS structure is key to improving the resolution and reliability of sequencing-based capsule typing methods and discovering novel capsule variants indistinguishable by conventional serotyping methods.

Aggregation of children in schools has been established to be a key driver of transmission of infectious diseases. Mathematical models of transmission used to predict the impact of control measures, such as vaccination and testing, commonly depend on self-reported contact data. However, the link between self-reported social contacts and pathogen transmission has not been well described. To address this, we used Staphylococcus aureus as a model organism to track transmission within two secondary schools in England and test for associations between self-reported social contacts, test positivity and the bacterial strain collected from the same students. Students filled out a social contact survey and their S. aureus colonization status was ascertained through self-administered swabs from which isolates were sequenced. Isolates from the local community were also sequenced to assess the representativeness of school isolates. A low frequency of genome-linked transmission precluded a formal analysis of links between genomic and social networks, suggesting that S. aureus transmission within schools is too rare to make it a viable tool for this purpose. Whilst we found no evidence that schools are an important route of transmission, increased colonization rates found within schools imply that school-age children may be an important source of community transmission.

The spread of carbapenemase-producing Enterobacterales (CPE) is of major public health concern. The transmission dynamics of CPE in hospitals, particularly at the national level, are not well understood. Here, we describe a retrospective nationwide genomic surveillance study of CPE in Ireland between 2012 and 2017. We sequenced 746 national surveillance CPE samples obtained between 2012 and 2017. After clustering the sequences, we used thresholds based on pairwise SNPs, and reported within-host diversity along with epidemiological data to infer recent putative transmissions. All clusters in circulating clones, derived from high-resolution phylogenies, of a species ( Klebsiella pneumoniae , Escherichia coli , Klebsiella oxytoca , Enterobacter cloacae , Enterobacter hormaechei and Citrobacter freundii ) were individually examined for evidence of transmission. Antimicrobial resistance trends over time were also assessed. We identified 352 putative transmission events in six species including widespread and frequent transmissions in three species. We detected putative outbreaks in 4/6 species with three hospitals experiencing prolonged outbreaks. The bla OXA-48 gene was the main cause of carbapenem resistance in Ireland in almost all species. An expansion in the number of sequence types carrying bla OXA-48 was an additional cause of the increasing prevalence of carbapenemase-producing K. pneumoniae and E. coli .

Summary: 10.6% patients were CRE positive. Only 27% patients were prescribed at least 1 antibiotic to which infecting pathogen was susceptible. Burn and ICU admission and antibiotics exposures facilitate CRE acquisition. Escherichia coli ST167 was the dominant CRE clone., Background: Given the high prevalence of multidrug resistance (MDR) across South Asian (SA) hospitals, we documented the epidemiology of carbapenem-resistant Enterobacterales (CRE) infections at Dhaka Medical College Hospital between October 2016 and September 2017., Methods: We enrolled patients and collected epidemiology and outcome data. All Enterobacterales were characterized phenotypically and by whole-genome sequencing. Risk assessment for the patients with CRE was performed compared with patients with carbapenem-susceptible Enterobacterales (CSE)., Results: 10.6% of all 1831 patients with a clinical specimen collected had CRE. In-hospital 30-day mortality was significantly higher with CRE [50/180 (27.8%)] than CSE [42/312 (13.5%)] (P = .001); however, for bloodstream infections, this was nonsignificant. Of 643 Enterobacterales isolated, 210 were CRE; blaNDM was present in 180 isolates, blaOXA-232 in 26, blaOXA-181 in 24, and blaKPC-2 in 5. Despite this, ceftriaxone was the most commonly prescribed empirical antibiotic and only 27% of patients were prescribed at least 1 antibiotic to which their infecting pathogen was susceptible. Significant risk factors for CRE isolation included burns unit and intensive care unit admission, and prior exposure to levofloxacin, amikacin, clindamycin, and meropenem. Escherichia coli ST167 was the dominant CRE clone. Clustering suggested clonal transmission of Klebsiella pneumoniae ST15 and the MDR hypervirulent clone, ST23. The major trajectories involved in horizontal gene transfer were IncFII and IncX3, IS26, and Tn3., Conclusions: This is the largest study from an SA public hospital combining outcome, microbiology, and genomics. The findings indicate the urgent implementation of targeted diagnostics, appropriate antibiotic use, and infection-control interventions in SA public institutions., Competing Interests: Potential conflicts of interest. The authors: No reported conflicts of interest. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest., (© The Author(s) 2022. Published by Oxford University Press on behalf of the Infectious Diseases Society of America.)

Mycobacterium bovis (M. bovis) is a causative agent of bovine tuberculosis, a significant source of morbidity and mortality in the global cattle industry. The Randomised Badger Culling Trial was a field experiment carried out between 1998 and 2005 in the South West of England. As part of this trial, M. bovis isolates were collected from contemporaneous and overlapping populations of badgers and cattle within ten defined trial areas. We combined whole genome sequences from 1,442 isolates with location and cattle movement data, identifying transmission clusters and inferred rates and routes of transmission of M. bovis. Most trial areas contained a single transmission cluster that had been established shortly before sampling, often contemporaneous with the expansion of bovine tuberculosis in the 1980s. The estimated rate of transmission from badger to cattle was approximately two times higher than from cattle to badger, and the rate of within-species transmission considerably exceeded these for both species. We identified long distance transmission events linked to cattle movement, recurrence of herd breakdown by infection within the same transmission clusters and superspreader events driven by cattle but not badgers. Overall, our data suggests that the transmission clusters in different parts of South West England that are still evident today were established by long-distance seeding events involving cattle movement, not by recrudescence from a long-established wildlife reservoir. Clusters are maintained primarily by within-species transmission, with less frequent spill-over both from badger to cattle and cattle to badger., Competing Interests: The authors have declared that no competing interests exist.

Viruses often encode proteins that mimic host proteins in order to facilitate infection. Little work has been done to understand the potential mimicry of the SARS-CoV-2, SARS-CoV, and MERS-CoV spike proteins, particularly the receptor-binding motifs, which could be important in determining tropism and druggability of the virus. Peptide and epitope motifs have been detected on coronavirus spike proteins using sequence homology approaches; however, comparing the three-dimensional shape of the protein has been shown as more informative in predicting mimicry than sequence-based comparisons. Here, we use structural bioinformatics software to characterize potential mimicry of the three coronavirus spike protein receptor-binding motifs. We utilize sequence-independent alignment tools to compare structurally known protein models with the receptor-binding motifs and verify potential mimicked interactions with protein docking simulations. Both human and non-human proteins were returned for all three receptor-binding motifs. For example, all three were similar to several proteins containing EGF-like domains: some of which are endogenous to humans, such as thrombomodulin, and others exogenous, such as Plasmodium falciparum MSP-1. Similarity to human proteins may reveal which pathways the spike protein is co-opting, while analogous non-human proteins may indicate shared host interaction partners and overlapping antibody cross-reactivity. These findings can help guide experimental efforts to further understand potential interactions between human and coronavirus proteins., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2021 The Authors.)

Bovine tuberculosis (bTB) is endemic in cattle in Ethiopia, a country that hosts the largest national cattle herd in Africa. The intensive dairy sector, most of which is peri-urban, has the highest prevalence of disease. Previous studies in Ethiopia have demonstrated that the main cause is Mycobacterium bovis , which has been investigated using conventional molecular tools including deletion typing, spoligotyping and Mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR). Here we use whole-genome sequencing to examine the population structure of M. bovis in Ethiopia. A total of 134 M . bovis isolates were sequenced including 128 genomes from 85 mainly dairy cattle and six genomes isolated from humans, originating from 12 study sites across Ethiopia. These genomes provided a good representation of the previously described population structure of M. bovis , based on spoligotyping and demonstrated that the population is dominated by the clonal complexes African 2 (Af2) and European 3 (Eu3). A range of within-host diversity was observed amongst the isolates and evidence was found for both short- and long-distance transmission. Detailed analysis of available genomes from the Eu3 clonal complex combined with previously published genomes revealed two distinct introductions of this clonal complex into Ethiopia between 1950 and 1987, likely from Europe. This work is important to help better understand bTB transmission in cattle in Ethiopia and can potentially inform national strategies for bTB control in Ethiopia and beyond.

Mycobacterium bovis is a causal agent of bovine tuberculosis (bTB), one of the most important diseases currently facing the cattle industry worldwide. Tracing the source of M. bovis infections of livestock is an important tool for understanding the epidemiology of bTB and defining control/eradication strategies. In this study, whole genome sequencing (WGS) of 74 M . bovis isolates sourced from naturally infected cattle in the State of Rio Grande do Sul (RS), southern Brazil, was used to evaluate the population structure of M. bovis in the region, identify potential transmission events and date the introduction of clonal complex (CC) European 2 (Eu2). In silico spoligotyping identified 11 distinct patterns including four new profiles and two CCs, European 1 (Eu1) and Eu2. The analyses revealed a high level of genetic diversity in the majority of herds and identified putative transmission clusters that suggested that within- and between-herd transmission is occurring in RS. In addition, a comparison with other published M. bovis isolates from Argentina, Brazil, Paraguay and Uruguay demonstrated some evidence for a possible cross-border transmission of CC Eu1 into RS from Uruguay or Argentina. An estimated date for the introduction of CC Eu2 into RS in the middle of the 19th century correlated with the historical introduction of cattle into RS to improve existing local breeds. These findings contribute to the understanding of the population structure of M. bovis in southern Brazil and highlight the potential of WGS in surveillance and helping to identify bTB transmission.

The aim of our study was to investigate phenotypic and genotypic features of streptococci misidentified (misID) as Streptococcus pneumoniae, obtained over 20 years from hospital patients in Poland. Sixty-three isolates demonstrating microbiological features typical for pneumococci (optochin susceptibility and/or bile solubility) were investigated by phenotypic tests, lytA and 16S rRNA gene polymorphism and whole-genome sequencing (WGS). All isolates had a 6-bp deletion in the lytA 3' terminus, characteristic for Mitis streptococc and all but two isolates lacked the pneumococcal signature cytosine at nucleotide position 203 in the 16S rRNA genes. The counterparts of psaA and ply were present in 100% and 81.0% of isolates, respectively; the spn9802 and spn9828 loci were characteristic for 49.2% and 38.1% of isolates, respectively. Phylogenetic trees and networks, based on the multilocus sequence analysis (MLSA) scheme, ribosomal multilocus sequence typing (rMLST) scheme and core-genome analysis, clearly separated investigated isolates from S. pneumoniae and demonstrated the polyclonal character of misID streptococci, associated with the Streptococcus pseudopneumoniae and Streptococcus mitis groups. While the S. pseudopneumoniae clade was relatively well defined in all three analyses, only the core-genome analysis revealed the presence of another cluster comprising a fraction of misID streptococci and a strain proposed elsewhere as a representative of a novel species in the Mitis group. Our findings point to complex phylogenetic and taxonomic relationships among S. mitis-like bacteria and support the notion that this group may in fact consist of several distinct species.

Streptococcus pneumoniae (pneumococcus) is a major human pathogen producing structurally diverse capsular polysaccharides. Widespread use of highly successful pneumococcal conjugate vaccines (PCVs) targeting pneumococcal capsules has greatly reduced infections by the vaccine types but increased infections by nonvaccine serotypes. Herein, we report a new and the 100th capsule type, named serotype 10D, by determining its unique chemical structure and biosynthetic roles of all capsule synthesis locus ( cps ) genes. The name 10D reflects its serologic cross-reaction with serotype 10A and appearance of cross-opsonic antibodies in response to immunization with 10A polysaccharide in a 23-valent pneumococcal vaccine. Genetic analysis showed that 10D cps has three large regions syntenic to and highly homologous with cps loci from serotype 6C, serotype 39, and an oral streptococcus strain ( S. mitis SK145). The 10D cps region syntenic to SK145 is about 6 kb and has a short gene fragment of wciN α at the 5' end. The presence of this nonfunctional wciN α fragment provides compelling evidence for a recent interspecies genetic transfer from oral streptococcus to pneumococcus. Since oral streptococci have a large repertoire of cps loci, widespread PCV usage could facilitate the appearance of novel serotypes through interspecies recombination. IMPORTANCE The polysaccharide capsule is essential for the pathogenicity of pneumococcus, which is responsible for millions of deaths worldwide each year. Currently available pneumococcal vaccines are designed to elicit antibodies to the capsule polysaccharides of the pneumococcal isolates commonly causing diseases, and the antibodies provide protection only against the pneumococcus expressing the vaccine-targeted capsules. Since pneumococci can produce different capsule polysaccharides and therefore reduce vaccine effectiveness, it is important to track the appearance of novel pneumococcal capsule types and how these new capsules are created. Herein, we describe a new and the 100th pneumococcal capsule type with unique chemical and serological properties. The capsule type was named 10D for its serologic similarity to 10A. Genetic studies provide strong evidence that pneumococcus created 10D capsule polysaccharide by capturing a large genetic fragment from an oral streptococcus. Such interspecies genetic exchanges could greatly increase diversity of pneumococcal capsules and complicate serotype shifts.

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 dissemination 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.

Objectives: We reported tet(S/M) in Streptococcus pneumoniae and investigated its temporal spread in relation to nationwide clinical interventions., Methods: We whole-genome sequenced 12 254 pneumococcal isolates from 29 countries on an Illumina HiSeq sequencer. Serotype, multilocus ST and antibiotic resistance were inferred from genomes. An SNP tree was built using Gubbins. Temporal spread was reconstructed using a birth-death model., Results: We identified tet(S/M) in 131 pneumococcal isolates and none carried other known tet genes. Tetracycline susceptibility testing results were available for 121 tet(S/M)-positive isolates and all were resistant. A majority (74%) of tet(S/M)-positive isolates were from South Africa and caused invasive diseases among young children (59% HIV positive, where HIV status was available). All but two tet(S/M)-positive isolates belonged to clonal complex (CC) 230. A global phylogeny of CC230 (n=389) revealed that tet(S/M)-positive isolates formed a sublineage predicted to exhibit resistance to penicillin, co-trimoxazole, erythromycin and tetracycline. The birth-death model detected an unrecognized outbreak of this sublineage in South Africa between 2000 and 2004 with expected secondary infections (effective reproductive number, R) of ∼2.5. R declined to ∼1.0 in 2005 and <1.0 in 2012. The declining epidemic could be related to improved access to ART in 2004 and introduction of pneumococcal conjugate vaccine (PCV) in 2009. Capsular switching from vaccine serotype 14 to non-vaccine serotype 23A was observed within the sublineage., Conclusions: The prevalence of tet(S/M) in pneumococci was low and its dissemination was due to an unrecognized outbreak of CC230 in South Africa. Capsular switching in this MDR sublineage highlighted its potential to continue to cause disease in the post-PCV13 era., (© The Author(s) 2019. Published by Oxford University Press on behalf of the British Society for Antimicrobial Chemotherapy.)

Background: Invasive pneumococcal disease remains an important health priority owing to increasing disease incidence caused by pneumococci expressing non-vaccine serotypes. We previously defined 621 Global Pneumococcal Sequence Clusters (GPSCs) by analysing 20 027 pneumococcal isolates collected worldwide and from previously published genomic data. In this study, we aimed to investigate the pneumococcal lineages behind the predominant serotypes, the mechanism of serotype replacement in disease, as well as the major pneumococcal lineages contributing to invasive pneumococcal disease in the post-vaccine era and their antibiotic resistant traits., Methods: We whole-genome sequenced 3233 invasive pneumococcal disease isolates from laboratory-based surveillance programmes in Hong Kong (n=78), Israel (n=701), Malawi (n=226), South Africa (n=1351), The Gambia (n=203), and the USA (n=674). The genomes represented pneumococci from before and after pneumococcal conjugate vaccine (PCV) introductions and were from children younger than 3 years. We identified predominant serotypes by prevalence and their major contributing lineages in each country, and assessed any serotype replacement by comparing the incidence rate between the pre-PCV and PCV periods for Israel, South Africa, and the USA. We defined the status of a lineage as vaccine-type GPSC (≥50% 13-valent PCV [PCV13] serotypes) or non-vaccine-type GPSC (>50% non-PCV13 serotypes) on the basis of its initial serotype composition detected in the earliest vaccine period to measure their individual contribution toward serotype replacement in each country. Major pneumococcal lineages in the PCV period were identified by pooled incidence rate using a random effects model., Findings: The five most prevalent serotypes in the PCV13 period varied between countries, with only serotypes 5, 12F, 15B/C, 19A, 33F, and 35B/D common to two or more countries. The five most prevalent serotypes in the PCV13 period varied between countries, with only serotypes 5, 12F, 15B/C, 19A, 33F, and 35B/D common to two or more countries. These serotypes were associated with more than one lineage, except for serotype 5 (GPSC8). Serotype replacement was mainly mediated by expansion of non-vaccine serotypes within vaccine-type GPSCs and, to a lesser extent, by increases in non-vaccine-type GPSCs. A globally spreading lineage, GPSC3, expressing invasive serotypes 8 in South Africa and 33F in the USA and Israel, was the most common lineage causing non-vaccine serotype invasive pneumococcal disease in the PCV13 period. We observed that same prevalent non-vaccine serotypes could be associated with distinctive lineages in different countries, which exhibited dissimilar antibiotic resistance profiles. In non-vaccine serotype isolates, we detected significant increases in the prevalence of resistance to penicillin (52 [21%] of 249 vs 169 [29%] of 575, p=0·0016) and erythromycin (three [1%] of 249 vs 65 [11%] of 575, p=0·0031) in the PCV13 period compared with the pre-PCV period., Interpretation: Globally spreading lineages expressing invasive serotypes have an important role in serotype replacement, and emerging non-vaccine serotypes associated with different pneumococcal lineages in different countries might be explained by local antibiotic-selective pressures. Continued genomic surveillance of the dynamics of the pneumococcal population with increased geographical representation in the post-vaccine period will generate further knowledge for optimising future vaccine design., Funding: Bill & Melinda Gates Foundation, Wellcome Sanger Institute, and the US Centers for Disease Control., (Copyright © 2019 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY 4.0 license. Published by Elsevier Ltd.. All rights reserved.)

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.

Background: Pneumococcal conjugate vaccines have reduced the incidence of invasive pneumococcal disease, caused by vaccine serotypes, but non-vaccine-serotypes remain a concern. We used whole genome sequencing to study pneumococcal serotype, antibiotic resistance and invasiveness, in the context of genetic background., Methods: Our dataset of 13,454 genomes, combined with four published genomic datasets, represented Africa (40%), Asia (25%), Europe (19%), North America (12%), and South America (5%). These 20,027 pneumococcal genomes were clustered into lineages using PopPUNK, and named Global Pneumococcal Sequence Clusters (GPSCs). From our dataset, we additionally derived serotype and sequence type, and predicted antibiotic sensitivity. We then measured invasiveness using odds ratios that relating prevalence in invasive pneumococcal disease to carriage., Findings: The combined collections (n = 20,027) were clustered into 621 GPSCs. Thirty-five GPSCs observed in our dataset were represented by >100 isolates, and subsequently classed as dominant-GPSCs. In 22/35 (63%) of dominant-GPSCs both non-vaccine serotypes and vaccine serotypes were observed in the years up until, and including, the first year of pneumococcal conjugate vaccine introduction. Penicillin and multidrug resistance were higher (p < .05) in a subset dominant-GPSCs (14/35, 9/35 respectively), and resistance to an increasing number of antibiotic classes was associated with increased recombination (R 2  = 0.27 p < .0001). In 28/35 dominant-GPSCs, the country of isolation was a significant predictor (p < .05) of its antibiogram (mean misclassification error 0.28, SD ± 0.13). We detected increased invasiveness of six genetic backgrounds, when compared to other genetic backgrounds expressing the same serotype. Up to 1.6-fold changes in invasiveness odds ratio were observed., Interpretation: We define GPSCs that can be assigned to any pneumococcal genomic dataset, to aid international comparisons. Existing non-vaccine-serotypes in most GPSCs preclude the removal of these lineages by pneumococcal conjugate vaccines; leaving potential for serotype replacement. A subset of GPSCs have increased resistance, and/or serotype-independent invasiveness., (Copyright © 2019 The Authors. Published by Elsevier B.V. All rights reserved.)

The introduction of pneumococcal conjugate vaccines (PCVs) into childhood vaccination programs has reduced carriage of vaccine serotypes and pneumococcal disease. The 10-valent PCV was introduced in Iceland in 2011. The aim of this study was to determine PCV impact on the prevalence of serotypes, genetic lineages, and antimicrobial-resistant pneumococci isolated from the lower respiratory tract (LRT) of adults. Pneumococci isolated between 2009 and 2017 at the Landspitali University Hospital were included ( n  = 797). The hospital serves almost three-quarters of the Icelandic population. Isolates were serotyped and tested for antimicrobial susceptibility, and the genome of every other isolate collected between 2009 and 2014 was sequenced ( n  = 275). Serotypes and multilocus sequence types (STs) were extracted from the genome data. Three study periods were defined, 2009 to 2011 (PreVac), 2012 to 2014 (PostVac-I), and 2015 to 2017 (PostVac-II). The total number of isolates and vaccine-type (VT) pneumococci decreased from PreVac to PostVac-II ( n  = 314 versus n  = 230 [ p  = 0.002] and n  = 170 versus n  = 33 [ p  < 0.001], respectively), but non-vaccine-type (NVT) pneumococci increased among adults 18 to 64 years old ( n  = 56 versus n  = 114 [ p  = 0.008]). Serotype 19F decreased in the PostVac-II period; these isolates were all multidrug resistant (MDR) and were members of the Taiwan 19F -14 PMEN lineage. Serotype 6A decreased among adults ≥65 years old in the PostVac-II period ( p  = 0.037), while serotype 6C increased ( p  = 0.021) and most serotype 6C isolates were MDR. Nonencapsulated Streptococcus pneumoniae (NESp) isolates increased among adults 18 to 64 years old in the PostVac-II period, and the majority were MDR ( p  = 0.028). An overall reduction in the number of LRT samples and pneumococcus-positive cultures and significant changes in the serotype distribution became evident within 4 years, thereby demonstrating a significant herd effect., (Copyright © 2019 Quirk et al.)

Background: Bovine tuberculosis (bTB) caused by Mycobacterium bovis is a re-emerging problem in both livestock and humans. The association of some M. bovis strains with hyper-virulence, MDR-TB and disseminated disease makes it imperative to understand the biology of the pathogen., Methods: Mycobacterium bovis (15) among 1755 M. tuberculosis complex (MTBC) isolated between 2012 and 2014 were characterized and analyzed for associated patient demography and other risk factors. Five of the M. bovis isolates were whole-genome sequenced and comparatively analyzed against a global collection of published M. bovis genomes., Results: Mycobacterium bovis was isolated from 3/560(0.5%) females and 12/1195(1.0%) males with pulmonary TB. The average age of M. bovis infected cases was 46.8 years (7-72years). TB patients from the Northern region of Ghana (1.9%;4/212) had a higher rate of infection with M. bovis (OR = 2.7,p = 0.0968) compared to those from the Greater Accra region (0.7%;11/1543). Among TB patients with available HIV status, the odds of isolating M. bovis from HIV patients (2/119) was 3.3 higher relative to non-HIV patients (4/774). Direct contact with livestock or their unpasteurized products was significantly associated with bTB (p<0.0001, OR = 124.4,95% CI = 30.1-508.3). Two (13.3%) of the M. bovis isolates were INH resistant due to the S315T mutation in katG whereas one (6.7%) was RIF resistant with Q432P and I1491S mutations in rpoB. M. bovis from Ghana resolved as mono-phyletic branch among mostly M. bovis from Africa irrespective of the host and were closest to the root of the global M. bovis phylogeny. M. bovis-specific amino acid mutations were detected among MTBC core genes such as mce1A, mmpL1, pks6, phoT, pstB, glgP and Rv2955c. Additional mutations P6T in chaA, G187E in mgtC, T35A in Rv1979c, S387A in narK1, L400F in fas and A563T in eccA1 were restricted to the 5 clinical M. bovis from Ghana., Conclusion: Our data indicate potential zoonotic transmission of bTB in Ghana and hence calls for intensified public education on bTB, especially among risk groups., Competing Interests: The authors have declared that no competing interests exist.

Understanding the structure of a bacterial population is essential in order to understand bacterial evolution. Estimating the core genome (those genes common to all, or nearly all, strains of a species) is a key component of such analyses. The size and composition of the core genome varies by dataset, but we hypothesized that the variation between different collections of the same bacterial species would be minimal. To investigate this, we analyzed the genome sequences of 3,118 pneumococci recovered from healthy individuals in Reykjavik (Iceland), Southampton (United Kingdom), Boston (United States), and Maela (Thailand). The analyses revealed a "supercore" genome (genes shared by all 3,118 pneumococci) of 558 genes, although an additional 354 core genes were shared by pneumococci from Reykjavik, Southampton, and Boston. Overall, the size and composition of the core and pan-genomes among pneumococci recovered in Reykjavik, Southampton, and Boston were similar. Maela pneumococci were distinctly different in that they had a smaller core genome and larger pan-genome. The pan-genome of Maela pneumococci contained several >25 Kb sequence regions (flanked by pneumococcal genes) that were homologous to genomic regions found in other bacterial species. Overall, our work revealed that some subsets of the global pneumococcal population are highly heterogeneous, and our hypothesis was rejected. This is an important finding in terms of understanding genetic variation among pneumococci and is also an essential point of consideration before generalizing the findings from a single dataset to the wider pneumococcal population.

Vaccination with pneumococcal conjugate vaccines (PCVs) disrupts the pneumococcal population. Our aim was to determine the impact of the 10-valent PCV on the serotypes, genetic lineages, and antimicrobial susceptibility of pneumococci isolated from children in Iceland. Pneumococci were collected between 2009 and 2017 from the nasopharynges of healthy children attending 15 day care centers and from the middle ears (MEs) of children with acute otitis media from the greater Reykjavik capital area. Isolates were serotyped and tested for antimicrobial susceptibility. Whole-genome sequencing (WGS) was performed on alternate isolates from 2009 to 2014, and serotypes and multilocus sequence types (STs) were extracted from the WGS data. Two study periods were defined: 2009 to 2011 (PreVac) and 2012 to 2017 (PostVac). The overall nasopharyngeal carriage rate was similar between the two periods (67.3% PreVac and 61.5% PostVac, P = 0.090). Vaccine-type (VT) pneumococci decreased and nonvaccine-type (NVT) pneumococci (serotypes 6C, 15A, 15B/C, 21, 22F, 23A, 23B, 35F, and 35B) significantly increased in different age strata post-PCV introduction. The total number of pneumococci recovered from ME samples significantly decreased as did the proportion that were VTs, although NVT pneumococci (6C, 15B/C, 23A, and 23B) increased significantly. Most serotype 6C pneumococci were multidrug resistant (MDR). Serotype 19F was the predominant serotype associated with MEs, and it significantly decreased post-PCV introduction: these isolates were predominantly MDR and of the Taiwan 19F -14 PMEN lineage. Overall, the nasopharyngeal carriage rate remained constant and the number of ME-associated pneumococci decreased significantly post-PCV introduction; however, there was a concomitant and statistically significant shift from VTs to NVTs in both collections of pneumococci., (Copyright © 2018 Quirk et al.)

Reactive oxygen species (ROS)-mediated oxidative stress and DNA damage have recently been recognized as contributing to the efficacy of most bactericidal antibiotics, irrespective of their primary macromolecular targets. Inhibitors of targets involved in both combating oxidative stress as well as being required for in vivo survival may exhibit powerful synergistic action. This study demonstrates that the de novo arginine biosynthetic pathway in Mycobacterium tuberculosis ( Mtb ) is up-regulated in the early response to the oxidative stress-elevating agent isoniazid or vitamin C. Arginine deprivation rapidly sterilizes the Mtb de novo arginine biosynthesis pathway mutants Δ argB and Δ argF without the emergence of suppressor mutants in vitro as well as in vivo. Transcriptomic and flow cytometry studies of arginine-deprived Mtb have indicated accumulation of ROS and extensive DNA damage. Metabolomics studies following arginine deprivation have revealed that these cells experienced depletion of antioxidant thiols and accumulation of the upstream metabolite substrate of ArgB or ArgF enzymes. Δ argB and Δ argF were unable to scavenge host arginine and were quickly cleared from both immunocompetent and immunocompromised mice. In summary, our investigation revealed in vivo essentiality of the de novo arginine biosynthesis pathway for Mtb and a promising drug target space for combating tuberculosis., Competing Interests: The authors declare no conflict of interest.

Competition among bacterial members of the same ecological niche is mediated by bacteriocins: antimicrobial peptides produced by bacterial species to kill other bacteria. Bacteriocins are also promising candidates for novel antimicrobials. Streptococcus pneumoniae (the "pneumococcus") is a leading cause of morbidity and mortality worldwide and a frequent colonizer of the human nasopharynx. Here, 14 newly discovered bacteriocin gene clusters were identified among >6,200 pneumococcal genomes. The molecular epidemiology of the bacteriocin clusters was investigated using a large global and historical pneumococcal dataset dating from 1916. These analyses revealed extraordinary bacteriocin diversity among pneumococci and the majority of bacteriocin clusters were also found in other streptococcal species. Genomic hotspots for the integration of different bacteriocin gene clusters were discovered. Experimentally, bacteriocin genes were transcriptionally active when the pneumococcus was under stress and when two strains were co-cultured in broth. These findings reveal much more diversity among bacterial defense mechanisms than previously appreciated, which fundamentally broaden our understanding of bacteriocins relative to intraspecies and interspecies nasopharyngeal competition and bacterial population structure.

Streptococcus pneumoniae is responsible for 240 000-460 000 deaths in children under 5 years of age each year. Accurate identification of pneumococcal serotypes is important for tracking the distribution and evolution of serotypes following the introduction of effective vaccines. Recent efforts have been made to infer serotypes directly from genomic data but current software approaches are limited and do not scale well. Here, we introduce a novel method, SeroBA, which uses a k-mer approach. We compare SeroBA against real and simulated data and present results on the concordance and computational performance against a validation dataset, the robustness and scalability when analysing a large dataset, and the impact of varying the depth of coverage on sequence-based serotyping. SeroBA can predict serotypes, by identifying the cps locus, directly from raw whole genome sequencing read data with 98 % concordance using a k-mer-based method, can process 10 000 samples in just over 1 day using a standard server and can call serotypes at a coverage as low as 15-21×. SeroBA is implemented in Python3 and is freely available under an open source GPLv3 licence from: https://github.com/sanger-pathogens/seroba.

A newly recognized pneumococcal serotype, 35D, which differs from the 35B polysaccharide in structure and serology by not binding to factor serum 35a, was recently reported. The genetic basis for this distinctive serology is due to the presence of an inactivating mutation in wciG , which encodes an O-acetyltransferase responsible for O-acetylation of a galactofuranose. Here, we assessed the genomic data of a worldwide pneumococcal collection to identify serotype 35D isolates and understand their geographical distribution, genetic background, and invasiveness potential. Of 21,980 pneumococcal isolates, 444 were originally typed as serotype 35B by PneumoCaT. Analysis of the wciG gene revealed 23 isolates from carriage ( n = 4) and disease ( n = 19) with partial or complete loss-of-function mutations, including mutations resulting in premature stop codons ( n = 22) and an in-frame mutation ( n = 1). These were selected for further analysis. The putative 35D isolates were geographically widespread, and 65.2% (15/23) of them was recovered after the introduction of pneumococcal conjugate vaccine 13 (PCV13). Compared with serotype 35B isolates, putative serotype 35D isolates have higher invasive disease potentials based on odds ratios (OR) (11.58; 95% confidence interval[CI], 1.42 to 94.19 versus 0.61; 95% CI, 0.40 to 0.92) and a higher prevalence of macrolide resistance mediated by mefA (26.1% versus 7.6%; P = 0.009). Using the Quellung reaction, 50% (10/20) of viable isolates were identified as serotype 35D, 25% (5/20) as serotype 35B, and 25% (5/20) as a mixture of 35B/35D. The discrepancy between phenotype and genotype requires further investigation. These findings illustrated a global distribution of an invasive serotype, 35D, among young children post-PCV13 introduction and underlined the invasive potential conferred by the loss of O-acetylation in the pneumococcal capsule., (Copyright © 2018 Lo et al.)

Mucosal-associated invariant T (MAIT) cells represent an innate T-cell population that can recognize ligands generated by the microbial riboflavin synthesis pathway, presented via the major histocompatibility complex class I-related molecule (MR1). Streptococcus pneumoniae is a major human pathogen that is also associated with commensal carriage; thus, host control at the mucosal interface is critical. The recognition of pneumococci by MAIT cells has not been defined nor have the genomics and transcriptomics of the riboflavin operon. We observed robust recognition of pneumococci by MAIT cells, using both MR1-dependent and MR1-independent pathways. The pathway used was dependent on the antigen-presenting cell. The riboflavin operon was highly conserved across a range of 571 pneumococci from 39 countries, dating back to 1916, and different versions of the riboflavin operon were also identified in related Streptococcus species. These data indicate an important functional relationship between MAIT cells and pneumococci.

Bacteriophages (phages) infect many bacterial species, but little is known about the diversity of phages among the pneumococcus, a leading global pathogen. The objectives of this study were to determine the prevalence, diversity and molecular epidemiology of prophages (phage DNA integrated within the bacterial genome) among pneumococci isolated over the past 90 years. Nearly 500 pneumococcal genomes were investigated and RNA sequencing was used to explore prophage gene expression. We revealed that every pneumococcal genome contained prophage DNA. 286 full-length/putatively full-length pneumococcal prophages were identified, of which 163 have not previously been reported. Full-length prophages clustered into four major groups and every group dated from the 1930-40 s onward. There was limited evidence for genes shared between prophage clusters. Prophages typically integrated in one of five different sites within the pneumococcal genome. 72% of prophages possessed the virulence genes pblA and/or pblB. Individual prophages and the host pneumococcal genetic lineage were strongly associated and some prophages persisted for many decades. RNA sequencing provided clear evidence of prophage gene expression. Overall, pneumococcal prophages were highly prevalent, demonstrated a structured population, possessed genes associated with virulence, and were expressed under experimental conditions. Pneumococcal prophages are likely to play a more important role in pneumococcal biology and evolution than previously recognised.

The pneumococcus is a leading global pathogen and a key virulence factor possessed by the majority of pneumococci is an antigenic polysaccharide capsule ('serotype'), which is encoded by the capsular ( cps ) locus. Approximately 100 different serotypes are known, but the extent of sequence diversity within the cps loci of individual serotypes is not well understood. Investigating serotype-specific sequence variation is crucial to the design of sequence-based serotyping methodology, understanding pneumococcal conjugate vaccine (PCV) effectiveness and the design of future PCVs. The availability of large genome datasets makes it possible to assess population-level variation among pneumococcal serotypes and in this study 5405 pneumococcal genomes were used to investigate cps locus diversity among 49 different serotypes. Pneumococci had been recovered between 1916 and 2014 from people of all ages living in 51 countries. Serotypes were deduced bioinformatically, cps locus sequences were extracted and variation was assessed within the cps locus, in the context of pneumococcal genetic lineages. Overall, cps locus sequence diversity varied markedly: low to moderate diversity was revealed among serogroups/types 1, 3, 7, 9, 11 and 22; whereas serogroups/types 6, 19, 23, 14, 15, 18, 33 and 35 displayed high diversity. Putative novel and/or hybrid cps loci were identified among all serogroups/types apart from 1, 3 and 9. This study demonstrated that cps locus sequence diversity varied widely between serogroups/types. Investigation of the biochemical structure of the polysaccharide capsule of major variants, particularly PCV-related serotypes and those that appear to be novel or hybrids, is warranted.

Background: One of the most important global pathogens infecting all age groups is Streptococcus pneumoniae (the 'pneumococcus'). Pneumococci reside in the paediatric nasopharynx, where they compete for space and resources, and one competition strategy is to produce a bacteriocin (antimicrobial peptide or protein) to attack other bacteria and an immunity protein to protect against self-destruction. We analysed a collection of 336 diverse pneumococcal genomes dating from 1916 onwards, identified bacteriocin cassettes, detailed their genetic composition and sequence diversity, and evaluated the data in the context of the pneumococcal population structure., Results: We found that all genomes maintained a blp bacteriocin cassette and we identified several novel blp cassettes and genes. The composition of the 'bacteriocin/immunity region' of the blp cassette was highly variable: one cassette possessed six bacteriocin genes and eight putative immunity genes, whereas another cassette had only one of each. Both widely-distributed and highly clonal blp cassettes were identified. Most surprisingly, one-third of pneumococcal genomes also possessed a cassette encoding a novel circular bacteriocin that we called pneumocyclicin, which shared a similar genetic organisation to well-characterised circular bacteriocin cassettes in other bacterial species. Pneumocyclicin cassettes were mainly of one genetic cluster and largely found among seven major pneumococcal clonal complexes., Conclusions: These detailed genomic analyses revealed a novel pneumocyclicin cassette and a wide variety of blp bacteriocin cassettes, suggesting that competition in the nasopharynx is a complex biological phenomenon.

The pneumococcus is a leading pathogen infecting children and adults. Safe, effective vaccines exist, and they work by inducing antibodies to the polysaccharide capsule (unique for each serotype) that surrounds the cell; however, current vaccines are limited by the fact that only a few of the nearly 100 antigenically distinct serotypes are included in the formulations. Within the serotypes, serogroup 6 pneumococci are a frequent cause of serious disease and common colonizers of the nasopharynx in children. Serotype 6E was first reported in 2004 but was thought to be rare; however, we and others have detected serotype 6E among recent pneumococcal collections. Therefore, we analyzed a diverse data set of ∼1,000 serogroup 6 genomes, assessed the prevalence and distribution of serotype 6E, analyzed the genetic diversity among serogroup 6 pneumococci, and investigated whether pneumococcal conjugate vaccine-induced serotype 6A and 6B antibodies mediate the killing of serotype 6E pneumococci. We found that 43% of all genomes were of serotype 6E, and they were recovered worldwide from healthy children and patients of all ages with pneumococcal disease. Four genetic lineages, three of which were multidrug resistant, described ∼90% of the serotype 6E pneumococci. Serological assays demonstrated that vaccine-induced serotype 6B antibodies were able to elicit killing of serotype 6E pneumococci. We also revealed three major genetic clusters of serotype 6A capsular sequences, discovered a new hybrid 6C/6E serotype, and identified 44 examples of serotype switching. Therefore, while vaccines appear to offer protection against serotype 6E, genetic variants may reduce vaccine efficacy in the longer term because of the emergence of serotypes that can evade vaccine-induced immunity., (Copyright © 2015, van Tonder et al.)

The bacterial core genome is of intense interest and the volume of whole genome sequence data in the public domain available to investigate it has increased dramatically. The aim of our study was to develop a model to estimate the bacterial core genome from next-generation whole genome sequencing data and use this model to identify novel genes associated with important biological functions. Five bacterial datasets were analysed, comprising 2096 genomes in total. We developed a Bayesian decision model to estimate the number of core genes, calculated pairwise evolutionary distances (p-distances) based on nucleotide sequence diversity, and plotted the median p-distance for each core gene relative to its genome location. We designed visually-informative genome diagrams to depict areas of interest in genomes. Case studies demonstrated how the model could identify areas for further study, e.g. 25% of the core genes with higher sequence diversity in the Campylobacter jejuni and Neisseria meningitidis genomes encoded hypothetical proteins. The core gene with the highest p-distance value in C. jejuni was annotated in the reference genome as a putative hydrolase, but further work revealed that it shared sequence homology with beta-lactamase/metallo-beta-lactamases (enzymes that provide resistance to a range of broad-spectrum antibiotics) and thioredoxin reductase genes (which reduce oxidative stress and are essential for DNA replication) in other C. jejuni genomes. Our Bayesian model of estimating the core genome is principled, easy to use and can be applied to large genome datasets. This study also highlighted the lack of knowledge currently available for many core genes in bacterial genomes of significant global public health importance.