Search

Abstract Respiratory syncytial virus (RSV) is the leading cause of hospitalisation for respiratory infection in young children. RSV disease severity is known to be age-dependent and highest in young infants, but other correlates of severity, particularly the presence of additional respiratory pathogens, are less well understood. In this study, nasopharyngeal swabs were collected from two cohorts of RSV-positive infants 100 pathogens, including all common respiratory viruses and bacteria, from samples collected from 433 infants, that burden of additional viruses is common (111/433, 26%) but only modestly correlates with RSV disease severity. In contrast, there is strong evidence in both cohorts and across age groups that presence of Haemophilus bacteria (194/433, 45%) is associated with higher severity, including much higher rates of hospitalisation (odds ratio 4.25, 95% CI 2.03–9.31). There is no evidence for association between higher severity and other detected bacteria, and no difference in severity between RSV genotypes. Our findings reveal the genomic diversity of additional pathogens during RSV infection in infants, and provide an evidence base for future causal investigations of the impact of co-infection on RSV disease severity.

ABSTRACT Recombination of short DNA fragments via horizontal gene transfer (HGT) can introduce beneficial alleles, create genomic disharmony through negative epistasis, and create adaptive gene combinations through positive epistasis. For non-core (accessory) genes, the negative epistatic cost is likely to be minimal because the incoming genes have not co-evolved with the recipient genome and are frequently observed as tightly linked cassettes with major effects. By contrast, interspecific recombination in the core genome is expected to be rare because disruptive allelic replacement is likely to introduce negative epistasis. Why then is homologous recombination common in the core of bacterial genomes? To understand this enigma, we take advantage of an exceptional model system, the common enteric pathogens Campylobacter jejuni and C. coli that are known for very high magnitude interspecies gene flow in the core genome. As expected, HGT does indeed disrupt co-adapted allele pairings, indirect evidence of negative epistasis. However, multiple HGT events enable recovery of the genome’s co-adaption between introgressing alleles, even in core metabolism genes (e.g., formate dehydrogenase). These findings demonstrate that, even for complex traits, genetic coalitions can be decoupled, transferred, and independently reinstated in a new genetic background—facilitating transition between fitness peaks. In this example, the two-step recombinational process is associated with C. coli that are adapted to the agricultural niche.IMPORTANCEGenetic exchange among bacteria shapes the microbial world. From the acquisition of antimicrobial resistance genes to fundamental questions about the nature of bacterial species, this powerful evolutionary force has preoccupied scientists for decades. However, the mixing of genes between species rests on a paradox: 0n one hand, promoting adaptation by conferring novel functionality; on the other, potentially introducing disharmonious gene combinations (negative epistasis) that will be selected against. Taking an interdisciplinary approach to analyze natural populations of the enteric bacteria Campylobacter, an ideal example of long-range admixture, we demonstrate that genes can independently transfer across species boundaries and rejoin in functional networks in a recipient genome. The positive impact of two-gene interactions appears to be adaptive by expanding metabolic capacity and facilitating niche shifts through interspecific hybridization. This challenges conventional ideas and highlights the possibility of multiple-step evolution of multi-gene traits by interspecific introgression.

ABSTRACTAmong genes present in all group A streptococci (GAS), those encoding M-fibril and T-pilus proteins display the highest levels of sequence diversity, giving rise to the two primary serological typing schemes historically used to define strain. A new genotyping scheme for the pilin adhesin and backbone genes is developed and, when combined with emm typing, provides an account of the global GAS strain population. Cluster analysis based on nucleotide sequence similarity assigns most T-serotypes to discrete pilin backbone sequence clusters, yet the established T-types correspond to only half the clusters. The major pilin adhesin and backbone sequence clusters yield 98 unique combinations, defined as “pilin types.” Numerous horizontal transfer events that involve pilin or emm genes generate extensive antigenic and functional diversity on the bacterial cell surface and lead to the emergence of new strains. Inferred pilin genotypes applied to a meta-analysis of global population-based collections of pharyngitis and impetigo isolates reveal highly significant associations between pilin genotypes and GAS infection at distinct ecological niches, consistent with a role for pilin gene products in adaptive evolution. Integration of emm and pilin typing into open-access online tools (pubmlst.org) ensures broad utility for end-users wanting to determine the architecture of M-fibril and T-pilus genes from genome assemblies.IMPORTANCEPrecision in defining the variant forms of infectious agents is critical to understanding their population biology and the epidemiology of associated diseases. Group A Streptococcus (GAS) is a global pathogen that causes a wide range of diseases and displays a highly diverse cell surface due to the antigenic heterogeneity of M-fibril and T-pilus proteins which also act as virulence factors of varied functions. emm genotyping is well-established and highly utilized, but there is no counterpart for pilin genes. A global GAS collection provides the basis for a comprehensive pilin typing scheme, and online tools for determining emm and pilin genotypes are developed. Application of these tools reveals the expansion of structural-functional diversity among GAS via horizontal gene transfer, as evidenced by unique combinations of surface protein genes. Pilin and emm genotype correlations with superficial throat vs skin infection provide new insights on the molecular determinants underlying key ecological and epidemiological trends.

IntroductionNon-baumannii Acinetobacter species are increasingly isolated in the clinical setting and the environment. The aim of the present study was to analyze a genome database of 837 Acinetobacter spp. isolates, which included 798 non-baumannii Acinetobacter genomes, in order to define the concordance of classification and discriminatory power of 7-gene MLST, 53-gene MLST, and single-nucleotide polymorphism (SNPs) phylogenies.MethodsPhylogenies were performed on Pasteur Multilocus Sequence Typing (MLST) or ribosomal Multilocus Sequence Typing (rMLST) concatenated alleles, or SNPs extracted from core genome alignment.ResultsThe Pasteur MLST scheme was able to identify and genotype 72 species in the Acinetobacter genus, with classification results concordant with the ribosomal MLST scheme. The discriminatory power and genotyping reliability of the Pasteur MLST scheme were assessed in comparison to genome-wide SNP phylogeny on 535 non-baumannii Acinetobacter genomes assigned to Acinetobacter pittii, Acinetobacter nosocomialis, Acinetobacter seifertii, and Acinetobacter lactucae (heterotypic synonym of Acinetobacter dijkshoorniae), which were the most clinically relevant non-baumannii species of the A. baumannii group. The Pasteur MLST and SNP phylogenies were congruent at Robinson-Fould and Matching cluster tests and grouped genomes into four and three clusters in A. pittii, respectively, and one each in A. seifertii. Furthermore, A. lactucae genomes were grouped into one cluster within A. pittii genomes. The SNP phylogeny of A. nosocomialis genomes showed a heterogeneous population and did not correspond to the Pasteur MLST phylogeny, which identified two recombinant clusters. The antimicrobial resistance genes belonging to at least three different antimicrobial classes were identified in 91 isolates assigned to 17 distinct species in the Acinetobacter genus. Moreover, the presence of a class D oxacillinase, which is a naturally occurring enzyme in several Acinetobacter species, was found in 503 isolates assigned to 35 Acinetobacter species.ConclusionIn conclusion, Pasteur MLST phylogeny of non-baumannii Acinetobacter isolates coupled with in silico detection of antimicrobial resistance makes it important to study the population structure and epidemiology of Acinetobacter spp. isolates.

Summary: Background: The Invasive Respiratory Infection Surveillance (IRIS) Consortium was established to assess the impact of the COVID-19 pandemic on invasive diseases caused by Streptococcus pneumoniae, Haemophilus influenzae, Neisseria meningitidis, and Streptococcus agalactiae. We aimed to analyse the incidence and distribution of these diseases during the first 2 years of the COVID-19 pandemic compared to the 2 years preceding the pandemic. Methods: For this prospective analysis, laboratories in 30 countries and territories representing five continents submitted surveillance data from Jan 1, 2018, to Jan 2, 2022, to private projects within databases in PubMLST. The impact of COVID-19 containment measures on the overall number of cases was analysed, and changes in disease distributions by patient age and serotype or group were examined. Interrupted time-series analyses were done to quantify the impact of pandemic response measures and their relaxation on disease rates, and autoregressive integrated moving average models were used to estimate effect sizes and forecast counterfactual trends by hemisphere. Findings: Overall, 116 841 cases were analysed: 76 481 in 2018–19, before the pandemic, and 40 360 in 2020–21, during the pandemic. During the pandemic there was a significant reduction in the risk of disease caused by S pneumoniae (risk ratio 0·47; 95% CI 0·40–0·55), H influenzae (0·51; 0·40–0·66) and N meningitidis (0·26; 0·21–0·31), while no significant changes were observed for S agalactiae (1·02; 0·75–1·40), which is not transmitted via the respiratory route. No major changes in the distribution of cases were observed when stratified by patient age or serotype or group. An estimated 36 289 (95% prediction interval 17 145–55 434) cases of invasive bacterial disease were averted during the first 2 years of the pandemic among IRIS-participating countries and territories. Interpretation: COVID-19 containment measures were associated with a sustained decrease in the incidence of invasive disease caused by S pneumoniae, H influenzae, and N meningitidis during the first 2 years of the pandemic, but cases began to increase in some countries towards the end of 2021 as pandemic restrictions were lifted. These IRIS data provide a better understanding of microbial transmission, will inform vaccine development and implementation, and can contribute to health-care service planning and provision of policies. Funding: Wellcome Trust, NIHR Oxford Biomedical Research Centre, Spanish Ministry of Science and Innovation, Korea Disease Control and Prevention Agency, Torsten Söderberg Foundation, Stockholm County Council, Swedish Research Council, German Federal Ministry of Health, Robert Koch Institute, Pfizer, Merck, and the Greek National Public Health Organization.

The genus Bordetella includes environmental bacteria as well as human pathogens. Here, the authors present a large database of environmental and clinical Bordetella isolates and genome sequences, and develop genotyping systems to facilitate evolutionary and epidemiological studies.

ABSTRACT Clostridioides difficile remains a key cause of healthcare-associated infection, with multidrug-resistant (MDR) lineages causing high-mortality (≥20%) outbreaks. Cephalosporin treatment is a long-established risk factor, and antimicrobial stewardship is a key control. A mechanism underlying raised cephalosporin MICs has not been identified in C. difficile, but among other species, this is often acquired via amino acid substitutions in cell wall transpeptidases (penicillin binding proteins [PBPs]). Here, we investigated five C. difficile transpeptidases (PBP1 to PBP5) for recent substitutions, associated cephalosporin MICs, and co-occurrence with fluoroquinolone resistance. Previously published genome assemblies (n = 7,096) were obtained, representing 16 geographically widespread lineages, including healthcare-associated ST1(027). Recent amino acid substitutions were found within PBP1 (n = 50) and PBP3 (n = 48), ranging from 1 to 10 substitutions per genome. β-Lactam MICs were measured for closely related pairs of wild-type and PBP-substituted isolates separated by 20 to 273 single nucleotide polymorphisms (SNPs). Recombination-corrected phylogenies were constructed to date substitution acquisition. Key substitutions such as PBP3 V497L and PBP1 T674I/N/V emerged independently across multiple lineages. They were associated with extremely high cephalosporin MICs; 1 to 4 doubling dilutions >wild-type, up to 1,506 μg/mL. Substitution patterns varied by lineage and clade, showed geographic structure, and occurred post-1990, coincident with the gyrA and/or gyrB substitutions conferring fluoroquinolone resistance. In conclusion, recent PBP1 and PBP3 substitutions are associated with raised cephalosporin MICs in C. difficile. Their co-occurrence with fluoroquinolone resistance hinders attempts to understand the relative importance of these drugs in the dissemination of epidemic lineages. Further controlled studies of cephalosporin and fluoroquinolone stewardship are needed to determine their relative effectiveness in outbreak control. IMPORTANCE Fluoroquinolone and cephalosporin use in healthcare settings has triggered outbreaks of high-mortality, multidrug-resistant C. difficile infection. Here, we identify a mechanism associated with raised cephalosporin MICs in C. difficile comprising amino acid substitutions in two cell wall transpeptidase enzymes (penicillin binding proteins). The higher the number of substitutions, the greater the impact on phenotype. Dated phylogenies revealed that substitutions associated with raised cephalosporin and fluoroquinolone MICs were co-acquired immediately before clinically important outbreak strains emerged. PBP substitutions were geographically structured within genetic lineages, suggesting adaptation to local antimicrobial prescribing. Antimicrobial stewardship of cephalosporins and fluoroquinolones is an effective means of C. difficile outbreak control. Genetic changes associated with raised MIC may impart a “fitness cost” after antibiotic withdrawal. Our study therefore identifies a mechanism that may explain the contribution of cephalosporin stewardship to resolving outbreak conditions. However, due to the co-occurrence of raised cephalosporin MICs and fluoroquinolone resistance, further work is needed to determine the relative importance of each.

Acinetobacter baumannii is an opportunistic pathogen responsible for healthcare-associated infections (HAIs) and outbreaks. Antimicrobial resistance mechanisms and virulence factors allow it to survive and spread in the hospital environment. However, the molecular mechanisms of these traits and their association with international clones are frequently unknown in low- and middle-income countries. Here, we analyze the phenotype and genotype of seventy-six HAIs and outbreak-causing A. baumannii isolates from a Mexican hospital over ten years, with special attention to the carbapenem resistome and biofilm formation. The isolates belonged to the global international clone (IC) 2 and the Latin America endemic IC5 and were predominantly extensively drug-resistant (XDR). Oxacillinases were identified as a common source of carbapenem resistance. We noted the presence of the blaOXA-143-like family (not previously described in Mexico), the blaOXA-72 and the blaOXA-398 found in both ICs. A low prevalence of efflux pump overexpression activity associated with carbapenem resistance was observed. Finally, strong biofilm formation was found, and significant biofilm-related genes were identified, including bfmRS, csuA/BABCDE, pgaABCD and ompA. This study provides a comprehensive profile of the carbapenem resistome of A. baumannii isolates belonging to the same pulse type, along with their significant biofilm formation capacity. Furthermore, it contributes to a better understanding of their role in the recurrence of infection and the endemicity of these isolates in a Mexican hospital.

Expansion of quinolone-resistant Neisseria meningitidis clone ChinaCC4821-R1-C/B from sequence type (ST) 4821 clonal complex (CC4821) caused a serogroup shift from serogroup A to serogroup C invasive meningococcal disease (IMD) in China. To determine the relationship among globally distributed CC4821 meningococci, we analyzed whole-genome sequence data from 173 CC4821 meningococci isolated from 4 continents during 1972–2019. These meningococci clustered into 4 sublineages (1–4); sublineage 1 primarily comprised of IMD isolates (41/50, 82%). Most isolates from outside China (40/49, 81.6%) formed a distinct sublineage, the Europe–USA cluster, with the typical strain designation B:P1.17-6,23:F3-36:ST-3200(CC4821), harboring mutations in penicillin-binding protein 2. These data show that the quinolone-resistant clone ChinaCC4821-R1-C/B has expanded to other countries. The increasing distribution worldwide of serogroup B CC4821 raises the concern that CC4821 has the potential to cause a pandemic that would be challenging to control, despite indirect evidence that the Trumenba vaccine might afford some protection.

It is unclear how gut-dwelling E. coli bacteria often emerge to cause systemic infection in chickens. Here, Mageiros et al. use population genomics and pangenome-wide association studies to identify genetic elements associated with pathogenicity in avian E. coli.

Horizontal gene transfer (HGT) can allow traits that have evolved in one bacterial species to transfer to another. This has potential to rapidly promote new adaptive trajectories such as zoonotic transfer or antimicrobial resistance. However, for this to occur requires gaps to align in barriers to recombination within a given time frame. Chief among these barriers is the physical separation of species with distinct ecologies in separate niches. Within the genus Campylobacter, there are species with divergent ecologies, from rarely isolated single-host specialists to multihost generalist species that are among the most common global causes of human bacterial gastroenteritis. Here, by characterizing these contrasting ecologies, we can quantify HGT among sympatric and allopatric species in natural populations. Analyzing recipient and donor population ancestry among genomes from 30 Campylobacter species, we show that cohabitation in the same host can lead to a six-fold increase in HGT between species. This accounts for up to 30% of all SNPs within a given species and identifies highly recombinogenic genes with functions including host adaptation and antimicrobial resistance. As described in some animal and plant species, ecological factors are a major evolutionary force for speciation in bacteria and changes to the host landscape can promote partial convergence of distinct species through HGT.

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.

ABSTRACT Clostridium perfringens is a spore-forming anaerobic pathogen responsible for a variety of histotoxic and intestinal infections in humans and animals. High-resolution genotyping aiming to identify bacteria at strain level has become increasingly important in modern microbiology to understand pathogen transmission pathways and to tackle infection sources. This study aimed at establishing a publicly available genome-wide multilocus sequence-typing (MLST) scheme for C. perfringens. A total of 1,431 highly conserved core genes (1.34 megabases; 50% of the reference genome genes) were indexed for a core genome-based MLST (cgMLST) scheme for C. perfringens. The scheme was applied to 282 ecologically and geographically diverse genomes, showing that the genotyping results of cgMLST were highly congruent with the core genome-based single-nucleotide-polymorphism typing in terms of resolution and tree topology. In addition, the cgMLST provided a greater discrimination than classical MLST methods for C. perfringens. The usability of the scheme for outbreak analysis was confirmed by reinvestigating published outbreaks of C. perfringens-associated infections in the United States and the United Kingdom. In summary, a publicly available scheme and an allele nomenclature database for genomic typing of C. perfringens have been established and can be used for broad-based and standardized epidemiological studies. IMPORTANCE Global epidemiological surveillance of bacterial pathogens is enhanced by the availability of standard tools and sharing of typing data. The use of whole-genome sequencing has opened the possibility for high-resolution characterization of bacterial strains down to the clonal and subclonal levels. Core genome multilocus sequence typing is a robust system that uses highly conserved core genes for deep genotyping. The method has been successfully and widely used to describe the epidemiology of various bacterial species. Nevertheless, a cgMLST typing scheme for Clostridium perfringens is currently not publicly available. In this study, we (i) developed a cgMLST typing scheme for C. perfringens, (ii) evaluated the performance of the scheme on different sets of C. perfringens genomes from different hosts and geographic regions as well as from different outbreak situations, and, finally, (iii) made this scheme publicly available supported by an allele nomenclature database for global and standard genomic typing.

Summary: Background: Streptococcus pneumoniae, Haemophilus influenzae, and Neisseria meningitidis, which are typically transmitted via respiratory droplets, are leading causes of invasive diseases, including bacteraemic pneumonia and meningitis, and of secondary infections subsequent to post-viral respiratory disease. The aim of this study was to investigate the incidence of invasive disease due to these pathogens during the early months of the COVID-19 pandemic. Methods: In this prospective analysis of surveillance data, laboratories in 26 countries and territories across six continents submitted data on cases of invasive disease due to S pneumoniae, H influenzae, and N meningitidis from Jan 1, 2018, to May, 31, 2020, as part of the Invasive Respiratory Infection Surveillance (IRIS) Initiative. Numbers of weekly cases in 2020 were compared with corresponding data for 2018 and 2019. Data for invasive disease due to Streptococcus agalactiae, a non-respiratory pathogen, were collected from nine laboratories for comparison. The stringency of COVID-19 containment measures was quantified using the Oxford COVID-19 Government Response Tracker. Changes in population movements were assessed using Google COVID-19 Community Mobility Reports. Interrupted time-series modelling quantified changes in the incidence of invasive disease due to S pneumoniae, H influenzae, and N meningitidis in 2020 relative to when containment measures were imposed. Findings: 27 laboratories from 26 countries and territories submitted data to the IRIS Initiative for S pneumoniae (62 837 total cases), 24 laboratories from 24 countries submitted data for H influenzae (7796 total cases), and 21 laboratories from 21 countries submitted data for N meningitidis (5877 total cases). All countries and territories had experienced a significant and sustained reduction in invasive diseases due to S pneumoniae, H influenzae, and N meningitidis in early 2020 (Jan 1 to May 31, 2020), coinciding with the introduction of COVID-19 containment measures in each country. By contrast, no significant changes in the incidence of invasive S agalactiae infections were observed. Similar trends were observed across most countries and territories despite differing stringency in COVID-19 control policies. The incidence of reported S pneumoniae infections decreased by 68% at 4 weeks (incidence rate ratio 0·32 [95% CI 0·27–0·37]) and 82% at 8 weeks (0·18 [0·14–0·23]) following the week in which significant changes in population movements were recorded. Interpretation: The introduction of COVID-19 containment policies and public information campaigns likely reduced transmission of S pneumoniae, H influenzae, and N meningitidis, leading to a significant reduction in life-threatening invasive diseases in many countries worldwide. Funding: Wellcome Trust (UK), Robert Koch Institute (Germany), Federal Ministry of Health (Germany), Pfizer, Merck, Health Protection Surveillance Centre (Ireland), SpID-Net project (Ireland), European Centre for Disease Prevention and Control (European Union), Horizon 2020 (European Commission), Ministry of Health (Poland), National Programme of Antibiotic Protection (Poland), Ministry of Science and Higher Education (Poland), Agencia de Salut Pública de Catalunya (Spain), Sant Joan de Deu Foundation (Spain), Knut and Alice Wallenberg Foundation (Sweden), Swedish Research Council (Sweden), Region Stockholm (Sweden), Federal Office of Public Health of Switzerland (Switzerland), and French Public Health Agency (France).

Campylobacter is the leading bacterial cause of gastroenteritis worldwide and its incidence is especially high in low- and middle-income countries (LMIC). Disease epidemiology in LMICs is different compared to high income countries like the USA or in Europe. Children in LMICs commonly have repeated and chronic infections even in the absence of symptoms, which can lead to deficits in early childhood development. In this study, we sequenced and characterized C. jejuni (n = 62) from a longitudinal cohort study of children under the age of 5 with and without diarrheal symptoms, and contextualized them within a global C. jejuni genome collection. Epidemiological differences in disease presentation were reflected in the genomes, specifically by the absence of some of the most common global disease-causing lineages. As in many other countries, poultry-associated strains were likely a major source of human infection but almost half of local disease cases (15 of 31) were attributable to genotypes that are rare outside of Peru. Asymptomatic infection was not limited to a single (or few) human adapted lineages but resulted from phylogenetically divergent strains suggesting an important role for host factors in the cryptic epidemiology of campylobacteriosis in LMICs.

Staphylococcus epidermidis is carried asymptomatically by virtually all humans but is also a major cause of nosocomial infection. Here, the authors study 141 isolates from healthy carriage and 274 isolates from clinical infections, and identify genes and genetic elements associated with pathogenicity.

Abstract Background Chlamydia trachomatis (Ct) plasmid has been shown to encode genes essential for infection. We evaluated the population structure of Ct using whole-genome sequence data (WGS). In particular, the relationship between the Ct genome, plasmid and disease was investigated. Results WGS data from 157 Ct isolates deposited in the Chlamydiales pubMLST database (http://pubMLST.org/chlamydiales/) were annotated with 902 genes including the core and accessory genome. Plasmid associated genes were annotated and a plasmid MLST scheme was defined allowing plasmid sequence types to be determined. Plasmid allelic variation was investigated. Phylogenetic relationships were examined using the Genome Comparator tool available in pubMLST. Phylogenetic analyses identified four distinct Ct core genome clusters and six plasmid clusters, with a strong association between the chromosomal genotype and plasmid. This in turn was linked to ompA genovars and disease phenotype. Horizontal genetic transfer of plasmids was observed for three urogenital-associated isolates, which possessed plasmids more commonly found in isolates resulting from ocular infections. The pgp3 gene was identified as the most polymorphic plasmid gene and pgp4 was the most conserved. Conclusion A strong association between chromosomal genome, plasmid type and disease was observed, consistent with previous studies. This suggests co-evolution of the Ct chromosome and their plasmids, but we confirmed that plasmid transfer can occur between isolates. These data provide a better understanding of the genetic diversity occurring across the Ct genome in association with the plasmid content.

A Correction to this paper has been published: https://doi.org/10.1038/s41467-021-22238-5

In Western Australia, Neisseria meningitidis serogroup W clonal complex 11 became the predominant cause of invasive meningococcal disease in 2016. We used core-genome analysis to show emergence of a penicillin-resistant clade that had the penA_253 allele. This new penicillin-resistant clade might affect treatment regimens for this disease.

Acinetobacter species assigned to the Acinetobacter calcoaceticus-baumannii (Acb) complex, are Gram-negative bacteria responsible for a large number of human infections. The population structure of Acb has been studied using two 7-gene MLST schemes, introduced by Bartual and coworkers (Oxford scheme) and by Diancourt and coworkers (Pasteur scheme). The schemes have three genes in common but underlie two coexisting nomenclatures of sequence types and clonal complexes, which complicates communication on A. baumannii genotypes. The aim of this study was to compare the characteristics of the two schemes to make a recommendation about their usage. Using genome sequences of 730 strains of the Acb complex, we evaluated the phylogenetic congruence of MLST schemes, the correspondence between sequence types, their discriminative power and genotyping reliability from genomic sequences. In silico ST re-assignments highlighted the presence of a second copy of the Oxford gdhB locus, present in 553/730 genomes that has led to the creation of artefactual profiles and STs. The reliability of the two MLST schemes was tested statistically comparing MLST-based phylogenies to two reference phylogenies (core-genome genes and genome-wide SNPs) using topology-based and likelihood-based tests. Additionally, each MLST gene fragment was evaluated by correlating the pairwise nucleotide distances between each pair of genomes calculated on the core-genome and on each single gene fragment. The Pasteur scheme appears to be less discriminant among closely related isolates, but less affected by homologous recombination and more appropriate for precise strain classification in clonal groups, which within this scheme are more often correctly monophyletic. Statistical tests evaluate the tree deriving from the Oxford scheme as more similar to the reference genome trees. Our results, together with previous work, indicate that the Oxford scheme has important issues: gdhB paralogy, recombination, primers sequences, position of the genes on the genome. While there is no complete agreement in all analyses, when considered as a whole the above results indicate that the Pasteur scheme is more appropriate for population biology and epidemiological studies of A. baumannii and related species and we propose that it should be the scheme of choice during the transition toward, and in parallel with, core genome MLST.

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.

Introduction and aimsSince 2013 MenC and MenW disease incidence and associated mortality rates have increased in the Republic of Ireland. From 2002/2003 to 2012/2013, the average annual MenC incidence was 0.08/100,000, which increased to 0.34/100,000 during 2013/2014 to 2017/18, peaking in 2016/17 (0.72/100,000) with an associated case fatality rate (CFR) of 14.7%. MenW disease incidence has increased each year from 0.02/100,000 in 2013/2014, to 0.29/100,000 in 2017/18, with an associated CFR of 28.6%. We aimed to characterise and relate recent MenC isolates to the previously prevalent MenC:cc11 ET-15 clones, and also characterise and relate recent MenW isolates to the novel 'Hajj' clones.MethodsUsing WGS we characterised invasive (n = 74, 1997/98 to 2016/17) and carried (n = 16, 2016/17) MenC isolates, and invasive (n = 18, 2010/11 to 2016/17) and carried (n = 15, 2016/17) MenW isolates. Genomes were assembled using VelvethOptimiser and stored on the PubMLST Neisseria Bacterial Isolate Genome Sequence Database. Isolates were compared using the cgMLST approach.ResultsMost MenC and MenW isolates identified were cc11. A single MenC:cc11 sub-lineage contained the majority (68%, n = 19/28) of recent MenC:cc11 disease isolates and all carried MenC:cc11 isolates, which were interspersed and distinct from the historically significant ET-15 clones. MenW:cc11 study isolates clustered among international examples of both the original UK 2009 MenW:cc11, and novel 2013 MenW:cc11clones.ConclusionsWe have shown that the majority of recent MenC disease incidence was caused by strain types distinct from the MenC:cc11 ET-15 clone of the late 1990s, which still circulate but have caused only sporadic disease in recent years. We have identified that the same aggressive MenW clone now established in several other European countries, is endemic in the RoI and responsible for the recent MenW incidence increases. This data informed the National immunisation Advisory Committee, who are currently deliberating a vaccine policy change to protect teenagers.

Single-strain outbreaks of Streptococcus pyogenes infections are common and often go undetected. In 2013, two clusters of invasive group A Streptococcus (iGAS) infection were identified in independent but closely located care homes in Oxfordshire, United Kingdom. Investigation included visits to each home, chart review, staff survey, microbiologic sampling, and genome sequencing. S. pyogenes emm type 1.0, the most common circulating type nationally, was identified from all cases yielding GAS isolates. A tailored whole-genome reference population comprising epidemiologically relevant contemporaneous isolates and published isolates was assembled. Data were analyzed independently using whole-genome multilocus sequencing and single-nucleotide polymorphism analyses. Six isolates from staff and residents of the homes formed a single cluster that was separated from the reference population by both analytical approaches. No further cases occurred after mass chemoprophylaxis and enhanced infection control. Our findings demonstrate the ability of 2 independent analytical approaches to enable robust conclusions from nonstandardized whole-genome analysis to support public health practice.

The PubMLST.org website hosts a collection of open-access, curated databases that integrate population sequence data with provenance and phenotype information for over 100 different microbial species and genera. Although the PubMLST website was conceived as part of the development of the first multi-locus sequence typing (MLST) scheme in 1998 the software it uses, the Bacterial Isolate Genome Sequence database (BIGSdb, published in 2010), enables PubMLST to include all levels of sequence data, from single gene sequences up to and including complete, finished genomes. Here we describe developments in the BIGSdb software made from publication to June 2018 and show how the platform realises microbial population genomics for a wide range of applications. The system is based on the gene-by-gene analysis of microbial genomes, with each deposited sequence annotated and curated to identify the genes present and systematically catalogue their variation. Originally intended as a means of characterising isolates with typing schemes, the synthesis of sequences and records of genetic variation with provenance and phenotype data permits highly scalable (whole genome sequence data for tens of thousands of isolates) means of addressing a wide range of functional questions, including: the prediction of antimicrobial resistance; likely cross-reactivity with vaccine antigens; and the functional activities of different variants that lead to key phenotypes. There are no limitations to the number of sequences, genetic loci, allelic variants or schemes (combinations of loci) that can be included, enabling each database to represent an expanding catalogue of the genetic variation of the population in question. In addition to providing web-accessible analyses and links to third-party analysis and visualisation tools, the BIGSdb software includes a RESTful application programming interface (API) that enables access to all the underlying data for third-party applications and data analysis pipelines.

ABSTRACT Neisseria meningitidis is a common cause of bacterial meningitis in children and young adults worldwide. The 4CMenB vaccine (Bexsero), developed to combat meningococcal serogroup B (MenB) disease, contains subcapsular antigens that may induce immunity against strains of N. meningitidis, regardless of serogroup. Owing to differential levels of expression and peptide diversity in vaccine antigens across meningococcal strains, the meningococcal antigen typing system (MATS) was developed to estimate the potential MenB strain coverage of 4CMenB. Prior to introducing the 4CMenB vaccine into routine use, we sought to estimate the potential 4CMenB coverage against invasive MenB strains isolated in the Republic of Ireland (RoI) over four consecutive epidemiological years. MATS was applied to a panel of 105 invasive MenB strains isolated during July 2009 to June 2013. Sequence data characterizing the multilocus sequence typing (MLST) alleles and the major 4CMenB target peptides were extracted from isolate genome sequence data, hosted in the Bacterial Isolate Sequencing database (BIGSdb). MATS data indicated that 4CMenB may induce protective immunity against 69.5% (95% confidence interval [CI95%], 64.8% to 84.8%) of circulating MenB strains. Estimated coverage was highest against the most prevalent disease-causing lineage, cc41/44, where the most frequently observed sequence types, ST-154 and ST-41 (21% of isolates, collectively), were typically covered by three antigens. No significant temporal trends were observed. Overall, these data provide a baseline of strain coverage prior to the introduction of 4CMenB and indicate that a decrease in invasive meningococcal disease (IMD) is predicted following the introduction of 4CMenB into the routine infant immunization schedule in the RoI. IMPORTANCE The meningococcal antigen typing system (MATS) is an enzyme-linked immunosorbent assay (ELISA) that measures both the levels of expression and the immune reactivity of the three recombinant 4CMenB antigens. Together with PorA variable-region sequence data, this system provides an estimation of how susceptible MenB isolates are to killing by 4CMenB vaccine-induced antibodies. Assays based on subcapsular antigen phenotype analyses, such as MATS, are important in situations where conventional vaccine coverage estimations are not possible. Subcapsular antigens are typically highly diverse across strains, and vaccine coverage estimations would require unfeasibly large efficacy trials and screening of an exhaustive strain panel for antibody functional activity. Here, MATS was applied to all invasive meningococcal serogroup B (MenB) strains isolated over four consecutive epidemiological years (n = 105) and predicted reasonably high 4CMenB vaccine coverage in the Republic of Ireland.

Dichelobacter nodosus (D. nodosus) is the causative pathogen of ovine footrot, a disease that has a significant welfare and financial impact on the global sheep industry. Previous studies into the phylogenetics of D. nodosus have focused on Australia and Scandinavia, meaning the current diversity in the United Kingdom (U.K.) population and its relationship globally, is poorly understood. Numerous epidemiological methods are available for bacterial typing; however, few account for whole genome diversity or provide the opportunity for future application of new computational techniques. Multilocus sequence typing (MLST) measures nucleotide variations within several loci with slow accumulation of variation to enable the designation of allele numbers to determine a sequence type. The usage of whole genome sequence data enables the application of MLST, but also core and whole genome MLST for higher levels of strain discrimination with a negligible increase in experimental cost. An MLST database was developed alongside a seven loci scheme using publically available whole genome data from the sequence read archive. Sequence type designation and strain discrimination was compared to previously published data to ensure reproducibility. Multiple D. nodosus isolates from U.K. farms were directly compared to populations from other countries. The U.K. isolates define new clades within the global population of D. nodosus and predominantly consist of serogroups A, B and H, however serogroups C, D, E, and I were also found. The scheme is publically available at https://pubmlst.org/dnodosus/.

Acinetobacter baumannii is a troublesome opportunistic pathogen with a high capacity for clonal dissemination. We announce the establishment of a database for the ampC locus in A. baumannii, in which novel ampC alleles are differentiated based on the occurrence of ≥ 1 nucleotide change, regardless of whether it is silent or missense. The database is openly accessible at the pubmlst platform for A. baumannii (http://pubmlst.org/abaumannii/). Forty-eight distinctive alleles of the ampC locus have so far been identified and deposited in the database. Isolates from clonal complex 1 (CC1), according to the Pasteur multilocus sequence typing scheme, had a variety of the ampC locus alleles, including alleles 1, 3, 4, 5, 6, 7, 8, 13, 14, 17, and 18. On the other hand, isolates from CC2 had the ampC alleles 2, 3, 19, 20, 21, 22, 23, 24, 26, 27, 28, and 46. Allele 3 was characteristic for sequence types ST3 or ST32. The ampC alleles 10, 16, and 25 were characteristic for CC10, ST16, and CC25, respectively. Our study points out that novel gene databases, in which alleles are numbered based on differences in their nucleotide identities, should replace traditional records that use amino acid substitutions to define new alleles.

Pathogenic Neisseria meningitidis isolates contain a polysaccharide capsule that is the main virulence determinant for this bacterium. Thirteen capsular polysaccharides have been described, and nuclear magnetic resonance spectroscopy has enabled determination of the structure of capsular polysaccharides responsible for serogroup specificity. Molecular mechanisms involved in N. meningitidis capsule biosynthesis have also been identified, and genes involved in this process and in cell surface translocation are clustered at a single chromosomal locus termed cps. The use of multiple names for some of the genes involved in capsule synthesis, combined with the need for rapid diagnosis of serogroups commonly associated with invasive meningococcal disease, prompted a requirement for a consistent approach to the nomenclature of capsule genes. In this report, a comprehensive description of all N. meningitidis serogroups is provided, along with a proposed nomenclature, which was presented at the 2012 XVIIIth International Pathogenic Neisseria Conference.

Bacterial species comprise related genotypes that can display divergent phenotypes with important clinical implications. Staphylococcus epidermidis is a common cause of nosocomial infections and, critical to its pathogenesis, is its ability to adhere and form biofilms on surfaces, thereby moderating the effect of the host's immune response and antibiotics. Commensal S. epidermidis populations are thought to differ from those associated with disease in factors involved in adhesion and biofilm accumulation. We quantified the differences in biofilm formation in 98 S. epidermidis isolates from various sources, and investigated population structure based on ribosomal multilocus typing (rMLST) and the presence/absence of genes involved in adhesion and biofilm formation. All isolates were able to adhere and form biofilms in in vitro growth assays and confocal microscopy allowed classification into 5 biofilm morphotypes based on their thickness, biovolume and roughness. Phylogenetic reconstruction grouped isolates into three separate clades, with the isolates in the main disease associated clade displaying diversity in morphotype. Of the biofilm morphology characteristics, only biofilm thickness had a significant association with clade distribution. The distribution of some known adhesion-associated genes (aap and sesE) among isolates showed a significant association with the species clonal frame. These data challenge the assumption that biofilm-associated genes, such as those on the ica operon, are genetic markers for less invasive S. epidermidis isolates, and suggest that phenotypic characteristics, such as adhesion and biofilm formation, are not fixed by clonal descent but are influenced by the presence of various genes that are mobile among lineages.